%Levinson
clear?all;???????????
N=128;??????????
P=60;????????????
f1=0.2;
f2=0.213;??????
k=1:1:N;?
sum1=0;
sum2=0;
snr=-15;

s=15;
d=1:1:s;
f1e=zeros（1s）;
f2e=zeros（1s）;
MSE1=zeros（1s）;
MSE2=zeros（1s）;


snr=-12;
x1=10*sin（2*pi*f1*k+pi/3）+5*sin（2*pi*f2*k+pi/4）;
xn=awgn（x1snr）;

%計算自相關函數?
Rx=[];?
for?m=1:N?????
????sum1=0;?
?for?n=1:N-（m-1）?????????
????????sum1=sum1+xn（n）*xn（n+m-1）;?
?end
????Rx（m）=sum1/N;?
end
?

%列文森遞推法
a（11）=-Rx（2）/Rx（1）;?
rou（1）=Rx（1）+a（11）*Rx（2）;?
for?p=2:P??????
????sum2=0;?
????for?k=1:p-1?????????
????????sum2=sum2+a（p-1k）*Rx（p-k+1）;
????end
?