

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、<p><b> 課程設計任務書</b></p><p> 學生姓名: xxx 專業(yè)班級: 信息xxx班 </p><p> 指導教師: xxx 工作單位:信息工程學院 </p><p> 題 目:MATLAB課程設計</p><p> 要求完成的主要任務:
2、 </p><p> 1. 在數(shù)字信號處理平臺上(PC機﹑MATLAB仿真軟件系統(tǒng))進行軟件仿真設計,并進行調試和數(shù)據(jù)分析。</p><p> 2. 利用MATLAB仿真軟件系統(tǒng)結合頻率取樣法設計一個數(shù)字高通FIR濾波器。</p><p><b> 課程設計的目的:</b></p><p><b> 1
3、.理論目的 </b></p><p> 課程設計的目的之一是為了鞏固課堂理論學習,并能用所學理論知識正確分析信號處理的基本問題和解釋信號處理的基本現(xiàn)象。</p><p><b> 2.實踐目的</b></p><p> 課程設計的目的之二是通過設計具體的圖像信號變換掌握圖像和信號處理的方法和步驟。</p><
4、;p><b> 時間安排:</b></p><p> 指導教師簽名: 年 月 日</p><p> 系主任(或責任教師)簽字: 年 月 日</p><p><b> 目 錄</b></p
5、><p><b> 1 引言1</b></p><p> 1.1MATLAB介紹1</p><p> 1.2MATLAB信號處理工具箱函數(shù)介紹1</p><p> 1.3濾波器的介紹2</p><p> 2 FIR數(shù)字濾波器設計原理3</p><p> 3
6、 FIR數(shù)字濾波器設計方法4</p><p><b> 3.1窗函數(shù)法4</b></p><p> 3.2頻率取樣法5</p><p> 4 頻率采樣法實際FIR高通濾波器7</p><p><b> 4.1設計原理7</b></p><p><b&g
7、t; 4.2設計步驟9</b></p><p> 5 MATLAB環(huán)境下設計FIR數(shù)字高通濾波器9</p><p><b> 5.1設計要求9</b></p><p> 5.2 FIR數(shù)字高通濾波器程序設計10</p><p> 5.3調試結果11</p><p>
8、 5.4 高通FIR數(shù)字濾波器的進一步設計12</p><p> 6 高通FIR數(shù)字濾波器性能測試14</p><p> 6.1高通FIR數(shù)字濾波器性能測試程序14</p><p> 6.2 性能測試結果15</p><p> 7 FDATOOL工具箱設計高通FIR濾波器16</p><p> 7
9、.1 FDATOOL工具箱16</p><p> 7.2 FIR濾波器參數(shù)設置17</p><p><b> 8心得體會19</b></p><p><b> 參考文獻20</b></p><p> 附件:MATLAB程序21</p><p><b&g
10、t; 摘 要</b></p><p> 數(shù)字濾波是語音和圖象處理、模式識別、頻譜分析等應用中的一個基本處理算法,數(shù)字濾波技術是信號消噪的基本方法。根據(jù)噪聲頻率分量的不同,可選用具有不同濾波特性的數(shù)字濾波器。當噪聲的頻率低于信號的頻率時,就應該選用高通濾波器。本設計采用了頻率抽樣法設計的FIR高通數(shù)字濾波器,其目的是為了讓中高頻率的信號通過,而且利用頻率抽樣法的優(yōu)點是可以在頻域直接設計,并且適合最
11、優(yōu)化設計。FIR濾波器為有限長沖激響應濾波器,因其在線性相位特性方面具有獨特的優(yōu)點,因此也越來越受到廣泛的重視。</p><p> MATLAB是一種面向科學和工程計算的語言,它集數(shù)值分析、矩陣運算、信號處理和圖形顯示于一體,具有編程效率高、調試手段豐富、擴展能力強等特點,MATLAB的信號處理工具箱具有強大的函數(shù)功能,它不僅可以用來設計數(shù)字濾波器,還可以使設計達到最優(yōu)化,是數(shù)字濾波器設計的強有力工具。<
12、/p><p> 關鍵詞: FIR數(shù)字濾波器 頻率抽樣 高通 MATLAB</p><p><b> Abstrct</b></p><p> The digital filter is the voice and image processing, pattern recognition, spectrum analysis in the
13、 application of a basic algorithm, digital filtering technology is the basic method of signal denoising. According to different frequency components of noise, can be used with different filtering characteristics of digit
14、al filter. When the noise of high frequency below the signal frequency, should choose a high-pass filter. This design uses the frequency sampling design method of FIR high pass digi</p><p> MATLAB is a scie
15、nce and engineerin-oriented computing language,which set numerical analysis,matrix xomputing,graphics and signal processing in one,the program has a highly efficient means of debugging a rich,strong,and other characteris
16、tics of the expansion,MATLAB signal processing tools Box features a powerful function,which not only can be used to design is a powerful tool.</p><p> Key word: FIR digital filter Frequency sampling Hig
17、h passt. MATLAB</p><p><b> 1 引言</b></p><p> 1.1MATLAB介紹</p><p> MATLAB是矩陣實驗室(Matrix Laboratory)的簡稱,是美國MathWorks公司出品的商業(yè)數(shù)學軟件,用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術計算語言和交互式環(huán)境,主要
18、包括MATLAB和Simulink兩大部分。 </p><p> MATLAB是由美國mathworks公司發(fā)布的主要面對科學計算、可視化以及交互式程序設計的高科技計算環(huán)境。它將數(shù)值分析、矩陣計算、科學數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中,為科學研究、工程設計以及必須進行有效數(shù)值計算的眾多科學領域提供了一種全面的解決方案,并在很大程度上擺脫了傳統(tǒng)非交互式程序設計語
19、言(如C、Fortran)的編輯模式,代表了當今國際科學計算軟件的先進水平。 </p><p> MATLAB和Mathematica、Maple并稱為三大數(shù)學軟件。它在數(shù)學類科技應用軟件中在數(shù)值計算方面首屈一指。MATLAB可以進行矩陣運算、繪制函數(shù)和數(shù)據(jù)、實現(xiàn)算法、創(chuàng)建用戶界面、連 接其他編程語言的程序等,主要應用于工程計算、控制設計、信號處理與通訊、圖像處理、信號檢測、金融建模設計與分析等領域。 M
20、ATLAB的基本數(shù)據(jù)單位是矩陣,它的指令表達式與數(shù)學、工程中常用的形式十分相似,故用MATLAB來解算問題要比用C,F(xiàn)ORTRAN等語言完成相同的事情簡捷得多,并且MATLAB也吸收了像Maple等軟件的優(yōu)點,使MATLAB成為一個強大的數(shù)學軟件。在新的版本中也加入了對C,F(xiàn)ORTRAN,C++ ,JAVA的支持??梢灾苯诱{用,用戶也可以將自己編寫的實用程序導入到MATLAB函數(shù)庫中方便自己以后調用,此外許多的MATLAB愛好者都編寫了
21、一些經典的程序,用戶可以直接進行下載就可以用。 </p><p> 1.2MATLAB信號處理工具箱函數(shù)介紹</p><p> 以下就本課程設計所用到的MATLAB函數(shù)進行簡要介紹。 </p><p> 1、freqs 求模擬濾波器Ha(s)的頻率響應函數(shù)。 </p><p> H=freqs(B,A,w) 計算由向量w(rad/
22、s)指定的頻率點上模擬濾波器Ha(s)的頻率響應Ha(jw),結果存于H向量中。向量B和A分別為模擬濾波器系統(tǒng)()aHs的分子和分母多項式系數(shù)。 </p><p> [H,w]=freqs(B,A,w) 計算出M個頻率點上的頻率響應存于H向量中,M個頻率存放在向量w 中。freqs函數(shù)自動將這M個頻點設置在適當?shù)念l率范圍。默認w和M時freqs自動選取200個頻率點計算。不帶左端輸出向量時,freqs函數(shù)將自
23、動繪出幅頻和相頻曲線。</p><p> 1.3濾波器的介紹 </p><p> 1.3.1濾波器功能及分類</p><p> 濾波器主要功能是對信號進行處理,保留信號中的有用成分,去除信號中的無用成分。其按處理的信號可分為數(shù)字濾波器(Digital Filter,DF)和模擬濾波器(AnalogFilter,AF),按頻域特性分為低通、高通、帶通、帶阻濾波器
24、,按時域特性可分為有限長沖激響應(FIR)濾波器和無限長沖激響應(IIR)濾波器。</p><p> 1.3.2模擬濾波器設計理論</p><p> 模擬濾波器的理論和設計方法已發(fā)展得相當成熟,且有若干典型的模擬低通濾波器的設計原型可供選擇,如巴特沃斯(Butterworth)濾波器、切比雪夫(Chebyshev)濾波器、橢圓(Ellips)濾波器、貝塞爾(Bessel)濾波器等。這些
25、濾波器各有特點,巴特沃斯濾波器具有通帶內最平坦且單調下降的幅頻特性;切比雪夫濾波器的幅頻特性在通帶或阻帶內有波動,可以提高選擇性;貝塞爾濾波器通帶內有較好的線性相位特性;而橢圓濾波器的選擇性相對前三種是最好的。 模擬低通濾波器的設計是最基本的,而高通、帶通、帶阻濾波器則可利用頻率轉換的方法由低通濾波器映射而得到。模擬濾波器的設計是根據(jù)一組設計規(guī)范來設計模擬系統(tǒng)函數(shù)Ha(s),使其逼近某個理想濾波器的特性。其中可以由幅度平方函數(shù)確定系統(tǒng)函
26、數(shù)。 </p><p> 下面介紹兩種常用的低通濾波器特性。一般以低通濾波器為基礎來討論逼近函數(shù),而高通、帶通、帶阻濾波器則可用變換方法有低通濾波器映射而得到。一種是巴特沃斯低通逼近,另一種是切比雪夫低通逼近。本設計中選用第一種方法巴特沃斯低通逼近。 </p><p> 由模擬低通濾波器系統(tǒng)函數(shù)確定模擬高通濾波器系統(tǒng)函數(shù)的方法如下: </p><p> (1)
27、定低通系統(tǒng)函數(shù)HLP(s),其參考角頻率?c(一般為截止頻率)由高 通參考角頻率(一般為截止頻率)?c選定,一般都選?c=1的歸一化原型低通濾波器;</p><p> (2)所得到的HLP(s)中代入變換關系式中,得到高通系統(tǒng)函數(shù)</p><p><b> ?。?)</b></p><p> 故模擬高通濾波器的實現(xiàn)可由模擬低通濾
28、波器的歸一化原型再經頻率變換得到。</p><p> 1.3.3數(shù)字濾波器介紹</p><p> 數(shù)字濾波器是在模擬濾波器的基礎上發(fā)展起來的,從結構上看,數(shù)字濾波器可以分內遞歸型(IIR)與非遞歸型(FIR)兩大類。與模擬濾波器之間存在著一些重要的差別。相比起來,數(shù)字濾波器具有精度高、穩(wěn)定性好、設計靈活、不存在阻抗匹配、便于大規(guī)模集成和可以實現(xiàn)多維濾波等優(yōu)點。在一般情況下,數(shù)字濾波器是
29、一個線性非移變系統(tǒng)。從頻域特性上看,它與模擬濾波器一樣,有低通、高通、帶通和帶阻之分。但在時域的實現(xiàn)方法與方式上,它們是完全不同的兩類系統(tǒng)。</p><p> 數(shù)字濾波器(digital filter)是由數(shù)字乘法器、加法器和延時單元組成的一種裝置。其功能是對輸入離散信號的數(shù)字代碼進行運算處理,以達到改變信號頻譜的目的。由于電子計算機技術和大規(guī)模集成電路的發(fā)展,數(shù)字濾波器已可用計算機軟件實現(xiàn),也可用大規(guī)模集成數(shù)
30、字硬件實時實現(xiàn)。數(shù)字濾波器廣泛用于數(shù)字信號處理中,如電視、VCD、音響等。 按照濾波電路的工作頻帶為其命名:設截止頻率為fp,頻率低于fp的信號可以通過,高于fp的信號被衰減的電路稱為低通濾波器,頻率高于fp的信號可以通過,低于fp的信號被衰減的電路稱為高通濾波器;而帶通嗎,就是頻率介于低頻段截止頻率和高頻段截止頻率的信號可以通過的電路。</p><p> 2 FIR數(shù)字濾波器設計原理</p>&
31、lt;p> FIR濾波器即有限長沖激響應濾波器,突出優(yōu)點是:系統(tǒng)總是穩(wěn)定的,易于實現(xiàn)線性相位、允許設計多通帶(或多阻帶)濾波器,但與IIR濾波器相比,在滿足同樣阻帶衰減的情況下需要的階數(shù)較高,濾波器的階數(shù)越高,占用的運算時間越多,因此在滿足指標要求的情況下應盡量減少濾波器的階數(shù)。</p><p> FIR濾波器的基本結構可以理解為一個分節(jié)的延時線,把每一節(jié)的輸出加權累加,可得到濾波器的輸出,F(xiàn)IR濾波器
32、的沖激響應h(n)是有限長的,數(shù)字上M階FIR濾波器可以表示為:</p><p><b> ?。?)</b></p><p><b> 其系統(tǒng)函數(shù)為:</b></p><p><b> ?。?)</b></p><p> 在自適應處理、數(shù)據(jù)通信等領域中往往要求信號在傳輸過程
33、中不能有明顯的相位失真,F(xiàn)IR濾波器可以做到線性相位滿足此要求。FIR濾波器實質上是一個分節(jié)的延遲線,把每一節(jié)的輸出加權累加,得到濾波器的輸出。對于FIR濾波器的單位脈沖響應h(i)只要滿足以下兩個條件之一,則為線性相位濾波器。線性相位的FIR濾波器具有中心對稱的特性,其對稱中心在N/2處。</p><p> h(i) = h(N-1-i) 偶對稱 (4)</p>
34、<p> h(i) = -h(N-1-i) 奇對稱 (5)</p><p> FIR濾波器的設計問題實質是確定所滿足要求的轉移序列或脈沖響應的常數(shù)的問題,設計方法主要有窗函數(shù)法、頻率取樣法和等波紋最佳逼近法等。</p><p> 3 FIR數(shù)字濾波器設計方法</p><p><b> 3.1窗函數(shù)法
35、</b></p><p> 窗函數(shù)法就是設計FIR數(shù)字濾波器的最簡單也是最常用的方法。它在設計FIR數(shù)字濾波器中有很重要的作用,正確的選擇窗函數(shù)可以提高設計數(shù)字濾波器的性能,或者在滿足設計要求的情況下,減少FIR濾波器的階次。</p><p> 各種窗函數(shù)的性能比較見表1,在設計FIR濾波器的過程中可以根據(jù)要求選擇合適的窗函數(shù):</p><p>
36、表1 各種窗函數(shù)的性能比較</p><p><b> 3.2頻率取樣法</b></p><p> 窗函數(shù)設計FIR數(shù)字濾波器是從時域出發(fā),把理想的濾波器的單位取樣響應Hd(n)用合適的窗函數(shù)截短成為有限長度的H(n),并使H(n)逼近理想的Hd(n),以實現(xiàn)所設計的濾波器的頻率響應Hd(ejw)逼近與理想濾波器的頻率響應Hd(ejw)。</p>&l
37、t;p> 一個有限長的序列,如果滿足頻率采樣定理,可以通過頻譜的有限個采樣點的值被準確的得以恢復。</p><p> 設理想濾波器的頻率響應是Hd(ejw),它是連續(xù)頻率w的周期函數(shù)。對其抽樣,使每一個周期有N個抽樣值,再對它做IDFT,可得到N點的單位抽樣序列Hd(n),即:</p><p> ,n=0,1,2,…………,N-1 (6)</p><p&
38、gt; 將h(n)作為所設計的濾波器的單位沖激響應,這樣就可以求出該濾波器的傳輸函數(shù),即:</p><p><b> ?。?)</b></p><p> 當然H(z)也可以用Hd(k)來表示,即:</p><p><b> ?。?)</b></p><p><b> 經推導,有:&l
39、t;/b></p><p><b> ?。?)</b></p><p> 由此得出的濾波器,其頻率響應在l=mk的抽樣點上嚴格的等于所希望的值Hd(k),而在l≠mk的點上,H(ejw)則由內插函數(shù)的插值決定。所以這種濾波器的設計方法被稱之為頻率抽樣法,該內插函數(shù)是:</p><p><b> ?。?0)</b>&
40、lt;/p><p> 在頻率抽樣法中指定Hd(k)要比窗函數(shù)法中指定Hd(ejw)復雜。</p><p> Hd(k)指定的原則是:</p><p> 在通帶內,可令|Hd(k)|=1,阻帶內|Hd(k)|=0,且在通帶內賦給Hd(k)相位函數(shù);</p><p> 指定的Hd(k)應保證 h(n)是實數(shù);</p><p
41、> 由抽樣序列h(n)求出的Hd(ejw)應具有線性相位。</p><p> 為保證所設計的濾波器具有線性相位,必須對頻率采樣值進行約束。</p><p><b> 有以下四種情況:</b></p><p> 抽樣序列h(n)偶對稱,N為奇數(shù)。</p><p> 抽樣序列h(n)偶對稱,N為偶數(shù)。<
42、/p><p> 抽樣序列h(n)奇對稱,N為奇數(shù)。</p><p> 抽樣序列h(n)奇對稱,N為偶數(shù)。</p><p> 對于理想濾波器的頻率響應Hd(ejw),應根據(jù)第一個采樣點的不同,可分為兩種采樣方法:</p><p> 方法1:第一個采樣點在w=0處</p><p><b> N為奇數(shù)<
43、/b></p><p> 若理想濾波器的頻率響應Hd(ejw)不存在相移時,即設計Ⅰ型FIR濾波器。</p><p> 若理想濾波器的頻率響應Hd(ejw)存在π/2相移時,即設計Ⅲ型FIR濾波器。</p><p><b> N為奇數(shù)</b></p><p> 若理想濾波器的頻率響應Hd(ejw)不存在相移
44、時,即設計Ⅱ型FIR濾波器。</p><p> 若理想濾波器的頻率響應Hd(ejw)存在π/2相移時,即設計Ⅳ型FIR濾波器。</p><p> 方法2:第一個采樣點在w=π/N處</p><p><b> N為奇數(shù)</b></p><p> 若理想濾波器的頻率響應Hd(ejw)不存在相移時,即設計Ⅰ型FIR濾波
45、器。</p><p> 若理想濾波器的頻率響應Hd(ejw)存在π/2相移時,即設計Ⅱ型FIR濾波器。</p><p><b> N為偶數(shù)</b></p><p> 若理想濾波器的頻率響應Hd(ejw)不存在相移時,即設計Ⅱ型FIR濾波器。</p><p> 若理想濾波器的頻率響應Hd(ejw)存在π/2相移時,
46、即設計Ⅳ型FIR濾波器。</p><p> 4 頻率采樣法實際FIR高通濾波器</p><p><b> 4.1設計原理</b></p><p> 4.1.1重構FIR的單位抽樣響應</p><p> 根據(jù)頻率抽樣定理,有N個頻率抽樣點可以唯一確定h(n),即對H(k)進行IDFT變換:</p>&
47、lt;p> ,n=0,1,2,……,N-1 (11)</p><p> 4.1.2重構系統(tǒng)函數(shù)</p><p><b> = </b></p><p><b> =</b></p><p><b> =</b></p><p
48、> = (12)</p><p> 4.1.3 FIR的頻率響應</p><p> 將z=ejw代入H(z)表達式可得:</p><p><b> =</b></p><p> = (13)</p><p>
49、 其中,為大家所知的內插函數(shù),分析可知,當,i=0,1,2,……,N-1時(采樣點)有:</p><p> i=0,1,2,……,N-1 (14)</p><p> 這說明,重構的頻率響應H(ejw)在采樣上嚴格等于H(k),而在采樣點之間,頻率響應則由加權的內插函數(shù)延伸疊加而成。</p><p> 4.1.4 線性相位的約束條件</p>
50、<p> 以h(n)為偶對稱,N為奇數(shù)的情況進行分析。</p><p> FIR濾波器的頻率響應具有線性相位的一般表達式</p><p> 當h(n)為偶對稱,N為奇數(shù)時,則</p><p><b> (15)</b></p><p> 而且幅度函數(shù)H(w)應為偶對稱,即</p>&
51、lt;p><b> ?。?6)</b></p><p> 采樣點H(k)具有線性相位的約束條件:</p><p><b> ?。?7)</b></p><p> 其中,表示采樣值的模(純標量),θk表示其相角。因此,在采樣點上具有線性相位的條件是:</p><p><b> (
52、18)</b></p><p> 而且,Hk必須滿足偶對稱,即</p><p><b> (19)</b></p><p><b> 4.2設計步驟</b></p><p> 1.根據(jù)指標要求,畫出頻率采樣序列的圖形;</p><p> 2.根據(jù)|Hk|
53、的對稱特點,可以使問題得以簡化;</p><p> 3.根據(jù)線性相位的約束條件,求出;</p><p> 4.將代入FIR濾波器的頻率響應表達式;</p><p> 5.有Hk的表達式畫出實際頻率響應。</p><p> 5 MATLAB環(huán)境下設計FIR數(shù)字高通濾波器</p><p><b> 5.
54、1設計要求</b></p><p> 1. 在數(shù)字信號處理平臺上(PC機﹑MATLAB仿真軟件系統(tǒng))進行軟件仿真設計,并進行調試和數(shù)據(jù)分析。</p><p> 2. 利用MATLAB仿真軟件系統(tǒng)結合頻率取樣法設計一個數(shù)字高通FIR濾波器。</p><p> 3. 性能指標如下:阻帶截止頻率0.5π,通帶截止頻率0.6π,阻帶衰減不小于30dB,通帶
55、衰減不大于1dB。</p><p> 5.2 FIR數(shù)字高通濾波器程序設計</p><p> 基于頻率采樣法的FIR數(shù)組高通濾波器程序如下:</p><p><b> N=33;%采樣點</b></p><p> wc=3*pi/5;%通帶截止頻率</p><p> N1=fix(wc/
56、(2*pi/N));%阻帶上采樣點</p><p> N2=N-2*N1;%通帶上采樣點</p><p> A=[zeros(1,N1),ones(1,N2),zeros(1,N1)];</p><p> theta=-pi*[0:N-1]*(N-1)/N;%線性相位約束條件</p><p> H=A.*exp(j*theta);%頻
57、率采樣</p><p> h=real(ifft(H));%傅里葉逆變換</p><p><b> v=1:N;</b></p><p> figure(1);</p><p> subplot(211)</p><p> plot(v,A,'*');</p>
58、<p> title('頻率采樣');</p><p> ylabel('H(k)');</p><p> axis([0,fix(N*1.1),-0.1,1.1]);</p><p> subplot(212)</p><p> stem(v ,h);title(‘脈沖響應’);yla
59、bel('h(n)');</p><p> axis([0,fix(N*1.1),min(h)*1.1,max(h)*1.1]); </p><p><b> M=5000;</b></p><p><b> nx=[1:N];</b></p><p> w=linspace
60、(0,pi,M);%0~π范圍等間隔</p><p> X=h*exp(-j*nx'*w);%求取頻率響應</p><p> figure(2);</p><p> plot(w./pi,abs(X),'k');xlabel('\omega/\pi');ylabel('|Hd(w)|');</p&g
61、t;<p> title('幅度響應');axis([0,1,-0.1,1.3]);</p><p> figure(3);</p><p> plot(w./pi,20*log10(abs(X)),'k');title(‘幅度特性’);</p><p> xlabel('\omega/\pi')
62、;ylabel('dB');axis([0,1,-80,10]);</p><p><b> 5.3調試結果</b></p><p> 圖1 理想高通濾波器頻率采樣情況及其脈沖響應</p><p> 分析:從頻率抽樣樣本及其脈沖響應來看,抽樣點的選擇滿足線性相位的條件,即在0~2π上等間隔采樣33個點,滿足對于高通和帶阻濾
63、波器,N只能取奇數(shù)的要求。同時,采樣點的選取9個0點,15個1點,9個0點,這樣在2π的區(qū)間上在π的附近就是通帶,滿足高頻濾波器頻帶處于π的奇數(shù)倍附近。</p><p> 圖2 高通FIR濾波器的幅頻響應曲線</p><p> 分析:從高通濾波器的幅頻響應曲線可以看出,大于0.6π區(qū)間上處于通帶,滿足高通濾波器阻帶截止頻率0.6π的要求,而阻帶截止頻率截止于0.5π同樣滿足設計的要求。
64、</p><p> 圖3 高通濾波器幅度特性</p><p> 分析:從高通濾波器幅度響應曲線可以看出,在通帶內高通濾波器幅度衰減滿足小于1dB的設計要求。而在阻帶內同樣滿足大于30dB的衰減要求。說明設計的高通FIR濾波器滿足設計的性能要求,達到了設計的標準。</p><p> 5.4 高通FIR數(shù)字濾波器的進一步設計</p><p>
65、; 5.4.1高通FIR濾波器改進措施分析</p><p> 以上設計的FIR數(shù)字濾波器滿足設計的要求,但是我們發(fā)現(xiàn)在阻帶內信號存在很大的震蕩,但是實際上我們總是希望信號在通帶內無衰減,而在阻帶內沒有信號通過。這也是理想濾波器的設計要求,但是實際上并不能達到理想的標準。但我們可以通過一定的手段減小阻帶內的信號震蕩,使之趨近于理想濾波器的特性要求。</p><p> 對于此次設計的高通
66、濾波器,我們可以通過增加頻率采樣點,即在頻率響應間斷點附近區(qū)間內插一個或幾個過渡采樣點,使不連續(xù)點變成緩慢過渡,這樣,雖然加大了過渡帶,但明顯增大了阻帶衰減。</p><p> 5.4.2 高通FIR濾波器改進設計程序</p><p><b> N=33;</b></p><p> wc=3*pi/5;T1=0.2;T2=0.7;%增加過
67、渡點值</p><p> N1=fix(wc/(2*pi/N));</p><p> N2=N-2*N1;</p><p> A=[zeros(1,N1),T1,T2,ones(1,N2-4),T2,T1,zeros(1,N1)];%增加過渡點</p><p> theta=-pi*[1:N]*(N-1)/N;</p>
68、<p><b> ………………</b></p><p> 在這個采樣頻率點的選取上,新增了四個過渡點在頻率響應間斷點上,即0.2,0.7,這樣的方法是提高阻帶衰減的最有效的方法。由于程序設計大部分不變,所以這部分只給出了部分程序。</p><p> 5.4.3 改進型高通FIR濾波器調試結果</p><p> 圖4 理想高通濾
69、波器增加頻率過度點</p><p> 圖5 改進型高通FIR濾波器幅度特性</p><p> 分析:在頻率響應間斷點附近區(qū)間內插了幾個過渡采樣點,是不連續(xù)的點變成緩慢過渡,我們可以看到,在改進后的高通FIR濾波器特性曲線中,可以看到明顯使得阻帶內的衰減幅度大大增加,達到了-50dB,有效減小了阻帶震蕩,使設計的高通FIR濾波器性能更加穩(wěn)定,更加趨近于理想高通FIR濾波器。</p&
70、gt;<p> 6 高通FIR數(shù)字濾波器性能測試</p><p> 6.1高通FIR數(shù)字濾波器性能測試程序</p><p> f1=200;f2=500;f3=1000;%輸入信號頻率</p><p> fs=3000;%采樣頻率</p><p> figure(1) </p><p> su
71、bplot(211) </p><p> t=0:1/fs:0.25;%定義時間范圍和步長</p><p> s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%濾波前信號</p><p> plot(t,s);%濾波前信號圖像 </p><p> xlabel('時間/秒'
72、); ylabel('幅度');title('信號濾波前時域圖'); </p><p> subplot(212) </p><p> Fs=fft(s,512);%將信號變換到頻域 </p><p> AFs=abs(Fs);%信號頻域圖的幅值</p><p> f=(0:255)*fs/512;%頻
73、率采樣 </p><p> plot(f,AFs(1:256));%濾波前信號頻域圖</p><p> xlabel('頻率/赫茲'); ylabel('頻譜'); title('信號濾波前頻譜圖'); </p><p> figure(2) </p><p> sf=filter(h,
74、1,s);</p><p> subplot(211) </p><p> plot(t,sf)%濾波后信號的圖像</p><p> xlabel('時間/秒'); ylabel('幅度'); title('信號濾波后時域圖'); </p><p> axis([0.2 0.25 -2
75、2]);%限定圖像坐標范圍</p><p> subplot(212) </p><p> Fsf=fft(sf,512);%濾波后信號頻域</p><p> AFsf=abs(Fsf);%信號頻域的幅值</p><p> f=(0:255)*fs/512;%頻率采樣</p><p> plot(f,AFsf
76、(1:256))%濾波后信號頻域圖</p><p> xlabel('時間/秒'); ylabel('幅度'); title('信號濾波后時域圖'); </p><p> 6.2 性能測試結果</p><p> 圖6輸入信號波形與頻譜圖</p><p> 圖7輸出信號波形與頻譜圖<
77、/p><p> 分析:從輸入輸出信號的波形和頻譜圖可以看出,輸入信號含有200Hz,500Hz,1000Hz的混合信號,經過高通FIR濾波器之后,只有高頻1000Hz的信號通過,低頻信號全部被過濾掉,檢驗了所設計的高通FIR濾波器的性能,說明所設計的高通FIR濾波器滿足設計的要求,達到了很好的濾波效果。</p><p> 7 FDATOOL工具箱設計高通FIR濾波器</p>
78、<p> 7.1 FDATOOL工具箱</p><p> fdatool(filter design & analysis tool)是matlab信號處理工具箱里專用的濾波器設計分析工具。fdatool可以設計幾乎所有的基本的常規(guī)濾波器,包括fir和iir的各種設計方法。它操作簡單,方便靈活。進入fdatool工具箱在matlab的命令行上輸入fdatool,回車,如圖所示:</p&
79、gt;<p> 圖8fdatool工作欄</p><p> 7.2 FIR濾波器參數(shù)設置</p><p> fdatool界面總共分兩大部分,一部分是design filter,在界面的下半部,用來設置濾波器的設計參數(shù),另一部分則是特性區(qū),在界面的上半部分,用來顯示濾波器的各種特性。design filter部分主要分為: filter type(濾波器類型)選項,包括
80、lowpass(低通)、highpass(高通)、bandpass(帶通)、bandstop(帶阻)和特殊的fir濾波器。 design method(設計方法)選項,包括iir濾波器的butterworth(巴特沃思)法、chebyshev type i(切比雪夫i型)法、 chebyshev type ii(切比雪夫ii型) 法、elliptic(橢圓濾波器)法和fir濾波器的equiripple法、least-squares(最小
81、乘方)法、window(窗函數(shù))法。 filter order(濾波器階數(shù))選項,定義濾波器的階數(shù),包括specify order(指定階數(shù))和minimum order(最小階數(shù))。在specify order中填入所要設計的濾波器的階數(shù)(n階濾波器,specify order=n-1),如果</p><p> 圖9 FDATOOL工具箱設計高通FIR濾波器 </p><p>&l
82、t;b> 8心得體會</b></p><p> 這學期我們學習了《數(shù)字信號處理》,這是我們的專業(yè)課,是理論方面的指示。正所謂“紙上談兵終覺淺,覺知此事要躬行”。學習任何知識,僅從理論上去求知,而不去實踐、探索是不夠的,所以在暨《數(shù)字信號處理》之后緊接著來一次MATLAB課程設計是很及時、很必要的。這樣不僅能加深我們對數(shù)字信號處理的認知,而且還及時、真正的做到了學以致用。</p>
83、<p> 在此過程中,我們通過查找大量資料,請教老師,以及不懈的努力,不僅培養(yǎng)了獨立思考、動手操作的能力,在各種其它能力上也都有了提高。更重要的是,在課設中,我們學會了很多學習的方法,而這也是日后最實用的。</p><p> 不管怎樣,這些都是一種鍛煉,一種知識的積累,能力的提高。完全可以把這個當作基礎東西,只有掌握了這些最基礎的,才可以更進一步,取得更好的成績。很少有人會一步登天,永不言棄才是最
84、重要的。而且,這對于我們的將來也有很大的幫助。</p><p><b> 參考文獻</b></p><p> 1.《數(shù)字信號處理及其MATLAB實現(xiàn)》. Vinay K.Ingle主編.電子工業(yè)出版社.1998年出版</p><p> 2.《Digital Signal Processing Laboratory Using Matlab
85、》. Sanjit K. Miltra編著. McGraw-Hill出版社.2000年出版</p><p> 3.《數(shù)字信號處理原理與實現(xiàn)》第二版.劉泉主編.電子工業(yè)出版社</p><p> 4.《數(shù)字信號處理教程—MATLAB釋義與實現(xiàn)》.陳懷琛主編.電子工業(yè)出版社.2004年版</p><p> 5.《信號與系統(tǒng)》.劉泉主編.高等教育出版社.2006年版&
86、lt;/p><p> 6.《MATLAB程序設計》.阮沈青主編.電子工業(yè)出版社.2004年版附件:MATLAB程序</p><p><b> N=33;</b></p><p> wc=3*pi/5;T1=0.2;T2=0.7;</p><p> N1=fix(wc/(2*pi/N));</p><
87、;p> N2=N-2*N1;</p><p> A=[zeros(1,N1),T1,T2,ones(1,N2-4),T2,T1,zeros(1,N1)];</p><p> theta=-pi*[1:N]*(N-1)/N;</p><p> H=A.*exp(j*theta);</p><p> h=real(ifft(H))
88、;</p><p><b> v=1:N;</b></p><p> figure(1);</p><p> plot(v,A,'*');</p><p> title('頻率采樣');</p><p> ylabel('H(k)');&l
89、t;/p><p> axis([0,fix(N*1.1),-0.1,1.1]);</p><p> figure(2);</p><p> stem(v ,h);title('脈沖響應'); ylabel('h(n)');</p><p> axis([0,fix(N*1.1),min(h)*1.1,max
90、(h)*1.1]);</p><p><b> M=5000;</b></p><p><b> nx=[1:N];</b></p><p> w=linspace(0,pi,M);</p><p> X=h*exp(-j*nx'*w);</p><p>
91、figure(3);</p><p> plot(w./pi,abs(X),'k');xlabel('\omega/\pi');ylabel('|Hd(w)|');</p><p> title('幅度響應');axis([0,1,-0.1,1.3]);</p><p> figure(4);&l
92、t;/p><p> plot(w./pi,20*log10(abs(X)),'k');title('幅度響應');</p><p> xlabel('\omega/\pi');ylabel('dB');axis([0,1,-80,10]);</p><p> n=0:N-1;figure(5);<
93、;/p><p> freqz(abs(h));</p><p> f1=200;f2=500;f3=1000;</p><p><b> fs=3000; </b></p><p> figure(6) </p><p> subplot(211) </p><p>
94、 t=0:1/fs:0.25; </p><p> s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t); </p><p> plot(t,s); </p><p> xlabel('時間/秒'); ylabel('幅度');title('信號濾波前時域圖'); <
95、;/p><p> subplot(212) </p><p> Fs=fft(s,512); </p><p> AFs=abs(Fs); </p><p> f=(0:255)*fs/512; </p><p> plot(f,AFs(1:256)); </p><p> xlabel
96、('頻率/赫茲'); ylabel('幅度'); title('信號濾波前頻域圖'); </p><p> figure(7) </p><p> sf=filter(h,1,s);</p><p> subplot(211) </p><p> plot(t,sf)</p>
97、<p> xlabel('時間/秒'); ylabel('幅度'); title('信號濾波后時域圖'); </p><p> axis([0.2 0.25 -2 2]);</p><p> subplot(212) </p><p> Fsf=fft(sf,512); </p>&
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- dsp課程設計fir數(shù)字濾波器
- 數(shù)字信號課程設計---fir數(shù)字濾波器
- fir數(shù)字濾波器設計及matlab實現(xiàn)
- fir數(shù)字濾波器的dsp課程設計
- 數(shù)字濾波器課程設計--數(shù)字濾波器的設計及實現(xiàn)
- 基于matlab的fir數(shù)字濾波器設計與優(yōu)化
- 基于labview的fir數(shù)字濾波器設計課程設計
- fir數(shù)字濾波器的設計
- 數(shù)字信號處理課程設計報告--fir數(shù)字濾波器設計及matlab實現(xiàn)
- 數(shù)字信號處理課程設計--fir數(shù)字濾波器的設計
- FIR數(shù)字濾波器的設計.pdf
- fir數(shù)字濾波器的設計方法
- 畢業(yè)設計----基于matlab的fir數(shù)字濾波器設計
- FIR數(shù)字濾波器的MATLAB設計與DSP實現(xiàn).pdf
- 基于matlab的數(shù)字濾波器的設計課程設計
- 基于matlab的iir數(shù)字濾波器的課程設計
- 三fir數(shù)字濾波器的設計
- iir數(shù)字濾波器課程設計
- 基于fir的數(shù)字濾波器設計研究
- 基于fir的數(shù)字濾波器設計研究
評論
0/150
提交評論