FIR數(shù)字濾波器設(shè)計(jì)

FIR數(shù)字濾波器設(shè)計(jì)第一頁(yè)，共103頁(yè)。第7章FIR數(shù)字濾波器設(shè)計(jì)7.1FIRDF設(shè)計(jì)的窗函數(shù)法7.2窗函數(shù)7.3FIRDF設(shè)計(jì)的頻率抽樣法7.4FIRDF設(shè)計(jì)的切比雪夫最佳一致逼近法7.5幾種簡(jiǎn)單形式的濾波器7.6簡(jiǎn)單整系數(shù)濾波器7.7差分濾波器第二頁(yè)，共103頁(yè)。IIR數(shù)字濾波器：有極點(diǎn)，也有零點(diǎn)，因此可以借用經(jīng)典的連續(xù)濾波器的設(shè)計(jì)方法，且取得非常好的效果，如好的衰減特性，準(zhǔn)確的邊緣頻率。由于FIR數(shù)字濾波器只有零點(diǎn)而沒(méi)有極點(diǎn)，所以沒(méi)辦法借用連續(xù)濾波器的設(shè)計(jì)方法。其思路是：直接從頻域出發(fā)，即以某種準(zhǔn)則逼近理想的頻率特性，且保證濾波器具有線性相位。第三頁(yè)，共103頁(yè)。7.1Fourier級(jí)數(shù)法（窗函數(shù)法）1.由理想的頻率響應(yīng)得到理的；2.由得到因果、有限長(zhǎng)的單位抽樣響應(yīng);3.對(duì)加窗得到較好的頻率響應(yīng)。理想頻率響應(yīng)一、思路與方法：第四頁(yè)，共103頁(yè)。設(shè)理想低通濾波器的幅頻為1，相頻為零：則：特點(diǎn)：無(wú)限長(zhǎng)非因果偶對(duì)稱解決方法：截短，移位保留即：隱含著使用了窗函數(shù)第五頁(yè)，共103頁(yè)。于是：注意：H(Z)是因果的，且是線性相位的，即？即事先給一線性相位為了省去每次的移位，可以令：

在通帶內(nèi)這樣：第六頁(yè)，共103頁(yè)。于是：使用了矩形窗上式的的表達(dá)式及設(shè)計(jì)的思路可推廣到高通、帶阻及帶通濾波器，也可推廣到其它特殊類型的濾波器。實(shí)際上，給定一個(gè) ，只要能積分得到，即可由截短、移位的方法得到因果的、且具有線性相位的FIR濾波器。第七頁(yè)，共103頁(yè)。高通：令：相當(dāng)于用一個(gè)截止頻率在處的低通濾波器（實(shí)際上是全通濾波器）減去一個(gè)截止頻率在處的低通濾波器。第八頁(yè)，共103頁(yè)。令：相當(dāng)于用一個(gè)截止頻率在處的低通濾波器減去一個(gè)截止頻率在處的低通濾波器。帶通：第九頁(yè)，共103頁(yè)。令：

：窗函數(shù)，自然截短即是矩形窗。 當(dāng)然也可以用其它形式的窗函數(shù)。相當(dāng)于？帶阻：第十頁(yè)，共103頁(yè)。例1.設(shè)計(jì)低通FIR

DF,令歸一化截止頻率

＝0.125,M＝10，20，40，用矩形窗截短。結(jié)果如右圖第十一頁(yè)，共103頁(yè)。接上例：M＝10分別用矩形窗和Hamming窗使用Hamming窗后，阻帶衰減變好，但過(guò)渡帶變寬。第十二頁(yè)，共103頁(yè)。例：理想差分器及其設(shè)計(jì)令：理想差分器的頻率特性：理想微分器的頻率特性：第十三頁(yè)，共103頁(yè)。奇對(duì)稱，純虛函數(shù)理想相頻特性第2類FIR濾波器第十四頁(yè)，共103頁(yè)。實(shí)際相頻特性有關(guān)各種差分器的性能，本章將繼續(xù)討論幅頻：1矩形窗2哈明窗第十五頁(yè)，共103頁(yè)。例：設(shè)計(jì)Hilbert變換器

