孔明昊 肖 龍 吳 雨 陳 亮

（1.中國(guó)艦船研究設(shè)計(jì)中心 武漢 430064）（2.電磁兼容性重點(diǎn)實(shí)驗(yàn)室 武漢 430064）

1 引言

近年來(lái)，微波無(wú)線傳能在生物醫(yī)學(xué)、太陽(yáng)能充電、衛(wèi)星導(dǎo)航等領(lǐng)域發(fā)揮了重要作用［1］。但是僅僅是傳輸能量還遠(yuǎn)遠(yuǎn)不夠，通過(guò)對(duì)電磁場(chǎng)的賦形，使得微波的能量聚集在我們所預(yù)期的位置，控制不同位置的場(chǎng)強(qiáng)也對(duì)提高我們利用電磁波能量的效率，實(shí)現(xiàn)更多的應(yīng)用有重要的意義。電磁場(chǎng)賦形隨著微波熱療、近場(chǎng)識(shí)別、醫(yī)療檢測(cè)等應(yīng)用的發(fā)展正在受到越來(lái)越多的關(guān)注。2008 年，Landy 等利用吸波材料實(shí)現(xiàn)了對(duì)入射電磁波幅度的調(diào)控［2］。2019年，J.W.W 等提出了一種綜合算法并應(yīng)用人工電磁材料實(shí)現(xiàn)了電磁場(chǎng)的賦形，得到了較好的賦形效果［3］。傳統(tǒng)的電磁場(chǎng)賦形集中在利用電磁超材料進(jìn)行傳播控制上，使用超材料的缺陷靈活性低，而且?guī)捿^窄和損耗較大。

時(shí)間反演（Time Reversal，TR）技術(shù)具有自適應(yīng)空時(shí)聚焦特性和環(huán)境自適應(yīng)性［4］，能實(shí)現(xiàn)任意位置的場(chǎng)聚焦，這種特性為復(fù)雜環(huán)境下任意形狀電磁場(chǎng)賦形提供了可能。利用這些特性有助于解決很多實(shí)際問(wèn)題，比如用于水下超聲探測(cè)和通信［5］、微波探測(cè)［6］、微波腫瘤治療［7］、微波能量的傳輸［8］等領(lǐng)域。2005 年，Devaney 使用從稀疏和非結(jié)構(gòu)化相控陣天線系統(tǒng)收集的多基地?cái)?shù)據(jù)，將時(shí)間反演方法應(yīng)用于雷達(dá)成像問(wèn)題［9］。2020年，何紫靜等通過(guò)仿真研究了相關(guān)系數(shù)的正交性與聚焦點(diǎn)距離和時(shí)間反演鏡數(shù)量等參數(shù)之間的關(guān)系，有助于提高時(shí)間反演聚焦場(chǎng)的質(zhì)量［10］。2021 年，丁帥等提出了一種基于時(shí)間反演理論的電磁（EM）元透鏡貝塞爾波束的合成和實(shí)現(xiàn)方法，可以在應(yīng)用微波通信領(lǐng)域，拓寬了非衍射光束的應(yīng)用范圍［11］。時(shí)間反演技術(shù)在接收電磁波的過(guò)程中，其回傳相位信息已經(jīng)包含了周圍環(huán)境的影響因素，因此回傳信號(hào)在多徑效應(yīng)明顯的環(huán)境中也能獲得良好的聚焦效果，電磁場(chǎng)能量可以在賦形場(chǎng)源處聚焦，形成一個(gè)預(yù)期形狀的電磁場(chǎng)分布，這樣賦形系統(tǒng)可以省去光電或者數(shù)據(jù)鏈等額外的對(duì)準(zhǔn)設(shè)備。

