畢業(yè)設(shè)計---ii數(shù)字低通濾波器的設(shè)計

1、<p>　　學(xué) 士 學(xué) 位 論 文</p><p>　　題 目 II數(shù)字低通濾波器的設(shè)計</p><p><b>　　學(xué) 生 </b></p><p><b>　　指導(dǎo)教師 </b></p><p><b>　　年 級 </b><

2、;/p><p>　　專 業(yè) 計算機科學(xué)與技術(shù)</p><p>　　系 別 信息工程系</p><p>　　學(xué) 院 計算機科學(xué)與信息工程</p><p><b>　　摘 要</b></p><p>　　摘 要：本論文首先介紹了濾波器的濾波原理以及模擬濾波器、數(shù)字濾波器的設(shè)

3、計方法。重點介紹了IIR數(shù)字濾波器的設(shè)計方法。即脈沖響應(yīng)不變法和雙線性變換法。在此基礎(chǔ)上，用DSP虛擬實現(xiàn)任意階IIR數(shù)字濾波器。此設(shè)計擴展性好，便于調(diào)節(jié)濾波器的性能，可以根據(jù)不同的要求在DSP上加以實現(xiàn)。</p><p>　　關(guān)鍵字： DSP 濾波、IIR數(shù)字低通濾波器</p><p><b>　　目 錄</b></p><p><b

4、>　　摘 要2</b></p><p><b>　　目 錄2</b></p><p>　　第1章 濾波器簡介3</p><p>　　1.1 濾波器的工作原理4</p><p>　　1.1.1 模擬濾波器的工作原理4</p><p>　　1.1.2 數(shù)字濾波器的工作

5、原理6</p><p>　　1.2 濾波器的基本特性7</p><p>　　1.2.1 模擬濾波器與數(shù)字濾波器的基本特性7</p><p>　　1.2.2 無限沖擊響應(yīng)IIR和有限沖擊響應(yīng)FIR濾波器9</p><p>　　1.3 濾波器的主要技術(shù)指標(biāo)10</p><p>　　第2章 模擬濾波器的設(shè)

6、計11</p><p>　　2.1 模擬濾波器的設(shè)計方法12</p><p>　　2.2 模擬原型濾波器及最小階數(shù)的選擇14</p><p>　　2.2.1 巴特沃斯濾波器及最小階數(shù)的選擇14</p><p>　　2.2.2低通原型濾波器的系統(tǒng)函數(shù)20</p><p>　　2.2.3 橢圓濾波器及最小階數(shù)

7、的選擇21</p><p>　　2.2.4貝塞爾濾波器22</p><p>　　第3章 IIR數(shù)字濾波器的設(shè)計23</p><p>　　3.1 IIR數(shù)字濾波器的設(shè)計方法23</p><p>　　3.2 IIR濾波器經(jīng)典設(shè)計24</p><p>　　3.3 IIR濾波器直接設(shè)計33</p>

8、<p>　　第4章 DSP仿真系統(tǒng)34</p><p>　　4.1 對低通模擬和數(shù)字濾波器的仿真34</p><p>　　4.1.1 模擬低通濾波器的仿真34</p><p>　　4.2.2 數(shù)字低通濾波器的仿真36</p><p>　　4.4.1 模擬帶通濾波器的仿真39</p><p&g

9、t;　　4.4.2 數(shù)字帶通濾波器的設(shè)計40</p><p>　　4.5 對帶阻模擬和數(shù)字濾波器的仿真42</p><p>　　4.5.1 模擬帶阻濾波器的設(shè)計42</p><p>　　4.5.2 數(shù)字帶阻濾波器的仿真43</p><p>　　第5章 總結(jié)與展望47</p><p>　　5.1 總 結(jié)47

10、</p><p>　　5.2 展 望47</p><p><b>　　參考文獻(xiàn)49</b></p><p>　　第1章 濾波器簡介</p><p>　　從廣義上講，任何對某些頻率（相對于其他頻率來說）進(jìn)行修正的系統(tǒng)稱為濾波器。嚴(yán)格地講，對輸入信號通過一定的處理得到輸出信號，這個處理通常是提取信號中某頻率范圍內(nèi)的信

11、號成分，把這種處理的過程稱為濾波。實現(xiàn)濾波處理的運算電路或設(shè)備稱為濾波器。</p><p>　　在許多科學(xué)技術(shù)領(lǐng)域中，廣泛應(yīng)用線性濾波和頻譜分析對信號進(jìn)行加工處理，模擬濾波是處理連續(xù)信號，數(shù)字濾波則是處理離散信號，而后者是在前者的基礎(chǔ)上發(fā)展起來的。我們知道，無源或有源模擬濾波器是分立元件構(gòu)成的線性網(wǎng)絡(luò)，他們的性能可以用線性微分方程來描述，而數(shù)字濾波器是個離散線性系統(tǒng)，要用差分方程來描述，并以離散變換方法來分析。這

12、些方程組可以用專用的或通用的數(shù)字計算機進(jìn)行數(shù)字運算來實現(xiàn)。因此，數(shù)字濾波器的濾波過程是一個計算過程，它將輸入信號的序列數(shù)字按照預(yù)定的要求轉(zhuǎn)換成輸出數(shù)列。</p><p>　　1.1 濾波器的工作原理</p><p>　　1.1.1 模擬濾波器的工作原理</p><p>　　我們知道，模擬濾波器是對模擬信號實行線性濾波的一種線性時不變系統(tǒng)，如圖1.1所示。在時域內(nèi)，它

13、的動態(tài)特性可以用系統(tǒng)的單位沖激函數(shù)的響應(yīng)來描述，也就是該濾波系統(tǒng)在任何時刻對輸入單位沖激信號=δ（t）的輸出響應(yīng)。這個函數(shù)從時域上反映了該濾波系統(tǒng)的傳輸特性。對于任意輸入信號，系統(tǒng)的輸出可以卷積表示：</p><p>　　= (1.1)</p><p>　　上式表明在對線性濾波器系統(tǒng)進(jìn)行時域分析時，采用了疊加原理，先將任意輸入信號波形分成不同時間的窄脈沖之和，

14、再分別求出各個脈沖通過濾波器之后的響應(yīng)，并進(jìn)行線性疊加從而得到總的輸出信號。</p><p>　　圖1.1模擬濾波器原理</p><p>　　在頻域分析時，線性濾波器的轉(zhuǎn)移函數(shù)等于系統(tǒng)的單位沖激函數(shù)的響應(yīng)的拉普拉斯變換：</p><p><b>　　（1.2）</b></p><p>　　很明顯，當(dāng)s=jω，上式就是傅立

