高 飛，劉建濤，王澤濤，王力寶，蔣國喜

(1.空軍預(yù)警學(xué)院三系， 武漢430019;2.駐上海航天局中心軍事代表室， 上海201109)

(3.國防科技大學(xué)電子科學(xué)與工程學(xué)院， 長沙410073;4.解放軍94535部隊， 江蘇徐州221005)

0 引言

機載預(yù)警雷達具有靈活機動、登高望遠的空中情報支援特點和較高的雷達測量精度，在進行遠程預(yù)警的同時還可以完成空中力量的指揮和引導(dǎo)作戰(zhàn)，其在現(xiàn)代信息化戰(zhàn)爭中的地位越來越受到重視，各國也在競相研究高性能的機載預(yù)警雷達。然而，以機載預(yù)警雷達為核心的相關(guān)技術(shù)研究是一個復(fù)雜而系統(tǒng)的工程，通常耗資巨大且技術(shù)難度高。20世紀90年代，美國根據(jù)機載預(yù)警雷達技術(shù)發(fā)展的需要，僅在機載雷達數(shù)據(jù)采集上就耗資上億美元，相繼實施了Mountaintop［1-3］和 MCARM［4］計劃。我國在此方面研究起步相對較晚，技術(shù)相對滯后，早期研究也缺乏相關(guān)數(shù)據(jù)的支撐，然而，隨著機載預(yù)警雷達數(shù)據(jù)和系統(tǒng)模型及相關(guān)仿真技術(shù)的發(fā)展，這一困境得到了明顯改善。高逼真一體化的機載預(yù)警雷達系統(tǒng)仿真不但可以靈活地為當前的研究提供各種不同外界環(huán)境(如:雜波環(huán)境、干擾環(huán)境和目標環(huán)境)下的數(shù)據(jù)需求，還可以滿足預(yù)先研究先進信號處理技術(shù)和數(shù)據(jù)處理技術(shù)的驗證需求。

文獻［5-7］從仿真實現(xiàn)技術(shù)的角度研究了雷達系統(tǒng)仿真引擎機制下雷達系統(tǒng)各執(zhí)程體間的通信問題，提高了雷達系統(tǒng)仿真過程的控制與監(jiān)視能力;文獻［8］在考慮雷達和目標的情況下，利用聯(lián)合建模與仿真環(huán)境，在信號層面對對空警戒雷達工作過程進行了仿真，但仿真系統(tǒng)中未考慮干擾，且目標模型中考慮的因素較少，還不夠逼真;文獻［9］研究了地基相控陣雷達相干視頻信號建模與仿真，給出了仿真系統(tǒng)的框架流程及關(guān)鍵信號模型;文獻［10］針對地基相控陣雷達系統(tǒng)，進行了功能級仿真和相干視頻信號級仿真，但信號級仿真中各模塊內(nèi)部的處理流程不夠明確;文獻［11］分析研究了相控陣雷達系統(tǒng)仿真的總體流程，構(gòu)建了仿真系統(tǒng)的模塊化結(jié)構(gòu);文獻［12］基于“組態(tài)可變、一機多用”的設(shè)計理念，提出了一種單脈沖雷達半實物仿真方法。

近年來隨著微波集成電路和高速/超高速數(shù)字電路技術(shù)的快速發(fā)展，數(shù)字雷達接收機的應(yīng)用越來越普遍，無論是射頻數(shù)字雷達接收機還是中頻數(shù)字雷達接收機，都可以形成零中頻的數(shù)字信號，為后續(xù)采取靈活的信號處理技術(shù)提供了可能。因此，雷達系統(tǒng)中從零中頻數(shù)據(jù)產(chǎn)生至后端的信號處理和數(shù)據(jù)處理的仿真研究受到普遍重視。本文針對相控陣機載預(yù)警雷達系統(tǒng)的特點，從雷達系統(tǒng)外部的零中頻數(shù)據(jù)源仿真(包括雜波數(shù)據(jù)、干擾數(shù)據(jù)和目標數(shù)據(jù))到系統(tǒng)內(nèi)部的信號處理(包括多個子模塊)和數(shù)據(jù)處理功能仿真進行了系統(tǒng)地研究，為了盡可能完整逼真地反映相控陣機載預(yù)警雷達系統(tǒng)探測目標的實際過程，仿真過程考慮了更為復(fù)雜的對空探測模式。即以一定的載機和目標機航線及電磁環(huán)境想定為背景，對各模塊節(jié)點信號進行輸出觀察，并將終端輸出的目標航跡與想定航線進行比較。

