一種雷達目標的動態(tài)雷達散射截面模擬

劉 哲 韓慧奇

(中國人民解放軍63726部隊，寧夏 銀川 750004)

0 引言

隨著隱身與反隱身技術(shù)的發(fā)展，研究目標的隱身與反隱身性能已經(jīng)成為測控領(lǐng)域必不可少的一項工作，而研究目標的雷達散射截面(radar cross section，RCS)，特別是目標的動態(tài)RCS就顯得尤為重要。

目前，目標RCS的獲取主要依靠理論計算、外場測試、暗室試驗以及電磁仿真軟件仿真等方法［1－2］。只有外場試驗能測量目標的動態(tài)RCS，但外場試驗所需條件苛刻，實施較為困難。文獻［1］在真實環(huán)境下測量目標的動態(tài)RCS，然后利用RCS起伏模型得到RCS的預處理結(jié)果。利用起伏模型計算具有一定的不確定性。文獻［2］基于暗室縮比模型條件下測得目標RCS的靜態(tài)值，然后根據(jù)目標航跡得到動態(tài)RCS，其中的縮比試驗需要在特定的暗室進行，并且一次試驗只能對一個目標進行，具有一定的局限性。

基于以上考慮，需要尋找一種較為通用、易于實現(xiàn)且具有一定精度的RCS預處理方法，同時需要實現(xiàn)目標動態(tài)RCS的顯示。雷達目標的動態(tài)RCS模擬對于研究和掌握目標的飛行特性起到了一定的作用。

1 幾種相關(guān)坐標系

1.1 測量坐標系和發(fā)射坐標系

在目標跟蹤的過程中，雷達工作頻率較高，目標波長一般遠大于雷達工作波長，在計算目標RCS時可以采用數(shù)值分析的高效算法。實現(xiàn)目標動態(tài)RCS模擬的關(guān)鍵是處理好測量坐標系和目標坐標系之間的轉(zhuǎn)換關(guān)系。在目標建模中合理地選取建模坐標系，以達到簡化坐標轉(zhuǎn)換的目的。

在測控領(lǐng)域中，測量坐標系是站心坐標系的一種。它是以測控設(shè)備中心為原點的坐標系，這里我們以O(shè)cl-xclyclzcl來表示測量坐標系，其中:Ocl位于天線回轉(zhuǎn)中心;ycl軸與過原點的地球參考橢球面方向一致，并指向橢球面外;xcl軸與ycl軸垂直，由原點指向大地北;zcl軸與xcl軸和ycl軸構(gòu)成右手坐標系。

發(fā)射坐標系也是站心坐標系的一種，用Ofs-xfsyfszfs來表示，其中:Ofs是目標質(zhì)心在發(fā)射平面的投影點;yfs軸與發(fā)射點的鉛垂線一致，指向向上;xfs軸與yfs軸垂直，由原點指向發(fā)射瞄準方向;zfs軸、xfs軸和yfs軸構(gòu)成右手坐標系。

1.2 目標坐標系

本文采用O-xMByMBzMB來表示目標坐標系，其中:O位于目標質(zhì)心;xMB軸與目標縱向重合，并指向目標頭部;yMB軸在目標縱對稱面內(nèi)，垂直于xMB軸，當目標平飛時指向上方;zMB軸、xMB和yMB軸構(gòu)成右手系。目標坐標系描述了空間點對目標的位置關(guān)系，除了采用直角坐標系(xMB、yMB、zMB)表示外，還可以用更加直觀的極坐標表示。

本文采用極坐標表示法，其示意圖如圖1所示。

圖1 目標坐標系的極坐標表示示意圖Fig.1 Polar coordinate representation of the target coordinates

圖1中，αB為目標徑向在O-yMBzMB平面的投影OP'與yMB軸的夾角，方向由 yMB軸向 zMB軸為正;βB為目標徑向與xMB軸的夾角，方向為從xMB軸向目標尾部轉(zhuǎn)動。

1.3 目標姿態(tài)的表示

