Elliptic濾波器設(shè)計(jì)詳解

本節(jié)將詳細(xì)講解橢圓函數(shù)濾波器(Elliptic Function Filter/Cauer Filter)的綜合設(shè)計(jì), 結(jié)合之前電路仿真的工作, 現(xiàn)將Elliptic濾波器集成到了濾波器設(shè)計(jì)App中, 提供一個(gè)試用版本, 并提供GitHub開(kāi)源鏈接。

橢圓函數(shù)濾波器特點(diǎn)和用途

橢圓函數(shù)濾波器(Elliptic Function Filter或Cauer Filter, 在微波射頻中有時(shí)也稱作Achieser-Zolotarev濾波器)是Cauer于1930-1931年首先使用橢圓有理函數(shù)(Elliptic Rational Function)逼近來(lái)進(jìn)行設(shè)計(jì)的濾波器, 其特點(diǎn)是通帶和阻帶內(nèi)都具有等紋波特性, 對(duì)比Chebyshev濾波器其更具陡峭的截止特性, 被廣泛用于實(shí)際工程應(yīng)用之中, 比如在DDS輸出級(jí)一般會(huì)加一個(gè)橢圓函數(shù)濾波器濾除鏡像頻率分量, 但是其缺點(diǎn)是對(duì)器件值要求非?？量?指精度和Q值), 很多濾波器理論上可行, 但是在實(shí)際中很難實(shí)現(xiàn)。

實(shí)際上橢圓函數(shù)濾波器的阻帶紋波為零則表現(xiàn)為切比雪夫?yàn)V波器(Chebyshev I Filter)特性, 當(dāng)通帶紋波為0則表現(xiàn)為逆切比雪夫?yàn)V波器(Chebyshev II Filter)特性, 當(dāng)通帶和阻帶紋波均為0則表現(xiàn)為巴特沃斯濾波器(Butterworth Filter)特性, 橢圓函數(shù)濾波器可以說(shuō)是濾波器的終極存在。

橢圓函數(shù)濾波器逼近

如何由通帶和阻帶等紋波推導(dǎo)出橢圓函數(shù)特征多項(xiàng)式(2)請(qǐng)參考《濾波器設(shè)計(jì)的逼近方法 - Butterworth, Chebyshev, Elliptic》。

橢圓函數(shù)濾波器傳遞函數(shù)

橢圓函數(shù)濾波器的幅頻響應(yīng)傳遞函數(shù)為:

橢圓函數(shù)濾波器傳遞函數(shù)是從其定義推導(dǎo)而來(lái), 即從通帶和阻帶都等紋波出發(fā), 推導(dǎo)出橢圓函數(shù)濾波器的特征函數(shù)表達(dá)式滿足如下關(guān)系：

式中是雅可比橢圓函數(shù),即通帶頻率和阻帶頻率的比值,即通帶紋波和阻帶紋波的比值(注意這里紋波是線性坐標(biāo)下的值,), 其中和由階方程(degree equation)聯(lián)系：

可以看到這個(gè)公式(2)和切比雪夫多項(xiàng)式

類比具有相同的形式, 實(shí)際上切比雪夫?yàn)V波器就是橢圓函數(shù)濾波器的特殊情況, 即當(dāng)阻帶紋波無(wú)窮大和阻帶頻率無(wú)窮大的特殊情況。

同樣的當(dāng)?shù)耐◣Ъy波為0, 且通帶頻率也為0, 即時(shí), 可以得到逆切比雪夫?yàn)V波器。

當(dāng)?shù)耐◣Ъy波為0, 阻帶紋波為無(wú)窮大, 且通帶頻率為0, 且阻帶頻率為無(wú)窮大, 也即時(shí), 可以得到逆切比雪夫?yàn)V波器。

下圖顯示了三種濾波器在相同階數(shù), 相同截止頻率下的動(dòng)圖展示：

橢圓到切比雪夫?yàn)V波器,從60dB變到1000dB動(dòng)圖:

橢圓到逆切比雪夫?yàn)V波器,從3dB變到0.000001dB動(dòng)圖:

