陳傳升 王秉中 王任

(電子科技大學(xué)應(yīng)用物理研究所, 成都 611731)

隨著電磁器件的集成化, 器件搭載的模塊、實(shí)現(xiàn)的功能愈發(fā)多樣.各模塊間的耦合難以忽略, 設(shè)計(jì)難度陡然增加, 傳統(tǒng)設(shè)計(jì)方法逐漸力不從心, 迫切需要尋找一種新的電磁綜合設(shè)計(jì)方法.本文利用時(shí)間反演電磁波的時(shí)空同步聚焦特性, 探索了將時(shí)間反演技術(shù)應(yīng)用于器件設(shè)計(jì)的可能性.首先, 基于通用的器件逆設(shè)計(jì)流程, 利用時(shí)間反演技術(shù)、并矢格林函數(shù)及電磁學(xué)的基本原理, 提出了將器件端口場(chǎng)分布轉(zhuǎn)換為內(nèi)部場(chǎng)分布的方法, 并證明由端口期望場(chǎng)的時(shí)間反演場(chǎng)在空間某一位置獲得的連續(xù)等效源的共軛分布可在端口處產(chǎn)生與期望場(chǎng)接近的場(chǎng)分布.且在單點(diǎn)頻逆設(shè)計(jì)過程中, 只需知道端口電場(chǎng)或磁場(chǎng)的切向分量即可完成端口場(chǎng)與內(nèi)部場(chǎng)的轉(zhuǎn)換.同時(shí), 借助格林函數(shù)的互易性對(duì)本文所提理論做適當(dāng)變換后, 進(jìn)行數(shù)值仿真驗(yàn)證, 分析討論了不同初始信息條件下該方法的適用性.仿真結(jié)果與理論相符, 證明了理論的正確性, 為將時(shí)間反演技術(shù)應(yīng)用于電磁器件的逆設(shè)計(jì)提供了可能.

1 引 言

隨著器件向多功能、高集成度發(fā)展, 各功能模塊之間耦合數(shù)量的增多導(dǎo)致傳統(tǒng)模塊設(shè)計(jì)方法的設(shè)計(jì)復(fù)雜度、設(shè)計(jì)難度陡然增加.尋找一種非模塊化的新型設(shè)計(jì)方法勢(shì)在必行, 正因如此逆設(shè)計(jì)問題在諸多領(lǐng)域都受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[1?3].在納米光學(xué)領(lǐng)域中有學(xué)者利用逆向方法設(shè)計(jì)了單入雙出的分頻器[4], 該分頻器與應(yīng)用傳統(tǒng)半解析方法設(shè)計(jì)的器件相比, 在結(jié)構(gòu)上截然不同.通過將結(jié)構(gòu)的介電常數(shù)參數(shù)化至二維圖像, 結(jié)合拓?fù)鋬?yōu)化[5]通過梯度下降法求出局部最優(yōu)值, 再利用伴隨方法[6]完成設(shè)計(jì).該課題組利用相同的方法也實(shí)現(xiàn)了一分三的分頻器[7], 證實(shí)了該方法的可移植性[8].拓?fù)鋬?yōu)化方法在微波波段同樣得到了應(yīng)用并取得了不錯(cuò)的效果[9?11].隨著拓?fù)鋬?yōu)化這類不依賴以及其他部分依賴先驗(yàn)知識(shí)的設(shè)計(jì)方法的提出[12], 學(xué)者們?cè)谀嬖O(shè)計(jì)的道路上邁出了第一步.這些方法擺脫先驗(yàn)知識(shí)約束的同時(shí), 在優(yōu)化過程中卻不可避免地引入了大量迭代操作.為了提高效率, 學(xué)者們針對(duì)逆設(shè)計(jì)問題中應(yīng)用的優(yōu)化算法也做出了諸多改進(jìn)[13,14].誠(chéng)然, 算法效率的提高是降低迭代次數(shù)的有效的手段, 但只能實(shí)現(xiàn)量的減少, 若能脫離迭代的桎梏, 不借助迭代優(yōu)化而實(shí)現(xiàn)器件設(shè)計(jì), 那么設(shè)計(jì)效率就可以有質(zhì)的提高[2,15].

