lms_nois_rejection.m

最小均方自适应滤波噪声抑制

%LMS算法在噪声抑制中的应用。在仿真中，原始信号选为sin((0.2*pi)*t),
%噪声信号采用标准白噪声，延迟D=1，收敛因子分别是u=0.001和u=0.3
%LMS_nois_rejection

clear;clc;
N=100;
t=0:N-1;
s=sin(2*pi*t/10); %原始信号
figure(1)
subplot(2,1,1)
plot(t,s)
title('原始信号')

x=s-0.3*randn(1,N);%原始信号和噪声叠加
subplot(2,1,2)
plot(t,x)
title('原始信号和噪声叠加')

D=1;
k=5;
y=zeros(1,N);
y(1:D+k)=x(1:D+k);
W=1/k*ones(1,k);
e=0;
u=0.001;%收敛因子u，可以选为0.001,或者0.3
%u=0.3
for i=(D+k+1):N
    X=x((i-D-k+1):(i-D));
    y(i)=W*X';
    e=x(i)-y(i);
    W=W+2*u*e*X;
end
figure(2);
plot(t,y)
title('采用LMS后的信号')