帶阻級(jí)電氣弘深班組題號(hào)設(shè)計(jì)切比雪夫數(shù)字濾波器

成績(jī)91按要求完成了切比雪夫帶阻數(shù)字濾波器設(shè)計(jì)，設(shè)計(jì)步驟完整，從理論推導(dǎo)和仿真結(jié)果看，其性能指標(biāo)不太能滿足設(shè)計(jì)要求，要求的阻帶為，而你們的結(jié)果，阻帶為，與設(shè)計(jì)值有偏差。同時(shí)，結(jié)果的分析太少?！稊?shù)字信號(hào)處理》課程設(shè)計(jì)題目：任務(wù)六：設(shè)計(jì)切比雪夫數(shù)字濾波器班級(jí)學(xué)號(hào)姓名任務(wù)分工工作量(%)11弘深電氣實(shí)驗(yàn)班0120114478王昌金問(wèn)題一分析與總結(jié)、論文撰寫(xiě)12.511弘深電氣實(shí)驗(yàn)班0120113825王宇森問(wèn)題二分析與總結(jié)、論文撰寫(xiě)12.511弘深電氣實(shí)驗(yàn)班0120114336王昱問(wèn)題三分析與總結(jié)、論文撰寫(xiě)12.511弘深電氣實(shí)驗(yàn)班0120114339馬達(dá)問(wèn)題四分析與總結(jié)、論文撰寫(xiě)12.511弘深電氣實(shí)驗(yàn)班0120114395萬(wàn)敏設(shè)計(jì)分析、論文撰寫(xiě)、結(jié)果分析12.511弘深電氣實(shí)驗(yàn)班0120114346曾宇航設(shè)計(jì)分析、仿真操作、結(jié)果分析12.511弘深電氣實(shí)驗(yàn)班0120114482胡昕睿設(shè)計(jì)分析、論文撰寫(xiě)、結(jié)果分析12.511弘深電氣實(shí)驗(yàn)班0120114484劉紅梅設(shè)計(jì)分析、論文綜合、結(jié)果分析12.52013年12月目錄TOC\o"1-3"\h\u296831.濾波器的性能特點(diǎn)及幅度響應(yīng)圖 2261851.1低通濾波器 2323301.2高通濾波器 2290341.3帶通濾波器 335991.4帶阻濾波器 381252.比較模擬巴特沃斯濾波器和切比雪夫?yàn)V波器的頻域特性 458412.1巴特沃斯濾波器 4175212.1.1巴特沃斯濾波器的簡(jiǎn)介 4321092.1.2巴特沃斯濾波器幅度特性 5239552.1.3巴特沃斯低通濾波器階次的確定 6270722.2切比雪夫低通濾波器 655132.2.1切比雪夫低通濾波器簡(jiǎn)介 6242372.2.2切比雪夫低通濾波器幅度特性 795912.3兩種低通模擬濾波器的比較 1049683.濾波器的設(shè)計(jì)方法 11183923.1巴特沃斯（Butterworth）濾波器 1182443.1.1巴特沃斯濾波器的幅度平方響應(yīng) 1194183.1.2設(shè)計(jì)巴特沃斯低通濾波器的指標(biāo) 12250753.1.3巴特沃斯低通濾波器的設(shè)計(jì)方法和步驟 13117523.2切比雪夫（Chebyshev）濾波器 15212953.2.1切比雪夫?yàn)V波器的幅度平方響應(yīng) 15198643.2.2設(shè)計(jì)切比雪夫?yàn)V波器的指標(biāo) 18201113.2.3切比雪夫Ⅰ型濾波器的設(shè)計(jì)方法和步驟 18118424.模擬低通濾波器轉(zhuǎn)換為數(shù)字濾波器的設(shè)計(jì)方法與流程 2135404.1三種設(shè)計(jì)方法 21171314.2設(shè)計(jì)步驟 22211505.切比雪夫模擬低通濾波器的設(shè)計(jì) 23223495.1設(shè)計(jì)技術(shù)指標(biāo)及要求 2338675.2設(shè)計(jì)方法選擇和原理簡(jiǎn)介 23129425.3設(shè)計(jì)流程圖 24210765.4設(shè)計(jì)過(guò)程 24216536.MATLAB濾波器仿真及結(jié)果分析 26269556.1MATLAB自編程設(shè)計(jì)濾波器設(shè)計(jì)步驟 26225736.2MATLAB仿真結(jié)果與分析 27179617.附件 301.濾波器的性能特點(diǎn)及幅度響應(yīng)圖1.1低通濾波器從頻率之間，幅頻特性平直，可以是信號(hào)中低于的成分幾乎不受衰減的通過(guò)，而高于的頻率成分受到極大地衰減。圖1-1低通濾波器1.2高通濾波器與低通濾波器相反，從頻率，其幅頻特性平直。使信號(hào)中高于的頻率成分幾乎不受衰減地通過(guò)，而低于的頻率成分將受到極大地衰減。圖1-2高通濾波器1.3帶通濾波器帶通濾波器的通頻帶在之間，使信號(hào)中高于而低于的頻率成分可以不受衰減地通過(guò)，而其它成分受到衰減。圖1-3帶同濾波器1.4帶阻濾波器與帶阻濾波器相反，阻帶在頻率之間。它使信號(hào)中高于而低于的頻率成分可以衰減，而其它成分的信號(hào)幾乎不受衰減地通過(guò)。圖1-4帶阻濾波器比較模擬巴特沃斯濾波器和切比雪夫?yàn)V波器的頻域特性2.1巴特沃斯濾波器2.1.1巴特沃斯濾波器的簡(jiǎn)介巴特沃斯濾波器是電子濾波器的一種。巴特沃斯濾波器的特點(diǎn)是通頻帶的頻率響應(yīng)曲線最平滑。這種濾波器最先由英國(guó)工程師斯提芬·巴特沃斯在1930年發(fā)表在英國(guó)《無(wú)線電工程》期刊的一篇論文中提出的。巴特沃斯濾波器的特點(diǎn)是通頻帶內(nèi)的頻率響應(yīng)曲線最大限度平坦，沒(méi)有起伏，而在阻頻帶則逐漸下降為零。在振幅的對(duì)數(shù)對(duì)角頻率的波特圖上，從某一邊界頻率開(kāi)始，振幅隨著角頻率的增加而逐漸減少，趨向負(fù)無(wú)窮大。見(jiàn)圖2-2。圖2-1一階巴特沃斯低通濾波器的幅頻特性2.1.2巴特沃斯濾波器幅度特性巴特沃斯低通濾波器可用如下幅度平方函數(shù)表示如下：其中，為濾波器的階數(shù),為截止頻率，即振幅下降為-3分貝時(shí)的頻率，

