超寬帶穿墻雷達成像的FDTD數(shù)值模擬

王芳芳 張業(yè)榮

(南京郵電大學(xué)電子科學(xué)與工程學(xué)院,江蘇南京210003)

超寬帶穿墻雷達成像的FDTD數(shù)值模擬

王芳芳 張業(yè)榮

(南京郵電大學(xué)電子科學(xué)與工程學(xué)院,江蘇南京210003)

研究了超寬帶雷達在穿墻探測中的保真度,通過FDTD建模產(chǎn)生精確的穿墻雷達數(shù)據(jù),利用后向投影算法對模擬得到的雷達數(shù)據(jù)進行成像。結(jié)果驗證:超寬帶雷達能準確探測和定位二維房間內(nèi)的物體。成像算法中必須考慮墻體作用,電磁波在穿墻傳播時發(fā)生折射、速度的改變、信號的衰減,這些將導(dǎo)致像的散焦、偏移,甚至可能產(chǎn)生虛像。運用迭代方法修正了成像算法中的時間因子,可有效地抑制像的偏移。

超寬帶雷達;穿墻;時域有限差分;后像投影算法

1.引 言

大規(guī)模城區(qū)巷戰(zhàn)的發(fā)展和武器隱身技術(shù)的進步對雷達的探測性能提出了更高的要求,能對葉簇、地表甚至墻壁、建筑、堡壘等非電磁透明介質(zhì)屏障進行穿透探測并分析隱蔽目標(biāo)的雷達已經(jīng)引起了越來越多的關(guān)注。穿墻探測雷達能夠穿透非金屬墻壁,并對墻壁后面的運動目標(biāo)進行探測、定位和跟蹤。隨著雷達技術(shù)的發(fā)展,穿墻探測在國防以及民用上得到了廣泛的應(yīng)用[1,4-5]。

超寬帶雷達有三種不同類型:頻率調(diào)制,頻率步進調(diào)制[7],短脈沖,我們應(yīng)用的是短脈沖雷達,超寬帶短脈沖雷達在穿墻成像方面有較大的優(yōu)勢,它的發(fā)射時寬為幾個納秒的超寬帶信號,相對帶寬大于25%,具有極高的距離向分辨率,同時超寬帶信號豐富的低頻分量保證其良好的傳輸特性和較大的作用距離。

通過數(shù)值計算得到的模擬雷達接收信號經(jīng)過后向投影算法可以進行成像,該算法在穿墻成像中同樣可以應(yīng)用,但是電磁波在穿透墻壁的同時會發(fā)生折射、傳播速度的改變、信號幅度的衰減。所以必須考慮墻體作用所引起信號傳播時間的變化。本文第二部分計算了該時間偏差并在成像算法中對其進行了補償。墻的結(jié)構(gòu)、厚度、介電常數(shù)、電導(dǎo)率、信號的入射角都會影響穿墻特性,如果在成像算法中沒有考慮這些,所得到的像會出現(xiàn)不同程度的偏差甚至畸變。本文提出了一種簡單的迭代方法來尋找折射點,從而找到了折射路線,正確修正了用于后向投影算法中信號從發(fā)射天線到接收天線的傳播時間。

文章首先建立了穿墻的物理模型,應(yīng)用折射定律得到正確的傳播路徑,第三部分介紹了BP成像算法,第四部分用FDT D數(shù)值模擬[2-4,6,8]驗證了該系統(tǒng)的正確性[9],比較了沒有考慮墻和考慮墻的成像情況。

2.電磁波傳播物理模型

電磁波在穿透墻壁的時候,在墻體和空氣的分界面上會發(fā)生折射,改變了電磁波的傳播路徑,為了正確地計算傳播時間,我們需要找到正確的傳播路徑。

信號從發(fā)射天線經(jīng)過墻體傳播到目標(biāo)發(fā)生了二次折射,簡單的幾何證明可知只要天線和目標(biāo)的位置確定下來,從發(fā)射天線到目標(biāo)的傳播距離就確定下來了,而與天線距離墻體的遠近沒有關(guān)系,所以我們考慮最特殊的情況,設(shè)天線緊貼墻壁放置。