1 系統(tǒng)仿真流程模塊

相控陣機載預(yù)警雷達系統(tǒng)中從零中頻數(shù)據(jù)產(chǎn)生至后端處理的仿真，按照仿真流程模塊可分為三部分:

1)數(shù)據(jù)源產(chǎn)生模塊。負責(zé)產(chǎn)生陣元級/通道級的多維數(shù)據(jù)。目前，機載預(yù)警雷達對空探測模式通常采用多個發(fā)射重頻的脈沖多普勒(Pulse Doppler，PD)體制，因此，在每個波位產(chǎn)生的雜波數(shù)據(jù)為多個重頻對應(yīng)的多個相干脈沖間隔(Coherent Processing Interval，CPI)內(nèi)的數(shù)據(jù)立方體;在有目標(單個點目標)的波位目標回波數(shù)據(jù)為多個重頻對應(yīng)的各自距離單元內(nèi)的空時二維矢量;在有干擾(人為施放的有源干擾)的波位干擾數(shù)據(jù)與目標數(shù)據(jù)形式相似，只是在距離上呈或離散或連續(xù)或離散連續(xù)混合的分布。

2)信號處理模塊。負責(zé)對數(shù)據(jù)源產(chǎn)生模塊產(chǎn)生的數(shù)據(jù)進行一系列處理，從而提取目標的相關(guān)參數(shù)?？紤]到相控陣機載預(yù)警雷達的體制特點，其信號處理中的雜波/干擾抑制除了可以采用常規(guī)PD處理，還可以靈活地采用先進的空時自適應(yīng)處理(Space-Time A-daptive Processing，STAP)技術(shù)，與常規(guī)地面對空警戒雷達目標檢測不同的是，機載預(yù)警雷達需要利用多個重頻的恒虛警率(Constant False Alarm Rate，CFAR)檢測結(jié)果進行聯(lián)合解模糊(解距離模糊和速度模糊)與目標檢測。

3)數(shù)據(jù)處理模塊。負責(zé)對信號處理模塊得到的目標參數(shù)(如距離、方位等)進行點跡處理和航跡處理，其處理過程與常規(guī)地面雷達的數(shù)據(jù)處理基本相同。

按照雷達系統(tǒng)的外部信號環(huán)境和內(nèi)部處理功能區(qū)分，我們將仿真系統(tǒng)中的數(shù)據(jù)源產(chǎn)生模塊劃為外部環(huán)境模塊，將信號處理和數(shù)據(jù)處理模塊劃為內(nèi)部處理模塊。圖1給出了機載預(yù)警雷達系統(tǒng)仿真流程模塊圖，其中仿真環(huán)境主要是對載機航線、目標機航線及電磁環(huán)境進行規(guī)劃設(shè)定，從而為整個仿真系統(tǒng)提供仿真參數(shù)，外部環(huán)境模塊和內(nèi)部處理模塊的內(nèi)部子功能將在下面具體分析。

圖1 機載預(yù)警雷達系統(tǒng)仿真流程模塊圖

1.1 外部環(huán)境模塊

外部環(huán)境模塊主要完成對雜波、干擾和目標信號的仿真，由于機載預(yù)警雷達采用的多個發(fā)射重頻之間在仿真上是獨立而對稱的，所以，此處只對一個重頻CPI內(nèi)的零中頻數(shù)據(jù)(后續(xù)的數(shù)據(jù)均指零中頻數(shù)據(jù))進行仿真，其他重頻的情況與此相同。

1.1.1 雜波仿真

