Yule-Walker法进行谱估计

Yule-Walker法的准则是是前向预测误差功率相对AR参数最小，该算法需要用到自相关函数序列，是所有AR参数求解方法中最简单的一种。对于长数据记录，因为可以得到较好的自相关估值，所以能够得到良好的谱估计。

clear all
N=128;
n=1:N;
p=55;
M=55;
x=sqrt（20）*sin（2*pi*0.2*n）+sqrt（2）*sin（2*pi*0.213*n）+randn（1N）;
subplot（211）;
plot（nx（n））;
axis（[0200-2020]）;
r0=0;
for n=1:N
    r0=r0+x（n）.^2;
end
r0=r0/N;
r=zeros（551）;
for m=1:p
    for n=1:N-m
        r（m）=r（m）+x（n）.*x（n+m）;
    end
    r（m）=r（m）/N;
end
q0=r0;
a（11）=-r（1）/q0;
k（1）=a（11