一、引 言

　　隨著(zhù)計算機技術(shù)的進(jìn)步，時(shí)域有限差分法（FiniteDifference TimeDomain）是一種求解電磁問(wèn)題的數值計算技術(shù)，由K.S.Yee于1966年提出。他的基本思想是根據時(shí)域麥克斯韋方程的場(chǎng)分量微分式，用差分替代微分式，進(jìn)行各場(chǎng)分量的迭代，但是這種方法隨著(zhù)頻率升高，計算網(wǎng)格將顯著(zhù)增加，PC機的性能將很難滿(mǎn)足需要，而且單純依靠計算機性能的提高也是不實(shí)際的。例如，在分析波導膜片濾波器時(shí)，為正確模擬全部膜片的幾何結構，FDTD柵網(wǎng)的網(wǎng)格尺寸選得非常小，從而導致描述整個(gè)波導濾波器的網(wǎng)格數量非常大。由于每?jì)蓚�(gè)膜片之間都是均勻波導傳輸線(xiàn)，使用與膜片相同的柵網(wǎng)顯然是不必要的。人們曾使用非均勻FDTD柵網(wǎng)的辦法解決這個(gè)問(wèn)題，當柵網(wǎng)的大小相差比較大時(shí)，不但收斂性不易控制，而且仍無(wú)法確保節省計算時(shí)間。將Diakoptics思想運用于微波電路的全波分析，通過(guò)將電路分割為若干獨立的部分，根據每部分的具體結構采用不同的網(wǎng)格，獨立地對各個(gè)部分進(jìn)行全波時(shí)域分析，由于每部分的網(wǎng)格是均勻的，因而容易保證算法的收斂性。

　　二、Diakoptics的概念

　　Diakoptics的概念來(lái)源于網(wǎng)絡(luò )理論。其定義為：將一個(gè)網(wǎng)絡(luò )分解為若干子網(wǎng)絡(luò )，對每個(gè)子網(wǎng)絡(luò )的沖擊響應單獨求解，最后通過(guò)一定的連接條件，由諸子網(wǎng)絡(luò )的沖擊響應求出網(wǎng)絡(luò )總的響應。連接條件按形式不同可分為串行連接及并行連接。串行連接是依照一定的順序由網(wǎng)絡(luò )的一端向另一端單向連接，見(jiàn)圖1,其優(yōu)點(diǎn)是簡(jiǎn)單，但最大的問(wèn)題是當其中一個(gè)子網(wǎng)絡(luò )的沖擊響應改變時(shí)，將對其后的網(wǎng)絡(luò )產(chǎn)生影響。并行連接可克服這個(gè)缺點(diǎn)。并行連接可在任意兩個(gè)相鄰的子網(wǎng)絡(luò )間進(jìn)行，且若干并行連接可同時(shí)獨立進(jìn)行，并行時(shí)域Diakoptics假設子網(wǎng)絡(luò )為M+N端口網(wǎng)絡(luò )，其中M個(gè)端口和前一級子網(wǎng)絡(luò )相連，N個(gè)端口和后一級子網(wǎng)絡(luò )相連。子網(wǎng)絡(luò )的離散格林函數為g（i,j,n′）即j（j=1,M+N）端口t=0時(shí)刻的激勵，在i（i=1,M+N）端口t=n′時(shí)刻的沖擊響應。

　　研究微波電路問(wèn)題時(shí)，若微波電路可以被等效為一個(gè)線(xiàn)性網(wǎng)絡(luò )的話(huà)，則可以設想描述微波電路特性的格林函數可對應于電路理論中的沖擊響應函數。從電磁場(chǎng)理論角度看，時(shí)域格林函數g（r,t;r0,t0）為位于r0點(diǎn)的點(diǎn)源t0時(shí)刻施加的單位沖擊信號在觀(guān)察點(diǎn)r及t時(shí)刻的場(chǎng)，且滿(mǎn)足方程

　　兩個(gè)微波子電路連接時(shí)，其連接參考面上存在著(zhù)復雜的耦合關(guān)系，這種耦合關(guān)系可以用電磁波在存在兩個(gè)不連續界面的媒質(zhì)中反射和透射現象來(lái)形象描述，如圖1所示。那么如何將Diakoptics算法應用于微波電路特性分析中呢？在介紹這一點(diǎn)之前，本文首先簡(jiǎn)要介紹Diakoptics算法的數學(xué)描述。

