%對語言信號進行處理，加噪和濾波
[yfsnbits]=wavread（‘voice.wav‘）;??????%讀取語音信號
y=80*y;?????????????????????????????????%若聲音太小，在此處調節放大倍數
sound（yfsnbits）;??????????????????????%回放錄制的語音信號
pause（12）;??????????????????????????????%延時程序

%對音頻信號加入噪聲
?L=length（y）;??????????????????????????%計算音頻信號的長度
?noise=0.08*randn（L1）;????????????????%加入高斯白噪聲
?y_z=y+noise;??????????????????????????%將兩個信號疊加成一個新的信號——加噪聲處理???
?sound（y_zfsnbits）;??????????????????%回放加噪后的語音信號
?pause（12）;????????????????????????????%延時程序

%對加噪的語音信號進行低通濾波處理程序如下：?
%低通濾波器設計階段
fp=1000;fc=1100;As=100;Ap=1;?????????????%低通濾波器的性能指標，fp為通帶截至頻率，fc為阻帶截至頻率，fc-fp為過渡帶的帶寬，AS為帶外抑制度，AP為通帶紋波
wc=2*fc/fs;?wp=2*fp/fs;????????????????%對頻率進行歸一化處理
[Nwn]=?buttord（wpwcApAs）;??????????%設計一個巴特沃斯濾波器，得到最低階數和截止頻率n為濾波器的最低階數，wn為3dB的截至頻率
b=fir1（Nwn‘low