常規(guī)機載預(yù)警雷達雜波仿真中，通常獨立地直接產(chǎn)生各個距離單元的雜波數(shù)據(jù)，沒有考慮各距離門前面若干個距離門回波的時延疊加效果和波形調(diào)制影響，因此，不能用于考察后續(xù)不同脈沖壓縮(Pulse Compression，PC)處理的效果。從PC前回波數(shù)據(jù)的實際情況和PC功能驗證的需求考慮，我們的仿真中反映了各距離門前面若干個距離門回波的時延疊加和波形調(diào)制影響。仿真的思路是，先獨立產(chǎn)生各個距離單元的雜波數(shù)據(jù)，然后，考慮波形調(diào)制情況下形成量化時延疊加的效果。其中，未考慮波形調(diào)制和時延疊加情況下獨立產(chǎn)生各個距離單元雜波數(shù)據(jù)的方法已經(jīng)在多個文獻中進行了研究［13-15］，本文不再研究。下面主要研究利用獨立產(chǎn)生的各距離單元雜波數(shù)據(jù)形成量化時延疊加效果的仿真方法。

假設(shè)雷達空域接收通道為N，發(fā)射K個相干脈沖，在每個脈沖發(fā)射結(jié)束的間歇期進行接收采樣，采樣起始時刻為t0，采樣周期為Ts，采樣得到的距離門個數(shù)為M，則不考慮波形調(diào)制和回波時延疊加時，第m個距離門的空時雜波數(shù)據(jù)矢量可表示為cm(m=1，2，…，M);若雷達采用的調(diào)制波形是脈寬為τ的線性調(diào)頻矩形脈沖信號，則該信號可表示為

式中:rect(·)為矩形函數(shù)，當0≤t＜τ時 rect(t，τ)=1，否則rect(t，τ)=0;μ為調(diào)頻斜率。以Ts對式(1)進行離散采樣，可得到第m個距離門的回波調(diào)制信號為

式中:n為離散時間變量，其對應(yīng)t0+(n-1)Ts時刻。實際上還可以證明，即使存在雜波距離模糊時，第m個距離門的回波調(diào)制信號仍為式(2)，此處不再推導(dǎo)。由于發(fā)射脈沖寬度為τ，某個距離門的回波采樣數(shù)據(jù)實際包含了該距離門前面τc/2(c為光速)范圍以內(nèi)所有距離門回波時延的疊加。令 Mτ=［τ/Ts］，「·」表示向上取整，考慮到起始采樣時刻(n=1)附近若干距離門回波時延疊加會受到cm邊界特性(cm取值在1≤m≤M范圍)的影響，那么在考慮波形調(diào)制和回波時延疊加情況下，第m個距離門的空時雜波數(shù)據(jù)矢量可表示為

對式(3)可以結(jié)合圖2給出的回波時延疊加的示意圖進行理解。

圖2 回波時延疊加示意圖

從式(3)可以看出，與cm相比，xm增加了波形幅度和相位調(diào)制以及回波時延疊加的影響，反映了PC前零中頻雜波數(shù)據(jù)的真實情況。

1.1.2 干擾仿真

目前，雷達干擾類型較為繁多，從干擾的意圖看，可分為有意干擾和無意干擾。對于預(yù)警雷達而言，無意干擾主要指氣象、鳥群、大型強散射地物及鄰近雷達、通信、電臺等電子設(shè)備無意中對雷達造成的電磁干擾。從廣義上講，雜波也可以看作一種無意干擾;有意干擾主要指為了干擾雷達正常工作而人為施放的各種有源和無源干擾。有意干擾從干擾功能表現(xiàn)上看，可分為壓制性干擾和欺騙性干擾，實際中施放的可能是二者的混合干擾或者同時具有壓制功能和欺騙功能的干擾，如靈巧干擾。各種類型的干擾模型都有較多研究，本文主要給出相控陣機載預(yù)警雷達中常見的有源壓制干擾和欺騙干擾的仿真方法。

1)有源壓制干擾模型

(1)射頻干擾模型。在仿真中均假設(shè)干擾載頻與雷達發(fā)射信號載頻匹配且不計載頻項(以下均如此假設(shè))，則射頻干擾信號模型為

式中:An(t)是均值為0，方差為的高斯白噪聲;φ(t)為［0，2π］上均勻分布且與An(t)獨立的隨機變量。

(2)噪聲調(diào)幅干擾模型。噪聲調(diào)幅干擾信號模型為

式中:An為載波電壓。

(3)噪聲調(diào)頻干擾模型。噪聲調(diào)頻干擾信號模型為

式中:KFM為調(diào)頻比例系數(shù)，表示單位調(diào)制信號強度所引起的頻率變化。

(4)噪聲調(diào)相干擾模型。噪聲調(diào)相干擾信號模型為

