基于Matlabatlab的FIRFIR數字濾波器設計

郭天賜

（河南大學歐亞國際學院，河南　開封　47755000011）

基于Matlabatlab的FIRFIR數字濾波器設計

郭天賜

（河南大學歐亞國際學院，河南開封47755000011）

FIR數字濾波器是數字信號處理的一個重要組成部分，由于FIR數字濾波器具有嚴格的線性相位，因此，在信息采集和處理過程中得到了廣泛的應用。本文介紹了FIR數字濾波器的概念，并在分析三種設計方法原理的基礎上，借助Matlab仿真軟件工具分別實現窗函數法、頻率采樣法和等波紋逼近法設計FIR濾波器。

數字濾波器；FIR濾波器；設計

現實中我們經常接觸的信號有很多種。對于信號的各個應用，濾波器是必不可少的一種工具。選用合適的濾波器對產品的質量有一定的提升作用。同樣，不適合的濾波器將影響產品的使用。

1　FIR濾波器及其設計方法

1.1FIR濾波器的基本概念

FIR濾波器即沖激響應是有限單位長的濾波器，是一個穩定的系統。因此，FIR濾波器在生活的各個方面都有廣泛的應用。

1.1.1FIR濾波器的優點

①對于嚴格的線性相位要求能簡單地完成，同時被處理的信號產生可減少相位失真；

②可得到多帶幅頻特性；

③極點全部在原點（永遠穩定），無穩定性問題；

④總是滿足因果性；

⑤無反饋運算，運算誤差小。

1.1.2FIR濾波器的缺點

①因為無極點，要獲得好的過渡帶特性，須以較高的階數為代價；

②需要計算機輔助設計以使程序完成。

1.2數字濾波器的設計方法

FIR數字濾波器有很多設計的方式，其中較為常用的是頻率采樣設計法、窗最優化設計法和函數設計法。本文主要論述利用窗函數法、頻率采樣法和等波紋切比雪夫逼近法來實現FIR濾波器的多種設計。

2　FIR數字濾波器設計

FIR數字濾波器設計過程一般包括以下基本問題：

①根據實際要求確定數字濾波器性能指標；

②用一個因果穩定的系統函數去逼近這個理想性能指標。

2.1窗函數法

窗函數法從時域出發，傳輸函數的系數由截取一段理想濾波器的單位脈沖響應來表示，被稱為設計FIR數字濾波器的最簡單的方法。由于窗函數法設計FIR數字濾波器是在時域進行的，且理想單位脈沖響應一般是非因果、無限長的，這就需要截取，截取的方法就是加窗。Matlab工具箱中加入了許多的窗函數，其中三角窗的最小阻帶值為25db、主瓣寬度為8π/M、精確過度帶寬6.1π/M，矩形窗最小阻帶值為21db、主瓣寬度為4π/M、精確過度帶寬1.8π/M，布拉克曼窗最小阻帶值為74db、主瓣寬度為12π/M、精確過度帶寬為11π/M，切比雪夫窗最小阻帶值為21db、主瓣寬度為4π/M、精確過度帶寬為1.1π/M，漢寧窗最小阻帶值為44db、主瓣寬度為8π/M、精確過度帶寬為6.6π/M，海明窗最小阻帶值為53db、主瓣寬度為4π/M、精確過度帶寬為1.8π/M，凱塞窗最小阻帶值為21db、主瓣寬度為4π/M、精確過度帶寬為1.8π/M［1-2］。由頻率相應Hd（ejw）求取沖擊響應，即分析出系統的穩定性和耦合性，然后通過在Matlab上編寫軟件程序，從而設計出一個近似的頻率響應，實現了硬件濾波器軟件的設計。設計過程如下：

為了把理想濾波器的無限長脈沖響應hd（n）乘以窗函數w（n）來產生一個被截斷的脈沖響應，人們使用了加窗的方法，即h（n）=hd（n）w（n），并且對頻率響應進行平滑。

窗函數主要用來減少序列因截斷而產生的Gibbs效應。為了克服這種現象，應該在設計窗函數時，濾波器的頻率特性最大程度地降低主瓣寬度的數值，且能量要求盡可能不被浪費；窗函數頻率特性中的ω越接近π其能量越接近零點［3］。

2.2頻率采樣法

頻率采樣法從頻域出發，理想濾波器的頻率響應為H（ejw），進行相同間隔的抽樣，將得出的理想頻率響應通過采樣的方式得到hd（k）：

再利用Hd（k）可求得FIR濾波器的系統函數H（Z）及頻率響應H（ejw），通過離散傅立葉反變換。而在各采樣點間的頻率響應則是其加權內插函數延伸疊加的結果。但在現實中信號大多是無限長的序列，插值逼近理論在實際應用中會有理想的頻率響應曲線和理想的頻率響應特性，內插函數的值越接近理想值，誤差越小，逐步轉變形狀的近似誤差，誤差的大小，為此設計時需要隨機地在不連續的區域添加不等量過度點。為了提高逼近的質量，可以通過在頻率相應的過渡帶內插入比較連續的采樣點，從而使得通帶和阻帶之間變換比較緩慢，以達到減少逼近誤差的目的。選取w∈［0，2π］內N個采樣點的約束條件為：

