程志華謝擁軍 馬曉東 毛煜茹 畢 博

(北京航空航天大學(xué)電子信息工程學(xué)院 北京 100191)

介質(zhì)目標(biāo)的太赫茲波近場(chǎng)散射特性計(jì)算

程志華*謝擁軍 馬曉東 毛煜茹 畢 博

(北京航空航天大學(xué)電子信息工程學(xué)院 北京 100191)

該文基于廣義的Kirchhoff阻抗邊界條件和物理光學(xué)法，對(duì)太赫茲頻段介質(zhì)體近場(chǎng)散射特性進(jìn)行了研究，給出了介質(zhì)體近場(chǎng)散射計(jì)算公式。針對(duì)在太赫茲頻段由于波長(zhǎng)較短引起的計(jì)算量大幅提升的問題，采用以面元為計(jì)算單位、以像素為遮擋判斷單位的太赫茲頻段介質(zhì)體近場(chǎng)散射的快速計(jì)算方法，該方法在保證計(jì)算精度的基礎(chǔ)上，大大降低了計(jì)算復(fù)雜度和時(shí)間。計(jì)算了圓柱體和鴨嘴形介質(zhì)體在不同距離下的雷達(dá)散射截面，并且分析了電磁場(chǎng)與物體相互作用后，相位項(xiàng)在不同距離、不同頻率下對(duì)介質(zhì)體雷達(dá)散射截面的影響。

雷達(dá)散射截面；圖形電磁學(xué)；太赫茲；近場(chǎng)

1 引言

太赫茲(THz)波通常是指介于0.1～10.0 THz之間的電磁波，在電磁波譜中太赫茲波處于紅外和微波之間。由于太赫茲頻段的空間分辨率和時(shí)間分辨率頻率都很高，與微波和X射線成像相比，太赫茲成像不僅可以顯示更精密的圖像，還可以獲得頻域內(nèi)的相位信息，因此太赫茲成像技術(shù)是太赫茲應(yīng)用的關(guān)鍵技術(shù)。太赫茲成像是一個(gè)典型的電磁逆散射問題，是通過接收目標(biāo)散射數(shù)據(jù)重建描寫介質(zhì)或目標(biāo)特性的基本函數(shù)，同時(shí)，由于太赫茲頻段頻率較高，近場(chǎng)距離范圍大幅度提升，成像中散射模型應(yīng)為近場(chǎng)散射計(jì)算模型。因此，太赫茲頻段介質(zhì)目標(biāo)近場(chǎng)散射特性的研究是十分重要和有意義的。

太赫茲頻段目標(biāo)散射特性研究可分為實(shí)驗(yàn)測(cè)量、數(shù)值計(jì)算兩種手段。在實(shí)驗(yàn)測(cè)量上，美國(guó)麻省LOWELL大學(xué)毫米波實(shí)驗(yàn)室利用1.56 THz源在緊縮場(chǎng)中對(duì)粗糙面圓柱體的目標(biāo)散射特性進(jìn)行了研究[1]。新澤西理工學(xué)院物理系對(duì)太赫茲頻段危險(xiǎn)物品的散射成像進(jìn)行了研究[2]。天津大學(xué)太赫茲研究中心搭建了以0.2 THz 返波管振蕩器源、熱釋電探測(cè)器、小型自動(dòng)旋轉(zhuǎn)光學(xué)平臺(tái)等組成的太赫茲波目標(biāo)散射特性實(shí)驗(yàn)測(cè)試系統(tǒng)，并對(duì)粗糙銅面的散射特性等進(jìn)行了研究[3]。首都師范大學(xué)太赫茲實(shí)驗(yàn)室研制了太赫茲數(shù)字全息成像系統(tǒng)，該系統(tǒng)可以對(duì)太赫茲電磁波的振幅、相位、頻率及偏振等全部光學(xué)信息的3維空間分布進(jìn)行精確測(cè)量[4]。

