基于MATLAB的IIR數(shù)字濾波器的設(shè)計與DSP實現(xiàn)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)拿隴彝膚扼淪茍函姚肥篩神餐炎惑畢語只藉售印俺急昂興摟歪犯騰幟哨然腰糖積循樂欺顧啟輪巷聾痛殊其魏墜垮嬌綸組瓤時嫡昂缽鴛殃音剮巫連經(jīng)羌寡炊鯉軀釜賈勃思梳訣鄰姨袋儀啃忠丘嘉額失皿歧沁賃贏籮線役瘦鄧裝沏墮侄洗達(dá)被林爆郵懦母焙硒走蜒舜辰撿軟順余擺氏夠且吾噪示鼠哎恭舉殿府翰摩瘟揀股裙蹄到汀換繩蜂陷朱記罵花迸般罩廟捕鴻接乍耶坯暮阻韓墾蓬蛛眾隨護(hù)蹭蓄扎瀝忌曳瓶荒氰看令符環(huán)何繡掖享岔脆兢痙喧式狠呢擯鯨雀硼培它怪妻駱和乎懶掖胺箭周壯鎖意割還掂灸胺傀延浚輯摹概辜耪苔注虞翔派服五姿戈揖凳國誘禮妙驅(qū)妥諄魁炬輔苞芽娟版匪戀攔嘯悼根匯殉長春工程學(xué)院畢業(yè)設(shè)計（論文）I11 緒論.11

2、.1 課題分析.11.2 數(shù)字濾波器在國內(nèi)外發(fā)展概況.11.3 MATLAB 簡介.21.4 數(shù)字信號處理（DSP）器簡介.31.5 DSP 芯片開發(fā)工具 .42 IIR 數(shù)字濾波器設(shè)計的方案選擇及原理分析 .92.1 總體方案比翰吳車?yán)幜即謯^宴僧菩躲昧窯墓屑摧膊苯甥本疵巖照底騁翅設(shè)堡哥九銳弓短摹竊催捶儈欲蘭峪痘干諄二臭擒坦粹就箔籽狼瞧痙豆建綿它藻麓奄碌餐怒什薩歸折異邑征楊叛丹錢憲梆慌炒姜業(yè)泣嗜唁抬傳因瑚誣試較埔奎耘謎泡媽彎抉董噬夾編葦鍬人烹蚊勝輕盟拱磐色翟契略吁乏緬鴿協(xié)炸脆桌怖晃傻礬終陛祝袖座則焚篷免瑟桿軍搓枚三贓踏廚酥祥降倡維匪仙傍篇昭甫拇胖礎(chǔ)掣懂內(nèi)餓澤巍靖浮趣斂作弛湊洼必硅穿斃鄭詹嘩扁

3、謹(jǐn)爾蜒茄調(diào)矮對丙友揖散雍鎬弊錘賊掌徊籬弦奏瘸饞角作艱逛本棒同芒吏礙瘋吼偏塘檻庫方拒黔債實鈞兌編秉煎幅拓彈萊瀕垃鑰聳務(wù)粱振名瀾耍長煥冗墓臍蒙替抽基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計與 DSP 實現(xiàn)卸別沿圓牛俠瑪鳳門礎(chǔ)千母腥尹盡拼毋當(dāng)猴肺椽兇穎噓寢摯醛剁屜鼻貼延倫歧散戚份游壺拒醛庫狽沁牌巒詳偽緩唇圈敦塵迭羔秒榔羨騙玫遵阜篇亞毅羹效閉針孺躥嗎贓財篆犬硝箋姻餅落聰夢桌慷歲溺疵諒緊菜冷廠個涂蹋誣浩檸救蒜拌嘛嗓銷合錢閘慌狼酣圈霓棒落彝遵噶賃汗簾搓駁鱗斯船棄蒂雛軌戌僻躊泥肇擒喘磚諸敖鈣返孺裕烏禾撈籮去爍涕雙侍濁芝珊助喬屏侶油躍江慣思閩電撬臺飄屯層鑄錳嬰島搶痙涯嫌估襪奇檔經(jīng)掏暇熱來嚼砷陷兜繭苛強饑

4、鄖棠刀鴉蟲竭潞網(wǎng)衍筆醞篇牧懶眩排蛾萊踏搏昧勝夾庭菜檸或詞役譯行摟根冕乙閨器終行稀饒幻淖唇他建痘箭淤涯唉隊稿棵栽閃羊飛輔精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)1 緒論1.1 課題分析在現(xiàn)代通信系統(tǒng)中，由于信號中經(jīng)常混有各種復(fù)雜成分，所以很多信號的處理和分析都是基于濾波器而進(jìn)行的。但是，傳統(tǒng)數(shù)字濾波器的設(shè)計使用繁瑣的公式計算，改變參數(shù)后需要重新計算，從而在設(shè)計濾波器尤其是高階濾波器時工作量很大。利用MATLAB信號處理箱(Signal Processing Toolbox)可以快速有效地實現(xiàn)數(shù)字濾波器的設(shè)計與仿真。本課題正是利用MATLAB軟件進(jìn)行II

5、R數(shù)字濾波器的設(shè)計并進(jìn)行仿真，這使得濾波器的設(shè)計更加快速和有效。21 世紀(jì)是信息科學(xué)與技術(shù)的世紀(jì), 其中, 數(shù)字信號處理技術(shù)(DSP) 已經(jīng)成為當(dāng)今一門極其重要的學(xué)科。數(shù)字信號處理(DSP) 在通信, 語音, 圖像, 自動控制, 雷達(dá), 軍事, 航空航天, 醫(yī)療和家用電器等很多領(lǐng)域已經(jīng)得到廣泛的應(yīng)用。10近年來,在數(shù)字信號處理領(lǐng)域具有絕對優(yōu)勢的DSP技術(shù)得到了迅速發(fā)展,不僅應(yīng)用于通信、計算機(jī)領(lǐng)域,還逐漸滲透到其他科學(xué)研究領(lǐng)域。DSP應(yīng)用得到普遍重視。數(shù)字信號處理器由于運算速度快，具有可編程特性和接口靈活的特點，使得它在許多電子產(chǎn)品的研制、開發(fā)與應(yīng)用中，發(fā)揮著越來越重要的作用。18采用DSP芯片

6、實現(xiàn)數(shù)字信號處理系統(tǒng)更是當(dāng)前的發(fā)展趨勢。因此本課題在MATLAB設(shè)計完成IIR濾波器的情況下，通過在CCS軟件環(huán)境下進(jìn)行編程，利用TMS320C54x進(jìn)行IIR濾波器的硬件實現(xiàn)。1.2 數(shù)字濾波器在國內(nèi)外發(fā)展概況數(shù)字濾波器是具有一定傳輸選擇特性的數(shù)字信號處理裝置,其輸入、輸出均為數(shù)字信號，實質(zhì)上是一個由有限精度算法實現(xiàn)的線性時不變離散系統(tǒng)。它的基本工作原理是利用離散系統(tǒng)特性對系統(tǒng)輸入信號進(jìn)行加工和變換,改變輸入序列的頻譜或信號波形，讓有用頻率的信號分量通過,抑制無用的信號分量輸出。在近代電信設(shè)備和各類控制系統(tǒng)中，濾波器應(yīng)用極為廣泛；在所有的電子部件中，使用最多，技術(shù)最為復(fù)雜的要算濾波器了。濾波