15、葉變換的表達(dá)式，它反映了濾波器的傳輸特性對各種頻率的響應(yīng)，也就是濾波器的頻率響應(yīng)函數(shù)，它決定著濾波特性。當(dāng)濾波器輸入信號與輸出信號的拉普拉斯變換，得</p><p><b>　　（1.3）</b></p><p>　　這表明兩信號卷積的變換等于各自變換的乘積。在頻譜關(guān)系上，一個輸入信號的頻譜，經(jīng)過濾波器的作用后，被變換成的頻譜。因此，根據(jù)不同的濾波要求來選定，就可以得

16、到不同類型的模擬濾波器。還可以看出，濾波器的濾波過程就是完成信號與它的單位沖激函數(shù)響應(yīng)之間的數(shù)學(xué)卷積運算過程。</p><p>　　1.1.2 數(shù)字濾波器的工作原理</p><p>　　在數(shù)字濾波中，我們主要討論離散時間序列。如圖1.2所示。設(shè)輸入序列為，離散或數(shù)字濾波器對單位抽樣序列的響應(yīng)為。因在時域離散信號和系統(tǒng)中所起的作用相當(dāng)于單位沖激函數(shù)在時域連續(xù)信號和系統(tǒng)中所起的作用。</

17、p><p>　　圖1.2 數(shù)字濾波器原理</p><p>　　數(shù)字濾波器的序列將是這兩個序列的離散卷積，即</p><p><b>　?。?.4） </b></p><p>　　同樣，兩個序列卷積的z變換等于個自z變換的乘積，即</p><p><b>　?。?.5）</b>

18、</p><p>　　用代入上式，其中T為抽樣周期，則得到</p><p><b>　?。?.6）</b></p><p>　　式中和 分別為數(shù)字濾波器輸入序列和輸出序列的頻譜，而為單位抽樣序列響應(yīng)的頻譜。由此可見，輸入序列的頻譜經(jīng)過濾波后，變?yōu)?，按照的特點和我們處理信號的目的，選取適當(dāng)?shù)氖沟臑V波后的符合我們的要求。</p>

19、<p>　　1.2 濾波器的基本特性</p><p>　　1.2.1 模擬濾波器與數(shù)字濾波器的基本特性</p><p>　　如利用模擬電路直接對模擬信號進(jìn)行處理則構(gòu)成模擬濾波器，它是一個連續(xù)時間系統(tǒng)。如果利用離散時間系統(tǒng)對數(shù)字信號（時間離散、幅度量化的信號）進(jìn)行濾波則構(gòu)成數(shù)字濾波器。</p><p>　　數(shù)字濾波器的差分方程表示為：</p>

20、<p><b>　　系統(tǒng)函數(shù)表示:</b></p><p>　　數(shù)字濾波器的特性通常用其頻率響應(yīng)函數(shù)來描述，</p><p>　　包括幅度特性和相位特性。</p><p>　　按信號通過系統(tǒng)時的特性（主要是幅頻特性）來分類：可以有低通、高通、帶通和帶阻四種基本類型。</p><p>　　低通數(shù)字濾波器：圖1.

21、3所示</p><p>　　圖1.3 低通數(shù)字濾波器的頻譜</p><p>　?。?） 高通數(shù)字濾波器：圖1.4所示</p><p>　　圖1.4 高通數(shù)字濾波器的頻譜</p><p>　?。?） 帶通數(shù)字濾波器：圖1.5所示</p><p>　　圖1.5 帶通數(shù)字濾波器的頻譜</p><p&

22、gt;　?。?）帶阻數(shù)字濾波器：圖1.6所示</p><p>　　圖1.6 帶阻數(shù)字濾波器的頻譜</p><p>　　其他較復(fù)雜的特性可以由基本濾波器組合。</p><p>　　1.2.2 無限沖擊響應(yīng)IIR和有限沖擊響應(yīng)FIR濾波器</p><p>　　按系統(tǒng)沖擊響應(yīng)（或差分方程）可以分成無限沖擊響應(yīng) IIR和有限沖擊響應(yīng)FIR濾波器兩

23、類。這兩種濾波器都可以現(xiàn)實各種頻率特性要求，但它們在計算流程、具體特性逼近等方面是有差別的。</p><p>　　(1) FIR濾波器(非遞歸型):</p><p>　　(2) IIR濾波器（遞歸型）</p><p>　　還有一些其他的分類方法，例如在特定場合使用的濾波器。</p><p>　　1.3 濾波器的主要技術(shù)指標(biāo)</p>

24、;<p>　　濾波器的主要技術(shù)指標(biāo)取決于具體的應(yīng)用或相互間的相互關(guān)系。具體的有最大通帶增益（即通帶允許起伏）；最大阻帶增益；通帶截止頻率；阻帶截止頻率。如圖1.7所示</p><p>　　第2章 模擬濾波器的設(shè)計</p><p>　　模擬濾波器的理論和設(shè)計方法已經(jīng)發(fā)展的相當(dāng)成熟，且有若干典型的模擬濾波器供我們選擇，如巴特沃斯（Butterworth濾波器.切比雪夫（Che

25、byshev）濾波器等。這些工作的理論分析和設(shè)計方法在20世紀(jì)30年代就完成，然而煩瑣.冗長的數(shù)字計算使它難以付諸實用。直到50年代，由于計算機技術(shù)的逐步成熟，求出大量設(shè)計參數(shù)和圖表，這種方法才得到廣泛應(yīng)用。這些典型的濾波器各有特點：巴特沃斯濾波器具有單調(diào)下降的幅頻特性；切比雪夫濾波器的幅頻特性在通帶或者阻帶有波動發(fā)，可以提高選擇性。這樣根據(jù)具體要求可以選擇不同類型的濾波器。</p><p>　　模擬濾波器按幅度

26、特征可以分成低通、高通、帶通和帶阻濾波器。它們的理想幅度特性如圖2.1所示，但我們設(shè)計濾波器時，總是先設(shè)計低通濾波器，再通過頻率變換將低通濾波器轉(zhuǎn)換成希望類型的濾波器</p><p>　　2.1 模擬濾波器的設(shè)計方法</p><p>　　利用頻率變換設(shè)計模擬濾波器的步驟為：</p><p>　　（1）給定模擬濾波器的性能指標(biāo)，如截止頻率或上、下邊界頻率等。<

27、/p><p>　?。?）確定濾波器階數(shù)</p><p>　?。?）設(shè)計模擬低通原型濾波器。</p><p>　?。?）按頻率變換設(shè)計模擬濾波器（低通、高通、帶通、帶阻）。</p><p>　　模擬低通濾波器的設(shè)計指標(biāo)有，和，其中和分別稱為通帶截止頻率和阻帶截止頻率。 是通帶Ω（=0—）中的最大衰減系數(shù)，是阻帶Ω≥的最小衰減系數(shù)，和一般用dB表示。