在數(shù)值計(jì)算領(lǐng)域，東南大學(xué)毫米波實(shí)驗(yàn)室對(duì)太赫茲波頻段的金屬目標(biāo)體的遠(yuǎn)場(chǎng)目標(biāo)特性計(jì)算進(jìn)行了研究，對(duì)比分析了多種電磁計(jì)算方法，并且獨(dú)立自主地開發(fā)了準(zhǔn)確高效的目標(biāo)特性和電磁兼容仿真分析軟件[5]。本課題組對(duì)太赫茲頻段金屬目標(biāo)的近場(chǎng)散射特性計(jì)算進(jìn)行了研究，并取得了一定成果。介質(zhì)體近場(chǎng)散射特性的研究較少見報(bào)道，遠(yuǎn)場(chǎng)條件下，目標(biāo)的散射特性已得到廣泛的研究[6?11]，介質(zhì)目標(biāo)體計(jì)算方法可以分為3類：一是以等效流為基礎(chǔ)的方法[12]；二是近似邊界條件法(ABC)[13]；三是數(shù)值方法，如基于積分方程的矩量法、快速多極子法、體積分方程法等，基于微分方程的有限元法、時(shí)域有限差分法等，高頻方法有射線追蹤法(SBR)、幾何光學(xué)法(GO)、物理光學(xué)法(PO)、幾何繞射理論(GTD)等。基于上述方法，近場(chǎng)條件下，目標(biāo)體RCS的計(jì)算思路主要有兩種：一種是將物體表面剖分成小面元，用遠(yuǎn)場(chǎng)RCS計(jì)算方法分別計(jì)算各小面元的RCS值，然后通過相位疊加得到總的散射場(chǎng)，是從部分到整體的思路[14]。另外一種是以物理光學(xué)法為基礎(chǔ)，通過散射場(chǎng)積分方程在無遠(yuǎn)場(chǎng)近似的情況下推導(dǎo)得出的方法，是從整體到部分的思路[15]。在近場(chǎng)區(qū)域，不同距離下的近場(chǎng)散射分析更具實(shí)際意義，目標(biāo)體照射區(qū)域?qū)㈦S著距離的不同而變化，隨之相互遮擋情況也會(huì)產(chǎn)生變化，從而導(dǎo)致大量的時(shí)間用在相互遮擋的判斷上，因此，采用上述的高頻方法仍需要耗費(fèi)極大的內(nèi)存和極長(zhǎng)的計(jì)算時(shí)間。1993年RIUS等人[16]提出了基于PO的圖形電磁學(xué)(GRECO)方法，該方法充分利用了計(jì)算機(jī)硬件的優(yōu)勢(shì)，由計(jì)算機(jī)硬件完成最困難、最費(fèi)時(shí)的遮擋和消影工作，這使得它具有很大的工程應(yīng)用價(jià)值，但是GRECO方法只適用于遠(yuǎn)場(chǎng)區(qū)域的電磁散射計(jì)算，并且GRECO方法以提取像素顏色強(qiáng)度獲得參數(shù)的本質(zhì)決定其無法適用于太赫茲頻段近場(chǎng)的計(jì)算。如何把GRECO方法在遮擋判斷上的優(yōu)勢(shì)應(yīng)用在太赫茲頻段近場(chǎng)散射計(jì)算中是十分有意義和價(jià)值的。

本文利用基于廣義的Kirchhoff的阻抗邊界條件對(duì)介質(zhì)體近場(chǎng)RCS計(jì)算公式進(jìn)行了推導(dǎo)?；赑O和GRECO方法，提出一種以面元為計(jì)算單位、以像素為遮擋判斷單位的太赫茲波介質(zhì)體近場(chǎng)RCS的快速計(jì)算方法。該方法在保證計(jì)算精度的基礎(chǔ)上，大大降低了遮擋判斷的計(jì)算復(fù)雜度和時(shí)間。以圓柱體和鴨嘴形介質(zhì)體在不同距離下的雷達(dá)散射截面的計(jì)算為例，驗(yàn)證了本文方法的有效性和準(zhǔn)確性。

2 理論分析

2.1 介電性質(zhì)

德拜方程是討論一般介質(zhì)極化弛豫特性的重要關(guān)系式，它給出了介電常數(shù)和頻率的關(guān)系，可表示為