針對(duì)本文討論了幾種時(shí)間反演實(shí)現(xiàn)方法并進(jìn)行了仿真分析，提出了一種面向電磁場(chǎng)賦形的時(shí)間反演異構(gòu)設(shè)計(jì)。基于此設(shè)計(jì)，構(gòu)建了墻壁環(huán)境的信道模型，并完成了特定圖案的電磁場(chǎng)賦形仿真。

2 時(shí)間反演技術(shù)

時(shí)間反演技術(shù)，即對(duì)接收到的電磁波信號(hào)先進(jìn)行時(shí)間反轉(zhuǎn)處理再發(fā)射出去，形成波矢方向反轉(zhuǎn)的時(shí)間反演電磁波。1959 年，Wigner 指出時(shí)間反演不是時(shí)間倒流，而是運(yùn)動(dòng)方向的倒轉(zhuǎn)，正反運(yùn)動(dòng)過(guò)程都必須遵循相同的因果性［12］。時(shí)間反演電磁波在傳播介質(zhì)中的存在和傳輸，這是時(shí)間反演應(yīng)用的前提，通過(guò)波動(dòng)方程可以推導(dǎo)均勻介質(zhì)和非均勻媒質(zhì)中的波動(dòng)方程同樣具備時(shí)間對(duì)稱性，而通過(guò)格林函數(shù)可以推導(dǎo)出三類邊界條件均有電磁場(chǎng)互易定理。對(duì)于均勻平面波在均勻介質(zhì)和非均勻媒質(zhì)中傳播，其回傳波滿足條件，就能實(shí)現(xiàn)電磁波的時(shí)間反演。

時(shí)間反演技術(shù)具有自適應(yīng)空時(shí)聚焦特性和環(huán)境自適應(yīng)性，利用時(shí)間反演的空間、時(shí)間雙重聚焦特性進(jìn)行電磁場(chǎng)賦形，可以自適應(yīng)匹配信道，有效抑制甚至利用多徑干涉效應(yīng)進(jìn)行傳播，增強(qiáng)聚焦效果，利用時(shí)間反演對(duì)單點(diǎn)和多點(diǎn)的聚焦效果仿真如圖1所示。

圖1 基于時(shí)間反演聚焦的焦平面能量分布圖

3 時(shí)間反演實(shí)現(xiàn)方法

3.1 啁啾變換原理

啁啾變換方法，即線性調(diào)頻。該技術(shù)發(fā)源于聲學(xué)領(lǐng)域，因?yàn)榫€性調(diào)頻信號(hào)類似于鳥(niǎo)聲，也叫做Chirp 信號(hào)，它的核心是對(duì)輸入信號(hào)在時(shí)域上進(jìn)行啁啾傅里葉變換（Chirp Fourier Transform，CFT）操作，簡(jiǎn)稱Chirp 變換。設(shè)r(t)為輸入信號(hào)，其Chirp變換原理如圖2所示。

圖2 Chirp變換原理框圖

圖2中，μ1掃頻本振信號(hào)的啁啾率，μ2為濾波器的啁啾率。色散濾波器的零狀態(tài)響應(yīng)為ejμ2t2，則系統(tǒng)輸出：

其中，Chirp變換結(jié)果y(t)為

當(dāng)μ1=-2μ2時(shí)，信號(hào)在時(shí)間序列上反轉(zhuǎn)，完整復(fù)原了r(t)，即完成了時(shí)間反演。對(duì)該過(guò)程可以進(jìn)行仿真分析，假設(shè)輸入一個(gè)三角波，時(shí)寬6ns，r=t；當(dāng)掃頻信號(hào)啁啾率μ1=100GHz/s；濾波器信號(hào)啁啾率μ2分別取μ1，-2μ1，可以得到如圖3 所示時(shí)間拉伸、時(shí)間壓縮的結(jié)果。圖4 展示了當(dāng)μ2=-0.5μ1時(shí)，信號(hào)載波、混頻后、時(shí)間反演后、解調(diào)后的數(shù)值仿真結(jié)果。