第2類FIR濾波器思考：能否用上一章的方法設(shè)計(jì)差分器和Hilbert變換器？第十六頁(yè)，共103頁(yè)。優(yōu)點(diǎn)：1.無(wú)穩(wěn)定性問(wèn)題；2.容易做到線性相位；3.可以設(shè)計(jì)各種特殊類型的濾波器；4.方法特別簡(jiǎn)單。缺點(diǎn)：1.不易控制邊緣頻率；2.幅頻性能不理想；3.較長(zhǎng)；二、FIRDF設(shè)計(jì)的窗函數(shù)法的特點(diǎn)：改進(jìn)：1.使用其它類型的窗函數(shù)；

2.改進(jìn)設(shè)計(jì)方法。第十七頁(yè)，共103頁(yè)。三、關(guān)于對(duì)截短的討論誤差曲線誤差能量第十八頁(yè)，共103頁(yè)。請(qǐng)自己導(dǎo)出此式什么情況下為最?。康谑彭?yè)，共103頁(yè)。最小所以，有限項(xiàng)傅立葉級(jí)數(shù)是在最小平方意義上對(duì)原信號(hào)的逼近。傅立葉級(jí)數(shù)是正交變換，這也體現(xiàn)了正交變換的性質(zhì)。第二十頁(yè)，共103頁(yè)。窗函數(shù)法周期信號(hào)展開為傅里葉級(jí)數(shù)傅里葉系數(shù)傅里葉級(jí)數(shù)法同一事情不同名稱第二十一頁(yè)，共103頁(yè)。7.2窗函數(shù)窗函數(shù)的使用在數(shù)字信號(hào)處理中是不可避免的。數(shù)據(jù)、頻譜、自相關(guān)函數(shù)等都需要截短。對(duì)窗函數(shù)提出那幾方面的要求？關(guān)鍵是要搞清楚使用窗函數(shù)后所產(chǎn)生的影響：一個(gè)域相乘，在另一個(gè)域是卷積。矩形窗第二十二頁(yè)，共103頁(yè)。2.邊瓣最大峰值 （dB）3.邊瓣譜峰衰減速度 （dB/oct）越小越好！越小越好！越大越好！對(duì)窗函數(shù)的技術(shù)要求：1.3dB帶寬：主瓣歸一化幅度降到－3dB時(shí)的帶寬；或直接用。令則的單位為；第二十三頁(yè)，共103頁(yè)。常用窗函數(shù)：1.矩形窗2.三角窗Bartlett窗3.漢寧窗Hanning4.漢明窗Hamming第二十四頁(yè)，共103頁(yè)。窗函數(shù)第二十五頁(yè)，共103頁(yè)。窗函數(shù)第二十六頁(yè)，共103頁(yè)。7.3FIRDF設(shè)計(jì)的頻率抽樣法窗函數(shù)法：給定連續(xù)的理想的，用得到因果的、具有線性相位的FIRDF逼近第二十七頁(yè)，共103頁(yè)。離散化直接賦值計(jì)算更容易思考：相等？濾波器已設(shè)計(jì)出第二十八頁(yè)，共103頁(yè)?？芍付ǎ喝绾沃付ǎ康诙彭?yè)，共103頁(yè)。轉(zhuǎn)移函數(shù)、頻率響應(yīng)和給定的 的關(guān)系：用DFT系數(shù)作為權(quán)函數(shù)來(lái)表示設(shè)計(jì)出的有何特點(diǎn)第三十頁(yè)，共103頁(yè)。原抽樣再抽樣關(guān)系？第三十一頁(yè)，共103頁(yè)。用插值的方法得到所要的濾波器：插值函數(shù)權(quán)重線性相位應(yīng)為實(shí)數(shù)所以：如何指定？第三十二頁(yè)，共103頁(yè)。

