基于矢量擬合的過孔等效電路提取方法

孔 繁 盛衛(wèi)星 韓玉兵 王 昊 馬曉峰

(南京理工大學(xué)電子工程與光電技術(shù)學(xué)院，江蘇 南京 210094)

引 言

隨著高速集成電路的快速發(fā)展，信號完整性問題[1]變得越來越復(fù)雜，單純的電路仿真不再能滿足精度的需要，而完整的電磁場分析建模又會使問題變得非常巨大，無論是對硬件中央處理器(Central Processing Unit,CPU)和內(nèi)存的需求上，還是電磁場問題本身的復(fù)雜度和困難度上，都會大幅度增加解決問題的成本.如何在效率和精確性上取得平衡，在既能滿足精度要求，又能快速有效地進行信號完整性仿真，甚至更好地兼容電路域的信號完整性仿真正逐步成為信號完整性分析方法研究的熱點.高速電路的不連續(xù)情況會嚴重影響信號傳播的質(zhì)量，而分析這些不連續(xù)情況則是電磁場建模中最為復(fù)雜的.例如印刷電路板上最常見的過孔結(jié)構(gòu)，它的形狀結(jié)構(gòu)復(fù)雜，穿越電路板平面層數(shù)多，為了得到精確的模型，其相對應(yīng)的建模也必須是精細的，尤其是帶有色散材質(zhì)的電磁分析，其計算量和分析過程是相當巨大的.但實際使用當中又常常希望能將這樣精細的建模的結(jié)果用于信號完整性分析中，甚至?xí)r域分析當中去.

現(xiàn)階段的很多信號完整性分析軟件(如ADS，HSPICE等等)都使用等效電路模型對高速互聯(lián)線進行快速的頻域和時域分析.使用電路模型仿真和使用全波的電磁模型仿真相比較有兩個優(yōu)勢：第一，可以有效地和傳統(tǒng)的電子器件模型兼容；第二，使用電路模型仿真對計算機的硬件資源要求比較低，節(jié)省對CPU和內(nèi)存的消耗.假如對整個電路板都采用全波建模仿真，其結(jié)果可以達到最精確，但是建模目標的相對物理尺寸跨度很大，建模的復(fù)雜度會非常高，仿真所需要的硬件資源會異常巨大，同時也會消耗很長的時間.在信號完整性分析的階段需要能夠在精確性和效率上做出平衡，在這樣的需求下，一種可行的方法是只針對感興趣的部分例如過孔結(jié)構(gòu)這樣的局部采用測量或者全波的電磁場分析，得到其特性參數(shù)，進而獲得其行為描述的宏模型或者稱之為黑箱模型，將其看作是兩端口或者是多端口的器件.這樣的宏模型建立要遵循以下的原則：首先是因果性條件，一個實際器件有兩端口，當信號從一端輸入，另一端口要滯后一段時間輸出信號.第二是穩(wěn)定性條件要求如果是帶限信號激勵，輸出的還是帶限信號.因此用一個有理函數(shù)來代替穩(wěn)定系統(tǒng)需要滿足極點在S平面的左半平面,有理函數(shù)不能在虛軸上有多重極點.第三是無源性，無源性條件要求無源電路不產(chǎn)生能量，用系統(tǒng)函數(shù)H(s)來解釋，就是H(s)不含右半平面的極點,在虛軸上沒有多重極點,系數(shù)是實數(shù),H(s)的實部對于所有的頻率都大于或等于零，即，Re(H(s))≥0.可以用來擬合宏模型的方法主要有基于Pade近似[2]的矩匹配法[3-4]、漸近波形估計法[3]、Krylov子空間變換法[6]、最小二乘估計法[7]、矢量擬合法[8-10].其中基于Pade近似的矩匹配法對矩的變化特別敏感容易不穩(wěn)定，不能很好滿足第二點的要求，而漸近波形估計法和Krylov子空間變換法無法解決第三點無源性的問題，最小二乘法在零極點個數(shù)比較多的情況下無法應(yīng)用.然而矢量擬合法都可以很好地滿足以上三個要求.矢量擬合法將網(wǎng)絡(luò)模型擬合為極點-留數(shù)型有理函數(shù)的形式，在擬合過程中保證極點都具有負的實部，滿足系統(tǒng)穩(wěn)定的條件.而且擬合成的零點-極點型有理函數(shù)還可以利用阻抗綜合技術(shù)得到直觀的電阻、電感、電容的電路形式，這是其他方法都無法做到的，方便了后續(xù)用于時域方法[11]的仿真，也方便加入到電路仿真軟件中仿真分析.但是，已有的采用矢量擬合法建立宏模型的方法通過有理函數(shù)擬合建立RL等效電路[12]來實現(xiàn)參數(shù)向電路的轉(zhuǎn)換有一定的局限性，僅是通過對二端口模型spice型電路[13]建模走線的等效電路方法進行電路模型的建立，不能解決差分對過孔這樣的復(fù)雜結(jié)構(gòu).我們在使用矢量擬合的基礎(chǔ)上提出了新的模型形式，從而對單個過孔或者差分對過孔都能很好地擬合出電路模型.