28、對于單調(diào)下降的幅度特性，可表示成：</p><p><b>　?。?.1）</b></p><p><b>　　（2.2）</b></p><p>　　如果Ω=0處幅度已歸一化為一，即，和表示為</p><p><b>　?。?.3）</b></p><p&

29、gt;<b>　?。?.4）</b></p><p>　　以上技術(shù)指標(biāo)用圖2.2表示，圖中稱為3dB 截止頻率，因,-20 </p><p>　　圖2.2 低通濾波器的幅度特性</p><p>　　濾波器的技術(shù)指標(biāo)給定以后，需要設(shè)計一個傳輸函數(shù)，希望其幅度平方函數(shù)滿足給定的指標(biāo)和，一般濾波器的單位沖激響應(yīng)為實數(shù)，因此</p>&l

30、t;p>　　= (2.5)</p><p>　　如果能由，，，求出，那么就可以求出所需的，對于上面介紹的典型濾波器，其幅度平方函數(shù)有自己的表達(dá)式，可以直接引用。這里要說明的是必須是穩(wěn)定的。因此極點必須落在s平面的左半平面，相應(yīng)的的極點落在右半平面。</p><p>　　2.2 模擬原型濾波器及最小階數(shù)的選擇</p><p>

31、;　　2.2.1 巴特沃斯濾波器及最小階數(shù)的選擇</p><p>　　巴特沃斯濾波器是最基本的逼近方法形式之一。它的幅頻特性模平方為</p><p><b>　　(2.6)</b></p><p>　　式中N是濾波器的階數(shù)。當(dāng)Ω=0時，；當(dāng)Ω=時，，是3dB截止頻率。</p><p>　　不同階數(shù)N的巴特沃斯濾波器特性如

32、圖2.3所示，這一幅頻特性具有下列特點：</p><p>　?。?）最大平坦性：可以證明：在Ω=0點，它的前（2N-1）階導(dǎo)數(shù)都等于0，這表明巴特沃斯濾波器在Ω=0附近一段范圍內(nèi)是非常平直的，它以原點的最大平坦性來逼近理想低通濾波器?！白钇巾憫?yīng)”即由此而來。</p><p>　?。?）通帶，阻帶下降的單調(diào)性。這種濾波器具有良好的相頻特性。</p><p>　　（3）

33、3dB的不變性：隨著N的增加，頻帶邊緣下降越陡峭，越接近理想特性，但不管N是多少，幅頻特性都通過-3dB點。當(dāng)Ω≥時，特性以20NdB/dec速度下降。</p><p>　　圖2.3 不同階數(shù)N的巴特沃斯濾波器特性</p><p>　　現(xiàn)根據(jù)式（2.6）求巴特沃斯濾波器的系統(tǒng)函數(shù)Ha（s）。令Ω=s/j，帶入式（2.6）</p><p><b>　　對應(yīng)的

34、極點：</b></p><p><b>　　(2.7)</b></p><p>　　即為的極點，此極點分布有下列特點：</p><p>　　(1)的2N個極點以π/N為間隔均勻分布在半徑為的圓周上，這個圓稱為巴特沃斯圓。</p><p>　　(2)所有極點以jΩ軸為對稱軸成對稱分布，jΩ軸上沒有極點。<

35、/p><p>　　(3)當(dāng)N為奇數(shù)時，有兩個極點分布在的實軸上；N為偶函數(shù)時，實軸上沒有極點。所有復(fù)數(shù)極點兩兩呈共軛對稱分布。圖2.4畫出了N=3時的極點分布。全部零點位于s=∞處。</p><p>　　圖2.4 N=3時Ha（s）Ha（-s）極點分布</p><p>　　為得到穩(wěn)定的，取全部左半平面的極點。</p><p><b>

36、;　　（2.8）</b></p><p><b>　　當(dāng)N為偶數(shù)時</b></p><p><b>　?。?.9）</b></p><p><b>　　當(dāng)N為奇數(shù)時</b></p><p><b>　　（2.10）</b></p>

37、<p>　　為使用方便把式（2.9）和式（2.10）對進(jìn)行歸一化處理，為此，分子分母各除以，并令，稱為歸一化復(fù)頻率：</p><p>　?。∟為偶數(shù)） （2.11）</p><p>　?。∟為奇數(shù)）（2.12）</p><p>　　用歸一化頻率表示的頻率特性稱為原型濾波特性（Ω?即歸一化復(fù)頻率s? 的虛部）。對式（2.6）所示的低通巴特沃斯特性用Ω?表

38、示得到：</p><p><b>　?。?.13）</b></p><p>　　稱為巴特沃斯低通原型濾波器幅頻特性。在低通原型濾波頻率特性上，截止頻率=1。</p><p>　　若給出模擬低通濾波器的設(shè)計性能指標(biāo)要求：通帶邊界頻率，阻帶邊界頻率，通帶波紋，阻帶衰減，要確定butterworth ,，低通濾波器最小階數(shù)N及截止頻率。，，，的意義如

39、圖所示。</p><p>　　當(dāng)=時， 即，以截至頻率（幅值下降3dB）為1，化為相對為相對的相對頻率由上式可寫為。</p><p>　　同理，當(dāng)=時，， 。</p><p>　　由此可見 N應(yīng)向上取整，再用MATLAB 編程計算濾波器最小階數(shù)N和截止頻率。</p><p>　　就是切比雪夫濾波器的極點，給定N，，ε即可求的2N個極點分布。

40、由式（2.22）實部與虛部的正弦和余弦函數(shù)平方約束關(guān)系可以看出，此極點分布滿足橢圓方程，其短軸和長軸分別為</p><p><b>　　（2.23）</b></p><p>　　圖2.7畫出了N=3時切比雪夫濾波器的極點分布。</p><p>　　極點所在的橢圓可以和半徑為a的圓和半徑為b的圓聯(lián)系起來，這兩個圓分別稱為巴特沃斯小圓和巴特沃斯大圓

41、。N階切比雪夫濾波器極點的縱坐標(biāo)，而橫坐標(biāo)等于N階巴特沃斯小圓極點的橫坐標(biāo)取左半平面的極點：</p><p>　　k=1,2…，N （2.24）</p><p>　　則切比雪夫濾波器的系統(tǒng)函數(shù)：</p><p><b>　?。?.25）</b></p><p>　　其中，常數(shù)A=。因而切比雪夫濾波器的系統(tǒng)函數(shù)

42、表示為：</p><p><b>　?。?.27）</b></p><p>　　切比雪夫濾波器的截止角頻率不是像巴特沃斯濾波器中所規(guī)定的（-3dB）處角頻率，而是通帶邊緣的頻率。若波紋參數(shù)滿足，可以求的-3dB處的角頻率為</p><p><b>　?。?.28）</b></p><p>　　將式（