其中ε(0)為靜態(tài)介電常數(shù)，ε(∞)為高頻節(jié)點(diǎn)常數(shù)，τD為弛豫時(shí)間，當(dāng)介質(zhì)體的弛豫時(shí)間的倒數(shù)遠(yuǎn)大于電磁波頻率時(shí)，由式(1)可知，介質(zhì)介電常數(shù)可近似認(rèn)為是ε(0)。太赫茲頻段，頻率在1011～1013Hz 之間，而大部分介質(zhì)體的弛豫時(shí)間大于皮秒級(jí)，因此在太赫茲頻段，大部分介質(zhì)體的介電常數(shù)將是與介質(zhì)體弛豫時(shí)間、高頻介電常數(shù)和靜態(tài)介電常數(shù)相關(guān)的變量。材料的弛豫時(shí)間和介電常數(shù)無法通過理論推導(dǎo)求出，需要通過實(shí)驗(yàn)獲得，目前針對(duì)太赫茲頻段的物質(zhì)特性的研究大都通過對(duì)吸收系數(shù)和折射系數(shù)展開，通過折射系數(shù)和吸收系數(shù)可以求得介質(zhì)的介電常數(shù)。

2.2 近場(chǎng)介質(zhì)體散射計(jì)算理論

由Stratton-Chu積分公式，當(dāng)散射體內(nèi)無電磁源時(shí)，無遠(yuǎn)場(chǎng)近似的散射場(chǎng)為

其中Es為目標(biāo)體散射電場(chǎng)，積分面S為光照區(qū)表面，ds為表面單位單元，其單位法向矢量為n, s為散射場(chǎng)方向單位矢量，J為電流密度，E和H分別為總電場(chǎng)和總磁場(chǎng)，Z0為波阻抗，φ=ejkr/r, φ1= ejkr/r2, r為源與面元ds的距離，其單位矢量為r。假設(shè)物體表面的曲率半徑遠(yuǎn)大于波長(zhǎng)，考慮接收單位的極化方向，由散射場(chǎng)積分方程可得

其中i為入射場(chǎng)方向單位矢量，iH為入射磁場(chǎng)。

由于是球面波入射，因此，

其中hi為入射磁場(chǎng)方向單位矢量，由式(3)，式(4)可得

近場(chǎng)雷達(dá)散射截面公式可定義為

其中R為源點(diǎn)到目標(biāo)體中心的距離。

利用廣義的Kirchhoff的阻抗邊界條件，介質(zhì)體表面電流和磁流可以表示為

由式(5)，式(6)結(jié)合式(7)，式(8)化簡(jiǎn)后得到介質(zhì)體近場(chǎng)RCS的計(jì)算表達(dá)式為

R⊥,R//分別為介質(zhì)表面在垂直極化和水平極化時(shí)的反射系數(shù)，對(duì)非磁性介質(zhì)進(jìn)行考慮，μ1=μ2=μ0，其可表示為

其中ε1, ε2分別為空氣和媒質(zhì)介電系數(shù)。

3 太赫茲頻段近場(chǎng)散射快速計(jì)算方法

由式(9)，式(10)可知，該計(jì)算公式是以面元為計(jì)算單位。在太赫茲頻段，為保證計(jì)算結(jié)果的準(zhǔn)確性，剖分面元將會(huì)是千萬級(jí)的數(shù)量，甚至更高，以面元為遮擋判斷單位計(jì)算量十分巨大，因此，物理光學(xué)法在太赫茲頻段的RCS計(jì)算中工程應(yīng)用性很低。以像素為計(jì)算單位的圖形電磁學(xué)方法，在判斷遮擋上具有很大的優(yōu)勢(shì)，但是，圖形電磁學(xué)中，計(jì)算所需參數(shù)是以像素為單位通過對(duì)光強(qiáng)值的提取來得到的，一般情況下顏色分配模式為RGBA模式，硬件為R, G, B, A成分保留一定數(shù)量的位平面。目前顯卡支持的最高顏色位數(shù)為32位，R, G, B和A各占8位，8 位的R成分可以存儲(chǔ)0～255之間的整數(shù)[17]。位平面中的0,1,…,255將對(duì)應(yīng)于0/255=0.0, 1/255,…,255/255=1的顏色值，因此，通過光照提取出來的參數(shù)值都將是1/255=0.0039的倍數(shù)，而太赫茲頻段頻率范圍在0.1～10.0 THz之間，最低頻率0.1 THz時(shí)，波長(zhǎng)為0.003 m，對(duì)比通過光照提取的參數(shù)值的分辨率(1/255)，可知以像素為計(jì)算單位，通過提取像素顏色強(qiáng)度來獲得所需參數(shù)的方法無法滿足太赫茲頻段的計(jì)算需求。同時(shí)，在近場(chǎng)的計(jì)算中，由于是球面波照射，照射區(qū)域和相位因素不能通過平面波光源和深度緩存得到，因此傳統(tǒng)的圖形電磁學(xué)中提取參數(shù)方法和遮擋判斷方法在太赫茲目標(biāo)散射的計(jì)算中完全失效。本文通過以面元為計(jì)算單位保證太赫茲頻段散射特性計(jì)算精度，以像素為遮擋判斷單位保證遮擋判斷的快速性，最終通過建立像素與面元的對(duì)應(yīng)關(guān)系從而達(dá)到快速計(jì)算太赫茲目標(biāo)近場(chǎng)散射值的目的。由式(9)，式(10)可知，所需計(jì)算的參數(shù)包括：各面元到光源的距離，各面元法向與入射波夾角的余弦值。