為通頻帶邊緣頻率，為通頻帶邊緣頻率處的幅度平方值。當(dāng)截止頻率時(shí)，可以將其歸一化為：在二維\o"復(fù)平面"復(fù)平面上，在點(diǎn)的數(shù)值,因此通過(guò)\o"解析延拓（頁(yè)面不存在）"解析延拓:上述函數(shù)的\o"極點(diǎn)"極點(diǎn)等距離地分布在\o"半徑"半徑為的\o"圓"圓上由于在左半平面的極點(diǎn)即是的極點(diǎn)：分析巴特沃思低通濾波器的幅頻特性，不難發(fā)現(xiàn)，無(wú)論在通帶與阻帶內(nèi)都隨頻率而單調(diào)變化，濾波特性簡(jiǎn)單。在通帶內(nèi)誤差分布不均勻，靠近頻帶邊緣誤差最大。當(dāng)濾波器階數(shù)N較小時(shí)，阻帶幅頻特性下降較慢，與理想濾波器的特性相差較遠(yuǎn)。若要求阻帶特性下降迅速，則需增加濾波器的階數(shù)，設(shè)計(jì)該濾波器時(shí)所用元器件數(shù)量增多，線路也趨于復(fù)雜。若將誤差均勻地分布在通帶內(nèi)，就可以設(shè)計(jì)出階數(shù)較低的濾波器。這種誤差均勻分布的辦法可通過(guò)選擇具有等波紋特性的逼近函數(shù)來(lái)完成。2.1.3巴特沃斯低通濾波器階次的確定令2.2切比雪夫低通濾波器2.2.1切比雪夫低通濾波器簡(jiǎn)介切比雪夫?yàn)V波器是由切比雪夫多項(xiàng)式的正交函數(shù)推導(dǎo)出來(lái)的，采用了在通帶內(nèi)等波動(dòng)，在通帶外衰耗單調(diào)遞增的準(zhǔn)則去逼近理想濾波器特性。在通帶內(nèi)是等波紋的，在阻帶內(nèi)則是單調(diào)下降的，稱為切比雪夫Ⅰ型。在通帶內(nèi)是單調(diào)的，在阻帶內(nèi)是等波紋的，稱為切比雪夫Ⅱ型。2.2.2切比雪夫低通濾波器幅度特性（1）切比雪夫Ⅰ型濾波器切比雪夫型濾波器幅度響應(yīng)式如下：其中，紋波參數(shù)，且，其值越大，紋波越大；若截止頻率為，濾波器在截止頻率的放大率如下：

