1、<p>　　IIR數字濾波器的計算機優(yōu)化與設計</p><p>　　[摘要] 數字濾波技術是數字信號處理的一個重要組成部分，濾波器的設計是信號處理的核心問題之一，本文介紹了IIR數字濾波器的特性以及應用，分析了IIR數字濾波器的基本原理、系統(tǒng)函數、差分方程、技術指標、基本結構、設計流程。分別采用脈沖響應不變法和雙線性變換法，設計出巴特沃斯、切比雪夫I型、切比雪夫2型以及橢圓濾波器，并對這些不同類型的濾波

2、器進行了比較，分析它們的濾波性能；同時還對比分析了不同階數情況下，同一類型低通數字濾波器的濾波性能，設計出滿足技術指標的最優(yōu)IIR數字濾波器。</p><p>　　[關鍵詞] 數字濾波，IIR濾波器， Matlab</p><p>　　Optimization and Designing of IIR Digital Filter</p><p>　　Abstrac

3、t: Digital filter technology is an important part of digital signal processing, and the design of filter is a core problems in signal processing. The author introduced the characteristics and application of the IIR digit

4、al filter, and the basic principle, system function, difference equation, technical index, basic structure, design process of IIR digital filter are analyzed in this paper. Using the impulse response and bilinear transfo

5、rmation method, the Butterworth, chebyshev 1, chebysh</p><p>　　Key words: digital filter; IIR filter; Matlab </p><p><b>　　目 錄</b></p><p><b>　　1引言1</b><

6、;/p><p>　　1.1數字濾波器的研究背景與意義1</p><p>　　1.2數字濾波器的應用現狀與發(fā)展趨勢1</p><p>　　1.3論文的主要框架2</p><p>　　2 IIR數字濾波器的概述3</p><p>　　2.1 IIR數字濾波器基本原理3</p><p>　　2.

7、2 IIR數字濾波器特點3</p><p>　　2.3 IIR數字濾波器的優(yōu)越性3</p><p>　　2.4數字濾波器的分類和技術指標4</p><p>　　2.5數字濾波器的設計方案5</p><p><b>　　2.6設計思路5</b></p><p>　　3 IIR濾波器的算法設

8、計6</p><p>　　3.1 基本結構6</p><p>　　3.2 技術指標和系統(tǒng)函數6</p><p>　　3.3 IIR濾波器的設計原理和思路7</p><p>　　3.3.1低通濾波器的設計過程7</p><p>　　3.3.2 高通、帶通和帶阻濾波器的設計過程7</p><

9、p>　　3.4 IIR數字濾波器的設計方法7</p><p>　　3.4.1脈沖響應不變法7</p><p>　　3.4.2 雙線性變換法8</p><p>　　3.5 IIR數字濾波器的計算機實現10</p><p>　　3.5.1 IIR數字濾波器的計算機設計步驟10</p><p>　　3.5.2

10、 MATLAB工具箱中巴特沃斯濾波器的設計函數11</p><p>　　4 IIR數字濾波器的Matlab實現12</p><p>　　4.1直接設計法設計巴特沃斯帶通數字濾波器12</p><p>　　4.2 不同階數的巴特沃斯低通濾波器的特性比較12</p><p>　　4.3參數相同時四種低通濾波器的特性比較13</p&

11、gt;<p>　　4.4當階數相同時四種類型模擬濾波器的比較14</p><p>　　4.5設計巴特沃斯模擬高通濾波器15</p><p>　　4.6 設計巴特沃斯模擬帶通和帶阻濾波器15</p><p>　　4.7 設計相同參數的橢圓和巴特沃斯低通模擬濾波器16</p><p>　　5脈沖響應不變法和雙線性變換法設計I

12、IR數字低通濾波器16</p><p>　　5.1脈沖響應不變法與雙線性變換法的比較17</p><p>　　5.2用雙線性變換法設計IIR數字低通濾波器17</p><p>　　5.3給不同濾波器中加入載波信號進行對比20</p><p>　　6濾波器的硬件設計21</p><p>　　6.1三階低通濾波器

13、的設計21</p><p>　　6.2舉例說明21</p><p><b>　　7 結語23</b></p><p><b>　　致謝23</b></p><p><b>　　參考文獻24</b></p><p><b>　　附錄2

14、5</b></p><p><b>　　1 引言</b></p><p>　　1.1 數字濾波器的研究背景與意義</p><p>　　當今，數字信號處理[1]DSP (Digtal Signal Processing)技術正飛速發(fā)展，它不但自成一門學科，更是以不同形式影響和滲透到其他學科。它與國民經濟息息相關，與國防建設緊密相連；它

15、影響或改變著我們的生產、生活方式，因此受到人們普遍的關注。</p><p>　　數字化、智能化和網絡化是當代信息技術發(fā)展的大趨勢，而數字化是智能化和網絡化的基礎，實際生活中遇到的信號多種多樣，例如廣播信號、電視信號、雷達信號、通信信號、導航信號、射電天文信號、生物醫(yī)學信號、控制信號、氣象信號、地震勘探信號、機械振動信號、遙感遙測信號等等。上述這些信號大部分是模擬信號，也有小部分是數字信號。模擬信號是自變量的連續(xù)函

16、數，自變量可以是一維的，也可以是二維或多維的。大多數情況下一維模擬信號的自變量是時間，經過時間上的離散化(采樣)和幅度上的離散化(量化)，這類模擬信號便成為一維數字信號。因此，數字信號實際上是用數字序列表示的信號，語音信號經采樣和量化后，得到的數字信號是一個一維離散時間序列；而圖像信號經采樣和量化后，得到的數字信號是一個二維離散空間序列。數字信號處理，就是用數值計算的方法對數字序列進行各種處理，把信號變換成符合需要的某種形式。例如，對數

17、字信號進行濾波以限制他的頻帶或濾除噪音和干擾，或將他們與其他信號進行分離；對信號進行頻譜分析或功率譜分析以了解信號的頻譜組成，進而對信號進行識別；對信號進行某種變換，使之更適合于傳輸，存儲和應用；對信號進行編碼以達到數據壓縮</p><p>　　數字濾波技術是數字信號分析、處理技術的重要分支。無論是信號的獲取、傳輸，還是信號的處理和交換都離不開濾波技術，它對信號安全可靠和有效靈活地傳輸是至關重要的。在所有的電子系

