周洪澄 王秉中 丁 帥 歐海燕

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

(2012年11月29日收到;2013年2月22日收到修改稿)

1 引言

時間反演技術(shù)最開始在聲波中得以提出和應(yīng)用，并于2004年被引入電磁學(xué)領(lǐng)域[1].電磁學(xué)中的時間反演是指：天線接收到一串電磁信號之后，將其在時域上進行反轉(zhuǎn)，并重新發(fā)射出去的過程.這類似于堆棧操作，最先接收到的信號最后發(fā)射出去，而最后接收到的信號則最先發(fā)射.2007年，F(xiàn)ink等人把有限個接收天線組成時間反演鏡[2](time reversalmirror，TRM)放置在具有豐富多徑的環(huán)境中進行實驗，驗證了時間反演電磁波在近場[3]和遠場[4]都具有超分辨率聚焦特性.研究發(fā)現(xiàn)，時間反演電磁波具有時空聚焦特性[5]，即：經(jīng)過時間反演處理的電磁波在時間域內(nèi)具有時間壓縮特性、在空間域內(nèi)具有空間聚焦特征[6,7].因此，時間反演技術(shù)可廣泛用于復(fù)雜媒質(zhì)中目標(biāo)的探測與成像、超寬帶無線通信等領(lǐng)域[8-11].

2010年，F(xiàn)usco等人利用相位共軛透鏡(phase conjugated mirror，PCM)，研究了相位共軛電磁波的近、遠場超分辨率聚焦特性，對頻域信號的相位共軛所具有的超分辨率特性進行了理論分析[12-14].研究表明，凋落波借助于散射體轉(zhuǎn)換為傳播波，激發(fā)物體表面的等離子激元，是相位共軛具有超分辨率聚焦特性的重要原因.由傅里葉變換關(guān)系可知，信號在時域進行時間反演等效于在頻域上進行相位共軛.與相位共軛鏡使用單頻信號相比，時間反演鏡使用的是寬帶或超寬帶脈沖，它允許對任何時間窗內(nèi)的信號進行時間反演處理.

自從Pendry預(yù)測負折射材料可以構(gòu)成一個完美透鏡后[15]，出現(xiàn)了多種可以用于亞波長天線成像的完美透鏡.Fink與Fusco等都曾指出，單純的時間反演技術(shù)并不能突破衍射極限、獲得電磁波的超分辨率聚焦特性，只有在目標(biāo)體的近場周圍加入散射體，才能獲得這種特性[16-18].本文在已有的研究基礎(chǔ)上，提出了利用一種規(guī)則排列的金屬絲陣列來實現(xiàn)時間反演的超分辨率聚焦成像的方法.這里說的金屬絲陣列，是一組平行的理想導(dǎo)線，即損耗可忽略不計，如圖1所示.這種結(jié)構(gòu)與左手材料不一樣，它沒有負折射，并且Belov等已經(jīng)證實該種金屬絲陣列可以為凋落波提供一種的傳播渠道[19]，因此它是一種不同于完美透鏡的亞波長成像透鏡.

文獻[19]利用金屬絲陣列能夠傳遞凋落波這一特性，成功地在源信號的異側(cè)實現(xiàn)了對源信號信息的提取.可是，一旦源信號停止激勵，則無法準(zhǔn)確地確定源信號的位置.2010年，Ding等人在時間反演技術(shù)的支持下，利用金屬絲陣列證實了時間反演的時空聚焦特性、實現(xiàn)了對源信號的提取，且在源信號停止激勵后仍舊能夠確定目標(biāo)天線的位置及目標(biāo)天線的形狀，得到了具有超分辨率特征的源像[20].本文在文獻[20]的基礎(chǔ)上，進一步證實利用金屬絲陣列的傳播特性和時間反演技術(shù)能夠確定源像的位置，首次提出利用金屬絲陣列在異地實現(xiàn)對源像恢復(fù)的方法，并在仿真軟件CST Microwave Studio中驗證了該方法.

2 線陣式時間反演鏡的聚焦成像

