【正文】 矩形窗漢寧窗布拉克曼窗頻譜圖比較濾波效果比較巴特沃斯圖 21 流程圖 1彩色圖片IIR 巴特沃斯低通濾波器FIR 布拉克曼低通濾波器圖片濾波效果比較圖 22 流程圖 2基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)6 數(shù)字濾波器的簡介作為電子濾波器中的一種——數(shù)字濾波器，它和完全工作在模擬域的模擬濾波器不一樣，數(shù)字濾波器工作在數(shù)字域,而且數(shù)字濾波器處理的信號是經(jīng)過 A/D 轉(zhuǎn)換后的數(shù)字信號。正是因為系統(tǒng)中反饋部分的存在，所以 IIR 濾波器的脈沖響應(yīng)就為無限長。因此，對離散時間系統(tǒng)的結(jié)構(gòu)有基本的認(rèn)識是有必要的。作為一個單位脈沖響應(yīng) h(z)是一個有限序列的 FIR 數(shù)字濾波器，它主要采用為濾波結(jié)果非遞歸結(jié)構(gòu)?，F(xiàn)在我們設(shè) FIR 數(shù)字濾波器是一個具有 N 點序列的單位脈沖響應(yīng) h(n)，其中0≤n≤N—1，則濾波器系統(tǒng)函數(shù)為：從函數(shù)是可以得出，在 Z 域中，有 N1 個零點在 Z 平面的任意位置，而在 Z=0 處，具有 N1 個極點。????10)()(NnnzhzH基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)8 數(shù)字濾波器的設(shè)計方法概述有很多種設(shè)計方法來設(shè)計數(shù)字濾波器，比如熟知的雙線性變換、窗函數(shù)、Chebyshev 逼近和差值逼近法等等。幅度指標(biāo)給出方式一般是兩種，絕對指標(biāo)（提供對幅度響應(yīng)函數(shù)的要求）和相對指標(biāo)（以分貝的形式給出要求） 。有兩種系統(tǒng)函數(shù)逼近所用系統(tǒng)函數(shù)，無限長單位沖激響應(yīng)（IIR）系統(tǒng)函數(shù)與有限長單位沖激響應(yīng)（FIR ） 。 數(shù)字濾波器的性能指標(biāo)模擬濾波器按頻率特性劃分有低通、高通、帶通和全阻等類型，數(shù)字濾波器和模擬濾波器類似，按照頻率響應(yīng)也分為四類。圖 24 非理想低通濾波器頻率響應(yīng)在通帶內(nèi)： 1 A P≤| H(ejω)| ≤1 |ω|≤ωc在阻帶中： |H(ejω)| ≤ Ast ωst ≤|ω|≤ωc基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)10ωc 為通帶截止頻率, ω st 為阻帶截止頻率， Ap 為通帶誤差, A st 為阻帶誤差。程序設(shè)計、圖像繪制、數(shù)值計算、文件管理以及輸入輸出等等各項操作用戶可以在此軟件下進行。因為模擬濾波器設(shè)計已經(jīng)具有一系列很容易掌握的方法，而且模擬濾波器有很多簡單現(xiàn)成的設(shè)計模擬濾波器的公式，并且設(shè)計參數(shù)有比較充分的圖標(biāo)來查詢。?cMatlab 信號處理工具箱為低通模擬巴特沃斯濾波器產(chǎn)生提供了 buttap 函數(shù)，它的調(diào)用格式：[z,p,k]=buttap(N),z 表示零點，p 表示極點,k 是表示增益，N 表示濾波器的階次，下面運用 Matlab 語言設(shè)計 Butterworth 低通模擬濾波器：基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)12圖 31 巴特沃斯程序圖 32 階數(shù) N 分別是 14 的巴特沃斯幅頻圖基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)13由上圖可以清楚地分析出，Butterworth 具有平滑的頻率響應(yīng)，而且在截止頻率外，頻率響應(yīng)單調(diào)下降，階數(shù) N 越大，過渡帶的陡峭程度越大，由此可以推出高階的Butterworth 濾波器頻率響應(yīng)近似于理想低通濾波器。但是在設(shè)計中只考慮到了幅度特性，沒有考慮到相位特性，所以設(shè)計的濾波器的相位一般是非線性的。其中對于窗函數(shù)法，本文主要討論矩形窗（Rectangle window） 、漢寧窗（Hanning window） 、布拉克曼窗（blackman window） 。由上圖的矩形窗，漢寧窗以及布拉克曼窗的比較可知矩形窗偶對稱，主瓣的寬度最窄，旁瓣的幅度最大；漢寧窗偶對稱，主瓣寬度中等，旁瓣幅度中等；布拉克曼窗偶對稱，主板寬度最大，旁瓣幅度最小?；?Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)16所以，在現(xiàn)實生活中對于 N 的設(shè)置不能太大，太大了不僅會增加我們的制作成本還會影響濾波后信號的延時。性能：IIR 濾波器的傳遞函數(shù)有零點和極點兩個因素，在單位圓內(nèi)是對極點的唯一限制，所以獲得較高的選擇性可以用比較低的階數(shù) N。相位：如果按照相位要求來看，IIR 需要增加對相位的糾正，因此結(jié)構(gòu)會復(fù)雜很多，而 FIR 濾波器本身就可以保證嚴(yán)格的線性相位。圖 37 原混合信號以及其幅頻圖（1） IIR 和 FIR 帶通濾波器濾波效果比較基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)18基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)19圖 38IIR 和 FIR 帶通濾波器的實現(xiàn)部分程序段結(jié)果：圖 39 混合正弦信號在巴特沃斯帶通濾波器和矩形窗帶通濾波器濾波后結(jié)果由圖 39 可知，IIR 的巴特沃斯帶通濾波器，在通帶臨界頻率 Wp 是 110 和 190，阻帶臨界頻率 Ws 是 80 和 220，通帶內(nèi)最大衰減 Rp 是 1，阻帶內(nèi)最小衰減 Rs 是 50 的情況下，能把 150HZ 的頻率保留下來，把 50HZ 和 300HZ 的頻率濾除，而 FIR 的矩形基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)20窗帶通濾波器在同樣的條件下也能達到預(yù)期的效果?；?Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)21圖 311 巴特沃斯、矩形窗、漢寧窗和布拉克曼窗在同樣條件下的帶通頻譜圖 312 巴特沃斯、矩形窗、漢寧窗和布拉克曼窗在同樣條件下的歸一化帶通頻譜有上圖 311 和圖 312 可知，巴特沃斯、矩形窗、漢寧窗和布拉克曼窗的幅頻圖像都關(guān)于 150HZ 的頻率幾乎偶對稱，并且在 150HZ 處為幅度的最大值，于是很顯然這四種方式都能把 150HZ 的頻率保留下來。從波形圖可看出，經(jīng)過巴特沃斯低通濾波后的波形與標(biāo)準(zhǔn)波形幾乎重合，只是相位發(fā)生了細微的平移。從信號波形圖來看，漢寧窗的低通濾波和布拉克曼的低通濾波效果都不是很好，漢寧床只是幅度和相位有適當(dāng)變化，而布拉克曼低通濾波后不僅幅度和相位變化，其波形的也發(fā)生了包絡(luò)變化。在圖像傳輸過程中，因為它的龐大數(shù)據(jù)量所以還需要很大的帶寬，這不僅是圖像傳輸過程中的特點也是難點，在計算機系統(tǒng)和網(wǎng)絡(luò)上存儲、處理以及傳輸對數(shù)量巨大的數(shù)據(jù)都是很難的。的分號去掉就可以 mand window 中得到 N和 Wn 的值基于 Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)27圖 41 在 mand window 中得出 N 和 Wn 的值在此相同參數(shù)分別利用 FIR 巴特沃斯低通濾波器和 FIR 的布拉克曼窗進行低通濾波。這種原因要從圖 42 巴特沃斯和布拉克曼窗的頻譜圖來解釋，在同樣的條件要求下，巴特沃斯的頻譜圖過渡帶窄，很快就可以讓幅值降低為零，這樣可以最大可能的濾除大量的高頻信息，而布拉克曼的頻譜圖，主瓣寬度太寬，這樣就保留了部分的高頻分量 [10,13]。剛開始選擇這個題目，是因為數(shù)字信號和 Matlab 都學(xué)過，所以有一定的基礎(chǔ)感覺會簡單一點，而且自己對于圖像處理方面挺感興趣。在學(xué)習(xí)中，熟練掌握一門語言也很重要，比如本論文是基于 Matlab 來進行數(shù)字濾波器的設(shè)計，因為 Matlab 是同類軟件中最出色的。雖然論文完成了，但是還是有很多不足之處，因為我只是簡單的用 Matlab 分析了IIR 和 FIR 濾波器，并且用程序去實現(xiàn)，并不能真正運用到現(xiàn)實生活中?；?Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)31致 謝在這忙碌的幾個月中，感謝譚談老師不辭勞苦的督促和教導(dǎo)，讓自己剛開始對論文一頭霧水開始有了清晰地思路與方向。其次，我要感謝在我碰到困難時，鼓勵我和幫助過我的同學(xué)，他們不僅給了我精神上的鼓勵還在論文的寫作上給予了很大幫助。為我以后的生活、工作和學(xué)習(xí)都提供了很大的幫助。因此，萬分感謝所有老師的付出。同時，感謝學(xué)校圖書館，給了我良好的學(xué)習(xí)環(huán)境和豐富的知識庫，讓我能盡情的在知識的海洋里面尋找屬于自己的那片天地，讓我忘卻外界的紛紛擾擾。接下來，自己認(rèn)為最好的設(shè)計方向是把各種濾波器的方式設(shè)計出來并組合在一起連成一個整體，在這樣的基礎(chǔ)上設(shè)計一個界面，來真正的做成一個比較完善整體的設(shè)計。在論文寫作過程中對比的重要性很突出，通過對 IIR 和 FIR 對比的過程，發(fā)現(xiàn)不同方式之間的相似點，不同點，然后得出各自的優(yōu)缺點，這些對現(xiàn)實中投入使用，根據(jù)外界條件選擇最優(yōu)選的方式提供了很好的基礎(chǔ)。所以我放棄直接著手論文，重新拿起書來把理論知識弄透徹，這樣基礎(chǔ)打的牢固，才有信心和能力進行論文的寫作?；?Matlab 的數(shù)字濾波器的設(shè)計與開發(fā)30總結(jié)一眨眼，幾個月的課程設(shè)計接近了尾聲。然后我們利用上面的巴特沃斯和布拉克曼