在分析過孔問題時，電路模型的建立[14]是根據(jù)過孔的高度、孔徑尺寸、焊盤半徑等參數(shù)估算出過孔的電感值和焊盤的電容值來進行粗糙的過孔建模，此類方法的缺點是對于某些隨頻率變化的參數(shù)，在模型中無法考慮，在高頻情況下，此類模型不準確.使用準靜的矩量法計算過孔的等效電路模型[15]，能在較低的頻段下滿足仿真需要，得到的是簡單的過孔準靜的電感和電容的計算公式，組成簡單的П型等效電路模型，針對這種簡單П型等效電路可更進一步改進過孔等效電路模型[16-17]，當加入了地平面間的耦合電容，使得仿真的頻率有所提高，但過孔電感和焊盤電容的計算方式，仍可使用低頻下的計算公式，其地平面間的耦合電容公式是基于特定尺寸下測量值擬合出來的，適用范圍會受到限制.

本文提出了先精細建模得到過孔頻域參數(shù)再將其通過矢量擬合建立成等效電路的方法來解決單個過孔和差分對過孔的信號完整性分析方法.該方法根據(jù)過孔的全波模型，推導(dǎo)了一般的單個過孔的電路等效模型提取方法，并且推廣到差分對過孔的等效電路，該等效電路模型不僅可以精確仿真過孔的行為，還具有實際電路的物理意義.通過與專業(yè)電磁場分析軟件的分析結(jié)果對比，驗證了本文提取的等效電路的正確性，為大規(guī)模仿真含有過孔結(jié)構(gòu)的信號完整性仿真提供了有效的方法.

1 理論分析

建模的主要思想是通過全波仿真軟件精細建模過孔的模型，第一步先得到過孔的導(dǎo)納參數(shù)或者散射參數(shù).第二步導(dǎo)納參數(shù)(Y參數(shù))經(jīng)過變換可以得到用于等效電路擬合的參數(shù)，對于散射參數(shù)(S參數(shù))的情況經(jīng)過變換為導(dǎo)納參數(shù)之后同樣可以推導(dǎo)出用于等效電路擬合的參數(shù).第三步將等效電路擬合之后的參數(shù)轉(zhuǎn)換為相應(yīng)的具有物理意義的等效電路形式.

圖1 單個過孔等效電路

通過全波分析軟件得到導(dǎo)納參數(shù)之后，建立等效電路的時候是沒有直接的對應(yīng)關(guān)系的.所以首先要建立從Y參數(shù)矩陣提取П型等效電路的方法，推導(dǎo)出用于擬合的參數(shù)形式，如圖1所示，對于單個過孔的結(jié)構(gòu)可以看作一個二端口網(wǎng)絡(luò).對于一般的二端口網(wǎng)絡(luò)Y參數(shù)為：

I1=Y11·U1+Y12·U2；

I2=Y21·U1+Y22·U2.

(1)

如圖1所示，可以得到相應(yīng)用于擬合的參數(shù)YA、YB、YC與電流、電壓的對應(yīng)關(guān)系如下所示：

I1=YA·U1+YB·(U1-U2)；

I2=YB·(U2-U1)+YC·U2.

(2)