43、2.27）表示的對歸一化，得到切比雪夫I型</p><p>　　2.2.2低通原型濾波器的系統(tǒng)函數(shù)</p><p><b>　?。?.29）</b></p><p>　　對不同的N，式（2.29）的分母多項式已制成表格，供設(shè)計參考。</p><p>　　和butterworth低通模擬濾波器設(shè)計一樣，若給定性能指標(biāo)要求：

44、，，，確定Chebyshev低通模擬濾波器最小階數(shù)N和截止頻率（-3dB頻率）。</p><p>　　2.2.2.1 Chbbyshev I型</p><p>　　由式可得 故階數(shù)N可由下式求得 式中，，截至頻率由上面兩式用Matlab 編程計算濾波器最小階數(shù)N和截止頻率 </p><p>　　2.2.2.2 Chbbyshev II型</p>

45、<p>　　Chbbyshev II型通帶內(nèi)是平滑的，而阻帶具有等波紋起伏特性。因此，在階數(shù)N的計算公式上是相同的，而-3dB截止頻率則不同。</p><p>　　2.2.3 橢圓濾波器及最小階數(shù)的選擇</p><p>　　橢圓的模擬低通濾波器圓形的平方幅值響應(yīng)函數(shù)為</p><p>　　式中，為小于1的正書，表示波紋情況；為截止頻率；）為橢圓函數(shù)，定義為

46、</p><p>　　當(dāng)N為偶數(shù)（N=2m）時,</p><p>　　當(dāng)N為奇數(shù)（N=2m+1）時, 其中</p><p>　　橢圓模擬濾波器特點是：在通帶和阻帶內(nèi)均具有等波紋起伏特性。何以上濾波器相比，相同的性能指標(biāo)所需要的階數(shù)最小。但頻率響應(yīng)應(yīng)具有明顯的非線性。由式</p><p>　　濾波器的階數(shù)可由下式確定， ， 式中

47、由上式計算濾波器的最小階數(shù)N和截止頻率。</p><p>　　2.2.4貝塞爾濾波器</p><p>　　貝塞爾模擬低通濾波器原型的特點是在零頻時具有最平坦的群延遲，并在整個通帶內(nèi)延遲幾乎不變。在零頻時的群延遲為。由于這一特點，貝塞爾模擬濾波器通帶內(nèi)保持信號形狀不變。濾波器傳遞漢書具有下面形式</p><p><b>　　]</b></p

48、><p>　　第3章 IIR數(shù)字濾波器的設(shè)計</p><p>　　3.1 IIR數(shù)字濾波器的設(shè)計方法</p><p>　　IIR濾波器是一種數(shù)字濾波器，濾波器的系統(tǒng)函數(shù)如式 （3.1）</p><p>　　由于它的脈沖響應(yīng)序列是無限長的，故稱無限沖激響應(yīng)濾波器。IIR濾波器的設(shè)計就是根據(jù)濾波器某些性能指標(biāo)要求，設(shè)計濾波器的分子和分母多項

49、式。它和FIR濾波器相比優(yōu)點是在滿足相同性能指標(biāo)要求條件下，IIR濾波器的階數(shù)要明顯低于FIR濾波器。但I(xiàn)IR濾波器的相位是非線性的。</p><p>　　IIR濾波器設(shè)計方法可分為三種：模擬濾波器變換（經(jīng)典設(shè)計法）、直接設(shè)計法、參數(shù)模型設(shè)計法、最大平滑濾波器設(shè)計。</p><p>　　3.2 IIR濾波器經(jīng)典設(shè)計</p><p>　　基于模擬濾波器變換原理，首先

50、是根據(jù)濾波器的技術(shù)指標(biāo)設(shè)計出相應(yīng)的模擬濾波器，然后再將設(shè)計好的模擬濾波器變換成滿足給定技術(shù)指標(biāo)的數(shù)字濾波器。這就是IIR濾波器設(shè)計的經(jīng)典法。經(jīng)典法可設(shè)計出低通、高通、帶通、帶阻等各種IIR濾波器。在MATLAB中，經(jīng)典法設(shè)計IIR數(shù)字濾波器采用下面的主要步驟：</p><p>　　由上可見，經(jīng)典設(shè)計法是利用模擬濾波器的設(shè)計成果。第二步完成后，一個達(dá)到期望性能指標(biāo)的模擬濾波器（低通、高通、帶通和帶阻）已經(jīng)設(shè)計出來。

51、第三步離散化主要任務(wù)就是把模擬濾波器變換成數(shù)字濾波器，即把模擬濾波器的系數(shù)映射成數(shù)字濾波器的系統(tǒng)函數(shù)。數(shù)字濾波器的設(shè)計工作就全部完成。</p><p>　　實現(xiàn)系統(tǒng)傳遞函數(shù)s域至z域映射有脈沖響應(yīng)不變法和雙線性映射兩種方法。</p><p><b>　?、倜}沖響應(yīng)不變法</b></p><p>　　利用模擬濾波器成熟的理論和設(shè)計方法來設(shè)計IIR

52、數(shù)字低通濾波器是經(jīng)常用的方法。設(shè)計過程是：按照技術(shù)要求設(shè)計一個模擬低通濾波器，得到模擬低通濾波器的傳輸函數(shù)轉(zhuǎn)換成數(shù)字低通濾波器的系統(tǒng)函數(shù)H(z)。這樣設(shè)計的關(guān)鍵問題就是找這種轉(zhuǎn)換關(guān)系，將s平面上的轉(zhuǎn)換成z平面上的H(z)。為了保證轉(zhuǎn)換后的H(z)穩(wěn)定且滿足技術(shù)要求，對轉(zhuǎn)換關(guān)系提出兩點要求：</p><p>　　（1）因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。我們知道，模擬濾波器因果穩(wěn)定要求其傳輸函數(shù)的

53、極點全部位于s平面的左半平面；數(shù)字濾波器因果則要求H(z)的極點全部在單位圓內(nèi)。因此，轉(zhuǎn)換關(guān)系應(yīng)是s平面的左半平面映射z平面的單位圓內(nèi)部。</p><p>　?。?）數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻響，s平面的虛軸映射z平面的單位圓，相應(yīng)的頻率至間成線性關(guān)系。</p><p>　　將傳輸函數(shù)從s平面?zhèn)鲹Q到z平面的方法有多種，但工程上常用的是脈沖響應(yīng)不變法和雙線性變換法。我們先研究脈沖

54、響應(yīng)不變法。</p><p>　　設(shè)模擬濾波器的傳輸函數(shù)為，相應(yīng)的單位沖激響應(yīng)是，=LT[]</p><p>　　LT[.]代表拉氏變換，對進(jìn)行等間隔采樣，采樣間隔為T，得到，將h(n)= 作為數(shù)字濾波器的單位取樣響應(yīng)，那么數(shù)字濾波器的系統(tǒng)函數(shù)H(z)便是h(n)的Z變換。因此脈沖響應(yīng)不變法是一種時域上的轉(zhuǎn)換方法，它是h(n)在采樣點上等于。</p><p>　　設(shè)

