用Levinson方法实现功率谱密度的估计。在N=1024条件下调试程序过程中用Levinson-Durbin快速递推法得出的图形不稳定，每次调试结果都不大相同，AR模型的阶次应选择在N/3<P<N/2之间，式中的N是记录的数据长度，这样可以得到高分辨率的估计谱，并且很少出现伪峰。-Levinson method used to achieve power spectral density estimate. Commissioning conditions under N = 10