對于不包含有源器件結(jié)構(gòu)的網(wǎng)絡(luò)，如過孔這樣的結(jié)構(gòu)是可以看作互易網(wǎng)絡(luò)的，而互易網(wǎng)絡(luò)有Y12=Y21，所以可以推導(dǎo)出:

YA=Y11+Y12;

YC=Y22+Y12;

YB=-Y12=-Y21.

(3)

對于差分對過孔，可以建立如圖2的等效電路結(jié)構(gòu)來映射差分對過孔的Y參數(shù)矩陣.

圖2 差分對過孔等效電路

差分過孔結(jié)構(gòu)包含了連接兩條走線的兩個過孔，每條走線的直通通路在1、3端口之間用一個П型等效電路來表示，相互靠近的近端端口(1、2端口)之間也使用一個П型電路來等效，對于遠端端口(1、4端口)之間為了表示出其耦合關(guān)系也用一個П型電路來等效.建立起圖2所示的等效電路后要推導(dǎo)出一般的四端口Y參數(shù)與之對應(yīng)的表達式.這樣的等效電路結(jié)構(gòu)包含了差分對所有的耦合信息，一般的差分對過孔的網(wǎng)絡(luò)電壓與電流的對應(yīng)關(guān)系如下：

I1=Y11·U1+Y12·U2+Y13·U3+Y14·U4；

I2=Y21·U1+Y22·U2+Y23·U3+Y24·U4；

I3=Y31·U1+Y32·U2+Y33·U3+Y34·U4；

I4=Y41·U1+Y42·U2+Y43·U3+Y44·U4.

(4)

對應(yīng)圖2的結(jié)構(gòu)可以推出其對應(yīng)端口用于擬合的參數(shù)的電壓與電流的表達式：

I1= (YA+YE+YI+YH)·U1+(-YE)·U2+

(-YH)·U3+(-YI)·U4，

I2=(-YE)·U1+(YB+YE+YF+YJ)·U2+

(-YJ)·U3+(-YF)·U4，

I3=(-YH)·U1+(-YJ)·U2+

(YD+YH+YG+YJ)·U3+(-YG)·U4，

I4=(-YI)·U1+(-YF)·U2+(-YG)·U3+

(YC+YF+YG+YI)·U4.

(5)

同樣對于差分對過孔這樣無源的結(jié)構(gòu)，可以將其看作互易網(wǎng)絡(luò)處理，同時對照式(4)與式(5)，可以推導(dǎo)出具體的用于擬合的參數(shù)表達式：

YA=Y11+Y12+Y13+Y14，

YB=Y21+Y22+Y23+Y24，

YC=Y31+Y32+Y33+Y34，

YD=Y41+Y42+Y43+Y44，

YE=-Y12=-Y21，

YF=-Y24=-Y42，

YG=-Y34=-Y43，

YH=-Y13=-Y31，

YI=-Y14=-Y41，

YJ=-Y23=-Y32.

(6)

由此已經(jīng)得到了全部的用于擬合參數(shù)的具體表達式，通過將其擬合成極點-留數(shù)型函數(shù)之后就可以對應(yīng)的找出電路元件表達的具體形式.矢量擬合可以采用文獻[8]的方法，將對應(yīng)的參數(shù)進行擬合，得到有理函數(shù)的形式,即

(7)

式中：s=jω，ω是角頻率； 常數(shù)項d和一次項e是實數(shù)；ak是極點，一般為實數(shù)或者共軛復(fù)數(shù)對；ck是第k個極點對應(yīng)的留數(shù)；N表示擬合的全部極點的個數(shù).對于常數(shù)項d和一次項e，可以對應(yīng)為電阻和電容，電阻的阻值為1/d，電容值是e.經(jīng)過擬合之后的N個有理函數(shù)的極點均在左半平面，經(jīng)過整理可以分為M個實極點和(N-M)/2對共軛復(fù)極點的情況.

對于實數(shù)極點的情況可以表示為式

(8)

通過圖3所示的電感與電阻串聯(lián)支路的形式可以對應(yīng)地找出電感與電阻的值.于是可以將Yr用導(dǎo)納的形式寫為

圖3 RL實數(shù)極點-留數(shù)對應(yīng)等效電路