55、模擬濾波器只有單階極點，且分母多項式的階次高于分子多項式的階次，將用部分分式表示：</p><p><b>　　（3.2）</b></p><p>　　式中為的單階極點。將逆拉氏變換得到：</p><p><b>　　（3.3）</b></p><p>　　式中u(t)是單位階躍函數(shù)。對進(jìn)行等間隔采

56、樣，采樣間隔為T,得到：</p><p><b>　?。?.4）</b></p><p>　　對上式進(jìn)行Z變換，得到數(shù)字濾波器的系統(tǒng)函數(shù)H(z)：</p><p><b>　　(3.5) </b></p><p>　　對比3.2、3.5式，的極點映射到z平面，其極點變成，系數(shù)不變化。下面我們分析從模

57、擬濾波器轉(zhuǎn)換到數(shù)字濾波器，s平面和z平面之間的映射關(guān)系，從而找到這種轉(zhuǎn)換方法的優(yōu)缺點。這里我們以采樣信號作為橋梁，推導(dǎo)其映射關(guān)系。</p><p>　　設(shè)的采樣信號用表示，</p><p><b>　　（3.6）</b></p><p>　　對進(jìn)行拉氏變換，得到：</p><p><b>　?。?.7）<

58、/b></p><p>　　式中是在采樣點t=n T時的幅度值，它與序列h(n)的幅度值相等，它與序列h(n)的幅度值相等，即h(n)= ，因此得到：</p><p>　　上式表示采樣信號的拉氏變換與相應(yīng)的序列的Z變換之間的映射關(guān)系可用下式表示：</p><p><b>　?。?.8） </b></p><p>　

59、　我們知道模擬信號的傅立葉變換和其采樣信號的傅立葉變換和其采樣信號的傅立葉變換之間的關(guān)系滿足式：</p><p><b>　?。?.9）</b></p><p><b>　　將代入上式，得</b></p><p><b>　?。?.10）</b></p><p>　　由 （3.

60、8）、（3.10）式得到：</p><p><b>　　（3.11）</b></p><p>　　上式表明將模擬信號的拉氏變換在s平面上沿虛軸按照周期延遲后，再按照（3.8）式映射關(guān)系，映射到z平面上，就得到H（z）。（3.8）式可稱為標(biāo)準(zhǔn)映射關(guān)系。下面進(jìn)一步分析這種映射關(guān)系。設(shè)</p><p>　　按照（3.8）式，得到：</p>

61、<p><b>　　因此得到：</b></p><p><b>　　（3.12）</b></p><p><b>　　那么 </b></p><p>　　上式關(guān)系式說明，s平面的虛軸（）映射z平面的單位圓（），s平面左平面（）映射z平面單位內(nèi)（），s平面右半平面映射ｚ平面單位圓外（）

62、。這說明如果因果穩(wěn)定，轉(zhuǎn)換后得到Ｈ（ｚ）仍是因果穩(wěn)定的。</p><p>　　另外，注意到z=是一個周期函數(shù)，可寫成</p><p>　　M為任意整數(shù)當(dāng)不變，模擬頻率變化2/T的整數(shù)倍時，映射值不變?；蛘哒f，將s平面沿著j軸分割成一條條寬為2/T的水平帶，每條水平面都按照前面分析的映射關(guān)系對應(yīng)著整個z平面。此時所在的s平面與所在的z平面的映射關(guān)系如圖所示。當(dāng)模擬頻率從變化到時，數(shù)字頻率則從

63、變化到，且按照（3.12）式，，即與之間成線性關(guān)系。</p><p>　　但是，從模擬信號到采樣信號，其拉氏變換要按照　　式。其為周期，沿虛軸方向進(jìn)行周期化。如果原模擬信號　　得頻帶不是限于之間，則會在得奇數(shù)倍附近產(chǎn)生頻率混疊，從而映射到ｚ平面上，在附近產(chǎn)生頻率混疊。脈沖相應(yīng)不變法的頻率混疊現(xiàn)象如圖所示。這種會使設(shè)計出的濾波器在附近的頻率特性，程度不同地偏離模擬濾波在　　　　附近的頻率特性，嚴(yán)重時使濾波器不滿足給

64、定的技術(shù)指標(biāo)。因此，希望設(shè)計的濾波器使帶限濾波器，如果不是帶限的，例如高通濾波器，帶阻濾波器，需要在高通帶阻濾波器之前加保護濾波器，濾除高于折疊頻率以上的頻帶，以避免產(chǎn)生頻率混疊現(xiàn)象。但這樣會增加系統(tǒng)的成本和復(fù)雜性，因此，高通與帶阻濾波器不適合用這種方法設(shè)計。</p><p>　　假設(shè)沒有頻率混疊現(xiàn)象，即滿足</p><p>　　將關(guān)系式代入，得到：</p><p>

65、;　　說明用脈沖響應(yīng)不變法設(shè)計的數(shù)字濾波器可以很好地重現(xiàn)原模擬濾波器的頻響。上式中，的幅度特性與采樣間隔成反比，這樣當(dāng)T較小時，就會有太高的增益。為避免這一現(xiàn)象，令那么,此時</p><p>　　一般的極點是一個復(fù)數(shù)，且以共軛成對的形式出現(xiàn)，在式中講一對復(fù)數(shù)共軛極點放在一起，形成一個二階基本節(jié)。如果模擬濾波器的二階基本節(jié)的形式為，極點為可以推導(dǎo)出相應(yīng)的數(shù)字濾波器二階基本節(jié)（只有實數(shù)乘法）的形式為如果模擬濾波器二階

66、基本節(jié)的形式為，極點為，則對應(yīng)的數(shù)字濾波器二階基本節(jié)的具體形式為利用以上這些變換關(guān)系，可以簡化設(shè)計。</p><p><b>　?、陔p線性變換法</b></p><p>　　這種變換方法，采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到±π/T之間，再用轉(zhuǎn)換到z平面上。設(shè)Ha（s），s=jΩ，經(jīng)過非線性頻率壓縮后用Ha（s1），=jΩ1表示，這里用正切

67、變換實現(xiàn)頻率壓縮：</p><p><b>　　（3.13）</b></p><p>　　式中T仍是采樣間隔，當(dāng)Ω1從-π/T經(jīng)過0變化到π/T時，Ω則由-∞經(jīng)過0變化到+∞，實現(xiàn)了s平面上整個虛軸完全壓縮到平面上虛軸的±π/T之間的轉(zhuǎn)換。這樣便有</p><p><b>　　（3.14）</b></p&g

