【正文】 Hz39。)。subplot(313)。y1=fft(x1)。x1=x+noise。 %繪制原始語(yǔ)音信號(hào)頻率響應(yīng)圖title(39。y1=fft(x1,1024)。subplot(211)。xlabel(39。title(39。參考文獻(xiàn)[1]王風(fēng)文，舒冬梅，趙宏才．?dāng)?shù)字信號(hào)處理（第2版） [M]．北京郵電大學(xué)出版社2006:3336,192193[2]劉波，文忠， [M].電子工業(yè)出版社2006:2933[3]高西全，丁玉美，——原理實(shí)現(xiàn)及應(yīng)用(第二版) [M].電子工業(yè)出版社2010:178192[4] [M].電子工業(yè)出版社 2005 年7月:153162[5]丁玉美,(第二版) [M].西安電子科技大學(xué)出版社2000:175193附錄源程序代碼：%語(yǔ)音信號(hào)采集程序設(shè)計(jì)fs=22050。圖33原始信號(hào)和加噪信號(hào)時(shí)域圖圖34原始信號(hào)和加噪后信號(hào)的頻譜圖通過(guò)對(duì)原始語(yǔ)音信號(hào)的頻譜圖與加噪后的語(yǔ)音信號(hào)頻譜圖的對(duì)比可得，加噪后語(yǔ)音信號(hào)頻率幅值發(fā)生了明顯的變化。常見(jiàn)的窗函數(shù)，可以分為以下主要類型：（1）冪窗采用時(shí)間變量某種冪次的函數(shù)，如矩形、三角形、梯形或其它時(shí)間t的高次冪；（2）三角函數(shù)窗應(yīng)用三角函數(shù)，即正弦或余弦函數(shù)等組合成復(fù)合函數(shù)，例如漢寧窗、海明窗、Kaiser窗等；（3）指數(shù)窗采用指數(shù)時(shí)間函數(shù)，例如高斯窗等。無(wú)限長(zhǎng)單位沖激響應(yīng)（IIR）濾波器有以下幾個(gè)特點(diǎn):（1）系統(tǒng)的單位沖激響應(yīng)h(n)是無(wú)限長(zhǎng)；（2）系統(tǒng)函數(shù)H(z)在有限z平面（0|z|）；（3）結(jié)構(gòu)上存在著輸出到輸入的反饋，也就是結(jié)構(gòu)上是遞歸型的。采樣頻率越高，即采樣的間隔時(shí)間越短，則在單位時(shí)間內(nèi)計(jì)算機(jī)得到的聲音樣本數(shù)據(jù)就越多，對(duì)聲音波形的表示也越精確。 1924年奈奎斯特(Nyquist)就推導(dǎo)出在理想低通信道的最高大碼元傳輸速率的公式: 理想低通信道的最高大碼元傳輸速率=2W*log2 N (其中W是理想低通信道的帶寬,N是電平強(qiáng)度)為什么把采樣頻率設(shè)為8kHz?在數(shù)字通信中，根據(jù)采樣定理, 最小采樣頻率為語(yǔ)音信號(hào)最高頻率的2倍頻帶為F的連續(xù)信號(hào) f(t)可用一系列離散的采樣值f(t1),f(t1177。Δt)，f(t1177。采樣頻率與聲音頻率之間有一定的關(guān)系，根據(jù)采樣定理，只有采樣頻率高于聲音信號(hào)最高頻率的兩倍時(shí)，才能把數(shù)字信號(hào)表示的聲音還原成為原來(lái)的聲音。IIR濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。其性能如表11所示：表11 常見(jiàn)窗函數(shù)性能表名稱濾波器過(guò)渡帶寬最小阻帶衰減名稱濾波器過(guò)渡帶寬最小阻帶衰減矩形21dBPARZENWIN56dB巴特利特25dBFLATTOPWIN108dB漢寧44dBGAUSSWIN60dB漢明51dBBARTHANNWIN40dB布萊克曼11π/M74dBBLACKMANHARRIS109dB本設(shè)計(jì)采用的是凱塞窗（Kaiser）設(shè)計(jì)的濾波器，凱塞窗（Kaiser）通過(guò)改變參數(shù)可以達(dá)到不同的性能，正由于其優(yōu)良的窗函數(shù)特性，廣泛應(yīng)用于高通、低通、帶通、帶阻等各種濾波器的設(shè)計(jì)。再通過(guò)對(duì)原始語(yǔ)音信號(hào)的回放效果與加噪后的語(yǔ)音信號(hào)回放的效果的對(duì)比，人耳可以明顯辨別出兩種語(yǔ)音信號(hào)不一樣了，加噪后的語(yǔ)音信號(hào)在聽(tīng)覺(jué)上比原始語(yǔ)音信號(hào)要渾濁很多，而且還有吱吱嘎嘎的混雜音。 %語(yǔ)音信號(hào)采集頻率為22050[x,fs,bits]=wavread(39。原始信號(hào)相位39。Hz39。plot(x)。 %對(duì)加噪信號(hào)做1024點(diǎn)FFT變換f=fs*(0:511)/512。加噪語(yǔ)音信號(hào)頻譜圖39。 %利用kaiser濾波器對(duì)語(yǔ)音信號(hào)濾波168。x2=fftfilt(hn,x1)。 plot(x2)。subplot(312)。)。)xlabel(39。fuzhi39。)。y=fft(x)。noise=*randn(L,2)。plot(f,abs(y1(1:512)))。figure(2)。 %對(duì)信號(hào)做1024點(diǎn)FFT變換figure(1)。)。plot(angX)。然后與不同信源相同濾波方法的同學(xué)比較各種信源的特點(diǎn)，與相同信源不同濾波方法的同學(xué)比較各種濾波方法性能的優(yōu)劣。原始語(yǔ)音信號(hào)的幅值和相位如圖31所示：圖31 原始信號(hào)幅值和相位原始信號(hào)時(shí)域和頻域圖如圖32所示：圖32 原始信