clear?all;
%-------------------------------讀入帶噪語音文件---------------------------
[filenamepathname]=uigetfile（‘*.wav‘‘請選擇純凈語音文件：‘）;
tidy=wavread（[pathname?filename]）‘;
[filenamepathname]=uigetfile（‘*.wav‘‘請選擇帶噪語音文件：‘）;
wavin=wavread（[pathname?filename]）‘;
%-------------------------------參數定義---------------------------------
frame_len=256;?%幀長
step_len=0.5*frame_len;?%分幀時的步長，相當于重疊50%
wav_length=length（wavin）;
R?=?step_len;
L?=?frame_len;?
f?=?（wav_length-mod（wav_lengthframe_len））/frame_len;
k?=?2*f-1;?%?幀數
h?=?sqrt（1/101.3434）*hamming（256）‘;?%?漢寧窗乘以系數的原因是使其復合條件要求；
wavin?=?wavin（1:f*L）;??%?帶噪語音與純凈語音長度對齊
tidy=?tidy（1:f*L）;
win?=?zeros（1f*L）;?%?設定初始值；
enspeech?=?zeros（1f*L）;?????????????????????????
%-------------------------------分幀-------------------------------------
for?r?=?1:k?
????y?=?wavin（1+（r-1）*R:L+（r-1）*R）;?%?對帶噪語音幀間重疊一半取值；
????y?=?y.*h;?%?對取得的每一幀都加窗處理；
????w?=?fft（y）;?%?對每一幀都作傅里葉變換；
????Y（1+（r-1）*L:r*L）?=?w（1:L）;?%?把傅里葉變換值放在Y中；
end
%-------------------------------估計噪聲-----------------------------------
???NOISE=?stationary_noise_evaluate（YLk）;?%噪聲最小值跟蹤算法
%?????NOISE=?non_stationary_noise_evaluate（YLk）;?%?基于統計信息的非平穩噪聲自適應算法
%-------------------------------譜減法-------------------------------------
for?????t?=?1:k?????
?????????X?=?abs（Y）.^2;???
?????????S?=?X（1+（t-1）*L:t*L）-NOISE（1+（t-1）*L:t*L）;?%?含噪語音功率譜減去噪聲功率譜；
?????????S?=?sqrt（S）;
?????????A?=?Y（1+（t-1）*L:t*L）./abs（Y（1+（t-1）*L:t*L））;?%?帶噪于語音的相位；
?????????S?=?S.*A;?%?因為人耳對相位的感覺不明顯，所以恢復時用的是帶噪語音的相位信息；
?????????s?=?ifft（S）;???
?????????s?=?real（s）;?%?取實部；
?????????enspeech（1+（t-1）*L/2:L+（t-1）*L/2）?=?enspeech（1+（t-1）*L/2:L+（t-1）*L/2）+s;?%?在實域疊接相加；
?????????win（1+（t-1）*L/2:L+（t-1）*L/2）?=?win（1+（t-1）*L/2:L+（t-1）*L/2）+h;?%?窗的疊接相加；
end
enspeech?=?enspeech./win;?%?去除加窗引起的增益得到增強的語音；
%-----------------------信噪比--------------------------------------------
SNR_before=SNR1（tidywavin）;
SNR_after=SNR2（tidyenspeech）;
%-------------------------------畫出波形----------------------------------------

?subplot（311）;plot（tidy）;title（‘（a）純凈語音‘）;xlabel（‘樣點數‘）;ylabel（‘幅度‘）;axis（[0?2.5*10^4?-0.3?0.3]）;
?subplot（312）;plot（wavin）;title（‘（b）帶噪語音（10dB白噪聲）‘）;xlabel（‘樣點數‘）;ylabel（‘幅度‘）;axis（[0?2.5*10^4?-0.3?0.3]）;
?subplot（313）;plot（enspeech）;title（‘（c）譜減法-增強語音‘）;xlabel（‘樣點數‘）;ylabel（‘幅度‘）;axis（[0?2.5*10^4?-0.3?0.3]）;
?axis（[0?7*10^4?-1?1]）;
?axis（[0?2.5*10^4?-0.3?0.3]）;
?wavwrite（enspeech‘spectruesub_enspeech.wav‘）;?%?寫出增強語音；

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-06-08?17:25??20160608語音教材MATLAB程序提交出版社\
?????目錄???????????0??2016-06-07?20:10??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\
?????文件???????19968??2016-06-08?17:53??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\Readme.doc
?????目錄???????????0??2016-06-07?20:08??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\
?????文件???????45102??2009-06-05?17:00??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\SNR0_sp01.wav
?????文件???????45102??2009-06-05?17:00??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\SNR10_sp01.wav
?????文件???????45102??2009-06-05?17:08??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\SNR15_sp01.wav
?????文件???????45102??2009-06-05?18:02??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\SNR5_sp01.wav
?????文件????????2720??2016-03-10?14:46??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\chapter10_1.m
?????文件????????2723??2016-03-07?11:59??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\chapter10_2.m
?????文件????????3899??2010-03-31?19:19??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\chapter10_3.m
?????文件????????4678??2015-12-19?22:41??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\chapter10_4.m
?????文件????????2023??2016-03-10?14:55??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\chapter10_5.m
?????文件???????45102??2005-07-12?08:58??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\sp01.wav
?????文件?????????987??2009-05-12?15:28??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第10章\stationary_noise_evaluate.m
?????目錄???????????0??2016-06-07?20:08??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\
?????文件?????????349??2016-03-08?15:19??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_1.m
?????文件?????????366??2016-03-08?15:20??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_2.m
?????文件?????????274??2016-03-08?15:21??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_3.m
?????文件?????????373??2016-03-08?15:22??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_4.m
?????文件?????????318??2016-03-08?15:23??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_5.m
?????文件?????????307??2016-03-08?15:24??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_6.m
?????文件?????????248??2016-03-08?15:25??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_7.m
?????文件????????2597??2016-03-08?15:25??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_8.m
?????文件????????1570??2016-03-09?11:49??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\chapter11_9.m
?????文件????????3562??2015-12-20?19:58??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\wu.wav
?????文件????????1650??2016-03-08?15:47??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第11章\yingshe.m
?????目錄???????????0??2016-06-07?20:08??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第12章\
?????目錄???????????0??2016-06-07?20:08??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第12章\RBF\
?????文件?????????967??2016-02-24?20:35??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第12章\RBF\RBF.asv
?????文件????????1445??2016-02-24?20:37??20160608語音教材MATLAB程序提交出版社\語音教材MATLAB程序提交出版社\第12章\RBF\RBF_Test.m
............此處省略607個文件信息