直接法求解循环谱密度，画出3D图，并与参考文献中的代码进行比较。