梁貴書,周安東

(華北電力大學(xué) 電氣與電子工程學(xué)院,河北 保定 071003)

0 引 言

電氣設(shè)備與器件的分?jǐn)?shù)階電路模型構(gòu)建離不開電路綜合理論。由于自然界中的材料固有的分?jǐn)?shù)階本質(zhì)[1],電氣設(shè)備與器件實(shí)際上也具有分?jǐn)?shù)階的特性。在電氣與電子工程領(lǐng)域中,分?jǐn)?shù)階理論已經(jīng)廣泛應(yīng)用于單相逆變器的分?jǐn)?shù)階建模[2]、分?jǐn)?shù)階控制[3-4]、分?jǐn)?shù)階混沌系統(tǒng)的研究[5]。分?jǐn)?shù)階電路綜合即利用分?jǐn)?shù)階元件與傳統(tǒng)整數(shù)階元件構(gòu)建電網(wǎng)絡(luò)來實(shí)現(xiàn)分?jǐn)?shù)階策動(dòng)點(diǎn)阻抗函數(shù)。其綜合方法通常存在兩種方式,一種是通過整數(shù)階電路以逼近擬合的方式來實(shí)現(xiàn)分?jǐn)?shù)階電路,另一種是直接通過分?jǐn)?shù)階元件構(gòu)建電路模型來實(shí)現(xiàn)。相比之下,后者得到的電路模型更為簡潔,因此可有效降低元件的冗余度,且精確度更高更易于仿真分析。隨著材料科學(xué)的發(fā)展與制造水平的提升,已經(jīng)可制造出分?jǐn)?shù)階元件[6],因此可直接使用分?jǐn)?shù)階元件來實(shí)現(xiàn)分?jǐn)?shù)階電路。電路實(shí)現(xiàn)的關(guān)鍵環(huán)節(jié)就在于選擇合適的電路綜合方法,考慮到無源網(wǎng)絡(luò)在仿真過程中的穩(wěn)定性要優(yōu)于有源網(wǎng)絡(luò),因此研究探索分?jǐn)?shù)階電網(wǎng)絡(luò)的無源綜合理論方法具有重要意義。

由于分?jǐn)?shù)階電路的策動(dòng)點(diǎn)阻抗函數(shù)復(fù)雜度較高,國內(nèi)外目前難以有一個(gè)通用的實(shí)現(xiàn)方法。文獻(xiàn)[7-8]研究了由整數(shù)階二端口RLC阻抗網(wǎng)絡(luò)端接分?jǐn)?shù)階電容與分?jǐn)?shù)階電感元件來實(shí)現(xiàn)特殊的分?jǐn)?shù)階電路,但該方法對超過兩個(gè)分?jǐn)?shù)階電抗元件的情形并不適用。基于雙變量達(dá)林頓綜合方法,文獻(xiàn)[9]解決了分?jǐn)?shù)階雙元次阻抗網(wǎng)絡(luò)的綜合問題,然而該綜合方法的缺點(diǎn)是必須借助多口變壓器,不利于分?jǐn)?shù)階電網(wǎng)絡(luò)的建模與分析。文獻(xiàn)[10]提出了一個(gè)三端口電阻網(wǎng)絡(luò)端接兩個(gè)分?jǐn)?shù)階電抗元件的的最少儲(chǔ)能元件綜合方法,但該方法可適用的范圍太小。

目前,尚未見分?jǐn)?shù)階三種元件電路綜合的相關(guān)報(bào)道。文章對分?jǐn)?shù)階RLαCβ三種元件電路的綜合問題進(jìn)行了探索。首先,基于特勒根能量函數(shù)的多變量域表達(dá)形式,推導(dǎo)出RLαCβ三種元件電路的Foster綜合方法以及Cauer綜合方法,其中,Foster綜合方法是基于部分分式的展開的思想來實(shí)現(xiàn)電路,Cauer綜合方法基于連分式的展開的思想來實(shí)現(xiàn)電路。進(jìn)一步,給出了分?jǐn)?shù)階RLαCβ三種元件電路阻抗函數(shù)的一般表達(dá)形式,最后,以一個(gè)具體的實(shí)例計(jì)算與仿真,驗(yàn)證了文中電路綜合方法。與以往綜合方法相比較,所提出的綜合方法不必使用多口變壓器來實(shí)現(xiàn)電路,計(jì)算方法且較為簡單實(shí)用,因此更有利于分?jǐn)?shù)階電路系統(tǒng)的建模與分析。