圖1　媒質(zhì)中反射和透射現象可以用來(lái)形象描述兩個(gè)微波子電路間的耦合關(guān)系

　　三、Diakoptics算法的數學(xué)描述

　　以?xún)蓚�(gè)二端口網(wǎng)絡(luò )的串、并行連接給出Diakoptics算法的數學(xué)描述。圖2假設兩個(gè)子電路的反射及透射波的沖擊響應函數分別為：gr1（t），gr2（t），gt1（t），gt2（t）和hr1（t），hr2（t），ht1（t），ht2（t），上標"r"表示反射波，"t" 表示傳輸波，下標1表示從輸入參考面對電路作激勵，下標2表示從輸出參考面對電路作激勵。設f為兩個(gè)子電路連接后電路的沖擊響應函數。使用串行算法，從f 網(wǎng)絡(luò )輸入參考面看入的沖擊響應為：

　　fr1（t）=gr1（t）+gt2（t）*hr1（t）*gt1（t）+gt2（t）*hr1（t）

　　*gr2（t）*hr1（t）*gt1（t）+…+gt2（t）*（hr1（t）

　　*gr2（t））n*hr1（t）*gt1（t）+…； （2）

　　使用并行算法，從f電路的輸入端口看入的沖擊響應函數fr1（t），ft2（t）以及從f電路的輸出端口看入的沖擊響應函數fr2（t），ft1（t）分別為：

　　fr1（t）=gr1（t）+gt2（t）*hr1（t）*gt1（t）+gt2（t）*hr1（t）

　　*gr2（t）*hr1（t）*gt1（t）+…+gt2（t）*（hr1（t）

　　*gr2（t））n*hr1（t）*gt1（t）+…

　　ft2（t）=gt2（t）*hr2（t）+gt2（t）*hr1（t）*gr2（t）*ht2（t）+…

　　+gr2（t）*（hr1（t）*gr2（t））n*hr2（t）+… （3）

　　fr2（t）=hr2（t）+ht1（t）*gr2（t）*ht2（t）+ht1（t）*gr2（t）

　　*hr1（t）*gt2（t）*ht2（t）+…+ht1（t）*（gr2（t）

　　*hr1（t））n*gr2（t）*ht2（t）+…

　　ft1（t）=ht1（t）*gt1（t）+ht1（t）*gr2（t）*hr1（t）*gt1（t）+…

　　+ht1（t）*（gr2（t）*hr1（t））n*gr1（t）+…

　　其中，*代表時(shí)域卷積，上下標的含義不變。

圖2　可說(shuō)明Diakoptics算法的兩個(gè)子電路連接示意圖

　　多端口子電路連接時(shí)，上述算法依然成立，只是式中各沖擊函數應換為相應的子矩陣。例如設g網(wǎng)絡(luò )為輸入端有M個(gè)、輸出端有N個(gè)端口的M+N端口網(wǎng)絡(luò )，h網(wǎng)絡(luò )為輸入端有N個(gè)、輸出端有L個(gè)端口的N+L端口網(wǎng)絡(luò )（g與h相鄰面的端口數目應相同），g網(wǎng)絡(luò )輸入參考面處的反射、傳輸子矩陣分別為：

和

　　式中下標代表參考面，i←j的意思為：i為響應所在參考面，j為激勵所在參考面；上標代表端口，m←n的意思為：n為輸入端口，m為輸出端口。同理，g網(wǎng)絡(luò )輸出參考面處的反射、傳輸子矩陣分別為：

和

　　h網(wǎng)絡(luò )相應子矩陣可用同樣方法求得。連接后網(wǎng)絡(luò )的沖擊響應函數[f]為：

　　[fr1（t）]=[gr1（t）]+[gt2（t）]*[hr1（t）]*[gt1（t）]+[gt2（t）]

　　*[hr1（t）]*[gr2（t）]*[hr1（t）]*[gt1（t）]+…

　　[ft2（t）]=[gt2（t）]*[ht2（t）]+[gt2（t）]*[hr1（t）]*[gr2（t）]*[ht2（t）]+…

　　[fr2（t）]=[hr2（t）]+[ht1（t）]*[gr2（t）]*[ht2（t）]+[ht1（t）]

　　*[gr2（t）]*[hr1（t）]*[gr2（t）]*[ht2（t）]+…

　　[ft1（t）]=[ht1（t）]*[gt1（t）]+[ht1（t）]*[gr2（t）]*[hr1（t）]*[gt1（t）]+… （4）