18、統(tǒng)中，使用最多技術最復雜的要算數字濾波器了。數字濾波器的優(yōu)劣直接決定產品的優(yōu)劣。</p><p>　　1.2 數字濾波器的應用現狀與發(fā)展趨勢</p><p>　　在信號處理過程中，所處理的信號往往混有噪音，從接收到的信號中消除或減弱噪音是信號傳輸和處理中十分重要的問題。根據有用信號和噪音的不同特性，提取有用信號的過程稱為濾波，實現濾波功能的系統(tǒng)稱為濾波器。在近代電信設備和各類控制系統(tǒng)中，數

19、字濾波器應用[2-5]極為廣泛，這里只列舉部分應用最成功的領域。</p><p><b>　　(1) 語音處理 </b></p><p>　　語音處理是最早應用數字濾波器的領域之一，也是最早推動數字信號處理理論發(fā)展的領域之一。該領域主要包括5個方面的內容：第一，語音信號分析，即對語音信號的波形特征、統(tǒng)計特性、模型參數等進行分析計算；第二，語音合成，即利用專用數字硬件或

20、在通用計算機上運行軟件來產生語音；第三，語音識別，即用專用硬件或計算機識別人講的話，或者識別說話的人；第四，語音增強，即從噪音或干擾中提取被掩蓋的語音信號。第五，語音編碼，主要用于語音數據壓縮，目前已經建立了一系列語音編碼的國際標準，大量用于通信和音頻處理。近年來，這5個方面都取得了不少研究成果，并且，在市場上已出現了一些相關的軟件和硬件產品，例如，盲人閱讀機、啞人語音合成器、口授打印機、語音應答機，各種會說話的儀器和玩具，以及通信和視

21、聽產品大量使用的音頻壓縮編碼技術。</p><p><b>　　(2) 圖像處理</b></p><p>　　數字濾波技術以成功地應用于靜止圖像和活動圖像的恢復和增強、數據壓縮、去噪音和干擾、圖像識別以及層析X射線攝影，還成功地應用于雷達、聲納、超聲波和紅外信號的可見圖像成像。</p><p><b>　　(3) 通信</b&g

22、t;</p><p>　　在現代通信技術領域內，幾乎沒有一個分支不受到數字濾波技術的影響。信源編碼、信道編碼、調制、多路復用、數據壓縮以及自適應信道均衡等，都廣泛地采用數字濾波器，特別是在數字通信、網絡通信、圖像通信、多媒體通信等應用中，離開了數字濾波器，幾乎是寸步難行。其中，被認為是通信技術未來發(fā)展方向的軟件無線電技術，更是以數字濾波技術為基礎。</p><p><b>　　(

23、4) 電視</b></p><p>　　數字電視取代模擬電視已是必然趨勢。高清晰度電視的普及指日可待，與之配套的視頻光盤技術已形成具有巨大市場的產業(yè)；可視電話和會議電視產品不斷更新?lián)Q代。視頻壓縮和音頻壓縮技術所取得的成就和標準化工作，促成了電視領域產業(yè)的蓬勃發(fā)展，而數字濾波器及其相關技術是視頻壓縮和音頻壓縮技術的重要基礎。</p><p><b>　　(5) 雷達&l

24、t;/b></p><p>　　雷達信號占有的頻帶非常寬，數據傳輸速率也非常高，因而壓縮數據量和降低數據傳輸速率是雷達信號數字處理面臨的首要問題。告訴數字器件的出現促進了雷達信號處理技術的進步。在現代雷達系統(tǒng)中，數字信號處理部分是不可缺少的，因為從信號的產生、濾波、加工到目標參數的估計和目標成像顯示都離不開數字濾波技術。雷達信號的數字濾波器是當今十分活躍的研究領域之一。</p><p&g

25、t;<b>　　(6) 聲納</b></p><p>　　聲納信號處理分為兩大類，即有源聲納信號處理和無源聲納信號處理，有源聲納系統(tǒng)涉及的許多理論和技術與雷達系統(tǒng)相同。例如，他們都要產生和發(fā)射脈沖式探測信號，他們的信號處理任務都主要是對微弱的目標回波進行檢測和分析，從而達到對目標進行探測、定位、跟蹤、導航、成像顯示等目的，他們要應用到的主要信號處理技術包括濾波、門限比較、譜估計等。</

26、p><p>　　(7) 生物醫(yī)學信號處理 </p><p>　　數字濾波器在醫(yī)學中的應用日益廣泛，如對腦電圖和心電圖的分析、層析X射線攝影的計算機輔助分析、胎兒心音的自適應檢測等。</p><p><b>　　(8) 音樂</b></p><p>　　數字濾波器為音樂領域開辟了一個新局面，在對音樂信號進行編輯、合成、以及在音

27、樂中加入交混回響、合聲等特殊效果方面，數字濾波技術都顯示出了強大的威力。數字濾波器還可用于作曲、錄音和播放，或對舊錄音帶的音質進行恢復等。</p><p><b>　　(9) 其他領域</b></p><p>　　數字濾波器的應用領域如此廣泛，以至于想完全列舉他們是根本不可能的，除了以上幾個領域外，還有很多其他的應用領域。例如，在軍事上被大量應用于導航、制導、電子對抗

28、、戰(zhàn)場偵察；在電力系統(tǒng)中被應用于能源分布規(guī)劃和自動檢測；在環(huán)境保護中被應用于對空氣污染和噪聲干擾的自動監(jiān)測，在經濟領域中被應用于股票市場預測和經濟效益分析等等。</p><p>　　1.3 論文的主要框架</p><p>　　本篇論文主要討論IIR濾波器的基本結構和特點，硬件設計有源三階低通濾波器。計算機設計不同參數的IIR數字濾波器，首先將要設計的數字濾波器的技術指標轉換成相應的模擬濾波

29、器的技術指標，接著按照轉換后的模擬濾波器的技術指標設計出相應的原型模擬濾波器，最后以一定的方法將設計出的模擬濾波器轉換為用戶所需的數字濾波器。用Matlab工具箱設計不同階數的巴特沃斯低通濾波器性能的對比，再用同種方法設計出巴特沃斯高通模擬濾波器，然后設計出巴特沃斯帶阻濾波器、帶通濾波器，設計出在同種參數下不同類型的數字濾波器的對比，最后用脈沖響應不變法和雙線性變化法設計出低通數字濾波器，找出最優(yōu)的設計方案。</p>&l

30、t;p>　　2 IIR數字濾波器的概述</p><p>　　2.1 IIR數字濾波器基本原理</p><p>　　實際應用中，多數情況利用數字濾波器來處理模擬信號。處理模擬信號的數字濾波器基本結構如圖2.1所示。</p><p>　　圖2.1中，輸入端接入一個低通濾波器H1(S)，其作用是對輸入信號的頻帶進行限制，以避免頻譜混疊；輸出端也接一個低通濾波器H2