為

階切比雪夫多項(xiàng)式，表達(dá)式如下：

展開(kāi)成多項(xiàng)式如下：下面給出切比雪夫型濾波器不同階數(shù)下的幅度特性圖：圖2-2切比雪夫型濾波器不同階數(shù)下的幅度特性圖當(dāng)時(shí)，在1與之間等幅波動(dòng)，愈小，波動(dòng)幅度愈小。所有曲線在時(shí)通過(guò)點(diǎn)。當(dāng)時(shí)，若N為奇數(shù)，則=1；若n為偶數(shù)，則；通帶內(nèi)誤差分布是均勻的。當(dāng)時(shí)，曲線單調(diào)下降，N值愈大，曲線下降愈快。綜上可知Ⅰ型切比雪夫?yàn)V波器有三個(gè)參數(shù)需要確定：波動(dòng)系數(shù)，通帶截止頻率和階數(shù)N。通帶截止頻率一般按照實(shí)際要求給定；表示通帶內(nèi)最大損耗，由容許的通帶最大衰減確定。（2）切比雪夫Ⅱ型濾波器第二類切比雪夫?yàn)V波器也稱倒數(shù)切比雪夫?yàn)V波器，較不常用，因?yàn)轭l率截止速度不如第一類快，也需要用更多的電子元件。第二類切比雪夫?yàn)V波器在通頻帶內(nèi)沒(méi)有幅度波動(dòng)，只在阻頻帶內(nèi)有幅度波動(dòng)。第二類切比雪夫?yàn)V波器的轉(zhuǎn)移函數(shù)為：參數(shù)與

\o"阻頻帶（頁(yè)面不存在）"阻帶的\o"衰減度（頁(yè)面不存在）"衰減度有如下關(guān)系：