圖3 基于Chirp變換的時(shí)間反演拉伸壓縮數(shù)值仿真圖

圖4 基于Chirp變換的時(shí)間反演調(diào)制解調(diào)器數(shù)值仿真圖

基于啁啾變換方法實(shí)現(xiàn)時(shí)間反演的需要設(shè)計(jì)特定啁啾率的啁啾色散濾波器，且要滿足μ1=-2μ2的條件，對(duì)混頻器的輸入頻率范圍要求很高，適合于窄帶輸入信號(hào)。

3.2 時(shí)間透鏡原理

時(shí)域上的相位調(diào)制過(guò)程與物理中的薄透鏡十分相似，通過(guò)色散和二次相位調(diào)制，可以完成成像系統(tǒng)的時(shí)域模擬，信號(hào)會(huì)發(fā)生拉伸、壓縮和反演這三種變換［13］。當(dāng)時(shí)間透鏡在特定條件下實(shí)現(xiàn)反演的作用時(shí)，這種透鏡即為T(mén)RM。時(shí)間透鏡方法是利用模擬電路技術(shù)實(shí)現(xiàn)電磁信號(hào)時(shí)間反演的經(jīng)典方法［14］。時(shí)間透鏡方法實(shí)現(xiàn)時(shí)間反演鏡與Chirp變換方法的時(shí)間反演系統(tǒng)有類似之處［15］，不同之處在于它多了一段啁啾色散延遲線。

其工作原理并不復(fù)雜，射頻信號(hào)先經(jīng)過(guò)一個(gè)啁啾率為μ1的色散器件，再與到啁啾率為μ2的掃頻本振混頻，然后再經(jīng)過(guò)一個(gè)啁啾率μ3為的色散器件，最后解調(diào)出包絡(luò)信號(hào)［16］。假設(shè)輸入信號(hào)為r(t)，則經(jīng)過(guò)圖5所示的系統(tǒng)后，解調(diào)之前的輸出：

圖5 基于時(shí)間透鏡原理的時(shí)間反演系統(tǒng)原理框圖

1989 年，美國(guó)Hewlett-Packard 實(shí)驗(yàn)室的Brian H Kolner 推導(dǎo)出來(lái)的窄帶脈沖通過(guò)時(shí)間透鏡后時(shí)域成像條件［13］，和物理透鏡的空間成像條件具有驚人的相似。推廣到放大特性上，物理透鏡放大系數(shù)M為透鏡到虛像距離與實(shí)體距離的比(-di/do)。時(shí)間透鏡的放大特性與輸出色散與輸入色散之比M有關(guān)，即

若要重建信號(hào)，還須滿足放大倍數(shù)M 為-1，即使用相同的色散延遲線。滿足μ3=μ1,μ2=-2μ1時(shí)，脈沖通過(guò)時(shí)間透鏡后完成反演，成為T(mén)RM，如式（5）。

對(duì)式（5）進(jìn)行解調(diào)，最終解調(diào)出的信號(hào)即為輸入信號(hào)的時(shí)間反演信號(hào)。對(duì)該過(guò)程可以進(jìn)行仿真分析，假設(shè)r(t) 為輸入中心頻率為1GHz，帶寬200M0Hz 的高斯信號(hào)，時(shí)寬10ns。取色散延遲線1的啁啾率μ1=100GHz/s，掃頻本振信號(hào)的啁啾率μ2=-20GHz/s，色散延遲線2 的啁啾率μ3=10GHz/s。時(shí)間反演鏡反演過(guò)程數(shù)值仿真圖如圖6 所示，數(shù)值仿真結(jié)果如圖7所示。

圖6 時(shí)間反演鏡反演過(guò)程數(shù)值仿真圖

圖7 時(shí)間反演鏡反演結(jié)果數(shù)值仿真圖