凋落波能夠攜帶亞波長信息，但在普通的導(dǎo)波結(jié)構(gòu)，甚至是自由空間中，凋落波的傳播通常會呈現(xiàn)指數(shù)型衰減，因此很難觀測到凋落波，也就很難獲得亞波長信息.傳統(tǒng)的成像系統(tǒng)只能傳遞源信號的傳播波，并受到衍射極限的限制.而本文利用的這種金屬絲陣列結(jié)構(gòu)，能夠為凋落波提供傳播渠道，使得源信號所成的像可以在透鏡的另一端被觀測到.因此可以推斷，該金屬絲媒質(zhì)可以用于對微波源進行亞波長成像.

本文成像所采用的結(jié)構(gòu)如圖1所示.其中，圖1(a)為放置在一個完全封閉的諧振腔中的一個31×31的正方形規(guī)則排列的理想金屬導(dǎo)線陣列，其中金屬絲半徑r=1 mm，相鄰金屬絲之間的中心間距a=3 mm，金屬絲的長度均為d=110 mm.如圖1(b)所示，這個金屬絲陣列可以等效為一種特殊介質(zhì)，它不僅能為傳播波提供豐富的多徑和散射，還能為凋落波到傳播波的轉(zhuǎn)換提供一種渠道.即使在低頻階段，該金屬絲陣列結(jié)構(gòu)仍有強烈的空間色散特性[21]，因此本文選擇該種金屬絲陣列作為亞波長成像的傳輸渠道，來研究時間反演電磁波的時空聚焦特性以及其超分辨率聚焦結(jié)果.在三維空間中，其該種金屬絲陣列的介電常數(shù)可用張量形式表示為

其中，x方向是金屬絲陣列的伸展方向，qx是波矢q在x方向上的分量，kp是對應(yīng)于等離子體頻率ωp的波數(shù).因此在此介質(zhì)的波長應(yīng)表示為

在這種結(jié)構(gòu)中，kp取決于金屬絲的間距a以及金屬絲的半徑r，公式如下[19]：

圖1為用于實現(xiàn)異地成像的結(jié)構(gòu)圖，由圖1可看出，目標(biāo)天線和時間反演鏡天線陣列分別位于金屬絲陣列的兩側(cè)，時間反演鏡呈線陣式分布(從右至左依次為天線2,3,···,8).這8個天線的結(jié)構(gòu)一樣，均為理想偶極子天線.這些偶極子天線的長度均為20 mm，且目標(biāo)天線距離金屬絲陣列的距離為2 mm，時間反演鏡中的所有天線距離金屬絲陣列的距離為5 mm.

圖1 成像結(jié)構(gòu)圖 (a)三維結(jié)構(gòu)圖;(b)平面圖

在時間反演技術(shù)下，利用該金屬絲陣列在不同位置恢復(fù)源像的具體步驟描述如下：

1)由目標(biāo)天線(1號天線)發(fā)射一個頻譜帶寬為0.1—6 GHz的調(diào)制高斯脈沖信號x(t)，信號波形如圖2所示.

2)時間反演鏡中的天線 i(i=2,3,···,8)接收到的信號 yi(t)(i=2,3,···,8)表示為

其中hi(t)是目標(biāo)天線與天線i之間的信道沖激響應(yīng).圖3(a)顯示的是接收天線2接收到的信號，由于金屬絲陣列的散射作用，復(fù)雜的多路徑效應(yīng)顯而易見.

3) 將 接 收 到 的 信 號 yi(t)(i=2,3,···,8) 分別在時域上進行時間反演操作，記為yi(t)(i=2,3,···,8).圖 3(b)表示的是天線2的接收信號經(jīng)時域反轉(zhuǎn)之后的信號.

5)取消目標(biāo)天線，利用CST仿真軟件，在目標(biāo)天線的原位置設(shè)置一個場觀測面，如圖5所示，其中場觀測面的頻率設(shè)置在3.5 GHz.

圖2 激勵信號

圖3 接收天線2的信號 (a)接收到的信號;(b)時域反轉(zhuǎn)后的信號

圖4 天線1接收到的信號

6)將時域反轉(zhuǎn)信號作為源信號，分別由對應(yīng)的7個接收天線同時發(fā)射出去，即天線i的激勵信號為 yi(-t)，i=2,3,···,8.在觀測面上顯示的最佳聚焦效果如圖6所示.

7)以時間反轉(zhuǎn)鏡中心位置的天線5為對稱點，將其它天線的時域反轉(zhuǎn)信號進行左右對稱交換之后再同時發(fā)射，即天線i的激勵信號為y5-(i-5)(-t)，i=2,3,···,8.這種時間反演激勵下，在觀測面上顯示的最佳聚焦效果如圖7所示.

