周期性非均勻介質波導特征模分析的七點頻域有限差分方法技術

本發明專利技術公開了一種周期性非均勻介質波導特征模分析的七點頻域有限差分方法，主要步驟是：首先進行有限差分的網格剖分，使非均勻介質的分界面與剖分單元吻合，計算節點為剖分單元頂點；進而利用周期性邊界條件，通過等效諧振腔理論，將介質中電磁場實頻域的場方程轉化為復頻域的場方程；通過計算節點相鄰的8個剖分單元分界面上場的切向連續性邊界條件，利用與計算節點相鄰的7個計算節點插值實現非均勻介質的場方程離散；通過離散的場方程，利用特征值分解技術得到特征模式的頻率及場分布；建立波導各個模式的傳播常數與工作頻率的函數關系，從而指導周期性非均勻介質波導的設計。本發明專利技術解決了現有技術無法實現沿傳播方向介質分布非均勻的問題。

【技術實現步驟摘要】
周期性非均勻介質波導特征模分析的七點頻域有限差分方法
本專利技術屬于毫米波集成電路
，涉及電磁場數值分析中的不均勻介質分析方法，是基于頻域有限差分方法以及周期性邊界條件對非均勻介質波導傳播模式的分析，結合介質分界面邊界條件以及三維空間離散建立特征方程，可用于周期性毫米波集成電路中非均勻介質波導的模式分析，從而指導其系統設計，本專利技術具體是一種周期性非均勻介質波導特征模分析的七點頻域有限差分方法。
技術介紹
毫米波技術、系統及應用是工業界、學術界以及商業界長期以來最為關注的領域之一。其工作頻段介于微波及光波之間，兼有二者的優點，同時又有其獨特的特點。相對于微波波段，毫米波具有極寬的帶寬，對于目前日趨緊張的頻率資源顯得尤為引人注目；與此同時，毫米波天線具有極窄的波束寬度，因此可以提供更高的空間分辨率；最為重要的是，毫米波器件體積極小，其系統更易小型化、集成化，同時也帶來了質量輕的特點，提高其便攜性。而相對于光波，毫米波穿透霧、煙、灰塵的能力強，具有全天候全天時的特點。正是因為這些特點，毫米波技術得到了廣泛的應用，諸如可以提供更為細節信息的安檢系統，可用于提高制藥、材料工藝等其它技術的高精度分譜技術，高安全性的地面通信系統，高效的星際通信系統，高精度汽車防撞雷達系統，等等其它與民生息息相關的領域；與此同時，毫米波技術在國防領域也有著重要的應用，諸如可與微波雷達互補的毫米波雷達，高抗干擾能力的制導系統等等。雖然毫米波電路具有小型化的天生優勢，但相對于已有完善理論體系的微波電路，毫米波集成電路的理論體系尚不完善，迫切需要能夠完成復雜電路工作模式分析的技術。工程技術人員常用的商業軟件雖然能夠實現對復雜毫米波電路的場分析，但大多只能提供給定激勵模式下沿傳輸線的場分布或者網絡參數，難以提供特征模式的分析。因此，有眾多學者和機構基于計算電磁學技術，開展了針對毫米波集成電路傳輸線結構的特征模式分析。目前已開展的研究多是基于均勻介質傳輸線，或者橫截面非均勻傳輸線，然而目前毫米波集成電路中廣泛應用的基片集成波導結構以及基片集成介質波導結構都是在傳播方向上存在介質不均勻，雖然有學者開展過此類結構的模式分析，但都基于主模分析，無法分析多模問題。Karlhein等人在Finite-DifferenceAnalysisofRectangularDielectricWaveguideStructures文章中提出采用五點差分技術實現分層介質的模式分析，但無法適用于延傳播方向非均勻的結構。XuFeng等人在EquivalentResonantcavitymodelofperiodicguide-wavestructuresanditsapplicationinfinitedifferencefrequencydomainalgorithm一文中提出采用等效諧振腔模型結合頻域有限差分方法分析周期性金屬介質混合非均勻結構，但只適用于分析該類結構的主模。隨著毫米波集成技術的發展，多模問題已經是毫米波集成電路亟需解決的關鍵問題之一，然而對于沿傳播方向存在復雜周期性的非均勻介質結構，目前尚缺少適用的特征模式分析技術。
技術實現思路
本專利技術的目的是克服上述現有技術中存在的問題，提供一種可以分析類似于基片集成介質波導結構特征模式的方法，使得既能夠分析延傳播方向的非均勻結構，也能夠分析周期性非均勻結構，能夠適用于周期性非均勻介質邊界。為此，本專利技術提供了一種周期性非均勻介質波導特征模分析的七點頻域有限差分方法，其技術方案是：周期性非均勻介質波導特征模分析的七點頻域有限差分方法，包含以下步驟：01)對周期性非均勻介質波導結構在計算機中進行電磁建模：將波導體放置在直角坐標系中,使波導壁與坐標平面平行,且波導中填充的介質沿z軸方向呈現周期性,周期為d，建立三維坐標數組記錄外形及非均勻介質分布的三維坐標參數,同時建立介電常數數組記錄該波導結構中介電常數的空間分布；02)基于三維坐標數組以及介電常數數組對該波導結構進行離散剖分；03)對傳播常數β進行賦初值；04)選取傳播常數β的最大值,判斷待分析的傳播常數是否小于選取的最大傳播常數，如果滿足條件，執行步驟05)，否則跳轉至步驟11)；05)利用Floquet定理在周期性單元沿z軸方向的兩個邊界面上施加周期性邊界條件：選取周期性非均勻介質波導結構的一個周期性單元作為分析目標，若周期單元長度為d，單元中傳播的電磁場角頻率為ω，復傳播常數為γ＝α+jβ，其中α為衰減因子，β為傳播常數；由周期性邊界條件可知，其中和分別為點(x,y,z)處電場強度矢量和磁場強度矢量，且均為周期函數，即06)提取等效諧振腔模型：定義復頻率ω'＝ω+jαvg，其中vg為群速，利用復頻率定義復振幅矢量,利用步驟02)的離散結果以及步驟05)在z＝0和z＝d面施加周期性邊界條件進一步在周期性單元四壁施加理想導體邊界條件εx(l,0,n)＝εx(l,My,n)＝0，εy(0,m,n)＝εy(Lx,m,n)＝0，即可得到將周期性單元等效為諧振腔模型，其內部電磁場滿足復頻域Maxwell方程組07)建立離散FDFD方程組，利用步驟06)結果，將周期性單元內部任意點P處滿足的復頻域Maxwell方程組，利用七點差分技術將方程組中的微分用差分表示，此時可得其中，下標F、B、L、R、U、D、P分別表示前、后、左、右、上、下及自身7個節點；遍取周期性單元的所有計算節點，可得矩陣方程其中[x]為計算節點處磁場強度矢量z及x方向的分量，[A]為上式形成的系數矩陣，顯然該矩陣方程為矩陣[A]的特征值方程；08)提取特征值及特征向量：對步驟07)形成的特征值方程進行特征值分解，得到矩陣[A]的復特征值及其對應的特征向量[x]；09)提取實特征頻率：利用步驟08)中的復特征值得到實特征頻率其中c為真空中的光速，即c＝3×108m/s；10)增加傳播常數：根據選定的傳播常數增量步進，增大傳播常數，增大傳播常數后轉至步驟04)；11)建立傳播常數隨特征頻率變化關系，利用步驟05)至09)，得到不同傳播常數下周期性單元所能激勵起的所有模式的特征頻率，從而得傳播常數與特征頻率變化關系的曲線。上述步驟02)中所述的基于三維坐標數組以及介電常數數組對該波導結構進行離散剖分，具體是將周期性非均勻介質波導結構剖分為正六面體單元的組合,若該結構工作在頻帶f1～f2之間，則采用邊長不超過λ2/20的正六面體單元進行剖分，其中λ2為f2所對應電磁波波長；在介電常數發生變化區域的剖分單元，必須保證頂點位于非均勻介質分界面，剖分單元的頂點即為計算節點，若該周期性單元在x方向剖分為Lx個單元，在y方向剖分為My個單元，在z方向剖分為Nz個單元，則節點序號從(0,0,0)開始，終止于(Lx,My,Nz)，其間的任意計算節點表示為(l,m,n)，l,m,n分別表示該計算節點在x、y、z方向上與起始計算節點相距l×d，m×d，n×d；此時可用周期性單元表面及內部所有計算節點處電場強度及磁場強度離散表示整個周期性單元的電磁場分布；上述步驟03)中，當分析傳播常數與頻率的關系時，取傳播常數的初值為0rad/s。上述步驟04)中，當分析傳播常數與頻率的關系時，取傳播常數的最大值為1000rad/s。本文檔來自技高網...