31、(S)，以便將D/A變換輸出的模擬量良好地恢復成連續(xù)時間信號[6]。</p><p>　　設計IIR數字濾波器一般采用間接法（脈沖響應不變法和雙線性變換法），應用最廣泛的是雙線性變換法?；驹O計過程是：（1）將給定的數字濾波器的指標轉換成過渡模擬濾波器的指標；（2）設計過渡模擬濾波器；（3）將過渡模擬濾波器系統(tǒng)函數轉換成數字濾波器的系統(tǒng)函數。</p><p>　　IIR數字濾波器的設計一般

32、是利用目前已經很成熟的模擬濾波器的設計方法進行設計，通常采用模擬濾波器原型有巴特沃斯函數、切比雪夫函數、橢圓濾波器函數等[1]。</p><p>　　2.2 IIR數字濾波器特點</p><p>　　(1)IIR數字濾波器的系統(tǒng)函數可以寫成封閉函數的形式。</p><p>　　(2)IIR數字濾波器采用遞歸型結構，即結構上帶有反饋環(huán)路。IIR數字濾波器運算結構通常由

33、延時、乘以系數和相加等基本運算組成，可以組合成直接型、級聯(lián)型、并聯(lián)型三種結構形式，都具有反饋回路。</p><p>　　(3)IIR數字濾波器在設計上可以借助成熟的模擬濾波器的設計方法進行設計，通常采用模擬濾波器原型有巴特沃斯函數、切比雪夫函數、橢圓濾波器函數等，由現成的設計數據或圖表可查，其設計工作量比較小，對計算工具的要求不高，在設計一個IIR數字濾波器時，我們根據指標寫出模擬濾波器的公式，然后通過一定的變換

34、，將模擬濾波器的公式轉換成數字濾波器的公式。</p><p>　　(4)IIR數字濾波器的相位特性不好控制，對相位要求較高時，需加相位校準網絡[7]。</p><p>　　2.3 IIR數字濾波器的優(yōu)越性</p><p>　　從性能上來說，IIR數字濾波器傳輸函數的極點可位于單位圓內的任何地方，因此可用較低的階數獲得較高的選擇性，所用的存貯單元少，信號延遲小，所以經

35、濟而效率高。目前已廣泛應用于通信、語音、雷達、聲納、地震、生物醫(yī)學、遙控遙測、地質勘探、自動控制、圖像處理、航空航天、故障檢測、自動儀表等領域。數字濾波是數字信號處理的一個重要分支，目前，數字信號濾波器的設計圖像處理、數據壓縮等方面的應用取得了令人矚目的進展和成就。它具有可靠性好、精度高、靈活、體積小、重量輕等優(yōu)點，它是通過數字濾波器來實現的。鑒于此，數字濾波器的設計就顯得尤為重要。</p><p>　　2.4數

36、字濾波器的分類和技術指標</p><p>　　按照不同的分類方法，數字濾波器有許多種類，但總起來可以分成兩大類：經典濾波器和現代濾波器。經典濾波器的特點是其輸入信號中有用的頻率成分和希望濾除的頻率成分占有不同的頻帶，通過一個合適的選頻濾波器濾除干擾，得到純凈信號，達到濾波的目的。但是，如果信號和干擾的頻譜相互重疊，則經典濾波器不能有效地濾除干擾，最大限度地恢復信號，這時就需要現代濾波器，例如維納濾波器、卡爾曼濾波

37、器、自適應濾波器等最佳濾波器?，F代濾波器是根據隨機信號的一些統(tǒng)計特性，在某種最佳準則下，最大限度地抑制干擾，同時最大限度地回復信號，從而達到最佳濾波的目的[1]。</p><p>　　數字濾波器從頻率特性可以分為：低通濾波器、高通濾波器、帶通濾波器、帶阻濾波器。</p><p>　　數字濾波器從實現的網絡結構上分類，可分為：無限脈沖響應濾波器和有限脈沖響應濾波器。</p>&

38、lt;p>　　濾波器的技術要求主要包括四個方面，即：</p><p>　　1. 低通濾波器的截止頻率主要包括通帶截止頻率和阻帶下限截止頻率；</p><p>　　2. 高通濾波器的截止頻率主要包括通帶截止頻率和阻帶上限截止頻率；</p><p>　　3. 帶通濾波器的截止頻率主要包括通帶下限截止頻率，通帶上限截止頻率；</p><p>

39、　　4. 帶祖濾波器的截止頻率與帶通濾波器一樣，也主要包括通帶下限截止頻率，通帶上限截止頻率，下阻帶截止頻率，以及上阻帶截止頻率。</p><p>　　濾波器中帶通帶阻的容限的具體技術指標，往往由允許的最大衰減 及阻帶應達到的最小衰減給出。通帶及阻帶的衰減，分別定義為：</p><p><b>　　（2.1）</b></p><p><b

40、>　　（2.2）</b></p><p>　　式中均假定已被歸一化為1。</p><p>　　由于在數字濾波器設計中是用弧度表示的，而實際上給出的頻率要求往往是實際頻率，單位是，因此在數字濾波器的設計中還應給出采樣頻率。</p><p>　　2.5數字濾波器的設計方案</p><p>　　數字濾波器的設計方法有兩種：直接法和

41、間接法。</p><p>　　IIR數字濾波器和FIR數字濾波器都可以采用直接法進行設計，對于IIR數字濾波器主要分為兩種設計方法：時域直接設計法和頻域直接設計法。圖2.2為直接法設計數字濾波器的過程[7]。</p><p>　　圖2.2直接法設計數字濾波器框圖</p><p>　　間接法的思路是先設計一個合適模擬濾波器，然后通過脈沖響應不變法和雙線性變換法來變換成

42、滿足預定指標的數字濾波器，用間接法設計IIR數字濾波器可以利用現成的模擬濾波器設計公式，是目前手工設計階數較低的IIR數字濾波器的主要方法。圖2.3是間接法設計數字濾波器的過程[7]。</p><p>　　圖2.3間接法設計數字濾波器框圖</p><p>　　由圖2.2和圖2.3比較可以看出：間接法設計思路清晰，步驟詳盡，可參閱公式等。但是由于其計算繁瑣，手工計算大多只能用來進行簡單的低階