(9)

對比式(8)與式(9)，于是可以推出電阻、電感與極點和留數(shù)的關(guān)系：

(10)

對于極點和留數(shù)是共軛復(fù)數(shù)對的情況，等效電路形式如圖4所示.

圖4 共軛復(fù)數(shù)對對應(yīng)的等效電路

為了找出RLC電路元件具體的數(shù)值，可以將共軛復(fù)數(shù)對的極點-留數(shù)有理函數(shù)的形式表示為

(11)

由圖4中的電路形式，其在復(fù)頻域的導(dǎo)納參數(shù)可以寫為

(12)

對比式(11)與(12)可以計算出所需要的具體的電路參數(shù)值:

(13)

至此已經(jīng)可以根據(jù)擬合出來的零極點情況構(gòu)成局部等效電路，再由局部等效電路構(gòu)成П型等效電路模型.

圖5 綜合的局部等效電路

前面推導(dǎo)出的都是導(dǎo)納參數(shù)，所以電路是并聯(lián)關(guān)系，將M個實數(shù)極點和(N-M)/2個共軛復(fù)極點的電路全部并聯(lián)起來組成擬合出來的局部等效形式，如圖5所示.例如二端口網(wǎng)絡(luò)的П等效電路有三個局部等效電路，需要進行三次擬合，形成圖5形式的三條支路.如果是差分對過孔的結(jié)構(gòu)則需要擬合出局部等效電路十個，將YA至YJ的所有支路構(gòu)成局部等效電路組成如圖2所示的電路形式.這樣就完成了過孔等效電路的建模，之后可以方便地加入信號完整性電路仿真軟件中進行仿真.

2 實驗結(jié)果分析

使用本文推導(dǎo)的單個過孔電路模型提取方法，先對單個過孔的情況進行了仿真.在HFSS中對單個過孔的建模如圖6所示.

圖6 單個過孔的仿真模型

過孔穿越地平面需要有反焊盤保護來防止短路，一般反焊盤的半徑會大于焊盤半徑，這里建模時選擇了焊盤半徑rpad為0.7 mm，反焊盤半徑ranti比焊盤半徑略大為0.9 mm，過孔的半徑rvia為0.6 mm.整個過孔的高度h是0.63 mm，連接上下兩段微帶傳輸線，為了仿真需要微帶線不可能無限長，截取長度L為5 mm，線寬ω設(shè)定為0.6 mm，微帶結(jié)構(gòu)填充介質(zhì)相對介電常數(shù)為4.4.先通過這樣的精細建模仿真得到Y(jié)參數(shù)矩陣，再通過第一節(jié)中推導(dǎo)的公式，計算出需要擬合的YA、YB、YC參數(shù)，通過矢量擬合將這些參數(shù)擬合為極點-留數(shù)型有理函數(shù).然后根據(jù)得到的極點和留數(shù)構(gòu)造出等效電路，最后使用電路模型進行S參數(shù)的仿真，電路域的仿真結(jié)果再同HFSS仿真的S參數(shù)結(jié)果比較，從而驗證本文單個過孔等效電路的正確性.

針對需要擬合的YA、YB等參數(shù)，圖7(a)、(b)給出了矢量擬合的結(jié)果，可以看到原始YA或者YB參數(shù)的實部和虛部曲線在整個帶寬內(nèi)都與擬合的有理函數(shù)穩(wěn)和非常好，這說明了提取出的極點-留數(shù)型有理函數(shù)可以很好地模擬相應(yīng)參數(shù)在頻段內(nèi)的行為特性.為了縮減篇幅省略了YC參數(shù)的擬合效果圖，實際計算的YC參數(shù)擬合效果與YA，YB的一樣，得到的有理函數(shù)都能很好地還原原始參數(shù)的行為特性，為進一步提取出等效電路提供了保證.對應(yīng)有理函數(shù)的極點和留數(shù)情況，使用式(10)、(13)可以計算出相應(yīng)的電路參數(shù)，構(gòu)成等效的電路，在電路域仿真S參數(shù)的結(jié)果可以在圖7(c)中看到，單個過孔等效電路的S參數(shù)與全波專業(yè)電磁仿真軟件HFSS的仿真結(jié)果一致，很好地證明了單個過孔等效電路的有效性.