68、t;<p>　　再通過 轉(zhuǎn)換到z平面上，得到：</p><p><b>　?。?.15）</b></p><p><b>　?。?.16）</b></p><p>　　式（3.15）或式（3.16）稱為雙線性變換。從s平面映射到平面，再從平面映射到z平面，其映射情況如圖3.1所示。由于從s平面到平面具有非線性

69、頻率壓縮的功能，因此不可能產(chǎn)生頻率混疊現(xiàn)象。另外，從平面轉(zhuǎn)換到z平面仍然采用標(biāo)準(zhǔn)轉(zhuǎn)換關(guān)系 ，平面的 ±π/T之間水平帶的左半部分映射z平面單位圓內(nèi)部，虛軸映射單位圓。這樣，Ha（s）因果穩(wěn)定，轉(zhuǎn)換成的H（z）也是因果穩(wěn)定的。下面分析模擬頻率Ω和數(shù)字頻率ω之間的關(guān)系。</p><p>　　令 ，并代入（3.15）式中，有</p><p><b>　　（3.17）</

70、b></p><p>　　上式說明，s平面上Ω與平面上的ω成非線性正切關(guān)系，如圖3.2所示。在ω=0附近接近線性關(guān)系；當(dāng)ω增加是，Ω增加得愈來愈快；當(dāng)ω趨近π時，Ω趨近于∞。正是因為這種非線性關(guān)系，消除了頻率混疊現(xiàn)象。</p><p>　　圖3.2 雙線性變換法的頻率變換</p><p>　　ω與Ω之間的非線性關(guān)系是雙線性變換法的缺點，直接影響數(shù)字濾波器頻響

71、逼真的模仿模擬濾波器的頻響，幅度特性和相位特性失真的情況如圖3.3所示。這種非線性影響的實質(zhì)問題是：如果Ω的刻度是均勻的，則影射到z平面ω的刻度不是均勻的，而是隨ω增加愈來愈密。</p><p>　　圖3.3 雙線性變換法幅度和相位特性的非線性映射</p><p>　　雙線性變換法可由簡單的代數(shù)公式（3.15）將Ha（s）直接轉(zhuǎn)換成H（z），這是該變換法的優(yōu)點。但當(dāng)階數(shù)稍高時，將H（z

72、）整理成需要的形式，也不是一件簡單的工作。</p><p>　　3.3 IIR濾波器直接設(shè)計</p><p>　　IIR數(shù)字濾波器的經(jīng)典設(shè)計法只限于幾種標(biāo)準(zhǔn)的低通、高通、帶通、帶阻濾波器，而對于具有形狀或多頻帶濾波器的設(shè)計是無能為力的。</p><p>　　如果所設(shè)計IIR濾波器幅頻特性比較復(fù)雜，可采用最小二乘法擬合給定的幅頻響應(yīng)，使設(shè)計的濾波器幅頻特性逼近期望的

73、頻率特性，這種方法稱為IIR濾波器的直接設(shè)計法。</p><p>　　MATLAB信號處理工具箱函數(shù)YULEWALK采用直接設(shè)計IIR數(shù)字濾波器。</p><p>　　函數(shù)YULEWALK采取下面步驟計算分子多項式：</p><p>　　計算與分子多項式相應(yīng)的幅值平方響應(yīng)的輔助式；</p><p>　　由輔助分子式和分母多項式計算完全的頻率響

74、應(yīng)；</p><p>　　計算濾波器的脈沖響應(yīng)；</p><p>　　采用最小二乘法擬合脈沖響應(yīng)最終求得濾波器的分子多項式系數(shù)。</p><p>　　函數(shù)YULEWALK允許我們自由定義濾波器的頻率向量f和幅值向量m，因此該函數(shù)可設(shè)計具有任意形狀的幅頻響應(yīng)的濾波器，包括多頻帶濾器。</p><p>　　第4章 DSP仿真系統(tǒng)</p

75、><p>　　4.1 對低通模擬和數(shù)字濾波器的仿真</p><p>　　4.1.1 模擬低通濾波器的仿真</p><p>　　第三章所討論的設(shè)計思想，在DSP對模擬濾波器進(jìn)行仿真的過程中依然適用，其具體步驟總結(jié)如下：</p><p>　　確定模擬濾波器的性能指標(biāo)，如截止頻率（對于低通和高通）或上、下邊界頻率,；波紋特性；帶阻衰減等。</p&

76、gt;<p><b>　　確定濾波器階數(shù)。</b></p><p>　　計模擬低通濾波原型濾波器。MATLAB信號處理工具箱的濾波器原型函數(shù)buttap,cheb1ap。</p><p>　　按頻率變換設(shè)計模擬濾波器（低通、高通、帶通、帶阻）。</p><p>　　DSP信號處理工具箱的頻率變換函數(shù)lp2lp,lp2hp,lp2b

77、p,lp2bs。但是，按照這種設(shè)計思想的編程較為麻煩。MATLAB信號處理工具箱還提供模擬濾波器的完全設(shè)計函數(shù)：butter,cheby1等。用戶只需調(diào)用一次設(shè)計函數(shù)就可自動完成全部設(shè)計過程，編程十分簡單。下面將以模擬低通Butterworth濾波器的設(shè)計為例，解釋這種設(shè)計方法。</p><p>　　設(shè)計指標(biāo)：通帶截止頻率=200π，阻帶截止頻率=300π，通帶衰減=1dB，阻帶衰減=16dB。仿真結(jié)果見圖4.1

78、</p><p>　　%Matlab program2.1</p><p>　　%Design a buttworth analog lowpass filter</p><p>　　ws=300*pi;</p><p>　　wp=200*pi;</p><p><b>　　Rp=1;</b>&l

79、t;/p><p><b>　　Rs=16;</b></p><p>　　%compute oder and cuttoff frequency</p><p>　　[N,Wn]=buttord(wp,ws,Rp,Rs,'s') </p><p>　　Fc=Wn/(2*pi) </p><p

80、>　　[b,a]=butter(N,Wn,'s'); </p><p><b>　　%output</b></p><p>　　w=linspace(1,3000,1000)*2*pi; </p><p>　　H=freqs(b,a,w); </p><p>　　magH=abs(H); </

81、p><p>　　phaH=unwrap(angle(H)); </p><p>　　plot(w/(2*pi),20*log10(magH)); </p><p>　　xlabel('Frequency(Hz)'); </p><p>　　ylabel('Magnidute(dB)'); </p>&

82、lt;p><b>　　grid on</b></p><p><b>　　N = 7</b></p><p>　　Wn =725.7292</p><p>　　Fc =115.5034</p><p>　　圖4.1 模擬低通濾波器</p><p>　　4.2.2 數(shù)字

