【正文】 頻率響應(yīng)曲線是否滿足指標(biāo)。（2）課程設(shè)計(jì)附件（主要是源程序）。由分析結(jié)果可知，濾波器后的語音信號與原始信號基本一致，即設(shè)計(jì)的并聯(lián)型IIR切比雪夫I型濾波器能夠去除信號中所加單頻噪聲，因此達(dá)到了課程設(shè)計(jì)目的。本課程設(shè)計(jì)的主要目的是通過設(shè)計(jì)一個濾波器對含噪的語音信號進(jìn)行濾波去噪處理，并在處理時采用雙線性變換法設(shè)計(jì)的并聯(lián)型切比雪夫I型濾波器。這一點(diǎn)與驗(yàn)證性的基本實(shí)驗(yàn)有本質(zhì)性的區(qū)別。（3）處理結(jié)果和分析結(jié)論應(yīng)該一致，而且應(yīng)符合理論。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如C、Fortran）的編輯模式，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先進(jìn)水平[1]。MATLAB已經(jīng)成為國際上最流行的科學(xué)與工程計(jì)算軟件工具，現(xiàn)在的MATLAB已經(jīng)不僅僅是一個“矩陣實(shí)驗(yàn)室”了，它已經(jīng)成為一種具有廣泛應(yīng)用前景的全新的計(jì)算機(jī)高級語言，在國內(nèi)外高校和研究部門中正扮演重要角色[1]。同一種系統(tǒng)函數(shù)H（z）可以有多種不同的結(jié)構(gòu)，基本網(wǎng)絡(luò)結(jié)構(gòu)有直接Ⅰ型、直接Ⅱ型、級聯(lián)型、并聯(lián)型四種，都具有反饋回路。這一套成熟、有效的設(shè)計(jì)方法，也可通過各種變換所得到。通過選擇一種具有等波紋特性而不是單調(diào)特性的逼近方法可以實(shí)現(xiàn)這一點(diǎn)。其幅度平方函數(shù)表示為： （21）上式中，ε為小于1的正數(shù)，表示通帶內(nèi)幅度波動的程度，ε愈大，波動幅度也愈大。 （22）是通帶邊界頻率，一般是預(yù)先給定的。N的值是根據(jù)阻帶的邊界條件來確定的，當(dāng)時，,即 （27）從而 （28）將此條件代入=cosh(N arcosh x),得到 （29）因此，要求阻帶邊界頻率處的衰減越大，也就是過渡帶內(nèi)幅頻特性越陡，則需用的階數(shù)N越高。通?？梢越o定的參數(shù)為：低通數(shù)字濾波器通帶邊界頻率、阻帶邊界頻率和對應(yīng)的通帶衰減函數(shù)、阻帶衰減函數(shù)。一般情況下，可以通過頻率的預(yù)畸變進(jìn)行校正。 錄制語音信號先在windows平臺下用麥克風(fēng)錄制一段語音，錄音機(jī)的路徑：開始—所有程序—附件—娛樂工具—錄音機(jī)，錄制時，可以通過點(diǎn)擊錄音機(jī)界面上的文件—屬性，在彈出的對話框中選擇馬上轉(zhuǎn)換，設(shè)置錄音格式為PCM編碼，屬性為8kHz，單聲道，7kB/s，這樣的設(shè)置在相同的時間下數(shù)據(jù)最少，處理時間也就最短。如“”。)。 % 計(jì)算信號x的長度fn=2200。 sound(y,fs,bits)。 Y=Y(1:N/2)。圖形如35所示：。fp=2000。 % 定義通帶波紋和阻帶衰減wp=fp/fs*2*pi。Omegas=2/T*tan(ws/2)。再使用plot等函數(shù)畫出濾波器的幅度，相位圖，然后觀察圖36，思考所設(shè)計(jì)的率波器是否達(dá)標(biāo)，能否較好的完成濾波功能，不足是在哪里，濾波器的參數(shù)是否適合，是否需要進(jìn)行改進(jìn)等，這些都需要我們從整幅圖中得出的信息。Y_fil=Y_fil(1:N/2)。通過以下語句來進(jìn)行語音信號回放比較：sound（x，fs，bits） %播放原始單音頻信號sound（y，fs，bits）wavwrite(y,fs,bits,39。39。來自這些節(jié)的輸出相加以形成濾波器的輸出。4出現(xiàn)的問題及解決方法在這次的課程設(shè)計(jì)中，由于理論知識的不踏實(shí)以及其他各種原因，我們遇到了不少問題。（5）在用wavread函數(shù)加入語音信號時，提示錯誤。（7）不明白原理，無從下手處理問題，此時就需要將學(xué)過的理論知識進(jìn)行溫習(xí)一遍，有不會的就像老師請教。在二周的課程設(shè)計(jì)時間里，我每天都過的很充實(shí)，查資料、讀程序，反復(fù)琢磨，學(xué)習(xí)的過程是艱辛的，但是同時也是快樂的。在整個課程設(shè)計(jì)過程中我再次認(rèn)識到堅(jiān)持、耐心、細(xì)心等品質(zhì)的重要，這對今后的學(xué)習(xí)和工作是有很大幫助的。這些在書本上是沒有的，如果書本基礎(chǔ)那么實(shí)踐就是提高，沒有基礎(chǔ)當(dāng)然不可能會有提高，但是沒有實(shí)際理論也不會有很多的價(jià)值。我在其中體會很深，開始的階段是會碰到很多的問題，對軟件的不熟練，對模塊的不認(rèn)識，對原理的不理解...這些使我感到了壓力和沮喪，但是困難只要去解決總是有辦法的，在老師和同學(xué)的幫助下，自己的努力中把碰到的問題一一化解，心中的愉快是不言而喻。但是在此要感謝指導(dǎo)過我的老師特別是胡雙紅老師，給予了我很大幫助，同時給了我深深的認(rèn)識和體會。)。 % 計(jì)算信號x的長度fn=2200。 sound(y,fs,bits)。)。ylabel(39。plot(t,y) % 繪制截取后的語音信號title(39。時間(單位：s)39。)。 Y=Y(1:N/2)。plot(f,abs(X)) % 繪制頻譜圖title(39。頻率（單位：Hz）39。) 。)。ylabel(39。fc=2150。ws=fc/fs*2*pi。 %截止頻率預(yù)畸變[cs,ds]=afd_chb1(Omegap,Omegas,Rp,As) % 計(jì)算濾波器系統(tǒng)函數(shù)分子分母系數(shù)[b,a]=bilinear(cs,ds,1/T) % 雙線性變換得到數(shù)字濾波器系統(tǒng)函數(shù)分子分母系數(shù)[db,mag,pha,grd,w]=freqz_m(b,a)。plot(w/pi,db,39。grid on。As,Rp,10,10]。r39。LineStyle39。濾波器幅度響應(yīng)圖39。)。subplot(2,2,2)。axis([0,1,0,])。)。ylabel(39。plot(w/pi,pha)。濾波器相位響應(yīng)圖39。)。% 繪制相位特性曲線圖subplot(2,2,4)。 axis([0,100,])。)。ylabel(39。 % 對信號y進(jìn)行濾波處理，其中b,a為上面設(shè)計(jì)好的濾波器參數(shù)Y_fil=abs(fft(y_fil))。plot(t,x)。原始語音信號時間x39。)。% 繪制原始錄制的語音信號時間波形圖subplot(3,2,3)。title(3