%例程7－7  自適應信號分離器
%自適應信號分離器

t=0:1/10:400;
s=sin(2*pi*t/10);                 %周期信號
x=awgn(s,15);                      
D=50;%延遲
r=[zeros(1,D),x];                 %信號延遲D
x=[x zeros(1,D)];

N=5;                          %r經LMS自適應濾波
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([1 200 -1.2 1.2]);
subplot(3,1,2); 
plot(t,y(1:(length(x)-D)));
title('周期信號');
axis([1 200 -1.2 1.2]);
subplot(3,1,3);
plot(t,e(1:(length(x)-D)));
title('寬帶信號');
axis([0 200 -1.2 1.2]);