雷達角閃爍的混沌特性研究

張 錦 苗俊剛 鄭 成 謝 衍

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

引 言

在對目標雷達特性的研究中，角閃爍和雷達散射截面（RCS）是同等重要的兩種物理量，是雷達目標的特征信號［1］。角閃爍的概念與擴展目標的概念緊密相連。復(fù)雜目標不同部位的散射強度和相位會造成回波相位的畸變，回波波前在接收天線口面上的傾斜與擺動產(chǎn)生了角閃爍［2］。凡是尺寸與波長能比擬、具有兩個或兩個以上等效散射中心的任何體目標都會產(chǎn)生角閃爍線偏差，這類體目標稱為擴展目標。角閃爍誤差對目標跟蹤和識別有較大影響，強角閃爍信號可能導(dǎo)致目標跟蹤失敗和誤判，因此，如何預(yù)測和抑制角閃爍成為散射研究領(lǐng)域中的重要議題［3］。相關(guān)研究對角閃爍的后處理及抑制方法的討論很多，但只有很少的文獻涉及到角閃爍的混沌特性的詳細判定及與標準模型的結(jié)果對照。因此，討論角閃爍的混沌特性是一個新穎的課題。

首先采用相位梯度法對共線兩散射中心和五散射中心模型的角閃爍進行了仿真計算，并用圖形電磁學(xué)算法（GRECO）計算出了一個典型目標的螺旋運動角閃爍數(shù)據(jù)。隨后設(shè)計了一個全面和細致的非線性和混沌驗證流程算法，并對照業(yè)內(nèi)公認的標準混沌模型Lorenz模型經(jīng)過同樣流程的仿真結(jié)果，對三種角閃爍計算結(jié)果進行了一系列混沌定性、定量判別，確定了角閃爍具有顯著的混沌特性，為抑制角閃爍和目標識別拓展了新思路。

1 角閃爍精確建模

1.1 角閃爍精確建模

相關(guān)領(lǐng)域研究人員對雷達角閃爍的建模提出了幾類方法，包括精確建模、統(tǒng)計建模等。為了驗證角閃爍的非線性及混沌特性，采用相位梯度法對同軸共線多散射中心模型進行角閃爍仿真計算［3－5］，后續(xù)內(nèi)容驗證其混沌特性。

假設(shè)目標由排列在同一條直線上多個散射源組成，如圖1所示，O是坐標原點，θ是水平面上觀測角，散射中心排列在Z軸上。P是雷達位置，ri是第i個目標到雷達的距離。

雷達距離：

單個理想導(dǎo)體球體散射中心的散射場為［6］

式中，a為球體直徑，直徑的大小代表了接收散射場的強度。設(shè)第i個散射中心的散射源幅度為Ai，相位為φi，則多個散射中心的后向散射總場為

圖1 共線多散射中心目標幾何模型

目標對觀察點P處散射場相位為

利用相位梯度法，可得俯仰面和方位面上的角閃爍計算方法如公式（5）所示。限于篇幅，本文只討論方位面上方位面角θ變化時的角閃爍精確計算。雙散射中心角閃爍建模采用兩個與中心線等距且直徑差距較小的散射中心作為建模模型，得到的結(jié)果具有較強的簡單目標普適性。仿真參數(shù)設(shè)置為：f＝10GHz，R＝300m，L1＝－0.3m，L2＝0.3m，a1＝0.1m，a2＝0.15m（符合遠場條件及散射中心體積與波長可比擬的條件），角閃爍量綱為米，仿真步長為0.01°，水平面角度掃描范圍為0°～180°.為計算其時間序列特性，并與實際雷達掃描時的情景相匹配，圖形坐標橫軸采用時間序列的方式顯示為雷達系統(tǒng)中的時域采樣點數(shù)量。角閃爍仿真計算結(jié)果請見后文混沌特性分析中圖5及圖6.可見角閃爍曲線偏向中心正向的距離，即散射中心直徑較大，產(chǎn)生較強角閃爍的方向，符合散射中心偏移理論。

