陳紹榮,劉郁林,朱行濤,徐 舜
(1.陸軍工程大學(xué)通信士官學(xué)校,重慶 400035;2.重慶市經(jīng)信委,重慶 400015)
在國(guó)內(nèi)外《數(shù)字信號(hào)處理》著作[1-3]中,介紹了歸一化巴特沃斯模擬低通濾波器的設(shè)計(jì)方法及歸一化切比雪夫Ⅰ型模擬低通濾波器的設(shè)計(jì)方法。巴特沃斯模擬低通濾波器的幅頻特性,無(wú)論在通帶和阻帶都是隨頻率單調(diào)遞減,若在通帶邊緣滿足指標(biāo)要求,則在通帶內(nèi)肯定會(huì)有富裕量,也就是會(huì)超過(guò)指標(biāo)的要求,因而并不經(jīng)濟(jì),故更有效的辦法是將指標(biāo)的精度要求均勻地分布在通帶內(nèi),或均勻地分布在阻帶內(nèi),或同時(shí)均勻地分布在通帶、阻帶內(nèi)。這時(shí)就可以設(shè)計(jì)出階數(shù)較低的濾波器。這種精度均勻地分布的辦法可通過(guò)選擇具有等波紋特性的逼近函數(shù)來(lái)完成。切比雪夫模擬低通濾波器的幅度特性就是在一個(gè)頻帶(通帶或阻帶)范圍內(nèi)具有這種等波紋特性。一種是通帶范圍內(nèi)為等波紋,阻帶范圍內(nèi)為單調(diào)遞減函數(shù),稱(chēng)為切比雪夫Ⅰ型,一種是阻帶范圍內(nèi)為等波紋,通帶范圍內(nèi)為單調(diào)遞減函數(shù),稱(chēng)為切比雪夫Ⅱ型。在應(yīng)用中,應(yīng)根據(jù)要求來(lái)確定采用哪種類(lèi)型的切比雪夫模擬低通濾波器。在《數(shù)字信號(hào)處理》著作[4]中,雖然介紹了歸一化切比雪夫Ⅱ型模擬低通濾波器的設(shè)計(jì)方法,但是推導(dǎo)過(guò)略,不便于理解,而且給出的確定歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)公式有誤,本文在《數(shù)字信號(hào)處理》著作[5]基礎(chǔ)上,詳細(xì)介紹歸一化切比雪夫Ⅱ型模擬低通濾波器的設(shè)計(jì)方法,基于歸一化切比雪夫Ⅱ型模擬低通濾波器來(lái)設(shè)計(jì)切比雪夫Ⅱ型模擬高帶通濾波器的原理、步驟及方法。
歸一化切比雪夫多項(xiàng)式定義為:
(1)當(dāng)|λ|≤1 時(shí),由式(1)可知:
①當(dāng)n=0 時(shí),由式(2)可得:
②當(dāng)n=1 時(shí),由式(2)可得:
③當(dāng)n ≥1 時(shí),在式(2)中,令λ=cosθ,可得:
由于:
考慮到式(5),由式(6)可得:
由式(7)可得遞推公式:
(2)當(dāng)|λ|≥1 時(shí),由式(1)可知:
①當(dāng)n=0 時(shí),由式(9)可得:
②當(dāng)n=1 時(shí),由式(9)可得:
③當(dāng)n ≥1 時(shí),在式(9)中,令λ=coshθ,可得:
由于:
考慮到式(12),由式(13)可得:
由式(14)可得與式(8)相同的遞推公式,即:
當(dāng)n ≥1 時(shí),由遞推式(8)或遞推式(15),并考慮到式(3)及式(4)或(10)及式(11),可以得到:
現(xiàn)將計(jì)算的結(jié)果制成表格,如表1 所示。
表1 歸一化切比雪夫多項(xiàng)式表
由表1 可知,歸一化切比雪夫多項(xiàng)式Cn(λ)是λ的n 次多項(xiàng)式,并且,首項(xiàng)為2n-1λn。當(dāng)|λ|≤1 時(shí),Cn(λ)在-1 至-1 之間波動(dòng)變化,當(dāng)|λ|>1 時(shí),Cn(λ)按雙曲余弦單調(diào)遞增。
切比雪夫Ⅱ型模擬低通濾波器的幅頻特性平方函數(shù)定義為:
式中,Ωp為通帶上截止角頻率,Ωst為阻帶下截止角頻率;n 為切比雪夫Ⅱ型模擬低通濾波器的階數(shù);ε 為常數(shù),并且0<ε<1。
當(dāng)Ω=Ωp時(shí),由式(20)可得:
令
考慮到式(23),則式(20)可寫(xiě)成:
由式(24)可得:
若對(duì)切比雪夫Ⅱ型模擬低通濾波器的角頻率作歸一化處理,由式(20)可知,則歸一化切比雪夫Ⅱ型模擬低通濾波器的幅頻特性平方函數(shù)可寫(xiě)成:
式中,λ=Ω/Ωp為歸一化角頻率,λs=Ωst/Ωp。
考慮到式(26),則有:
由式(27)可知,HL(p)HL(-p)除了有2n個(gè)極點(diǎn)外,還有2n 個(gè)零點(diǎn)。
令:
即:
亦即:
由式(30)可得HL(p)HL(-p)的零點(diǎn),即:
式中,r=0,1,2,…,2n-1。
式(31)表明,HL(p)HL(-p)的零點(diǎn)分布在在p平面的虛軸上,若n 為奇數(shù),當(dāng)2r+1=n 時(shí),則有pr=j∞,即HL(p)HL(-p)在p 平面虛軸上無(wú)窮遠(yuǎn)處有一個(gè)零點(diǎn)。
令:
由式(32)可得:
定義φ=arccos(jλs/p),則有:
顯然,φ 應(yīng)是復(fù)數(shù)。為此,令φ=φ1+jφ2,則有:
考慮到式(36),則有:
由式(37)可知cos(nφ1)=0,于是:
將式(39)代入式(38),則有:
由式(39)及式(41)分別可以得到φ1及φ2,將φ1及φ2代入式(35),得到HL(p)HL(-p)的極點(diǎn),即:
式中,l=1,2,3,…,2n。
如果令pl=σl+jλl,考慮到式(42),則有:
顯然cosh φ2>sinh φ2,因此橢圓的焦點(diǎn)在p 平面的虛軸上。
式(43)表明,歸一化切比雪夫Ⅱ型模擬低通濾波器極點(diǎn)pl的實(shí)部和虛部滿足橢圓方程,即極點(diǎn)pl落在橢圓上。
由式(42)求出2n 個(gè)極點(diǎn)pl,一半屬于HL(p),一半屬于HL(-p),為保證所設(shè)計(jì)的歸一化切比雪夫Ⅱ型模擬低通濾波器是穩(wěn)定系統(tǒng),應(yīng)將p 平面左半面的極點(diǎn)賦于HL(p)。
若規(guī)定φ2>0,l=1,2,…,n,則由式(42)可以得 到HL(p)的n 個(gè) 極 點(diǎn)pl(l=1,2,…,n)。由式(27)可知,HL(p)HL(-p)的分子和分母分別為和由 于是 常 數(shù),因 此,由表1 的歸一化切比雪夫多項(xiàng)式可知,HL(p)HL(-p)的分子和分母均是jλs/p 的2n 次多項(xiàng)式,并且最高冪次項(xiàng)均為[2n-1(jλs/p)n]2。換言之,由于描述HL(p)HL(-p)的分子和分母的2n 次多項(xiàng)式的變量為jλs/p,而且最高冪次(jλs/p)2n的系數(shù)相同,因此,歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)可表示成:
式中,pr(r=0,1,2,…,n-1)及pl(l=1,2,…,n)分別為歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)的零點(diǎn)和極點(diǎn)。
(1)將角頻率作歸一化處理,得到歸一化切比雪夫Ⅱ型模擬低通濾波器幅頻特性平方函數(shù)。
(2)確定歸一化切比雪夫Ⅱ型模擬低通濾波器涉及的常數(shù)ε 及階數(shù)n。
考慮到式(26),則衰減函數(shù)α(λ)可表示成:
設(shè)Ω=Ωp,即λ=λp=Ωp/Ωp=1 時(shí),通帶允許的最大衰減為αp,考慮到式(45),則有:
設(shè)Ω=Ωst,即λ=λs時(shí),阻帶應(yīng)達(dá)到的最小衰減為αs,考慮到式(45),則有:
考慮到式(46),由式(48)可得:
由式(49)可得:
這樣,首先利用式(46)求出常數(shù)ε,再利用式(49)計(jì)算出a2值,最后利用式(50)就可求出歸一化切比雪夫Ⅱ型模擬低通濾波器的階數(shù)n。
(3)利用式(31)、式(42)及式(44),可求出歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)HL(p)。
(4)作逆歸一化處理,即令p=s/Ωp,可從歸一化的切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)HL(p)得到切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)H(s),即:
可以通過(guò)頻率變換的方法,基于歸一化切比雪夫Ⅱ型模擬低通濾波器來(lái)完成切比雪夫Ⅱ型模擬高通濾波器的設(shè)計(jì),其設(shè)計(jì)過(guò)程,如圖1 所示。
圖1 切比雪夫Ⅱ型模擬高通濾波器的設(shè)計(jì)過(guò)程
歸一化切比雪夫Ⅱ型模擬高通濾波器及模擬低通濾波器的幅頻特性,分別如圖2 及圖3 所示。
圖2 歸一化切比雪夫Ⅱ型模擬高通濾波器的幅頻特性
圖3 歸一化切比雪夫Ⅱ型模擬低通濾波器的幅頻特性
為了保證頻率變換后,兩個(gè)歸一化的幅頻特性|HG(jη)|和|HL(jλ)|,在通帶內(nèi)和阻帶內(nèi)的衰減dB 數(shù)相同,各頻率點(diǎn)的對(duì)應(yīng)關(guān)系,如表2 所示。
表2 歸一化切比雪夫Ⅱ型模擬高通與低通濾波器幅頻特性各頻率點(diǎn)的對(duì)應(yīng)關(guān)系
由表2可知,歸一化角頻率λ 與η 滿足下述關(guān)系:
考慮到式(52),則|HL(jλ)|中對(duì)應(yīng)的復(fù)變量p=jλ 與|HG(jη)|中對(duì)應(yīng)的復(fù)變量q=jη 滿足下述關(guān)系:
于是,作逆復(fù)頻率變換,可得歸一化切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù),即:
由于:
基于式(55),作逆歸一化處理,并注意到式(54),可得切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù),即:
式(56)表明,可將逆復(fù)頻率變換和逆歸一化處理兩個(gè)步驟,合為一步完成便可。
(1)將待設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的通帶下截止角頻率和阻帶上截止角頻率作歸一化處理,可得:
(2)利用式(52)作頻率變換,并考慮到頻率點(diǎn)的對(duì)應(yīng)關(guān)系,則有:
因此,λp=1 可直接給出,而αp、αs保持不變;其中,αp、αs分別為待設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器通帶允許的最大衰減及阻帶應(yīng)達(dá)到的最小衰減。
(3)按要求設(shè)計(jì)歸一化切比雪夫Ⅱ型模擬低通濾波器,得到HL(p)。
(4)利用式(54)作逆復(fù)頻率變換,得到歸一化切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù)HG(q)。
(5)利用式(56)作逆歸一化處理,得到切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù)H(s)。
例:切比雪夫Ⅱ型模擬高通濾波器的技術(shù)指標(biāo)為:通帶下截止角頻率Ωp=5×106rad/s,通帶最大衰減αp=101g(8 194/4 225)dB,阻帶上截止角頻率Ωst=4×106rad/s,阻帶最小衰減αs=12 dB。
試按下述要求,完成該模擬高通濾波器的設(shè)計(jì)。
(1)設(shè)計(jì)歸一化切比雪夫Ⅱ型模擬低通濾波器。
(2)若取Ωp=4×106rad/s,對(duì)(1)的結(jié)果作逆歸一化處理,即令p=s/Ωp,寫(xiě)出相應(yīng)的切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù),并畫(huà)出實(shí)現(xiàn)電路。
(3)基于(1)的結(jié)果,寫(xiě)出所設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù),并畫(huà)出實(shí)現(xiàn)電路。
解:(1)設(shè)計(jì)歸一化切比雪夫Ⅱ型模擬低通濾波器
①將待設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的通帶下截止角頻率和阻帶上截止角頻率作歸一化處理
因?yàn)棣竝=5×106rad/s,Ωst=4×106rad/s,所以ηp=1,ηs=Ωst/Ωp=0.8。
②作頻率變換,考慮到式(57)及式(58),則有:λp=1/ηp=1,λs=1/ηs=5/4
并且,αp=101g(8 194/4 225)dB;αs=12 dB。
③確定歸一化切比雪夫Ⅱ型模擬低通濾波器涉及的常數(shù)ε 及階數(shù)n
由式(46)可得:
由式(61)可得:
將αs=12 dB 及式(61)代入式(49),可得:
由式(63)可得:
設(shè):
由式(65)可得:
由式(66)可得一元二次方程,即:
解得:
考慮到式(68),則有:
將式(64)代入式(69)可得:
同理,當(dāng)λs=5/4 時(shí),由式(69)可得:
將式(70)及式(71)代入式(50),可得:
由式(72)可知,?。?/p>
④確定歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)HL(p)
將λs=5/4,式(62)及式(73)代入到式(41),可得:
考慮到:
將式(74)代入式(75)可得:
考慮到式(69),由式(76)可得:
考慮到式(77),則?。?/p>
考慮到式(78),則有:
將λs=5/4,式(73)、式(79)、式(80)及式(81)代入式(42),則有:
由式(82)可得極點(diǎn):
將λs=5/4,式(73)代入式(31),則有:
由式(86)可得零點(diǎn):
將式(83)、式(84)、式(85)、式(87)、式(88)及式(89)代入式(44),可得歸一化切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù)HL(p),即:
(2)考慮到式(90),由式(51)可得切比雪夫Ⅱ型模擬低通濾波器的轉(zhuǎn)移函數(shù),即:
考慮到式(91),則描述切比雪夫Ⅱ型模擬低通濾波器的信號(hào)流圖,如圖4 所示。
圖4 切比雪夫Ⅱ型模擬低通濾波器的信號(hào)流圖
依據(jù)圖4 所示的信號(hào)流圖,則利用電阻、電容及運(yùn)算放大器來(lái)實(shí)現(xiàn)切比雪夫Ⅱ型模擬低通濾波器的電路圖,如圖5 所示,其中,R1=100 kΩ,R2=R4=200 kΩ,R3=80 kΩ,R5=R6=R7=100 kΩ,R8=50 kΩ,Rp2=2.1 kΩ,Rp3=2 kΩ,Rp4=1.5 kΩ。
圖5 實(shí)現(xiàn)切比雪夫Ⅱ型模擬低通濾波器的電路圖
(3)考慮到式(56),可得所設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的轉(zhuǎn)移函數(shù),即:
下面檢驗(yàn)所設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的通帶最大衰減和阻帶最小衰減是否滿足技術(shù)指標(biāo)要求。
考慮到式(92),則切比雪夫Ⅱ型模擬高通濾波器的頻率特性為:
由式(93)可得:
由式(94)可得:
由式(93)可得:
由式(96)可得:
因此,所設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的通帶最大衰減和阻帶最小衰減滿足技術(shù)指標(biāo)要求。
其實(shí),式(92)可寫(xiě)成:
考慮到式(98),則所設(shè)計(jì)的切比雪夫Ⅱ型模擬高通濾波器的信號(hào)流圖,如圖6 所示。
依據(jù)圖6 所示的信號(hào)流圖,則利用電阻、電容及運(yùn)算放大器來(lái)實(shí)現(xiàn)切比雪夫Ⅱ型模擬高通濾波器的電路圖,如圖7 所示,其中,R1=30 kΩ,R2=150 kΩ,R3=100 kΩ,R4=50 kΩ,R5=300 kΩ,R6=R8=50 kΩ,R7=100 kΩ,Rp2=5 kΩ,Rp3=Rp4=10 kΩ。
圖6 切比雪夫Ⅱ型模擬高通濾波器的信號(hào)流圖
圖7 實(shí)現(xiàn)切比雪夫Ⅱ型模擬高通濾波器的電路圖
本文介紹了歸一化切比雪夫多項(xiàng)式及其遞推公式,歸一化切比雪夫Ⅱ型模擬低通濾波器的極點(diǎn)分布的特點(diǎn),確定歸一化切比雪夫Ⅱ型模擬低通濾波器的零極點(diǎn)和轉(zhuǎn)移函數(shù)的公式,設(shè)計(jì)歸一化切比雪夫Ⅱ型模擬低通濾波器的步驟;再介紹了利用頻率變換來(lái)設(shè)計(jì)切比雪夫Ⅱ型模擬高通濾波器的方法;最后給出了基于歸一化切比雪夫Ⅱ型模擬低通濾波器來(lái)設(shè)計(jì)切比雪夫Ⅱ型模擬高通濾波器的實(shí)例。