jiaochalms.m

来自「传统双话筒降噪应用中」· M 代码 · 共 36 行

sample_n=60000;
[n1,n2,n3,n4,n5]=dataread5('D:\chen\科研\matlab算法汇总\2006年MATLAB语音增强算法\5ch8k20cm02.dat',sample_n);%噪声；
if max(n1)>max(n2)
    dd=max(n1);
else
    dd=max(n2);
end
x=n1/dd;y=n2/dd;


Eout=zeros(1,sample_n);
yout=zeros(1,sample_n);
N=32;
%---------LMS--------
W=zeros(1,N);
u=0.02;           %0.00000001;
for i=N+1:sample_n
   yout(1,i)=W(1,1:N)*x(1,i-N+1:i)';
   Eout(1,i)=y(1,i)-yout(1,i);
   if i<10000
   W(1,1:N)=W(1,1:N)+2*u*Eout(1,i)*x(1,i-N+1:i);
   end
end
W(1,1:N)=zeros(1,N);
for i=N+1:sample_n
   yout2(1,i)=W(1,1:N)*Eout(1,i-N+1:i)';
   Eout2(1,i)=y(1,i)-yout2(1,i);
   W(1,1:N)=W(1,1:N)+2*u*Eout2(1,i)*Eout(1,i-N+1:i);  
end

E=Eout2;
plot(x(N+1:sample_n))
hold on
plot(yout2(N+1:sample_n),'r')