課程設(shè)計---利用fdatool設(shè)計iir數(shù)字濾波器

1、<p>　　電子電路CAD及仿真技術(shù)專題訓(xùn)練說明書</p><p>　　題 目：利用FDATool設(shè)計IIR數(shù)字濾波器 </p><p>　　系 部： </p><p>　　專 業(yè)： </p><p>　　班 級： </p>

2、;<p>　　學(xué)生姓名: 學(xué) 號: </p><p>　　指導(dǎo)教師: </p><p>　　2012年 05 月24 日</p><p><b>　　目 錄</b></p><p><b>　　目 錄2</b&

3、gt;</p><p>　　1 實訓(xùn)任務(wù)與要求1</p><p>　　2 實訓(xùn)方法與內(nèi)容1</p><p>　　2.1. 數(shù)字濾波器簡介1</p><p>　　2.1.1 IIR濾波器的基本結(jié)構(gòu)1</p><p>　　2.1.2 數(shù)字濾波器的設(shè)計2</p><p>　　2.2 F

4、DATooL簡介2</p><p>　　2.3 利用FDATool設(shè)計IIR數(shù)字濾波器3</p><p>　　2.3.1 FDATooL程序界面3</p><p>　　2.3.2 濾波器各項指標(biāo)的設(shè)置4</p><p>　　2.3.3 濾波器的幅頻響應(yīng)4</p><p>　　2.3.4 濾波器的相頻

5、響應(yīng)5</p><p>　　2.3.5 濾波器的脈沖相應(yīng)5</p><p>　　2.3.6 濾波器的零極點分布6</p><p>　　2.3.7 濾波器系數(shù)6</p><p>　　2.3.8 濾波器的系統(tǒng)函數(shù)7</p><p><b>　　3 結(jié)論7</b></p>

6、;<p><b>　　4 參考文獻8</b></p><p><b>　　1 實訓(xùn)任務(wù)與要求</b></p><p>　　利用FDATool設(shè)計IIR數(shù)字濾波器。</p><p>　　輸出內(nèi)容包括濾波器指標(biāo)、幅頻響應(yīng)、相頻響應(yīng)、脈沖響應(yīng)、零極點分布、系數(shù)，寫出相應(yīng)的系統(tǒng)函數(shù)。</p><

7、;p>　　通過實驗加深對數(shù)字濾波器各項指標(biāo)的理解。</p><p><b>　　2 實訓(xùn)方法與內(nèi)容</b></p><p>　　2.1. 數(shù)字濾波器簡介</p><p>　　作為線形時不變系統(tǒng)的數(shù)字濾波器可以用系統(tǒng)函數(shù)來表示，而實現(xiàn)一個系統(tǒng)函數(shù)表達式所表示的系統(tǒng)可以用兩種方法：一種方法是采用計算機軟件實現(xiàn)；另一種方法是用加法器、乘法器、和

8、延遲器等元件設(shè)計出專用的數(shù)字硬件系統(tǒng)，即硬件實現(xiàn)。不論軟件實現(xiàn)還是硬件實現(xiàn)，在濾波器設(shè)計過程中，由同一系統(tǒng)函數(shù)可以構(gòu)成很多不同的運算結(jié)構(gòu)。對于無限精度的系數(shù)和變量，不同結(jié)構(gòu)可能是等效的，與其輸入和輸出特性無關(guān)；但是在系數(shù)和變量精度是有限的情況下，不同運算結(jié)構(gòu)的性能就有很大的差異。</p><p>　　2.1.1 IIR濾波器的基本結(jié)構(gòu)</p><p>　　一個數(shù)字濾波器可以用系統(tǒng)函數(shù)表示

9、為： </p><p>　　(2-1) </p><p>　　由這樣的系統(tǒng)函數(shù)可以得到表示系統(tǒng)輸入與輸出關(guān)系的常系數(shù)線形差分程為： </p><p><b>　　(2-2) </b></p><p>　　可見數(shù)字濾波器的功能就是把輸入序列x(n)通過一定的運算變換成輸出序列y(n)。不同的運算處