依相同方法可計算五散射中心目標的方位面面角閃爍：L1＝－0.3m，L2＝0.3m，L3＝0.6m，L4＝－0.5m，L5＝1m，a1＝0.1m，a2＝0.15m，a3＝0.1 m，a4＝0.3m，a5＝0.8m，其余參數(shù)不變，仿真結(jié)果如圖3所示。按相關(guān)文獻論述，五散射中心角閃爍已屬復(fù)雜散射體［7］。

1.2 典型目標GRECO角閃爍計算

GRECO算法是對角閃爍精確建模的拓展，被電磁散射理論界公認為是與實際測量最為接近的仿真方法之一［6，8］，是相位梯度法在動目標圖形化實時仿真系統(tǒng)中的應(yīng)用，綜合了物理光學(xué)法、物理繞射法、幾何繞射法等常用電大尺寸目標的散射場計算方法［9－10］，將每個像素視為是一個散射中心，當目標在屏幕上顯示后，所有像素點組成的散射中心模型即代表該復(fù)雜目標，通過獲得每個散射中心的回波場強和位置信息求取復(fù)雜目標的角閃爍。

復(fù)雜目標由n個統(tǒng)計獨立的散射中心組成，假定觀察距離遠大于散射中心與雷達間距和波長，則入射波可以近似看成是平面波，略去時諧因子后各散射中心被雷達接收的總回波信號為

式中：rN表示第N個散射中心到接收雷達處的距離；，δn分別是第n個散射中心的幅度和初始相位。由此得到回波信號的總相位為

式中：

針對目標的每個實時姿態(tài)，由相位梯度法，可計算出水平面和俯仰面的角閃爍線偏差為

GRECO計算的目標模型采用公認的典型目標［8］，目標運動方式［10－11］選為螺旋運動（符合目標運動規(guī)律），目標和坐標設(shè)置如圖2所示。

角閃爍計算參數(shù)如表1所示。

圖2 典型復(fù)雜目標模型［8］

表1 典型目標角閃爍計算參數(shù)

由GRECO算法計算得到的角閃爍波形請見后文混沌分析中圖7.

2 角閃爍時間序列的混沌驗證

圖3 混沌判定流程圖

非線性科學(xué)與混沌理論是近年學(xué)術(shù)界熱點，很多看似用常規(guī)方法無法解決的問題從非線性科學(xué)中都獲得了嶄新的解決思路。在經(jīng)濟學(xué)、地理、氣象數(shù)據(jù)的處理和預(yù)測中，非線性與混沌理論已經(jīng)獲得了廣泛的應(yīng)用。本文參考國外非線性科學(xué)學(xué)者的權(quán)威成果［12－14］，利用混沌理論設(shè)計了一個全面和細致的驗證流程算法（圖3），并對照公認的混沌模型Lorenz吸引子的仿真結(jié)果，驗證了角閃爍的混沌特性。

本流程中混沌特性的判定可依定性［12，14］和定量［15－16］分析算法分別驗證，三個定性判據(jù)包括吸引子、龐加萊截面和主分量分析，兩個定量判據(jù)包括關(guān)聯(lián)維數(shù)和最大Lyapunov指數(shù)。后續(xù)分析將Lorenz吸引子和三種角閃爍序列均通過此流程驗證，得到一系列定性和定量的判定結(jié)果。以下對此驗證流程進行簡要敘述。

1）互信息函數(shù)確定最佳時間延遲。取多個時間延遲進行平均信息量的迭代計算，得到的向量是延遲后以τ為自變量的時間序列重構(gòu)；取其第一個最小值點對應(yīng)的時間點作為延遲的估計值［12－13］。