1 分?jǐn)?shù)階電抗元件

分?jǐn)?shù)階電抗元件主要是指分?jǐn)?shù)階電感與分?jǐn)?shù)階電容。分?jǐn)?shù)階電感的時(shí)域特性方程為[11]:

(1)

分?jǐn)?shù)階電感的電感值與階次一般表示為(Lα,α),電感值單位為H/s1-α,其中,0≤α≤1。

分?jǐn)?shù)階電容的時(shí)域特性方程為[11]:

(2)

分?jǐn)?shù)階電容的電容值與階次一般表示為(Cβ,β),電容值單位為F/s1-β。其中,0≤β≤1。對式(1)與式(2)取拉式變換,得到分?jǐn)?shù)階電感與分?jǐn)?shù)階電容的頻域特性方程為:

U(s)=LαsαI(s)

(3)

(4)

圖1給出了分?jǐn)?shù)階電感與分?jǐn)?shù)階電容的電路符號。其中u(t)表示分?jǐn)?shù)階電感或分?jǐn)?shù)階電容的端電壓,i(t)表示流過分?jǐn)?shù)階電感或分?jǐn)?shù)階電容的電流。且u(t)與i(t)為關(guān)聯(lián)參考方向。

圖1 分?jǐn)?shù)階電抗元件的電路符號

2 分?jǐn)?shù)階RLαCβ電路多變量域綜合

2.1 RLαCβ電路阻抗函數(shù)的多變量域Foster電路實(shí)現(xiàn)

包含RLαCβ三種元件的一般支路如圖2所示。基于特勒根定理及其能量函數(shù),得到圖2所示的單端口策動(dòng)點(diǎn)阻抗函數(shù)為:

圖2 分?jǐn)?shù)階RLαCβ的一般支路

(5)

基于變量代換sα=p1,sβ=p2,式(5)變?yōu)?

(6)

進(jìn)一步對式(6)整理得到:

(7)

根據(jù)式(6)與式(7)可知Rk=Rk1+Rk2,對于式(7),由于Rk,Lk,Ck均為實(shí)常數(shù),令:

(8)

進(jìn)一步,令:

(9)

將式(8)與式(9)帶入式(7),從而可使得式(7)可簡化整理為:

(10)

顯然,在式(10)中,Z(q)是一個(gè)關(guān)于變量q的奇函數(shù),其表達(dá)式為:

(11)

根據(jù)文獻(xiàn)[12]的結(jié)論,具有式(11)形式的Z(q)為q域的電抗函數(shù),其可寫為式(12)所示的部分分式展開的形式為:

(12)

Z(q)可實(shí)現(xiàn)為梯形電抗網(wǎng)絡(luò),在式(12)中:k∞,k0是q平面上Z(q)在無窮遠(yuǎn)與原點(diǎn)處的留數(shù);ki表示q平面Z(q)在虛軸上極點(diǎn)留數(shù)之和,k∞,k0,ki均為非負(fù)實(shí)數(shù)。根據(jù)式(8)～式(12),得到RLαCβ三種元件電路分?jǐn)?shù)階阻抗函數(shù)多變量域的展開式為:

(13)

在式(13)中,k∞,k0,ki為:

(14)

顯然,根據(jù)式(13)可知,第一項(xiàng)可以實(shí)現(xiàn)為電感串聯(lián)電阻的形式,第二項(xiàng)可以實(shí)現(xiàn)為電容串聯(lián)電阻的形式,剩余的其他項(xiàng)可以實(shí)現(xiàn)為電感串聯(lián)電阻之后整體并聯(lián)電容串聯(lián)電阻的形式,之后將每一項(xiàng)所實(shí)現(xiàn)的部分依次串聯(lián)起來,最終整體阻抗函數(shù)所實(shí)現(xiàn)的Foster形式的電路如圖3所示。

圖3 RLαCβ阻抗函數(shù)的Foster電路

此外,可以完全對偶地推出RLαCβ電路的導(dǎo)納函數(shù)多變量域部分分式展開式為:

(15)

2.2 RLαCβ電路阻抗函數(shù)多變量域Cauer梯形實(shí)現(xiàn)

RLαCβ電路多變量域阻抗函數(shù)為Z(p1,p2),假定其在p1→∞處有一極點(diǎn),則可以得到:

Z(p1,p2)=k∞1(p1+a)+Z2(p1,p2)=k∞1(p1+a)+

(16)

移去Z(p1,p2)在p1→∞處的極點(diǎn),根據(jù)式(16)可知,k∞1(p1+a)可以實(shí)現(xiàn)為電感k∞1p1串聯(lián)電阻k∞1a的串聯(lián)支路形式,且余函數(shù)Z2(p1,p2)仍然為RLαCβ電路的阻抗函數(shù),進(jìn)一步,根據(jù)式(13)可知,Z2(p1,p2)在p2=-b-1處,有:

(17)

根據(jù)式(17)可知Z2(p1,p2)在p2=-b-1處有一零點(diǎn),所以Y2(p1,p2)在p2=-b-1處有一極點(diǎn),根據(jù)式(15),得到:

(18)

(19)

(20)

重復(fù)進(jìn)行式(16)～式(20)的極點(diǎn)提取過程,每提取一次p1→∞處的極點(diǎn)便得到一個(gè)電感串電阻的串聯(lián)支路,每提取一次p2=-b-1處的極點(diǎn)一個(gè)電容串電阻的并聯(lián)支路,最終阻抗函數(shù)Z(p1,p2)的Cauer梯形電路的實(shí)現(xiàn)如圖4所示。

圖4 RLαCβ阻抗函數(shù)的Cauer電路

需要說明的是,以上推導(dǎo)過程首先假定在p1→∞處為阻抗函數(shù)Z(p1,p2)的極點(diǎn),若在p1→∞時(shí),RLαCβ阻抗函數(shù)Z(p1,p2)→0,此時(shí)首先對Y(p1,p2)進(jìn)行極點(diǎn)的移出,其電路實(shí)現(xiàn)步驟與式(16)～式(20)的過程是完全相同的。

根據(jù)式(16)～式(20),我們可以得到圖4所對應(yīng)的阻抗函數(shù)Z(p1,p2)的連分展開式為:

Z(p1,p2)=k∞1(p1+a)+

(21)

2.3 分?jǐn)?shù)階RLαCβ電路阻抗函數(shù)s域一般表達(dá)式

根據(jù)式(12),Z(q)為q域的電抗函數(shù),所以Z(q)為奇函數(shù),分兩種情況討論并最終得到分?jǐn)?shù)階RLαCβ阻抗函數(shù)s域的一般表達(dá)形式:

(1)當(dāng)Z(q)在原點(diǎn)處有零點(diǎn)時(shí),由于q域網(wǎng)絡(luò)僅由電感和電容組成,所以,Z(q)可表示為[12]:

(22)

式中ci,di(i=0,1,2,…),k均為正實(shí)數(shù),且要求有[12]di

Z(s)=(sα+a)·

(23)

式中A0,…,AM,B0,…,BN>0。

(2)當(dāng)Z(q)在原點(diǎn)處有極點(diǎn)時(shí),由于q域網(wǎng)絡(luò)僅由電感和電容組成,所以,Z(q)可表示為[12]:

(24)

式中c′i,d′i(i=0,1,2,…),k均為正實(shí)數(shù),且要求有[12]c′i

(25)

同理可知,A0,…,AM,B0,…,BN>0。

綜上所述,分?jǐn)?shù)階RLαCβ電路的阻抗函數(shù)s域表達(dá)形式必然為式(23)或式(25)的形式。

3 算例驗(yàn)證與仿真

如式(24)所示的分?jǐn)?shù)階s域的阻抗函數(shù):

(26)

對式(26)重新整理,得到:

(27)

將式(27)與式(25)對比可知,二者形式一致。因此式(27)一定可以實(shí)現(xiàn)為圖3與圖4所示的RLαCβ電路。對式(26)變量代換s0.1=p1,s0.2=p2,進(jìn)而得到式(26)的多變量域表達(dá)形式為:

(28)