時(shí)間反演技術(shù), 是指對(duì)接收到的電磁波信號(hào)先進(jìn)行時(shí)間反轉(zhuǎn)處理, 再發(fā)射出去, 形成波矢方向反轉(zhuǎn)的時(shí)間反演電磁波.時(shí)間反演電磁波并非時(shí)間倒流[16], 但是卻伴有很多有趣的現(xiàn)象[17], 其中波的自適應(yīng)空時(shí)同步聚焦特性[18]更是被廣泛地應(yīng)用于無(wú)線通信、微波探測(cè)等諸多方面[19?22].本文正是基于這種聚焦特性, 根據(jù)逆設(shè)計(jì)的流程結(jié)合并矢格林函數(shù)提出了將器件端口場(chǎng)分布轉(zhuǎn)換為器件內(nèi)部區(qū)域場(chǎng)分布的理論方法, 并給出證明.借助該理論的互易性, 分別在自由空間和類波導(dǎo)腔體中完成了數(shù)值仿真驗(yàn)證, 仿真結(jié)果與理論結(jié)果較為一致.考慮到在器件設(shè)計(jì)過程中, 由已知參數(shù)信息往往不能得到完整的端口場(chǎng)分布, 又在已知部分端口場(chǎng)分布的條件完成了數(shù)值仿真, 證明了在某些條件下, 僅已知部分端口場(chǎng)分布并不影響該方法的應(yīng)用.

2 理論分析

電磁器件的逆設(shè)計(jì)的流程可概括為以下三步:第一步, 由目標(biāo)的端口參數(shù)出發(fā), 得到器件的端口場(chǎng)分布, 這一步是逆設(shè)計(jì)過程中初始條件的準(zhǔn)備工作.例如已知端口S參數(shù), 可以根據(jù)端口導(dǎo)波結(jié)構(gòu)(矩形波導(dǎo)、微帶線、或其他)的特性確定其定義S參數(shù)的工作模式的場(chǎng)的空間分布, 并由S參數(shù)導(dǎo)出端面場(chǎng)時(shí)域/頻域分布, 進(jìn)而推導(dǎo)出含有時(shí)空信息完整的端口場(chǎng)分布供逆設(shè)計(jì)過程使用.第二步:由器件的端口場(chǎng)分布推導(dǎo)、計(jì)算出器件內(nèi)部局域場(chǎng)分布, 這一步涉及逆設(shè)計(jì)的實(shí)質(zhì)是逆設(shè)計(jì)的核心步驟, 亦是本文的研究所在.第三步: 將器件內(nèi)部局域場(chǎng)分布轉(zhuǎn)換為器件內(nèi)部對(duì)應(yīng)的結(jié)構(gòu)分布完成設(shè)計(jì), 探尋場(chǎng)與結(jié)構(gòu)的關(guān)系, 這一步也是逆設(shè)計(jì)方法未來(lái)研究工作的重難點(diǎn)之一, 目前僅有部分文章對(duì)個(gè)別特例提出了場(chǎng)到結(jié)構(gòu)的轉(zhuǎn)換方法, 且方法存在一定的局限性.例如, 文獻(xiàn)[3]在不考慮磁電耦合的情況下, 由源分布計(jì)算得到對(duì)應(yīng)的電磁介質(zhì)參數(shù)分布, 并利用FDTD-GTSC 求解器[23]和3D EFIE求解器分別在二維、三維模型上完成了正向驗(yàn)證,結(jié)果符合預(yù)期.在文末卻也指出, 這套方法雖然成功地將源分布轉(zhuǎn)化為電磁介質(zhì)參數(shù)分布, 但如何在物理上實(shí)現(xiàn)這種電磁介質(zhì)參數(shù)分布仍是亟待解決的問題.又如, 文獻(xiàn)[24]通過求解金屬小球陣列的等效偶極子矢量間的耦合方程組也定量地建立了場(chǎng)與結(jié)構(gòu)的關(guān)聯(lián).