43、選頻濾波器的設計，而直接設計方法，可應用MATLAB軟件信號處理工具箱中所提供的一整套函數，通過原型變換法直接設計各種典型的濾波器。</p><p><b>　　2.6設計思路</b></p><p>　　數字濾波器的設計就是用一個因果穩(wěn)定的離散線性時不變系統(tǒng)（LTI）去逼近給定的性能指標要求。系數的求解在數學上是一個逼近問題，如果在s平面上去逼近，就得到模擬濾波器；

44、如果在z平面上去逼近，則是數字濾波器的設計。</p><p>　　無論是IIR數字濾波器設計，還是FIR數字濾波器設計，數字濾波器的大致由以下步驟完成。</p><p>　?。?）根據要求，確定濾波器的性能指標。</p><p>　　（2）用因果穩(wěn)定離散LTI系統(tǒng)的系統(tǒng)函數去逼近這一性能指標。</p><p>　　（3）用有限精度算法來實現這

45、個系統(tǒng)函數（其中包括濾波器結構選擇、有限字長效應的處理方法等）。</p><p>　?。?）利用計算機軟件或專用數字濾波器硬件的實際技術實現。</p><p>　　3 IIR濾波器的算法設計</p><p><b>　　3.1 基本結構</b></p><p>　　無限長單位脈沖響應網絡，存在輸出對輸入的反饋支路，IIR

46、系統(tǒng)的基本網絡結構有三種，直接型、級聯(lián)型、并聯(lián)型。實際應用中大量采用級聯(lián)和并聯(lián)兩種結構[1]。</p><p>　　3.2 技術指標和系統(tǒng)函數</p><p>　　設計濾波器時，總是先設計低通濾波器，再通過頻率變換將低通濾波器轉換成希望類型的濾波器。低通濾波器的技術指標是：</p><p>　?。和◣Ы刂诡l率（又稱通帶上限頻率）</p><p&g

47、t;<b>　　：阻帶下限截止頻率</b></p><p>　?。和◣г试S的最大衰減</p><p>　?。鹤鑾г试S的最小衰減（，的單位dB）</p><p><b>　?。和◣舷藿穷l率 </b></p><p><b>　?。鹤鑾孪藿穷l率</b></p>

48、<p>　　對于單調下降的幅度特性，和可表示成：</p><p><b>　　（3.1）</b></p><p><b>　?。?.2）</b></p><p><b>　　（3.3）</b></p><p><b>　?。?.4）</b><

49、;/p><p><b>　?。?.5）</b></p><p>　　如果能有，，和求出可求出所需要的Ha(s)。Ha(s)必須是因果穩(wěn)定的，因此極點必須落在s平面的左半平面，相應的Ha(-s)極點必然落在右半平面，Z域N階IIR數字濾波器的系統(tǒng)函數為公式：</p><p><b>　?。?.6）</b></p>

50、<p>　　IIR 數字濾波器可用一個n階差分方程表示：</p><p><b>　?。?.7）</b></p><p>　　公式(3.6)中、為濾波系數。當全為零時，該濾波器為FIR數字濾波器；當不全為零時，則為IIR濾波器。對照模擬濾波器的傳遞函數：</p><p><b>　?。?.8）</b></p

51、><p>　　數字濾波器與模擬濾波器的設計思路相仿，其設計實質也是尋找一組系數{b,a}，去逼近所要求的頻率響應，使其在性能上滿足預定的技術要求；不同的是模擬濾波器的設計是在S平面上，用數學逼近法去尋找近似的所需特性H(S)，而數字濾波器則是在Z平面尋找合適的H(Z)。IIR數字濾波器的單位響應是無限長的，而模擬濾波器一般都具有無限長的單位脈沖響應，因此與模擬濾波器相匹配。由于模擬濾波器的設計在理論上已十分成熟，因此

52、數字濾波器設計的關鍵是將H(S)→H(Z)，即：利用復值映射將模擬濾波器離散化。已經證明，脈沖響應不變法和雙線性變換法能較好地擔當此任，則在此基礎上，數字濾波器的設計就可首先歸結為模擬濾波器的設計了 [3]。</p><p>　　3.3 IIR濾波器的設計原理和思路</p><p>　　對于IIR數字濾波器的設計，通常采用模擬濾波器設計技術來實現IIR數字濾波器的設計。這是因為模擬濾波器的

53、設計已經有很多簡單而又有現成的設計公式可循，并且設計參數已經表格化，設計起來準確且方便，這種方法可使數字濾波器設計變得簡單有序、容易實現。通常模擬濾波器設計IIR數字濾波器有以下幾步完成：</p><p>　　(1)根據給定的數字濾波器指標要求，確定模擬濾波器的技術指標。</p><p>　　(2)按照模擬濾波器的性能指標要求，設計模擬濾波器的系統(tǒng)函數。模擬濾波器的種類和設計方法很多，其中

54、較為典型的、性能優(yōu)良的設計方法有：巴特沃斯（Butterworth）濾波器設計法、切比雪夫(Chebyshev)濾波器設計法和橢圓(Elliptic)濾波器設計法等。</p><p>　　(3)用適當的數字化方法（如脈沖響應不變法、雙線性變換法等），將模擬濾波器的系統(tǒng)函數轉換成數字濾波器的系統(tǒng)函數。</p><p>　　3.3.1低通濾波器的設計過程</p><p>

55、;　　(1)由得到象限對稱的s平面函數。</p><p>　　(2)將因式分解，得到零點和極點。</p><p>　　(3)由求出的的零點、極點及增益常數，可完全確定系統(tǒng)函數。</p><p>　　3.3.2 高通、帶通和帶阻濾波器的設計過程</p><p>　　高通、帶通和帶阻濾波器設計的一般過程是：</p><p>

56、;　　(1)通過頻率變換公式，現將需要設計的濾波器指標轉換為相應的低通濾波器的指標；</p><p>　　(2)設計相應的低通系統(tǒng)函數Q(p)；</p><p>　　(3)對Q(p)進行頻率變換，得到希望設計的濾波器系統(tǒng)函數H(S)。</p><p>　　3.4 IIR數字濾波器的設計方法 </p><p>　　3.4.1脈沖響應不變法<

57、;/p><p><b>　　1．變換原理</b></p><p>　　脈沖響應不變法的基本原理是將模擬濾波器的單位沖擊響應加以等間隔的抽樣（抽樣周期為T），使數字濾波器的單位沖激響應正好等于的抽樣乘值T ，即</p><p><b>　?。?.9）</b></p><p>　　沖激響應不變法的實現流程如

