LMS自適應濾波

clear all;
t=0:0.0001:1.9999;
N=20000;

%生成有用訊號s
s=3*sin(2*pi*100*t);
figure(1);
subplot(3,1,1);
plot(t,real(s));
title('有用訊號s');
xlabel('n');
ylabel('s');

%生成兩個相關的噪聲n1,n2
n1=wgn(1,N,3);
n2=zeros(1,N);
n2(1:18000)=n1(1:18000);
n3=wgn(1,N,3);
n2(18001:N)=n3(18001:N);
figure(1);
subplot(3,1,2);
plot(t,real(n2));
title('噪聲參考輸入');
xlabel('n');
ylabel('n2');

%輸入端期望訊號d=s+n1
d=s+n1;
subplot(3,1,3);
plot(t,real(d));
title('含有噪聲的訊號d');
xlabel('n');
ylabel('d');

%初始化各參量
M=8;  %濾波器階數
w=zeros(M,1); %濾波器權係數
mu=0.01;  %更新步長
iter=2000; %迭代次數
x=zeros(M,1); %輸入向量

%LMS迴圈演算法
for i=1:iter
    i
   x(1)=n2(i);   
   y(i)=dot(w,x);
   e=d(i)-y(i);
   w=w+mu*e*x;
   
    for j=M:-1:2
        x(j)=x(j-1);    
    end

    error(i)=e; %儲存輸出誤差
   cost(i)=e*e;  
end

   

figure(2);
subplot(1,1,1);
plot(error);
title('除噪後的訊號');
xlabel('n');
ylabel('error');