(中國人民解放軍92941部隊93分隊,遼寧葫蘆島125001)

0 引言

雷達系統(tǒng)模擬的核心是建立雷達目標回波信號及各種雜波信號散射、傳播特性的模型,從而對雷達的性能和指標進行評估和測試。雜波的散射特性直接影響雷達對目標的檢測和跟蹤性能,例如雜波的功率譜特性與雷達的運動目標顯示及空時自適應(yīng)信號處理(STAP)濾波器的設(shè)計密切相關(guān),雜波的幅度統(tǒng)計特性是雷達恒虛警檢測器設(shè)計的重要依據(jù)[1-2]。對于彈載雷達導引頭來說,面臨一個很重要的背景環(huán)境——海雜波。

目前對海雜波的建模主要集中在物理散射機理建模和統(tǒng)計特性建模上。物理散射機理建模主要研究電磁波與海面的相互作用以及海面雜波單元的散射機理,能夠從本質(zhì)上對雜波散射特性進行解釋,缺點是模型復雜,計算量大,不易形成快速算法。統(tǒng)計特性建模建立在大量的實測數(shù)據(jù)的基礎(chǔ)之上,是一種經(jīng)驗?zāi)Ｐ蚚3]。雜波仿真的實質(zhì)是產(chǎn)生服從特定概率密度分布的相關(guān)序列,常用的方法有零記憶非線性法(ZMNL)和球不變的隨機過程法(SIRP)。

本文結(jié)合反艦導彈不同的作戰(zhàn)背景,提出了一種基于K分布的海雜波回波信號模擬的方法。首先采用SIRP法對海雜波的幅度分布進行建模,功率譜模型采用高斯譜,然后將產(chǎn)生的K分布雜波序列疊加到PD雷達的發(fā)射波形中,合成海雜波的回波信號,最后給出了仿真結(jié)果和實驗分析。

1 PD雷達導引頭海雜波信號建模

1.1 海雜波回波信號的數(shù)學模型

常用的PD雷達導引頭發(fā)射信號有線性調(diào)頻信號,相參脈沖串信號以及二進制相位編碼信號。本文采用的是線性調(diào)頻信號,它的時域表達式為

式中,f0為發(fā)射信號中心頻率,表示脈寬為τ的矩形脈沖函數(shù),μ為線性調(diào)頻系數(shù),其值為μ=B/τ。

雷達回波信號產(chǎn)生的實質(zhì)是雷達的發(fā)射信號作用于目標并發(fā)生后向散射,可以認為回波信號是發(fā)射信號乘以幅度衰減因子A(t)并附加一定的延遲τ0和隨機相位?(t)而得到。數(shù)學表達式為

由式(1)和式(2)可得

地、海雜波可以認為是面雜波,云雨、氣象雜波則是體雜波。無論是面雜波還是體雜波,都是基于點目標來模擬的。單個雷達分辨單元內(nèi)的目標可以看作點目標來處理。單個點目標的回波信號可以表示為

式中,A ij可以通過雷達方程得到,它的具體值為

式中,r ij為雜波單元到雷達之間的距離,σ0S ij為雜波單元的雷達橫截面積(σ0是雜波的后向散射系數(shù),S ij是雜波單元面積),G ij為雜波單元對應(yīng)的天線增益,n ij為雜波的調(diào)制系數(shù)。

1.2 導彈軌跡方程的確定

比例導引法是現(xiàn)代導彈制導中較常見的一種方法,它是指導彈在攻擊目標的過程中,導彈速度矢量的旋轉(zhuǎn)角速度與目標線的旋轉(zhuǎn)角速度成固定比例的一種導引方法。它克服了直接追蹤法彈道彎曲和平行接近法對控制系統(tǒng)要求嚴格的缺點,在工程上容易實現(xiàn)。通過目標飛行的特點選擇合適的導引系數(shù)k,就可以實現(xiàn)較高的追蹤效率。對于不同機動特性的目標適應(yīng)能力也比較強,因此廣泛應(yīng)用于各種類型的導彈上,它的原理結(jié)構(gòu)圖如圖1所示。

