%正弦干擾陷波LMS算法
clear all

%正弦信號(hào)產(chǎn)生 噪聲s 正弦信號(hào)v
%輸入信號(hào)x=v+s
[s,pv]=rands(8,pi/3);
v=randv(pv);
x=v+s;

%步長(zhǎng)u，期望d
N=10;
d(1)=0;
for i=1:100
    d(i)=x(i);
end
u=0.00028;%u應(yīng)再計(jì)算  收斂范圍
%LMS算法
[w,y,EE,e,r]=LMS(x,d,u,N);


%畫(huà)時(shí)域波形
%正弦信號(hào)s(n)
figure;
subplot(4,1,1);
plot(s);
title('正弦信號(hào)s(n)');

%噪聲v(n)
subplot(4,1,2);
plot(v);
axis([0 100 -10 10]);
title('噪聲v(n)');

%輸入信號(hào)x(n)=s(n)+v(n)
subplot(4,1,3);
plot(x);
title('輸入信號(hào)x(n)=s(n)+v(n)');


%畫(huà)輸出信號(hào)e(n)
subplot(4,1,4);
plot(e);
title('輸出信號(hào)e(n)');


%畫(huà)均方誤差曲線
figure;
subplot(2,1,1);
plot(EE);
title('誤差曲線');


X=fft(x,256);
E=fft(e,256);
% H=EE/X;
for i=1:256
    H(i)=abs(E(i)/X(i));
end

%畫(huà)|H(w)|
subplot(2,1,2);
plot(H);
title('|H(w)|');