時(shí)間透鏡方法相比于Chirp 變換解調(diào)效果更佳，但是依舊需要設(shè)計(jì)特定啁啾率的啁啾色散延遲線，且時(shí)寬和帶寬要滿足的條件較為嚴(yán)苛，對(duì)混頻器的輸入頻率范圍要求很高，適合于窄帶輸入信號(hào)。

3.3 方向回溯方法

從時(shí)域上看，時(shí)間反演處理是指天線接收到一串電磁信號(hào)f(t)后，將其在時(shí)域上翻轉(zhuǎn)為f(-t)后再次發(fā)射出去，這種過(guò)程接近于“方向回溯”。1964年，Y.C.PON 提出了基于外差混頻的方向回溯陣列［17］。按照實(shí)現(xiàn)方式，可以分為射頻混頻和中頻混頻，如圖8所示。

圖8 外差混頻陣列

射頻混頻是采用一個(gè)頻率為射頻（RF）兩倍的本振信號(hào)（LO）和接收到的射頻信號(hào)混頻，三次諧波分量被低通濾波器濾掉，從而得到的中頻IF 頻率和RF頻率相等，相位共軛［18］，具體理論推導(dǎo)如式所示。采用這種混頻方式可以以較為簡(jiǎn)單的電路實(shí)現(xiàn)方向回溯，但由于RF 和IF 相等，因此實(shí)際使用中，要考慮RF泄露的問(wèn)題。

在上述的推導(dǎo)中只是對(duì)于單路天線單元而言。在實(shí)際電路實(shí)現(xiàn)中，共軛電路總會(huì)對(duì)信號(hào)產(chǎn)生額外相移。由于假定每一路天線以及電路都是完全相同的，所以這一部分額外相移對(duì)于每一個(gè)天線單元都是相同的，不會(huì)影響波束指向。實(shí)際上每一路天線單元和電路總會(huì)有或多或少的差異，這一部分差異會(huì)使每個(gè)天線單元之間有額外的相位差，從而影響方向圖。相對(duì)于VanAtta 陣列，外差混頻方向回溯陣采用每個(gè)單元內(nèi)部獨(dú)立實(shí)現(xiàn)相位共軛，入射波可以不用限制為平面波，其應(yīng)用范圍較為廣泛［19］。

圖9 是導(dǎo)航信號(hào)進(jìn)行外差混頻的仿真圖，與式結(jié)果一致，實(shí)際混頻器會(huì)內(nèi)置低通濾波，該信號(hào)再通過(guò)混頻器就可以得到時(shí)間反演信號(hào)。該方式適用于單頻點(diǎn)信號(hào)，要考慮混頻器的RF 泄露的問(wèn)題。

圖9 導(dǎo)航信號(hào)進(jìn)行外差混頻的仿真圖

3.4 數(shù)字信號(hào)處理方法

時(shí)間反演的實(shí)現(xiàn)過(guò)程中可以理解為時(shí)域上翻轉(zhuǎn)操作，或者頻域上共軛處理再回發(fā)，即可實(shí)現(xiàn)了時(shí)間反演時(shí)空聚焦的效果。數(shù)字處理方法實(shí)現(xiàn)時(shí)間反演處理就是基于此推論。首先通過(guò)數(shù)字取樣將模擬信號(hào)轉(zhuǎn)變?yōu)閿?shù)字信號(hào)，即可對(duì)數(shù)字信號(hào)進(jìn)行傅里葉變換，做完共軛處理之后進(jìn)行及逆變換，最終經(jīng)過(guò)DA轉(zhuǎn)換，即可獲得時(shí)域的時(shí)間反演信號(hào)，其實(shí)現(xiàn)框圖如圖10所示。該方式較為靈活，限制較少，難點(diǎn)在于切換如何低成本地準(zhǔn)確采集相位信息。

圖10 數(shù)字信號(hào)頻域時(shí)間反演處理原理框圖

4 系統(tǒng)框架設(shè)計(jì)

