【正文】 序流程圖 音頻信號部分程序 （ 1）寫入聲音信號 [y,Fs,bits]=wavread(39。198。 最后 謝謝陳華老師上課的細致講解 ，使我這次課設得以順利地完成。)。xlabel(39。as1=65。%求出模擬濾波器階數(shù) N 和截止頻率 [B,A]=butter(N,wc,39。幅度 39。Y1=fft(y1,n)。)。grid on。xlabel(39。 subplot(3,2,5)。as=70。 figure(2)。Hz39。xlabel(39。 subplot(3,2,6)。 sound(y1,Fs,bits)。title(39。Y1=Y1(1:n/2,:)。)。s39。 %高通濾波器設計 Fs=16000。title(39。 Y=Y(1:n/2,:)。)。 subplot(3,2,4)。)。 xlabel(39。wsh=[wsh1,wsh2]。fp2=4000。)。 subplot(3,2,1)?？傊估碚撀?lián)系了實際，鞏固并深化了對課本基本知識的認識和理解，使理論得以升華。11\181。 [Bz,Az]=bilinear(B,A,) 程序流程圖 首先確定性能指標，把頻率轉化為數(shù)字角頻率，進而在進行頻率預畸變，用Ω =tan(w/2)對高 （帶） 通模擬數(shù)字濾波器 H(z)的數(shù)字邊界頻率預畸變 ,得到高（帶） 通模擬濾波器 H(s)的邊界頻率主要是通帶截止頻率 Wp。fs2=5000。 ws=2*pi*fs/Fs (3) 預畸變校正求出相應模擬濾波器模擬角頻率指標 用Ω =tan(w/2)對模擬濾波器 濾波器 的數(shù)字邊界頻率預畸變 ,得到模擬濾波器的邊界頻率主要是通帶截止頻率 Wp。S 平面整個 jΩ 軸單值地對應于 Z 平面單位圓一周，即頻率軸是單值變換關系。 IIR 數(shù)字濾波器的特征是，具有無限持續(xù)時間沖激響應，需要用遞歸模型 來實現(xiàn)，其差分方程為： ? ? ? ? ? ????? ????Ni iNi i inyinxany 10 b 系統(tǒng)函數(shù)為： ? ????????? NkkkMrrrzazbzH101 設計 IIR濾波器的任務就是尋求一個物理上可實現(xiàn)的系統(tǒng)函數(shù) H(z)，使其頻率響應 H(z)滿足所希望得到的頻域指標，即符合給定的通帶截止頻率、阻帶截止頻率、通帶衰減系數(shù)和阻帶衰減系數(shù)。 IIR 數(shù)字濾波器設計原理 IIR數(shù)字濾波器是一種離散時間系統(tǒng)，其系統(tǒng)函數(shù)為 ? ? ? ?? ?zX zYzazNkkkkk?????????1M0k1bzH 假設 M≤ N，當 M＞ N時 ,系統(tǒng)函數(shù)可以看作一個 IIR的子系統(tǒng)和一個 (MN)的 FIR子系統(tǒng)的級聯(lián)。這個關系重寫如下： 上式表明， S 平面上 Ω 與 Z 平面的 ω 成非線性的正切關系，如圖 14所示。阻帶截止頻率 Ws 的轉換。Fs=16000。阻帶截止頻率 Ws的轉換。247。 在做本次課程設計的過程中，我深深地感受到了自己所學到知識的有限，明白了只學好課本上的知識是不夠的，要通過圖書館和互聯(lián)網(wǎng)等各種渠道來擴充自己的知識。plot(t,y)。grid on。fs2=5000。%帶通指標 [N,wc]=buttord(wph,wsh,ap,as,39。Hz39。grid on。plot(w*Fs/n,abs(Y1))。title(39。%取一半頻譜 figure(1)。濾波前音樂信號頻譜 39。fp=6000。)。 subplot(2,1,2)。 figure(1)。濾波后音樂信號頻譜 39。 figure(3)。grid on。plot(t,y1)。xlabel(39。%雙線性法把模擬高通轉化為數(shù)字高通 [H,w]=freqz(Bz,Az)。rp=1。grid on。plot(t,y)。)。Hz39。 y1=filter(Bz,Az,y)。ylabel(39。)。ap=1。plot(k*Fs/n,angle(Y))。秒 39。但是從中我們學習到了如何對待遇到的困難，進一步 培養(yǎng)了我們一絲不茍的科學態(tài)度和 不厭其煩的耐心。195。 上述準備工作做好之后，最后利用雙線性變換法將模擬高 （帶） 通濾波器 Ha(s)轉換成數(shù)字高 （帶） 通濾波器 H(z)。as1=65。 通帶截止頻率 wph=tan(wp/2)。 11112????? zzTs sTsTsTsTz?????? 222121???????????? ?? jTjeeTs jj 2tan2112 ??????????jTjTz??22222222||???????? ????????? ????TTz?－ ? ?o? ?2ta n2 ?T? ＝???????? 2tan2 ?T 5 圖