通常環(huán)境下，時間反演聚焦特性受到瑞利判據(jù)的限制，聚焦分辨率是λ/2.而在本文所構(gòu)建的仿真模型，能夠?qū)崿F(xiàn)電磁波超分辨率的聚焦特性.圖6所示為t=181.86 ns時(聚焦效果最好的時刻)的觀測面上呈現(xiàn)出的聚焦效果，從圖中可以清晰地看到偶極子天線的形狀.時間反演鏡天線反演之后發(fā)射出的信號，能夠在目標(biāo)天線原位置很好實現(xiàn)時空聚焦，而且成像尺寸為金屬絲的三個結(jié)構(gòu)周期，約為λ/14，突破衍射極限，實現(xiàn)了亞波長成像.

圖5 設(shè)置場觀測面

圖6 原地成像效果

圖7同樣為t=181.86 ns時(聚焦效果最好的時刻)的觀測面上呈現(xiàn)出的聚焦效果，由圖7可看出，將時間反演鏡時域反轉(zhuǎn)后的信號以天線5為中心左右對稱交換之后再發(fā)射，同樣能夠突破衍射極限，并在異于源像的位置出現(xiàn)了清晰的成像效果.

文獻[20]中指出，如果不使用時間反演技術(shù)，將時間反演鏡接收到的信號直接發(fā)射出去，即使是穿過金屬絲陣列，也只能在目標(biāo)天線對應(yīng)的位置上，出現(xiàn)時間反演鏡天線的模糊成像，無法在觀測面上得到清晰的聚焦成像結(jié)果.但由圖6和圖7可以知道，在時間反演技術(shù)的支持下，只要接收天線陣列接收到了目標(biāo)天線發(fā)射出的信號，并經(jīng)過時間反演技術(shù)處理，即使取消了目標(biāo)天線，同樣能夠在目標(biāo)天線的觀測面上得到良好的聚焦效果.與此同時，根據(jù)圖6和圖7可知，原地和異地的成像能量分布一致、成像位置對稱，并且兩者在同一時刻得到空間上的最佳聚焦效果.因此，利用時間反演技術(shù)，改變信號對時間反演鏡天線的激勵方式，能夠在除了源像原位置之外的其他位置實現(xiàn)對源信號的提取和分析.圖7所示的異地成像效果滿足遠場成像條件(r＞2D2/λ)，即利用該方法能夠在遠場實現(xiàn)超分辨率聚焦成像.

圖7 異地成像效果

3 面陣式時間反演鏡的聚焦成像

上一小節(jié)證實了利用金屬絲陣列結(jié)合時間反演技術(shù)能夠在兩個不同的位置恢復(fù)目標(biāo)天線的信號，這一節(jié)將介紹一種能夠在4個不同的位置實現(xiàn)對源信號恢復(fù)的情況.采用的結(jié)構(gòu)設(shè)置與上一節(jié)相似，金屬絲陣列保持不變，僅改變目標(biāo)天線和時間反演鏡的分布和數(shù)目.其中，1號天線為目標(biāo)天線，2—6號天線組成面陣式分布的時間反演鏡.天線2處于金屬絲陣列的中心位置，天線3—6以天線2為中心呈對稱分布，具體分布如下圖8所示.

圖8 四個成像結(jié)果的結(jié)構(gòu)設(shè)置圖

具體成像步驟與上一節(jié)的類似，首先由目標(biāo)天線發(fā)射一個如圖2所示的時域信號x(t)，然后時間反演鏡接收到的信號記為 ym(t)(m=2,3,···,6)，其時域反轉(zhuǎn)之后的信號記為 ym(-t)(m=2,3,···,6)，接著取消目標(biāo)天線，利用CST仿真軟件，在目標(biāo)天線的原位置設(shè)置一個場觀測面，場觀測面的頻率設(shè)置在3.5 GHz，最后分別按照如表1所示的四種方式(S1—S4)對時間反演鏡天線進行激勵.

表1 時域反轉(zhuǎn)信號對時間反演鏡的激勵方式

