基于自适应语音增强的算法研究与MATLAB仿真(20)

只适应滤波算法

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

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); %算法收敛曲线

t=1:N-k;

plot(t,bi,'r');

title('收敛曲线');

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