本文主要針對(duì)上述流程的第二步提出解決方法, 即如何完成由端口場(chǎng)到內(nèi)部場(chǎng)之間的轉(zhuǎn)換.本文提出的理論表明在利用端口參數(shù)得到端口場(chǎng)分布后, 對(duì)端口場(chǎng)分布施加時(shí)間反演操作獲得端口反演源, 端口反演源在器件內(nèi)部形成一定的場(chǎng)分布,只要能實(shí)現(xiàn)該場(chǎng)分布的時(shí)間反演分布, 便可在端口處獲得期望的端口參數(shù).本節(jié)立足于等效原理和并矢格林函數(shù)從數(shù)理模型和理論證明兩個(gè)方面闡述該理論.本文所提方法與拓?fù)鋬?yōu)化等逆設(shè)計(jì)方法相比, 最大的優(yōu)勢(shì)在于避免了冗長(zhǎng)的迭代優(yōu)化過程.有必要說明, 點(diǎn)頻情況下的時(shí)間反演操作等價(jià)于頻域的共軛操作, 最終解為穩(wěn)態(tài)場(chǎng), 因此本方法可適用于器件逆設(shè)計(jì)中端口場(chǎng)與內(nèi)部場(chǎng)的轉(zhuǎn)換.

2.1 數(shù)理模型

考慮如圖1 所示的空間拓?fù)潢P(guān)系, 與逆設(shè)計(jì)的流程相對(duì)應(yīng).表面a作為設(shè)計(jì)域, 代表著器件內(nèi)部某一部分的場(chǎng)分布, 對(duì)應(yīng)器件內(nèi)部的具體結(jié)構(gòu), 是待求解量; 表面b作為目標(biāo)域, 代表著設(shè)計(jì)過程中器件的端口場(chǎng)分布, 在器件逆設(shè)計(jì)過程中這是是已知的, 其中上角標(biāo)的正負(fù)符號(hào)表示緊貼該表面的一層虛擬平面.

圖1 空間拓?fù)潢P(guān)系Fig.1.Topological relation of space.

若在表面a上存在著連續(xù)的面電流、磁流源Ja,Ma, 根據(jù)電磁場(chǎng)格林函數(shù)和電磁場(chǎng)的疊加原理, 可將Ja,Ma在表面b上形成的場(chǎng)可以表示為(1)式, 簡(jiǎn)寫為(2)式:

其中GPQ(r,r′) 為并矢格林函數(shù), 表示位于r′處的Q源與位于r處的P場(chǎng)之間的聯(lián)系;〈〉表示兩個(gè)函數(shù)在公共區(qū)域上的乘積, 用“, ”隔開, 寫作“;”表示乘積的形式為內(nèi)積.

器件設(shè)計(jì)的最終目的是在設(shè)計(jì)域a中合理地排布結(jié)構(gòu), 使其等效電流源、磁流源在表面b上形成的場(chǎng)符合特定需求.顯然, 這是一個(gè)逆問題, 即尋找一組源分布使其在b面處產(chǎn)生的場(chǎng)最大程度上與已知的期望場(chǎng)相同; 即求解使得代價(jià)函數(shù)式 最 小,ρa(bǔ),b表示a,b兩組數(shù)據(jù)的相關(guān)系數(shù).

2.2 理論證明

本小節(jié)借助并矢格林函數(shù)和等效原理, 理論證明了利用端口期望場(chǎng)的時(shí)間反演場(chǎng)可在空間某一位置獲得的連續(xù)等效源, 該連續(xù)等效源可在端口處產(chǎn)生與期望場(chǎng)接近的場(chǎng)分布的結(jié)論, 一定程度上說明了將時(shí)間反演理論用于器件設(shè)計(jì)中是可行的.

在目標(biāo)域b待實(shí)現(xiàn)特定場(chǎng)分布是已知的, 其時(shí)間反演場(chǎng)在a–處產(chǎn)生場(chǎng)借助等效原理將其看作由位于表面b處的等效源在a–處產(chǎn)生的場(chǎng), 故寫作(6)式, 其時(shí)間反演場(chǎng)寫作(8)式,其中由等效原理寫作(7)式.(8)式中,已知, 故可 求,再利用等效原理得(9)式即為所求.

