Levinson快速递推法估计功率谱

%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