2）Cao算法確定最小嵌入維數(shù)［17］。①采用上步求得的最佳時間延遲τ和d＋1（為待定嵌入維數(shù)）個嵌入維數(shù)作為輸入?yún)⒘?，采用偽臨近算法定義范數(shù)商參量，并將a參量加權(quán)求和且歸一化得到參量并繪制其圖形；②以a的第一個明顯拐點作為最小嵌入維；③以τ和min－ed重構(gòu)吸引子［15］并觀察其運動特性。通常情況下，混沌的最小嵌入維數(shù)為2或者3，方能成功重構(gòu)吸引子動力學(xué)特性。由此得到的吸引子圖形為第一個混沌定性判據(jù)：吸引子在有限空間內(nèi)不斷伸長和折疊，構(gòu)成回復(fù)性永不相交的非周期運動，具有良好的幾何構(gòu)型，則序列具有混沌特性。

3）相空間重構(gòu)并生成二維龐加萊截面圖［13－14，16］，即在相空間中適當選取一截面（要有利于觀察系統(tǒng)的運動特征和變化，如截面不能和軌跡相切，更不能包含軌跡），稱此截面為龐加萊截面，相空間的連續(xù)軌跡與龐加萊截面的交點稱為截點。設(shè)記錄得到的龐加萊點為：B0，B1，…，Bn，…。在龐加萊截面上使相空間軌跡連續(xù)運動，降為低維的離散點之間的映射，T稱為龐加萊映射［16］。

由此得到的龐加萊截面圖成為第二個混沌定性判據(jù)［17－18］，即：1）當龐加萊截面上有且僅有一個不動點或少數(shù)離散點時，序列是周期的；2）當龐加萊截面是一條封閉曲線時，序列是準周期的；3）當龐加萊截面上是一些成片的具有分形結(jié)構(gòu)或良好幾何構(gòu)型的形體時，序列是混沌的。

4）依據(jù)關(guān)聯(lián)積分G－P算法［18］選擇最佳嵌入維數(shù)重構(gòu)相空間［16］，并求得關(guān)聯(lián)維數(shù)。步驟如下：1）按照前述計算得到最小嵌入維，并以此重構(gòu)相空間，然后計算關(guān)聯(lián)積分；2）不斷增大嵌入維數(shù)并重新重構(gòu)相空間，觀察關(guān)聯(lián)積分曲線圖形，選擇除去斜率為0和無窮的曲線中最接近直線的曲線，求得該曲線的斜率即為關(guān)聯(lián)維數(shù)D；3）考察關(guān)聯(lián)積分曲線的收斂性，若曲線在某嵌入維數(shù)以上呈發(fā)散狀態(tài)并密集聚攏，則該嵌入維數(shù)即為最佳嵌入維數(shù)opt－ed.由此得到關(guān)聯(lián)維數(shù)的定量混沌判據(jù)：如果關(guān)聯(lián)維數(shù)是整數(shù)，則序列是準周期或隨機的；如果關(guān)聯(lián)維數(shù)不是整數(shù)，則序列是混沌的。

5）以最佳嵌入維數(shù)opt－ed重構(gòu)相空間，進行主分量譜圖分析［16，19］。由時間延遲和最佳嵌入維數(shù)重構(gòu)相空間，求出軌線矩陣；計算協(xié)方差矩陣，并計算出其特征值和特征向量；將特征值按由大到小排列，則特征值和對應(yīng)的特征向量稱為該時間序列的主分量。求出特征值的和γ，以γ為縱坐標，嵌入維數(shù)1至opt－ed為橫坐標繪制出的圖形即為主分量譜圖。由此得到主分量譜圖判據(jù)：隨機噪聲序列的主分量譜圖是一條與橫軸接近平行的直線，而混沌序列則表現(xiàn)為一條沿負斜率下降的曲線，且在低維迅速下降。

6）Wolf法計算最大李氏（Lyapunov）指數(shù)［20］。該指數(shù)的計算方法為最小臨近點算法。將重構(gòu)相空間后的序列取一初始點，找到該點的最鄰近點，追蹤兩點的時間演化，直到二者間距大于某門限值，保留該點并再尋找下一點，使距離小于該門限值，并且前兩點與原始點連線的夾角盡可能小，對序列中每點重復(fù)上述過程直到序列終點，再用對數(shù)距離商的加權(quán)平均求得最大李氏指數(shù)。其定量判據(jù)為：若該指數(shù)大于零，則序列是混沌的。

