【導(dǎo)讀】近年來，雙負(fù)媒質(zhì)材料的研究成為研究合成材料的熱點。鑒于其異于傳統(tǒng)媒質(zhì)的典。型特性，其應(yīng)用十分廣泛。在信息化的今天，快速計算雙負(fù)媒質(zhì)散射體的電磁散射特性，對實現(xiàn)雙負(fù)媒質(zhì)的寬頻帶的電磁散射特性分析很有現(xiàn)實意義和價值。然后介紹PMCHWT方程推導(dǎo)過程和矩量法的基本理論，闡述漸近波。負(fù)媒質(zhì)的PMCHWT方程，從而解出單頻率入射波下的電流和雷達(dá)散射截面。興趣的頻段內(nèi)能很好逼近矩量法逐點求解的結(jié)果，而且在計算效率上也大大提高。

　　

【正文】 的電流，通過 Pad233。 逼近拓展有效計算帶寬 0000()()()Lllnln Mmmnma k kIkb k k??????? （ 323） LM、 是 Pad233。 逼近階數(shù)。聯(lián)立式 (321)與 (323)，整理對比 0kk?（ ） 同類項系數(shù)可得系數(shù) lmnnab、 ，從而可以計算出任意頻點的極化電流和 RCS 。 從上面的推導(dǎo)可知， AWE 技術(shù)矩陣的求逆次數(shù)只有一次，特別對雙負(fù)媒質(zhì)而言，由于阻抗矩陣元素的復(fù)雜性，避免了矩量法多次求逆的弊端，將大幅提高計算效率。 二維散射體阻抗矩陣及阻抗元素求導(dǎo) 從雙負(fù)（ DNG）媒質(zhì)的 PMCHWT 方程的推導(dǎo)過程可以清楚看出，雙負(fù)媒質(zhì)的阻抗元素在外等效時和普通介質(zhì)目標(biāo)沒有區(qū)別，因此對外等效阻抗矩陣 00AB、 不作討論，以下主要研究內(nèi)等效時的阻抗矩陣 00CD、 。內(nèi)等效時阻抗元素的具體表達(dá)式為： ? ?( 1 ) 21 10 1 0 1001 110( ) ( r ) 4( , )21 l n , 44k C H k j H k r m nC m nk k CjC m ne??? ???? ??? ? ?????? ?? ? ?? ? ? ?? ? ? ?? ? ? ???? ? ? ? ????（ ）（ ） ， （ 331）( 1 ) 211 1 1 10( ) ( ) , 4( , ) 1 , 2rjk C H k r jH k r n e m nD m nmn? ??? ? ? ? ????? ?? ???（ ）（ ） （ 332） 其中 1? 為雙負(fù)媒質(zhì)內(nèi)波阻抗， 0? 為真空中波阻抗， 1k 為雙負(fù)媒質(zhì)中的波數(shù)，Δ C 為離散步長， r 為場點與源點距離， n 為電流元的單位外法向矢量，re為場點至源點的單位矢量。 由于 ??、 都是關(guān)于頻率的復(fù)雜函數(shù)，如果將雙負(fù)媒質(zhì)的本構(gòu)關(guān)系式代入到阻抗元素表 22 達(dá)式后，直接對阻抗元素進行求導(dǎo)，求導(dǎo)工作將非常繁瑣，高階導(dǎo)數(shù)將 很難實現(xiàn)。采用如下處理，令 ? ?? ? ? ?110010( 1 ) ( 2 )0 0 1 0 1421 l n4kgCk d ljhep H k r jH k r?????? ?????????????????? ?????????? 則 000 00 , ( , ) , g p m nC m n g h m n??? ? ?? (333) 然后利用組合的思想對阻抗元素進行求導(dǎo)，從而簡化了求導(dǎo)的復(fù)雜過程。得出阻抗元素及其各階導(dǎo)數(shù)的矩陣表達(dá)如下： 76543210010 1 20 1 2 30 1 2 3 40 1 2 3 4 50 1 2 3 4 5 60 1 2 3 4 5 6 700 0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 20 0 0 0 3 30 0 0 4 6 40 0 5 10 10 50 6 15 20 15 67 21 35 35 21 7gpgppgp p pgp p p pp p p p p gp p p p p p gp p p p p p pgp p p p p p p pg???? ???? ?????????? ?????????? ?????? ????? ?? ?? ?? ?? ?? ?? ?? ?01234567, , ,C m nC m nC m nC m nmnC m nC m nC m nC m n???????????????????????????? (334) 76543210010 1 20 1 2 30 1 2 3 40 1 2 3 4 50 1 2 3 4 5 60 1 2 3 4 5 6 700 0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 20 0 0 0 3 30 0 0 4 6 40 0 5 10 10 50 6 15 20 15 67 21 35 35 21 7ghghhgh h hgh h h hh h h h h gh h h h h h gh h h h h h hgh h h h h h h hg???? ?????????????? ???????? ?????? ????? ?? ?? ?? ?? ?? ?? ?? ?01234567, , ,C m nC m nC m nC m nmnC m nC m nC m nC m n???????????????????????????? (335) 其中 ? ?1 2 3 4 5 6 7i i i ih g C p i ?、 、 、 、 、 、 、 、 、分別是 0 0 0 0h g C p、 、 、 的 i 階導(dǎo)數(shù)。同理可對磁流阻抗? ?0 ,D mn 進行求導(dǎo)。 23 計算實例與分析 在本文算例中，泰勒級數(shù)取 7階， Pad233。逼近取 (L=4， M=3)，矩量法計算步長 取 ,AWE計算步長 取 ,離散步長 取 10C ??? 。 算例 1:半 經(jīng) a=3mm 的 DNG 媒質(zhì)的無限長圓柱 ， 在 TM 平面波垂直照射下， DNG 媒質(zhì)參數(shù)分 別 取 0emHz? ?? ? ， 926 10 /em rad s??? ? ?，計算頻帶 為 18 22GH z GH z， 中心頻率分別 20GHz ， 解析解 、 矩量法和 AWE 計算的 RCS 對比如圖 32所示 。 算例 2： 半 經(jīng) a=3mm 的 DNG 媒質(zhì)的無限長圓柱 ， 在 TM 平面波垂直照射下， DNG 媒質(zhì)參數(shù)分 別 取 10em Hz? ? ? ? ?, 926 10 /em rad s??? ? ?，計算頻帶為 8 12GHz GHz，中心頻率為 10GHz ，解析解、矩量法和 AWE技術(shù) 計算的 RCS 對比如圖 33所示。 24 圖 32 雙負(fù)媒質(zhì)圓柱 RCS比較（ 0emHz? ?? ? ， 926 10 /em rad s??? ? ?） 25 圖 33雙負(fù)媒質(zhì)圓柱 RCS比較（ 10em Hz? ? ? ? ?， 926 10 /em rad s??? ? ?） 算例 3： 邊長 a=1mm 雙負(fù)媒質(zhì)二維方柱，在 TM 高斯脈沖平面波垂直照射 下，矩量法和 AWE技術(shù)計算的 RCS對比。 DNG媒質(zhì)參數(shù)分 別 取 0emHz? ?? ? ， 926 10 /em rad s??? ? ?和 10em Hz? ? ? ? ?， 926 10 /em rad s??? ? ?。 計算頻帶在 8 12GHz GHz， 中心頻率 為 10GHz。 解析解、矩量法和 AWE 技術(shù)計算的 RCS 對比如圖 34 和 35所示。 26 圖 34 雙負(fù)媒質(zhì)方柱 RCS比較（ 0emHz? ?? ? ， 926 10 /em rad s??? ? ?） 27 圖 35 雙負(fù)媒質(zhì)方柱 RCS比較（ 10em Hz? ? ? ? ?， 926 10 /em rad s??? ? ?） 從圖 32— 圖 35 的計算結(jié)果可以清晰看出： AWE 技術(shù)能很好逼近雙負(fù)媒質(zhì)的 MOM 逐點計算解； 泰勒級數(shù)展開和 Pad233。 逼近都在一定的頻帶內(nèi)和矩量法解 相吻合，但 Pad233。逼近比泰勒展開近似頻帶更寬，與理論分析相一致，達(dá)到了預(yù)期的效果。表 31 表明 AWE 能更有效計算雙負(fù)媒質(zhì)的寬頻 RCS 頻率響應(yīng)，計算效率比矩量法逐點計算大大提高。 表 31 MOM 和 AWE技術(shù)計算效率對照表 圖 32 圖 33 圖 34 圖 35 MOM 耗時 /s AWE 耗時 /s 提高倍數(shù) 28 4 總結(jié)與展望 寬頻帶的雷達(dá)電磁散射問題是計算電磁學(xué)領(lǐng)域一個極其重要而且快速發(fā)展 的課題。隨著科學(xué)技術(shù)的不斷發(fā)展，對復(fù)雜電磁環(huán)境下各種研究的展開，例如反隱身技術(shù)和隱身技術(shù)、寬帶雷達(dá)和超 寬帶雷達(dá)的技術(shù)的研究等， 都迫切需 對散射體目標(biāo)進行寬帶電磁散射特性以及電磁理論的研究與分析。在以往，分析寬頻帶的電磁散射特性都是利用傳統(tǒng)的矩量法進行逐個頻率點求解，在信息化的今天，特別是數(shù)據(jù)的快速傳輸和編碼，這種 傳統(tǒng)的做法是不可行的，極大的浪費 時間和計算機資源。因此，在計算機資源和 時間有限的情況下，急需實現(xiàn)更精確、 有效的寬頻計算方法。 作為一種新型的合成材料 — 雙負(fù)媒質(zhì)材料 ，分析其寬頻 電磁散射特性是及其重要 的。本文主要工作概括 為以下幾個方面 ： 1. 介紹了 雷達(dá)及 雷達(dá)散射截面的基本概念 ，闡述了雙負(fù)媒質(zhì)的研究背景及發(fā)展趨勢， 說明寬頻電磁計算的重要意義和研究價值。 2. 介紹了 PMCHWT 方程推導(dǎo)過程，作為最基本的理論 — 矩量法，本文 給予了重點的介紹， 詳細(xì)地闡述 了漸進波形估計技術(shù)的基本原理； 3. 研究了雙負(fù)媒質(zhì)的一般特性和基本模型 的建立， 針對雙負(fù) 媒質(zhì)的求導(dǎo)過程比較復(fù)雜，通過數(shù)學(xué)處理實現(xiàn)了雙負(fù)媒質(zhì)阻抗矩陣的求導(dǎo)， 并將 AWE 技術(shù)應(yīng)用于 雙負(fù)媒質(zhì)體散射的電磁散射特性分析 。 由于時間和自身等 原因， 本文只對形狀簡單的二維雙負(fù)媒質(zhì) 目標(biāo)進行了分析和方法的驗證，對于形狀復(fù)雜的介質(zhì)目標(biāo)還需進一步研究和探討，特別是對具有廣泛實際應(yīng)用的三維目標(biāo)同樣需要進行深入的研究?？偟膩碚f，今后尚需在以下幾方面開展更進一步的研究： 1. 有理函數(shù)逼近技術(shù)方面， 本文中所討論的僅是基于有理多項式分式函數(shù)思想的 Pad233。逼近，還可以考慮其余多項式的擬合來逼近， 如傅立葉葉級 數(shù)、切比雪夫多項式、賣 克勞林級數(shù)、勒讓德多項式等 ； 2. 散射體阻抗矩陣 求逆方面，本文 僅采用 MATLAB 語言求逆，其效果不太好， 所求出的阻抗矩陣都是特殊的矩陣，應(yīng)該 考慮一些特殊的預(yù)處理方法，提高阻抗矩陣求逆運算的效率。 3 可以進一步研究雙負(fù)媒質(zhì)其他方面特性及應(yīng)用。 29 參考文獻 [1] Veselago . The electrodynamics of substances with simultaneously negative values of ε and μ [J].Soviet PhysUspekhi, 1968, 10(4): 509514. [2] Pendry J B, Holden A J, Robbins D J, et al. Magism from Conductors and Enhanced Nonlinear Phenomena. IEEE Trans Microwave Theory Tech.,1999 ,47(11):20752084. [3] Shelby B A, Smith D R, Schultz S. Experimental verification of a negative index of refraction Science , 2001,292(6):7779. [4] Weiland T, Cchuhmann R, Greegor R B, et al. Ab initio numerical simulation of left _handed meta materials : Comparison of calculations and Appl Phys,2001,90:54195424. [5] Markos P, Soukoulis C M