2.3等波紋逼近法

2.3.1等波紋切比雪夫逼近準則

對于通帶和阻帶的誤差數值在濾波器的設計中有不同的要求。等波紋切比雪夫逼近準則通過對通帶應用加權函數A，對阻帶使用應用加權函數B，以此來實現相同的在不同頻段的（通常指的是通帶和阻帶）加權誤差最大值，從而實現濾波器最大誤差在滿足性能指標的條件下達到最小值。

窗函數法與頻率采樣法在FIR數字濾波器的設計中被人們廣泛地應用，但兩者均無法成為最優化的設計，被廣泛應用的原因是因為考慮到成本與使用的問題。通常情況下線性相位濾波在不同的頻帶內逼近的最大容許誤差要求不同。即使得Hd（ejw）和H（ejw）之間的最大絕對誤差最小［4］。

等波紋切比雪夫逼近是采用加權逼近誤差E（ejw），它可以表示為：

其中，W（ejw）為逼近誤差加權函數在誤差要求高的頻段上，可以取較大的加權值，否則，應當取較小的加權值。

雖然按照FIR數字濾波器單位取樣響應h（n）的對稱性和采樣點數的奇偶性，FIR濾波器可被分為不同類型的幾種濾波器，但濾波器的頻率響應的公式卻是可以寫成同樣的式子：

其中，k∈｛0，1｝，H（ω）為幅度函數，且是一個純實數，表達式也可以寫成統一的形式：

其中，Q（ω）為ω的固定函數，P（ω）為M個余弦函數的線性組合［5］。

2.3.2仿真函數

利用remezord和remez函數可以實現FIRDF的最優化設計，而且remezord和remez函數存在于數字信號處理工具箱中：

（1）n，f0，a0，weight=remezord（f，a，dev）

是用來代替remezord中的N，歸一化頻帶邊緣、頻帶的振幅頻帶振幅中存在的功能和在每個頻帶之內權重的加權系數的輸入參數，它們代表不同的含義。

（2）h=reme（n，f0，a0，weight，′ftype′）

功能：利用remez函數可以得到最優化設計的FIR DF的h（n）系數，輸入參數中濾波器的階數由第一個數據n來表示。為了區分濾波器的類型，人們在最后加入ftype值來區分不同類型的濾波器。remezord函數的階數n提升的高估或低估。因此在通過remezoed函數求出濾波器的系數后，對于阻帶最小衰減的檢測觀察其是否滿足設計要求［6］。如果此時的阻帶最小衰減不能滿足設計要求，那么就必須提高濾波器的階數。

3　結果分析

在間斷點處增加一個過渡點后的情況。可以看出濾波器的帶外衰減指標有了明顯的改善，但同時增加了濾波器的過渡帶寬。

采樣點上2/N點集中在指定的通帶和阻帶截止頻率的整數倍的頻率，這種方法受到一定的限制，相對而言，不夠靈活。充分地提高N，對于任何給定的頻率可以接近。但計算量和復雜性增加。頻率采樣法明顯偏離設計指標，且相對來說更不易操作，只有過渡帶樣點值與實際相對出入較小時，衰耗特性才會取得令人滿意的結果。

4　結論

在設計中，要使設計的濾波器的沖激響應同時滿足預先設定的全部參數是十分困難的，如果該濾波器的特性不滿足要求，但是原有參數必須做適當調整。這在程序中很容易實現，只需對參數M、δ1、δ2事先給定，允許ωp與ωs隨著要求的改變而變化，就可以得到新條件下濾波器的特性。從而達到節約資源的效果。這樣使得設計出來的濾波器更為簡單經濟。

相對于前兩種濾波器的設計方法，等波紋逼近法有著無法忽視的優勢。其不僅程序簡單而且在實現中的應用層面也相對較寬廣。當一個產品需要壓縮程序空間時，對于濾波器設計模塊最優化設計給出的濾波器必定是我們的首選方式。

［1］郭士劍，賀志國.Matlab7.x數字信號處理［M］.北京：人民郵電出版社，2008.

［2］程耕國，陳華麗.信號與系統實驗教程（Matlab）版［M］.北京：機械工業出版社，2009.

［3］謝仕宏.Matlab R2008控制系統動態仿真實例教程［M］.北京：化學工業出版社，2005.

［4］徐昕，伯曉晨，等.Matlab工具箱應用指南——控制工程篇［M］.北京：電子工業出版社，2000.

［5］程耕國.信號與系統：上冊［M］.北京：機械工業出版社，2009.

［6］程耕國.信號與系統：下冊［M］.北京：機械工業出版社，2009.

Design of FIR Digital Filter based on Matlab

Guo Tianci
（Eurasian International College，Henan University，Kaifeng Henan 475001）

FIR digital filter is an important part of digital signal processing，due to its strict linear phase，so in the process of information acquisition and processing，it has been widely used.The concept of FIR digital filter is introduced in this paper，and based on the analysis of three design principles and methods，with the tool of MATLAB simulation software，FIR digital filter is designed by the window function method and frequency sampling method and the ripple approximation method respectively.

digital filter；FIR filter；design

TN713.7

A

1003-5168（2015）10-0042-3

2015-9-20

郭天賜（1995-），男，本科，研究方向：通信工程。