10、理方法決定了濾波器實現(xiàn)結(jié)構(gòu)的不同。無限沖激響應(yīng)濾波器的單位抽樣響應(yīng)h(n)是無限長的，其差分方程如(2-2)式所示，是遞歸式的，即結(jié)構(gòu)上存在著輸出信號到輸入信號的反饋，其系統(tǒng)函數(shù)具有(2-1)式的形式，因此在z平面的有限區(qū)間(0<︱z︱<∞)有極點存在。</p><p>　　前面已經(jīng)說明，對于一個給定的線形時不變系統(tǒng)的系統(tǒng)函數(shù)，有著各種不同的等效差分方程或網(wǎng)絡(luò)結(jié)構(gòu)。由于乘法是一種耗時運算，而每個延遲單

11、元都要有一個存儲寄存器，因此采用最少常熟乘法器和最少延遲支路的網(wǎng)絡(luò)結(jié)構(gòu)是通常的選擇，以便提高運算速度和減少存儲器。然而，當(dāng)需要考慮有限寄存器長度的影響時，往往也采用并非最少乘法器和延遲單元的結(jié)構(gòu)。IIR濾波器實現(xiàn)的基本結(jié)構(gòu)有：</p><p>　　(1)IIR濾波器的直接型結(jié)構(gòu)</p><p>　　(2)IIR濾波器的級聯(lián)型結(jié)構(gòu)；</p><p>　　(3)IIR濾

12、波器的并聯(lián)型結(jié)構(gòu)。</p><p>　　a、直接型 b、并聯(lián)型</p><p><b>　　c、串聯(lián)型</b></p><p>　　2.1.2 數(shù)字濾波器的設(shè)計</p><p>　　數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時域特性，可分為兩種，即無限長沖激響應(yīng)(I

13、IR)濾波器和有限長沖激響應(yīng)(FIR)濾波器。IIR濾波器的特征是，具有無限持續(xù)時間沖激響應(yīng)。這種濾波器一般需要用遞歸模型來實現(xiàn)，因而有時也稱之為遞歸濾波器。FIR濾波器的沖激響應(yīng)只能延續(xù)一定時間，在工程實際中可以采用遞歸的方式實現(xiàn)，也可以采用非遞歸的方式實現(xiàn)。數(shù)字濾波器的設(shè)計方法有多種，如雙線性變換法、窗函數(shù)設(shè)計法、插值逼近法和Chebyshev逼近法等等。隨著MATLAB軟件尤其是MATLAB的信號處理工作箱的不斷完善，不僅數(shù)字濾波

14、器的計算機輔助設(shè)計有了可能，而且還可以使設(shè)計達到最優(yōu)。本次試驗中利用分析設(shè)計工具箱FDATOOL設(shè)計IIR數(shù)字濾波器</p><p>　　2.2 FDATooL簡介</p><p>　　FDATool（Filter Design & Analysis Tool）是MATLAB信號處理工具箱里專用的濾波器設(shè)計分析工具，MATLAB6.0以上的版本還專門增加了濾波器設(shè)計工具箱（Fil

15、ter Design Toolbox）。FDATool可以設(shè)計幾乎所有的基本的常規(guī)濾波器，包括FIR和IIR的各種設(shè)計方法。它操作簡單，方便靈活。</p><p>　　FDATooL(Filter Design & Analysis Tool)是MATLAB信號處理工具箱里專用的濾波器設(shè)計分析工具，MATLAB 6.0以上的版本還專門增加了濾波器設(shè)計工具箱(Filter Design Toolbox) 。

16、FDATooL可以設(shè)計幾乎所有的常規(guī)濾波器，包括FIR和IIR的各種設(shè)計方法[3]。它操作簡單，方便靈活。</p><p>　　FDATooL界面總共分兩大部分，一部分是Design Filter，在界面的下半部，用來設(shè)置濾波器的設(shè)計參數(shù)；另一部分則是特性區(qū)，在界面的上半部分，用來顯示濾波器的各種特性。</p><p>　　Design Filter部分主要分為:Filter Type(濾

17、波器類型) 選項，包括Lowpass(低通) 、Highpass(高通) 、Bandpass (帶通) 、Bandstop(帶阻)和特殊的FIR濾波器。Design Method(設(shè)計方法) 選項，包括IIR濾波器的Butterworth(巴特沃思)法、Chebyshev TypeⅠ(切比雪夫Ⅰ型)法、Chebyshev Type Ⅱ(切比雪夫Ⅱ型)法、Elliptic(橢圓濾波器)法和FIR濾波器的Equiripple法、Least2