對(duì)上述幾種方法進(jìn)行總結(jié)整理如表1所示。

表1 時(shí)間反演實(shí)現(xiàn)方法對(duì)比

本文針對(duì)啁啾變換、時(shí)間透鏡、方向回溯等模擬方法其啁啾期間參數(shù)固定，無(wú)法調(diào)整頻率，而數(shù)字處理方法實(shí)時(shí)性差，工作頻率受限，不能處理非周期信號(hào)等痛點(diǎn)，綜合考慮了反演質(zhì)量和成本，結(jié)合了模擬和數(shù)字方法，提出了一種面向電磁場(chǎng)賦形的時(shí)間反演綜合設(shè)計(jì)方法，做到高效實(shí)時(shí)的時(shí)間反演處理。

圖11 所示為該電磁場(chǎng)賦形系統(tǒng)時(shí)間反演處理部分的框架設(shè)計(jì)圖，每個(gè)時(shí)間反演天線陣元和射頻芯片的單刀雙擲開(kāi)關(guān)相連，對(duì)應(yīng)發(fā)射和接收，分別用在時(shí)間反演過(guò)程的導(dǎo)引階段和回傳階段。矢量網(wǎng)絡(luò)分析儀（Vector Network Analyzer，VNA）使用本實(shí)驗(yàn)室自研發(fā)的基于FPGA 的便攜式矢量網(wǎng)絡(luò)分析儀，利用該儀器對(duì)某C 波段串口多功能芯片進(jìn)行驅(qū)動(dòng)控制，程控完成放大（兩次）、開(kāi)關(guān)（兩次）、移相和衰減等操作。賦形過(guò)程中，導(dǎo)引端發(fā)送的導(dǎo)引信號(hào)，使時(shí)間反演天線陣列接受導(dǎo)航信號(hào)；FPGA控制射頻芯片里的開(kāi)關(guān)進(jìn)入前端接收鏈路，F(xiàn)PGA 控制每個(gè)射頻芯片里接收鏈路的開(kāi)關(guān)，讓矢網(wǎng)能夠分時(shí)采樣各通道的相位信息，由FPGA 計(jì)算相位補(bǔ)償信息?；貍麟A段，F(xiàn)PGA 控制射頻芯片里移相器和放大器，完成進(jìn)行移相和放大。對(duì)于每一路信號(hào)，其傳輸和處理流程如圖12所示。

圖11 電磁場(chǎng)賦形系統(tǒng)時(shí)間反演處理部分框架設(shè)計(jì)圖

圖12 陣列信號(hào)時(shí)間反演處理流程

5 時(shí)間反演賦形場(chǎng)仿真

本文構(gòu)建的信道模型如圖13（a）所示，在墻壁多徑環(huán)境下，使用4×4 的時(shí)間反演面天線陣列作為接收端，兩個(gè)天線作為導(dǎo)引端，構(gòu)建多目標(biāo)電磁場(chǎng)賦形電磁仿真模型。多目標(biāo)墻壁多徑環(huán)境賦形仿真模型如圖13（a）所示，為了模擬墻壁多徑環(huán)境，在模型的六面分別放置墻壁材料，墻壁相對(duì)介電常數(shù)取6.25［20］，墻壁厚度為20mm。賦形天線和時(shí)間反演天線均使用諧振頻率5.8GHz的微帶天線。16只時(shí)間反演天線圍成一個(gè)陣列大小200mm×200mm的天線陣，天線陣元的距離50mm，構(gòu)成時(shí)間反演面天線陣列。賦形天線AN17 和AN18 位于時(shí)間反演面天線陣列上方150mm 處，坐標(biāo)為（-50，-50，150）和（50，50，150）。

圖13 墻壁環(huán)境信道模型

可以做出如圖14 所示的電磁場(chǎng)賦形數(shù)字。具體步驟如下：

圖14 時(shí)間反演進(jìn)行數(shù)字賦形仿真圖