目標姿態(tài)由目標坐標系相對于發(fā)射坐標系的轉(zhuǎn)動角描述，包括俯仰角φ、偏航角ψ和滾動角γ。俯仰角φ是指目標坐標軸縱軸在Ofs-xfsyfs平面上的投影與xfs軸的夾角，方向由xfs轉(zhuǎn)向yfs;偏航角ψ是指xMB與Ofsxfsyfs平面夾角，方向從Ofs-xfsyfs平面逆時針旋轉(zhuǎn);滾動角γ是指O-xMByMB面繞xMB軸的轉(zhuǎn)角，方向順xMB軸方向右為正。其中，下標MB代表目標坐標系，下標fs表示發(fā)射坐標系。

為便于直觀理解，兩個坐標系和俯仰角φ、偏航角ψ、滾動角γ之間的關(guān)系如圖2所示。

圖2 目標坐標系的姿態(tài)角定義Fig.2 Definition of posture angle in target coordinates

2 坐標轉(zhuǎn)換

2.1 坐標系的轉(zhuǎn)換

為了計算測量坐標系在目標坐標系中的位置，先將其轉(zhuǎn)換到發(fā)射坐標系下［3］，即:

其中，i分別代表下標f(發(fā)射坐標系)或c(測量坐標系);a為參考橢球長半軸;e為第一偏心率;Lf、Bf、hf、Lc、Bc、hc分別為發(fā)射坐標系原點和測量坐標系原點的大地經(jīng)緯度和高程;Af為目標姿態(tài)的發(fā)射方位角。

2.2 直角坐標表示

將發(fā)射坐標表示轉(zhuǎn)換為直角坐標表示的計算公式為:

2.3 極坐標表示

測量坐標系的極坐標示意圖如圖3所示。

圖3 測量坐標系的極坐標表示示意圖Fig.3 Polar coordinate representation of the measurement coordinates

將直角坐標表示轉(zhuǎn)換為極坐標表示的計算公式為:

3 快速多極子算法原理

快速多極子方法是將積分方程中的格林函數(shù)進行加法處理，在角譜空間中用平面波進行算子對角化，將計算問題轉(zhuǎn)化為稀疏陣與矢量的相乘。算法的核心思想是分組。

任意兩個子散射體間的互耦根據(jù)所在組的位置關(guān)系，采用不同的算法來減少計算量。相鄰組直接進行數(shù)值計算，非相鄰組采用聚合、轉(zhuǎn)移、配置的步驟計算。

對于場點組，先將各個非相鄰組的所有子散射體產(chǎn)生的貢獻聚合到組中心，再將這些貢獻轉(zhuǎn)移到給定場點組的中心，然后將散射體產(chǎn)生的貢獻聚合到組中心，接著將這些貢獻轉(zhuǎn)移到給定場點組的中心，最后將由中心配置到組內(nèi)各個子散射體。

若散射體表面有N個子散射體，當直接計算互耦時計算量為O(N2)，而采用快速多極子計算量只為O(N1.5)，大大減少了計算量［4］。

矩陣與矢量相乘的快速多極子表達公式為［5］:

式中:j∈Gm;αi為第i個源子散射體的電流幅度;NG為附近組的貢獻;FG為非附近組的貢獻;Vsm'i(k^)、Vfmj(k^)αmm'為聚合因子、轉(zhuǎn)移配置因子。

對于電大尺寸目標的散射問題，其未知量數(shù)目N＞＞1，應用多層快速多極子法(multilevel fast multipole method，MLFMM)可以獲得更高的效率。基于多層級結(jié)構(gòu)的思想，提出了多層快速多極子方法。對于N體互耦，多層快速多極子方法采用多層分區(qū)計算，附近區(qū)強耦合量直接計算，非附近區(qū)耦合量用多層快速多極子方法實現(xiàn)。

多層分區(qū)計算的特點是逐層聚合、逐層轉(zhuǎn)移、逐層配置、嵌套遞推。從實質(zhì)上講，多層快速多極子方法是快速多極子方法的多層擴展，多層快速多極子方法除了聚合、轉(zhuǎn)移、配置的相似步驟外，還有父層、子層間的遞推計算［6－7］。

