直接解Yule-Walker方程法估计功率谱

%Yule Walker
N=128;
n=0:1:N-1;
P=80;

snr=1;
w=2*pi*128/128;
x1=10*sin（2*pi*0.2*n+pi/3）+2*sin（2*pi*0.213*n+pi/4）;
xn=awgn（x1snr）;

%自相关函数
sum1=0;
Rx=zeros（1N+1）;
for m=1:N+1;
sum1=0;
for n=1:N-m+1;
sum1=sum1+xn（n）.*xn（n+m-1）;
end
Rx（m）=sum1/N;
end

%构建P*P矩阵
R=zeros（PP）;
for p=1:1:P;
    R（pp:P）=Rx（1:P-p+1）;
end