u=0.00025;
%算法的收斂性由收斂因子u決定，只有u在理論要求的范圍內時，算法才收斂，否則發散。
%在u值使算法收斂的范圍內，u值越大，收斂速度越快，u值越小，收斂速度越慢。
N=200;fo=100;f=1000;fai=pi/6;C=10;noise=1*randn（1N）;
pow1=（noise.*noise）;????????????????%pow1相當于randn（12*N）
power1=0;
for?k=1:N?
????power1=power1+pow1（1k）;
end
power1=power1/N;
for?k=1:N?
????signal（k）=2*cos（2*pi*k*fo/f）;
end
power2=sqrt（2）/2*2;?????????????????%有效值
INPUT_SNR=10*log（power2/power1）?????%信噪比
d=noise+signal;?????????????????????%合成信號
for?k=1:N???????????????????????????%參考信號
????x1（k）=C*cos（2*pi*k*fo/f+fai）;
????x2（k）=C*cos（2*pi*k*fo/f+fai+pi/2）;
end
W1=0.5*ones（1N+1）;
W2=0.5*ones（1N+1）;
y=zeros（1N）;
e=zeros（1N）;
for?j=1:N???????????????????????????%LMS（線形均方估計）算法，即最小均方估計
????y（j）=W1（j）*x1（j）+W2（j）*x2