(9)式的意義在于通過引入格林函數(shù), 利用取切向分量和頻域共軛的操作, 由已知場(chǎng)可求解出一組位于表面a處源分布且該源分布的求解只與已知場(chǎng)和a,b平面之間的格林函數(shù)有關(guān).不難證明,這組位于a處的源分布, 在b處產(chǎn)生的場(chǎng)與期望場(chǎng)一致.理想情況下,代價(jià)函數(shù)為0, 代價(jià)函數(shù)為0 的物理意義是由端口期望場(chǎng)的時(shí)間反演場(chǎng)可在器件內(nèi)部形成的場(chǎng)構(gòu)建出的連續(xù)等效源, 該連續(xù)等效源的時(shí)間反演場(chǎng)在端口產(chǎn)生的場(chǎng)與端口期望場(chǎng)是一致的.由對(duì)偶原理式(10)易證, 其中c1,c2為常數(shù), 代價(jià)函數(shù)的理論值為0.

3 數(shù)值仿真流程

圖2 用于數(shù)值驗(yàn)證的問題轉(zhuǎn)化示意圖 (a)數(shù)理模型對(duì)應(yīng)的數(shù)值驗(yàn)證模型; (b)利用互易原理轉(zhuǎn)換后的數(shù)值驗(yàn)證模型Fig.2.A problem transformation diagram for numerical validation: (a) Numerical verification model corresponding to mathematical model; (b) numerical verification model converted by reciprocity principle.

根據(jù)格林函數(shù)的互易性, 將上一節(jié)的理論推導(dǎo)中涉及的上角標(biāo)a,b對(duì)換推導(dǎo)仍然成立.事實(shí)上,上角標(biāo)a,b的調(diào)換, 可視作已知源和未知源的調(diào)換.如圖2 所示, 第二節(jié)中理論對(duì)應(yīng)的數(shù)值驗(yàn)證模型為圖2(a), 借助互易原理轉(zhuǎn)換為圖2(b)所示的情形, 即原有的已知源視作未知源, 原有的未知源視作已知源.數(shù)值驗(yàn)證的目的是驗(yàn)證本文所提出理論的正確性, 為了得到更為直觀的觀察效果, 這里選用圖2(b)所示的情況進(jìn)行驗(yàn)證.即對(duì)本文所提理論的數(shù)值驗(yàn)證問題轉(zhuǎn)化為當(dāng)表面a處為已知源,驗(yàn)證up, down 兩個(gè)表面的反演源在表面a處生成的重建源與已知的表面a處的初始源的一致性.根據(jù)上一節(jié)的理論推導(dǎo), up, down 兩個(gè)表面的反演源經(jīng)激勵(lì)會(huì)在表面a處生成與初始源趨勢(shì)一致的重建源, 兩者的模值存在一定的比例關(guān)系, 代價(jià)函數(shù)應(yīng)趨近0.

數(shù)值仿真具體可分為以下三個(gè)步驟:

1) 表面a上已知初始源的獲得;

2) 第一次仿真, 初始源激勵(lì), 由up, down 雙平面記錄場(chǎng)并獲得反演源;

3) 第二次仿真, 由不同反演源激勵(lì), 比較表面a處的重建源與已知初始源的異同, 并計(jì)算其代價(jià)函數(shù).

借助達(dá)索公司的商用電磁仿真軟件CST 微波仿真實(shí)驗(yàn)室和python 完成數(shù)值仿真及后續(xù)的數(shù)據(jù)分析處理.數(shù)值仿真實(shí)驗(yàn)的示意圖所圖3 所示, 對(duì)于圖3(a), 自由空間中存在著up, down 和表面a三個(gè)有限大平面, 其平面方程為(11)式, D 為該平面方程對(duì)應(yīng)的定義域,λ為真空中10 GHz 電磁波所對(duì)應(yīng)的波長(zhǎng); 圖3(b)與圖3(a)尺寸相同, 唯一的區(qū)別是圖(b)四周是與up, down 兩個(gè)平面嚙合的理想電導(dǎo)體(perfect electric conductor, PEC)邊界.各平面均設(shè)有場(chǎng)監(jiān)視器用以記錄場(chǎng)分布, 各平面的場(chǎng)分布記為其中field ∈{Ex,Ey,Ez,Hx,Hy,Hz},surface ∈{a,up,down}.