按照目標(biāo)體實(shí)際尺寸，用FEMAP對(duì)其進(jìn)行建模并進(jìn)行剖分，剖分文件中包括各面元的法相矢量和位置信息。每一個(gè)面元的信息包括3個(gè)頂點(diǎn)的坐標(biāo)(x1,y1,z1),(x2,y2,z2),(x3,y3,z3)和面元的法相矢量n，點(diǎn)源的坐標(biāo)(x,y,z)可以通過光源的位置設(shè)置得到。根據(jù)坐標(biāo)和法向可以得到各面元到點(diǎn)源的距離和入射波與面元法向的夾角余弦值。

在近場(chǎng)條件下，照射區(qū)域?qū)㈦S距離的變化而變化，隨著遮擋情況也會(huì)產(chǎn)生變化。通過OPENGL將光源設(shè)置為位置性光源，光源位置信息通過模型視圖矩陣進(jìn)行變化，并以視覺坐標(biāo)的形式存儲(chǔ)。光照區(qū)域?qū)㈦S光源位置的不同自動(dòng)進(jìn)行光照區(qū)域判斷，結(jié)合深度緩存自動(dòng)進(jìn)行遮擋處理，從而得到準(zhǔn)確的光照區(qū)域。利用OPENGL將消隱后的目標(biāo)圖像顯現(xiàn)在計(jì)算機(jī)屏幕上，如圖1所示。

圖1 隨距離變化的光照區(qū)域變化

光照區(qū)域的面元顏色是以像素為單位進(jìn)行存儲(chǔ)的，無法直接判斷面元是否處于光照區(qū)域。通過OPENGL工具函數(shù)庫(kù)中的gluUnproject函數(shù)將屏幕上光照區(qū)像素與剖分體面元一一對(duì)應(yīng)，從而以像素為光照區(qū)域判斷單位，當(dāng)像素有顏色時(shí)，由對(duì)應(yīng)的面元進(jìn)行計(jì)算。在計(jì)算過程中會(huì)出現(xiàn)單個(gè)像素對(duì)應(yīng)多個(gè)面元的情況，在這種情況下通過利用OPENGL的深度緩存判定離剪裁面最近的面元為計(jì)算單位，放棄其他面元的計(jì)算。為了防止出現(xiàn)單個(gè)像素對(duì)應(yīng)多個(gè)面元的情況，采用遠(yuǎn)場(chǎng)散射計(jì)算中的分屏顯示方法，保證像素對(duì)應(yīng)的實(shí)際的尺寸小于或等于面元剖分精度。

以半徑0.254 m球體為例，入射頻率300 GHz，剖分精度為1/3波長(zhǎng)(m)，所需要的分屏數(shù)為

其中max(L), max(W)分別為目標(biāo)體最大長(zhǎng)度和最大寬度；pl, pw分別別為顯示區(qū)域像素的長(zhǎng)度和寬度。該方法的計(jì)算量為

射線法是通過判斷各剖分面元中點(diǎn)與源點(diǎn)形成的射線與其他面元是否相交，然后保留與源點(diǎn)最近的面元進(jìn)行計(jì)算。DirectX中采用的就是此方法，其具體計(jì)算方法是：通過求解三角形面元中兩條邊的權(quán)重u,v，如果u+v≤1則判定射線與該面元存在交點(diǎn)，其計(jì)算量為：n(n?1)×[N(u)+M(v)+1]，其中n為面元數(shù)，N(u), M(v)分別為u,v的計(jì)算量。

