7-8.m

来自「经典数字信号处理滤波器的源代码 重点是利用巴特沃斯模拟滤波器转而设计其它数字滤」· M 代码 · 共 35 行

%例程7－8  自适应窄带干扰抑制

t=0:1/10:400;
s=sin(2*pi*t/10);             %窄带干扰
x=awgn(s,15);               %叠加宽带噪声
D=100;                     %延迟
r=[zeros(1,D),x];
x=[x zeros(1,D)];

%r经LMS自适应滤波
N=5;
u=0.02;
M=length(r);
y=zeros(1,M);
w=zeros(1,N);
for  n=N:M
    x1=r(n:-1:n-N+1);
    y(n)=w*x1';
    e(n)=x(n)-y(n);
    w=w+u.*e(n).*x1;
end

subplot(3,1,1);
plot(t,x(1:(length(x)-D)));
title('输入信号');
axis([50 200 -1.2 1.2]);
subplot(3,1,2); 
plot(t,e(1:(length(x)-D)));
title('预期信号');
axis([50 200 -1.2 1.2]);
subplot(3,1,3);
plot(t,s(1:(length(x)-D)));
title('窄带干扰');
axis([50 200 -1.2 1.2]);