定義一種變換, 記作Tm,n(x) , 該操作表示將連續(xù)量x(x∈R2)離散采樣為m×n的矩陣, 通過這一變換, 將理論推導(dǎo)中的連續(xù)源轉(zhuǎn)換為具體的矩陣, 象征離散源, 進(jìn)而方便后續(xù)數(shù)值仿真處理與計(jì)算, 例如連續(xù)的面電流源離散采樣如(12)式所示.

圖3 數(shù)值仿真實(shí)驗(yàn)示意圖 (a)自由空間驗(yàn)證示意圖;(b)四周為理想電導(dǎo)體邊界的驗(yàn)證示意圖Fig.3.Schematic diagram of numerical simulation experiment: (a) Schematic diagram of free space verification;(b) schematic diagram of verification of PEC boundaries around.

3.1 數(shù)值仿真初始源的獲得

的定義域?yàn)?Da, 而p(x,y) 的定義域?yàn)?Dp,又因?yàn)?Dp→Da之間為單射, 則有K(p(x,y))→或K表示某種映射關(guān)系.在數(shù)值仿真中, 設(shè)置了頻率均為10 GHz 點(diǎn)頻的兩種初始源, 第一種源, 其電場(chǎng)幅度等于1, 電場(chǎng)相位分布為peaks函數(shù); 第二種源的電場(chǎng)相位同相分布, 幅度分布為peaks函數(shù),兩種源分別記為簡(jiǎn)記為如(15)式所示.需要說明的是, 在仿真過程中, 采用作為第一次數(shù)值仿真中使用的激勵(lì)源.

3.2 第一次仿真, 獲得端口場(chǎng)反演源

3.3 不同反演源激勵(lì)情況下, 重建源與初始源的異同

為了量化重建源與初始源的異同, 選取合適的代價(jià)函數(shù)是很有必要的.

選取相關(guān)距離作為代價(jià)函數(shù)的原因是, 由(10)式可知重建場(chǎng)與初始場(chǎng)只存在模值上的縮小, 其變化趨勢(shì)應(yīng)當(dāng)一致, 而相關(guān)距離恰好是研究變量之間線性相關(guān)程度的量.因此, 將(12)式代入(16)式得(17)式, 系離散情況下衡量重建場(chǎng)與期望場(chǎng)之間的代價(jià)函數(shù), 用以衡量重建場(chǎng)與原場(chǎng)的相似程度.在自由空間、四周為理想電導(dǎo)體邊界兩種條件下對(duì)兩種不同的初始源, 采用多種不同的信息進(jìn)行重建, 并計(jì)算其代價(jià)函數(shù).源的選用及計(jì)算結(jié)果總結(jié)在表1 中, 縮略詞F P/A CF (free space phase/amplitude rebuild cost function), 意為在圖3(a)所示的自由空間條件下, 重建源與初始源或代價(jià)函數(shù)的計(jì)算結(jié)果; P P/A CF(space surrounded by PEC phase/amplitude rebuild cost function)則表示在圖3(b) 所示的四周為理想電導(dǎo)體邊界條件下, 重建源和初始源的代價(jià)函數(shù)計(jì)算結(jié)果.在3.2 節(jié)中雖然獲得了反演場(chǎng)的所有場(chǎng)信息, 但在重建過程中卻并非要使用所有的反演源信息, 既可以利用兩個(gè)反演源的部分信息,也可以選用其中一個(gè)反演源的全部信息, 甚至只利用部分反演源的部分信息.對(duì)和使用重建時(shí), 反演源的選用有以下幾種可能, 其代價(jià)函數(shù)計(jì)算結(jié)果總結(jié)在表1 中, 代價(jià)函數(shù)曲線繪制在圖4 中.

表1 數(shù)值仿真結(jié)果表Table 1.Table of numerical simulation results.

圖4 重建源代價(jià)函數(shù)計(jì)算結(jié)果Fig.4.The calculation results of the cost function of the rebuild source.

3.4 重建結(jié)果及分析