如圖1所示,墻體厚度為dw,相對介電常數(shù)為εr,在直角坐標(biāo)系x-o-y中,發(fā)射天線P(px,py),目標(biāo)Q(qx,qy),折射點W(wx,wy)。設(shè)矢量F為點W到Q的矢量,如圖建立坐標(biāo)系A(chǔ)-P-B,在該坐標(biāo)系下,Q的坐標(biāo)為V(va,vb),W 的坐標(biāo)為(wa,wb),F的坐標(biāo)為(fa,fb),由幾何關(guān)系可得

圖1 折射模型

考慮兩種極限情況,假設(shè)墻體為空氣時,即介電常數(shù)為εr=1,傳播路徑為直線PW2V,即假設(shè)不存在折射現(xiàn)象,交點為W2(w2,wb),另外再假設(shè)墻體介電參數(shù)為無窮大時,電磁波在墻內(nèi)可近似垂直于墻面?zhèn)鞑?傳播路徑為 PW1V,交點為 W1(w1, wb),實際上,電磁波在兩種介質(zhì)分界面上必須滿足折射定律,所以可知w1＜wa＜w2,有幾何關(guān)系如下

式中:κ為折射率;入射角為θi;折射角為θr.

通過迭代法求得實際折射點W(wa,wb),設(shè)wa初始值為(w1,w2)區(qū)間中點w0,通過幾何關(guān)系以及(1)(2)(3)式求得κ,將其與墻體介電常數(shù)εr進行比較,如果κ＞εr,則wa取(w0,w2)區(qū)間中點,如果κ＜εr,則wa取(w1,w0)區(qū)間中點,以此迭代下去直到在容許誤差范圍內(nèi)κ=εr成立而結(jié)束,最后求得折射點W(wa,wb).

設(shè)墻內(nèi)傳播距離為l1,非墻體內(nèi)傳播距離為l2

def f為有效傳播距離,為電磁波在墻內(nèi)的傳播速度,則,

3.成像算法

3.1 BP算法

FDTD模擬的接收信號數(shù)據(jù)是離散時間函數(shù)的電場值,每個抽樣時間稱為一個時間位,后向投影算法首先計算了從發(fā)射天線到成像區(qū)域每個網(wǎng)格(像素點)再到接收天線總的往返時間,在模擬數(shù)據(jù)中找到對應(yīng)時間的電場幅度,然后對所有信道電場值進行疊加,在目標(biāo)位置處,信號相干疊加,而在非目標(biāo)處,信號非相干疊加,最后可得目標(biāo)成像圖。BP算法的具體實現(xiàn)步驟如下:

◆ 把整個成像區(qū)域劃分成網(wǎng)格得到像素點;

◆ 對每一個像素點計算從發(fā)射天線到該像素點再從本像素點到接收天線的往返時間;

◆ 對應(yīng)某個像素點根據(jù)往返時間記錄所有接收天線的電場幅度;

◆ 對應(yīng)該像素點對所有記錄值進行幅度累加;

◆ 對所有的像素點重復(fù)b,c,d三個步驟;

設(shè)任意像素點的坐標(biāo)為(xi,yj),像素點的像素為S(xi,yj)

式中:c為自由空間的光速;(xT,yT)為發(fā)射天線的坐標(biāo);(xR(k),yR(k))為第k個接收天線坐標(biāo)。

3.2 考慮時間修正的BP算法

在BP算法中,對往返時間的正確計算是極其重要的一步,在穿墻系統(tǒng)中,由于墻的作用,電磁波的傳播發(fā)生了折射,從而改變了傳播路徑,電磁波在墻內(nèi)的速度也不再是光速。為了正確地應(yīng)用BP算法,必須對總的往返時間進行修正,用實際物理模型計算出有效路徑deff,從而計算出往返時間tij,具體見第2章的式(4)。

4.FDT D建模仿真結(jié)果

這部分用FDT D建模產(chǎn)生了穿墻成像中的模擬雷達數(shù)據(jù),房間模型如圖2所示。