(a) YA擬合結(jié)果

(b) YH擬合結(jié)果

(c) 等效電路與HFSS仿真的S參數(shù)比對結(jié)果圖7 單個過孔等效電路仿真結(jié)果驗證

差分對過孔(如圖8)的走線線寬為0.1 mm，線間距Sp為0.1 mm，兩個過孔之間的間距Vs為1.5 mm，由于過孔有焊盤和反焊盤，走線到連接過孔的位置就必須采用分叉的連接增加了線間距，過孔的半徑rvia為0.15 mm，焊盤半徑rpad為0.2 mm，反焊盤半徑ranti為0.4 mm，過孔連接上下兩段微帶線，過孔的整體高度h為0.45 mm，截取的走線長度L為2.6 mm，微帶線填充介質(zhì)相對介電常數(shù)4.4.

圖8 差分對過孔的仿真模型

(a) YA參數(shù)擬合結(jié)果

(b) YH參數(shù)擬合結(jié)果

(c) YI參數(shù)擬合結(jié)果

(d) 等效電路與HFSS仿真的S參數(shù)比對結(jié)果圖9 差分對過孔等效電路仿真結(jié)果驗證

由圖9(a)～(c)部分的參數(shù)擬合結(jié)果可以看出，使用矢量擬合針對公式(6)中的YA、YH、YI參數(shù)的擬合效果很好.將擬合的有理函數(shù)按照之前計算等效電路的方法建立起等效電路，再仿真差分對過孔的S參數(shù)，由圖9(d)可以看出HFSS仿真結(jié)果和等效電路建立的模型的S參數(shù)仿真結(jié)果一致，因為在HFSS中針對局部的建模是精細的建模，因此差分過孔電路模型很好地保存了這些物理特性，在行為上完全能夠模擬過孔的傳輸特性，同時等效電路可以方便地加入任何電路仿真軟件中和其他電路元器件聯(lián)合仿真，適用于大規(guī)模的信號完整性仿真.

3 結(jié) 論

通過對二端口和四端口網(wǎng)絡(luò)參數(shù)的推導(dǎo)，導(dǎo)出了適用于建立單個過孔和差分對過孔的用于擬合的參數(shù).再通過矢量擬合將相應(yīng)的參數(shù)擬合成極點-留數(shù)型有理函數(shù)的形式，將抽取出來的極點和留數(shù)映射為相應(yīng)的RLC局部等效電路，最后把局部等效電路組合成完整的等效電路模型.該模型保留了精細建模過孔的物理特性，在行為仿真上能夠替代全波模型用于電路仿真器中.文中提出的等效電路提取方法是針對使用精細建模方法或者實驗測量的基礎(chǔ)上，對已經(jīng)得到的單個過孔或者差分過孔的導(dǎo)納參數(shù)進行擬合，使用該方法的頻率限制取決于精細建模方法所能達到精確計算結(jié)果的限制或者實驗測量儀器的帶寬限制.所以使用該方法進行正確建模的前提是得到精確的導(dǎo)納參數(shù)，該導(dǎo)納參數(shù)的頻帶范圍應(yīng)該根據(jù)具體問題而分析.在信號完整性分析中一般認為信號的最高頻分量為上升時間或者下降時間的80%，相應(yīng)的仿真所需精確的導(dǎo)納參數(shù)的帶寬應(yīng)為最高頻率分量的3到5倍.例如本文中的10 GHz的過孔模型分析，適用于分析的信號上升或者下降時間不超過0.625 ns.相應(yīng)的如果要進行信號完整性分析中的信號速度越快，上升/下降時間越短，則需要對更高頻率的導(dǎo)納參數(shù)進行計算或者測量.通過仿真驗證，過孔的等效電路模型仿真結(jié)果與HFSS仿真結(jié)果一致，證明了等效電路模型的有效性和正確性，該模型可以用于大規(guī)模的信號完整性分析.

[1] GREEN L. Understanding the importance of signal integrity[J]. IEEE Circuits and Devices Magazine, 1999, 15(6): 7-10.