3 仿真結(jié)果與分析

3.1 Lorenz模型的建立

Lorenz是典型的混沌系統(tǒng)，其迭代微分方程組為［15－16］：

為清晰地分析角閃爍的混沌特性，以標準Lorenz模型作為對照模型，給出該模型和角閃爍模型在定性和定量特性上的相似點。參數(shù)設(shè)置為：序列長度與角閃爍同為18 000點，步長為0.025，σ，r，b的取值分別為0.1，－0.1，－0.02，考慮x分量序列，如圖4左上角所示。在混沌定性驗證分析中，所有軌跡均為歸一化長度，故運動軌跡無量綱。

3.2 混沌定性分析結(jié)果比對與分析

如圖4～7所示，混沌定性分析產(chǎn)生三個定性判據(jù)：吸引子圖形，主分量分析和龐加萊截面。

從圖4～7右上角吸引子圖形可得，Lorenz模型與三個角閃爍序列吸引子結(jié)果均在有限空間內(nèi)不斷伸長和折疊，構(gòu)成回復(fù)性永不相交的非周期運動，且具有良好的幾何構(gòu)型，從運動形態(tài)上均滿足混沌序列的吸引子特性。

圖7 典型目標角閃爍混沌定性分析

從圖4～7左下角各圖可得，四者的二維龐加萊截面圖上均為成片具有分形結(jié)構(gòu)的密集點，并且密集區(qū)域與非密集點區(qū)域有明顯的界限，符合混沌序列的龐加萊截面特征。

從圖4～7右下角各圖可得，四者的主分量均具有沿負斜率下降的特性，并且在低嵌入維時主分量迅速下降（圖7的主分量分析明顯表現(xiàn)了該特征），故符合混沌的主分量譜圖特性。

故由定性分析三個判據(jù)及判定標準，三組角閃爍序列具有明確的混沌特性。

3.3 混沌定量分析結(jié)果比對與分析

根據(jù)圖3所示計算流程得到兩種定量判據(jù)：關(guān)聯(lián)維數(shù)和最大李氏指數(shù)的結(jié)果如表2所示，四者關(guān)聯(lián)維數(shù)均為分數(shù)，最大李氏指數(shù)均大于零，且三組角閃爍序列的最大李氏指數(shù)與Lorenz系統(tǒng)的最大李氏指數(shù)處于同一數(shù)量級。故由定量分析的兩個判據(jù)及判定標準，三種角閃爍序列均表現(xiàn)出明顯的混沌特性。

表2 混沌定量分析結(jié)果

4 結(jié) 論

為獲得角閃爍原始數(shù)據(jù)，首先進行了同一直線上多散射中心目標的角閃爍仿真，并對典型目標采用GRECO算法進行典型目標動態(tài)角閃爍仿真，隨后設(shè)計了詳細的混沌驗證流程，以Lorenz系統(tǒng)作為對照，從定性和定量兩個方面驗證了三個角閃爍模型的混沌特性。混沌相關(guān)研究表明：混沌具有短期可預(yù)測性，且可以依據(jù)相關(guān)算法對角閃爍進行抑制，從而達到減少丟失目標概率和提高目標識別能力的目的。因此，繼續(xù)研究角閃爍的混沌特性及混沌抑制措施具有很高的理論和應(yīng)用價值。

［1］HOWARD D D.Radar target glint in tracking and guidance system based on echo signal phase distortion［C］／／Proc of NEC，1959，5：840－849.

［2］王 濤，王雪松，肖順平.一種極化測量雷達的角閃爍抑制方法［J］.電波科學(xué)學(xué)報，2004，19（6）：702－707.

WANG Tao，WANG Xuesong，XIAO Shunping.Suppressing the angular glint of polarimetric radar ［J］.Journal of Radio Science，2004，19（6）：702－707.（in Chinese）

［3］黃培康.雷達目標特征信號［M］.北京：宇航出版社，1993.

［4］SHI Xiaoli，LIU Yongcai.Angular glint of aircraft formation and its applications［C］／／IEEE ICMA.Harbin，5－8August，2007：1334－1339.

