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