[2] NATORI K. A design method of time-delay systems with communication disturbance observer by using Pade approximation[C]//12th IEEE International Workshop on Advanced Motion Control (AMC). Sarajevo, March 25-27, 2012: 1-6.

[3] 秦建軍, 殷兆偉, 曹祥玉, 等. 漸近波形估計技術(shù)在線天線分析中的應(yīng)用[J]. 電波科學(xué)學(xué)報, 2005, 20(6): 824-828.

QIN Jianjun, YIN Zhaowei, CAO Xiangyu, et al. Applications of asymptotic waveform evaluation technique in the analysis of wire antenna[J]. Chinese Journal of Radio Science, 2005, 20(6): 824-828.(in Chinese)

[4] IONESCU T C, ASTOLFI A. Moment matching for linear port Hamiltonian systems[C]// 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC). Orlando, December 12-15, 2011: 7164-7169.

[5] WU Jiang. Simultaneous computation of wide-angle and wide-band base on asymptotic waveform evaluation[C]// 2012 International Conference on Computer Science and Electronics Engineering (ICCSEE), March 2012, 2: 359-362

[6] JAGELS C, REICHEL L. Recursion relations for the extended Krylov subspace method[J]. Original Research Article Linear Algebra and its Applications, 2011, 434(7): 1716-1732.

[7] 朱凱然, 何學(xué)輝, 鄭小保, 等. 基于總體最小二乘的核四極矩共振參數(shù)估計[J]. 電波科學(xué)學(xué)報, 2012, 27(1): 134-140.

ZHU Kairan, HE Xuehui, ZHENG Xiaobao, et al. The estimation of NQR parameters based on total least square[J]. Chinese Journal of Radio Science, 2012, 27(1): 134-140.(in Chinese)

[8] GUSTAVSEN B, SEMLYEN A. rational approximation of frequency domain responses by vector fitting[J].IEEE Trans Power Delivery, 1999, 14(3): 1052-1061.

[9] GUSTAVSEN B. Improving the pole relocating properties of vector fitting[J]. IEEE Trans Power Delivery, 2006, 21(3): 1587-1592.

[10] DESCHRIJVER D, MROZOWSKI M, DHAENE T, et al. Macromodeling of? multiport systems using a fast implementation of the vector fitting method[J].IEEE Microwave and Wireless Components Letters, 2008, 18(6): 383-385.

[11] HOLLOWAY C L, MCKENNA P, JOHNK R T. The effects of gaps in ferrite on both absorber and chamber performance[C]// Proc IEEE 1999 Int Symp Electromagnetic Compatibility. Seattle, August 2-6, 1999: 239-244.

[12] SARTO M S, SCARLATTI A, HOLLOWAY C. On the use of fitting models for the time-domain analysis of problems with frequency-dependent parameters[C]// Proc IEEE 2001 Int Symp Electromagnetic Compatibility. Montreal, August 13-17, 2001: 588-593.

[13]ANTONINI G. SPICE equivalent circuits of frequency-domain responses[J]. IEEE Transactions on Electromagnetic Compatibility, 2003, 45(3): 502-512.

[14] 劉學(xué)觀, 江盼盼, 郭輝萍, 等. 高速數(shù)字信號在PCB 中的傳輸特性分析[J]. 電波科學(xué)學(xué)報, 2009, 24(3): 498-517

LIU Xueguan, JIANG Panpan, GUO Huiping, et al. Transmission characteristics of high-speed digital signal in PCB[J]. Chinese Journal of Radio Science, 2009, 24(3): 498-517.(in Chinese)

[15] JOHNSON H, GRAHAM M. High-Speed Digital Design-A Handbook of Black Magic[M]. Saddle River: Prentice Hall, 1993.

[16] ANTONINI G, SCOGNA A C, ORLANDI A. Full-wave modeling of via holes and equivalent circuit extraction for signal integrity analysis[C]//Proceedings of International Symposium on Electromagnetic Compatibility (EMC Europe 2002). Sorrento, September 9-13, 2002.

[17] NGUYEN T V, MORALES A, AGILI S. Characterization of differential via holes using equivalent circuit extraction technique[J]. International Journal of RF and Microwave Computer-Aided Engineering, 2007, 17(6): 552-559