(1)Foster電路實(shí)現(xiàn)

根據(jù)式(13)與式(14)得到式(28)的部分分式展開式為:

(29)

根據(jù)式(29),得到其Foster電路如圖5所示。

圖5 阻抗函數(shù)的Foster電路

(2)Cauer電路實(shí)現(xiàn)

根據(jù)式(16)～式(20)的過程,移去式(28)在p1→∞處與p2=-1處的極點(diǎn),得到式(28)的連分式展開:

(30)

根據(jù)式(30),得到其Cauer電路如圖6所示。

圖6 阻抗函數(shù)的Cauer電路

為驗(yàn)證圖5與圖6所得到的分?jǐn)?shù)階電路的正確性,我們進(jìn)行頻域仿真驗(yàn)證,將圖7所示的正弦穩(wěn)態(tài)電壓激勵(lì)與階躍暫態(tài)電壓激勵(lì)應(yīng)用于式(26)的分?jǐn)?shù)階阻抗函數(shù)與圖5與圖6的分?jǐn)?shù)階阻抗網(wǎng)絡(luò),進(jìn)行數(shù)學(xué)計(jì)算與電路仿真兩方面的互相印證,從而得到圖8所示的端口電流響應(yīng)。其中,數(shù)學(xué)計(jì)算是基于端口阻抗頻域表達(dá)式U(s)=Z(s)·I(s),計(jì)算穩(wěn)態(tài)與暫態(tài)電壓激勵(lì)下的電流響應(yīng),之后通過快速傅里葉變換得到端口電流時(shí)域響應(yīng)曲線;電路計(jì)算是通過借助分?jǐn)?shù)階電抗元件的分?jǐn)?shù)階微積分定義以及L1插值法得到分?jǐn)?shù)階電抗元件的離散化模型,進(jìn)一步,通過改進(jìn)節(jié)點(diǎn)法編寫程序軟件對得到的分?jǐn)?shù)階電路進(jìn)行離散仿真。由圖8可知,其電路仿真曲線與數(shù)學(xué)計(jì)算曲線是相吻合的,因此圖5與圖6所得到的分?jǐn)?shù)階電路實(shí)現(xiàn)是正確的。

圖7 電壓激勵(lì)

圖8 電流響應(yīng)

與傳統(tǒng)綜合方法相比,文中所提出的分?jǐn)?shù)階電路綜合方法有以下顯著的特點(diǎn)與優(yōu)勢:

(1)與傳統(tǒng)的雙變量達(dá)林頓電路綜合法相比,文中所提出的電路綜合方法大大降低了計(jì)算量與計(jì)算難度。傳統(tǒng)雙變量達(dá)林頓綜合法必須借助雙變量矩陣的譜分解理論來實(shí)現(xiàn)電路,計(jì)算難度較高,因此不利于計(jì)算機(jī)的編程分析;

(2)與傳統(tǒng)的基于阻抗換標(biāo)思想的電路綜合方法相比,文中所提出的綜合方法適用于分?jǐn)?shù)階RLαCβ三種元件電路,因此具有更廣的適用范圍;

(3)在電路結(jié)構(gòu)的實(shí)現(xiàn)方面,傳統(tǒng)的分?jǐn)?shù)階電路綜合方法必須使用多口變壓器來實(shí)現(xiàn)電路,得到的電路模型復(fù)雜度較高,因此不利于對電氣設(shè)備的分?jǐn)?shù)階電路模型進(jìn)行仿真分析。文中所提出的綜合方法克服了這一缺點(diǎn)。

4 結(jié)束語

基于特勒根定理與RLαCβ單端口網(wǎng)絡(luò)的能量函數(shù),推導(dǎo)出了分?jǐn)?shù)階RLαCβ阻抗函數(shù)的部分分式展開形式的Foster電路綜合方法,以及連分式展開形式的Cauer電路綜合方法。通過具體算例的計(jì)算與仿真,驗(yàn)證了所提出的電路綜合方法。文中提出的綜合方法計(jì)算方法簡單且實(shí)用性更強(qiáng),更有利于分?jǐn)?shù)階電路系統(tǒng)的建模與分析,同時(shí)進(jìn)一步完善了分?jǐn)?shù)階電路綜合的理論體系。