Y-W方程的L-D递推算法和伯格递推算法的比较.pptx

Y-W方程的L-D递推算法和伯格递推算法 的比较AR模型参数与自相关函数的关系Y-W 方程的矩阵表示例: p=3 时的Y-W 方程 若已知Ry[n] ，由Y-W方程解出各参数a3(1), a3 (2), a3 (3), s2。则可由AR模型参数获得功率谱Py(?)的估计值。Y-W方程的L-D递推算法 一阶Y-W方程的解解此方程得 Y-W方程的L-D递推算法 二阶Y-W方程的解Y-W方程的L-D递推算法 p阶Y-W方程的递推解Y-W方程的L-D递推算法 L-D算法估计功率谱的步骤(1) 计算自相关函数的估计值(2) 求解一阶模型参数关函数的估计值Y-W方程的L-D递推算法 L-D算法估计功率谱的步骤(3)由递推算法求解p阶模型参数(4) 求出功率谱估计利用L-D算法估计频谱的MATLAB函数[Pxx,f] = pyulear(x,p,NFFT,Fs)x：进行功率谱估计的输入有限长序列；p:模型的阶数NFFT：DFT的点数；Fs ：绘制功率谱曲线的抽样频率，默认值为1；Pxx：功率谱估计值； f：Pxx值所对应的频率点例：利用L-D算法进行谱估计一序列含有白噪声和两个频率的余弦信号，利用L-D算法估计该序列的功率谱。N=128;p=40;NFFT=2048;Fs=2; n=0:N-1;randn(state,0);x=cos(0.3*pi*n)+cos(0.4*pi*n)+randn(size(n));[P,f]=periodogram(x,[],NFFT,2) ;[Py,fy]=pyulear(x,p,NFFT,2);subplot(211);plot(f,10*log(P));grid;title( Periodogram );axis([0 1 -60 60]);subplot(212);plot(fy,10*log(Py));grid;title( Yule);axis([0 1 -60 60]); Periodogram 6040200-20-40-6000.10.20.30.40.50.60.70.80.91 Yule6040200-20-40-6000.10.20.30.40.50.60.70.80.91谱估计结果——p=40,N=128伯格(Burg)递推算法 L-D算法缺点： 在计算相关函数估计时，对N个观测数据以外的数据作零的假设，故谱估计误差较大。 伯格(Burg)递推算法基本思想： 直接从观测的数据利用线性预测器的前向和后向预测的总均方误差之和为最小的准则来估计反射系数，进而通过L-D算法的递推公式求出AR模型的优化参数。伯格(Burg)递推算法前向预测误差的递推公式 2阶前向预测误差 1阶后向预测误差 L-D算法的递推公式 K2=a2(2) 2阶预测器的反射系数 伯格(Burg)递推算法预测误差的递推公式 一般地同理可得后向预测误差的递推公式 Kp=ap(p)为p阶预测器的反射系数。 伯格(Burg)递推算法预测误差滤波器的格形结构 伯格(Burg)递推算法反射系数Kp 的确定前向和后向均方预测误差的总和为 由可得伯格(Burg)递推算法 伯格(Burg)递推算法步骤(1) 确定初始条件(2) 从p=1开始迭代计算:计算AR模型参数递推p阶均方误差伯格(Burg)递推算法 伯格(Burg)递推算法步骤(3) 递推高一阶前、后向预测误差(4) 若阶数小于p，则阶数加1，回到步骤(2)进行下一次迭代，直到达到预定阶数p。(5) 估计功率谱Burg算法估计频谱的MATLAB函数[Pxx,f] = pburg (x,p,NFFT,Fs)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]);L-D6040200-2000.10.20.30.40.50.60.70.80.91Burg6040200-2000.10.20.30.40.50.60.70.80.91AR模型阶数p=50 的谱估计结果L-D6