【技術保護點】
周期性非均勻介質波導特征模分析的七點頻域有限差分方法，其特征在于，包含以下步驟：01)對周期性非均勻介質波導結構在計算機中進行電磁建模：將波導體放置在直角坐標系中,使波導壁與坐標平面平行,且波導中填充的介質沿z軸方向呈現周期性,周期為d，建立三維坐標數組記錄外形及非均勻介質分布的三維坐標參數,同時建立介電常數數組記錄該波導結構中介電常數的空間分布；02)基于三維坐標數組以及介電常數數組對該波導結構進行離散剖分；03)對傳播常數β進行賦初值；04)選取傳播常數β的最大值,判斷待分析的傳播常數是否小于選取的最大傳播常數，如果滿足條件，執行步驟05)，否則跳轉至步驟11)；05)利用Floquet定理在周期性單元沿z軸方向的兩個邊界面上施加周期性邊界條件：選取周期性非均勻介質波導結構的一個周期性單元作為分析目標，若周期單元長度為d，單元中傳播的電磁場角頻率為ω，復傳播常數為γ＝α+jβ，其中α為衰減因子，β為傳播常數；由周期性邊界條件可知，其中和分別為點(x,y,z)處電場強度矢量和磁場強度矢量，且均為周期函數，即06)提取等效諧振腔模型：定義復頻率ω'＝ω+jαvg，其中vg為群速，利用復頻率定義復振幅矢量,利用步驟02)的離散結果以及步驟05)在z＝0和z＝d面施加周期性邊界條件進一步在周期性單元四壁施加理想導體邊界條件εx(l,0,n)＝εx(l,My,n)＝0，εy(0,m,n)＝εy(Lx,m,n)＝0，即可得到將周期性單元等效為諧振腔模型，其內部電磁場滿足復頻域Maxwell方程組07)建立離散FDFD方程組，利用步驟06)結果，將周期性單元內部任意點P處滿足的復頻域Maxwell方程組，利用七點差分技術將方程組中的微分用差分表示，此時可得其中，下標F、B、L、R、U、D、P分別表示前、后、左、右、上、下及自身7個節點；遍取周期性單元的所有計算節點，可得矩陣方程其中[x]為計算節點處磁場強度矢量z及x方向的分量，[A]為上式形成的系數矩陣，顯然該矩陣方程為矩陣[A]的特征值方程；08)提取特征值及特征向量：對步驟07)形成的特征值方程進行特征值分解，得到矩陣[A]的復特征值及其對應的特征向量[x]；09)提取實特征頻率：利用步驟08)中的復特征值得到實特征頻率其中c為真空中的光速，即c＝3×108m/s；10)增加傳播常數：根據選定的傳播常數增量步進，增大傳播常數，增大傳播常數后轉至步驟04)；11)建立傳播常數隨特征頻率變化關系，利用步驟05)至09)，得到不同傳播常數下周期性單元所能激勵起的所有模式的特征頻率，從而得傳播常數與特征頻率變化關系的曲線。...

