半无限大空间中无耗二维电磁目标的微波成像新方法.docVIP

半无限大空间中无耗二维电磁目标的微波成像新方法第21卷第1期1999年1月电子科学学刊JOURNALOFELECTRONICSV0l21No1Jan.1999半无限大空间中无耗二维电磁目标的微波成像新方法,/,f/摘要本文提出利用模拟得到的时域散射信息对埋于半无限大空间中的无耗二维电磁目标进行微渡成像的新方法一时域变分迭代法(time-domainvariationaliterativemethod,TVIM).这是一种基于电磁体等效原理,变分原理和傅立叶变换的迭代方法,每次迭代各计算正问题和逆问题一次,正问题采用时域有限差分(finite—differencetime—domain)方法处理.列举了一些典型的数值反演实例,并与有关文献结果作了比较,考察了TVIM的收敛性能,反演复杂目标能力,抗随机噪声等反演性能.并从理论上简要分析了形成TVIM良好反演性能的原因.关键词微波成像,中圈号西r●■●■___--一半无限大空间,无耗二维电磁目标工言电磁逆散射(微波成像)是当今电磁学研究中应用性很强的前沿性研究课题【.地下管线等目标的探测即是一个很重要的应用场合,这一问题的物理模型可视为如图1所示的半无限大空间中二维电磁F/I标的微波成像,相比自由空间中目标的微波成像而言,这类问题的求解要困难得多.不久前,我们提出了一种方法【,初步探讨了上述问题.该方法仪能有效处理低对比度,小尺寸的弱散射F/I标,且收敛速度太慢,迭代耗时太多,在较少的迭代次数内反演精度有限.TTT丁_|0-一.R--..R—+R…R一R-0-0-0一__R0-一2基本原理一●--—+…一……●率腮自由空间0.口斛面物理模型如图1所示,其中冗表示接收点,e和e(z,Y)分别为背景媒质和未知目标的舟电常数,表示z轴方向的时域线电流源,它激励出满足如下时域标量波动方程的电场E(p,t),X7(p,t)一0E(p)a.E(pO/ot=uooJ~(o,)/.(1)1997-03?03收到.1998—05—25定稿国家自然科学基金资助项目—■■_■一:藿辩1期卿安永等;半无限太空间中无耗二维电磁目标的微波成像新方法105由体等效原理并利用E=ESC+E可得如下积分方程组删pOp,t,t景:(,t),【2),.p,(∥)蒜p,∥,(3)其中Pr与P.分别为位于自由空间中和未知目标内的场点,∥(,Y)为激励源位置,o(p)=【(p)一E1】/E0】为对比度函数,gl(,P.,亡J)和g2(P.,p,,亡J)为时域标量格林函数对(2)式和(3)式求关于o(p)的变分,有.=删()景足(p~pett)州翕()rdl=,,(4)(州,∥)景)州袅幔()(5)对(4)式和(5)式作傅立叶变换,得0=或(p.,)+.吣//G1(p.,P,)x(p)丘(p,p.)+Q(p)(p,p,w)]dp,(6)占壶:砒(pr,)=一eo//G(,pl,)×【(p)息(p,∥,u)+Q(p)最(p,p,)】dp,(7)其中G.(u)和G2(p,pl,)为频域格林函数,它们可由平面渡展开法[81推得.令0(p)=∑loR(p),其中R(p)为脉冲函数.采用点匹配法.则(,u)+..∑【Q,P,)+Qn或(p,P,)】×//G1(,P,w)dp=0.(8)(p,)=一.吣∑Qn丘(p,P,)+Qn(p,P,u)]×//G2(,pl,w)dp.(9)将(8)式和(9)式改写为矩阵方程形式为A+fJ+B).=0,(10)J=c0+面.,(11)106电子科学学刊21卷其中J为单位阵消去(10)式和(11)式中含色的项可得={c—D(J+B)一A}Sa垒埘a,(12)(12)式即为TVIM的反演方程其解为a=[Re(.]一[Re(.砒)其中上标+表示矩阵共轭转置,Re表示取实部3数值结果(13)定义如下反映每次迭代后测量得到的散射场与经迭代得到的散射场之间相对误差的函数F(0),用以控制反演进程,当以3×10时终止迭代.R(.)=ll《一碟ll./llz㈨,I~e()其中和为列向量,分别储存测量或经n次迭代后得到的时域散射场,范数E譬atI{=(c.atc,at)=(cat)cat图2收敛性能=2.Oe0.迭代初值有三种选取方法,第一种和第二种方法分别将初值选为0或,而第三种方法则是根据一些先验知识确定迭代初值,本文采用第一种方法选取迭代初值.首先考察收敛性能,反演相对介电常数为2.5,4.0和8.0的三个点目标.经4,5和l5次迭代后得到反演结果为2.5000,40000和8.0o01,在IBM350-P75上每次迭代耗时108s.相对误差值随迭代次数n的变化趋势如图2所示.显然,TVIM迭代稳定,收敛迅速,反演结果准确.下面看看TVIM反演任意复杂目标的能力.我们反演了高对比度方柱,多柱,洞及大尺寸柱等典型目