畢業(yè)設(shè)計基于MATLAB的數(shù)字濾波器的設(shè)計與仿真分析學(xué)院名稱新華學(xué)院專業(yè)名稱電子信息工程年級2007級學(xué)生學(xué)號學(xué)生姓名指導(dǎo)教師申請日期2011年5月10日摘要本文研究了IIR數(shù)字濾波器和FIR數(shù)字濾波器在MATLAB環(huán)境下的設(shè)計方法及仿真分析。在無限脈沖響應(yīng)（IIR）數(shù)字濾波器設(shè)計中，先采用脈沖響應(yīng)不變法及雙線性Z變化法設(shè)計數(shù)字濾波器，最后進行濾波器的頻帶轉(zhuǎn)換，即進行模擬濾波器的設(shè)計，然后進行模擬數(shù)字濾波器轉(zhuǎn)換。在有限脈沖響應(yīng)（FIR）數(shù)字濾波器設(shè)計中，研究了FIR線性相位濾波器的特點和用窗函數(shù)法設(shè)計FIR濾波器兩個問題。這兩類濾波器全部設(shè)計過程都是由理論分析、編程設(shè)計、具體實現(xiàn)的步驟進行的。數(shù)字濾波器的單位沖激響應(yīng)是有限長的，可以用一個因果系統(tǒng)來實現(xiàn)，因而數(shù)字濾波器可以做成即是因果又是穩(wěn)定的系統(tǒng)。關(guān)鍵詞數(shù)字濾波器；MATLAB；無限脈沖響應(yīng)；有限脈沖響應(yīng)ABSTRACTINTHISPAPER,DESIGNSOFTHEINFINITEIMPULSERESPONSEDIGITALFILTERIIRANDFINITEIMPULSERESPONSEDIGITALFILTERFIRUNDERMATLABARESTUDIEDTHEDESIGNOFIIRFILTERSCANBEACHIEVEDTHROUGHTHREESTEPSFIRSTLY,THEDESIGNOFANALOGLOWPASSFILTERSECONDLY,ITISANALOGTODIGITALFILTERCONVERSIONLASTLY,ITISTHECONVERSIONOFFILTERFREQUENCYBANDINDESIGNOFFIRFILTERS,TWOQUESTIONSAREDISCUSSEDTHECHARACTERISTICSOFFIRLINEARPHASEFILTERANDREASONINGOFRELATEDFORMULASTHEOTHERISABOUTTHEDESIGNOFTHEFIRFILTERSBYMEANSOFWINDOWFUNCTIONSTHEDESIGNOFFIRANDIIRFOLLOWSTHEPROCEDURESOFTHEORETICALANALYSIS,PROGRAMMINGDESIGNANDREALIZATIONTHEUNITIMPULSERESPONSEOFDIGITALFILTERISLIMITED,CANBEACHIEVEDBYACAUSALSYSTEM,SODIGITALFILTERCANMAKEITTHATISCAUSALANDSTABLESYSTEMKEYWORDSDIGITALFILTERMATLABIIRFIR目錄第1章緒論111課題的背景及研究意義112國內(nèi)外的研究綜述113研究內(nèi)容和預(yù)期目標(biāo)214研究方法3第2章數(shù)字濾波器的基本概念421數(shù)字濾波器的定義422數(shù)字濾波器的分類423IIR數(shù)字濾波器設(shè)計方法4231模擬濾波器設(shè)計4232模擬數(shù)字濾波器設(shè)計5233頻帶變換724FIR數(shù)字濾波器設(shè)計方法8241線性相位條件8242線性相位FIR濾波器幅度特性的特點9243線性相位FIR濾波器零點分布特點10245各種窗函數(shù)的實現(xiàn)與比較1225本章小結(jié)12第3章基于MATLAB的IIR數(shù)字濾波器設(shè)計與仿真分析1331基于MATLAB的模擬濾波器設(shè)計與仿真分析13311巴特沃斯低通濾波器設(shè)計與仿真分析13312切比雪夫低通濾波器設(shè)計與仿真14313橢圓低通濾波器設(shè)計與仿真1732基于MATLAB的模擬數(shù)字濾波器變換仿真19321脈沖響應(yīng)不變法19322雙線性Z變換法2133基于MATLAB的頻帶變換仿真23331模擬低通濾波器轉(zhuǎn)換成數(shù)字高通濾波器23332模擬低通濾波器轉(zhuǎn)換成數(shù)字帶阻濾波器25333模擬低通濾波器轉(zhuǎn)換成數(shù)字帶通濾波器26第4章基于MATLAB的FIR數(shù)字濾波器設(shè)計與仿真分析2841基于MATLAB的窗函數(shù)的FIR數(shù)字濾波器設(shè)計與仿真分析28411海明窗設(shè)計與仿真數(shù)字低通濾波器28412漢寧窗設(shè)計與仿真數(shù)字高通濾波器29413漢寧窗設(shè)計與仿真數(shù)字帶阻濾波器33414布拉克曼窗及三角窗設(shè)計與仿真數(shù)字帶通濾波器35第5章工程實例4051IIR數(shù)字濾波器工程實例4052FIR數(shù)字濾波器工程實例41第6章結(jié)論43附錄44參考文獻47致謝48第1章緒論11課題的背景及研究意義信號與信號處理是信息科學(xué)中近幾十年來發(fā)展最為迅速的學(xué)科之一。長期以來，信號處理技術(shù)一直用于轉(zhuǎn)換、產(chǎn)生模擬或數(shù)字信號，其中最為頻繁應(yīng)用的領(lǐng)域就是信號的濾波。數(shù)字濾波是語音、圖像處理、模式識別和譜分析等應(yīng)用中的一個基本處理部件，它可以滿足濾波器對幅度和相位特性的嚴(yán)格要求，避免模擬濾波器無法克服的電壓漂移、溫度漂移和噪聲等問題。我的畢業(yè)設(shè)計的課題是基于MATLAB的數(shù)字濾波器的設(shè)計與仿真分析，其主要目的是通過此次課程設(shè)計進一步學(xué)習(xí)和鞏固數(shù)字信號處理及其相關(guān)知識，并學(xué)會利用所學(xué)的知識能力，在設(shè)計過程中能綜合運用所學(xué)知識內(nèi)容，進一步熟悉和掌握MATLAB的使用方法對數(shù)字濾波器的原理有較深的了解為即將進入社會參加工作打下堅實的基礎(chǔ)掌握收集資料、消化資料和綜合資料的能力等等。從事電子通信業(yè)而不能熟練操作使用MATLAB電子線路設(shè)計軟件，在工作和學(xué)習(xí)中將是寸步難行的。在數(shù)學(xué)、電子、金融等行業(yè)，使用MATLAB等計算機軟件對產(chǎn)品進行設(shè)計、仿真在很早以前就已經(jīng)成為了一種趨勢，這類軟件的問世也極大地提高了設(shè)計人員在通信、電子等行業(yè)的產(chǎn)品設(shè)計質(zhì)量與效率。眾所周知，實際過程中信號傳輸都要經(jīng)過調(diào)制與解調(diào)這一過程，由于消息傳過來的原始信號即調(diào)制信號具有頻譜較低的頻譜分量，這種信號在許多信道中不宜傳輸。因而，在通信系統(tǒng)的發(fā)送端通常需要有調(diào)制過程，反之在接收端則需要有解調(diào)過程。12國內(nèi)外的研究綜述數(shù)字濾波在DSP（數(shù)字信號處理）中占有重要地位。數(shù)字濾波器按實現(xiàn)的網(wǎng)絡(luò)結(jié)構(gòu)或者從單位脈沖響應(yīng)，分為IIR無限脈沖響應(yīng)）和FIR（有限脈沖響應(yīng)）濾波器。如果IIR濾波器和FIR濾波器具有相同的性能，那么通常IIR濾波器可以用較低的階數(shù)獲得高的選擇性，執(zhí)行速度更快，所有的儲存單元更少，所有既經(jīng)濟又高效。數(shù)字濾波器精確度高、使用靈活、可靠性高，具有模擬設(shè)備所沒有的許多優(yōu)點，已廣泛地應(yīng)用于各個科學(xué)技術(shù)領(lǐng)域，例如數(shù)字電視、語音、通信、雷達、聲納、遙感、圖像、生物醫(yī)學(xué)以及許多工程應(yīng)用領(lǐng)域。隨著信息時代數(shù)字時代的到來，數(shù)字濾波技術(shù)已經(jīng)成為一門及其重要的科學(xué)和技術(shù)領(lǐng)域。以往的濾波器大多采用模擬電路技術(shù)，但是，模擬電路技術(shù)存在很多難以解決的問題，例如，模擬電路元件對溫度的敏感性，等等。而采用數(shù)字技術(shù)則避免很多類似的難題，當(dāng)然數(shù)字濾波器在其他方面也有很多突出的優(yōu)點都是模擬技術(shù)所不能及的，所以采用數(shù)字濾波器對信號進行處理是目前的發(fā)展方向。MATLAB是“矩陣實驗室”MATRIXLABORATOY的縮寫，它是由美國MATHWORKS公司于1984年正式推出的一種面向科學(xué)與工程的計算軟件，是一種以矩陣運算為基礎(chǔ)的交互式程序語言，專門針對科學(xué)、工程計算及繪圖的需求。隨著版本的不斷升級，內(nèi)容不斷擴充，功能更加強大，從而被廣泛應(yīng)用于仿真技術(shù)、自動控制和數(shù)字信號處理等領(lǐng)域。全球現(xiàn)有超過五十萬的企業(yè)用戶和上千萬的個人用戶，廣泛分布在航空航天、金融財務(wù)、機械化工、電信、教育等各個行業(yè)。通過MATLAB和相關(guān)工具箱，工程師、科研人員、數(shù)學(xué)家和教育工作者可以在統(tǒng)一的平臺下完成相應(yīng)的科學(xué)計算工作。MATLAB本身包含了600余個用于數(shù)學(xué)計算、統(tǒng)計和工程處理的函數(shù)，這樣，就可以迅速完成科學(xué)計算任務(wù)而不必進行額外的開發(fā)。業(yè)內(nèi)領(lǐng)先的工具箱算法極大的擴展了MATLAB的應(yīng)用領(lǐng)域，所以MATLAB自推出以來就受到廣泛的關(guān)注，信號處理工具箱就是其中之一，在信號處理工具箱中，MATLAB提供了濾波器分析、濾波器實現(xiàn)、FIR濾波器實現(xiàn)、IIR數(shù)字濾波器設(shè)計、IIR數(shù)字濾波器階次估計等方面的函數(shù)命令。由于我們所面對的工程問題越來越復(fù)雜，過去所依賴分析的技術(shù)已逐漸不敷使用。利用電腦來分析及解決工程問題已是當(dāng)今工程師的必要工具。使用MATLAB軟件進行科學(xué)計算，能夠極大加快科研人員進行研究開發(fā)的進度，減少在編寫程序和開發(fā)算法方面所消耗的時間和有限的經(jīng)費，從而獲得最大的效能。13研究內(nèi)容和預(yù)期目標(biāo)IIR數(shù)字濾波器技術(shù)設(shè)計最通用的方法是借助于模擬濾波器的設(shè)計方法。模擬濾波器設(shè)計已經(jīng)有了一套相當(dāng)成熟的方法，它不但有完整的設(shè)計公式，而且還有較為完整的圖表可供查詢，因此充分利用這些已有的資源將會給數(shù)字濾波器的設(shè)計帶來很大方便。濾波器設(shè)計最重要的是尋找一個穩(wěn)定、因果的系統(tǒng)函數(shù)去逼近濾波器的技術(shù)指標(biāo)，因此模擬濾波器的設(shè)計十分重要。實際中，有三種廣泛應(yīng)用的濾波器，即巴特沃斯低通濾波器、切比雪夫低通濾波器，橢圓低通濾波器。在設(shè)計了模擬低通濾波器后，就可以把它們變成數(shù)字濾波器了。這些變換均是復(fù)值映射，根據(jù)數(shù)字濾波器所保持的模擬濾波器的不同特性，研究出不同的變換技術(shù)。其中，最重要的有兩種脈沖響應(yīng)不變法（保持脈沖響應(yīng)不變，又叫沖激響應(yīng)不變法）和雙線性Z變換法（保持系統(tǒng)函數(shù)不變）。除了低通數(shù)字濾波器之外，實際中還常常需要高通、帶通、及帶阻數(shù)字濾波器。即基于濾波器轉(zhuǎn)換法的IIR數(shù)字濾波器設(shè)計直接由模擬低通濾波器轉(zhuǎn)換成數(shù)字高通、帶通或帶阻濾波器。IIR數(shù)字濾波器的設(shè)計方法是利用模擬濾波器成熟的理論及設(shè)計圖表進行的，因而保留了一些典型模擬濾波器優(yōu)良的幅度特性，但設(shè)計中只考慮到了幅度特性，沒考慮到相位特性，所設(shè)計的濾波器相位特性一般是非線性的。為得到線性相位特性，必須增加相位校正網(wǎng)絡(luò)，使濾波器設(shè)計變得復(fù)雜。而FIR濾波器在保證幅度特性滿足技術(shù)要求的同時，很容易做到有嚴(yán)格的線性相位特性。設(shè)計FIR數(shù)字濾波器的基本方法有窗函數(shù)法，頻率抽樣法，等波紋切比雪夫逼近法等，窗函數(shù)在設(shè)計FIR數(shù)字濾波器中有很重要的作用，正確的選擇窗函數(shù)可以提高所設(shè)計的數(shù)字濾波器的性能，或者在滿足設(shè)計要求的情況下，減小FIR數(shù)字濾波器的階數(shù)。窗函數(shù)法的FIR數(shù)字濾波器設(shè)計海明窗設(shè)計數(shù)字低通濾波器、漢寧窗設(shè)計數(shù)字高通濾波器、布拉克曼窗及三角窗設(shè)計數(shù)字帶通濾波器和漢寧窗設(shè)計數(shù)字帶阻濾波器。14研究方法通過數(shù)學(xué)方法對數(shù)字濾波器進行分析、建模、編寫文件，應(yīng)用MATLAB軟件進行仿真及結(jié)果分析。第2章數(shù)字濾波器的基本概念21數(shù)字濾波器的定義濾波器FILTER，是一種用來消除干擾雜訊的器件，將輸入或輸出經(jīng)過過濾而得到純凈的直流電。對特定頻率的頻點或該頻點以外的頻率進行有效濾除的電路，就是濾波器，其功能就是得到一個特定頻率或消除一個特定頻率。數(shù)字濾波器由數(shù)字乘法器、加法器和延時單元組成的一種算法或裝置。數(shù)字濾波器的功能是對輸入離散信號的數(shù)字代碼進行運算處理，以達到改變信號頻譜的目的。22數(shù)字濾波器的分類根據(jù)沖激響應(yīng)函數(shù)的時間特性，數(shù)字濾波器分為IIRINFINITEIMPULSERESPONSE數(shù)字濾波器和FIRFINITEIMPULSERESPONSE數(shù)字濾波器。IIR數(shù)字濾波器存在反饋支路，即信號流圖中存在環(huán)路，其單位沖激響應(yīng)為無限長；FIR數(shù)字濾波器不存在反饋支路，其單位沖激響應(yīng)為有限長。23IIR數(shù)字濾波器設(shè)計方法231模擬濾波器設(shè)計IIR數(shù)字濾波器設(shè)計通常是借助于模擬濾波器的設(shè)計方法進行的。模擬濾波器的設(shè)計方法已經(jīng)很成熟，不僅有完整的設(shè)計公式，還有完善的圖表和曲線供查閱，因此充分利用這些已有的資源將會給數(shù)字濾波器的設(shè)計帶來很大方便。1、巴特沃斯低通濾波器設(shè)計巴特沃斯的特點是同頻帶內(nèi)的頻率響應(yīng)曲線最為平坦，沒有起伏，在截止頻率以外，頻率響應(yīng)單調(diào)下降。在通帶中是理想的單位響應(yīng)，在阻帶中響應(yīng)為0。在截止頻率處有3DB的衰減。巴特沃斯濾波器過渡帶的陡峭程度正比于濾波器的階數(shù)。而高階巴特沃斯濾波器的頻率響應(yīng)近似于理想低通濾波器。巴特沃斯BUTTERWORTH低通濾波器的幅度平方函數(shù)用下式表示|2JH21NCJH221|式中，N表示濾波器的階數(shù)。當(dāng)0時，1；當(dāng)時，1/，|JC|J2為3DB截止頻率。C2、切比雪夫低通濾波器設(shè)計巴特沃斯濾波器的頻率特性曲線，無論在通帶和阻帶內(nèi)部都是頻率的單調(diào)減函數(shù)。因此當(dāng)通帶的邊緣能夠滿足指標(biāo)要求時，在通帶的內(nèi)部肯定會有較大富余量，造成濾波器的階次N比較高。因此，更有效地設(shè)計方法是將逼近精確度均勻地分布在整個通帶內(nèi)，或者均勻分布在整個阻帶內(nèi)，或者同時分布在兩者之內(nèi)。這樣，就可以使濾波器階數(shù)大大降低。這可通過選擇具有等波紋特性的逼近函數(shù)來達到。切比雪夫濾波器的振幅特性就具有這種等波紋特性。它有兩種形式振幅特性在通帶內(nèi)是等波紋的、在阻帶內(nèi)是單調(diào)的切比雪夫I型濾波器；振幅特性在通帶內(nèi)是單調(diào)下降、在阻帶內(nèi)是等波紋的切比雪夫II型濾波器。采用何種形式的切比雪夫濾波器取決于實際用途。切比雪夫I型濾波器幅度平方函數(shù)用表示|2JH22CJN221|式中，為小于1的正數(shù)，表示通帶內(nèi)幅度波動的程度，越大，波動幅度也越大；是的切比雪夫多項式，它的定義為XCN231|,ARCOSXNHXN而切比雪夫II型濾波器幅度平方函數(shù)用表示|2JH24CNCJ2121|式中，為小于1的正數(shù)，表示通帶內(nèi)幅度波動的程度，為截止頻率；N為濾波器C的階次。3、橢圓低通濾波器設(shè)計橢圓ELLIPTIC濾波器在通帶和阻帶內(nèi)都具有等波紋幅頻響應(yīng)特性。由于其極點位置與經(jīng)典場論中的橢圓函數(shù)有關(guān)，所以由此取名為橢圓濾波器。又因為在1931年，考爾CAUER首先提出了采樣有限零點設(shè)計的濾波器，所以其另一個通用名字為考爾CARER濾波器。橢圓低通濾波器幅度平方函數(shù)用表示|2JH25UJN221|式中是雅可比橢圓函數(shù)，是與通帶衰減有關(guān)的函數(shù)。濾波器階次N等于通帶和UN2阻帶內(nèi)最大點和最小點的和。232模擬數(shù)字濾波器設(shè)計利用模擬濾波器成熟的理論及其設(shè)計方法來設(shè)計IIR數(shù)字低通濾波器是常用的方法。設(shè)計過程是按照數(shù)字濾波器技術(shù)指標(biāo)要求設(shè)計一個過渡模擬低通濾波器，SHA再按照一定的轉(zhuǎn)換關(guān)系將轉(zhuǎn)換成數(shù)字低通濾波器的系統(tǒng)函數(shù)。由此可見，設(shè)SHAZ計的關(guān)鍵問題就是找到這種轉(zhuǎn)換關(guān)系，將S平面上的轉(zhuǎn)換成Z平面上的。為SHAZ了保證轉(zhuǎn)換后的穩(wěn)定且滿足技術(shù)指標(biāo)要求，對轉(zhuǎn)換關(guān)系提出兩點要求Z1因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。我們知道，模擬濾波器因果穩(wěn)定的條件是其系統(tǒng)函數(shù)的極點全部位于S平面的左半平面；數(shù)字濾SHA波器因果穩(wěn)定的條件時的極點全部在單位圓內(nèi)。因此，轉(zhuǎn)換關(guān)系應(yīng)使S平面的左Z半平面映射到Z平面的單位圓內(nèi)部。2數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻響特性，S平面的虛軸映射為Z平面的單位圓，相應(yīng)的頻率之間呈線性關(guān)系。將系統(tǒng)函數(shù)從S平面轉(zhuǎn)換到Z平面的方法有很多種，但工程上常用的是脈沖HA響應(yīng)不變法和雙線性變換法。1、脈沖響應(yīng)不變法設(shè)模擬濾波器的系統(tǒng)函數(shù)為，相應(yīng)的單位沖激響應(yīng)為,并將脈沖響應(yīng)SATH進行等間隔采樣，使得數(shù)字濾波器的單位抽樣響應(yīng)剛好等于的采樣值，即THN0SNSNTTNTHTTHHS（26）其中為采樣周期。ST是模擬濾波器的系統(tǒng)傳遞函數(shù)，又令是數(shù)字濾波器的系統(tǒng)傳遞函數(shù)。采HAZH樣信號的拉式變換與相應(yīng)的采樣序列Z變換的映射關(guān)系為EZST（27）所以系統(tǒng)函數(shù)和的關(guān)系為SAZH28KSASEZJKHTS1|式28的物理意義為首先將模擬濾波器的系統(tǒng)函數(shù)作周期的延拓，在經(jīng)過A式27的映射變換，映射到Z平面上，從而得到數(shù)字濾波器的系統(tǒng)函數(shù)。且模擬ZH和數(shù)字頻率滿足下列關(guān)系。經(jīng)過式27的映射，S平面的左半平面映射為Z平面的單位圓內(nèi)，因此，一個因果的和穩(wěn)定的模擬濾波器映射成因果的和穩(wěn)定的數(shù)字濾波器。經(jīng)過以上分析，按照脈沖響應(yīng)不變法，通過模擬濾波器的系統(tǒng)傳遞函數(shù)，可SA直接求得數(shù)字濾波器系統(tǒng)函數(shù),其設(shè)計具體步驟歸納如下ZH1利用可由關(guān)系式推出，將數(shù)字濾波器指標(biāo)，轉(zhuǎn)換為模擬ESTPS濾波器指標(biāo)，PS2根據(jù)指標(biāo)，來設(shè)計模擬濾波器SA3利用部分分式展開法，把展成SA29NKKAPSAH14最后把模擬極點轉(zhuǎn)換為數(shù)字極點，得到數(shù)字濾波器PKETS210NKTSKZZ112、雙線性Z變換法為了克服脈沖響應(yīng)不變法設(shè)計數(shù)字濾波器時，由于的頻率關(guān)系是根據(jù)推導(dǎo)的，所以是軸每隔便映射到單位圓上一周，引起的頻域混疊，人們EZSTJ/2找到了另一種映射關(guān)系2111ZTS此關(guān)系稱為雙線性Z變換法。233頻帶變換1、模擬低通濾波器轉(zhuǎn)換成數(shù)字高通濾波器設(shè)模擬低通濾波器的系統(tǒng)傳遞函數(shù)為，則模擬低通濾波器S平面到數(shù)字高通SHA濾波器的Z平面的變換公式為212ZTS12頻率變換公式為2132COTS得到數(shù)字高通濾波器的傳遞函數(shù)表達式為214|12ZTSASHZ2、模擬低通濾波器轉(zhuǎn)換成數(shù)字帶阻濾波器設(shè)模擬低通濾波器的系統(tǒng)傳遞函數(shù)為，則模擬低通濾波器S平面到數(shù)字帶阻SA濾波器的Z平面的變換公式為215ZS201COS2頻率變換公式為216SIN0其中2172COS0LH式中為下截止頻率，為上截止頻率。LH從而可得數(shù)字帶阻濾波器的傳遞函數(shù)表達式為218|201COS2ZSASHZ3、模擬低通濾波器轉(zhuǎn)換成數(shù)字帶通濾波器設(shè)模擬低通濾波器的系統(tǒng)傳遞函數(shù)為，則模擬低通濾波器S平面到數(shù)字帶通A濾波器的Z平面的變換公式為219ZS1COS220頻率變換公式為220SINC0其中2212COS0LH式中為下截止頻率，為上截止頻率。LH從而可得數(shù)字帶通濾波器的傳遞函數(shù)表達式為222|1COS220ZSAHZ24FIR數(shù)字濾波器設(shè)計方法穩(wěn)定和線性相位特性是FIR濾波器突出的特點。設(shè)FIR濾波器單位脈沖響應(yīng)長度為N，其系統(tǒng)函數(shù)為NHZ22310NNNZHH是的N1次多項式，它在Z平面上有N1個零點，原點是N1階重極ZH10Z點。因此，永遠穩(wěn)定。ZFIR濾波器的設(shè)計方法與IIR數(shù)字濾波器的設(shè)計方法有很大的不同。FIR濾波器的設(shè)計任務(wù)是選擇有限長度的，使傳遞函數(shù)滿足技術(shù)要求。NHEHJ241線性相位條件對于長度為N的,傳遞函數(shù)為22410NNNJJEHE225HJIJ式中，稱為幅度特性，稱為相位特性。的線性相位是指是的HEJ線性函數(shù)，即，為常數(shù)226如果滿足下式，是起始相位22700以上兩種情況都滿足群延遲是一個常數(shù)，即D滿足226為第一類線性相位；滿足227為第二類線性相位。第一類線性相位特性是是實序列且對偶對稱，即。第二類線性相位特性是NH2/1N1NNH是實序列且對奇對稱，即。具有線性相位的FIR數(shù)字濾波/NH器既有恒定的延遲群，又有恒定的相延遲。242線性相位FIR濾波器幅度特性的特點的長度N取奇數(shù)還是偶數(shù)影響的特性。所以，這兩類線性相位，有四種NHH情況討論其幅度特性特點。1為偶對稱，且N為奇數(shù)型濾波器當(dāng)為偶對稱，且N為奇數(shù)時，濾波器的幅頻函數(shù)可以表示為（228）1/20COSNNAN其中，。此時，對呈102AH1,ANHH0,2偶對稱。2為偶對稱，且N為偶數(shù)型濾波器N當(dāng)為偶對稱，且N為偶數(shù)時，濾波器的幅頻函數(shù)可以表示為H（229）/211COS2NHBN其中。此時對呈奇對稱。但是，故2,1,BNHN，0H高通濾波器不能用這種方法實現(xiàn)。3為奇對稱，且N為奇數(shù)型濾波器當(dāng)HN為奇對稱，且為奇數(shù)時，濾波器的幅頻函數(shù)可以表示為（230）1/2SINNNHC其中。此時對呈奇對稱。但是當(dāng)2,12,CNHN02，時，所以低通、高通濾波器不能采用這種形式。02，0H4為奇對稱，且N為偶數(shù)型濾波器NH當(dāng)為奇對稱，且N為偶數(shù)時，濾波器的幅頻函數(shù)可表示為（231）/211SIN2NND其中。此時對呈奇對稱，對呈偶對2,1,DNHNH0，稱。但是當(dāng)時，所以低通濾波器不能采用這種形式。0,0243線性相位FIR濾波器零點分布特點第一類和第二類線性相位的系統(tǒng)函數(shù)綜合起來表示為232ZHZN1上式表明，如是的零點，其倒數(shù)也必然是其零點；又因為是實序列，ZIHI1NH的零點必定共軛成對，因此和也是其零點。這樣，線性相位FIR濾波器零點ZHZII分布特點是零點必須是互為倒數(shù)的共軛對，確定其中一個，另外三個零點也就確定了。244常用窗函數(shù)介紹1、矩形窗RECTANGULARWINDOW定義為233NRN頻率響應(yīng)為23412SINNJJRWE主瓣寬度為，第一副瓣比主瓣低13DB。EJR/4實現(xiàn)函數(shù)為RECTWIN，調(diào)用格式為WRECTWINNN是窗函數(shù)的長度，返回值W是一個N階的向量，它由窗函數(shù)的值組成。2、三角形窗TRIANGULARWINDOW定義為23512,01NNNR其頻率響應(yīng)為2362121SIN4NJJRWE其主瓣寬度為，第一副瓣比主瓣低26DB。/8實現(xiàn)函數(shù)為TRIANG，其調(diào)用格式為WTRIANGN3、漢寧窗HANNINGWINDOW漢寧窗函數(shù)又稱升余弦函數(shù)，其時域表達式為237NRNNR12COS150其頻率響應(yīng)為238EWJRJR21其最大旁瓣值比主瓣值低31DB，但是主瓣寬度比矩形窗函數(shù)的主瓣寬度增加了一倍，為。N/8實現(xiàn)函數(shù)為HANN，其調(diào)用格式為WHANNN4、海明窗HAMMINGWINDOW海明窗函數(shù)是一種改進的升余弦函數(shù)，其時域表達式為239NRNNNR12COS4605其幅度頻率特性為240121340WWRRRHM其主瓣寬度為，能量更集中在主瓣之中，主瓣的能量約占9996，第一旁瓣的N/8峰值比主瓣小40DB。實現(xiàn)函數(shù)為HAMMING，其調(diào)用格式為WHAMMINGN5、布拉克曼窗BLACKMANWINDOW定義為241110425COS08COS4KKWKNN其中。其幅度頻率特性為1,2KN241111RRRRRWWWN（242）為矩形窗函數(shù)的幅度頻率特性函數(shù)。布拉克曼的主瓣寬度是矩形窗主瓣寬度的R三倍，為。它的最大旁瓣值比主瓣值低57DB。N/12實現(xiàn)函數(shù)為BLACKMAN，其調(diào)用格式為WBLACKMANN245各種窗函數(shù)的實現(xiàn)與比較在設(shè)計FIR數(shù)字濾波器時，要選擇合適的窗函數(shù)來進行處理，調(diào)整窗函數(shù)形狀可使帶內(nèi)波動減少以及加大阻帶衰減；調(diào)整窗口長度N可以有效地控制過渡帶寬度。在設(shè)計時應(yīng)當(dāng)根據(jù)實際需要和技術(shù)要求來選取窗函數(shù)進行設(shè)計。如果找到的窗函數(shù)形狀，使其譜函數(shù)的主瓣包含更多的能量，則相應(yīng)的旁瓣幅度就能減小，而旁瓣的減小可使通帶阻帶波動減小，加大阻帶衰減，但這樣總是以加寬過渡帶為代價的。25本章小結(jié)這一章，IIR濾波器的設(shè)計步驟分為三步，即模擬低通濾波器設(shè)計，模擬數(shù)字濾波器變換，濾波器的頻帶變換。模擬低通濾波器的設(shè)計中，主要研究三種設(shè)計方法在模擬數(shù)字濾波器變換中，討論了兩種變換方法，為脈沖響應(yīng)不變法和雙線性Z變換法；在頻帶變換的實現(xiàn)中，以巴特沃斯濾波器為例進行分析設(shè)計。FIR濾波器的設(shè)計包括三部分線性相位FIR數(shù)字濾波器的條件和特點，常用窗函數(shù)及其MATLAB仿真，基于窗函數(shù)的FIR數(shù)字濾波器設(shè)計及其MATLAB仿真。在設(shè)計任意幅度頻率特性濾波器的同時，保證精確、嚴(yán)格的線性相位特性。FIR數(shù)字濾波器的單位沖激響應(yīng)是有限長的，可以用一個因果系統(tǒng)來實現(xiàn)，因而FIR數(shù)字濾波器可以做成即是因果又是穩(wěn)定的系統(tǒng)。在常用窗函數(shù)及MATLAB仿真中，討論了五種窗函數(shù)及其各自特性和相互比較，并結(jié)合MATLAB進行仿真；在基于窗函數(shù)的FIR數(shù)字濾波器設(shè)計中，根據(jù)設(shè)計的不同技術(shù)要求，選擇了不同的窗函數(shù)，并利用MATLAB提供的相關(guān)函數(shù)來設(shè)計實現(xiàn)FIR數(shù)字濾波器。正確地選擇窗函數(shù)可以提高所設(shè)計的數(shù)字濾波器的性能，減小FIR數(shù)字濾波器的階次。同時它也是數(shù)字濾波器設(shè)計中最簡單的方法。第3章基于MATLAB的IIR數(shù)字濾波器設(shè)計與仿真分析31基于MATLAB的模擬濾波器設(shè)計與仿真分析311巴特沃斯低通濾波器設(shè)計與仿真分析MATLAB信號處理工具箱函數(shù)BUTTAP是巴特沃斯低通模擬濾波器的設(shè)計函數(shù)。其調(diào)用的格式為Z,P,KBUTTAPN，其中，Z表示零點,P表示極點,K表示增益，N表示階次。取N階巴特沃斯濾波器（N1,4,10,22），進行比較分析N00013FORI14SWITCHICASE1N1CASE2N4CASE3N10CASE4N22ENDZ,P,KBUTTAPN設(shè)計巴特沃斯低通濾波器B,AZP2TFZ,P,K零極點增益模型轉(zhuǎn)換為傳遞函數(shù)模型H,WFREQSB,A,N求解模擬濾波器頻率響應(yīng)MAGH2ABSH2取模值函數(shù)HOLDON控制是否保持當(dāng)前圖形PLOTW,MAGH2畫二維線性圖AXIS0301控制坐標(biāo)軸比例和外觀ENDXLABELW/WCYLABEL|HJW|2TITLE巴特沃斯低通模擬濾波器張晉川GRIDON運行結(jié)果如圖31所示圖31巴特沃斯濾波器的幅頻特性由圖31可知，對N1階濾波器和N22階濾波器的平方幅度響應(yīng)進行比較后，證明了高階巴特沃斯濾波器有著更好的幅度特性，更接近理想低通濾波器。巴特沃斯濾波器擁有平滑的頻率響應(yīng)，在截止頻率以外，頻率響應(yīng)單調(diào)下降。其過渡帶的陡峭程度正比于濾波器的階數(shù)，高階巴特沃斯濾波器的頻率響應(yīng)近似于理想低通濾波器。312切比雪夫低通濾波器設(shè)計與仿真MATLAB信號處理工具箱函數(shù)CHEB1AP及CHEB2AP是切比雪夫濾波器設(shè)計函數(shù)。其調(diào)用格式如下Z,P,KCHEB1APN,RPZ,P,KCHEB2APN,RS其中，Z表示零點,P表示極點,K表示增益，N表示階次，RP為通帶波紋DB,RS為阻帶波紋DB。取N階切比雪夫I型濾波器（N5,6），進行比較分析切比雪夫I型低通模擬濾波器N00012FORI12SWITCHICASE1N5CASE2N6ENDRP1通帶紋波系數(shù)為1DBZ,P,KCHEB1APN,RP設(shè)計切比雪夫I型低通濾波器B,AZP2TFZ,P,KH,WFREQSB,A,N求解模擬濾波器頻率響應(yīng)MAGH2ABSH2POSPLOT12NUM2STRI輸出圖形SUBPLOTPOSPLOT創(chuàng)建子圖圖區(qū)PLOTW,MAGH2AXIS0201XLABELW/WC張晉川YLABEL|HJW|2TITLENNUM2STRNEND運行結(jié)果如圖32所示圖32切比雪夫I型濾波器的幅頻特性取N階切比雪夫II型濾波器（N5,6），進行比較分析切比雪夫II型低通模擬濾波器N00012FORI12SWITCHICASE1N5CASE2N6ENDRS10阻帶文波系數(shù)為10DBZ,P,KCHEB2APN,RS函數(shù)CHEB2設(shè)計切比雪夫II型低通濾波器B,AZP2TFZ,P,KH,WFREQSB,A,NMAGH2ABSH2POSPLOT12NUM2STRISUBPLOTPOSPLOTPLOTW,MAGH2AXIS02011XLABELW/WC張晉川YLABEL|HJW|2TITLENNUM2STRNEND運行結(jié)果如圖33所示圖33切比雪夫II型濾波器的幅頻特性由圖32與圖33分別比較可得出，切比雪夫濾波器有效地將技術(shù)指標(biāo)的精確度均勻分布，切比雪夫I型濾波器的振幅特性在通帶內(nèi)是等波紋的、在阻帶內(nèi)是單調(diào)的；切比雪夫II型濾波器的振幅特性在通帶內(nèi)是單調(diào)的、在阻帶內(nèi)是等波紋的。313橢圓低通濾波器設(shè)計與仿真MATLAB信號處理工具箱函數(shù)ELLIPAP是橢圓低通模擬濾波器設(shè)計函數(shù)。其調(diào)用格式如下Z,P,KELLIPAPN,RP,RS，其中，Z表示零點,P表示極點,K表示增益，N表示階次,RP為通帶波紋DB,RS為阻帶波紋DB。取N階橢圓濾波器（N2,3,4,5），進行比較分析N00012FORI14SWITCHICASE1N2CASE2N3CASE3N4CASE4N5ENDRP1通帶紋波系數(shù)為1DBRS10阻帶紋波系數(shù)為10DBZ,P,KELLIPAPN,RP,RS設(shè)計橢圓低通濾波器B,AZP2TFZ,P,KH,WFREQSB,A,NMAGH2ABSH2POSPLOT22NUM2STRISUBPLOTPOSPLOTPLOTW,MAGH2AXIS02011XLABELW/WC張晉川YLABEL|HJW|2TITLENNUM2STRNGRIDONEND運行結(jié)果如圖34所示圖34橢圓濾波器的幅頻特性由圖34可得出，橢圓濾波器在通帶及阻帶內(nèi)均成等波紋振蕩，并具有最窄的過渡帶，其性質(zhì)比較優(yōu)越，但橢圓濾波器的轉(zhuǎn)移函數(shù)不但有極點，而且在軸上還有零J點，所以橢圓濾波器的設(shè)計比較復(fù)雜。32基于MATLAB的模擬數(shù)字濾波器變換仿真321脈沖響應(yīng)不變法在MATLAB環(huán)境下用函數(shù)實現(xiàn)脈沖響應(yīng)不變法設(shè)計數(shù)字低通濾波器。其函數(shù)為B,AIMPINVARC,D,T其中，B表示數(shù)字濾波器自變量為的分子多項式，A表示數(shù)字濾波器自變量為Z1的分母多項式，C表示模擬濾波器自變量為S的分子多項式,D表示模擬濾波器自變Z1量為S的分母多項式，T表示采樣變換參數(shù)。根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個IIR數(shù)字低通濾波器，02WP4S3DBRP20DAS并繪出所設(shè)計的濾波器的幅度響應(yīng)。WP02PIWS04PIRP1AS15T1采樣周期RIP10RP/20ATN10AS/20OMGPWPTOMGSWSTN,OMGCBUTTORDOMGP,OMGS,RP,AS,S選取模擬濾波器的階數(shù)CS,DSBUTTERN,OMGC,S設(shè)計出所需的模擬低通濾波器B,AIMPINVARCS,DS,T應(yīng)用脈沖響應(yīng)不變法進行轉(zhuǎn)換DB,MAG,PHA,GRD,WFREQZ_MB,A求得相對，絕對頻響及相位、群延遲響應(yīng)SUBPLOT2,2,1PLOTW/PI,MAGTITLE幅頻特性張晉川XLABELW/PIYLABEL|HJW|AXIS0,1,0,11SETGCA,XTICKMODE,MANUAL,XTICK,00203051SETGCA,YTICKMODE,MANUAL,YTICK,0ATNRIP1GRIDSUBPLOT2,2,2PLOTW/PI,DBTITLE幅頻特性DB張晉川XLABELW/PIYLABELDBAXIS0,1,40,5SETGCA,XTICKMODE,MANUAL,XTICK,00203051SETGCA,YTICKMODE,MANUAL,YTICK,40ASRP0GRIDSUBPLOT2,2,3PLOTW/PI,PHA/PITITLE相頻特性張晉川XLABELW/PIYLABELPHA/PIAXIS0,1,1,1SETGCA,XTICKMODE,MANUAL,XTICK,00203051GRIDSUBPLOT2,2,4PLOTW/PI,GRDTITLE群延遲張晉川XLABELW/PIYLABELSAMPLEAXIS0,1,0,12SETGCA,XTICKMODE,MANUAL,XTICK,00203051GRIDON運行結(jié)果如圖35所示圖35脈沖響應(yīng)不變法設(shè)計數(shù)字低通濾波器的特性脈沖響應(yīng)不變法的優(yōu)點是頻率坐標(biāo)變換是線性的，即，如不考慮頻率混疊現(xiàn)象，用這種方法設(shè)計數(shù)字濾波器會很好的重現(xiàn)原模擬濾波器的頻率響應(yīng)。另外一個優(yōu)點是數(shù)字濾波器的單位脈沖響應(yīng)完全模仿模擬濾波器的單位沖激響應(yīng)，時域逼近好。但其也具有很大的缺點，若抽樣頻率不高或其它原因?qū)a(chǎn)生混疊失真，不能重現(xiàn)原模擬濾波器頻率響應(yīng)。所以，脈沖響應(yīng)不變法不適合高通、帶阻濾波器的設(shè)計，適合低通、帶通濾波器設(shè)計。322雙線性Z變換法雙線性Z變換可通過BILINEAR函數(shù)實現(xiàn)。其調(diào)用格式如下BZ,AZBILINEARB,A,FS其中B,A為模擬濾波器的傳遞函數(shù)的分子分母多項式系數(shù)分量，而BZ,AZ為SHA數(shù)字濾波器的傳遞函數(shù)HZ的分子分母多項式的系數(shù)分量。根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個IIR數(shù)字低通濾波器通帶邊界頻率FP100HZ,阻帶邊界頻率FS300HZ,采樣頻率FS1KHZ,通帶最大衰減RP3DB，阻帶最小衰減RS20DB并繪出所設(shè)計的濾波器的幅度響應(yīng)。FP100FS300FS1000RP3RS20WP2PIFP/FSWS2PIFS/FSFSFS/FSWAPTANWP/2WASTANWS/2N,WNBUTTORDWAP,WAS,RP,RS,SZ,P,KBUTTAPNBP,APZP2TFZ,P,KBS,ASLP2LPBP,AP,WAPBZ,AZBILINEARBS,AS,FS/2H,WFREQZBZ,AZ,256,FS1000PLOTW,ABSHTITLE雙線性Z變換設(shè)計數(shù)字低通濾波器張晉川GRIDON運行結(jié)果FS1N2WN04363BP001AP100001414210000BS01056AS100000459501056BZ006750134900675AZ100001143004128運行結(jié)果如圖36所示圖36雙線性Z變換設(shè)計數(shù)字低通濾波器的特性與脈沖響應(yīng)不變法相比，雙線性變換的主要優(yōu)點靠頻率的嚴(yán)重非線性關(guān)系得到S平面與Z平面的單值一一對應(yīng)關(guān)系，整個軸單值對應(yīng)于單位圓一周，這個關(guān)系就J是式所表示的，其中和為非線性關(guān)系。如圖圖中看到，在零頻率附近，接近于線性關(guān)系，進一步增加時，增長變得緩慢，終止于折疊頻率處，所以雙線性變換不會出現(xiàn)由于高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象。33基于MATLAB的頻帶變換仿真331模擬低通濾波器轉(zhuǎn)換成數(shù)字高通濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個IIR數(shù)字高通濾波器，05WP3S1DBAP5DS并繪出所設(shè)計的濾波器的幅度響應(yīng)。WP05PIWS03PIAP1AS15N,WNBUTTORDWP/PI,WS/PI,AP,AS計算巴特沃斯濾波器階次和截至頻率B,ABUTTERN,WN,HIGH頻率變換法設(shè)計巴特沃斯高通濾波器B0,B,ADIR2CASB,A數(shù)字高通濾波器級聯(lián)型DB,MAG,PHA,GRD,WFREQZ_MB,A數(shù)字濾波器響應(yīng)SUBPLOT211PLOTW/PI,MAGTITLE高通數(shù)字濾波器幅頻響|HEJOMEGA|張晉川SUBPLOT212PLOTW/PI,DBTITLE高通數(shù)字濾波器幅頻響DB張晉川運行結(jié)果B001482B100002000010000100002000010000A100002000010000100002000010000運行結(jié)果如圖37所示圖37巴特沃斯模擬濾波器設(shè)計數(shù)字高通濾波器的幅頻響應(yīng)曲線由此結(jié)果可知，所設(shè)計的高通濾波器滿足設(shè)計要求。332模擬低通濾波器轉(zhuǎn)換成數(shù)字帶阻濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個IIR數(shù)字帶阻濾波器，,073WP,082S1DBAP30DS并繪出所設(shè)計的濾波器的幅度響應(yīng)。WP03PI,07PIWS02PI,08PIAP1AS30N,WNBUTTORDWP/PI,WS/PI,AP,AS計算巴特沃斯濾波器階次和截至頻率B,ABUTTERN,WN,STOP頻率變換法設(shè)計巴特沃斯帶阻濾波器B0,B,ADIR2CASB,A數(shù)字帶阻濾波器級聯(lián)型DB,MAG,PHA,GRD,WFREQZ_MB,A數(shù)字濾波器響應(yīng)SUBPLOT211PLOTW/PI,MAGTITLE數(shù)字帶阻濾波器幅頻響|HEJOMEGA|張晉川SUBPLOT212PLOTW/PI,DBTITLE數(shù)字帶阻濾波器幅頻響DB張晉川運行結(jié)果B000294B100000009309978100000007510059100000004109915100000000010096100000004109915100000007510059100000009309978A100000009309978100000007510059100000004109915100000000010096100000004109915100000007510059100000009309978運行結(jié)果如圖38所示圖38基于巴特沃斯模擬濾波器設(shè)計的數(shù)字帶阻濾波器的幅頻響應(yīng)曲線由此結(jié)果可知，所設(shè)計的帶阻波器滿足設(shè)計要求。333模擬低通濾波器轉(zhuǎn)換成數(shù)字帶通濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個IIR數(shù)字帶通濾波器，,，,043WP061S5DBAP20DS并繪出所設(shè)計的濾波器的幅度響應(yīng)。WP03PI,04PIWS01PI,06PIAP5AS20N,WNBUTTORDWP/PI,WS/PI,AP,AS計算巴特沃斯濾波器階次和截至頻率B,ABUTTERN,WN,BANDPASS頻率變換法設(shè)計巴特沃斯帶通濾波器B0,B,ADIR2CASB,A數(shù)字帶通濾波器級聯(lián)型DB,MAG,PHA,GRD,WFREQZ_MB,A數(shù)字濾波器響應(yīng)SUBPLOT211PLOTW/PI,MAGTITLE帶通數(shù)字濾波器幅頻響|HEJOMEGA|張晉川SUBPLOT212PLOTW/PI,DBTITLE帶通數(shù)字濾波器幅頻響DB張晉川運行結(jié)果B000460B100002000010000100002000010000A100002000010000100002000010000運行結(jié)果如圖39所示圖39基于巴特沃斯模擬濾波器設(shè)計的數(shù)字帶通濾波器的幅頻響應(yīng)曲線由此結(jié)果可知，所設(shè)計的帶通濾波器滿足設(shè)計要求。第4章基于MATLAB的FIR數(shù)字濾波器設(shè)計與仿真分析41基于MATLAB的窗函數(shù)的FIR數(shù)字濾波器設(shè)計與仿真分析FIR數(shù)字濾波器的基本設(shè)計方法有窗函數(shù)法，頻率抽樣法，等波紋切比雪夫逼近法等，在此主要研究基于窗函數(shù)法設(shè)計FIR數(shù)字濾波器及其MATLAB實現(xiàn)的問題。411海明窗設(shè)計與仿真數(shù)字低通濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個FIR數(shù)字低通濾波器，02WP04S25DBAP0S選擇一個適當(dāng)?shù)拇昂瘮?shù)，確定單位沖激響應(yīng)，繪出所設(shè)計的濾波器的幅度響應(yīng)。WP02PIWS04PITR_WIDTHWSWP過渡帶寬度NCEIL66PI/TR_WIDTH1濾波器長度N01N1WCWSWP/2理想低通濾波器的截至頻率HDIDEAL_LPWC,N理想低通濾波器的單位沖激響應(yīng)W_HAMHAMMINGN海明窗HHDW_HAM截取得到實際的單位脈沖響應(yīng)DB,MAG,PHA,WFREQZ_M2H,1計算實際濾波器的幅度響應(yīng)DELTA_W2PI/1000APMINDB11WP/DELTA_W1實際通帶紋波ASROUNDMAXDBWS/DELTA_W11501實際阻帶紋波SUBPLOT221STEMN,HD火柴棍圖TITLE理想單位脈沖響應(yīng)HDN張晉川SUBPLOT222STEMN,W_HAMTITLE海明窗WN張晉川SUBPLOT223STEMN,HTITLE實際單位脈沖響應(yīng)HN張晉川SUBPLOT224PLOTW/PI,DBTITLE幅度響應(yīng)DB張晉川AXIS0,1,100,10運行結(jié)果N34AP00477AS52運行結(jié)果如圖41所示圖41II型FIR數(shù)字低通濾波器脈沖響應(yīng)與幅度響應(yīng)由此結(jié)果可知，所設(shè)計的低通濾波器為型濾波器，它的通帶紋波和阻帶紋波均滿足設(shè)計要求。412漢寧窗設(shè)計與仿真數(shù)字高通濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個FIR數(shù)字高通濾波器，06WP04S25DBAP0S選擇一個適當(dāng)?shù)拇昂瘮?shù)，確定單位沖激響應(yīng)，繪出所設(shè)計的濾波器的幅度響應(yīng)。WP06PIWS04PITR_WIDTHWPWS過渡帶寬度NCEIL62PI/TR_WIDTH濾波器長度N01N1WCWSWP/2理想高通濾波器的截至頻率HDIDEAL_HP1WC,N理想高通濾波器的單位沖激響應(yīng)W_HANHANNINGN漢寧窗HHDW_HAN截取得到實際的單位脈沖響應(yīng)DB,MAG,PHA,WFREQZ_M2H,1計算實際濾波器的幅度響應(yīng)DELTA_W2PI/1000APMINDBWP/DELTA_W11501實際通帶紋波ASROUNDMAXDB11WS/DELTA_W1實際阻帶紋波SUBPLOT221STEMN,HD火柴棍圖TITLE理想單位脈沖響應(yīng)HDN張晉川SUBPLOT222STEMN,W_HANTITLE漢寧窗WN張晉川SUBPLOT223STEMN,HTITLE實際單位脈沖響應(yīng)HN張晉川SUBPLOT224PLOTW/PI,DBTITLEFIR數(shù)字高通濾波器幅度響應(yīng)DB張晉川AXIS0,1,100,10運行結(jié)果N31AP00887AS44運行結(jié)果如圖42所示圖42I型FIR數(shù)字高通濾波器脈沖響應(yīng)與幅度響應(yīng)由此結(jié)果可知，所設(shè)計的高通濾波器為型濾波器，它的通帶紋波和阻帶紋波均滿足設(shè)計要求。根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個具有/2相移的FIR數(shù)字高通濾波器，06WP04S25DBAP0S選擇一個適當(dāng)?shù)拇昂瘮?shù)，確定單位沖激響應(yīng)，繪出所設(shè)計的濾波器的幅度響應(yīng)。WP06PIWS04PITR_WIDTHWPWS過渡帶寬度NCEIL62PI/TR_WIDTH1濾波器長度N01N1WCWSWP/2理想低通濾波器的截止頻率HDIDEAL_HP2WC,N理想低通濾波器的單位沖激響應(yīng)W_HANHANNINGN漢寧窗HHDW_HAN截取得到實際的單位脈沖響應(yīng)DB,MAG,PHA,WFREQZ_M2H,1計算實際濾波器的幅度響應(yīng)DELTA_W2PI/1000APMINDBWP/DELTA_W11501實際通帶紋波ASROUNDMAXDB11WS/DELTA_W1實際阻帶紋波SUBPLOT221STEMN,HDTITLE理想單位脈沖響應(yīng)HDN張晉川SUBPLOT222STEMN,W_HANTITLE漢寧窗WN張晉川SUBPLOT223STEMN,HTITLE實際單位脈沖響應(yīng)HDN張晉川SUBPLOT224PLOTW/PI,DBTITLE幅度響應(yīng)DB張晉川AXIS0,1,100,10運行結(jié)果N32AP00711AS44運行結(jié)果如圖43所示圖43IV型FIR數(shù)字高通濾波器脈沖響應(yīng)與幅度響應(yīng)由此結(jié)果可知，所設(shè)計的高通濾波器為型濾波器，它的通帶紋波和阻帶紋波均滿足設(shè)計要求。413漢寧窗設(shè)計與仿真數(shù)字帶阻濾波器根據(jù)下列技術(shù)指標(biāo)，設(shè)計一個FIR數(shù)字帶阻濾波器，,02WPL08PH04WSL06SH，1DBALDBLA選擇一個適當(dāng)?shù)拇昂瘮?shù)，確定單位沖激響應(yīng)，繪出所設(shè)計的濾波器的幅度響應(yīng)。WPL02PIWPH08PIWSL04PIWSH06PITR_WIDTHMINWSLWPL,WPHWSH過渡帶寬度NCEIL62PI/TR_WIDTH濾波器長度N01N1WCLWSLWPL/2理想帶阻濾波器的下截至頻率WCHWSHWPH/2理想帶阻濾波器的上截至頻率HDIDEAL_BSWCL,WCH,N理想帶阻濾波器的單位沖激響應(yīng)W_HANHANNINGN漢寧窗HHDW_HAN截取得到實際的單位脈沖響應(yīng)DB,MAG,PHA,WFREQZ_M2H,1計