【技術特征摘要】
1.周期性非均勻介質波導特征模分析的七點頻域有限差分方法，其特征在于，包含以下步驟：01)對周期性非均勻介質波導結構在計算機中進行電磁建模：將波導體放置在直角坐標系中,使波導壁與坐標平面平行,且波導中填充的介質沿z軸方向呈現周期性,周期為d，建立三維坐標數組記錄外形及非均勻介質分布的三維坐標參數,同時建立介電常數數組記錄該波導結構中介電常數的空間分布；02)基于三維坐標數組以及介電常數數組對該波導結構進行離散剖分；所述的基于三維坐標數組以及介電常數數組對該波導結構進行離散剖分，具體是將周期性非均勻介質波導結構剖分為正六面體單元的組合,若該結構工作在頻帶f1～f2之間，則采用邊長不超過λ2/20的正六面體單元進行剖分，其中λ2為f2所對應電磁波波長；在介電常數發生變化區域的剖分單元，必須保證頂點位于非均勻介質分界面，剖分單元的頂點即為計算節點，若該周期性單元在x方向剖分為Lx個單元，在y方向剖分為My個單元，在z方向剖分為Nz個單元，則節點序號從(0,0,0)開始，終止于(Lx,My,Nz)，其間的任意計算節點表示為(l,m,n)，l,m,n分別表示該計算節點在x、y、z方向上與起始計算節點相距l×d，m×d，n×d；此時可用周期性單元表面及內部所有計算節點處電場強度及磁場強度離散表示整個周期性單元的電磁場分布；03)對傳播常數β進行賦初值；04)選取傳播常數β的最大值,判斷待分析的傳播常數是否小于選取的最大傳播常數，如果滿足條件，執行步驟05)，否則跳轉至步驟11)；05)利用Floquet定理在周期性單元沿z軸方向的兩個邊界面上施加周期性邊界條件：選取周期性非均勻介質波導結構的一個周期性單元作為分析目標，若周期單元長度為d，單元中傳播的電磁場角頻率為ω，復傳播常數為γ＝α+jβ，其中α為衰減因子，β為傳播常數；由周期性邊界條件可知，其中和分別為點(x,y,z)處電場強度矢量和磁場強度矢量，且均為周期函數，即06)提取等效諧振腔模型：定義復頻率ω'＝ω+jαvg，其中vg為群速，利用復頻率定義復振幅矢量,利用步驟...

【專利技術屬性】
技術研發人員：徐樂，李蕊，史小衛，王俊，仁學施，
申請(專利權)人：西安電子科技大學，
類型：發明
國別省市：陜西;61