　　其中[fr1（t）]、[ft1（t）]、[fr2（t）]和[ft2（t）]分別為M×M、L×M、L×L和M×L階子矩陣。下面以[gt2（t）]*[ht2（t）]為例說(shuō)明如何計算矩陣卷積，并以[gt2（t）]*[ht2（t）]的第一個(gè)元素為例，說(shuō)明其物理意義：

　　g1←11←2*h1←11←2:h子網(wǎng)絡(luò )輸出參考面上第一個(gè)端口的輸入通過(guò)gh連接面第1個(gè)端口的耦合在g子網(wǎng)絡(luò )輸入參考面上端口1產(chǎn)生的輸出；g1←21←2*h2←11←2:h子網(wǎng)絡(luò )輸出參考面上第一個(gè)端口的輸入通過(guò)gh交界面第2個(gè)端口的耦合在g子網(wǎng)絡(luò )輸入參考面上端口1產(chǎn)生的輸出；g1←N1←2*hN←11←2:h子網(wǎng)絡(luò )輸出參考面上第一個(gè)端口的輸入通過(guò)gh交界面第N個(gè)端口的耦合，在g子網(wǎng)絡(luò )輸入參考面上端口1產(chǎn)生的輸出；所以[gt2（t）]*[ht2（t）]的第一個(gè)元素描述了h網(wǎng)絡(luò )輸出參考面上第一個(gè)端口上的輸入耦合到g網(wǎng)絡(luò )輸入參考面第一個(gè)端口的輸出。

　　四、Diakoptics算法在微波電路分析中的實(shí)現

　　Diakoptics源于網(wǎng)絡(luò )理論，為將其應用于微波電路的分析中，首先需要建立適于使用Diakoptics方法的微波電路的等效電路模型。

　　1.微波電路的等效時(shí)域網(wǎng)絡(luò )模型

　　建立微波電路等效時(shí)域網(wǎng)絡(luò )模型的基本方法是：利用基函數技術(shù)（或稱(chēng)時(shí)域模技術(shù)）將參考面處的場(chǎng)表示為選定的正交基函數的線(xiàn)性組合，將一個(gè)微波網(wǎng)絡(luò )等效為一個(gè)多模電路，進(jìn)而再將多模電路等效為多端口網(wǎng)絡(luò )的方法。

　　選定的基函數滿(mǎn)足下述條件：只是空間坐標的函數；與時(shí)間無(wú)關(guān)；構成一個(gè)完備正交集。且對于給定的微波電路，選定的基函數應能夠有效地描述電路中電磁場(chǎng)的分布規律。假設：X-Y平面為電路橫截面所在平面，Z為傳播方向，電路在Dirac-δ函數激勵下在z=z0處的電場(chǎng)分布為 Ei（x,y,z0,t），{φmn（x,y）}為基函數族，用φmm（x,y）可將微波電路中t=t0,z=z0處的場(chǎng)表示為：

　　其中amn（z0,t0）為第（m,n）次基函數的系數，即幅度，這樣從參考面z=z0看入的微波電路可等效為一個(gè)基于基函數的等效時(shí)域多模電路。 基函數的函數形式既可以是適用于一般電路的正交函數形式，也可以是特別適用于某類(lèi)電路的特殊正交函數。一般說(shuō)來(lái)，當電路幾何結構比較復雜，不易根據電路特性選取特殊的正交函數作為基函數時(shí)，可以選取矩形脈沖函數（取網(wǎng)格結點(diǎn)的值作為整個(gè)網(wǎng)格的平均值，故脈沖寬度為一個(gè)網(wǎng)格的寬度）。但因脈沖函數描述的只是系統的局部信息，因此要達到足夠的精度，函數的展開(kāi)項數較多。當正交函數可以有效表述電路的全局信息時(shí)，通常只需幾項或十幾項，就可以達到所需的精度。例如，對于均勻填充的矩形波導問(wèn)題，如根據波導內的場(chǎng)的分布特性，把基函數選為{sin,cos}正交函數集，通常只需5項就可以滿(mǎn)足要求。相比較之下，至少需要60個(gè)脈沖即60個(gè)結點(diǎn)方可較準確地描述波導系統橫截面上的空間場(chǎng)分布。

　　基函數的正交性使得每一個(gè)基函數可以被視為一個(gè)獨立的端口，因此上述基于基函數的等效時(shí)域多模電路就可以進(jìn)一步被視為一個(gè)多端口網(wǎng)絡(luò )。