式中:KPM為調(diào)相比例系數(shù)，表示單位調(diào)制信號強度所引起的相位變化。

文中制作的高通濾波器由上圖所示的級聯(lián)電容與接地電感組成。為了節(jié)省材料，文中電容輸入與輸出端和饋電接口由小段傳輸線連接，亦可直接由電容的一端直接連接饋電端口。C2與C3亦可由一塊金屬面版來替代兩塊由一段傳輸線級聯(lián)的電容。

可以看出，有源壓制干擾J(t)為一隨機變量，設(shè)其產(chǎn)生的K×1維隨機時域干擾矢量為Jt=［J1，J2，…，JK］T對應(yīng)的空域?qū)蚴噶繛镴s(由干擾相對天線的來向決定)，則對應(yīng)該有源壓制干擾，第m個距離門的空時干擾數(shù)據(jù)矢量可表示為

式中:?為Kronecker積。

2)欺騙干擾模型

按照目標參數(shù)信息，欺騙干擾可分為距離欺騙干擾、速度欺騙干擾、角度欺騙干擾和多參數(shù)欺騙干擾，本文主要給出常見的距離欺騙干擾和速度欺騙干擾模型。距離欺騙干擾和速度欺騙干擾通常是通過距離波門拖引和速度波門拖引技術(shù)實現(xiàn)的，主要是通過對目標的正常回波信號實施一定的延時和多普勒頻移來實現(xiàn)。距離/速度欺騙干擾信號模型可表示為

式中:AJ為干擾信號幅度;c為光速;fd和R0為正常目標的多普勒頻率和距離;Δt和Δfd分別為實施距離欺騙和速度欺騙時的拖引量;f'd和Δt'為假目標的多普勒頻率和距離對應(yīng)的時延。設(shè)假目標在第m'個距離門，即Δt'=t0+(m'-1)Ts，可得到該假目標的回波調(diào)制信號為 Sτ(n，m')=rect［(n-m')Ts，τ］·ejπμ［(n-m')Ts］2，由于t=t0+(n-1)Ts，由式(9)可知快時間上的多普勒信號相位項為 ud(n)=ej2πfd［t0+(n-1)Ts］。設(shè)由 J(t)產(chǎn)生的K×1維時域干擾矢量為Jt=［J(0)，J(Tr)，…，J［(K-2)Tr］，J［(K-1)Tr］］T，其中 Tr為發(fā)射脈沖重復(fù)周期。那么對應(yīng)該欺騙干擾，第m個距離門的空時干擾數(shù)據(jù)矢量可表示為

1.1.3 目標仿真

真實目標信號模型與欺騙干擾信號模型較類似，可表示為

1.2 內(nèi)部處理模塊

1.2.1 信號處理仿真

機載預(yù)警雷達信號處理流程示意圖，如圖3所示。其中，各模塊功能的實現(xiàn)算法已經(jīng)比較成熟［11，16］，此處不再贅述。值得說明的是，作為信號處理核心功能之一的雜波/干擾抑制，可以采用常規(guī)非自適應(yīng)處理，如PD處理、動目標顯示或者動目標檢測等，對于空域多通道體制的機載雷達還可以采用STAP技術(shù);而對于采取多個發(fā)射重頻的機載預(yù)警雷達，此處的CFAR處理除了要在每個重頻內(nèi)進行，還包括利用多個重頻CFAR結(jié)果進行聯(lián)合解模糊與目標檢測，其意味著單個重頻內(nèi)實現(xiàn)恒虛警率檢測的同時也實現(xiàn)了多個重頻聯(lián)合的恒虛警率檢測，保證了系統(tǒng)總的恒虛警率要求。

圖3 機載預(yù)警雷達信號處理流程示意圖

1.2.2 數(shù)據(jù)處理仿真

圖4給出了機載預(yù)警雷達數(shù)據(jù)處理流程示意圖，其原理與地面對空警戒雷達基本相同，只是由于載機運動及姿態(tài)變化，當需要知道不同坐標系下目標狀態(tài)時，通常需要對雷達的直接量測結(jié)果在點跡處理之后進行運動補償和坐標轉(zhuǎn)換。雷達數(shù)據(jù)處理的主要任務(wù)可分為點跡處理和航跡處理，其中點跡處理是對錄取的點跡數(shù)據(jù)進行野值剔除、濾波和凝聚處理，而航跡處理是將同一目標的點跡連成航跡的處理過程，一般包括航跡起始、點/航跡相關(guān)(數(shù)據(jù)關(guān)聯(lián))、航跡濾波與預(yù)測及航跡起始/終止。雷達數(shù)據(jù)處理中各個環(huán)節(jié)的具體原理和仿真方法較為繁多，文獻［16］和文獻［17］分別進行了詳細的介紹，可進行參考。

