denoising.m

来自「噪声环境下的予以增强程序」· M 代码 · 共 27 行

clear
clc
%在噪声环境下语音信号的增强
%语音信号为读入的声音文件
%噪声为正态随机噪声
sound=wavread('sound.wav');
count=length(sound);
noise=0.1*randn(1,count);
y=sound'+noise;
%用小波函数'db6'对信号进行3层分解
[c,l]=wavedec(y,3,'db6');
%估计尺度l的噪声标准偏差
sigma=wnoisest(c,l,1);
alpha=2;
%获取消噪过程中的阈值
thr=wbmpen(c,l,sigma,alpha);
keepapp=1;
%对信号进行消噪
yd=wdencmp('gbl',c,l,'db6',3,thr,'s',keepapp);
subplot(121);
plot(sound);
title('原始语音信号');
subplot(122);
plot(yd);
title('去噪后的语音信号');
wavwrite(yd,11025,8,'sounddenoising.wav');
wavwrite(y,11025,8,'soundnoising.wav');