圖1中,A為導彈的位置,M為目標的位置,vm為導彈的速度,vt為目標的速度。比例導引法制導的相對運動方程組為

式中,r為目標與導彈的距離,ηt,η分別為目標和導彈速度矢量的前置角,k為比例系數(shù),σt為航向角,σ為導彈導引角,θ為目標線角。

圖1 比例導引示意圖

1.3 仿真場景的設(shè)置

本文所設(shè)置的仿真場景是彈載導引頭在向海面目標飛行的過程中,海雜波對導引頭上雷達回波信號的影響。

1.3.1 波束單元照射面積的劃分

對面雜波的仿真,通常是將雷達波束照射單元劃分成許多小的距離單元,每一個小的距離單元反射的電磁波強度、方位均不相同,對這類波束照射區(qū)域為面雜波時,采用經(jīng)典而又高效的處理方法——網(wǎng)格映像法[4]。網(wǎng)格映像法是根據(jù)雷達的距離分辨率和方位分辨率來將照射區(qū)域劃分成許多獨立的小網(wǎng)格,每一個網(wǎng)格可以認為是點目標,計算每一個網(wǎng)格的回波信號,最后再將所有網(wǎng)格內(nèi)的回波信號進行相干疊加,得到總的回波信號。原理圖如圖2(a)、(b)所示。

圖2(a)、(b)中,VR是雷達平臺的速度,φ是網(wǎng)格單元的擦地角,θ是方位角,Δr是距離分辨率,Δθ是方位角的分辨率,R是雷達到網(wǎng)格單元的距離,H是雷達平臺的高度。根據(jù)文獻[5]可知,網(wǎng)格單元的距離分辨率和方位分辨率分別為

圖2 網(wǎng)格映像法原理圖

式中,τ為雷達發(fā)射信號的脈寬,θmax為雜波散射單元的最大方位角,Δf為多普勒分辨率,Δf=,其中fr是脈沖重復頻率,N是雷達發(fā)射脈沖串個數(shù)。

1.4 雜波功率譜特性和幅度統(tǒng)計特性的計算

目前,對雜波統(tǒng)計建模的研究主要為統(tǒng)計特性和功率譜特性。雜波是隨機起伏信號,其功率譜特性可以用功率譜密度或自相關(guān)函數(shù)來表示。常用的雜波功率譜模型有高斯模型、Cauchy模型,以及全極點模型。本文采用高斯譜來進行仿真,它的表達式為

式中,f0為頻譜最大點的位置,通常令f0=0;f3dB為兩個半功率點的波束寬度,它由海情及雷達的波長所決定;a的取值滿足S(f3dB/2)=0.5,解得

雜波的幅度統(tǒng)計特性與雷達的分辨率、波束入射角等因素有關(guān)。常用的雜波幅度模型有瑞利分布、對數(shù)正態(tài)分布、Weibull分布和K分布等。這些幅度統(tǒng)計模型與雷達平臺參數(shù)、海情、入射角和極化等因素有關(guān),適用范圍如表1所示。

表1 雜波幅度模型

產(chǎn)生雜波序列的方法一般有兩種:零記憶非線性法(ZMNL)和球不變隨機過程法(SIRP)。由于SIRP受所求序列的階數(shù)及自相關(guān)函數(shù)的限制,同時這種方法的計算量非常大,不易形成快速算法,而ZMNL方法克服了這種缺陷。但ZMNL中輸入序列和輸出序列存在著復雜的非線性關(guān)系,文獻[6]對此進行了很好的總結(jié),使得這種非線性關(guān)系已經(jīng)明確。本文采用ZMNL產(chǎn)生相干相關(guān)的K分布,它的原理圖如圖3所示。

圖3 ZMNL產(chǎn)生K分布示意圖

z i的幅度分布特性由非線性變換(ZMNL)獲得,功率譜特性通過濾波器H(w)獲得。使用ZMNL法時應(yīng)先求得z i和w i相關(guān)函數(shù)之間的非線性關(guān)系,然后再求得v i,即通過z i的相關(guān)系數(shù)s ij求得v i的相關(guān)系數(shù)ρij,然后再確定濾波器H(w),這里