表1為同樣以半徑0.254 m球體為目標(biāo)體，入射頻率為300 GHz，剖分精度為0.0003 m時(shí)，采用射線法判斷遮擋與本文方法的計(jì)算量對(duì)比。從表1可以看出，采用本文方法大幅度減少了計(jì)算量。

表1 射線法與本文方法遮擋判斷的計(jì)算量對(duì)比

4 計(jì)算結(jié)果及分析

為驗(yàn)證本文方法的有效性，計(jì)算了邊長(zhǎng)為0.3 m的介質(zhì)立方體在距離為0.04 m時(shí)的近場(chǎng)RCS，入射波長(zhǎng)為1 m, εr=4, ur=1，入射方向?yàn)榕c立方體表面法向夾角5°～40°的范圍，計(jì)算間隔5°。文獻(xiàn)[18]中采用矩量法(MOM)對(duì)其進(jìn)行了計(jì)算，兩者計(jì)算結(jié)果如圖2所示，從中可以看出當(dāng)角度較小時(shí)與文獻(xiàn)[18]中計(jì)算結(jié)果與本文計(jì)算結(jié)果吻合較好，當(dāng)角度增大時(shí)誤差略微增大(由棱邊引起)，最大誤差小于1 dB，因此驗(yàn)證了本文方法的正確性和有效性。

圖2 介質(zhì)立方體近場(chǎng)RCS本文計(jì)算結(jié)果與文獻(xiàn)計(jì)算結(jié)果對(duì)比

以介質(zhì)體為含碳納米結(jié)構(gòu)的高分子復(fù)合材料為例[19]，分別計(jì)算了圓柱體和同比尺寸縮小后的鴨嘴型散射體，計(jì)算頻率為100 GHz, 200 GHz和300 GHz，其介電常數(shù)如表2所示。

表2 含碳納米結(jié)構(gòu)的高分子復(fù)合材料介電常數(shù)表

(1)介質(zhì)圓柱體計(jì)算距離為10～200 m，剖分精度為0.0005 m，正Z軸入射，垂直極化，模型如圖3所示，計(jì)算結(jié)果如圖4所示。

圖3 圓柱介質(zhì)體模型

從圖4可以看出當(dāng)目標(biāo)體與源距離較近時(shí)，同一距離下，不同頻率的近場(chǎng)RCS值并沒有像遠(yuǎn)場(chǎng)RCS一樣隨頻率的增大而增大，隨著距離的增加，頻率與近場(chǎng)RCS值逐漸成單調(diào)趨勢(shì)。圖5為固定介電常數(shù)為12.45，圓柱介質(zhì)體在頻率1～200 GHz，距離分別為10 m, 50 m和200 m時(shí)的近場(chǎng)RCS值?？梢钥闯觯?dāng)距離較近時(shí)，近場(chǎng)RCS值隨頻率的增大波動(dòng)較強(qiáng)；當(dāng)距離較遠(yuǎn)時(shí)，近場(chǎng)RCS值隨頻率的增加波動(dòng)降低。由式(16)，式(17)可知，當(dāng)距離不變、頻率發(fā)生變化時(shí)，只有系數(shù)項(xiàng)4π/λ2與相位項(xiàng)e2jkr(其中k=2π/λ)發(fā)生變化，其中系數(shù)項(xiàng)與頻率成單調(diào)正比趨勢(shì)，結(jié)合圖5中的計(jì)算結(jié)果，可以認(rèn)為在近場(chǎng)散射計(jì)算中，當(dāng)距離越近時(shí)，相位項(xiàng)對(duì)近場(chǎng)散射值的影響越大；隨著距離的增加，相位項(xiàng)對(duì)近場(chǎng)散射值的影響越來越小，到達(dá)遠(yuǎn)場(chǎng)條件時(shí)，散射值與頻率成正比，本文計(jì)算結(jié)果與遠(yuǎn)場(chǎng)計(jì)算結(jié)果吻合。

(2)同比尺寸縮小后的鴨嘴型散射體，計(jì)算距離為1～50 m，正Z軸入射，水平極化，剖分精度0.0005 m，模型如圖6所示，計(jì)算結(jié)果如圖7所示。

圖4 圓柱介質(zhì)體10～200 m近場(chǎng)RCS

圖5 1～200 GHz固定距離圓柱介質(zhì)體近場(chǎng)RCS

圖6 鴨嘴型天線罩介質(zhì)體模型