　　2.等效多端口網(wǎng)絡(luò )特性的計算

　　沖擊函數的頻譜是無(wú)限寬的，因此不能直接使用FDTD算法求解系統的沖擊響應函數。FDTD-Diakoptics使用高斯脈沖調制波作為激勵，通過(guò)加窗Fourier變換技術(shù)，求得有限帶寬微波電路的沖擊響應函數。其中，高斯脈沖激勵的調制頻率為電路工作頻帶的中心頻率，脈沖寬度和脈沖時(shí)間采樣間隔取決于頻率分辨率和帶寬。盡管激勵脈沖具有有限帶寬導致FDTD-Diakoptics求得的沖擊響應函數中包含了加窗帶來(lái)的影響（稱(chēng)此時(shí)的沖擊響應函數為準沖擊響應函數），但是只要滿(mǎn)足下述條件：使用FDTD-Diakoptics分析整個(gè)電路特性時(shí)，各個(gè)子電路使用具有相同頻譜特性的激勵脈沖，計入加窗對時(shí)域脈沖的展寬效應，最終得到的沖擊響應函數的頻域響應是足夠準確的。

　　五、FDTD-Diakoptics應用實(shí)例及討論

　　時(shí)域有限差分方法（Finite Difference Time Domain），簡(jiǎn)稱(chēng)FDTD.FDTD方法是把 Maxwell 方程式在時(shí)間和空間領(lǐng)域上進(jìn)行差分化。利用蛙跳式（Leap frog algorithm）--空間領(lǐng)域內的電場(chǎng)和磁場(chǎng)進(jìn)行交替計算，通過(guò)時(shí)間領(lǐng)域上更新來(lái)模仿電磁場(chǎng)的變化，達到數值計算的目的。用該方法分析問(wèn)題的時(shí)候要考慮研究對象的幾何參數，材料參數，計算精度，計算復雜度，計算穩定性等多方面的問(wèn)題。其優(yōu)點(diǎn)是能夠直接模擬場(chǎng)的分布，精度比較高，是目前使用比較多的數值模擬的方法之一。

　　本文以一個(gè)波導帶通濾波器的特性分析為例說(shuō)明該算法的應用。圖3為一個(gè)具有5個(gè)膜片的矩形波導帶通濾波器（WR34）。按照本方法首先將該濾波器分為5個(gè)部分，使用FDTD對其進(jìn)行計算，求出在所有連接參考面處（圖中虛線(xiàn)所示）的場(chǎng)分布。FDTD計算中，高斯脈沖調制函數為：f（t）=AmaxA（x,y）exp[-（（t- t0）/T）2].sin（2πf0t），其中調制頻率f0為WR34-TE10模單模工作頻帶的中心頻率；A（x,y）為激勵幅度空間分布，Diakoptics算法中需計算所有可能存在的基函數單一激勵時(shí)的響應，所以A（x,y）依次選為每一個(gè)基函數。激勵函數幅度Amax取決于其沿傳播方向的衰減及計算精度，基本原則是達到不連續性處和觀(guān)察面處的場(chǎng)仍具有足夠大的幅度。T的取值要保證在激勵脈沖的頻譜上截止頻率點(diǎn)處的能量具有足夠小的值。本例中，WR34的單模工作頻帶為：fTE10=17.369GHz,fTE20=34.738GHz,f0=26.0535GHz,T=200（ps），t0=3T,Amax=10,基函數為φn（x）=sin,相應的系數an（z0,t）如圖4所示（由于文章篇幅原因，只給出一個(gè)結果）。圖5為用本文方法得到的濾波器頻率特性，圖中可見(jiàn)該結果與FDTD結果吻合很好。

圖3　五膜片WR34波導帶通濾波器示意圖

　　　　　　　　　

圖4　本文方法得到的圖3中第一個(gè)子電路反射波基函數的系數

圖5　圖3所示W(wǎng)R34波導濾波器的頻率特性

　　六、結 論

　　本文介紹了一種分析復雜微波電路的新方法：FDTD-Diakoptics方法，它可克服傳統的FDTD方法需要大內存、長(cháng)計算時(shí)間的弊病，并可充分發(fā)揮FDTD可易于研究復雜幾何結構電路的優(yōu)勢，經(jīng)作者的若干分析設計實(shí)例證明，該方法不但比較靈活，且具有較高的精度，是一種比較有效的微波電路仿真分析方法。