橢圓到巴特沃斯濾波器,從60dB變到1000dB且從3dB變到0.000001dB動(dòng)圖:

橢圓函數(shù)濾波器的零極點(diǎn)

橢圓函數(shù)濾波器的極點(diǎn)可以通過(guò)求解傳遞函數(shù)分母為0求解:

其中，是濾波器階數(shù)。橢圓函數(shù)濾波器的零點(diǎn)和的極點(diǎn)相同, 這也比較容易理解, 設(shè), 代入頻率響應(yīng)表達(dá)式, 通分后得到其零點(diǎn)和的一樣, 在復(fù)平面中找到對(duì)應(yīng)的極點(diǎn)位置得到:

濾波器零極點(diǎn)分析工具可以對(duì)各種不同濾波器進(jìn)行零極點(diǎn)分析, 源碼見(jiàn), 下圖是橢圓函數(shù)濾波器參數(shù)變化時(shí)濾波器的零極點(diǎn)和頻響變化仿真結(jié)果:

下圖是不同類型濾波器零極點(diǎn)對(duì)比圖:

可見(jiàn)橢圓函數(shù)濾波器的極點(diǎn)形狀不是一個(gè)橢圓, 所以之前的猜想并不正確, 橢圓函數(shù)濾波器得名并不是由于其極點(diǎn)形狀是一個(gè)橢圓, 而是于其衰減特性可以由雅可比橢圓函數(shù)來(lái)描述。濾波器名稱由特征函數(shù)決定, 可以看到以后討論的貝塞爾濾波器、高斯濾波器、勒讓德濾波器均是依據(jù)特征函數(shù)或傳遞函數(shù)來(lái)得名的。

下圖是濾波器在-平面的曲線圖，從中可以看到4階橢圓函數(shù)濾波器的零極點(diǎn)的分布情況(注意這個(gè)并不能直接通過(guò)電路綜合，需要莫比烏斯變換后才能綜合，見(jiàn)下文)：

由衰減確定濾波器階數(shù)

橢圓函數(shù)濾波器的參數(shù)由如下圖所示(,)：

橢圓函數(shù)濾波器階數(shù)和濾波器的衰減參數(shù)由階方程決定, 階方程將橢圓函數(shù)濾波器的階數(shù), 紋波以及通帶阻帶頻率聯(lián)系起來(lái)了, 其底層原理是由于在解橢圓積分方程中需要讓虛部周期相同,實(shí)部周期呈整數(shù)倍,這樣才能滿足通帶和阻帶等紋波。階方程的數(shù)值計(jì)算用到了蘭登變換(Landen Transform，鏈接提供了多種計(jì)算方法), 現(xiàn)在有現(xiàn)成的迭代方法可以求解, Matlab中可以使用ellipdeg函數(shù)求解。

上式是已知通帶阻帶頻率, 通帶阻帶紋波求階數(shù), 其值可能是個(gè)小數(shù), 需要向上取整, 取整后還需要再使用一次階方程求解出對(duì)應(yīng)的值。

這里也給出階方程的直接求解方法(也是雅可比所提出的模方程, Matlab中也使用了此算法)：

這里給出一個(gè)匹配型3階橢圓函數(shù)濾波器設(shè)計(jì)實(shí)例：

非匹配型橢圓函數(shù)濾波器設(shè)計(jì)

非匹配型橢圓函數(shù)濾波器傳遞函數(shù)為:

其傳遞函數(shù)的零極點(diǎn)和匹配型濾波器的相同，不同的在于特征函數(shù)零點(diǎn)，依據(jù)特征函數(shù)定義，并令其等于0可以求得零點(diǎn)為:

其中為虛數(shù),為雅可比橢圓反正弦函數(shù),為通帶紋波,為濾波器階數(shù),為增益因子, 和之前有效衰減一樣:

通過(guò)以上求解可以得到計(jì)算出濾波器輸入阻抗，同逆切比雪夫?yàn)V波器設(shè)計(jì)一樣，由零點(diǎn)移除技術(shù)可以綜合出所需要的電路網(wǎng)絡(luò)。

