clc
clear?all
clf

N=128;??
n=1:N;???????
w（n）=randn（1N）;?????%產(chǎn)生零均值方差為1的高斯白噪聲
xn=sqrt（20）*sin（2*pi*0.2*n）+sqrt（2）*sin（2*pi*0.213*n）+w（n）;
a=zeros（N-1N-1）;
p0=sum（xn.^2）/N;
f0=xn;???%%%
g0=xn;???%%%
Km=0;
P=zeros（1N-1）;
f=zeros（N-1N-1）;
g=zeros（N-1N-1）;
zi=0;
mu=0;
%?w=0:2*pi;
m=1;
for?n=m+1:N
????????????zi=zi+f0（n）*g0（n-1）;
????????????mu=mu+f0（n）*f0（n）+g0（n-1）*g0（n-1）;
end
????????Km=-2*zi/mu;
%?????????a（10）=1;
????????a（m1）=Km;
????????P（m）=（1-Km^2）*p0;??????%%%
????????for?n=2:N
????????f（1n）=f0（n）+Km*g0（n-1）;
????????g（1n）=conj（Km）*f0（n）+g0（n-1）;
????????end

fo