clear?all;
%-------------------------------讀入帶噪語音文件---------------------------
load?yi.mat
tidy=load（‘yi.mat‘‘rt‘）?;
[filenamepathname]=uigetfile（‘*.cnt‘‘請選擇帶噪語音文件：‘）;
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）;?%?基于統計信息的非平穩噪聲自適應算法
%-------------------------------譜減法-------