為偶數(shù)：

為奇數(shù)：其它賦值方法見(jiàn)書。當(dāng)然，阻帶內(nèi)應(yīng)指定為零。另外，為了得到好的幅頻響應(yīng)，在1和0之間加過(guò)渡點(diǎn)，如0.5。第三十三頁(yè)，共103頁(yè)。7.4用Chebyshev最佳一致逼近設(shè)計(jì) FIRDF7.4.1最佳一致逼近定理7.4.2利用最佳一致逼近理論設(shè)計(jì)FIRDF

7.4.3關(guān)于誤差函數(shù)的極值特性7.4.4FIRDF的四種表示形式7.4.5設(shè)計(jì)舉例7.4.6濾波器階次估計(jì)第三十四頁(yè)，共103頁(yè)。上述兩種方法（窗函數(shù)法和頻率抽樣法）設(shè)計(jì)的FIRDF的頻率響應(yīng)都不理想，即通帶不夠平，阻帶衰減不夠大，過(guò)渡帶過(guò)寬，頻率邊緣不能精確指定。因此我們要尋找新的設(shè)計(jì)方法。此方法即是Chebyshev最佳一致逼近法。該方法在數(shù)字信號(hào)處理中占有重要的定位，是設(shè)計(jì)

FIRDF最理想的方法。但是，該方法的原理稍為復(fù)雜。第三十五頁(yè)，共103頁(yè)。給定理想的，設(shè)計(jì)，使是對(duì)的“最佳”逼近。1.最小平方逼近：著眼積分區(qū)間內(nèi)整體誤差最小。傅立葉級(jí)數(shù)法即是這一種。對(duì)函數(shù)f(x)逼近的方法：目標(biāo)：？第三十六頁(yè)，共103頁(yè)。插值法：尋找階多項(xiàng)式，使其在個(gè)點(diǎn)上滿足：頻率抽樣方法3.最佳一致逼近：尋找，使誤差在區(qū)間[a,b]均勻一致，并使誤差的極大值

達(dá)到最小。第三十七頁(yè)，共103頁(yè)。Chebyshev最佳一致逼近理論解決了的存在性、唯一性及構(gòu)造方法等問(wèn)題。將最佳一致逼近理論應(yīng)用于FIRDF的設(shè)計(jì)，是數(shù)學(xué)和信號(hào)處理理論相結(jié)合的又一典型范例。該方法可以設(shè)計(jì)出性能優(yōu)良的FIRDF，是FIR設(shè)計(jì)的主要方法。該方法又稱McClellan----Parks方法第三十八頁(yè)，共103頁(yè)。一、切比雪夫最佳一致逼近定理在階多項(xiàng)式的集合中，尋找多項(xiàng)式使其相對(duì)其它所有的多項(xiàng)式對(duì)的偏差為最?。鹤钚∽畲笤淼谌彭?yè)，共103頁(yè)。交錯(cuò)點(diǎn)組原理：令：誤差最大值誤差曲線是最佳一致逼近的充要條件是，在上至少存在個(gè)交錯(cuò)點(diǎn)使得：極值點(diǎn)交錯(cuò)點(diǎn)第四十頁(yè)，共103頁(yè)。所以：是的極值點(diǎn)，它們構(gòu)成了一個(gè)“交錯(cuò)點(diǎn)組”什么樣的函數(shù)（或多項(xiàng)式）可以充當(dāng)誤差曲線？

