【文章內容簡介】 Kc? ?0? ??j 圖 21 理想低通濾波器頻率特性 圖 22 理性濾波器的單位 沖激 響應 假定信號頻率成分為 ： c??? ，噪音頻率成分為 c??? 。則完成濾波的理想低通濾波器特性是 ??1 ： | |()() () 0 | |j ccKeYjHj Xj ???? ? ? ? ??? ? ? ?? ? ? ?? （ 29） 即： | || ( ) | 0 | | ccKHj ? ? ???? ? ? ? ?? （ 210） arg( ( ))Hj ?? ? ?? （ 211） 系統(tǒng)的單位 沖激 響應為 ??1 ： ? ? ? ?? ?? ??? ?? ? ? ? ? ?????? ???? ??? ttKdeKeth cjj s i n2 1 （ 212） 理性低通濾波器的頻率特性如圖 21 所示，單位 沖激 響應的波形如圖 22 所示。 h(t) t 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 7 理想濾波器具有非因果、無限長的單位脈沖響應 和不連續(xù)的頻率特性，要用穩(wěn)定的線性時不變（ LTI）系統(tǒng)來實現(xiàn)這樣的特性是不可能的。工程上是用 沖激響應為有限長的、因果的、穩(wěn)定的線性時不變系統(tǒng)或具有連續(xù)頻率特性的線性時不變系統(tǒng)來逼近理想特性。在滿足一定的誤差要求的情況下來實現(xiàn)理想濾波特性。因此實際的濾波器的頻率特性如圖 23 所示。 0通 帶 阻 帶 過 渡 帶 )( ?jeH11 ? 11 ??2?c?s?? 圖 23 實際濾波器頻率特性 其中： c?—— 截止頻率 s?—— 阻帶起始頻率 cs ???—— 過渡帶寬 在通帶內幅度響應以 1?? 的誤差接近于 1，即 ??1 ： 111 ( ) 1j cHe ?? ? ? ?? ? ? ? ? （ 213） s?為阻帶起始頻率 ， 在阻帶內幅度響應以小于 2? 的誤差接近于零 ， 即 ??1 ： 2()j sHe ? ? ? ? ?? ? ? （ 214） 為了使逼近理想低通濾波器的方法成為可能 ， 還必須提供一帶寬為cs ???的不為零的過渡帶。在這個頻帶內，幅度響應從通帶平滑 地 下落到阻帶。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 8 MATLAB 概述 MATLAB是美國 MathWorks公司開發(fā)的一種功能極其強大的高技術計算語言和內容極其豐富的軟件庫 ， 集數(shù)值計算 、 矩陣運算和信號處理與顯示于一身 。該軟件最初是由美國教授 Cleve Moler創(chuàng)立的 。 1980年前后 ， 他在教線性代數(shù)課程時，發(fā)現(xiàn)用其他高級語言編程時極不方便 ， 便構思開發(fā)了 MATLAB， 即矩陣實驗室（ Matrix Laboratory） 。 該軟件利用了當時代表數(shù)值線性代數(shù)領域最高水 平的 EISPACK和 LINPACK兩大軟件包 ， 并且利用 Fortran語言編寫了最初的一套交互式軟件系統(tǒng)， MATLAB的最初版本便由此 誕生 了 。 最初的 MATLAB由于語言單一 ， 只能進行矩陣的運算 ， 繪圖也只能用原始的描點法 ， 內部函數(shù)只有幾十個 ， 因此功能十分簡單 。 1984年該公司推出了第一個 MATLAB的商業(yè)版 ， 并用 C語言作出了全部改寫?，F(xiàn)在的 MATLAB程序是MathWorks公司用 C語言開發(fā)的 ， 第一版由 Steve Bangert主持開發(fā)編譯解釋程序，Steve Kleiman完成圖形功能的設計， John Little和 Cleve Moler主持開發(fā)了各類數(shù)學分分析的子模塊，撰寫用戶指南和大部分的 M文件。接著又添加了豐富的圖形圖像處理、多媒體功能、符號運算和與其它流行軟件的接口功能，使 MATLAB的功能越來越強大 ??6 。 MATLAB系統(tǒng)主要由以下五個部分組成： （ 1） MATALB語言體系 。 MATLAB是高層次的矩陣／數(shù)組語言，具有條件控制、函數(shù)調用、數(shù)據結構、輸入輸出、面向對象等程序語言特性。利用它既可以進行小規(guī)模端程，完成算法設計和算法實驗的基本任務，也可以進行大規(guī)模編程，開發(fā)復雜的應用程序。 （ 2） MATLAB工作環(huán)境 。 這是對 MATLAB提供給用戶使用的管理功能的總稱。包括管理工作空間中的變量據輸入輸出的方式和方法，以及開發(fā) 、 調試 、 管理 M文件的各種工具。 （ 3） 圖形句相系統(tǒng) 。 這是 MATLAB圖形系統(tǒng)的基礎，包括完成 2D和 3D數(shù)據圖示 、 圖像處理 、 動畫生成 、 圖形顯示等功能的高層 MATLAB命令 ， 也包括用戶對圖形圖像等對象進行特性控制的低層 MATLAB命令，以及開發(fā) GUI應用程序的各種工具。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 9 （ 4） MATLAB數(shù)學函數(shù)庫 。 這是對 MATLAB使用的各種數(shù)學算法的總稱 。包括各種初等函數(shù)的算法 ， 也包括矩陣運算 、 矩陣分析等高層次數(shù)學算法。 （ 5） MATLAB應用程序接口 （ API）。 這是 MATLAB為用戶提供的一個函數(shù)庫 ， 使得用戶能夠在 MATLAB環(huán)境中使用 C程序或 FORTRAN程序 ， 包括從MATLAB中調用于程序 （動態(tài)鏈接）， 讀寫 MAT文件的功能 ??4 。 除此之外， MATLAB系統(tǒng)還具有如下特點： （ 1） 具有易學易用的語言體系； （ 2） 具有交互式的工作環(huán)境 ； （ 3） 具有多層面的圖像處理系統(tǒng) ； （ 4） 具有豐富高效的 MATLAB工具箱 ； （ 5） 具有便利的程序接口（ API） ； （ 6） 應用領域廣泛 。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 10 第 3章 FIR濾波器的設計方法研究 FIR 數(shù)字濾波器的數(shù)學模型和基本結構 1. 直接型：如圖 31 可以看出直接型結構共需要 N 個乘法器，若系數(shù)不對稱則不能設計線性相 位 ??8 。 ? ?nX1?Z 1?Z1?Z ? ?nY? ?0h ? ?1h ? ?2h ? ?2Nh ? ?1Nh圖 31 FIR濾波器的 直接 型結構 2. 級聯(lián)型：這種結構的每一節(jié)控制一對共軛極點，因此調整傳輸零點方便，但是這種結構所需的系數(shù)和所需的乘法運算比直接型多，所以這種結構使用的比較少 ??8 。 ? ?nX 1?Z 1?Z ? ?nY1?Z 1?Z 1?Z 1?Z01?11?21?02?12?22?2N0?2N1?2N2? 圖 32 FIR濾波器的級聯(lián)型結 構 3. 頻率抽樣型：把一個有限長序列（長度為 N 點）的 z 變換 ??zH 在單位圓上作 N 等分抽樣，就得到 ??kH ，其主值序列就等于 h（ n）的離散傅里葉變換 ??kH 。 諧振器的極點正好與梳狀濾波器的零點相抵消，保證了網絡的穩(wěn)定性。 N 個并聯(lián)諧振器與梳狀濾波器級聯(lián)后，得到圖 33 的頻率抽樣結構 ??8 。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 11 ? ?nX 1 ?ZNZ ? ? ?nY? ?0H1 ?Z1 ?Z? ?1H? ?1NH 圖 33 FIR濾波器的頻率抽樣型結構 窗函數(shù)設計法 窗函數(shù)法的設計流程如圖 34所示 ??8 ： 開 始H ( ej w)是 否 滿 足 要 求得 到 所 需 F I R 濾 波 器 的 H ( z ) 或 者 h ( n )結 束重 新 選 擇窗 口 函 數(shù)改 變 窗 口長 度N確 定 逼 近 理 想 濾 波 器 頻 率 響 應 函 數(shù)H d ( ej w)求 理 想 濾 波 器 的 單 位 抽 樣 函 數(shù)h d ( n )選 擇 窗 函 數(shù)w ( n )和 窗 口 長 度N加 窗 ：h ( n ) = hd( n ) w ( n )， 并 求 出H ( ej w)是否 圖 34 窗函數(shù)法設計流程圖 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 12 常用的窗函數(shù)有：矩形窗、漢寧窗、海明窗、布萊克曼窗、凱塞窗、三角窗等。在窗函數(shù)中，凱塞窗是比較靈活的一種窗函數(shù)， 調整凱塞窗中的參數(shù) ? 的大小，我們可以得到不同性能的濾波器。窗函數(shù)設計法是一種常用的設計數(shù)字濾波器的方法，不過在設計精度和性能方面并不是十分理想。 幾種窗函數(shù)的性能比較如表 31 所示： 表 31 各種窗函數(shù)的參數(shù)表 窗函數(shù) 過渡帶寬度（ P/N） 阻帶最小衰減（ dB） 矩形窗 4 21 三角窗 8 25 漢寧窗 8 44 海明窗 8 53 凱塞窗（ ? =） 60 布萊克曼窗 12 74 1. 用窗函數(shù)法設計低通濾波器，滿足如下指標： ?? ?p rad， ?? ?s rad， p? =3dB， s? =50dB 程序如下： %阻帶衰減 50dB，選用海明窗 clear。close all。 wd=*pi。N=133。M=(N1)/2。 nn=M:M。 n=nn+eps。 hd=sin(wd*n)./(pi*n)。 % 理想沖激響應 w=hamming(N)39。 % 海明窗 h=hd.*w。 % 實際沖激響應 H=20*log10(abs(fft(h,1024)))。 % 實際濾波器的幅頻特性 HH=[H(513:1024) H(1:512)]。 subplot(221),stem(nn,hd,39。k39。)。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 13 xlabel(39。n39。)。title(39。理想沖激響應 39。)。axis([70 70 ])。 subplot(222),stem(nn,w,39。k39。)。axis([70 70 ])。 title(39。海明窗 39。)。xlabel(39。n39。)。 subplot(223),stem(nn,h,39。k39。)。 axis([70 70 ])。xlabel(39。n39。)。title(39。實際沖激響應 39。)。 w=(512:511)/511。 subplot(224),plot(w,HH,39。k39。)。 axis([ 140 20])。xlabel(39。\omega/\pi39。)。title(39。濾波器幅頻特性 39。)。 set(gcf,39。color39。,39。w39。)。 運行后結果： 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 14 圖 35 用海明窗設計 FIR低通濾波器的運行結果 該低通濾波器的設計參數(shù)為通帶頻率 ?? ?p ，阻帶頻率 ?? ?s ，過渡帶寬度為 ? ，由運行結果（圖 35）可以看出，設計出的濾波器實際特性過渡帶為 ? ，所以該設計滿足要求。 2. 用窗函數(shù)法設計高通濾波器，滿足如下指標： ?? ?p rad， ?? ?s rad， p? =， s? =50dB 程序如下： %阻帶衰減 50dB，選擇海明窗 clear。close all。 wd=*pi。N=33。M=(N1)/2。 nn=M:M。 西北工業(yè)大學明德學院 本科畢業(yè)設計論文 15 n=nn+eps。 hd=2*((1).^n).*sin(wd*n)./(pi*n)。 % 理想沖激響應 w=hamming(N)39。 % 海明窗 h=hd.*w。 % 實際沖激響應 H=20*log10(abs(fft(h,1024)))。 % 實際濾波器的幅頻特性 HH=[H(513:1024) H(1:512)]。 subplot(221),stem(nn,hd,39。k39。)。