Levinson-Durbin迭代算法求AR模型参数

clear all;
clc;
close all;

fs=1;
f1=0.15;
f2=0.1700;
f3=0.2600;
t=0:499;
N=500;

v=（randn（1N）+j*rand（1N））/sqrt（2）;
%长500的复高斯白噪声序列
A1=sqrt（10^3）;
A2=sqrt（10^3）;
A3=sqrt（10^2.7）;
%产生观测信号长500 

s1=A1*exp（i*2*pi*f1*t）;
s2=A2*exp（i*2*pi*f2*t）;
s3=A3*exp（i*2*pi*f3*t）;
u=s1+s2+s3+v;

%-----------3.17（3）------------------
u_LD=u（1:256）;
p=16;
r0=xcorr（u_LD16‘biased‘）;
r=r0（p+1:2*p+1）;


a=zeros（1616）;
var=zeros（116）;
a（11）=-r（2）/r（1）;
var（1）=r（1）-abs（r（2））*abs（r（2））/r（1）;

k=zeros（116）;
for m=2:16
    temp=0;
    for i=1:m-1
        temp=temp