巴特沃斯低通濾波器的設(shè)計與仿真分析

0   引言

巴特沃斯濾波器（Butterworth Filter）也被稱作最大平坦濾波器。其特點(diǎn)是在通頻帶內(nèi)，其頻率響應(yīng)曲線最大限度平坦、單調(diào)遞減無波紋產(chǎn)生；而在阻頻帶內(nèi)，其頻率響應(yīng)曲線逐漸下降為零。這些優(yōu)點(diǎn)使其在信號處理領(lǐng)域有著廣泛的應(yīng)用。對于濾波器的性能一般考慮其幅頻特性，但對于更深層次的信號處理應(yīng)用方面，濾波器的相頻特性也是其重要特性之一。尤其是在多點(diǎn)激勵、載荷建立以及傳遞路徑識別等方面問題的研究中具有重要作用。

1   巴特沃斯濾波器傳遞函數(shù)分析

連續(xù)時間巴特沃斯低通濾波器可用式（1）表示。

式中：B(jω) 為連續(xù)時間巴特沃斯傳遞函數(shù)，ω 為頻率，ω_C 為濾波器截止頻率，N 為階數(shù)，該濾波器Bode 圖如圖1 所示。

將式（1）轉(zhuǎn)換為Laplace 域分析，即令s = jω，可得到式（2）。

圖1 巴特沃斯濾波器Bode圖

求解式（2）的極點(diǎn)，可得式（3）。

對極點(diǎn)作歸一化處理，即令：

作歸一化巴特沃斯濾波器的拉普拉斯域平面，可以得到2N 個以虛軸Im 對稱的極點(diǎn)，以N = 3 為例其極點(diǎn)分布如圖2 所示，所有極點(diǎn)都均分在以原點(diǎn)為中心的單位圓上。

為形成因果穩(wěn)定系統(tǒng)，取左半平面極點(diǎn)構(gòu)建式（5）的傳遞函數(shù)。

圖2 極點(diǎn)分布圖

其中多項式系數(shù)b_i可構(gòu)為建歸一化N 階巴特沃斯濾波器的傳遞函數(shù)系數(shù)，不同階數(shù)系數(shù)如表1 所示。

對應(yīng)不同階數(shù)N， 繪制歸一化巴特沃斯濾波器Bode圖如圖3 所示。

圖3 歸一化巴特沃斯濾波器Bode圖

2   數(shù)字濾波器設(shè)計實例

現(xiàn)以表2 所述濾波器性能參數(shù)，實例設(shè)計巴特沃斯型數(shù)字濾波器。

針對通帶截止頻率與阻帶截止頻率設(shè)計要求，建立式（6）與（7）的濾波器設(shè)計約束條件。

通過約束條件，可通過式（8）得到理論設(shè)計階數(shù)。

由于階數(shù)需為整數(shù)，將理論階數(shù)N* 向高階取整以得到N，并代入式（9），即可得到截止頻率ω_c。

經(jīng)計算，N = 4，ω_c = 5.275 kHz，通過歸一化濾波器系數(shù)，即可表得到濾波器的傳遞函數(shù)為式（10）。

采用雙線性變換法，即通過式（11）完成s 域到z域的映射。由此即可得到在采樣頻率下的數(shù)字濾波器傳遞函數(shù)。

3   濾波器的仿真分析

仿真過程中，通過構(gòu)建兩路不同頻率正弦波信號x1和x2進(jìn)行疊加作為輸入，其中x1為300Hz，x2為15kHz。按圖4 所示原理圖，將疊加信號通過采樣器采樣后輸入數(shù)字濾波器，通過示波器分別顯示原始疊加信號、低頻信號與濾波信號的波形。

圖4 仿真系統(tǒng)原理

通過圖5的波形分析可以看出，該濾波器可以很好的將高頻信號濾除，濾波信號的波形基本能夠還原輸入的低頻信號。只是在幅值上有0.028%的衰減，相位上存在127.433μs的延遲。

圖5 波形分析

4   結(jié)束語

本文從巴特沃斯濾波器的設(shè)計原理著手，通過極點(diǎn)分析推導(dǎo)了巴特沃斯系數(shù)表的構(gòu)建方法，并分析了設(shè)計階數(shù)對濾波器的幅頻特性與相頻特性所產(chǎn)生的影響。通過結(jié)合實例指標(biāo)參數(shù)要求，使用雙線性變換法完成了巴特沃斯型數(shù)字濾波器的設(shè)計。最后通過仿真分析，驗證了該數(shù)字濾波器的濾波效果，結(jié)果表明在通頻段信號的幅值有微弱的降低，相位有少量滯后，但總體還原程度較高，效果理想。

參考文獻(xiàn)：

[1] [美]ALAN V O,等.信號與系統(tǒng)[M].3版.劉樹棠,譯.北京:電子工業(yè)出版社,2013.

[2] 龔作豪,沈君鳳.巴特沃斯低通濾波器的仿真設(shè)計[J].信息通信,2014(139):40-41.

[3] 仲帥.基于改進(jìn)粒子群算法的數(shù)字濾波器設(shè)計及應(yīng)用[D].長春:吉林大學(xué),2016.

（本文來源于《電子產(chǎn)品世界》雜志2022年12月期）