jiangzao1.m

这是一个有关语音信号中的噪声降噪的原代码

[y,fs]=wavread('C:\Documents and Settings\feipo\桌面\数字信号实验选作题目\LAB07 语音降噪实验\speech_noisy_8k.wav');
for i=1:25 
     nw=y(1+200*(i-1):200*i);
     NW(i)=var(nw);
end
sum=0;
for i=1:25
    sum=sum+NW(i);
end
sum=sum/25;
for k=1:144
   yi=y(1+200*(k-1):200*k);
   YI=fft(yi);
   Y2=abs(YI).^2;
   SW2=Y2-sum;
  for m=1:length(SW2)
    if SW2(m)<0
      SW2(m)=0;
    end
  end
  SW=SW2.^0.5;
  ae=exp(angle(YI)*j);
  SW1=SW.*ae;
  sw=ifft(SW1);
  ya(1+200*(k-1):200*k)=real(sw);
end
plot(ya);