自适应滤波器MATLAB仿真(3)

clear

t=0:1/10000:1-0.0001;%设置迭代次数 s=sin(2*pi*t);%标准正弦信号

n=randn(size(t));%与时间t等长随机信号 x=s+n;%加噪信号

w=[0,0.5];%初始2阶加权系数 u=0.00026;%最佳参数

for i=1:9999;%自适应算法 y(i+1)=n(i:i+1)*w'; e(i+1)=x(i+1)-y(i+1); w=w+2*u*e(i+1)*n(i:i+1); end;

%画图程序 figure(1)

subplot(4,1,1) plot(t,n);

title('Noise signal');xlabel('t');ylabel('n(t)'); subplot(4,1,2) plot(t,s);

title('Sinusoidal signal');xlabel('t');ylabel('s(t)'); subplot(4,1,3) plot(t,x);

title('Sinusoidal signal with noise');xlabel('t');ylabel('x(t)'); subplot(4,1,4) plot(t,e);

title('the result of filtering');xlabel('t');ylabel('e(t)');