北京交通大学数字信号处理方法研究生课程-现代功率谱估计.pptVIP

(3) 递推高一阶前、后向预测误差 (4) 若阶数小于p，则阶数加1，回到步骤(2)进行下一次迭代，直到达到预定阶数p。 (5) 估计功率谱 伯格(Burg)递推算法步骤 伯格(Burg)递推算法 [Pxx,f] = pburg (x,p,NFFT,Fs) Burg算法估计频谱的MATLAB函数 x：进行功率谱估计的输入有限长序列； p: 模型的阶数 NFFT：DFT的点数； Fs ：绘制功率谱曲线的抽样频率，默认值为1； Pxx：功率谱估计值； f：Pxx值所对应的频率点 利用Burg法进行谱估计程序 N=512;NFFT=1024;Fs=2;p=40; n=0:N-1;randn(state,0); x=cos(0.3*pi*n)+cos(0.32*pi*n)+randn(size(n)); [P,f]=pyulear(x,p,NFFT,2); [Pw, f2]=pburg(x,p,NFFT,2); subplot(211);plot(f,10*log(P));grid;title(L-D); axis([0 1 -30 60]); subplot(212);plot(f2,10*log(Pw));grid;title(Burg); axis([0 1 -30 60]); AR模型阶数p=50 的谱估计结果 0 0.1 0.2 0.3 0.4 0.5 0.6 0.