%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