急求懂LMS算法的MATLAB高手帮忙!!!万分感激

来源：百度知道 编辑：UC知道 时间：2024/06/03 20:01:12
急求懂LMS算法的MATLAB高手帮忙!!!万分感激
最近正在做毕设,题目为 基于LMS算法的语音增强
请哪位高手帮忙给个改进LMS算法用于语音增强的程序吧...
真的很急,万分万分的感激,太谢谢啦

clear all；
g=100;
N=256;
k=16;
pp=zeros(g,N-k);
u=0.01;
for q=1:g,
t=1:N;
a=1;
s=a*sin(0.05*pi*t);
figure(1);
subplot(3,1,1)
plot(t,real(s));
title('信号s时域波形');
xlabel('n');
ylabel('s');
axis([0,N,-a-1,a+1]);
xn=awgn(s,5);
y=zeros(1,N);
y(1:k)=xn(1:k);
w=zeros(1,k);
e=zeros(1,N);
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(3,1,2)
plot(t,real(xn));
title('信号s加噪声后的时域波形');
subplot(3,1,3)
plot(t,real(y));
title('自适应滤波后的输出时时域波形');
for c=1:N-k;
bi(c)=sum(pp(:,c))/g;
end;
figure(2);
T=1:N-k;
plot(T,bi,'b');
hold on
plot(T,bi,'b');