2.2.2切比雪夫低通濾波器階次的確定 切比雪夫?yàn)V波器的衰減函數(shù)定義為：其中為通帶最大衰減，又稱為通帶波紋，定義為：波動(dòng)系數(shù)為：由濾波器的通帶截止頻率及通帶內(nèi)允許的最大衰減和阻帶下限截止頻率及阻帶內(nèi)允許的最小衰減，可以確定濾波器所需的階數(shù)N。阻帶內(nèi)允許的最小衰減為濾波器的階次為：2.3兩種低通模擬濾波器的比較下圖為三種同階低通濾波器的幅度特性圖，可見(jiàn)巴特沃斯濾波器通帶和阻帶幅度特性均平坦，切比雪夫則是通帶平坦，另外阻帶有波動(dòng)，或者阻帶平坦，通帶有波動(dòng)。圖2-3三種同階低通濾波器幅度特性圖從幅頻特性上看，巴特沃斯濾波器在全頻段具有單調(diào)下降的幅度特性，切比雪夫Ⅰ型在通帶中呈等波紋形的，在阻帶則是單調(diào)下降的，切比雪夫Ⅱ型在通帶是單調(diào)下降的，在阻帶是呈等波紋形的，衰減特性依次為切比雪夫、巴特沃斯；通帶波紋依次為巴特沃斯、切比雪夫；巴特沃斯濾波器具有最平具有“最平幅度”特性。從過(guò)渡帶寬比較：具有相同階數(shù)（極點(diǎn)數(shù)）的各類濾波器中，切比雪夫的過(guò)渡帶要窄。在滿足同樣的通帶、阻帶衰減指標(biāo)的情況下，巴特沃斯的過(guò)渡帶要比切比雪夫的寬。從階次N比較：若濾波器具有相同的幅度特性指標(biāo)，則所需階次N依次為巴特沃斯型最大，切比雪夫次之。從濾波器對(duì)參數(shù)量化的靈敏度比較，量化靈敏度巴特沃斯型最低，切比雪夫Ⅰ、Ⅱ型較高。相位響應(yīng)比較：巴特沃斯在部分通帶中有線性相位，切比雪夫Ⅰ、Ⅱ型次之。巴特沃斯、切比雪夫?yàn)V波器是從濾波器的幅頻特性上考慮的，濾波器的相位特性較差，切比雪夫?yàn)V波器相位特性比巴特沃斯濾波器差。3.濾波器的設(shè)計(jì)方法3.1巴特沃斯（Butterworth）濾波器3.1.1巴特沃斯濾波器的幅度平方響應(yīng)巴特沃斯低通濾波器可用如下振幅的平方對(duì)頻率的公式表示：注：濾波器的階數(shù)；為一常數(shù)參數(shù)；為歸一化低通截止頻率,。N階巴特沃斯低通濾波器的傳遞函數(shù)的逼近函數(shù)為：N為電路階數(shù)，為截止頻率，為通帶增益。令，則歸一化幅頻響應(yīng)為：3.1.2設(shè)計(jì)巴特沃斯低通濾波器的指標(biāo)：通帶截止頻率；：通帶衰減，單位：；：阻帶起始頻率；：阻帶衰減，單位：。11巴特沃斯低通濾波器指標(biāo)圖3-1巴特沃斯模擬濾波器幅度特性說(shuō)明：衰減在這里以分貝（dB）為單位；即（2）當(dāng)時(shí)為通常意義上的截止頻率。（3）在濾波器設(shè)計(jì)中常選用歸一化的頻率，即3.1.3巴特沃斯低通濾波器的設(shè)計(jì)方法和步驟模擬濾波器（AF）的設(shè)計(jì)方法都?xì)w結(jié)于先設(shè)計(jì)一個(gè)“樣本”的歸一化原型低通濾波器，然后通過(guò)模擬頻帶變換得到巴特沃斯模擬濾波器.巴特沃斯低通濾波器的設(shè)計(jì)方法如下：計(jì)算歸一化頻率；計(jì)算巴特沃斯濾波器的參數(shù)和階次；利用查表獲得歸一化巴特沃斯低通原型濾波器的系統(tǒng)函數(shù)；將轉(zhuǎn)換為巴特沃斯低通濾波器的系統(tǒng)函數(shù)。具體步驟如下：（1）計(jì)算歸一化頻率將實(shí)際頻率歸一化得到：（2）根據(jù)已知的幅度平方函數(shù)確定和。由帶入得即因?yàn)?，所以由兩邊取?duì)數(shù)得：注意：當(dāng)時(shí)，即，此時(shí)巴特沃斯濾波器只剩下一個(gè)參數(shù)。（3）確定歸一化的巴特沃斯低通濾波器的傳遞函數(shù)。由于由解得極點(diǎn)為：將左半平面的極點(diǎn)賦予即其中也可以通過(guò)查階次表的方法得到：表3-1歸一化巴特沃斯模擬低通濾波器系統(tǒng)函數(shù)表N多項(xiàng)式因子Bn(s)1(s+1)2s2+1.414s+13(s+1)(s2+s+1)4(s2+0.7654s+1)(s2+1.8478s+1)5(s+1)(s2+0.6180s+1)(s2+1.6180s+1)6(s2+0.5176s+1)(s2+1.414s+1)(s2+1.9318s+1)7(s+1)(s2+0.4450s+1)(s2+1.247s+1)(s2+1.8022s+1)8(s2+0.3986s+1)(s2+1.111s+1)(s2+1.6630s+1)(s2+1.9622s+1)去掉歸一化影響上面設(shè)計(jì)中采用歸一化的頻率即，而實(shí)際中截止頻率為，所以要進(jìn)行如下的變量代換：即即可得到巴特沃斯低通濾波器的系統(tǒng)函數(shù)。3.2切比雪夫（Chebyshev）濾波器3.2.1切比雪夫?yàn)V波器的幅度平方響應(yīng)切比雪夫歸一化濾波器的幅度平方函數(shù)為注：通帶波紋參數(shù)，的正數(shù)，越大，波紋也越大；為歸一化截止頻率,；截止頻率，也是濾波器的通帶帶寬；（注：切比雪夫?yàn)V波器的通帶帶寬并不一定是3dB帶寬）為階切比雪夫多項(xiàng)式，定義為其中為反余弦函數(shù)；為雙曲余弦函數(shù)，；為反雙曲余弦函數(shù)，；將階切比雪夫多項(xiàng)式展開(kāi)為多項(xiàng)式的形式如表5-2所示：表3-2切比雪夫多項(xiàng)式展開(kāi)式N01x234567由表5-2可歸納出各階切比雪夫多項(xiàng)式的遞推公式為下圖為時(shí)切比雪夫多項(xiàng)式的特性。11-11x圖3-2N=0，4，5時(shí)切比雪夫多項(xiàng)式的特性可得：1.切比雪夫多項(xiàng)式的零值在的間隔內(nèi)。2.當(dāng)時(shí)，，且具有等波紋幅度特性。根據(jù)以上所述，在，在接近1處振蕩，其最大值為1，最小值為。在此范圍之外，隨著增大，，則很快接近于零。振幅特性的起伏為1～，因，所以在時(shí)，，所以切比雪夫=1\*ROMANI型濾波器的截止頻率并不對(duì)應(yīng)3dB的衰減。3.2.2設(shè)計(jì)切比雪夫?yàn)V波器的指標(biāo)：通帶截止頻率；：通帶衰減，單位：；：阻帶起始頻率；：阻帶衰減，單位：。3.2.3切比雪夫Ⅰ型濾波器的設(shè)計(jì)方法和步驟模擬濾波器（AF）的設(shè)計(jì)方法都?xì)w結(jié)于先設(shè)計(jì)一個(gè)“樣本”的歸一化原型低通濾波器，然后通過(guò)模擬頻帶變換得到切比雪夫?yàn)V波器.計(jì)算歸一化頻率；根據(jù)通帶波紋（通帶衰減），計(jì)算；根據(jù)阻帶起始頻率，阻帶衰減和。計(jì)算濾波器的階數(shù)；根據(jù)階數(shù)，查表得歸一化原型切比雪夫?yàn)V波器系統(tǒng)函數(shù)；根據(jù)的低頻特性求出待定系數(shù)；去掉歸一化影響，計(jì)算切比雪夫?yàn)V波器的系統(tǒng)函數(shù)；具體方法如下：（1）將實(shí)際頻率歸一化得（2）計(jì)算確定和。定義通帶波紋(即通帶衰減)（以分貝為單位）為：代入得即因?yàn)樗詣t其中這樣可以求出和，其中.（3）在已知的情況下，就可以根據(jù)幅度平方函數(shù)求出濾波器的零點(diǎn)和極點(diǎn)，從而確定濾波器的系統(tǒng)函數(shù)。也可以通過(guò)查表得到歸一化原型切比雪夫低通濾波器的系統(tǒng)函數(shù)。表3-3歸一化切比雪夫原型濾波器分母多項(xiàng)式設(shè)計(jì)系數(shù)n波紋（）12.862775221.51620261.425624530.71569381.53489541.252913040.37905061.02545531.71686621.197385650.17892340.75251811.30957471.93736751.1724909波紋（）11.965226721.10251031.097734330.49130671.23840920.988341240.27562760.74261941.45392480.952811450.12282670.58053420.97439611.68881600.9368201波紋（）11.307560320.63676810.803816430.32689011.02219030.737821640.20576510.51679811.25648190.716215050.08172250.45934910.69347701.49954330.7064606（4）去歸一化影響，系統(tǒng)函數(shù)為：對(duì)于切比雪夫Ⅱ型濾波器與Ⅰ型類似，考慮篇幅這里不做介紹了模擬低通濾波器轉(zhuǎn)換為數(shù)字濾波器的設(shè)計(jì)方法與流程由于模擬濾波器設(shè)計(jì)方法成熟，且有完整的表格數(shù)據(jù)可供使用，因而一般利用模擬濾波器，對(duì)其進(jìn)行數(shù)字化，即是從s平面映射到z平面，來(lái)間接得到所需的數(shù)字濾波器。4.1三種設(shè)計(jì)方法實(shí)際應(yīng)用中的數(shù)字濾波器有低通、高通、帶通、帶阻等類型。設(shè)計(jì)各類數(shù)字濾波器從歸一化原型的模擬低通濾波出發(fā)主要有以下三種方法：4.1.1方法1第一種方法分兩步實(shí)現(xiàn)：第一步是將設(shè)計(jì)出的歸一化樣本模擬低通濾波器經(jīng)模擬—模擬頻帶變換成各種模擬濾波器；第二步是數(shù)字化（采用沖激響應(yīng)不變法或者雙線性變換法）成各相應(yīng)頻帶的數(shù)字濾波器。設(shè)計(jì)流程框圖如下：4.1.2方法2直接導(dǎo)出由設(shè)計(jì)出的樣本模擬低通濾波器變換成各種通帶數(shù)字濾波器，則是把頻帶變換和數(shù)字化結(jié)合起來(lái)，一步變換完成設(shè)計(jì)。設(shè)計(jì)流程框圖如下：4.1.3方法3將設(shè)計(jì)出的歸一化的樣本模擬低通濾波器先數(shù)字化成低通數(shù)字濾波器，然后用數(shù)字—數(shù)字頻帶變換法設(shè)計(jì)出各種通帶的數(shù)字濾波器。設(shè)計(jì)流程框圖如下：4.2設(shè)計(jì)步驟利用上面的設(shè)計(jì)方法將模擬低通濾波器設(shè)計(jì)成IIR數(shù)字濾波器的設(shè)計(jì)步驟如下：將給定的數(shù)字濾波器的性能指標(biāo)，按某一變換(映射)規(guī)則轉(zhuǎn)換成相應(yīng)的模擬流波器的性能指標(biāo)。即是把平面映射到平面，使模擬系統(tǒng)函數(shù)變換成所需的數(shù)宇濾波器的系統(tǒng)函數(shù)；用所得到的模擬低通濾波器的性能指標(biāo)利用某種模擬濾波器逼近方法設(shè)計(jì)并查表求得此模擬低通濾波器的系統(tǒng)函數(shù)，以它作為設(shè)計(jì)數(shù)字濾波器的“樣本”。將作為“樣本”的模擬原型低通濾波器的系統(tǒng)函數(shù)最終變換成所需的數(shù)字各型濾波器的系統(tǒng)函數(shù)。注：步驟一種這種由復(fù)變量到復(fù)變量之間的映射(變換)關(guān)系，必須滿足以下兩條基本要求:的頻率響應(yīng)須有可能模仿的頻率響應(yīng)，即平面的虛軸必須映射到平面的單位圓上，也就是頻率軸要對(duì)應(yīng)。因果穩(wěn)定的應(yīng)能映射成因果穩(wěn)定的。也就是平面的左半平面必須映射到平面單位圓的內(nèi)部。5.切比雪夫模擬低通濾波器的設(shè)計(jì)5.1設(shè)計(jì)技術(shù)指標(biāo)及要求下通帶截止頻率為，上通帶截止頻率為；阻帶下截止頻率，阻帶上截止頻率為；通帶最大衰減率：；阻帶最小衰減率：；抽樣頻率；要求：采用從低通濾波器直接變換到帶阻濾波器的辦法設(shè)計(jì)5.2設(shè)計(jì)方法選擇和原理簡(jiǎn)介本題采用的設(shè)計(jì)方法是從模擬低通濾波器直接變換到數(shù)字帶阻濾波器。該種方法的設(shè)計(jì)原理是把一個(gè)歸一化原型摸擬低通濾波器經(jīng)模擬頻帶變換成所需要類型(包括高通、帶通、帶阻與另一截止頻率的低通)的模擬濾波器，然后再通過(guò)沖激響應(yīng)不變法或雙線性變換法數(shù)字化為所需類型的字濾波器。實(shí)際設(shè)計(jì)時(shí)是把這一方法中的兩步合成一步來(lái)實(shí)現(xiàn)，即把模擬歸一化低通原型變換到模擬低通、高通、帶通、帶阻等濾波器的公式與用雙線性變換得到相應(yīng)數(shù)字濾波器的公式合并，就可直接從模擬低通歸一化原型通過(guò)一定的頻率變換關(guān)系，一步完成各類型數(shù)字濾波器的設(shè)計(jì)，因沖激響應(yīng)不變法有頻率響應(yīng)混疊失真效應(yīng)，只適用于設(shè)計(jì)嚴(yán)格限帶的數(shù)字低通、帶通濾波器;而對(duì)于數(shù)字高通、帶阻濾波器，則不能直接應(yīng)用。在本題中要求設(shè)計(jì)的是帶阻濾波器，所以這里采用的是雙線性變換法。5.3設(shè)計(jì)流程圖5.4設(shè)計(jì)過(guò)程按參數(shù)要求計(jì)算出階數(shù)、截止頻率或通帶紋波系數(shù)、階數(shù)等，寫(xiě)出模擬低通濾波器的表達(dá)式：（1）各對(duì)應(yīng)的數(shù)字域臨界頻率（不必預(yù)畸）為求此帶阻數(shù)字低通濾波器的中心頻率（3）雙線性變換法需要將各數(shù)字域臨界頻率預(yù)畸為各模擬域臨界頻率，求樣本低通AF的可得取注：在較小頻率處滿足衰減要求，則較大頻率處衰減一定會(huì)更大，超過(guò)要求，所以此處選擇（4）求樣本低通AF的階次，利用公式由于可求出，擇取取階數(shù)。查詢切貝雪夫?yàn)V波器分母多項(xiàng)式的系數(shù)的表格，當(dāng)，，得切貝雪夫?yàn)V波器分母多項(xiàng)式的系數(shù)：則歸一化原型切貝雪夫型低通模擬濾波器的為（6）由于是切貝雪夫型濾波器，故，可直接利用“去歸一化”得到樣本低通模擬低通濾波器的系統(tǒng)函數(shù)：其中利用雙線性變換關(guān)系式求所需帶阻數(shù)字濾波器的系統(tǒng)函數(shù)將模擬低通變換為數(shù)字帶阻時(shí)，域到域的映射關(guān)系應(yīng)該是映射為映射為滿足這一變換要求的雙線性變換為先設(shè)其中的，則所求的數(shù)字濾波器的表達(dá)式為MATLAB濾波器仿真及結(jié)果分析6.1MATLAB自編程設(shè)計(jì)濾波器設(shè)計(jì)步驟運(yùn)用自編程序設(shè)計(jì)與仿真的步驟如下：給定技術(shù)指標(biāo)：通帶截止頻率,阻帶截止頻率,通帶衰減,阻帶衰減。求濾波器階次及設(shè)計(jì)用的通帶截止頻率利用工具箱中，求歸一化原型低通濾波器。求非歸一化的低通濾波器。6.2MATLAB仿真結(jié)果與分析1.首先都得切比雪夫Ⅰ型數(shù)字帶通幅度響應(yīng)圖像如下；圖6-1切比雪夫帶阻幅度響應(yīng)圖從幅度響應(yīng)圖可以明顯看出此濾波器類型為帶阻濾波器，下面我們就用不同頻率的信號(hào)來(lái)檢測(cè)此濾波器是否為帶阻濾波器。2.選擇如下低頻、中頻、高頻信號(hào)三組信號(hào)圖6-2低頻、中頻、高頻信號(hào)波形3.將低頻信號(hào)與中頻信號(hào)疊加作為輸入信號(hào)，觀測(cè)輸出信號(hào)如下：圖6-3低頻、中頻疊加信號(hào)輸入與輸出波形疊加信號(hào)在通過(guò)濾波器過(guò)后中頻段信號(hào)基本被過(guò)濾，通過(guò)濾波器的為低頻信號(hào)。4.將中頻信號(hào)與高頻信號(hào)疊加作為輸入信號(hào)，觀測(cè)輸出信號(hào)如下：圖6-4中頻、高頻疊加信號(hào)輸入與輸出波形疊加信號(hào)在通過(guò)濾波器過(guò)后中頻段信號(hào)基本被過(guò)濾，通過(guò)濾波器的為高頻信號(hào)。根據(jù)以上的不同頻率信號(hào)疊加后的輸出分析可知，此濾波器滿足帶阻濾波器的要求。根據(jù)程序運(yùn)行可得以下濾波器參數(shù)：N=2b=0.2489-0.22020.5465-0.22020.2489a=1.0000-0.4127-0.0734-0.14190.3880C=0.6414B=-0.2455-0.4234-0.14700.2718A=1.00000.96230.57531.0000-1.37500.6745dbHx=2.000016.599545.04582.0000與計(jì)算結(jié)果相比存在一定差異，帶阻頻率范圍相較于計(jì)算值增大，主要體現(xiàn)在上通帶截止頻率增大，同時(shí)導(dǎo)致中心截至頻率增大。下圖為模擬仿真結(jié)果：圖6-5切比雪夫Ⅰ型模擬帶阻幅度響應(yīng)與數(shù)字仿真結(jié)果相比較更為接近計(jì)算值，同時(shí)發(fā)現(xiàn)在適當(dāng)提高數(shù)字濾波器的抽樣頻率后可以使數(shù)字濾波器濾波器的仿真結(jié)果更加接近于實(shí)際計(jì)算值。附件主程序：Fs=2*10^4;wp1=2*2000/Fs;wp2=2*7000/Fs;%切比雪夫Ⅰ型數(shù)字帶阻技術(shù)指標(biāo)ws1=2*3000/Fs;ws2=2*4000/Fs;%適當(dāng)提高帶阻下截止頻率，避免與帶通下截止頻率重合wp=[wp1,wp2];ws=[ws1,ws2];Rp=2;As=20;[N,wc]=cheb1ord(wp,ws,Rp,As);%階次[b,a]=cheby1(N,Rp,wc,'stop');[C,B,A]=tf2par(b,a);%直接型轉(zhuǎn)化為并聯(lián)型w0=[wp1*pi,ws1*pi,ws2*pi,wp2*pi];%檢驗(yàn)各臨界頻率的衰減Hx=freqz(b,a,w0);[H,w]=freqz(b,a);dbHx=-20*log10(abs(Hx)/max(abs(H)));[db,mag,pha,grd,w]=freqz_m(b,a);%求頻率響應(yīng)各項(xiàng)信息figure;plot(w/pi,db);title('切比雪夫Ⅰ型數(shù)字帶阻幅度響應(yīng)（db）');xlabel('\omega/\pi');ylabel('db');axis([0,1,-70,10]);set(gca,'xtickmode','manual','xtick',0:0.1:1);set(gca,'ytickmode','manual','ytick',-70:5:10);grid;n=0:200;x1=cos(0.1*pi*n);x2=cos(0.4*pi*n);x3=cos(pi*n);%x1為低頻信號(hào)，x2為中頻信號(hào)，x3為高頻信號(hào)figure;%分別繪制出三種頻段的抽樣信號(hào)圖像subplot(3,1,1);plot(n,x1);title('cos(0.1*pi*n)');subplot(3,1,2);plot(n,x2);title('cos(0.4*pi*n)');subplot(3,1,3);plot(n,x3);title('cos(pi*n)');figure;%繪制低頻信號(hào)與中評(píng)信號(hào)疊加的輸入輸出圖形y1=filter(b,a,x1+x2);subplot(2,1,1);plot(n,x1+x2);title('cos(0.1*pi*n)+cos(0.4*pi*n)');subplot(2,1,2);plot(n,y1);title('輸出信號(hào)');figure;%繪制高頻信號(hào)與中評(píng)信號(hào)疊加的輸入輸出圖形y2=filter(b,a,x2+x3);subplot(2,1,1);plot(n,x2+x3);title('cos(0.4*pi*n)+cos(pi*n)');subplot(2,1,2);plot(n,y2);title('輸出信號(hào)');NbaCBAdbHx輸出：N=2b=0.2489-0.22020.5465-0.22020.2489a=1.0000-0.4127-0.0734-0.14190.3880C=0.6414B=-0.2455-0.4234-0.14700.2718A=1.00000.96230.57531.0000-1.37500.6745dbHx=2.000016.599545.04582.0000子程序-tf2par.mfunction[C,B,A]=tf2par(b,a)M=length(b);N=length(a);[r1,p1,C]=residuez(b,a);%先求系統(tǒng)的單根p1,對(duì)應(yīng)的留數(shù)r1及直接項(xiàng)Cp=cplxpair(p1,1e-9);%先配對(duì)函數(shù)cplxpair由p1找共軛復(fù)根p，le-9的誤差I(lǐng)=cplxcomp(p1,p);%找p1變?yōu)閜的排序變化r=r1(I);%變?yōu)?階子系統(tǒng)K=floor(N/2);B=zeros(K,2);A=zeros(K,3);ifK*2==N;fori=1:2:N-2pi=p(i:i+1,:);%取出一對(duì)極點(diǎn)ri=r(i:i+1,:);%取出一對(duì)對(duì)應(yīng)留數(shù)[Bi,Ai]=residuez(