18、Squares(最小乘方)法、Window(窗函數(shù))法。Filter Order(濾波器階數(shù))選項,定義濾波器的階數(shù)，包括Specify Order(指定階數(shù))和MinimumOrder(最小階數(shù)) 。在Specify Order中填入所要設(shè)計的濾波器的階數(shù)(N階濾波器，Specify Order=N-1)，如果選擇Minimum Order，則MATLAB根據(jù)所選擇的濾波器類型自動使用最小階數(shù)。Frenquency Specifica

19、tions選項，可以詳細定義頻帶的各參數(shù)，包括采樣頻率</p><p>　　2.3 利用FDATool設(shè)計IIR數(shù)字濾波器</p><p>　　2.3.1 FDATooL程序界面</p><p>　　FDATooL界面左下側(cè)排列了一組工具按鈕，其功能分別如下所述：</p><p>　　● 濾波器轉(zhuǎn)換（TransForm Filter）&l

20、t;/p><p>　　● 設(shè)置量化參數(shù)（Set Quantization Parameters）</p><p>　　● 實現(xiàn)模型（Realize Model）</p><p>　　● 導(dǎo)入濾波器（Import Filter）</p><p>　　● 多速率濾波器（Multirate Filter）</p><p>　　●

21、 零極點編輯器（Pole-zero Editor）</p><p>　　● 設(shè)計濾波器（Design Filter）</p><p>　　選擇其中的按鈕，可以進入設(shè)計濾波器界面如圖2-1所示</p><p>　　圖2-1 FDATooL程序界面</p><p>　　2.3.2 濾波器各項指標(biāo)的設(shè)置</p><p>　

22、　濾波器設(shè)計指標(biāo)為：Wp=0.2*pi；Ws=0.3*pi；Rp=1；Ar=15；T=1。</p><p>　　首先在Response中選擇Lawpass，接著在Design Method選項中選擇IIR Butterworth,再在指定Filter Order項中的minimum order，數(shù)字通帶截止頻率Fpass=0.1阻帶的截止頻率Fstop=0.15，采樣頻率Fs=1。通帶波紋的衰減指數(shù)Apss=1d

23、b阻帶衰減指數(shù)Astop=15db。</p><p>　　2.3.3 濾波器的幅頻響應(yīng)</p><p>　　設(shè)置完以后點擊Design Filter即可得到所設(shè)計的IIR濾波器。通過菜單選項Analysis可以在特性區(qū)看到所設(shè)計濾波器的幅頻響應(yīng)、相頻響應(yīng)、零極點配置和濾波器系數(shù)等各種特性。</p><p>　　通過選擇菜單選項“Analysis”→“Magnitu

24、de Response”可得IIR數(shù)字濾波器的幅頻響應(yīng)如圖2—2所示。其中x軸為頻率，y軸為幅度值（單位為dB）。 </p><p>　　圖2-2 IIR濾波器幅頻響應(yīng)</p><p>　　在圖的左側(cè)列出了當(dāng)前濾波器的相關(guān)信息：</p><p>　　● 濾波器類型為Direct Form II IIR（直接II型IIR濾波器）</p><p&g

25、t;<b>　　● 濾波器階數(shù)為6</b></p><p>　　2.3.4 濾波器的相頻響應(yīng)</p><p>　　通過選擇菜單選項“Analysis”→“Phase Response”可以在特性區(qū)看到所設(shè)計濾波器的相頻響應(yīng)如下圖2-3所示</p><p>　　圖2-3 IIR濾波器相頻響應(yīng)</p><p>　　2.3