58、下：</p><p><b>　　(3.10)</b></p><p>　　抽樣序列的z變換與模擬信號的拉普拉斯變換的關系：</p><p><b>　?。?.11）</b></p><p>　　沖激響應不變法將模擬濾波器的s平面變換映射成數字濾波器的z平面，其變換關系式為：</p>

59、<p><b>　　(3.12)</b></p><p>　　2．模擬濾波器的數字化方法</p><p>　　從（3.11）式可知，模擬濾波器轉換成數字濾波器是通過以下公式得來的：</p><p>　　可將展成部分分式： （3.13）</p>&

60、lt;p>　　該模擬濾波器的單位脈沖響應為： （3.14）</p><p>　　其中u(t)是連續(xù)時間的單位脈沖響應。</p><p>　　沖激響應不變法要求h(n)等于的抽樣，由（3.10）和（3.15）可得：</p><p><b>　　（3.15）</b></p><p>　

61、　對h(n)求z變換可得數字濾波器的系統(tǒng)函數：</p><p><b>　?。?.16）</b></p><p>　　若時，（4.17）式可改寫為：</p><p>　?。?.17） 3．脈沖響應不變法的實現步驟</p><p>　　綜上所述，在實際設計過程

62、中，脈沖響應不變法設計IIR數字濾波器的步驟如下：</p><p>　　(1)確定數字濾波器的性能指標和。</p><p>　　(2)利用，將數字濾波器的性能指標變換成模擬濾波器的性能指標和。</p><p>　　(3)利用成熟的模擬濾波設計方法設計模擬濾波器，求出。</p><p>　　(4)若能寫成的形式，就可直接通過式將模擬濾波器轉變成

63、數字濾波器的系統(tǒng)函數H(z)。否則，需對進行拉氏逆變換求的，然后對按式進行抽樣，才能獲得數字濾波器的單位沖激響應h(n)。</p><p>　　3.4.2 雙線性變換法</p><p>　　脈沖響應不變法是使數字濾波器在時域上模仿模擬濾波器，其缺點是產生頻率響應的混疊失真。造成此缺陷的主要原因是，從s平面到z平面是多值得映射關系。為了辟免沖激響應不變法的缺點，對于模擬濾波器的數字化，人們又

64、提出了雙線性變換法。</p><p><b>　　1.變換原理</b></p><p>　　為了克服沖激響應不變法多只映射這一缺點，雙線性變換法則將整個s平面映射到平面中的一個帶寬為橫帶上（該帶域的下限頻率為，帶域的上限頻率為），然后通過變換關系把平面映射到整個z平面上。由于整個s平面與帶域為到的平面與z平面上的點也是一一對應的，由此可以消除多值映射所帶來的頻譜混疊現

65、象。</p><p>　　2.模擬濾波器數字化方法</p><p>　　首先，為了將s平面整個地變換到平面的橫帶到范圍，取頻率變換關系： （3.18）</p><p>　　其中，c是待定常數；Ω是s平面的縱坐標，即s平面的頻率變量；是平面的縱坐標，即平面的頻率變量。根據上式，當從到變化時，Ω將從變化，且一一對應。

66、將上式改寫為：</p><p><b>　　(3.19)</b></p><p>　　解析延拓到整個s平面和平面，令，，則得到s平面到平面變換的關系：</p><p>　　(3.20)然后，將平面通過以下標準變換關系映射到z平面：</p><p>　　Z= (

67、3.21)</p><p>　　將(3.22)代到（3.21）中，從而得到s平面與z平面的單值映射關系：</p><p><b>　?。?.22）</b></p><p><b>　?。?.23）</b></p><p>　　（3.22）和（3.23）則稱為雙線性變換公式。</p>&

68、lt;p>　　當c取值隨頻率不同而不同，分為兩種：（1）時，；（2）使數字濾波器的某一頻率與模擬濾波器的一個特定頻率嚴格相對應，則。</p><p>　　3.雙線性變換法的實現步驟</p><p>　?。?）確定數字濾波器的技術指標（通帶頻率、通帶最大衰減、阻帶頻率、阻最小衰減）。</p><p>　?。?）采用頻率轉換關系，將數字濾波器的技術指標轉換成模擬濾

69、波器的技術指標，即和轉換成和，對于通帶最大衰減和阻帶最小衰減。</p><p>　?。?）按照模擬濾波器的技術指標設計模擬濾波器，獲得模擬濾波器的系統(tǒng)函數。</p><p>　?。?）利用式（4.4.5），將s平面轉換到z平面，得到數字濾波器的系統(tǒng)函數，即</p><p><b>　?。?.24）</b></p><p>

70、;　　通常情況下，c一般取，T可任意取，一般情況下取1[8]。</p><p>　　3.5 IIR數字濾波器的計算機實現</p><p>　　3.5.1 IIR數字濾波器的計算機設計步驟</p><p>　　(1)指標歸一化處理。如果采用雙線性變換法，還需進行預畸變。</p><p>　　(2)確定最小階數N和頻率參數?？晒┻x用的階數擇函數有

71、：buttord，cheblord，cheb2ord，ellipord等。</p><p>　　(3)運用最小階數N設計模擬低通濾波器原型。模擬低通濾波器的創(chuàng)建函數有：buttap,cheblap, cheb2ap和ellipap，這些函數輸出的是零極點式形式，還要用zp2tf函數轉換成分子分母多項式形式。如果想根據最小階數直接設計模擬低通濾波器原型，可用butter，chebyl，cheby2，ellip等函數

72、，只是注意要將函數中的設為1。</p><p>　　(4)根據第2步的頻率參數，模擬低通濾波原型轉換模擬低通、高通、帶通、帶阻濾波器，可用函數分別是：lp21p,lp2hp,lp2bp,lp2bs。</p><p>　　(5)運用脈沖響應不變法或雙線性變法把模濾波器轉數字濾波器，調用的函數是impinvar和bilinear。脈沖響應不變法適用于采樣頻率大于4倍截止頻率的銳截止低通帶通濾波

73、器，而雙線性變換法適合于相位特性要求不高的各型濾波器。</p><p>　　具體設計框圖[8]如圖3.1所示：</p><p>　　圖3.1利用MATLAB設計IIR的框圖</p><p>　　由圖3.1可知，利用MATLAB設計IIR數字濾波器可分為以下幾步來實現：</p><p>　　按一定規(guī)則將數字濾波器的技術指標轉換為模擬低通濾波器技