83、低通濾波器的仿真</p><p>　　IIR數(shù)字濾波器經(jīng)典設(shè)計法的一般步驟是：</p><p>　　（1）根據(jù)給定的性能指標(biāo)和方法不同，首先對設(shè)計性能指標(biāo)中的頻率指標(biāo)進(jìn)行轉(zhuǎn)換，轉(zhuǎn)換后的頻率指標(biāo)作為模擬濾波器原型設(shè)計性能指標(biāo)。</p><p>　?。?）估計模擬低通濾波器最小階數(shù)和邊界頻率，利用MATLAB工具函數(shù)buttord、cheb1ord等。</p>

84、;<p>　?。?）設(shè)計模擬低通濾波器原型，利用MATLAB工具函數(shù)buttap、cheb1ap等。</p><p>　?。?）由模擬低通原型經(jīng)頻率變換獲得模擬濾波器（低通、高通、帶通、帶阻），利用MATLAB工具函數(shù)lp2lp、lp2hp、lp2bp、lp2bs。</p><p>　?。?）將模擬濾波器離散化獲得IIR數(shù)字濾波器，利用MATLAB工具函數(shù)bilinear。&

85、lt;/p><p>　　設(shè)計IIR濾波器時，給出的性能指標(biāo)通常分為數(shù)字指標(biāo)和模擬指標(biāo)兩種。</p><p>　　數(shù)字性能指標(biāo)給出通帶截止頻率，阻帶截止頻率，通帶衰減，阻帶衰減等。數(shù)字頻率和的取值范圍為0～π，單位：弧度，而MATLAB工具函數(shù)常采用標(biāo)準(zhǔn)化頻率，和的取值范圍為0～1。</p><p>　　模擬性能指標(biāo)給出通帶截止頻率，阻帶起始頻率，通帶衰減，阻帶衰減等。模

86、擬頻率和的單位均為弧度/秒。下面是一個利用雙線性變換設(shè)計一個Butterworth低通濾波器的示例：</p><p>　　技術(shù)指標(biāo)：通帶截止頻率=2π，阻帶截止頻率，通帶衰減小于3dB，阻帶衰減大于15dB,采樣頻率。仿真結(jié)果見圖4.2</p><p>　　Wp= 2000*2*pi; </p><p>　　Ws=3000*2*pi; </p><

87、;p><b>　　Rp=3; </b></p><p><b>　　Rs=15; </b></p><p>　　Fs=10000; </p><p><b>　　Ts=1/Fs; </b></p><p><b>　　Nn=128;</b></

