answer1.asv

正弦干扰陷波LMS算法

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

%正弦信号产生 噪声s 正弦信号v
%输入信号x=v+s
[s,pv]=rands(8,pi/3);
v=randv(pv);
x=v+s;

%步长u，期望d
N=10;
d(1)=0;
for i=1:100
    d(i)=x(i);
end
u=0.00028;%u应再计算  收敛范围
%LMS算法
[w,y,EE,e,r]=LMS(x,d,u,N);


%画时域波形
%正弦信号s(n)
figure;
subplot(4,1,1);
plot(s);
title('正弦信号s(n)');

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

%输入信号x(n)=s(n)+v(n)
subplot(4,1,3);
plot(x);
title('输入信号x(n)=s(n)+v(n)');


%画输出信号e(n)
subplot(4,1,4);
plot(e);
title('输出信号e(n)');


%画均方误差曲线
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

%画|H(w)|
subplot(2,1,2);
plot(H);
title('|H(w)|');