谱减法降噪的matlab代码实现

最新推荐文章于 2024-08-02 22:15:09 发布

phymat.nico

最新推荐文章于 2024-08-02 22:15:09 发布

阅读量7.8k

点赞数 4

分类专栏：信号处理

本文链接：https://blog.csdn.net/shixin_0125/article/details/100518511

版权

本文详细介绍了如何使用MATLAB进行谱减法降噪的代码实现，包括噪声估计、频谱减法原理及具体步骤，适用于信号处理和噪声抑制领域的研究。

摘要由CSDN通过智能技术生成

clc;clear;close all;
winsize=256;%窗长
n=0.2;%噪声水平
a=2;
b=5;
c=0.01;
fs=44100;
[speech,fs,nbits]=wavread('file_1.wav') ;%读入wav文件
sound(speech,fs) ;
pause(5);
size=length(speech);%语音长度
numofwin=floor(size/winsize);%窗数
%定义汉明窗
ham=hamming(winsize)';
hamwin=zeros(1,size);
enhanced=zeros(1,size);
improved=zeros(1,size);
%生成噪声信号

y=speech(:,1)'+n*randn(1,size)
sound(y,fs) ;
pause(5);
%噪声处理
noisy=n*randn(1,winsize);
Q=fft(speech);
N=fft(noisy);
npow=abs(N);
for q=1:2*numofwin-1
yframe=y(1+(q-1)*winsize/2:winsize+(q-1)*winsize/2);%分帧
hamwin(1+(q-1)*winsize/2:winsize+(q-1)*winsize/2)=hamwin(1+(q-1)*winsize/2:winsize+(q-1)*winsize/2)+ham;%
%加噪信号FFT
y1=fft(yframe.*ham);
ypow=abs(y1);%加噪信号幅度
yangle=angle(y1);%相位
%计算功率谱密度
Py=ypow.^2;
Pn=npow.^2;
Pyy=ypow.^a;
Pnn=npow.^a;
%基本谱减
for i=1:winsize
if Py(i)-Pn(i)>0
Ps(i