74、術指標。</p><p>　　根據轉換后的技術指標，使用濾波器階數函數，確定濾波器的最小階數N和截止頻率Wc。</p><p>　　利用最小階數N產生模擬低通濾波原型。</p><p>　　利用截止頻率Wc把模擬濾波器原型轉換成模擬低通、高通、帶通或帶阻濾波器。</p><p>　　利用脈沖響應不變法或雙線性變換法把模擬濾波器轉換成數字濾波器

75、。</p><p>　　3.5.2 MATLAB工具箱中巴特沃斯濾波器的設計函數</p><p>　　[Z,P,K]=buttap(N)，用于計算N階巴特沃斯歸一化模擬低通原型濾波器系統(tǒng)函數的零、極點和增益因子，其中K表示濾波器增益。</p><p>　　[N,wc]=buttord(,,RP,As)，用于計算巴特沃斯濾波器的階數N和3dB截止頻率。</p&g

76、t;<p>　　[N,wc]=buttord(,,RP,As, ‘s’)，用于計算巴特沃斯濾波器的階數N和3dB截止頻率。</p><p>　　[B,A]=butter(N, ,‘ftype’)，用于計算N階巴特沃斯模擬濾波器的系統(tǒng)函數分子和分母多項式的系數向量B和A，有系統(tǒng)向量B和A可以寫出數字濾波器的系統(tǒng)函數： </p><p><b>　?。?.25） <

77、;/b></p><p>　　式（3.9）中，B(N)和A(N)分別為向量B和A的第N個元素。</p><p>　　[B,A]=butter（N,wc,‘ftype’,‘s’），用于計算N階巴特沃斯模擬濾波器的系統(tǒng)函數分子和分母多項式的系數向量B和A，有系統(tǒng)向量B和A可以寫出模擬濾波器的系統(tǒng)函數：</p><p><b>　?。?.26）</b

78、></p><p>　　用參數ftype來區(qū)分低通還是高通濾波器。當ftype=high，設計3dB截止頻率為的高通濾波器。缺省ftype時默認設計低通濾波器。ftype=stop時，設計3dB截止頻率為的帶阻濾波器，此時為二元向量[,],和分別為帶阻濾波器的通帶3dB下截止頻率和上截止頻率[2]。</p><p>　　注意：設計的帶通和帶阻濾波器的系統(tǒng)函數是2N階的。這是因為帶通濾

79、波器相當于N階低通濾波器與N階高通濾波器級聯(lián)[7-14]。</p><p>　　4 IIR數字濾波器的Matlab實現</p><p>　　4.1直接設計法設計巴特沃斯帶通數字濾波器</p><p>　　圖4.1 波特沃斯數字濾波器的幅頻響應</p><p>　　圖4.1中給出的是巴特沃斯型模擬帶通幅頻特性與數字帶通幅頻特性的比較，其中給定

80、的參數為：fp=[480Hz,520Hz]，fs=[450Hz,550Hz]，wp=2*pi*[480Hz,520Hz]，ws=2*pi*[450 Hz,550 Hz]， rp=3dB，rs=20dB。計算程序見附錄A。</p><p>　　幅頻特性對比由圖4.1可看出：模擬濾波器在500Hz時達到最高幅度，而數字濾波器在1.6Hz達到最高幅度；相頻特性對比由圖4.1可看出：模擬濾波器的通帶為500Hz， 數字濾

81、波器的通帶為1.6Hz。由此可得出結論：當給定數字濾波器指標時，采樣間隔T的取值對頻譜混疊程度影響很小，而對模擬濾波器的影響較大；數字濾波器的幅頻響應較好地逼近了模擬濾波器的幅度響應，但相位響應卻出現了較大失真。</p><p>　　4.2 不同階數的巴特沃斯低通濾波器的特性比較</p><p>　　圖4.2給出的是5階、4階、3階、2階巴特沃斯低通濾波器的幅頻特性曲線，給定參數為：通帶截

82、止頻率fp=5khz，fp1=4khz，fp2=3khz，fp3=1khz，通帶最大衰減ap=2dB，阻帶截止頻率fs=12khz，阻帶最小衰減as=30dB。計算程序見附錄B。</p><p>　　將B和A代入公式3.9寫出系統(tǒng)函數為：</p><p>　　濾波器的階數N選取了5階、4階、3階、2階，由圖4.2可以看出N值越大，通帶和阻帶的近似就越好，過渡帶的特性越陡，響應曲線在通帶內越

83、平緩，在阻帶內衰減的速度越大，實際應用中，根據需要選擇合適的濾波器階數。幅度平方函數特性如圖4.2所示。這種函數具有以下特點：通帶內具有最大平坦幅度特性，在正頻率范圍內 ，隨頻率升高而單調下降；階數越高，過渡帶越窄；沒有零點。</p><p>　　圖4.2 5階、4階、3階、2階低通濾波器的幅頻特性 </p><p>　　4.3參數相同時四種低通濾波器的特性比較</p>

84、<p>　　給定相同的濾波器參數：wp=20*2*pi，ws=50*2*pi，Rp=1，Rs=30，Fs=200，Ts=1/Fs，分別利用巴特沃斯濾波器、切比雪夫I型濾波器、切比雪夫II型濾波器、橢圓濾波器進行低通濾波器的設計。設計程序見附錄C。</p><p>　　圖4.3 相同參數下不同濾波器的低通濾波器的幅頻特性</p><p>　　調用MATLAB中相應的函數，運行

85、程序得出圖4.3：</p><p>　　由圖4.3運行結果很容易驗證，當給定相同濾波器參數時，通帶最大衰減ap和阻帶最小衰減as參數相同時，巴特沃斯濾波器具有單調下降的幅頻特性，過渡帶最寬。優(yōu)點是它的通帶最平坦，趨向阻帶時衰減單調增大，缺點是從通帶到阻帶的過渡帶最寬，對于帶外干擾信號的衰減作用最弱，過渡帶不夠陡峭。</p><p>　　兩種類型的切比雪夫濾波器相同點是：它們的過渡帶比巴特沃

86、斯濾波器的過渡帶窄，但比橢圓濾波器的過渡帶寬；切比雪夫1型和II型濾波器都具有等波紋幅度特性。不同之處是：切比雪夫I型濾波器的過渡帶比切比雪夫II型濾波器的過渡帶寬；切比雪夫I型濾波器的通帶幅頻響應幾乎與巴特沃思濾波器相同，通帶是等波紋特性；切比雪夫II型濾波器的通帶頻率比其他濾波器延遲大約15Hz左右，阻帶是等波紋幅頻特性。</p><p>　　橢圓濾波器的過渡帶最窄，通帶和阻帶均是等波紋幅頻特性。</p