2 系統(tǒng)各模塊聯(lián)合仿真實驗

仿真實驗中，想定設(shè)計為:某P波段機載預(yù)警雷達，以7 600 m飛行高度水平向正北(定為極坐標零度方向)飛行，速度為 132 m/s，雷達工作在高重頻(HPRF)PD模式，采用2/4檢測準則，天線為方位機械掃描的4×18陣列，并按列子陣合成18個接收通道，掃描周期為12 s;兩個飛機目標RCS分別為10 m2和5 m2且起伏模型為Sweilling I型，目標1以300 m/s速度向南偏西20°方向水平飛行，目標2以300 m/s速度向南偏東40°方向水平飛行，兩目標起始位置(斜距，方位，高度)分別為(100 km，30°，4 000 m)和(200 km，-10°，5 000 m);干擾為兩個射頻噪聲干擾，干噪比分別為20 dB和30 dB，來向(方位，俯仰)分別為(66°，-2°)和(108°，-3°);雷達探測目標累計38幀，時長7.6 min。

圖4 機載預(yù)警雷達數(shù)據(jù)處理流程示意圖

2.1 回波仿真實驗

圖5和圖6分別給出了某重頻波束方位指向為30.88°(對應(yīng)天線陣面斜側(cè)角為-59.12°)時，21.33 km處距離門內(nèi)無干擾和存在干擾(兩干擾相對陣面軸向的空間錐角分別為125.1°和167.78°)情況下，參數(shù)計算與回波仿真的空時二維雜波譜(錐角-多普勒譜)分布。

圖5 無干擾情況下參數(shù)計算與回波仿真的空時二維雜波譜分布

從圖5和圖6中參數(shù)計算與回波仿真結(jié)果比較可以看出，回波仿真得到的雜波/干擾空時二維譜分布與通過參數(shù)理論計算得到的結(jié)果非常吻合，只是回波仿真中考慮了雷達波束寬度、天線轉(zhuǎn)動及雜波起伏導(dǎo)致的譜展寬與距離模糊等因素，因此，仿真的結(jié)果存在一定的譜展寬和高重頻導(dǎo)致的距離模糊。仿真結(jié)果反映了單個距離門內(nèi)雜波/干擾仿真的逼真性。

2.2 信號處理仿真實驗

圖7a)～圖7c)分別給出與圖5和圖6相同重頻和波束指向情況下，無干擾時的距離-多普勒譜分布的參數(shù)計算結(jié)果及回波仿真脈壓前后結(jié)果。從圖中可以看出，回波仿真的距離-多普勒譜分布與參數(shù)計算結(jié)果比較吻合，由于仿真中參數(shù)及譜展寬影響，近程雜波隨距離的非均勻彎曲不太明顯，從目標和高度線雜波可以看出，脈壓前后回波信號在距離-多普勒譜上的變化與理論實際相符。

圖6 存在干擾情況下參數(shù)計算與回波仿真的空時二維雜波譜分布

圖7 無干擾時參數(shù)計算與回波仿真的距離-多普勒譜分布

圖8給出了與圖7相同重頻和波束指向情況下，主雜波跟蹤后常規(guī)PD處理方法和自適應(yīng)3DT方法的雜波抑制效果(距離-多普勒譜)比較。從圖中可以看出，相對常規(guī)PD處理，3DT處理的雜波抑制性能在低速區(qū)(主雜波附近)具有顯著優(yōu)勢，可以大大提高低速區(qū)信雜噪比。圖9給出了圖8兩種雜波抑制算法處理后采用十字CFAR處理的結(jié)果。雖然兩種雜波抑制算法后的CFAR結(jié)果基本相同，但由于雜波抑制性能的差異，3DT方法低速區(qū)雜波遮擋要比常規(guī)PD處理小很多，因而CFAR處理時低速區(qū)舍棄的檢測范圍要比常規(guī)PD處理小很多(可以從過門限的高度線雜波在多普勒上的寬度看得出)，這一點非常利于低速小目標的檢測。