采用多層快速多極子法，對文獻［7］中的某椎柱幾何模型進行建模仿真，模型坐標系z軸正向指向椎尖，圓柱半徑為80 mm、柱高260 mm、椎半角為11°、入射頻率為5 GHz、入射角度 φ=90°，得到單站RCS的變化曲線如圖4所示。

圖4 RCS變化曲線Fig.4 RCS variation curve

由圖4可以看出，RCS曲線與文獻［7］中驗證的結(jié)果基本吻合，說明此方法的有效性和正確性，滿足工程實際需要。

4 動態(tài)RCS仿真程序流程

目標動態(tài)RCS仿真程序流程如圖5所示。

圖5 目標動態(tài)RCS仿真程序流程圖Fig.5 Flowchart of the target dynamic RCS simulation program

根據(jù)以上分析，本文中實現(xiàn)目標動態(tài)RCS模擬的具體步驟如下。

①利用幾種坐標變換原理，編程實現(xiàn)初始航跡對應的目標坐標系下的(αc，βc)數(shù)組對，并形成數(shù)據(jù)庫等待調(diào)用;

②根據(jù)多層快速多極子算法，對目標在建模的基礎(chǔ)上進行遠場仿真，建模坐標系的建立與目標系一致，由此得到與(αc，βc)數(shù)組對具有一一對應關(guān)系的RCS數(shù)據(jù)庫文件;

③按照目標航跡的飛行時間將目標仿真庫的RCS值以動態(tài)時間序列顯示，模擬目標動態(tài)RCS。

根據(jù)以上步驟得到文獻［7］中的某椎柱幾何模型單站動態(tài)RCS模擬曲線如圖6所示。

圖6 目標單站動態(tài)RCS模擬曲線圖Fig.6 Curve of target single station dynamic RCS simulation

從圖6可以看出，400序列RCS值達到最大，525序列以后有幾個波動較大的峰值，并且RCS值在這段序列內(nèi)變化劇烈。這也恰好說明這段時間目標姿態(tài)調(diào)整頻繁，與航跡數(shù)據(jù)分析基本一致。從工程角度講，這種模擬具有大致掌握目標飛行特性的實際意義［8］。

5 結(jié)束語

本文采用多層快速多極子算法，通過對坐標系的恰當選取和坐標系轉(zhuǎn)換，編程實現(xiàn)了對目標動態(tài)RCS的模擬，并通過實例說明了這種模擬對于實際工程的意義。目前，采用這種方法模擬顯示目標動態(tài)RCS的程序已經(jīng)應用在某設(shè)備中，并能提供一定意義下的目標飛行特性信息。

［1］劉英芝，曾勇虎，汪連棟，等.動態(tài)雷達目標仿真中目標姿態(tài)角的計算［J］.電光與控制，2007，14(1):38 －41.

［2］張居鳳，馮德軍，王雪松，等.雷達目標動態(tài)RCS仿真研究［J］.系統(tǒng)仿真學報，2005，17(4):834 －837.

［3］黃學德.導彈測控系統(tǒng)［M］.北京:國防工業(yè)出版社，2000.

［4］John H M ，Kurtis D F.數(shù)值方法［M］.周璐，陳渝，錢方，等，譯.北京:電子工業(yè)出版社，2002.

［5］聶在平.目標與環(huán)境電磁散射特性建模——理論、方法與實現(xiàn)(基礎(chǔ)篇)［M］.北京:國防工業(yè)出版社，2009.

［6］胡俊，聶在平，王軍，等.基于多層快速多極子方法的三維目標RCS高效數(shù)值求解技術(shù)［J］.電波科學學報，2004，19(10):92 －94.

［7］ 王超.NURBS曲面 RCS算法研究［D］.西安:西北工業(yè)大學，2004.

［8］吉雅嫻.ESD與TPS之間的數(shù)據(jù)通信［J］.石油化工自動化，2008(3):54－56.