87、><p>　　4.4當階數相同時四種類型模擬濾波器的比較</p><p>　　圖4.4 同階數不同類型的濾波器的對比圖</p><p>　　圖4.4是巴特沃斯濾波器（左上）和同階第一類切比雪夫濾波器（右上）、第二類切比雪夫濾波器（左下）、橢圓函數濾波器（右下）的頻率響應圖。</p><p>　　調用MATLAB濾波器設計函數，由4.4圖很容易驗

88、證，當階數相同時，對相同的通帶最大衰減ap和阻帶最小衰減as，巴特沃斯濾波器具有單調下降的幅頻特性，過渡帶最寬。它的通帶內沒有波紋，在靠近零頻處，有最平坦通帶，趨向阻帶時衰減單調增大，缺點是從通帶到阻帶的過渡帶最寬，對于帶外干擾信號的衰減作用最弱，過渡帶不夠陡峭。</p><p>　　兩種類型的切比雪夫濾波器的過渡帶寬度相等，比巴特沃斯濾波器的過渡帶窄，但比橢圓濾波器的過渡帶寬。切比雪夫1型濾波器在通帶具有等波紋

89、幅度特性，過渡帶和阻帶是單調下降的幅頻特性。切比雪夫2型濾波器的通帶幅頻響應幾乎與巴特沃思濾波器相同，阻帶是等波紋特性。</p><p>　　橢圓濾波器的過渡帶最窄，通帶和阻帶均是等波紋幅頻特性。</p><p>　　由上述比較可見，四種濾波器各具特點。工程中實際中選擇哪種濾波器取決于對濾波器階數和相位特性的具體要求。例如，在滿足所給相同指標的條件下希望濾波器階數最低時，就應當選擇橢圓濾波

90、器。</p><p>　　4.5設計巴特沃斯模擬高通濾波器</p><p>　　圖4.5 Q(P)損耗函數曲線圖 圖4.6 由低通轉換高通后的濾波器 圖4.7 直接設計的高通巴特沃斯濾波器</p><p>　　圖4.5給出的是5階巴特沃斯高通濾波器的幅頻特性曲線，計算中已知參數為通帶邊界頻率為4KHz,阻帶邊界頻率為1 KHz，通帶最大衰減為0.1dB

91、，阻帶最小衰減為40dB。先設計低通濾波器，由低通濾波器轉換為高通濾波器，計算程序見附錄D，圖4.6是轉化為低通濾波器的幅頻曲線圖。圖4.7是調用函數buttord和butter直接設計巴特沃斯高通濾波器，計算程序見附錄E。</p><p>　　4.6 設計巴特沃斯模擬帶通和帶阻濾波器</p><p>　　圖4.8 5階巴特沃斯帶通濾波器損耗函數曲線圖 圖4.9

92、 5階巴特沃斯帶阻濾波器幅頻圖</p><p>　　圖4.8的運行結果：</p><p><b>　　N = 5</b></p><p>　　圖4.9的運行結果：</p><p><b>　　N= 5</b></p><p>　　圖4.8給出的是5階巴特沃斯模擬帶通濾波器的損

93、耗函數曲線圖，計算中已知參數為通帶上、下邊界頻率分別為4KHZ和7KHZ，阻帶上、下邊界頻率分別為2KHZ和9KHZ，通帶最大衰減為1dB,阻帶最小衰減為20dB。計算程序見附錄F。圖4.9給出的是5階巴特沃斯模擬帶阻濾波器的幅頻函數圖，計算中已知參數是阻帶上、下邊界頻率分別為4KHZ和7KHZ，通帶上、下邊界頻率分別為2KHZ和9KHZ，通帶最大衰減為1dB，阻帶最小衰減為20dB，計算程序見附錄G。</p><p

94、>　　4.7 設計相同參數的橢圓和巴特沃斯低通模擬濾波器</p><p>　　圖4.10 相同參數不同類型的模擬擬低通濾波器損耗圖</p><p>　　運行結果：Nb = 7， Ne = 4</p><p>　　圖4.10給出的是7階巴特沃斯低通濾波器和4階橢圓模擬低通濾波器，計算中已知參數通帶截止頻率fp=3khz，通帶最大衰減ap=0.1dB，阻帶截止頻

95、率fs=12khz，阻帶最小衰減as=60dB。計算中程序見附錄H。由圖4.11可知，同種參數用不同類型設計的濾波器，所需的階數是不同的，巴特沃斯需要更高的階數（N=7），這就延長了計算時間。而橢圓濾波器的階數最低（N=4），階數差別較大。所以，就滿足濾波器幅頻響應指標而言，橢圓濾波器的性能價格比最高，應用將會更廣泛。</p><p>　　5脈沖響應不變法和雙線性變換法設計IIR數字低通濾波器</p>

96、<p>　　5.1 脈沖響應不變法設計數字低通濾波器</p><p>　　5.2 雙線性變換法設計數字低通濾波器</p><p><b>　　圖5.1運行結果：</b></p><p><b>　　N = 4</b></p><p><b>　　圖5.2運行結果：

97、</b></p><p><b>　　N = 2</b></p><p>　　圖5.1為脈沖響應不變法設計數字低通濾波器的幅頻特性曲線，指標參數如下：fp=2100，fs=8000，Fs=20000，Rp=0.5，Rs=30，T=1/Fs，計算中程序見附錄I。圖5.2為雙線性變換法設計低通數字濾波器的幅頻函數圖，已知參數頻率低于2100Hz時，容許

98、幅度誤差在0.5dB以內；在頻率8000Hz到10000Hz之間的阻帶衰減大于30dB。指定模擬濾波器采用巴特沃斯低通濾波器。計算程序見附錄J。通過運行結果得出以下結論：當給定相同參數來設計數字濾波器時，沖擊響應不變法的階數N=4，幅頻和相頻特性的過渡帶頻率為6800Hz；雙線性變換法的階數N=2，幅頻相頻特性過渡帶頻率為5500Hz。通過比較可知，對于給定相同參數來設計濾波器的方法中，雙線性變換法在硬件實現上性價比高，同時濾波效果好。

99、</p><p>　　5.1脈沖響應不變法與雙線性變換法的比較</p><p>　　脈沖響應不變法的優(yōu)點是：模擬頻率和數字頻率之間的轉換關系是線性的，并保持了模擬濾波器的時域瞬態(tài)特性。缺點：當模擬濾波器頻率響應不是嚴格限帶時，則用沖擊響應不變法設計出的數字濾波器在頻域出現混疊現象。脈沖響應不變法會造成頻響混疊，不宜用來設計高通、帶阻濾波器，適用于基本上是限帶的濾波器，如低通和帶通濾波器。設