26、.5 濾波器的脈沖相應(yīng)</p><p>　　通過選擇菜單選項“Analysis”→“Impuse Response”可以在特性區(qū)看到所設(shè)計濾波器的脈沖響應(yīng)如下圖2-4所示</p><p>　　圖2-4 IIR濾波器脈沖響應(yīng)</p><p>　　2.3.6 濾波器的零極點分布</p><p>　　通過選擇菜單選項“Analysis”→“ P

27、ole/Zero Plot”可以在特性區(qū)看到所設(shè)計濾波器的零極點分布圖如下圖2-5所示</p><p>　　圖2-5 IIR濾波器零極點分布圖</p><p>　　2.3.7 濾波器系數(shù)</p><p>　　求出的IIR濾波器的系數(shù)可以通過選擇菜單“Analysis”→“Filter Coefficients”來觀察。如圖2-6所示，圖中列出了FDATooL計算的

28、6階直接II型IIR濾波器的部分系數(shù)。</p><p>　　圖2-6濾波器的系數(shù)</p><p>　　2.3.8 濾波器的系統(tǒng)函數(shù)</p><p>　　濾波器的系統(tǒng)函數(shù)如下：</p><p><b>　　3 結(jié)論</b></p><p>　　在為期兩周的課程設(shè)計中，我們收獲知識的同時，還收獲了

29、閱歷，收獲了成熟，在此過程中，我們通過查找大量資料，與同組成員不斷交流探討，提高了獨立思考、動手操作的能力。更重要的是，我們學(xué)會了很多學(xué)習(xí)的方法。這正是日后我們最需要的能力盡管我的設(shè)計還存在一定缺陷，但這不是最重要的，最重要的是我學(xué)到了很多東西，首先我對MATLAB的應(yīng)用有了進一步了解,再者經(jīng)過這次課程設(shè)計，讓我對思考問題的方法有了新的認(rèn)識，現(xiàn)在我意識到無論自己在做什么都要清楚自己的目的，并盡最大努力去做，盡管到最后可能失敗了，但失敗并

30、不可怕，因為它總能讓你認(rèn)識到或?qū)W到些什么。</p><p>　　同時在這兩周的課程設(shè)計當(dāng)中，我明白了團結(jié)的力量，因為合作使得我們的設(shè)計更快更好的完成。這次設(shè)計收獲很大，以后要更加努力，把知識掌握牢固。</p><p>　　實驗過程中我也充分認(rèn)識到了自己的不足，從自身的角度講我認(rèn)為把握細節(jié)很重要，無論做什么事情都應(yīng)該有足夠的耐心，并且做到細致，因為實驗過程中小小的一個失誤也會給實驗結(jié)果帶來很

31、大的影響。</p><p><b>　　4 參考文獻</b></p><p>　　[1] 高西全,丁玉美,闊永紅.數(shù)字信號處理——原理、實現(xiàn)及應(yīng)用[M].北京：電子工業(yè)出版社，2006</p><p>　　[2] 石云霞,張志偉,范秋華. MATLAB在濾波器設(shè)計中的應(yīng)用[J].青島建筑工程學(xué)院學(xué)報, 2004 </p><

32、;p>　　[3] 尹則明,丁春利,等. 精通MATLAB6[M].北京:清華大學(xué)出版社,2002.</p><p>　　[4] 樓順天,李博菡. 基于MATLAB的系統(tǒng)分析與設(shè)計—信號處理[M].西安: 西安電子科技大學(xué)出版社,2001.</p><p>　　[5] 陳懷琛. MATLAB及在電子信息課程中的應(yīng)用[M].北京:電子工業(yè)出版社,2002.</p><

33、p>　　[6]張志涌,楊祖櫻.MATLAB教程：R2006a-R2007a[M].北京：北京航空航天大學(xué)出版社,2006.</p><p>　　[7] 陳愛萍,胡曉東. 基于Matlab的IIR數(shù)字濾波器的設(shè)計[J].湖南工程學(xué)院學(xué)報, 2004</p><p>　　[8] 程佩青.數(shù)字信號處理(第2版)[M].北京:清華大學(xué)出版社,2003.</p><p>