圖8 兩種雜波抑制算法處理效果比較

圖9 兩種雜波抑制算法處理后的CFAR結(jié)果

2.3 數(shù)據(jù)處理仿真實驗

實際中考慮到載機橫滾和俯仰等姿態(tài)變化及坐標系的選取，需要在數(shù)據(jù)處理的點跡處理之后航跡處理之前進行運動補償和坐標變換。本節(jié)中載機如上述參數(shù)飛行，并采用載機坐標系(以載機為原點)，航跡處理中采用聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)濾波(JPDAF)算法以及3/4航跡起始準則。圖10給出了38幀數(shù)據(jù)處理輸出與想定比較的類似PPI顯示結(jié)果。從圖中比較可以看出，數(shù)據(jù)處理輸出的兩目標航跡與想定設(shè)計較為吻合，只是由于雷達量測和跟蹤誤差，輸出航跡存在些許抖動。

圖10 數(shù)據(jù)處理輸出與想定比較

數(shù)據(jù)處理過程中目標徑向速度的量測值與真實值比較，如圖11所示?？梢钥闯?，在目標徑向速度較高的大多時候，其量測值與真實值基本相符(包含一定的誤差)，而當目標徑向速度較低(接近主雜波)時，由于雜波遮擋，直接的點跡量測丟失(存在外推的填充值)，根據(jù)量測點跡計算的目標量測速度亦丟失(圖中未給出利用外推填充值計算的結(jié)果)。

圖11 數(shù)據(jù)處理過程中目標徑向速度的量測值與真實值比較

3 結(jié)束語

基于相控陣機載預(yù)警雷達技術(shù)發(fā)展中對數(shù)據(jù)來源和技術(shù)驗證的需求，本文從高逼真一體化仿真的角度，將雷達系統(tǒng)中從零中頻數(shù)據(jù)產(chǎn)生至后端的信號處理和數(shù)據(jù)處理整套流程的仿真，分為外部環(huán)境模塊和內(nèi)部處理模塊進行研究，對體現(xiàn)高逼真特性的主要模塊-外部環(huán)境模塊內(nèi)部各功能(產(chǎn)生各種數(shù)據(jù)源)進行了詳細研究，并在內(nèi)部處理模塊中利用現(xiàn)有的信號處理算法和數(shù)據(jù)處理算法對外部環(huán)境模塊產(chǎn)生的數(shù)據(jù)進行處理，實現(xiàn)一體化仿真。為了盡可能驗證仿真的高逼真和一體化特性，文中設(shè)計了仿真想定并參照實際雷達系統(tǒng)進行過程仿真，比較了模塊中各節(jié)點輸出信號與參數(shù)的理論計算結(jié)果，以及最終的數(shù)據(jù)處理輸出與想定設(shè)計。結(jié)果表明:該仿真過程保證了相控陣機載預(yù)警雷達回波數(shù)據(jù)較高的逼真性，反映了其數(shù)據(jù)采樣至處理的一體化過程。

［1］ Berin M O，Haimovich A M.Signal cancellation effects in adaptive radar mountaintop data-set［C］//IEEE International Conference on Acoustics，Speech，and Signal Processing.Newark，USA:IEEE Press，1996:2614-2617.

［2］ Seliktar Y，Williams D B，McClellan J H.Evaluation of partially adaptive STAP algorithms on the mountain top data set［C］//IEEE International Conference on Acoustics，Speech，and Signal Processing.Atbnta，USA:IEEE Press，1996:1169-1172.

［3］ Titi G W，Marshall D F.The A/N mountaintop program:adaptive signal processing for airborne early warning radar［C］//IEEE InternationalConference on Acoustics，Speech，and Signal Processing.Lexington，USA:IEEE Press，1996:1165-1168.

［4］ Little M O，Berry W P.Real-time multichannel airborne radar measurements［C］//IEEE National Radar Conference.［S.l.］:IEEE Press，1997:138-142.

［5］ 王 磊，陳明燕，張 偉，等.基于消息驅(qū)動的雷達仿真引擎控制算法［J］.系統(tǒng)工程與電子技術(shù)，2012，34(9):1807-1812.Wang Lei，Chen Mingyan，Zhang Wei，et al.Control algorithm of radar simulation engine based on message driven［J］.Systems Engineering and Electronics，2012，34(9):1807-1812.