表 1中的 S1，S2，S3，S4四種情況，分別表示將時間反演鏡接收信號做時域反轉(zhuǎn)后的信號用于直接激勵、左右對稱交換后激勵、上下對稱交換后激勵以及對角交換后激勵，由此分別得到如圖9(a)—(d)所示的四個成像結(jié)果.圖9展示的是在同一時刻點，不同激勵方式下最佳的聚焦成像效果.由圖9可以看出，通過改變信號的激勵方式，能夠在不同的位置實現(xiàn)對源像的恢復(fù)，每個成像效果都突破衍射極限，成像尺寸約為λ/8，實現(xiàn)了超分辨率成像，并且在除了目標(biāo)天線原位置之外的三個遠場位置實現(xiàn)了對源像的恢復(fù).

圖9 成像結(jié)果圖 (a)原地成像;(b)左右交換;(c)上下交換;(d)對角交換

本文還分析了如表2所示的兩種方式，將時間反演鏡接收信號做時域反轉(zhuǎn)后的信號用于激勵.表2表示的是時間反演鏡天線處理激勵信號的過程中，只是對一部分信號進行了對稱交換，分別得到如圖10(a)和(b)所示的兩種結(jié)果.由圖10可以看出，如果只對稱交換一部分接收天線時域反轉(zhuǎn)后的信號，將只能在觀測面上得到模糊的成像效果，我們無法通過成像結(jié)果來分析和提取源信號的信息.因此在信號交換的過程中，須對稱地交換所有的信號，才能在不同的位置實現(xiàn)對源像的恢復(fù).

表2 只對稱交換部分信號的激勵方式

圖10 對稱交換部分信號的成像效果 (a)只交換3號和5號天線;(b)只交換4號和6號天線

4 目標(biāo)天線不同位置的聚焦成像

接下來討論在時間反演鏡陣列一定的情況下，通過改變目標(biāo)天線的位置，能夠在另外4個位置得到相應(yīng)的成像效果.此時的結(jié)構(gòu)設(shè)置與圖8相似，唯一的變化就是將目標(biāo)天線的位置向下移動了15 mm.激勵方式與表1所示一致，分別為：直接激勵、左右交換后激勵、上下交換后激勵以及對角交換后激勵.得到的四個成像效果分別如下圖11(a)—(d)所示.

由圖11可知，即使改變目標(biāo)天線的位置，通過時間反演技術(shù)處理之后，同樣能夠在成像平面上的4個不同位置得到超分辨率聚焦成像效果(成像尺寸約為λ/7)，實現(xiàn)對源像的恢復(fù).可以推斷，如果金屬絲介質(zhì)的某一部分出現(xiàn)問題，則會在成像平面上出現(xiàn)不同的能量分布和成像效果[13].通過改變目標(biāo)天線的位置和觀察成像效果圖，能夠?qū)崿F(xiàn)對結(jié)構(gòu)的檢測和掃描.

綜合以上仿真結(jié)果圖，可以發(fā)現(xiàn)，在目標(biāo)天線位置和初始激勵信號都不改變的時候，不同的時間反演鏡陣列將會在不同的地方實現(xiàn)對源像的恢復(fù)，并且成像分辨率也不一樣;此外，在時間反演鏡陣列不改變的情況下，改變目標(biāo)天線的位置，也將會得到不同效果的成像精度.這是因為，改變時間反演鏡排布或是改變目標(biāo)天線的位置之后，將會改變目標(biāo)天線與時間反演鏡之間的多徑分布，使得時間反演鏡接收到不同的信號，最終導(dǎo)致了在觀測面上得到不同的成像效果.綜合分析，成像效果的主要影響因素有：1)近遠場轉(zhuǎn)換效率.轉(zhuǎn)換效率將主要取決于目標(biāo)天線與時間反演鏡陣列之間的多徑分布，因此金屬絲的排布方式將會影響近遠場轉(zhuǎn)換效率.2)傳播過程中自由空間的波譜損失.因為凋落波與傳播波之間的互換過程發(fā)生在金屬絲陣列中，但并不是所有的波(凋落波/傳播波)都能夠到達時間反演鏡平面，因此盡可能地將目標(biāo)天線和接收天線排布得離金屬絲近一些.3)傳播過程中的色散效應(yīng).在整個實驗過程中，金屬絲的電尺寸、信號在金屬絲陣列中接收和發(fā)送形式的差異，都將帶來圖像的不同擴展[13].