房間長寬分別為359 cm,236 cm,墻壁的介電常數(shù)為2.75,電導(dǎo)率為0.003,在(-30 cm,95 cm)處放置邊長為10 cm的導(dǎo)體箱子。發(fā)射天線位于墻壁前面5 cm處,接收天線距離墻壁4 cm,在縱向上每隔設(shè)置一個接收天線,總的孔徑長度為240 cm。發(fā)射信號為調(diào)制高斯脈沖,時寬為0.6 ns,中心頻率為4 GHz,時域、頻域波形如圖3所示

圖2 房間模型

圖3. 發(fā)射信號時域,頻域波形圖

FDT D的空間步為dx=dy=1 cm,劃分網(wǎng)格數(shù)為300×425,時間步為dt=19.25×10-12s,整個計算區(qū)域用PML(最佳匹配層)吸收邊界條件截斷。

用BP算法對房間內(nèi)物體成像的時候,實際上有用信號是物體散射的電磁波,為了防止房間墻壁的影響,我們需要把背景信號去掉。用以下方法得到物體的散射信號,先對房間內(nèi)沒有物體的情況進行仿真,記錄接收天線的信號,再把物體放回房間,所有的參數(shù)設(shè)置與房間內(nèi)無物體情況設(shè)置相同,記錄接收天線的信號,有物體在房間內(nèi)情況下的模擬接收信號減去沒有物體在房間內(nèi)情況下的模擬接收信號,此時該信號即為物體的散射信號。下面列出了第41個抽樣位置處接收天線的的信號,即在該位置處天線接收到的物體散射信號,如圖4所示,圖中信號1表示物體本身的散射信號,信號2表示由于后面墻壁的作用,經(jīng)物體散射的信號作用在后墻上再次引起的散射信號,等效為二次散射信號,該信號會在仿真空間外形成對應(yīng)的虛像。

圖4 第41個接收天線位置的接收信號

圖5～圖7分別為在成像算法中沒有考慮墻的影響,在成像算法中對墻體作用進行了補償,在仿真中沒有設(shè)置墻壁的成像情況。圖7情況也即單純的物體散射問題。其中圖6(a)、(b)分別為相應(yīng)的二維、三維圖,(c)、(d)分別為在方位向和距離向上的幅度分布圖。

圖中交叉點表示發(fā)射天線,墻前面的離散點表示采樣位置處放置的接收天線,虛線外框,內(nèi)框分別表示外墻和內(nèi)墻,方框為金屬物體。比較圖5和圖6可以看出在對時間因子修正前后像的位置發(fā)生了變化,未考慮墻體作用時像與實際位置存在一定的偏差,而考慮墻體的折射作用后,即對墻壁作用進行補償后像在金屬物體實際位置處。從而可以看出墻壁的影響不可忽略,在成像算法中必須考慮墻體的作用。由圖6(c)、(d)可以看出峰值點在物體的實際位置(-0.3 m,095 m)附近,圖中像出現(xiàn)了一定程度的散焦,特別是在方位向上,但是從波瓣寬度來看,應(yīng)用超寬帶雷達技術(shù)可以達到很好的距離向和方位向分辨率。

圖5 沒有考慮墻體作用的成像結(jié)果

圖7 沒有房間的物體的成像結(jié)果

5.結(jié) 論

用FDT D數(shù)值方法模擬了超寬帶雷達在穿墻成像中的應(yīng)用,從文中結(jié)果可以看出,超寬帶穿墻雷達能夠有效地對墻后目標(biāo)進行成像,并且具有很好的距離向和方位向分辨率,可以更有效地探測墻后的目標(biāo)。對成像結(jié)果進行比較,可以發(fā)現(xiàn)在穿墻成像中,考慮墻體作用后,像的偏移得到了有效的抑制。被探測物體的形狀以及如何提取物體的特性參數(shù),比如:介電常數(shù),電導(dǎo)率等,將是下一步工作的研究焦點。

[1] FERRIS D D,CURRIE N C.A survey of current technologies for through-the-wall surveillance(TWS) [J].Proc SPIE,1999,3577(62).

[2] T AFLOVE A,HAGNESS S C.Computational Electrodynamics:T he Finite-Difference T ime-Domain[M]. Boston:Artech House,2000.