88、p><p>　　[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s') </p><p>　　[b,a]=butter(N,Wn,'s'); </p><p>　　[z,p,k]=buttap(N); </p><p>　　[Bap,Aap]=zp2tf(z,p,k); </p><p&

89、gt;　　[b,a]=lp2lp(Bap,Aap,Wn); </p><p>　　[bz,az]=bilinear(b,a,Fs) </p><p>　　freqz(bz,az,Nn,Fs) </p><p><b>　　N =5</b></p><p>　　Wn =1.3387e+004</p><

90、p>　　bz =0.0171 0.0854 0.1708 0.1708 0.0854 0.0171</p><p><b>　　az =</b></p><p>　　1.0000 -1.2271 1.1622 -0.5176 0.1450 -0.0159</p><p>　　4.4.1

91、 模擬帶通濾波器的仿真</p><p>　　函數(shù)BUTTER用于Butterworth濾波器設(shè)計，調(diào)用格式：</p><p>　　[b,a]=butter(n,,’s’) </p><p>　　[b,a]=butter(n, ,’ftype’,’s’) </p><p>　　其中，n為濾波器階數(shù)；為濾波器截止頻率，‘s’為模擬濾波器，確省時

92、為數(shù)字濾波器。</p><p>　　‘ftype’濾波器類型：</p><p>　　‘high’為高通濾波器，截止頻率；</p><p>　　‘stop’為帶阻濾波器，=（）；</p><p>　　‘ftype’缺省時為低通或帶通濾波器。</p><p>　　以下設(shè)計一個Butterworth模擬帶通濾波器，設(shè)計指標(biāo)為

93、：逼近頻率1000—2000Hz,兩側(cè)過渡帶寬500Hz，通帶衰減1dB，阻帶衰減大于100dB。結(jié)果見圖4.5。</p><p>　　%Matlab program4.1</p><p>　　%Design a butterworth analog bandpass filter</p><p>　　Wp=[1000 2000]*2*pi; </p>

94、<p>　　Ws=[500 2500]*2*pi; </p><p><b>　　Rp=1; </b></p><p><b>　　Rs=100; </b></p><p>　　[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s') </p><p>　　[b,a

95、]=butter(N,Wn,'s'); </p><p><b>　　N =23</b></p><p>　　Wn =1.0e+004 *</p><p>　　0.6220 1.2695</p><p>　　4.4.2 數(shù)字帶通濾波器的設(shè)計</p><p>　　設(shè)計一個帶通c

96、hebyshev I 型數(shù)字濾波器，通帶為100Hz—200Hz，過渡帶寬均為50Hz，通帶衰減小于1dB，阻帶衰減30Hz，采樣頻率1000Hz。仿真結(jié)果見圖4.6</p><p>　　%Matlab program4.2</p><p>　　%Design a chebyshev I digital bandpass filter</p><p><b&g

97、t;　　Fs=1000; </b></p><p>　　ws=[100 200]*2/Fs; </p><p>　　wp=[50 250]*2/Fs; </p><p><b>　　Rp=1; </b></p><p><b>　　Rs=30; </b></p><p

98、><b>　　Nn=128; </b></p><p>　　%compute oder and cuttoff frequency</p><p>　　[N,Wn]=cheb1ord(wp,ws,Rp,Rs,'s')</p><p>　　[b,a]=cheby1(N,Rp,Wn); </p><p>

99、;<b>　　%output</b></p><p>　　freqz(b,a,Nn,Fs); </p><p><b>　　N =5</b></p><p>　　Wn =0.1600 0.5000</p><p>　　4.5 對帶阻模擬和數(shù)字濾波器的仿真</p><p>

100、;　　4.5.1 模擬帶阻濾波器的設(shè)計</p><p>　　這里設(shè)計一個chebyshev I 型模擬帶阻濾波器，設(shè)計指標(biāo)為：阻帶頻率1000Hz—2000Hz，兩側(cè)過渡帶寬500Hz，通帶衰減1dB，阻帶衰減大于50dB。仿真結(jié)果見圖4.7。</p><p>　　%Matlab program5.1</p><p>　　%Design a chebyshev I

101、analog bandstop filter</p><p>　　ws=[1000 2000]*2*pi; </p><p>　　wp=[500 2500]*2*pi; </p><p><b>　　Rp=1; </b></p><p><b>　　Rs=50; </b></p>&l

102、t;p>　　%compute oder and cuttoff frequency</p><p>　　[N,Wn]=cheb1ord(wp,ws,Rp,Rs,'s') </p><p>　　[b,a]=cheby1(N,Rp,Wn,'stop','s'); </p><p><b>　　%output

103、</b></p><p>　　w=linspace(1,3000,1000)*2*pi; </p><p>　　H=freqs(b,a,w); </p><p>　　magH=abs(H); </p><p>　　phaH=unwrap(angle(H)); </p><p>　　plot(w/(2*pi)

104、,20*log10(magH)); </p><p>　　xlabel('Frequency(Hz)'); </p><p>　　ylabel('Magnidute(dB)'); </p><p><b>　　grid on</b></p><p>　　4.5.2 數(shù)字帶阻濾波器的仿真&l

105、t;/p><p>　　設(shè)計一個Butterworth數(shù)字帶阻濾波器，設(shè)計指標(biāo)為：阻帶頻率為100Hz—200Hz,兩側(cè)過渡帶寬500Hz，通帶衰減1dB,阻帶衰減大于30dB。仿真結(jié)果見圖4.8</p><p><b>　　Fs=1000; </b></p><p>　　ws=[100 200]*2/Fs; </p><p>

106、;　　wp=[50 250]*2/Fs; </p><p><b>　　Rp=1; </b></p><p><b>　　Rs=30; </b></p><p><b>　　Nn=128; </b></p><p>　　[N,Wn]=buttord(wp,ws,Rp,Rs) &

107、lt;/p><p>　　[b,a]=butter(N,Wn,'stop'); </p><p>　　freqz(b,a,Nn,Fs); </p><p><b>　　N =7</b></p><p>　　Wn =0.1606 0.4720</p><p>　　同樣的技術(shù)指標(biāo)，設(shè)計

108、一個chebyshev I 型帶阻濾波器，仿真結(jié)果見圖4.9</p><p>　　%Matlab program5.2.1</p><p>　　%Design a chebyshev I digital bandstop filter</p><p><b>　　Fs=1000;</b></p><p>　　ws=[10

109、0 200]*2/Fs;</p><p>　　wp=[50 250]*2/Fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rs=30;</b></p><p><b>　　Nn=128;</b></p><p>　　%

110、compute oder and cuttoff frequency</p><p>　　[N,Wn]=cheb1ord(wp,ws,Rp,Rs)</p><p>　　[b,a]=cheby1(N,Rp,Wn,'stop');</p><p><b>　　%output</b></p><p>　　fre

111、qz(b,a,Nn,Fs);</p><p><b>　　N = 4</b></p><p>　　Wn =0.1476 0.5000</p><p>　　前面已經(jīng)介紹，巴特沃斯濾波器在Ω=0附近一段范圍內(nèi)是非常平直的，它以原點的最大平坦性來逼近理想低通濾波器。而切比雪夫等波紋濾波器是在通帶內(nèi)以最大誤差最小化對理想低通濾波器的最佳一致逼近，

112、通帶內(nèi)誤差分布是均勻的。對比Butterworth濾波器和Chebshev I 型濾波器的設(shè)計結(jié)果，可以看出，巴特沃斯濾波器通帶，阻帶單調(diào)下降，這種濾波器具有良好的相頻特性。切比雪夫濾波器通帶內(nèi)的相頻特性有相應(yīng)的起伏波動，即相位是非線性的，這給信號傳輸時帶來線性畸變，所以在要求群時延為常數(shù)時不宜采用這種濾波器。但是，在同樣的技術(shù)指標(biāo)情況下，Chebshev I 濾波器的最小階數(shù)N比Butterworth濾波器的要小，而且濾波幅度特性更明

113、顯，更有利與濾除低通頻率。</p><p><b>　　第5章 總結(jié)與展望</b></p><p><b>　　5.1 總 結(jié)</b></p><p>　　幾乎在所有的工程技術(shù)領(lǐng)域中都會涉及到信號處理問題 ，濾波器作為信號處理的重要組成部分，已發(fā)展的相當(dāng)成熟。本論文主要是針對無限脈沖響應(yīng)（IIR）數(shù)字濾波器的設(shè)計與仿真。

114、在以上五章的論述中，系統(tǒng)分析了無限沖擊響應(yīng)數(shù)字濾波器的網(wǎng)絡(luò)結(jié)構(gòu)，研究了模擬濾波器（包括巴特沃斯濾波器和切比雪夫濾波器）的設(shè)計原理和方法，并在此基礎(chǔ)上論述了利用雙線性變換法完成對無限脈沖響應(yīng)（IIR）數(shù)字濾波器（包括低通、高通、帶通、帶阻）的設(shè)計。最后，利用MATLAB對所設(shè)計的模擬、數(shù)字濾波器進(jìn)行了仿真。</p><p><b>　　5.2 展 望</b></p><p

115、>　　信號處理已有一段悠久而豐富的歷史了, ,這一領(lǐng)域總是得益于它的理論,應(yīng)用與實現(xiàn)信號處理系統(tǒng)的技術(shù)之間的緊密結(jié)合, 日見增長的應(yīng)用范圍和日益增長的高級算法的需求總是與現(xiàn)實信號處理系統(tǒng)的器件技術(shù)的快速發(fā)展齊頭并進(jìn).在許多方面都清楚的表明,信號處理的重要性和地位在迅速提高和擴大。八十年代以來,DSP芯片的復(fù)雜性和容量一直成指數(shù)的增長著,并且沒有放慢的跡象.隨著整片集成技術(shù)的迅速發(fā)展,價廉、超微型和低功耗的很復(fù)雜的數(shù)字信號處理系統(tǒng)也

116、將會實現(xiàn)。因此，數(shù)字信號處理的重要性無疑仍會與日俱增，而濾波器作為數(shù)字信號處理的重要組成部分也必將迎來革命性的變革。</p><p>　　雖然數(shù)字信號處理是一個不斷更新和飛速發(fā)展的領(lǐng)域，但是它的基礎(chǔ)已經(jīng)日臻完善。本設(shè)計的目的就是為了在數(shù)字信號處理濾波理論方面給出一種條理清晰的論述。</p><p><b>　　參考文獻(xiàn)</b></p><p>

117、　　【1】王世一.數(shù)字信號處理.北京理工大學(xué)出版社. 1997【2】Bernard Gold, Lawrence R. Rabiner. "Theory and Application of DigitalSignal Processing". Prentice-Hall. Inc., 1975【3】 Dimitris G. Manolakis, John G. Proakis. "Digit

118、al Signal Processing: Principles, Algorithms, and Applications". Macmillan Publishing Company. 1992</p><p>　　【4】丁玉美 高西全 數(shù)字信號處理 西安電子科技大學(xué)</p><p>　　【5】王金龍 沈良 無線通信系統(tǒng) 人民郵電出版社</p><p>