圖11 改變目標(biāo)天線位置后的成像結(jié)果圖 (a)原地成像;(b)左右交換;(c)上下交換;(d)對角交換

5 結(jié)論

本文利用時間反演技術(shù)，對規(guī)則排列金屬絲陣列的超分辨率聚焦成像特性進行了研究.利用一種已有的規(guī)則排列的金屬絲陣列，對時間反演電磁波的超分辨率聚焦成像特性展開了研究.文中證實，在不同的信號激勵方式下，除了能在原地實現(xiàn)對源信號的提取和突破衍射極限聚焦之外，還能在同一時刻的不同位置實現(xiàn)對源像的恢復(fù)，并得到了多種超分辨率聚焦成像的仿真結(jié)果，最佳聚焦效果達到了λ/14.在時間反演鏡位置一定的情況下，改變目標(biāo)天線的位置依然能實現(xiàn)聚焦成像.但是，若只交換部分時間反演鏡中時域反轉(zhuǎn)后的信號，將無法在觀測面上得到清晰地超分辨率成像.這些結(jié)果表明，我們可以在傳統(tǒng)的材料和設(shè)備下，利用時間反演技術(shù)在源像的遠處實現(xiàn)對微細特征的提取，這種成像結(jié)構(gòu)對異地恢復(fù)源信號、遠場超分辨率成像以及結(jié)構(gòu)探測都有著重要的意義.

[1]Lerosey G，Rosny J，Tourin A，Derode A，MontaldoG，F(xiàn)ink M 2004 Phys.Rev.Lett.92 1

[2]Ding S，Wang B Z，Ge G D，Wang D 2011 Acta Phys.Sin.60 104101(in Chinese)[丁帥，王秉中，葛廣頂，王多2011物理學(xué)報60 104101]

[3]Rosny J，F(xiàn)ink M 2007 Phys.Rev.A 76 1

[4]Lerosey G，Rosny J，Tourin A，F(xiàn)ink M 2007 Science 315 1120

[5]Ge G D，Wang B Z，Huang H Y，Zheng G 2009 Acta Phys.Sin.58 8249(in Chinese)[葛廣頂，王秉中，黃海燕，鄭罡2009物理學(xué)報58 8249]

[6]Chen Y M，Wang B Z，Ge G D 2012 Acta Phys.Sin.61 024101(in Chinese)[陳英明，王秉中，葛廣頂2012物理學(xué)報61 024101]

[7]Chen Y M，Wang B Z 2012 Chin.Phys.B 21

[8]Zhang Z M，Wang B Z，Ge G D 2012 Acta Phys.Sin.61 098401(in Chinese)[章志敏，王秉中，葛廣頂2012物理學(xué)報61 098401]

[9]Liu X F，Wang B Z，XiaoS，Deng J 2008 Progress in Electromagnetics Research,PIER 79 401

[10]Liu X F，Wang B Z，XiaoS，LaiS 2010 IEEE Trans.Antennas Propag.58 1731

[11]Liu X F，Wang B Z，LiL 2009 IEEE Antennas and Wireless Propag.Letters 9 1429

[12]Malyuskin O，F(xiàn)uscoV 2010 IEEE Trans.Antennas Propagat.58 2884

[13]Malyuskin O，F(xiàn)uscoV 2010 IEEE Trans.Antennas Propagat.58 459

[14]Malyuskin O，F(xiàn)uscoV 2010 Proc.URSI InternationalSymposium on Electromagnetic Theory.Berlin 2010 p235

[15]Pendry J B 2010 Phys.Rev.Lett.85 3966

[16]Ge G D，Wang B Z，Wang D，ZhaoD，Ding S 2011 IEEE Trans.Antennas Propagat.59 4345

[17]Ge G D，Wang D，Wang B Z 2011 Progress In Electromagnetics Research 114 429

[18]Ge G D，Zang R，Wang D，Ding S，Wang B Z 2011 IET Electronics Letters 47 901

[19]Belov P A，Silveirinha M G 2006 Phys.Rev.73 056607

[20]Ding S，Wang B Z，Ge G D 2010 Microwave and Millimeter Wave Technology(ICMMT),InternationalConference on p1160

[21]Belov P A，Marques R，MaslovskiS，Nefedov I，Silverinha M，SimovskiC，Tretyakov S 2003 Phys.Rev.B 67 113103