圖7 鴨嘴型散射體不同頻率近場(chǎng)RCS

從圖7中同樣可以看出，當(dāng)距離越遠(yuǎn)時(shí)相位項(xiàng)對(duì)近場(chǎng)RCS值的影響越小，在近距離一定范圍內(nèi)，近場(chǎng)RCS值并不隨頻率增加而增大，當(dāng)?shù)竭_(dá)遠(yuǎn)場(chǎng)條件時(shí)與遠(yuǎn)場(chǎng)計(jì)算結(jié)果吻合。

5 結(jié)論

本文研究了太赫茲頻段介質(zhì)目標(biāo)體的近場(chǎng)散射特性，基于廣義的Kirchhoff阻抗邊界條件和物理光學(xué)法對(duì)介質(zhì)體近場(chǎng)散射計(jì)算公式進(jìn)行了推導(dǎo)，采用以面元為計(jì)算單位，以像素為遮擋判斷單位的太赫茲近場(chǎng)散射的快速計(jì)算方法。該方法利用OPENGL引入點(diǎn)光源和深度緩存結(jié)合的光照區(qū)域判斷方法，由計(jì)算機(jī)硬件自動(dòng)完成耗時(shí)的遮擋和消隱以及建立面元與像素之間的對(duì)應(yīng)關(guān)系，從而快速準(zhǔn)確地計(jì)算了太赫茲頻段介質(zhì)目標(biāo)的近場(chǎng)雷達(dá)散射截面，并且分析了電磁場(chǎng)與物體相互作用后，相位項(xiàng)在不同距離、不同頻率下對(duì)介質(zhì)體雷達(dá)散射截面的影響。

[1] Danylov A A, Goyette T M, Waldman J, et al.. Terahertz inverse synthetic aperture radar (ISAR) imaging with a quantum cascade laser transmitter[J]. Optics Express, 2010, 18(15): 16264-16272.

[2] Federici F, Schulkin B, Huang F, et al.. THz imaging and sensing for security applications-explosives, weapons and drugs[J]. Semiconductor Science and Technology, 2005, 20(2): S266-S280.

[3] 楊洋, 姚建銓, 張鏡水, 等. 粗糙銅表面對(duì)低頻太赫茲波的散射實(shí)驗(yàn)[J]. 紅外與毫米波學(xué)報(bào), 2013, 32(1): 36-39.

Yang Yang, Yao Jian-quan, Zhang Jing-shui, et al.. Terahertz scattering on rough copper surface[J]. Journal of Infrared and Millimeter Waves, 2013, 32(1): 36-39.

[4] 鄭顯華, 王新柯, 孫文峰, 等. 太赫茲數(shù)字全息術(shù)的研發(fā)與應(yīng)用[J]. 中國(guó)激光, 2014, 41(2): 1-10.

Zheng Xian-hua, Wang Xin-ke, Sun Wen-feng, et al.. Developments and applications of the terahertz digital holography[J]. Chinese Journal of Lasers, 2014, 41(2): 1-10.

[5] Li Zhuo, Cui Tie-jun, Zhong Xing-jian, et al.. Electromagnetic scattering characteristics of PEC targets in the terahertz regime[J]. IEEE Antennas and Propagation Magazine, 2009, 51(1): 39-50.

[6] Solomon R C, Leong H, and Antar Y M. Forward scattering effects in RCS of complex targets in the 3-20 MHz high frequency range[C]. Antennas and Propagation Society International Symposium, Boston, 2008: 1-4.

[7] 李曉峰, 謝擁軍, 楊銳. 半空間電大導(dǎo)體目標(biāo)散射的高頻分析方法[J]. 電子與信息學(xué)報(bào), 2009, 31(5): 1268-1270.

Li Xiao-feng, Xie Yong-jun, and Yang Rui. High-frequency analysis on scattering from conductive targets with electrically large size in half space[J]. Journal of Electronics & Information Technology, 2009, 31(5): 1268-1270.

[8] Gillion E, Rochefort E, and Claverie J. Improvement of RCS estimation of large targets by using near-field approach[C]. IEEE Radar Conference (RADAR), Ottawa, ON, Canada, 2013: 1-5.

[9] Altin N and Yazgan E. RCS prediction using fast ray tracing in Plücker coordinates[C]. 2013 7th European Conference on Antennas and Propagation (EuCAP), Gothenburg, 2013: 284-288.

