%例程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]);