clear?all
close?all

%語音增強程序
%噪聲為加性高斯白噪聲

%讀入聲音到數(shù)組sound中。
[soundfsnbits]=wavread（‘sound.wav‘）;
sound=sound（:1）;
count=length（sound）;
t=（0:count-1）/fs;

%完成聲音的加噪處理
noise=0.05*randn（1count）;
y=sound‘+noise;
wavwrite（yfs‘compound.wav‘）;

%用小波函數(shù)db6對信號進行三層分解
[cl]=wavedec（y3‘db6‘）;

%估計尺度為1的噪聲標準偏差
sigma=wnoisest（cl