[3] 周璧華,陳 彬,高 成,等.鋼筋網(wǎng)及鋼筋混凝土電磁脈沖屏蔽效能研究[J].電波科學(xué)學(xué),2000,15(3):251-259. ZHOU Bihua,CHEN Bin,GAO Cheng,et al.Study on EMP shielding effectiveness of wire-mesh reinforcement and reinforced-concrete[J].Chinese Journal of Radio Science,2000,15(3):251-259.(in Chinese)

[4] 張建華,張國俊,凌永順.脈沖波反射折射特性的分析計算[J].電波科學(xué)學(xué)報,2001,16(2):200-203. ZHANG Jianhua,ZHANG Guojun,LING Yongshun. Analysis of pulsed reflection and refraction of electromagnetic waves[J].Chinese Journal of Radio Science, 2001,16(2):200-203.(in Chinese)

[5] NAG S,FLUHLER,BARNES M.Preliminary interferometric images of moving targets obtained using a time-modulated ultra-wide band through-wall penetra-tion radar[C]//Proceedings of IEEE Radar Conference,Atlanta,GA,2001,64-69.

[6] CHAMMA W,KASHYAP S."Detection of Targets Behind Walls using U ltra Wide Band Short Pulse",in book U ltra wide Band Short Pulse Electromagnetic [M].New York:Plenum Press,2003.

[7] 趙永輝,吳健生,萬明浩.不同地下介質(zhì)條件下探地雷達的探測深度問題分析[J].電波科學(xué)學(xué)報,2003(02): 220-224. ZHAO Yonghui,WU Jiansheng,WAN M inghao.T he analysis of detectable range of ground-penetrating-radar system in different underground medium[J].Ch-i nese Journal of Radio Science,2003(02):220-224.(in Chinese)

[8] AHMAD F,AMIN M G,KASSAM S A.A beamforming approach to stepped-frequency synthetic aperture through-the-wall radarimaging,[C]//Proc. First IEEE International Workshop on Computational Advances in Mult-i sensor Adaptive Processing(CAMSAP05),Puerto Vallarta,Mexico,2005.

[9] 包永芳,呂英華,張洪欣.基于FDT D方法的集成材料特性分析[J].電波科學(xué)學(xué)報,2006(4):576-581. BAO Yongfang,LV Yinghua,ZHANG Hongxin. Characteristic analysis of composite material using FDT D method[J].Chinese Journal of Radio Science, 2006(4):576-581.(in Chinese)

[10] 陳 潔,方廣有,李 芳.時域波束形成在超寬帶穿墻成像雷達中的應(yīng)用[J].電子與信息學(xué)報,2008,30 (6):1341-1344. CHEN Jie,FANG Guangyou,LI Fang.The application of time domain beamforming to ultra-wideband through-wall imaging radar[J].Journal of Electronics &Information T echnology,2008,30(6):1341-1344. (in Chinese)

Through-wall imaging ultra-wideband radar:numerical simulation

WANG Fang-fang ZHANG Ye-rong

(College of Electronics Science and Technology,N anj ing University of Post and Telecom,N anj ing Jiangsu 210003,China)

T his paper studies the capability of ultra-wideband radar to provide surveillance through concrete walls.FDT D simulator is used to generate data of high fidelity through-wall radar.The raw data are transformed into radar imaging using back projection algorithm.It is shown that ultra-wideband radar can track and pos-i tion targets inside a room with wood wall.Various effects of presence of the wall, such as refraction,change in speed,and attenuation,can defocus target images, displace targets from their true positions and possibly produce false targets.It is effective to restrain the derivation of image modifing the time in imaging algorithm by interation.

ultra-wideband radar;through-wall;finite-difference time-domain (FDT D);back projection(BP)

TN011

A

1005-0388(2010)03-0569-05

王芳芳 (1985-),女,江蘇人,南京郵電大學(xué)在讀博士生。主要研究方向為計算電磁學(xué),電磁散射,微波成像技術(shù)等。

張業(yè)榮 (1963-),男,安徽人,南京郵電大學(xué)教授、博士生導(dǎo)師。主要研究方向為移動通信系統(tǒng)與設(shè)計、電磁場的數(shù)值計算、UWB信道等。