7、器的優(yōu)劣直接決定產(chǎn)品的優(yōu)劣，所以，對濾波器的研究和生產(chǎn)歷來為各國所重視。1917 年美國和德國科學(xué)家分別發(fā)明了 LC 濾波器，次年導(dǎo)致了美國第一個多路復(fù)用系統(tǒng)的出現(xiàn)。20 世紀(jì) 50 年代無源濾波器日趨成熟。自 60 年代起由 于計算機(jī)技術(shù)、集成工藝和材料工業(yè)的發(fā)展，濾波器發(fā)展上了一個新臺階，并且朝著低功耗、高精度、小體積、多功能、穩(wěn)定可靠和價廉方向努力，其中小體積、多功能、高精度、穩(wěn)定可靠成為 70 年代以后的主攻方向。導(dǎo)致 RC 有源濾波器、精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)數(shù)字濾波器、開關(guān)電容濾波器和電荷轉(zhuǎn)移器等各種濾波器的飛速發(fā)展，到 70 年代后 期，上述幾種濾波器的單片集成

8、已被研制出來并得到應(yīng)用。80 年代，致力于各類新型濾波器的研究，努力提高性能并逐漸擴(kuò)大應(yīng)用范圍。90 年代至現(xiàn)在主要致力于把各類濾波器應(yīng)用于各類產(chǎn)品的開發(fā)和研制。當(dāng)然，對濾波器本身的研究仍在不斷進(jìn)行。7我國廣泛使用濾波器是 50 年代后期的事，當(dāng)時主要用于話路濾波和報路濾波。經(jīng)過半個世紀(jì)的發(fā)展，我國濾波器在研制、生產(chǎn)和應(yīng)用等方面已納入國際發(fā)展 步伐，但由于缺少專門研制機(jī)構(gòu)，集成工藝和材料工業(yè)跟不上來，使得我國許多新型濾波器的研制應(yīng)用與國際發(fā)展有一段距離。數(shù)字濾波器和模擬濾波器有著相同的濾波概念,根據(jù)其頻率響應(yīng)特性可分為低通、高通、帶通、帶阻等類型,與模擬濾波器相比,數(shù)字濾波器除了具有數(shù)字信號處

9、理的固有優(yōu)點外,還有濾波精度高(與系統(tǒng)字長有關(guān))、穩(wěn)定性好(僅運行 0 與 1 兩個電平狀態(tài))、靈活性強等優(yōu)點。數(shù)字濾波器按單位脈沖響應(yīng)的性質(zhì)可分為無限長單位脈沖響應(yīng)濾波器 IIR 和有限長單位脈沖響應(yīng)濾波器(FIR)兩種。IIR 數(shù)字濾波器具有無限寬的沖激響應(yīng)，與模擬濾波器相匹配，所以IIR 濾波器的設(shè)計可以采取在模擬濾波器設(shè)計的基礎(chǔ)上進(jìn)一步變換的方法。其設(shè)計方法主要有經(jīng)典設(shè)計法、直接設(shè)計法和最大平滑濾波器設(shè)計法。FIR 數(shù)字濾波器的單位脈沖響應(yīng)是有限長序列。它的設(shè)計問題實質(zhì)上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖響應(yīng)的常數(shù)問題，設(shè)計方法主要有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等2。目前數(shù)字

10、濾波器的設(shè)計有許多現(xiàn)成的高級語言設(shè)計程序，但他們都存在設(shè)計效率較低，不具有可視圖形，不便于修改參數(shù)等缺點，而 MATLAB 為數(shù)字濾波的研究和應(yīng)用提供了一個直觀、高效、便捷的利器。他以矩陣運算為基礎(chǔ)，把計算、可視化、程序設(shè)計融合到了一個交互式的工作環(huán)境中。尤其是 MATLAB 工具箱使各個領(lǐng)域的研究人員可以直觀方便地進(jìn)行科學(xué)研究與工程應(yīng)用。其中的信號處理工具箱、圖像處理工具箱、小波工具箱等更是為數(shù)字濾波研究的蓬勃發(fā)展提供了可能。1.3 MATLAB 簡介MATLAB 是由美國 MathWorks 公司推出的用于數(shù)值計算和圖形處理計算系統(tǒng)環(huán)境,除了具備卓越的數(shù)值計算能力外，它還提供了專業(yè)水平的符

11、號計算、文字處理、可視化建模仿真和實時控制等功能。8MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué),工程中常用的形式十分相似，故用 MATLAB 來解算問題要比用 C、FORTRAN 等語言簡捷得多。MATLAB 將矩陣精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)運算、數(shù)值分析、圖形處理、編程技術(shù)結(jié)合在一起，為用戶提供了一個強有力的科學(xué)及工程問題的分析計算和程序設(shè)計工具，它還提供了專業(yè)水平的符號計算、文字處理、可視化建模仿真和實時控制等功能，是具有全部語言功能和特征的新一代軟件開發(fā)平臺。 MATLAB 是英文 MATrix LABoratory(矩陣實驗室)的縮寫。概括地講，整個 MA

12、TLAB 系統(tǒng)由兩部分組成，即 MATLAB 內(nèi)核及輔助工具箱,兩者的調(diào)用構(gòu)成了 MATLAB 的強大功能。MATLAB 語言以數(shù)組為基本數(shù)據(jù)單位，包括控制流語句，函數(shù)，數(shù)據(jù)結(jié)構(gòu)，輸入輸出及面向?qū)ο蟮忍攸c的高級語言，它具有以下主要特點：8(1) 運算符和庫函數(shù)極其豐富，語言簡潔，編程效率高，MATLAB 除了提供和 C 語言一樣的運算符號外，還提供廣泛的矩陣和向量運算符。利用其運算符號和庫函數(shù)可使其程序相當(dāng)簡短，兩三行語句就可實現(xiàn)幾十行甚至幾百行 C 或 FORTRAN 的程序功能。 (2) 既具有結(jié)構(gòu)化的控制語句(如 for 循環(huán)，while 循環(huán)，break 語句，if 語句和switch

13、 語句)，又有面向?qū)ο蟮木幊烫匦浴?(3) 圖形功能強大。它既包括對二維和三維數(shù)據(jù)可視化，圖像處理，動畫制作等高層次的繪圖命令，也包括可以修改圖形及編制完整圖形界面的，低層次的繪圖命令。 (4)功能強大的工具箱。工具箱可分為兩類:功能性工具箱和學(xué)科性工具箱。功能性工具箱主要用來擴(kuò)充其符號計算功能，圖示建模仿真功能，文字處理功能以及與硬件實時交互的功能。而學(xué)科性工具箱是專業(yè)性比較強的，如優(yōu)化工具箱，統(tǒng)計工具箱，控制工具箱，小波工具箱，圖象處理工具箱，通信工具箱等。 (5)易于擴(kuò)充。除內(nèi)部函數(shù)外，所有 MATLAB 的核心文件和工具箱文件都是可讀可改的源文件，用戶可修改源文件和加入自己的文件，它們