觀察表1 及圖4 不難發(fā)現(xiàn), 不論使用何種信息, 最終獲得的重建源與初始源都十分接近的, 且在相同空間條件下, 對(duì)相位源的重建結(jié)果總好于幅度源.使用兩個(gè)平面的全部場(chǎng)信息能獲得較好的重建效果, 與本文的理論分析一致; 使用兩個(gè)平面的部分場(chǎng)信息仍能實(shí)現(xiàn)較好的重建效果, 是因?yàn)楦鶕?jù)穩(wěn)態(tài)電磁場(chǎng)中的唯一性原理“在閉合面S包圍的區(qū)域V中, 當(dāng)邊界面S上的電場(chǎng)強(qiáng)度的切向分量E或者磁場(chǎng)強(qiáng)度H給定時(shí), 體積V中任一點(diǎn)的電磁場(chǎng)由Maxwell 方程唯一的確定.”序號(hào)為3,6, 9 的仿真實(shí)驗(yàn)的成功重建說明了這一點(diǎn).對(duì)于圖3(b)所示的條件, 根據(jù)唯一性定理只需要知道的切向分量即可唯一地確定surfacea的場(chǎng)分布, 因此field={Ex,Ey}.對(duì)于圖3(a)所示的空間示意圖, 其四周為開放空間, 根據(jù)唯一性原理若僅知道端口平面up, down 的切向分量在理論上是不能完整重建初始源的, 但數(shù)值仿真實(shí)驗(yàn)仍然實(shí)現(xiàn)了較好的重建效果, 分析原因可能是因?yàn)橹亟▓?chǎng)的極化較為單一.

考慮仿真實(shí)驗(yàn)的結(jié)果較多, 因而只選取部分結(jié)果作展示.對(duì)不同的仿真實(shí)驗(yàn)情況, 在實(shí)驗(yàn)序號(hào)為1—6 的實(shí)驗(yàn)選取一個(gè)結(jié)果, 在實(shí)驗(yàn)序號(hào)為7—9 的實(shí)驗(yàn)選取一個(gè)結(jié)果, 選取原則為選取重建結(jié)果代價(jià)函數(shù)最大的仿真實(shí)驗(yàn)代號(hào)作為展示.代價(jià)函數(shù)越大代表其復(fù)原效果越差, 但即使是最差的復(fù)原效果仍得到了較完好的重建效果, 充分證明了理論的正確性.

圖6 展示了在圖3(a)所示自由空間條件下對(duì)圖5 所示初始源選用不同反演信息激勵(lì)時(shí)得到的重建源幅相分布,其重建結(jié)果與圖5 十分接近.圖6(a)和圖6(b)是對(duì)的重建結(jié)果展示圖, 對(duì)應(yīng)仿真實(shí)驗(yàn)序號(hào)為5, 8, 代表了仿真實(shí)驗(yàn)中獲得的最差重建結(jié)果, 其中主圖為重建源的相位分布, 左上角為幅度分布.圖6(c)和圖6(d)是對(duì)的重建結(jié)果, 對(duì)應(yīng)仿真實(shí)驗(yàn)序號(hào)為2, 8, 其中主圖為重建源的幅度分布, 左上角為相位分布.圖6(d)的重建幅度約為圖6(c)的兩倍是因?yàn)槭褂脙蓚€(gè)反演源激勵(lì)時(shí)其輸入功率是單個(gè)反演源激勵(lì)的2 倍, 除幅度差異外, 圖6(c)和圖6(d)兩圖的趨勢(shì)一致且其重建的相位分布均為同向分布.不難發(fā)現(xiàn), 不論是選何種反演源信息重建何種源, 最終獲得的重建源與初始源的幅相分布都是十分接近的, 其代價(jià)函數(shù)計(jì)算結(jié)果均小于0.1, 代表著重建源與初始源十分接近.

圖5 初始源幅相分布圖 (a) E x 相位分布為特定函數(shù); (b) E x 幅度分布為特定函數(shù)Fig.5.The amplitude-phase distribution of the original source: (a) The phase distribution is a special function; (b) the amplitude distribution is a special function.