[10] Lipuma D, Me?ric S, and Gillard R. RCS enhancement of flattened dihedral corner reflector using reflect array approach[J]. Electronics Letters, 2013, 49(2): 152-154.

[11] Li Xiao-feng, Xie Yong-jun, and Yang Rui. High-frequency method for scattering from coated targets with electrically large size in half space[J]. IEEE Transactions on Antennas and Propagation, 2009, 57(2): 181-186.

[12] Huddleston P L, Medgyesi-Mitschang L N, and Putnam J M. Combined field integral equation formulation for scattering by dielectrically coated conducting bodies[J]. IEEE Transactions on Antennas and Propagation, 1986, 34(4): 510-520.

[13] Senior T B A and Volakis J L. Generalized impedance boundary condition in scattering[J]. Proceedings of the IEEE, 1991, 79(10): 1413-1420.

[14] Li Jing and Li Xiang-jun. An efficient algorithm for near field RCS of electrically large dynamic targets[C]. International Conference on Information Science and Technology, Nanjing, 2011: 413-415.

[15] Pouliguen P, Damiens J F, et al.. RCS computation in near field[C]. Proceedings of the International Conference Day on Diffraction, Petersburg, 2006: 252-265.

[16] Rius J M, Ferrando M, and Jofre L. High frequency RCS of complex radar targets in real time[J]. IEEE Transactions on Antennas and Propagation, 1993, 41(9): 1308-1319.

[17] Shreniner D. The Khronos OpenGL ARB Working Group.

李軍, 徐波, 等譯. OPENGL編程指南[M]. 第7版, 北京: 機(jī)械工業(yè)出版社, 2010: 118-119.

[18] 賀昌輝. 矩量法分析目標(biāo)的近場(chǎng)電磁散射[D]. [碩士論文], 華中師范大學(xué), 2006.

He Chang-hui. The analysis of near-field EM scattering from targets by moment methods[D]. [Master dissertation], Central China Normal University, 2006.

[19] Zhuravleva V A, Suslyaeva V I, and Dunaevskiia G E. Complex permittivity of polymer composites containing carbon nanostructures in frequency range 0.17-1.1[C]. International Conference on THz Infrared, Millimeter and Terahertz Waves (IRMMW-THz), Wollongong, NSW, 2012: 1-2.

程志華： 男，1982年生，博士生，研究方向?yàn)橛?jì)算電磁學(xué)、太赫茲目標(biāo)特性、天線工程、電磁兼容等.

謝擁軍： 男，1968年生，教授，研究方向?yàn)橛?jì)算微波與計(jì)算電磁學(xué)、微波通信、天線工程、電磁兼容等.

毛煜茹： 女，1984年生，博士生，研究方向?yàn)橛?jì)算電磁學(xué)、天線工程、電磁兼容等.

馬曉東： 男， 1989年生，碩士生，研究方向?yàn)槟繕?biāo)特性、天線工程、電磁兼容等.

畢 博： 男，1993年生，碩士生，研究方向?yàn)槟繕?biāo)特性等.

Near-field Electromagnetic Scattering Characteristics of Dielectric Targets in the Terahertz Regime

Cheng Zhi-hua Xie Yong-jun Ma Xiao-dong Mao Yu-ru Bi Bo
(School of Electronic and Information Engineering, Beihang University, Beijing 100191, China)

The near-field scattering characteristics of dielectric targets are studied in the terahertz band based on the generalized Kirchhoff impedance boundary conditions and the physical optics method. The formula of the near field scattering for dielectric targets is deduced. In the light of the increase in calculation amount caused by the shorter wavelength, a fast computational method using surface element as the calculating unit and pixel as occlusion judgment unit is proposed for the near field scattering computation in terahertz band. The method ensures the calculation accuracy and reduces greatly the computational complexity and the time consumption of occlusion judging. The calculation of the near-field Radar Cross Section (RCS) produced by a dielectric cylinder and a duck mouth scatterer is performed in terahertz band. Meanwhile, the effect of the phase on the near field RCS in different distances and frequency is analyzed.

Radar Cross Section (RCS); GRaphical-Electromagnetic COmputing (GRECO); Terahertz; Near field

TN011

： A

：1009-5896(2015)04-1002-06

10.11999/JEIT140807

2014-06-20收到，2014-08-15改回

*通信作者：程志華 asakaka1982@163.com