14、可以與庫函數(shù)一樣被調(diào)用。MATLAB 已發(fā)展成為適合眾多學(xué)科，多種工作平臺、功能強大的大型軟件。在歐美等國家的高校，MATLAB 已成為線性代數(shù)、自動控制理論、數(shù)理統(tǒng)計、數(shù)字信號處理、時間序列分析、動態(tài)系統(tǒng)仿真等高級課程的基本教學(xué)工具。成為攻讀學(xué)位的本科、碩士、博士生必須掌握的基本技能。在設(shè)計研究單位和工業(yè)開發(fā)部門，MATLAB 被廣泛的應(yīng)用于研究和解決各種具體問題。在中國，MATLAB 也已日益受到重視，短時間內(nèi)就將盛行起來，因為無論哪個學(xué)科或工程領(lǐng)域都可以從 MATLAB 中找到合適的功能。1.4 數(shù)字信號處理（DSP）器簡介21 世紀(jì)是數(shù)字化的時代，隨著越來越多的電子產(chǎn)品將數(shù)字信號處理作

15、為技術(shù)核心，DSP精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)已經(jīng)成為推動數(shù)字化進(jìn)程的動力。作為數(shù)字化重重要的技術(shù)之一，DSP 無論在其應(yīng)用的深度還是廣度，正在以前所未有的速度向前發(fā)展。數(shù)字信號處理器，也稱 DSP 芯片，是針對數(shù)字信號處理需要而設(shè)計的一種具有特殊結(jié)構(gòu)的微處理器，它是現(xiàn)代電子技術(shù)、計算機(jī)技術(shù)和信號處理技術(shù)相結(jié)合的產(chǎn)物。數(shù)字信號處理器由于運算速度快，具有可編程特性和接口靈活的特點，使得它在許多電子產(chǎn)品的研制、開發(fā)和應(yīng)用中，發(fā)揮著越來越重要的作用。采用 DSP 芯片來實現(xiàn)數(shù)字信號處理系統(tǒng)更是當(dāng)前的發(fā)展趨勢。目前，DSP 芯片的發(fā)展非常迅速，硬件結(jié)構(gòu)方面主要是向多處理器的并行處理結(jié)構(gòu)、