Chebyshev多項(xiàng)式：在區(qū)間[-1,1]上存在個(gè)點(diǎn)：第四十一頁(yè)，共103頁(yè)。輪流使取極值＋1，－1。是的階多項(xiàng)式，最高項(xiàng)系數(shù)是，在所有階多項(xiàng)式的集合中，和0的偏差為最小。因此，可用為誤差多項(xiàng)式。：個(gè)極值點(diǎn)交錯(cuò)點(diǎn)組定 理要求：個(gè)極值點(diǎn)？第四十二頁(yè)，共103頁(yè)。二、利用最佳一致逼近理論設(shè)計(jì)FIRDF理想濾波器要設(shè)計(jì)的濾波器第四十三頁(yè)，共103頁(yè)。四種情況下的“濾波器增益”都是實(shí)函數(shù)，也有四種表示形式。其一是：線性相位FIR濾波器有四種形式：我們用逼近理想濾波器。顯然，若能求出，則濾波器也就設(shè)計(jì)出來(lái)了。第四十四頁(yè)，共103頁(yè)。用一致逼近目標(biāo)使誤差曲線尋找的最大值為最小，并呈現(xiàn)交錯(cuò)。第四十五頁(yè)，共103頁(yè)。定義加權(quán)函數(shù)：在設(shè)計(jì)濾波器時(shí)，對(duì)通帶和阻帶往往有不同的要求，如通帶要求特別平，這是需要犧牲阻帶；反之，要想阻帶衰減特別大，則需要犧牲通帶。實(shí)現(xiàn)方法：給以不同的加權(quán)。需要離散化第四十六頁(yè)，共103頁(yè)。由交錯(cuò)點(diǎn)組定理：注意，將頻率分成了個(gè)離散的點(diǎn)。分點(diǎn)在通帶和阻帶上，過(guò)渡帶不考慮。目的是取得個(gè)極值點(diǎn)。寫成矩陣形式第四十七頁(yè)，共103頁(yè)。方陣，可唯一地求出然而，該方程的求解異常困難！第四十八頁(yè)，共103頁(yè)。McClellan.J.H&Parks.T.W等于70年代初提出用數(shù)值分析中的Remez算法，靠一次次的迭代來(lái)求解最優(yōu)的系數(shù)及。從而達(dá)到濾波器設(shè)計(jì)的目的。該方法不但可以用來(lái)設(shè)計(jì)低通、高通、帶通、帶阻等經(jīng)典濾波器，而且可以用來(lái)設(shè)計(jì)差分濾波器，Hilbert變換器。不但可以給出好的幅頻特性、線性相位，而且可以給出較為準(zhǔn)確定邊緣頻率。數(shù)字信號(hào)處理中最有名的算法之一！第四十九頁(yè)，共103頁(yè)。Step1.先在通帶、阻帶頻率軸上等間隔取M＋2個(gè)頻率點(diǎn)，計(jì)算出。它是相對(duì)第一次指定的交錯(cuò)點(diǎn)組產(chǎn)生的誤差A(yù).第五十頁(yè)，共103頁(yè)。求出后，利用插值公式，在不知 的情況下求出。B.第五十一頁(yè)，共103頁(yè)。當(dāng)然，初次求出的肯定不是最優(yōu)的！將求出的代入C.可求出誤差函數(shù)。如果第一次迭代即是最優(yōu)，那么應(yīng)是 的極值點(diǎn)。當(dāng)然，一次迭代是不夠的。完成第一次迭代！第五十二頁(yè)，共103頁(yè)。Step2.檢查是否有的頻率點(diǎn)（肯定有）。將出現(xiàn)這種情況的頻率點(diǎn)和原來(lái)指定的頻率點(diǎn)中相距最近的點(diǎn)相交換（注意：這樣的點(diǎn)可能不止一個(gè)），這樣，就得到一組新的頻率點(diǎn)組，當(dāng)然，它們不再是原頻率區(qū)間的等分。第五十三頁(yè)，共103頁(yè)。Step3.將新的頻率點(diǎn)組再重復(fù)步驟2，又可得到一組新的交錯(cuò)點(diǎn)組：代入公式A.求出新的再代入公式求出新的B.再代入公式求出新的C.重復(fù)了步驟1第五十四頁(yè)，共103頁(yè)。如此重復(fù)迭代，每一次都是把新的局部極值點(diǎn)當(dāng)作新的交錯(cuò)點(diǎn)組，所以，每一次的都是遞增的，最后收斂到自己的上限。再迭代一次，也不會(huì)再增加，頻率點(diǎn)組也不會(huì)再移動(dòng)，這時(shí)的即是對(duì)的最佳一致逼近。Step4.將最優(yōu)的配上線性相位，作傅立葉反變換，即可得所設(shè)計(jì)濾波器的。第五十五頁(yè)，共103頁(yè)。通帶內(nèi)的峰值偏差最佳一致逼近是在通帶與阻帶內(nèi)進(jìn)行的，過(guò)渡帶沒(méi)有考慮。迭代步驟是阻帶峰值偏差；第五十六頁(yè)，共103頁(yè)。三、關(guān)于誤差函數(shù)的極值特性（見(jiàn)書）四、FIRDF的四種表示形式第五十七頁(yè)，共103頁(yè)。把上述四種形式稍作改造，得到如下的統(tǒng)一形式，目的是便于編程：第五十八頁(yè)，共103頁(yè)。例1：設(shè)計(jì)低通FIRDF:調(diào)整通帶、阻帶的加權(quán)及濾波器的長(zhǎng)度。設(shè)計(jì)結(jié)果五、設(shè)計(jì)舉例第五十九頁(yè)，共103頁(yè)。參數(shù)調(diào)整對(duì)濾波器性能的影響：第六十頁(yè)，共103頁(yè)。例2：設(shè)計(jì)多帶濾波器，抽樣頻率500Hz,在 50Hz、100Hz及150Hz處陷波。通帶加權(quán)為8，阻帶為1-17dB通帶、阻帶加權(quán)都是1-25dB第六十一頁(yè)，共103頁(yè)。六、階次估計(jì)設(shè)計(jì)濾波器之前，濾波器的長(zhǎng)度（即階次）是未知道。顯然，要求：通帶越平，阻帶衰減越大，過(guò)渡帶越窄，濾波器的階次越高。：通帶紋波：阻帶紋波經(jīng)驗(yàn)公式第六十二頁(yè)，共103頁(yè)。另一估計(jì)公式：估計(jì)出的階次稍低例如，對(duì)例1的第一種情況：求出：和原來(lái)給定的相同第六十三頁(yè)，共103頁(yè)。7.5幾種簡(jiǎn)單形式的濾波器一、平均濾波器二、平滑濾波器三、梳狀濾波器這一類濾波器性能不是很好，但濾波器簡(jiǎn)單，有時(shí)很實(shí)用，有的具有一些特殊的用途。第六十四頁(yè)，共103頁(yè)。信噪比(SNR)與噪聲減少比（NRR）信噪比：觀察信號(hào)信號(hào)噪聲為了減少噪聲，將通過(guò)一個(gè)濾波器第六十五頁(yè)，共103頁(yè)。噪聲減少比（NoiseReductionRation,NRR):越小越好！可以證明：第六十六頁(yè)，共103頁(yè)。一、平均濾波器點(diǎn)平均器IIR系統(tǒng)第六十七頁(yè)，共103頁(yè)。第六十八頁(yè)，共103頁(yè)。第六十九頁(yè)，共103頁(yè)?？梢郧蟪觯嚎梢?jiàn)N足夠大，即可就可以獲得足夠小的NRR。但是，N過(guò)大會(huì)使濾波器具有過(guò)大的延遲：群延遲=（N－1）/2而且會(huì)使其主瓣的單邊的帶寬大大降低，這就有可能在濾波時(shí)使有用的信號(hào)s(n)也受到損失。因此，在平均器中，N不宜取得過(guò)大。第七十頁(yè)，共103頁(yè)。二、平滑濾波器Savitzky－Golay平滑器：基于多項(xiàng)式擬合的方法，具體推導(dǎo)過(guò)程見(jiàn)教材。5點(diǎn)2次（拋物線）擬合：7點(diǎn)3次擬合：在NRR和階次N之間取得折中。MATLAB文件： sgolay.m第七十一頁(yè)，共103頁(yè)。三、梳狀濾波器作用：去除周期性的噪聲，或是增強(qiáng)周期 性的信號(hào)分量。