100、計性能要求又較高時，則不宜使用這種方法。</p><p>　　雙線性變換法克服了頻率響應的混疊響應的混疊現象，但頻率變換關系產生了非線性。對具有分段常數幅頻特性的選頻濾波器來說，頻率非線性失真問題可以用預畸變方法解決。并且雙線性變換法同脈沖響應不變法相比，它具有計算簡單和易于實現的特點。因此，實際工作中廣泛采用雙線性變換法來設計IIR數字濾波器。</p><p>　　5.2用雙線性變換法設

101、計IIR數字低通濾波器</p><p><b>　　1.實現步驟： </b></p><p>　?。?）確定數字低通濾波器的技術指標：通帶截止頻率ωp、通帶衰減p、阻帶截止頻率ωs、阻帶衰減s。</p><p>　?。?）將數字低通濾波器的技術指標轉換成模擬低通濾波器的技術指標。如果采用脈沖響不變法，邊界頻率的轉換關系為。</p>

102、<p>　　如果采用雙線性變換法，邊界頻率的轉換關系為。</p><p>　　（3）按照模擬低通濾波器的技術指標設計模擬低通濾波器。</p><p>　?。?）利用雙線性變換法將模擬濾波器Ha(s)，從s平面轉換到z平面，得到數字低通濾波器系統(tǒng)函數H(z)。</p><p>　?。?）數字低通技術指標為：ωp=0.4πrad p=1dB ωs=0.

103、5πrad s=40dbB</p><p>　　（6）模擬低通的技術指標為：</p><p>　　歸一化截止角頻率為：wp=2pi*Fs/Ft， ws=2pi*Fs/Ft</p><p>　?。?）利用模擬切比雪夫濾波器設計數字濾波器。通帶截止頻率為：wp=0.4*pi； 阻帶截止頻率為：ws=0.5*pi；通帶最大衰減為：Rp=1；阻帶最大衰減為：As=15

104、；設定周期為1s；模擬低通濾波器的生成：[b,a]=cheby1(n,1,Wn,'low','s')； 滿足設計指標的最小階數n和截止頻率為：Wn[n,Wn]=cheb1ord(OmegaP,OmegaS,1,40,'s')。</p><p>　　上面的設計步驟可簡化為圖5.3：</p><p>　　圖5.3 濾波器設計流程方框圖</

105、p><p><b>　　2.舉例說明</b></p><p>　?。?）技術指標通帶邊界頻率為400Hz，阻帶邊界頻率為500Hz，通帶最大衰減1dB，阻帶最小衰減40dB，抽樣頻率為2000Hz。希望采用雙線性變換法設計切比雪夫I型數字濾波器。程序見附錄K。</p><p>　　圖5.4 切比雪夫1型低通濾波器及其傳輸零極點圖</p>

106、;<p>　　由圖5.4可知：極點全在單位元內部，因而該濾波器是穩(wěn)定的。其中濾波器傳輸函數的零極點對系統(tǒng)頻率響應的影響是：在原點（z=0）處的極點或零點至單位圓的距離大小不變，其值為，故對幅度影響不起作用；單位圓附近的零極點將對系統(tǒng)的幅度響應的位置和深度有明顯的影響；單位圓內且靠近單位圓附近的極點將對系統(tǒng)幅度響應凸峰的位置和峰度有明顯的影響。因此，適當的控制系統(tǒng)函數的極點、零點分布、就能改變數字濾波器的頻率響應特性，達到預

107、期的要求。</p><p>　　（2）當給（1）中加入信號，其中300Hz，600Hz，求濾波器的輸出，并給出的圖形。</p><p>　　圖5.5 濾波器輸入x、x1、x2與輸出y的波形</p><p>　　由5.5圖可知，輸出信號y的波形與給出的x1信號和x2信號的頻率有關，其中x信號是x1信號和x2 信號的疊加，x幅頻特性受x1和x2信號的影響，輸出信號y和

108、輸入信號x波形相同，頻域延遲了5Hz左右，由此可知，數字濾波器的作用就是通過一定運算關系改變輸入信號所含頻率成分的比例或者濾除某些頻率成分。</p><p>　　5.3給不同濾波器中加入載波信號進行對比</p><p>　　圖5.6 切比雪夫1型數字帶通濾波器輸入一個方波信號圖</p><p>　　圖5.7 巴特沃斯型數字帶通濾波器輸入一個方波信號</p&

109、gt;<p>　　圖5.6和圖5.7進行觀察對比可得：當給定指標參數為fp=[480,520]，fs=[450,550]，rp=3，rs=20，wp=2*pi*[480,520]， ws=2*pi*[450,550]，方波頻率為f=100Hz。用切比雪夫1和巴特沃斯設計的帶通濾波器來過濾該方波信號，通過觀察輸出函數波形，可以看出方波信號受到了濾波器的影響；同時將方波進行卷積作為輸入信號，觀察輸出信號可知，信號先增大最后逐漸

110、減小直到衰減到0dB，巴特沃斯濾波器輸出卷積波形時在150Hz以后的波動比較小，而切比雪夫1新數字濾波器在150Hz以后的波形還有所抖動。由此可知，巴特沃斯濾波器的通帶與阻帶的波紋比較平滑。在MATLAB中實現程序時，輸入信號頻率的改變會對相應數值的計算時間和波形圖有一定的影響，當頻率越高時，計算時間越長，仿真圖形越不逼真。所以得出以下結論，設置合理的輸入頻率會使結果更加清楚，同時節(jié)約了計算時間。計算中程序見附錄L和程序M。</p

111、><p><b>　　6濾波器的硬件設計</b></p><p>　　6.1三階低通濾波器的設計</p><p>　　在高于截止頻率f0時，其增益以-60dB/10倍頻程的速率下降。為了簡化設計，令R1=R2=R3，并引入計算電容的常數，如表6.1所示：</p><p>　　表6.1 階低通濾波器計算電容常數（k4、k5、k

112、6）</p><p>　　基本關系式為 AVD>>C2/(2C3),f0處AVD>100 (6.1)</p><p>　　R1=R2=R3=R (6.2)</p><p>　　m4=k/,m5=k/,m6=k

113、/ (6.3)</p><p>　　C1= m4/R,C2= m5/R,C3= m6/R (6.4)</p><p>　　式中，AVD為集成運放開環(huán)放大倍數；m4,m5,m6是截止頻率為f0時的計算系數[15-16]。</p><p><b>　　6.2舉例說明</b&g