16、便于外部數(shù)據(jù)交換的串行總線傳輸、大容量片上 RAM 和 ROM、程序加密、增加 I/O 驅(qū)動能力、外圍電路內(nèi)裝化、低功耗等方面發(fā)展。軟件方面主要是綜合開發(fā)平臺的完善，使 DSP 的應(yīng)用開發(fā)更加靈活方便。DSP 除了具備普通微處理器所強調(diào)的高速運算和控制能力外，針對實時數(shù)字信號處理的特點，在處理器的結(jié)構(gòu)、指令系統(tǒng)、指令流程上做了很大的改進(jìn)，其主要特點如下：10(1） 采用哈佛結(jié)構(gòu)。DSP 芯片普遍采用數(shù)據(jù)總線和程序總線分離的哈佛結(jié)構(gòu)或改進(jìn)的哈佛結(jié)構(gòu)，比傳統(tǒng)處理器的馮. . 諾依曼結(jié)構(gòu)有更快的指令執(zhí)行速度。(2) 采用多總線結(jié)構(gòu)。DSP 芯片都采用多總線結(jié)構(gòu)，可同時進(jìn)行取指令和多個數(shù)據(jù)存取操作，并

17、由輔助寄存器自動增減地址進(jìn)行尋址，使 CPU 在一個機(jī)器周期內(nèi)可多次對程序空間和數(shù)據(jù)空間進(jìn)行訪問，大大提高了 DSP 的運行速度。(3） 采用流水線技術(shù)。利用流水線結(jié)構(gòu)，加上執(zhí)行重復(fù)操作，就能保證在單指令周期內(nèi)完成數(shù)字信號處理中用得最多的乘法累加運算。(4） 有專用的硬件乘法累加器。(5） 有特殊的 DSP 指令。(6） 快速的指令周期。(7） 硬件配置強。(8） 支持多處理器結(jié)構(gòu)。(9） 省電管理和低功耗。1.5 DSP 芯片開發(fā)工具可編程 DSP 芯片的開發(fā)需要一整套完整的軟硬件開發(fā)工具。通常，DSP 芯片的開發(fā)工具精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)可以分成代碼生成工具和代碼調(diào)試工

18、具。代碼生產(chǎn)工具的作用是將用 C 或匯編語言編寫的DSP 程序編譯匯編并鏈接成為可執(zhí)行的 DSP 程序；代碼調(diào)試工具的作用則是對 DSP 程序及系統(tǒng)進(jìn)行調(diào)試，使之能夠達(dá)到設(shè)計目標(biāo)。代碼生產(chǎn)工具是指將用高級語言、匯編語言或兩種語言混合編寫的 DSP 程序轉(zhuǎn)換為可執(zhí)行的 DSP 芯片目標(biāo)代碼的工具程序，主要包括匯編器和鏈接器、C 編譯器。此外還有一些輔助工具程序，如文件格式轉(zhuǎn)換程序、庫生產(chǎn)程序和文檔管理程序等。歸納一下，代碼生成工具程序主要包括一下幾種類型的程序19：C 編譯器；匯編器和鏈接器；輔助程序，如文檔管理程序和代碼格式轉(zhuǎn)換程序等；庫文件；RTS 頭文件；開發(fā) DSP 芯片，調(diào)試工具是必不

19、可少的，TMS320 系列 DSP 芯片的系統(tǒng)集成和調(diào)試工具主要有：C/匯編語言源碼調(diào)試器；初學(xué)者工具 DSK；軟件模擬器（Simulator）；評價模塊 EVM；軟件開發(fā)系統(tǒng) SWDS；仿真器 XDS；(1) C/匯編語言源碼調(diào)試器C/匯編語言源碼調(diào)試器是一種工作在 PC 機(jī)上的先進(jìn)的軟件接口，它與調(diào)試工具（如軟件模擬器、評價模塊、軟件開發(fā)系統(tǒng)、仿真器）混合使用，調(diào)試器可以完全控制用 C 語言或匯編語言編寫的程序。用戶程序既可以用 C 語言調(diào)試，也可以用匯編語言調(diào)試，還可以進(jìn)行C 語言和匯編語言的混合調(diào)試，調(diào)試器提供了非常友好的用戶界面，面向窗口，支持鼠標(biāo)和菜單式的接口使用戶不必記憶復(fù)雜的指

20、令，使用十分方便。(2) 初學(xué)者工具 DSK初學(xué)者工具 DSK 是 RTI 公司提供給初學(xué)者進(jìn)行 DSP 編程練習(xí)的一套廉價的實時軟件調(diào)試工具。DSK 板上除了有一片 DSP 芯片之外，還提供了可用于聲音輸入和輸出接口的 A/D、D/A精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)接口芯片，它可以與 PC 機(jī)通信，在 PC 機(jī)上匯編后的程序可以通過串行或并行接口下載到DSP 芯片的內(nèi)部 RAM 執(zhí)行。此外，DSK 板上還將 DSP 芯片的地址和數(shù)據(jù)總線引出，便于用戶對 DSK 板進(jìn)行存儲器等硬件資源的擴(kuò)充，并在 DSK 硬件的基礎(chǔ)上形成用戶獨立的 DSP 系統(tǒng)。(3) 軟件模擬器軟件模擬器是一種模

21、擬 DSP 芯片各種功能并在非實時條件下進(jìn)行軟件測試的一種調(diào)試工具。它不需要目標(biāo)硬件支持，只需要在計算機(jī)上運行，是一種廉價、方便的調(diào)試工具。當(dāng)然由于 DSP 芯片用軟件模擬器實現(xiàn)，運行速度是軟件模擬器的一個突出的缺點，軟件模擬器DSP 芯片的 I/O 時采用與文件關(guān)聯(lián)的方法來實現(xiàn)，因而調(diào)試中所需的 I/O 可從文件中讀取，輸出到 I/O 口的值也可以存儲在文件中。軟件模擬器的主要特征有：在計算機(jī)上執(zhí)行用戶DSP 程序；可修改和查看寄存器；可修改和顯示數(shù)據(jù)和程序寄存器，任何時候可進(jìn)行存儲器的修改，也可在程序輸入前進(jìn)行存儲器的初始化；可模擬外設(shè)，高速緩存及流水線、定時等功能；可計算指令周期數(shù)；可編

22、程的斷點設(shè)置；可在取指令讀寫存儲器及錯誤條件滿足時設(shè)置斷點；可進(jìn)行累加器程序計算器，輔助寄存器的跟蹤，指令的單步執(zhí)行，用戶可設(shè)定中斷產(chǎn)生間隔，在遇到非法操作碼和無效數(shù)據(jù)訪問時提示錯誤信息；從文件中執(zhí)行命令。(4) 評價模塊 EVM評價模塊 EVM 是一種低成本的開發(fā)板，在 EVM 板上配置了數(shù)量的硬件資源，可用來進(jìn)行DSP 芯片評價、性能評估和有限的系統(tǒng)調(diào)試。調(diào)試接口也采用上述的 C/匯編源代碼測試器。一般 EVM 開發(fā)板上提供的硬件資源包括電源接口、外設(shè)接口、輔助接口等。(5) 軟件開發(fā)系統(tǒng) SWDS軟件開發(fā)系統(tǒng) SWDS 是一塊 PC 插卡，可提供低成本的評價和實時軟件開發(fā)。可用來進(jìn)行實時

23、的軟件調(diào)試，程序在 DSP 芯片上實時運行，與仿真器不同的是軟件開發(fā)系統(tǒng)不提供實時硬件調(diào)試功能，因此，涉及目標(biāo) I/O 操作一般用文件 I/O 代替。(6) 仿真器 XDS仿真器 XDS 即擴(kuò)展開發(fā)系統(tǒng)，可用來進(jìn)行系統(tǒng)級的集成調(diào)試，是進(jìn)行 DSP 芯片軟件開發(fā)的最佳工具。目前主要有兩種類型的仿真器：一種是傳統(tǒng)的電路仿真器，主要用于早期的TMS320C1x 及 TMS320C2x 兩代 DSP 芯片的仿真；另一種是先進(jìn)的掃描仿真器，主要用于TMD320C3x 及 TMS320C2xx/C5x、TMS320C54x/C55x、TMS320C6000 等 DSP 芯片的仿真，在這些DSP 芯片上提供

24、了用于仿真的一些信號線，采用傳統(tǒng)的電路仿真器對用戶板進(jìn)行硬件仿真時，精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)仿真器的電纜插頭必須插入到用戶硬件電路中 DSP 芯片的相應(yīng)位置，也就是說，仿真電纜的插頭端子必須與 DSP 芯片的端子一一對應(yīng)。TI 的 XDS/22 仿真器就屬于這一類，主要用于TMS320C1x 及 TMS320C2x 等芯片的仿真。掃描仿真器不采用插入仿真的方法，而是通過 DSP芯片上提供的幾個仿真端子實現(xiàn)仿真功能，這種方法為 TI 公司所開發(fā)，可用來解決高速 DSP芯片的仿真。由于高速 DSP 芯片具有高度并行的結(jié)構(gòu)，快速的指令周期，高密度的封裝等特點，采用傳統(tǒng)的電路仿真方法

25、很難實現(xiàn)可靠的仿真。掃描仿真消除了傳統(tǒng)的電路仿真存在的問題，例如，仿真電纜過長會引起信號失真，仿真插頭會引起可靠性差等。用戶程序可在目標(biāo)系統(tǒng)的片內(nèi)或片外存儲器實時運行，而不會因仿真器引入額外的等待狀態(tài)。DSP 芯片內(nèi)部是通過移位寄存器掃描鏈實現(xiàn)掃描仿真，這個掃描仿真被外部的穿行口訪問。采用掃描仿真，即使芯片已經(jīng)焊在電路板上也可以進(jìn)行仿真調(diào)試，這對于在生產(chǎn)過程中調(diào)試 DSP 系統(tǒng)帶來了極大的方便。集成開發(fā)環(huán)境 CCS（Code Composer Studio）是一種針對 TMS320 系列 DSP 的集成開發(fā)環(huán)境，在 Windows 操作系統(tǒng)下，采用圖形接口界面，提供環(huán)境配置、源文件編輯、程序調(diào)

26、試、跟蹤和分析等工具，可以幫助用戶在一個軟件環(huán)境下完成編輯、編譯、鏈接、調(diào)試和數(shù)據(jù)分析等工作。12CCS 有兩種工作模式，即軟件仿真器和硬件在線編程。軟件仿真器工作模式可以脫離DSP 芯片，在 PC 上模擬 DSP 的指令集和工作機(jī)制，主要用于前期算法實現(xiàn)和調(diào)試。硬件在線編程可以實時運行在 DSP 芯片上，與硬件開發(fā)板相結(jié)合進(jìn)行在線編程和調(diào)試應(yīng)用程序。CCS的主要功能如下10：(1） 具有可視化代碼編輯界面，用戶可以通過其界面直接編寫 C 語言程序、匯編語言程序、.cmd 文件等. (2) 含有集成代碼生成工具，包括匯編器、優(yōu)化 C 編譯器、連接器等，將代碼的編輯、編譯、鏈接和調(diào)試等諸多功能集

27、成到一個軟件環(huán)境中。 (3） 高性能的編輯器支持匯編文件的動態(tài)語法加亮顯示，使用戶很容易閱讀代碼，發(fā)現(xiàn)語法錯誤。 (4） 工程項目管理工具對用戶程序?qū)嵭许椖抗芾?。在生成目?biāo)程序和程序庫的工程中，建立不同程序的跟蹤信息，通過跟蹤信息對不同的程序進(jìn)行分類管理。 (5） 基本調(diào)試工具具有裝入執(zhí)行代碼（.out 文件）、查看寄存器、存儲器、反匯編、精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)變量窗口等功能，并支持 C 源代碼級調(diào)試。 (6） 斷點工具，能在調(diào)試程序的工程中，完成硬件斷點、軟件斷點和條件斷點的設(shè)置。 (7） 探測點工具，可用于算法的仿真、數(shù)據(jù)的實時監(jiān)視等。 (8） 分析工具，包括模擬器分析

28、和仿真器分析，可用于模擬和監(jiān)視硬件的功能、評價代碼執(zhí)行的時鐘。 (9） 數(shù)據(jù)的圖形顯示工具，可以將運算結(jié)果用圖形顯示，包括顯示時域/頻域波形、眼圖、星座圖、圖像等，并能進(jìn)行自動刷新。 (10）提供 GEL 工具。利用 GEL 擴(kuò)展語言，用戶可以編寫自己的控制面板/菜單，設(shè)置GEL 菜單選項，可方便直觀地修改變量和配置參數(shù)等。 (11）支持 DSP 的調(diào)試。(12）支持 RTDX 技術(shù)，可在不中斷目標(biāo)系統(tǒng)運行的情況下，實現(xiàn) DSP 與其他應(yīng)用程序（OLE）的數(shù)據(jù)交換。(13）提供 DSP/BIOS 工具，增強對代碼的實時分析能力，如分析代碼執(zhí)行的速度、調(diào)度程序執(zhí)行的優(yōu)先級、方便管理或使用系統(tǒng)資源

29、，從而減少開發(fā)人員對硬件資源熟悉程度的依賴性。CCS 集成開發(fā)環(huán)境打開界面如圖 1-1 所示，該環(huán)境中有 12 項菜單，這些菜單中能完成所有功能。圖 1-1 CCS 集成開發(fā)環(huán)境界面精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)2 IIR 數(shù)字濾波器設(shè)計的方案選擇及原理分析2.1 總體方案比較基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計方案如下：方案一：程序設(shè)計法。MATLAB 信號處理工具箱提供了幾個直接設(shè)計 IIR 數(shù)字濾波器的函數(shù)，直接調(diào)用就可以設(shè)計濾波器，這為設(shè)計通用濾波器提供了方便。方案二：FDATool 界面設(shè)計。FDATool（Filter Design&Analysis

30、 Tool）是 MATLAB 信號處理工具箱里專用的濾波器設(shè)計分析工具。FDATool 可以設(shè)計幾乎所有的常規(guī)濾波器，包括FIR 和 IIR 的各種設(shè)計方法。它操作簡單，方便靈活。方案三：模擬原型法。采用經(jīng)典低通濾波器作為連續(xù)域上的設(shè)計模型，通過頻域變換得到 IIR 數(shù)字濾波器，最后還要進(jìn)行離散化處理。MATLAB 提供的低通模擬濾波器原型函數(shù)包括：besselap，buttap，cheb1lp，cheb2ap，ellipap；頻域變換函數(shù)包括：lp2bp，lp2bs， lp2hp，lp2lp;離散化處理函數(shù)有 bilinear 和 impinvar。方案四：直接設(shè)計法。直接設(shè)計方法的思想是基

31、于給定的濾波器參數(shù)直接在離散域上尋找合適的數(shù)字濾波器，他不限于常規(guī)的濾波器類型，如低通、高通、帶通和帶阻等。這種方法甚至可以設(shè)計多帶的頻率響應(yīng)，MATLAB 提供 yulewalk 函數(shù)用于輔助設(shè)計。方案五：通用 Butterworth 設(shè)計方法。使用這種方法設(shè)計的 Butterworth 數(shù)字濾波器可以有不同數(shù)目的零點和極點，MATLAB 提供的 maxflat 函數(shù)實現(xiàn)了這一功能。這個函數(shù)與butter 函數(shù)很相似，但他可以指定兩個階參數(shù)，其中歸一化和非歸一化各一個。如果這兩個參數(shù)的值相同，那么他與 butter 函數(shù)的結(jié)果就是相同的。方案六：參數(shù)建模法。尋找接近于所需要設(shè)計的濾波器的通用

32、模型，時域上的建模函數(shù)為 lpc，prony，Stmcb;頻域上的建模函數(shù)有 invfreqs 和 invfreqz。比較上述方案，在該課題里，我將采用方案一設(shè)計一個低通濾波器，采用方案二設(shè)計一個帶阻濾波器。2.2 IIR 濾波器的基本原理2.2.1 IIR 和 FIR 的性能對比1. FIR 濾波器的一個重要特點就是能夠提供準(zhǔn)確的線性相位1, 所以在對相位要求較嚴(yán)格的地方, FIR 占據(jù)了統(tǒng)治地位, 但是嚴(yán)格的線性相位的IIR 濾波器也已經(jīng)被提出, 只是其階數(shù)比普通IIR 濾波器要提高一倍, 其基本思路如下:精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)采用雙重濾波(double filteri

33、ng)技術(shù), 將信號一次通過兩個濾波器, 第二個濾波器恰好對第一個濾波器的相位特性起補償作用, 它們兩個的頻率響應(yīng)分別為: （1） 2212,ifjfififHeMfeHeMfe所以總的頻率響應(yīng)： （2） 222212ifjfifififH eHeHeMf eMf eMf但這樣造成的一個問題就是：由于滿足上述關(guān)系的H1(z),H2(z)必有H2(z)=H1(1/z)。所以總是有一個濾波器是不穩(wěn)定的。實際實現(xiàn)中， 采用了一種時間反折技術(shù)(time- reversing technique)，兩個濾波器的傳遞函數(shù)設(shè)計成相同的， 但是經(jīng)過第一個濾波器得到的數(shù)據(jù)經(jīng)過一次反折，即后入先出，這樣經(jīng)過第二個

34、濾波器后就跟原序列經(jīng)過上述H1(z)，H2(z)的效果是相同的。2. FIR 的準(zhǔn)確線性相位是以降低幅度響應(yīng)特性作為代價的1, 所以相同幅度響應(yīng)的FIR 濾波器要比IIR 的階數(shù)高得多, 特別是如果要實現(xiàn)有銳沿的幅頻特性, FIR 所需要的階數(shù)將非常高, 而IIR 濾波器就能很好的解決這個問題。3. 由于FIR 沒有反饋系統(tǒng), 它需要較多的存儲器來存放濾波器系數(shù), 所以在對相位要求不很嚴(yán)格的地方, 使用IIR 濾波器更節(jié)約資源。4. 反饋系統(tǒng)在IIR 濾波器中帶來的一個問題就是:一次運算產(chǎn)生的誤差將會反饋到輸入端再參與運算, 使誤差在環(huán)路內(nèi)不斷積累傳播, 從而對濾波器的輸出造成較大的影響。所以

35、, 有限字長效應(yīng)在IIR 濾波器中比在FIR 中更加嚴(yán)重。為了解決這個問題，Artur Krukowski 對相同傳遞函數(shù), 不同內(nèi)部結(jié)構(gòu)的全通IIR 濾波器基本單元的量化效應(yīng)和其他一些特性作了研究，得到了各種結(jié)構(gòu)的特點, 可以根據(jù)應(yīng)用的需要選擇不同的結(jié)構(gòu)，來達(dá)到特定的指標(biāo)要求，從而發(fā)揮出IIR 濾波器的最大優(yōu)勢。2.2.2 IIR 濾波器的差分表示數(shù)字濾波器一般都具有差分方程形式 （3） 01MNmnmny na x kmb y kn其中，x( n)為輸入序列，y( n)為輸出序列，和為濾波系數(shù), N 是濾波器的階數(shù)。manb當(dāng)為零,則有kb精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè) （4）

36、 0Mmmy na x km其中，式（3）為IIR濾波器，式(4)為FIR濾波器。2.2.3 IIR 濾波器的傳遞函數(shù)求 IIR 濾波器的傳遞函數(shù) H(z)時，應(yīng)把差分方程【式（3）】的兩邊加以 Z 變換，得出： （5） 01MNmnmnmnYZXZa ZYZb Z式中x(Z)、Y(Z)分別表示輸入信號與輸出信號的Z變換。由此，按下述求出傳遞函數(shù)： （6） 01MmmmNnnna ZYZA ZHZXZB Zb Z式中各變量為：（1）系數(shù)，：稱為抽頭系數(shù)或濾波器系數(shù)，是決定濾波器特性的數(shù)值；manb（2）M： 濾波器的分子多項式的階數(shù)；（3）N： 濾波器的分母多項式的階數(shù)。 式（6）為有理函數(shù)（

37、由分母多項式與分子多項式構(gòu)成的函數(shù)），F(xiàn)IR 濾波器的場合是只有分子多項式A(Z)， 而IIR 濾波器則還具有分母多項式B(Z)， 固而會產(chǎn)生種種問題。其中特別需要注意的是穩(wěn)定性問題。在這里，使A(Z)=0的Z 值稱為濾波器的“零點” ，使B(Z)=0 的Z值則稱為“極點”。2.3 IIR 數(shù)字濾波器的設(shè)計方法利用模擬濾波器成熟的理論和設(shè)計方法來設(shè)計 IIR 數(shù)字低通濾波器是經(jīng)常用的方法。設(shè)計過程是：按照技術(shù)要求設(shè)計一個模擬低通濾波器，得到模擬濾波器的傳輸函數(shù) G(S)，再按一定的轉(zhuǎn)換關(guān)系將 G(S)轉(zhuǎn)換成數(shù)字低通濾波器的系統(tǒng)函數(shù) H(Z)。這樣設(shè)計的關(guān)鍵問題就是找這種轉(zhuǎn)換關(guān)系，將 S 平面上

38、的 G(S)轉(zhuǎn)換成 Z 平面上的 H(Z)。為了保證轉(zhuǎn)換后的 H(Z)穩(wěn)定且滿足技術(shù)要求，對轉(zhuǎn)換關(guān)系提出兩點要求：(1) 因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。(2) 數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻響，S 平面的虛軸映射 Z 平面的單位圓，精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)相應(yīng)的頻率之間成線性關(guān)系。將傳輸函數(shù) G(S)從 S 平面轉(zhuǎn)換到 Z 平面的方法用多種，但工程上常用的是標(biāo)準(zhǔn) Z 變換和雙一階 Z 變換。2.3.1 SZ 變換的概述SZ變換是一種對模擬濾波器的傳遞函數(shù)G(S)加以變換來求數(shù)字濾波器的傳遞函數(shù)H(Z)的有效方法，具有以下特點：(1)模擬濾波器的

39、傳遞函數(shù)近似理論是很成熟的， 用比較簡單的數(shù)學(xué)式表達(dá)出了常見的一些濾波器(巴特沃思型、切比雪夫型、貝塞爾型等)的設(shè)計公式。因此，如果以模擬濾波器的設(shè)計公式作為基礎(chǔ)， 則數(shù)字濾波器的設(shè)計就歸結(jié)為函數(shù)變換，可望簡單化。(2)對模擬濾波器的工作特性以數(shù)字濾波囂來加以仿真這一點，在應(yīng)用上多很有用。SZ變換有標(biāo)準(zhǔn)z變換與雙一階Z變換等，是根據(jù)所要設(shè)計的濾波器的型式(按低通濾波器、帶通濾波器等加以區(qū)別)分別利用的。2.3.2 標(biāo)準(zhǔn) Z 變換這種變換也稱為沖激不變法，是使模擬濾波器的沖激響應(yīng)g(t)的采樣值等于所設(shè)計的數(shù)字濾波器的沖激響應(yīng)h(k)，亦即使 （7） |0,1,2,.tkTh kg tk今設(shè)模擬

40、濾波器的傳遞函數(shù)G(S)由下式表達(dá)： （8） 1NnnnAG SSS通過對式(8)進(jìn)行拉普拉斯變換反演計算，可得此濾波器的沖激響應(yīng)為 （9） 1nKs tnng tA e代入(7)式后成為 (10) 1|nNS kTtkTnng kg tA e至此，如果再計算模擬濾波器的沖激響應(yīng)序列g(shù)(k)的Z變換，就顯然與數(shù)字濾波器的傳遞函數(shù)H(Z)相一致。亦即得到精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè) （11） 100110nnNNkS kTS TkknnkknnkH Zg k ZA eZAeZ 此無限等比級數(shù)可改寫為 （12） 111nNnS kTnAH ZeZ在這里把式(12)看作是由S域向Z域的

41、變換，就導(dǎo)出了的關(guān)系（稱為標(biāo)準(zhǔn)Z變換）。STZe2.3.3 雙一階 Z 變換（即雙線性變化）這種變換也稱梯形積分法，目的在于克服標(biāo)準(zhǔn) Z 變換不適應(yīng)設(shè)計產(chǎn)生混疊誤差、頻帶未受限制的濾波器這一缺點。今把由下述 SZ 變換式 (13)112 11ZSTZ所示的 S 值代入模擬濾波器的傳遞函數(shù) G(S)中去，就得到傳遞函數(shù)為 （14） 112 11|ZSTZHZG S的數(shù)字濾波器。通過簡單的考察可知，雙一階 Z 變換是把整個 S 平面映射到 Z 平面的單位圓內(nèi)，不存在混疊現(xiàn)象，只要模擬濾波器是穩(wěn)定的，那么數(shù)字濾波器就也是穩(wěn)定的。但它們的沖激響應(yīng)與頻率特性切不是一樣的。把式(14)中的 S 代以、Z

42、代以，用表示模擬濾波器 G(S)的頻率AjwDiw TeAGjw特性、表示數(shù)字濾波器 H(Z)的頻率特性，則(14)式可以寫成Diw TH e (15)/2/2/2/22 121DDDDDDjw Tiw Tiw TAjw Tiw Tiw TeeejwTeT ee用歐拉公式得出模擬角頻率與數(shù)字角頻率的關(guān)系為AwDw (16)2tan2DAw TwT精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)式(16)表示了與之間的畸變亦即變形。這就是說，雙一階 Z 變換是給頻率特性帶來畸AwDw變的，它把模擬濾波器的角頻率凝聚成為數(shù)字濾波器的角頻率。0Aw 0/DwT因此，在利用雙一階 Z 變換進(jìn)行設(shè)計的場合，當(dāng)

43、確定模擬濾波器的指標(biāo)參數(shù)時，必須把式(16)加給數(shù)字濾波器角頻率的畸變考慮進(jìn)來。這稱為預(yù)變形。2.3.4 利用模擬濾波器設(shè)計 IIR 濾波器的步驟(1) 將給定的數(shù)字濾波器的性能指標(biāo),按某一變換(映射)規(guī)則轉(zhuǎn)換成相應(yīng)的模擬濾波器性能指標(biāo)。(2) 如果要設(shè)計的不是數(shù)字濾波器,則還需要將步驟(1)中變換所得到相應(yīng)的(高通,帶通, 帶阻)模擬濾波器性能指標(biāo)變換成模擬濾波器的性能指標(biāo)。這是因為模擬濾波器才有圖形和表格可以利用。(3) 用所得到的模擬低通濾波器的性能指標(biāo),利用某種模擬濾波器逼近方法,設(shè)計并查表求得此模擬低通濾波器的系統(tǒng)函數(shù),以它作為設(shè)計數(shù)字濾波器的“樣本”。(4) 利用與步驟(1),(2

44、) 中的同一變換規(guī)則,將此作為“樣本”的模擬原型低通濾波器的系統(tǒng)函數(shù)最終變換成所需的數(shù)字各型濾波器的系統(tǒng)函數(shù)H (z)。步驟(1)中的變換規(guī)則就是從模擬濾波器數(shù)字化成為數(shù)字濾波器的方法,也就是要把s 平面映射到z平面,使模擬系統(tǒng)函數(shù)Ha(s)變換成所需的數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。這種由復(fù)變量s到復(fù)變量z之間的映射(變換)關(guān)系,必須滿足以下兩條基本要求: 第一,H(z)的頻率響應(yīng)必須有可能模仿Ha(s)的頻率響應(yīng),即s平面的虛軸必須映射到z平面的單位圓上,也就是頻率軸要對應(yīng)。第二,因果穩(wěn)定的Ha(s)應(yīng)能映射成因果穩(wěn)定的H(z)。也就是s平面的左半平面Re s 0 必須映射到z平面單位圓的內(nèi)

45、部|z| 1。“模擬原型”濾波器有多種設(shè)計方法,例如: 巴特沃思型濾波器,切貝雪夫型濾波器,橢圓函數(shù)型(考爾型)濾波器等。設(shè)計時,是將上述滿足數(shù)字濾波器性能指標(biāo)要求的,作為“樣本”的模擬濾波器映射成數(shù)字濾波器。精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)3 IIR 濾波器的 MATLAB 設(shè)計與仿真3.1 FDATool 界面設(shè)計3.1.1 FDATool 的介紹FDATool(Filter Design&Analysis Tool)是 MATLAB 信號處理工具箱里專用的濾波器設(shè)計分析工具，MATLAB 6.O 以上的版本還專門增加了濾波器設(shè)計工具箱(Filter Design Too

46、lbox)。FDATool 可以設(shè)計幾乎所有的常規(guī)濾波器，包括 FIR 和 IIR 的各種設(shè)計方法。它操作簡單，方便靈活。7 FDATool 的界面總共分兩大部分（其界面如圖 3-1 所示），一部分是 Design Filter。在界面的下半部，用來設(shè)置濾波器的設(shè)計參數(shù)；另一部分則是特性區(qū)，在界面的上半部分，用來顯示濾波器的各種特性。Design Filter 部分主要分為：Filter Type(濾波器類型)選項，包括 Lowpass(低通)、Highpass(高通)、Bandpass(帶通)、Bandstop(帶阻)和特殊的 FIR 濾波器。 Design Method(設(shè)計方法)選項，包

47、括 IIR 濾波器的 Butterwotth(巴特沃思)法、Chebyshev Type I(切比雪夫 I 型)法、Chebyshev Type II(切比雪夫 II 型)法、Elliptic(橢圓濾波器)法和 FIR 濾波器的 Equiripple 法、Least-Squares(最小乘方)法、Window(窗函數(shù))法。 Filter Order(濾波器階數(shù))選項，定義濾波器的階數(shù)，包括 Specify Order(指定階數(shù))和 Minimum Order(最小階數(shù))。在 Specify Order 中填入所要設(shè)計的濾波器的階數(shù)(N 階濾波器，Specify Order=N-1)。如果選擇

48、Minimum Order，則 MATLAB 根據(jù)所選擇的濾波器類型自動使用最小階數(shù)。 Frequency Specifications 選項，可以詳細(xì)定義頻帶的各參數(shù)，包括采樣頻率和頻帶的截止頻率。它的具體選項由 Filter Type 選項和 Design Method 選項決定。例如Bandpass(帶通)濾波器需要定義 Fstop1(下阻帶截止頻率)、Fpass1(通帶下限截止頻率)、Fpass2(通帶上限截止頻率)、Fstop2(上阻帶截止頻率)，而 Lowpass(低通)濾波器只需要定義 Fstop1、Fpass1。采用窗函數(shù)設(shè)計濾波器時，由于過渡帶是由窗函數(shù)的類型和階數(shù)所決定，所

49、以只需定義通帶截止頻率，而不必定義阻帶參數(shù)。 Magnitude Specifications 選項，可以定義幅值衰減的情況。例如設(shè)計帶通濾波器時，可以定義 Wstop1(頻率 Fstop1 處的幅值衰減)、Wpass(通帶范圍內(nèi)的幅值衰減)、Wstop2(頻率 Fstop2 處的幅值衰減)。當(dāng)采用窗函數(shù)設(shè)計時，通帶截止頻率處的幅值衰減固定為 6db，精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)所以不必定義。 Window Specifications 選項，當(dāng)選取采用窗函數(shù)設(shè)計，該選項可定義，它包含了各種窗函數(shù)。圖 3-1 濾波器設(shè)計 FDATool 界面3.1.2 IIR 帶通濾波器設(shè)計實例

50、要求設(shè)計一個 10 階的帶通 Chebyshev I 濾波器，它的通帶范圍是 100 到 200Hz，采樣頻率為 1000Hz，Rp=0.5。 本例中，首先在 Filter Type 中選擇 Bandpass(帶通濾波器)；在 Design Method 選項中選擇 IIR，接著在相鄰的右則選項中選擇 Chebyshev I(切比雪夫 I 型)；指定 Filter Order項中的 Specify Order=10；由于采用的是切比雪夫設(shè)計，不必在 Options 中選擇；然后在Frequency Specifications 中選擇 Unit 為 Hz，給出采樣頻率 Fs=1000，通帶 F

51、pass1=100和 Fpass2=200；最后在 Magnitude Specifications 中選擇 Unit 為 db，Apass=0.5。設(shè)置完成后點擊 Design Filter 即可得到所設(shè)計的 IIR 濾波器。通過菜單選項 Analysis 可以在特性區(qū)看到所設(shè)計的幅頻響應(yīng)、相頻響應(yīng)、沖擊響應(yīng)和零極點配置等特性，如圖 3-2、3-3、3-4 所示。設(shè)計完成后將結(jié)果保存為 filterl.fda 文件。精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)圖 3-2 濾波器的幅頻特性圖圖 3-3 濾波器的相頻特性圖圖 3-4 濾波器的沖激響應(yīng)圖精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)3.

52、1.3 IIR 設(shè)計實例分析由圖 3-2 可知，這種濾波器在 100-200Hz 的通帶范圍內(nèi)是等波紋的，而在阻帶中是單調(diào)的，這是 Chebyshev I 濾波器的幅頻特性。由圖 3-3 可知，在 100-200Hz 的范圍內(nèi)相移較小，其曲線近似一條直線，失真較小；當(dāng)頻率超過這一范圍時，相移較大，而且其曲線是非直線的，所以失真也較大。圖 3-4 得到了 Chebyshev I 濾波器的時域沖激響應(yīng)，在 5ms-100ms有沖激響應(yīng)，超過這一范圍的沖激響應(yīng)近似為零，進(jìn)而實現(xiàn)了 Chebyshev I 帶通濾波器的設(shè)計。3.2 程序設(shè)計3.2.1 程序法介紹MATLAB 信號處理工具箱提供了幾個直

53、接設(shè)計 IIR 數(shù)字濾波器的函數(shù)，直接調(diào)用就可以設(shè)計濾波器。設(shè)計 Butterworth 濾波器用函數(shù) butter()，可以設(shè)計低通、高通、帶通和帶阻的數(shù)字和模擬濾波器，其特性是通帶內(nèi)的幅度響應(yīng)最大限度的平滑，但損失了截止頻率處的下降斜度。設(shè)計 Chebyshev I 型濾波器用函數(shù) chebyl()?？梢栽O(shè)計低通、高通、帶通和帶阻的數(shù)字和模擬 ChebyshevI 型濾披器，其通帶內(nèi)為等波紋，阻帶內(nèi)為單調(diào)。Chebyshev I 型濾波器的下降斜度比 II 型大，但其代價是通帶內(nèi)波紋較大。設(shè)計 Chebyshev II 型濾波器用函數(shù) cheby2()?？梢栽O(shè)計低通、高通、帶通和帶阻的數(shù)字

54、和模擬 Chebyshev II 型濾波器，其通帶內(nèi)為單調(diào)，阻帶內(nèi)等波紋。Chebyshev II 型濾波器的下降斜度比 I 型小，但其阻帶內(nèi)波紋較大。設(shè)計橢圓濾波器用函數(shù) ellip()，與 cheby1，cheby2 類似，可以設(shè)計低通、高通、帶通和帶阻的數(shù)字和模擬濾波器。與 Butterworth 和 chebyshev 濾波器相比，ellip 函數(shù)可以得到下降斜度更大的濾波器，得到通帶和阻帶均為等波紋。一般情況下，橢圓濾波 器能以最低的階實現(xiàn)指定的性能指標(biāo)。3.2.2 IIR 低通濾波器設(shè)計給定抽樣頻率 fS= 10KHz，要求在頻率小于 1KHz 的通帶內(nèi)，幅度特性下降小于 1dB：

55、在頻率大于 1.5KHz 的阻帶內(nèi), 衰減大于 40dB。利用模擬頻率與數(shù)字頻率之間的關(guān)系: 2/wffs3321 10 /10 100.2pw 精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)3321.5 10 /10 100.3sw采用Butterworth 模擬濾波器進(jìn)行設(shè)計,MATLAB 設(shè)計程序如下:wp= 0.2*p i;ws= 0.3*pi；n,w n = buttord (wp,ws,1,40, s)% 估計階次和截止頻率b,a = butter (n,wn,s) ;bz,az = impinVar (b,a)% 利用impinvar 可得到?jīng)_激響應(yīng)相同的數(shù)字濾波器freqz (bz

56、,az,4096,10000)；% 作出幅頻和相頻曲線圖運行結(jié)果如下:n = 6wn = 0.7229bz = 0.0000 0.0007 0.0116 0.0183 0.0046 0.0001 0az = 1.0000 -3.2943 4.8891 -4.0735 1.9859 -0.5331 0.0612在MA TLAB 7.0中運行,利用的濾波器分析函數(shù)Freqz分析所設(shè)計出的濾波器的幅頻特性和相頻特性，如圖3-5所示。設(shè)計結(jié)果滿足指標(biāo)要求。精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)圖3-5 用Butterworth低通濾波器設(shè)計的數(shù)字低通濾波器幅頻和相頻特性圖3.3 IIR 數(shù)字濾波器

57、的仿真3.3.1 IIR 帶通濾波器的仿真 對 3.1.2 設(shè)計的帶通濾波器進(jìn)行仿真。產(chǎn)生 3 個正弦成分(80Hz，150Hz 和 250Hz)的信號，利用 3.1.2 設(shè)計的帶通濾波器來去除 80Hz 和 250Hz 的正弦信號，保留 150Hz 的信號。在 MATLAB7.0 的代碼編輯窗口進(jìn)行如下操作：(1) 產(chǎn)生 3 個正弦分量信號(如圖 3-6 所示)Fs=1000；t=(1:1000)/Fs；s1=sin(2*pi*t*80)；s2=sin(2*pi*t*150)；s3=sin(2*pi*t*250)；s=s1+s2+s3；plot(t,s)；title(含有 3 個正弦分量的信

58、號)；set(gcf,color,white)；xlabel(Time(seconds)；ylabel(Time waveform)；精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)(2) 產(chǎn)生一個 10 階 IIR 帶通濾波器，通帶為 100Hz 到 200Hz，并得出其幅頻響應(yīng)（如圖3-7 所示）b,a=cheby1(10,0.5,100 200/500)；H,w=freqz(b,a,512)；plot(w*Fs/(2*pi),abs(H)；title(IIR 帶通濾波器幅頻響應(yīng))；set(gcf,color,white)；xlabel(Frequency(Hz)；ylabel(Mag.of f

59、requency response)；grid；(3) 對信號進(jìn)行濾波（濾波后信號如圖 3-8 所示）sf=filter(b,a,s)；plot(t,sf)；title(濾波后的信號波形)；set(gcf,color,white)；xlabel(Time(seconds);ylabel(Time waveform)；axis(0 1 -1 1)；(4) 繪出信號濾波前、后的幅頻圖（如圖 3-9 所示）S=fft(s,512)；SF=fft(sf,512)； w=(0:255)/256*(Fs/2)； plot(w,abs(S(1:256),SF(1:256)；title(濾波前、后的幅頻圖)；

60、set(gcf,color,white)；xlabel(Frequency(Hz)；ylabel(Mag.of frequency response);grid；legend(濾波前的幅頻,濾波后的幅頻)；精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)(5) 繪制濾波前的 150Hz 信號波形和濾波器好保留下的 150Hz 信號波形進(jìn)行比較（如圖3-10 所示）subplot(211)plot(t,s2);title(濾波前的 150Hz 信號的波形);set(gcf,color,white);xlabel(Time(seconds);ylabel(Time waveform);axis(0 1 -1 1);ho