第七十二頁(yè)，共103頁(yè)。第七十三頁(yè)，共103頁(yè)。第七十四頁(yè)，共103頁(yè)。7.6建立在極零抵消基礎(chǔ)上的 簡(jiǎn)單整系數(shù)的濾波器對(duì)信號(hào)作實(shí)時(shí)濾波處理時(shí)，有時(shí)對(duì)濾波器的性能要求并不很高，但要求計(jì)算速度快，濾波器的設(shè)計(jì)也應(yīng)簡(jiǎn)單易行，因而希望濾波器的系數(shù)為整數(shù)。特別是當(dāng)用匯編語(yǔ)言編寫程序時(shí)，更希望如此。采用極零抵消的方法，可以設(shè)計(jì)出簡(jiǎn)單整系數(shù)的低通、高通、帶通和帶阻濾波器。第七十五頁(yè)，共103頁(yè)。1.低通單位圓上均勻分布M個(gè)零點(diǎn)設(shè)置一極點(diǎn)，抵消掉z=1處零點(diǎn)M點(diǎn)平均器第七十六頁(yè)，共103頁(yè)。2.高通單位圓上均勻分布M個(gè)零點(diǎn)設(shè)置一極點(diǎn)，抵消掉z=－1處零點(diǎn)第七十七頁(yè)，共103頁(yè)。上述低通和高通濾波器的系數(shù)都是整系數(shù)（系數(shù)1/N可最后單獨(dú)處理），如果認(rèn)為幅頻響應(yīng)不滿意，可以取濾波器系數(shù)仍為整數(shù)第七十八頁(yè)，共103頁(yè)。3.帶通實(shí)際應(yīng)用第七十九頁(yè)，共103頁(yè)。為保證分母取整數(shù)，要求取整數(shù)因此：在要求整系數(shù)的情況下，對(duì)帶通濾波器，其通帶的中心頻率收到限制。第八十頁(yè)，共103頁(yè)。4.帶阻設(shè)計(jì)方法幅頻：全通幅頻－帶通幅頻相頻：配置相頻例令，設(shè)計(jì)50Hz陷波器，中心頻率范圍在解：取第八十一頁(yè)，共103頁(yè)。由于因此增加一對(duì)共軛極點(diǎn)：150Hz現(xiàn)在需要確定M:?第八十二頁(yè)，共103頁(yè)。具有相同相位第八十三頁(yè)，共103頁(yè)。第八十四頁(yè)，共103頁(yè)。7.7低階低通差分濾波器？理想微分器：理想差分器：第八十五頁(yè)，共103頁(yè)。為了防止在高頻端將噪聲放大，取：低通差分器：差分器的一般形式：第八十六頁(yè)，共103頁(yè)。差分器的抽樣響應(yīng)：所以，差分器是奇對(duì)稱的。現(xiàn)在的任務(wù)是確定系數(shù)兩點(diǎn)中心差分：第八十七頁(yè)，共103頁(yè)?！白罴选辈罘制鞅平`差：得到最佳系數(shù)得到最佳通帶第八十八頁(yè)，共103頁(yè)。最佳通帶：可求出：M＝2可求出：M＝3第八十九頁(yè)，共103頁(yè)。M＝2和3時(shí)“最佳”差分器的幅頻特性：但是，上述“最佳”差分器的系數(shù)全是小數(shù)，我們希望得到整系數(shù)。實(shí)際上，人們從不同的角度，已給出了不同形式的整系數(shù)差分器。后來(lái)，人們還導(dǎo)出了“次最佳”的整系數(shù)差分器。第九十頁(yè)，共103頁(yè)。單純M次差分；牛頓－柯斯特差分；Lanczos差分（多項(xiàng)式擬合）；平滑化差分；最佳差分；整系數(shù)比較參數(shù)：第九十一頁(yè)，共103頁(yè)。7.8濾波器設(shè)計(jì)小結(jié)IIR濾波器的優(yōu)點(diǎn)：1.好的通帶與阻帶衰減;準(zhǔn)確的通帶與阻帶邊緣頻率;2.濾波時(shí)需要的計(jì)算量較少缺點(diǎn):不具有線性相位，有可能存在穩(wěn)定性問(wèn)題。FIR濾波器的優(yōu)點(diǎn)：1.可取得線性相位;2.無(wú)穩(wěn)定性問(wèn)題;缺點(diǎn)：濾波時(shí)需要的計(jì)算量較少第九十二頁(yè)，共103頁(yè)。FIR窗函數(shù)法頻率抽樣法一致逼近法簡(jiǎn)單平均簡(jiǎn)單平滑設(shè)計(jì)方法簡(jiǎn)單，性能不夠好性能非常好簡(jiǎn)單，實(shí)用，性能不夠好IIR梳狀濾波器極零抵消濾波器特殊用途，周期性簡(jiǎn)單實(shí)用，速度快第九十三頁(yè)，共103頁(yè)。與本章內(nèi)容有關(guān)的MATLAB文件：產(chǎn)生窗函數(shù)的文件有八個(gè)：bartlett（三角窗）;2.blackman（布萊克曼窗）;3.boxcar（矩形窗）;4.hamming（哈明窗）;5.hanning（漢寧窗）;6.triang（三角窗）;7.chebwin（切比雪夫窗）;8.kaiser（凱賽窗）;兩端為零兩端不為零調(diào)用方式都非常簡(jiǎn)單請(qǐng)見(jiàn)help文件稍為復(fù)雜第九十四頁(yè)，共103頁(yè)。9．fir1.m用“窗函數(shù)法”設(shè)計(jì)FIRDF。調(diào)用格式：（1）b=fir1(N,Wn);(2)b=fir1(N,Wn,‘high’);(3)b=fir1(N,Wn,‘stop’);N：階次，濾波器長(zhǎng)度為N＋1；Wn：通帶截止頻率，其值在0～1之間，1對(duì)應(yīng)Fs/2b:濾波器系數(shù)。CUST第九十五頁(yè)，共103頁(yè)。對(duì)格式（1），若Wn為標(biāo)量，則設(shè)計(jì)低通濾波器，若Wn是1×2的向量，則用來(lái)設(shè)計(jì)帶通濾波器，若Wn是1×L的向量，則可用來(lái)設(shè)計(jì)L帶濾波器。這時(shí)，格式（1）要改為:b=fir1(N,Wn,'DC-1'),或b=fir1(N,Wn,'DC-0')前者保證第一個(gè)帶為通帶，后者保證第一個(gè)帶為阻帶。格式（2）用來(lái)設(shè)計(jì)高通濾波器，格式（3）用來(lái)設(shè)計(jì)帶阻濾波器。在上述所有格式中，若不指定窗函數(shù)的類型，fir1自

動(dòng)選擇Hamming窗。第九十六頁(yè)，共103頁(yè)。10．fir2.m本文件采用“窗函數(shù)法”設(shè)計(jì)具有任意幅頻相應(yīng)的FIR數(shù)字濾波器。其調(diào)用格式是：

b=fir1(N,F,M);F是頻率向量，其值在0～1之間，M是和F相對(duì)應(yīng)的所希望的幅頻相應(yīng)。如同fir1,缺省時(shí)自動(dòng)選用Hamming窗。例：設(shè)計(jì)一多帶濾波器，要求頻率在0.2~0.3, 0.6~0.8之間為1，其余處為零。設(shè)計(jì)結(jié)果如下：第九十七頁(yè)，共103頁(yè)。N=30,90時(shí)幅頻響應(yīng)響應(yīng)及理想幅頻響應(yīng)；N=30N=90第九十八頁(yè)，共