［5］王國玉，汪連棟.雷達目標角閃爍的建模與仿真.雷達火控技術(shù)，2000，29（2）：1－5.

WANG Guoyu，WANG Liandong.Modeling and simulations for radar angular glint［J］.Radar and Weapon Control，2000，29（2）：1－5.（in Chinese）

［6］方 寧，寧煥生，王寶發(fā).復(fù)雜目標可視化角閃爍偏差計算［J］.北京航空航天大學(xué)學(xué)報，2006，32（2）：186－189.

FANG Ning，NING Huansheng，WANG Baofa.Visual computation of glint error for complex targets［J］.Journal of Beijing University of Aeronautics and Astronautics，2006，32（2）：186－189.（in Chinese）

［7］DELANO R H.A theory of target glint or angular scintillation in radar tracking［C］／／Proc IRE，1952（12）：1778－1734.

［8］RIUS J M，F(xiàn)ERRANDO M，JOFRE L.GRECO：graphical electromagnetics computing for RCS prediction in real time［J］.IEEE Antennas and Propagation Magazine，1993， 35 （2）： 7－17.doi：10.1109／74.207645.

［9］SUI Miao，XU Xiaojian.NFIPO based technique for near－field angular glint calculation of radar targets［J］.IEEE Transactions on Aerospace and Electronic Systems，2008，44（1）：372－380.

［10］曾永虎，王國玉，陳永光，等.動態(tài)雷達目標RCS的統(tǒng)計分析［J］.電波科學(xué)學(xué)報，2007，22（4）：610－613.

ZENG Yonghu，WANG Guoyu，CHEN Yongguang，et al.Statistical analysis for RCS of dynamic radar target［J］.Chinese Journal of Radio Science，2007，22（4）：610－613.（in Chinese）

［11］陳博韜，雷振亞，謝擁軍，等.基于改進物理光學(xué)法的電大目標雙站RCS的預(yù)估［J］.電波科學(xué)學(xué)報，2010，25（5）：960－965.

CHEN Botao，LEI Zhenya，XIE Yongjun，et al.Bistatic RCS prediction with improved PO method for electrically large targets［J］.Journal of Radio Science，2010，25（5）：960－965.（in Chinese）

［12］WOLF A，BESSOIR T.Diagnosing chaos in the space circle［J］.Physica D：Nonlinear Phenomena，1991，50（2）：239－258.

［13］MEISS J D.Asked questions about nonlinear science［R／OL］.［2011－10－13］http：／／amath.colorado.edu／pub／dynamics／papers／sci.nonlinearFAQ.pdf.

［14］CAMPBELL D K.Nonlinear Science from Paradigms and Practicalities［J］.Los Alamos Science，1987，15：218－262.

［15］李曉輝，李 湘，郭桂榮.擴展目標角閃爍非線性檢驗及其抑制技術(shù)［J］.系統(tǒng)工程與電子技術(shù)，2004，26（5）：592－595.

LI Xiaohui，LI Xiang，GUO Guirong.Nonlinearity detection and suppression technique for angle glint of extended targets［J］.System Engineering and Electronics，2004，26（5）：592－595.（in Chinese）

［16］MERKWIRTH C，PARLITZ U.TSTOOL User Manual Version 1.11［M／OL］.［2011－10－13］http：／／www.dpi.physik.uni－goettingen.de／tstool／manual.pdf.

［17］呂金虎，陸君安.混沌時間序列分析及其應(yīng)用［M］.武漢：武漢大學(xué)出版社，2002.

［18］CAO Liangyue.Practical method for determining the minimum embedding dimension of a scalar time series［J］.Physica D：Nonlinear Phenomena，1997，110（1／2）：43－50.

［19］GRASSBERGER P，PROCACCIA I.Characterization of stranger attractors［J］.Phys Rec Lett，1983，50（5）：346－349.

［20］WOLF A，SWIFT J B，SWINNEY L，et al.Determining Lyapunov exponents from a time series［J］.Physica D：Nonlinear Phenomena，1985，16（3）：285－317.