1）根據(jù)預(yù)期形狀，以導(dǎo)引天線為像素點(diǎn)，構(gòu)建各導(dǎo)引天線的位置；

2）各導(dǎo)引天線發(fā)出導(dǎo)引信號(hào)，時(shí)間反演天線采集相位信息，記錄并保存時(shí)間反演面天陣列線第1，…，16 個(gè)陣元采集到的信號(hào)相位φ1，φ2，…，φ16，存入數(shù)據(jù)庫(kù)；

3）能量回傳過(guò)程撤去導(dǎo)引天線的激勵(lì)，同時(shí)撤去導(dǎo)引天線，然后將幅度為1，相位為-φ1，-φ2，…，-φ16的回傳信號(hào)依次導(dǎo)入時(shí)間反演陣列第1，…，16個(gè)陣元；

4）記錄焦平面的能量分布，即為賦形電磁場(chǎng)；

5）再次需要該賦形電磁場(chǎng)時(shí)，直接導(dǎo)出數(shù)據(jù)庫(kù)中該形狀的共軛相位信息，無(wú)需重復(fù)步驟1）和2）。

要得到更加復(fù)雜的賦形電磁場(chǎng)，則需要更多的像素點(diǎn)，即放置更多的導(dǎo)引天線。實(shí)際仿真中發(fā)現(xiàn)，過(guò)多的導(dǎo)引天線會(huì)影響時(shí)間反演賦形電磁場(chǎng)的分辨效果，這是因?yàn)椴荒車?yán)格滿足信道互易性條件，導(dǎo)引天線的撤去對(duì)整個(gè)信道的影響已經(jīng)不能忽略了，所以需要更大規(guī)模的時(shí)間反演天線陣列，才能得到更加精細(xì)的賦形電磁場(chǎng)。賦形微波場(chǎng)在EMC 測(cè)試、電磁偽裝、賦形輻照源、區(qū)域干擾、腫瘤化療、醫(yī)學(xué)診斷、電磁成像等領(lǐng)域有廣泛應(yīng)用［21］。

6 結(jié)語(yǔ)

時(shí)間反演電磁波具有良好的空時(shí)聚焦特性，可用于電磁場(chǎng)賦形。本文探究了啁啾變換方法、時(shí)間透鏡方法、方向回溯方法和數(shù)字處理方法這幾種時(shí)間反演的實(shí)現(xiàn)方法，并進(jìn)行了數(shù)值仿真與分析，提出了一種面向電磁場(chǎng)賦形的時(shí)間反演異構(gòu)設(shè)計(jì)。本文基于此設(shè)計(jì)構(gòu)建了墻壁環(huán)境的信道模型，并基于時(shí)間反演完成了特定圖案的電磁場(chǎng)賦形。本文的工作對(duì)時(shí)間反演在醫(yī)學(xué)診療、工業(yè)加工、電磁干擾、無(wú)線傳能等領(lǐng)域的實(shí)際應(yīng)用具有一定的參考意義。

然而，有些方面的研究還有待進(jìn)一步深入，具體如下：1）本文仿真偏重于對(duì)通過(guò)仿真研究賦形，而在實(shí)際環(huán)境下的效果，還有待實(shí)驗(yàn)驗(yàn)證；2）本文所述的時(shí)間反演賦形研究在實(shí)際中應(yīng)用前景廣泛，但由于算力限制，只基于4×4 天線陣列設(shè)計(jì)一個(gè)仿真展示。后期如果有足夠的算力支持，可以設(shè)計(jì)更大規(guī)模的天線陣列，實(shí)現(xiàn)更高分辨率的磁場(chǎng)賦形仿真，并進(jìn)行實(shí)物驗(yàn)證。更高分辨率的磁場(chǎng)賦形在癌癥診療、電磁成像、無(wú)線傳能等領(lǐng)域具有較高的應(yīng)用價(jià)值。