［6］ 王 磊，盧顯良，陳明燕，等.雷達系統(tǒng)仿真中基于消息驅(qū)動的通信服務(wù)模型之研究［J］.計算機應(yīng)用研究，2011，28(10):3767-3774.Wang Lei，Lu Xianliang，Chen Mingyan，et al.Study on communication model based on message driven in radar system simulation［J］.Application Research of Computers，2011，28(10):3767-3774.

［7］ 王 磊，盧顯良，陳明燕，等.一種基于分布式結(jié)構(gòu)的雷達系統(tǒng)仿真引擎機制研究［J］.計算機應(yīng)用研究，2012，29(1):220-223.Wang Lei，Lu Xianliang，Chen Mingyan，et al.Study on engine based mechanism of radar system simulation on distributed architecture［J］.Application Research of Computers，2012，29(1):220-223.

［8］ 欒 磊，王全民，李 革.基于KD-JMASE的信號級雷達系統(tǒng)仿真［J］.計算機仿真，2011，28(9):18-21.Luan Lei，Wang Quanmin，Li Ge.Radar system simulation in signal level based on KD-JMASE［J］.Computer Simulation，2011，28(9):18-21.

［9］ 王 象，李 盾，畢 莉，等.相控陣雷達相干視頻信號建模與仿真方法研究［J］.系統(tǒng)仿真學(xué)報，2010，22(3):741-747.Wang Xiang，Li Dun，Bi Li，et al.Coherent video modeling and simulation method of phased array radar［J］.Journal of System Simulation，2010，22(3):741-747.

［10］ 王 路，邢清華.改進型功能級相控陣雷達系統(tǒng)仿真［J］.現(xiàn)代防御技術(shù)，2013，41(1):164-169.Wang Lu，Xing Qinghua.Improved functional system simulation of phased array radar［J］.Modern Defence Technology，2013，41(1):164-169.

［11］ 童維健，蔡桂友.相控陣雷達仿真技術(shù)研究［J］.現(xiàn)代雷達，2008，30(4):45-49.Tong Weijian，Cai Guiyou.Research on simulation technique of phased array radar［J］.Modern Radar，2008，30(4):45-49.

［12］ 張殿宗，趙嚴冰.基于組態(tài)可變的單脈沖雷達半實物仿真方法［J］.現(xiàn)代雷達，2013，35(10):76-80.Zhang Dianzong，Zhao Yanbing.A hardware-in-loop simulation method for mono-pulse radar based on changing configuration［J］.Modern Radar，2013，35(10):76-80.

［13］ 王永良，彭應(yīng)寧.空時自適應(yīng)信號處理［M］.北京:清華大學(xué)出版社，2000.Wang Yongliang，Peng Yingning.Space-time adaptive processing［M］.Beijing:Tsinghua University Press，2000.

［14］ 杜鵬飛.機載預(yù)警雷達恒虛警率檢測方法研究［D］.長沙:國防科技大學(xué)，2003.Du Pengfei.Study on constant false alarm rate detection approaches for airborne early warning radar［D］.Changsha:National Unversity of Defense Technology，2003.

［15］ 高 飛，謝文沖，段克清，等.共形陣機載相控陣雷達統(tǒng)一雜波建模與分析［J］.國防科技大學(xué)學(xué)報，2008，30(6):94-100.Gao Fei，Xie Wenchong，Duan Keqing，et al.Modeling and analysis of conformal array airborne phased radar clutter［J］.Journal of National University of Defense Technology，2008，30(6):94-100.

［16］ 吳順君，梅曉春，等.雷達信號處理和數(shù)據(jù)處理技術(shù)［M］.北京:電子工業(yè)出版社，2008.Wu Shunjun，Mei Xiaochun，et al.Radar signal processing and data processing technique［M］.Beijing:Publishing House of Electronics Industry，2008.

［17］ 何 友，修建娟，張晶煒，等.雷達數(shù)據(jù)處理及應(yīng)用［M］.北京:電子工業(yè)出版社，2006.He You，Xiu Jianjuan，Zhang Jingwei，et al.Radar data processing and applications［M］.Beijing:Publishing House of Electronics Industry，2006.