這里給出一個(gè)非匹配型3階橢圓函數(shù)濾波器設(shè)計(jì)實(shí)例：

偶數(shù)階橢圓函數(shù)濾波器設(shè)計(jì)

偶數(shù)階橢圓函數(shù)濾波器為什么和奇數(shù)階濾波器設(shè)計(jì)獨(dú)立開(kāi)來(lái)的原因與之前Chebyshev濾波器設(shè)計(jì)一樣, 也是由于0頻率處衰減有個(gè)和紋波一樣的衰減, 以及無(wú)窮頻率處衰減不為0所導(dǎo)致的, 解決此問(wèn)題的辦法就是使用頻率變換, 即雙線性變換(也稱為莫比烏斯變換, M?bius transformation):

基本原則是:

當(dāng)時(shí),,是偶數(shù)階橢圓函數(shù)濾波器的第一個(gè)紋波最高點(diǎn)位置所對(duì)應(yīng)的頻率。

當(dāng)時(shí),。

當(dāng)時(shí),,是偶數(shù)階橢圓函數(shù)濾波器的最后一個(gè)零點(diǎn)位置所對(duì)應(yīng)的頻率。

通過(guò)以上條件得到:

其中：

這里給出一個(gè)匹配型4階橢圓函數(shù)濾波器設(shè)計(jì)實(shí)例：

不同濾波器通帶類型之間的轉(zhuǎn)換

只要有了低通原型, 其他濾波器通帶類型之間的轉(zhuǎn)換同Butterworth濾波器。

橢圓函數(shù)濾波器設(shè)計(jì)軟件

基于Matlab的appdesign工具開(kāi)發(fā)了一套橢圓函數(shù)濾波器設(shè)計(jì)軟件,主要特點(diǎn)是：

支持橢圓函數(shù)濾波器(Elliptic Function Filter)、支持切比雪夫?yàn)V波器(Chebyshev I)、逆切比雪夫?yàn)V波器(Chebyshev II, Inverse Chebyshev)、巴特沃斯濾波器(Butterworth)設(shè)計(jì)

支持4種不同濾波器通帶類型(LPF,HPF,BPF,BRF)設(shè)計(jì)

T型和PI型結(jié)構(gòu)濾波器隨意切換

可以設(shè)置阻帶衰減決定濾波器階數(shù)

可以設(shè)置通帶衰減來(lái)綜合濾波器

可以隨意配置負(fù)載和終端阻抗, 并支持一端接載(源端電阻短路, 源端電流源, 終端開(kāi)路, 終端短路)設(shè)計(jì)

可以幅頻響應(yīng)分析、零極點(diǎn)分析、瞬態(tài)分析

可以顯示理想頻率響應(yīng)、零極點(diǎn)和實(shí)際仿真的的頻率響應(yīng)、零極點(diǎn)

可以支持實(shí)際標(biāo)準(zhǔn)器件逼近設(shè)計(jì)

Elliptic LPF設(shè)計(jì)舉例

設(shè)計(jì)一款-3dB截止頻率為1GHz,阻帶紋波為60dB，通帶紋波1dB的7階低通Elliptic濾波器，輸入輸出阻抗為50歐姆，設(shè)計(jì)過(guò)程如下：

最終設(shè)計(jì)參數(shù)如下：

零極點(diǎn)仿真結(jié)果：

Elliptic BPF設(shè)計(jì)舉例

設(shè)計(jì)6階帶通Elliptic濾波器, 通帶紋波為0.1dB, 阻帶衰減為70dB, 中心頻率為1GHz，帶寬為500MHz，50歐姆輸入，輸出阻抗為高阻，最后進(jìn)行瞬態(tài)仿真，設(shè)計(jì)過(guò)程如下：

最終設(shè)計(jì)參數(shù)如下：

瞬態(tài)仿真結(jié)果：

程序的Matalb源碼已經(jīng)上傳GitHub中(https://github.com/etools361/MatlabEllipticFilterDesignApp)，有興趣的同學(xué)可以下載試用體驗(yàn)，當(dāng)然也歡迎技術(shù)交流。

審核編輯：湯梓紅