新疆大学科学技术学院毕业论文（设计）

plot(s); %信号s时域波形 title('信号s时域波形'); xlabel('n');

axis([0,N,-a-1,a+1]);

xn=awgn(s,5); %加入均值为零的高斯白噪声 %设置初值

y=zeros(1,N); %输出信号y

y(1:k)=xn(1:k); %将输入信号xn的前k个值作为输出y的前k个值 w=zeros(1,k); %设置抽头加权初值 e=zeros(1,N); %误差信号 %用LMS算法迭代滤波 for i=(k+1):N

XN=xn((i-k+1):(i)); y(i)=w*XN'; e(i)=s(i)-y(i); w=w+u*e(i)*XN;

end

pp(q,:)=(e((k+1):N)).^2; end subplot(312)

plot(xn); %信号s时域波形 title('信号s加噪声后的时域波形'); subplot(313)

plot(y); %信号s时域波形 title('自适应滤波后的输出时域波形'); for b=1:N-k

bi(b)=sum(pp(:,b))/g; %求误差统计平均 end

figure(2); %算法收敛曲线

13

新疆大学科学技术学院毕业论文（设计）

t=1:N-k; plot(t,bi,'r'); title('收敛曲线');

hold on %将每次循环的图形显示结果保存下来

信号s时域波形20-20100200300500600700n信号s加噪声后的时域波形400800900100050-501002003004005006007008009001000自适应滤波后的输出时域波形20-201002003004005006007008009001000

图4.1 μ=0.00026系统输入、噪声、输出图

收敛曲线10.5001002003004005006007008009001000图4.1 μ=0.00026收敛曲线

14

新疆大学科学技术学院毕业论文（设计）

信号s时域波形20-202004006008001000n信号s加噪声后的时域波形50-502004006008001000自适应滤波后的输出时域波形50-502004006008001000

图4.2 μ=0.00058系统输入、噪声、输出图

收敛曲线10.5002004006008001000

图4.2 μ=0.00058收敛曲线

15

新疆大学科学技术学院毕业论文（设计）

信号s时域波形20-202004006008001000n信号s加噪声后的时域波形50-502004006008001000自适应滤波后的输出时域波形50-502004006008001000

图4.3μ=0.00116系统输入、噪声、输出图

收敛曲线10.5002004006008001000图4.3 μ=0.00116收敛曲线

16