H(w)可由下式得到:

H(w)得到后可以通過頻域法得到{v i}。

2 仿真實驗

海雜波回波信號產(chǎn)生的流程如圖4所示。

圖4 海雜波回波信號產(chǎn)生的流程

具體的參數(shù)值如表2所示。

表2 仿真參數(shù)設(shè)置

圖5(a)、(b)分別為線性調(diào)頻信號的時域波形和頻譜圖,可以看出線性調(diào)頻信號的頻譜類似與矩形窗函數(shù),這是因為頻率的線性變化,使得能量平均分布在各個頻點上。

圖6為仿真的K分布與理論K分布的對比,K分布的形狀參數(shù)為v=2,尺度參數(shù)a=0.5;可知,生成的海雜波PDF與理論值大致相吻合,說明了用球不變的隨機過程法實現(xiàn)雜波模擬的有效性。圖7為仿真雜波的功率譜與理論值對比,功率譜為高斯譜,f3dB=40 Hz,模擬產(chǎn)生雜波的功率譜采用Burg法估計,可知仿真的結(jié)果與理論值大致相吻合。

圖8 為產(chǎn)生的K分布雜波序列的實部和虛部,雜波序列長度為4 096點。

仿真產(chǎn)生的海雜波回波信號及其頻譜如圖9(a)、(b)所示,圖9(a)仿真的是50μs內(nèi)海雜波的回波信號,圖9(b)是回波信號的頻譜,可知,幅度在-55～-70 dB之間。

圖5 線性調(diào)頻信號時域與頻譜圖

圖6 幅度分布與理論分布對比

圖7 功率譜與理論對比

圖8 仿真產(chǎn)生雜波序列的實部和虛部

3 結(jié)束語

本文對PD雷達導引頭所面臨的背景環(huán)境——海雜波進行了仿真建模和頻譜分析,對海雜波的一些散射特性,如后向散射系數(shù)、幅度統(tǒng)計特性和頻譜特性進行了仿真分析,仿真實驗得出的結(jié)果與理論值相符。最后對海雜波的回波信號進行了仿真分析,本文所設(shè)置的發(fā)射信號比較簡單,雷達發(fā)射信號的參數(shù),如脈沖重復頻率、帶寬和調(diào)頻系數(shù)等對海雜波回波信號影響沒有進行深入討論,這也是后續(xù)工作所要研究討論的地方。本文的仿真模型和實驗結(jié)果可以作為海雜波中目標信號檢測的基礎(chǔ)。

圖9 海雜波回波信號及其頻譜圖

[1]WARD K D,WATTS S.Use of Sea Clutter Models in Radar Design and Development[J].IET Radar,Sonar and Navigation,2010,4(2):146-157.

[2]姜斌.地、海雜波建模及目標檢測技術(shù)研究[D].長沙:國防科學技術(shù)大學,2004.

[3]楊俊嶺,呂韶昱,萬建偉.一種新的相干K分布模型及其在海雜波仿真中的應(yīng)用[J].系統(tǒng)仿真學報,2007,19(2):250-253,260.YANG Jun-ling,LV Shao-yu,WAN Jian-wei.A New Coherent K-Distributed Model and Its Application in Sea-Clutter Simulation[J].Journal of System Simulation,2007,19(2):250-253,260.(in Chinese)

[4]李芾.基于實測海雜波數(shù)據(jù)的統(tǒng)計特性分析[D].西安:西安電子科技大學,2013.

[5]謝洪森,鄒鯤,周鵬.低掠射角海雜波的統(tǒng)計特性分析[J].雷達科學與技術(shù),2011,9(2):172-179.XIE Hong-sen,ZOU Kun,ZHOU Peng.Statistical Analysis of Sea Clutter at Low Grazing Angle[J].Radar Science and Technology,2011,9(2):172-179.(in Chinese)

[6]MOREAU I.Study and Simulation of Sea Clutter[C]∥IEEE Signal Processing Workshop on Higher-Order Statistics,South Lake Tahoe,CA,USA:IEEE,1993:178-181.