圖6 自由空間條件下重建源的幅相分布 (a)序號(hào)為5 的重建的幅相分布; (b) 序號(hào)為8 的重建 的幅相分布; (c)序號(hào)為2 的重建 的幅相分布; (d)序號(hào)為8 的重建 的幅相分布Fig.6.The amplitude-phase distribution of the reconstructed source in free space: (a) The amplitude-phase distribution of the reconstructed with experimental number 5; (b) the amplitude-phase distribution of the reconstructed with experimental number 8; (c) the amplitude-phase distribution of the reconstructed with experimental number 2; (d) the amplitude-phase distribution of the reconstructed with experimental number 8.

圖7 四周為理想電導(dǎo)體條件下重建源的幅相分布 (a)序號(hào)為4 的重建 的幅相分布; (b) 序號(hào)為7的重建的幅相分布; (c)序號(hào)為3 的重建 的幅相分布; (d)序號(hào)為9 的重建 的幅相分布Fig.7.The amplitude-phase distribution of the reconstructed source in PEC space: (a) The amplitude-phase distribution of the reconstructed with experimental number 4; (b) the amplitude-phase distribution of the reconstructed with experimental number 7; (c) the amplitude-phase distribution of the reconstructed with experimental number 3; (d) the amplitude-phase distribution of the reconstructed with experimental number 9.

圖7 與圖6 的不同之處僅在于其所處的仿真實(shí)驗(yàn)環(huán)境為圖5(b)所示的四周為理想電導(dǎo)體條件的環(huán)境, 重建結(jié)果展示的選取原則與圖6 一致.圖7(a)和圖7(b)是對(duì)的重建結(jié)果展示圖, 對(duì)應(yīng)仿真實(shí)驗(yàn)序號(hào)4, 7; 圖7(c)和圖7(d)是對(duì)的重建結(jié)果, 對(duì)應(yīng)仿真實(shí)驗(yàn)序號(hào)為3, 9.與圖5 所示的初始源對(duì)比不難發(fā)現(xiàn), 圖7重建結(jié)果雖然與圖5 較為相似, 但多出了許多規(guī)律性斑點(diǎn), 且斑點(diǎn)間距接近一個(gè)波長(zhǎng).分析原因可能是反演源的激勵(lì)不僅形成了重建源, 同時(shí)也激勵(lì)出腔體自身的本征模式.事實(shí)上, 仿真實(shí)驗(yàn)序號(hào)為1, 4, 7 的重建源斑點(diǎn)分布較其他仿真實(shí)驗(yàn)代號(hào)更為明顯, 分析原因可能是數(shù)值仿真中存在一定的仿真誤差, 而磁場(chǎng)反演源的誤差引入恰好激勵(lì)出了更強(qiáng)的腔體本征模式.在圖7 中同樣可以觀察到圖7(d)圖的重建幅度約為圖7(c)的兩倍, 與圖6(c)和圖6(d)兩圖存在差異的原因相同.

4 結(jié) 論

本文借助時(shí)間反演技術(shù)提出了將器件端口場(chǎng)分布轉(zhuǎn)換為器件內(nèi)部區(qū)域場(chǎng)分布的理論方法.利用理論的互易性實(shí)施了系列數(shù)值仿真實(shí)驗(yàn), 仿真實(shí)驗(yàn)結(jié)果與理論相符合, 說明了本文提出理論的正確性.即利用端口期望場(chǎng)的時(shí)間反演場(chǎng)可在空間某一位置獲得連續(xù)等效源, 該連續(xù)等效源的共軛分布可在端口處產(chǎn)生與期望場(chǎng)接近的場(chǎng)分布.從唯一性原理出發(fā), 在點(diǎn)頻逆設(shè)計(jì)由端口參數(shù)得到端口場(chǎng)分布的過程中, 不需獲得全部的端口場(chǎng)分布, 只需知道端口的電場(chǎng)或磁場(chǎng)的切向分量即可.此外, 本方法應(yīng)用于四周為理想電導(dǎo)體條件的腔體時(shí)應(yīng)當(dāng)充分考慮腔體本征模式的影響, 如何解決腔體本征模式帶來(lái)的影響, 也是后續(xù)亟待解決的問題之一.