時空相關(guān)海雜波序列的仿真設(shè)計

趙嘉斐，吳敬沄

（1.北京理工大學(xué)北京100081；2.北京十一學(xué)校北京100039）

海面環(huán)境對雷達(dá)信號散射的電磁波稱為海雜波，海雜波對雷達(dá)的目標(biāo)信號產(chǎn)生干擾，因而在雷達(dá)系統(tǒng)的建?；蚶走_(dá)特性的分析階段都需要產(chǎn)生海雜波序列用于模擬雷達(dá)的真實回波[1-2]。對海雜波序列的研究主要包括海雜波的建模與海雜波的仿真。海雜波的建模研究又包括對雜波幅度均值特性、幅度分布特性與相關(guān)性的建模[3]。本文基于海雜波的幅度均值特性、幅度分布特性與相關(guān)性經(jīng)驗?zāi)Ｐ蛯ｋs波進(jìn)行仿真設(shè)計，生成了4種服從不同分布的相關(guān)海雜波序列。

1 雜波幅度均值特性模型

常用的海雜波幅度均值特性模型[4]包括：GIT模型、HYB模型與TSC模型。海雜波的幅度均值特性與頻率、極化、擦地角、風(fēng)速、風(fēng)向及浪高等參數(shù)有關(guān)，用于表征海雜波的平均回波強(qiáng)度。模型的獲取基于海雜波實測數(shù)據(jù)。

1）GIT模型

20世紀(jì)70年代，GIT模型由美國喬治理工學(xué)院工作人員提出。GIT模型覆蓋了10～100 GHz的雷達(dá)頻率范圍，不同頻段對應(yīng)模型系數(shù)不同，本文只給出10 GHz以上的模型參數(shù)。GIT模型在擦地角0.1°～10°范圍、風(fēng)向視角0°～180°的范圍內(nèi)有效。GIT模型把海面平均雜波散射系數(shù)建模為擦地角、風(fēng)速、平均波浪高度、風(fēng)向視角、頻率以及極化方式的函數(shù)。GIT模型的一個重要特征是使用了波浪高度和風(fēng)速參數(shù)來綜合描述海面。該模型認(rèn)為，海雜波是多路徑因子、風(fēng)向因子和風(fēng)速因子3個因子作用的結(jié)果，其中多路徑因子是由高斯分布波浪高度經(jīng)過理論推導(dǎo)得出，其余2個因子則由測量數(shù)據(jù)用經(jīng)驗公式擬合而得出。

GIT模型HH極化和VV極化的歸一化后向散射系數(shù)σ0分別滿足以下表達(dá)式：

式中，?為雷達(dá)擦地角，以弧度制為單位；Al、Au、Aw分別為多路徑因子、風(fēng)向因子和風(fēng)速因子；hav為平均浪高，單位為m。

多路徑因子Al為用于估計多路徑或干涉的參數(shù)，取值如下：

λ為雷達(dá)波長，單位為m。

風(fēng)向因子Au為與風(fēng)向相關(guān)的參數(shù)，取值為：

θ為雷達(dá)觀測方向與逆風(fēng)方向的夾角。

風(fēng)速因子Aw為與風(fēng)速有關(guān)的因子，滿足以下表達(dá)式：

對于充分發(fā)展的海浪，風(fēng)速與平均浪高之間滿足以下關(guān)系：

2）TSC模型

TSC模型在結(jié)構(gòu)上與GIT模型類似，但TSC模型還包含了大氣波導(dǎo)傳播下采集的數(shù)據(jù)，更為全面。TSC模型頻率范圍覆蓋0.5～35 GHz，不同頻段對應(yīng)模型系數(shù)不同，本文只給出10 GHz以上的模型參數(shù)。TSC模型的風(fēng)向視角覆蓋了0°～180°的范圍，擦地角范圍覆蓋0°～90°。與GIT模型不同的是，TSC模型采用Douglas海況參數(shù)來描述海面，可以利用Douglas海況參數(shù)來計算波浪高度和風(fēng)速參數(shù)，但不能獨立輸入波浪高度和風(fēng)速參數(shù)。該模型不像GIT模型那樣平均散射系數(shù)隨著距離的增加（即擦地角的減少）而迅速衰減，原因是在模型推導(dǎo)時對實測數(shù)據(jù)進(jìn)行了修正，同時把實測數(shù)據(jù)看作全方位視角回波值的平均，即把它近似看作橫向視角回波值。

TSC模型的輸入是雷達(dá)頻率、極化方式、Douglas海況等級、擦地角和風(fēng)向角，而風(fēng)速和浪高是海況的函數(shù)，并不作為獨立的變量輸入。需要注意的是國際標(biāo)準(zhǔn)海況等級與Douglas海況等級有較大差別，實際應(yīng)用中不應(yīng)將兩者混淆。Douglas海況對應(yīng)參數(shù)如表1所示。

表1 Douglas海況等級

TSC模型認(rèn)為海面后向散射系數(shù)與擦地角因子、風(fēng)速因子和風(fēng)向因子3項因子有關(guān)。模型具體表示如下：

式中，?為雷達(dá)擦地角，以弧度制為單位；Gl、Gu、Gw分別為擦地角因子、風(fēng)向因子和風(fēng)速因子；λ為雷達(dá)波長，單位為m。

擦地角因子Gl為與擦地角、路徑有關(guān)的參數(shù)，取值如下：

式中，σz為海面高度標(biāo)準(zhǔn)差，單位為m；S為Douglas海況參數(shù)。

風(fēng)速因子Gw，滿足以下表達(dá)式：

式中，U為風(fēng)速，單位為m/s。參數(shù)A滿足以下關(guān)系式：

風(fēng)向因子Gu是與風(fēng)向相關(guān)的參數(shù)，取值為：

θ為雷達(dá)觀測方向與逆風(fēng)方向的夾角，?為擦地角，二者均以弧度表示。

3）HYB模型

HYB模型同樣適用于0.5～35 GHz頻率范圍，且不同頻段對應(yīng)模型系數(shù)不同，本文只給出10 GHz以上的模型參數(shù)。HYB模型的風(fēng)向視角也在0°～180°范圍，但HYB模型的擦地角范圍為0.1°～30°。HYB模型也利用Douglas海況參數(shù)來描述海面。不同的是在HYB模型中引入了海況為5、擦地角為0.1°、VV極化以及逆風(fēng)視角情況下的參考散射系數(shù)來計算其他任意情況下的實際散射系數(shù)。HYB模型的散射系數(shù)計算式為：

參考擦地角?ref=0.1°，臨界角為：

式中，?ref為參考擦地角，單位為度；?t為臨界角，單位為度；λ為波長，單位為m；f為雷達(dá)頻率，單位為GHz；σh=0.031S2，為浪高的均方根值，單位為m。

若?t≥?ref，則

若?t＜?ref，則

海況參數(shù)：

極化因子：

對于VV極化，有

對于HH極化，有

式中，hav=0.08S2，為平均浪高，單位為m；?為擦地角，單位為弧度。

風(fēng)向因子：

式中，θ為風(fēng)向視角。

2 雜波幅度分布模型

常用的海雜波幅度分布模型[5]包括：瑞利分布模型、對數(shù)-正態(tài)分布模型、韋布爾分布模型及K分布模型。其中，瑞利分布適用于低分辨率雷達(dá)、高擦地角的情況。對數(shù)-正態(tài)分布模型、韋布爾分布模型及K分布模型適用于高分辨率雷達(dá)、低擦地角、高海情的情況。本文討論了海雜波的經(jīng)驗?zāi)Ｐ图捌溥m用范圍，并在此基礎(chǔ)上通過仿真設(shè)計生成海雜波序列。

1）瑞利分布

瑞利分布[6]是雷達(dá)雜波中最常用也是用得最早的一種統(tǒng)計模型。在雷達(dá)可分辨范圍內(nèi)，當(dāng)散射體的數(shù)目很多時，根據(jù)散射體反射信號振幅和相位的隨機(jī)特性，認(rèn)為所有散射體合成回波的包絡(luò)振幅服從瑞利分布。假設(shè)x表示瑞利分布雜波回波的包絡(luò)振幅，其概率密度函數(shù)表示為：

式中，w決定了瑞利分布的形狀參數(shù)，w取不同值時瑞利分布的概率密度函數(shù)曲線尖銳程度不同，表現(xiàn)在w越小，PDF曲線越尖銳。瑞利分布模型的適用范圍為低分辨率雷達(dá)及高擦地角的情況。

2）韋布爾分布

韋布爾分布[7]含有兩個參數(shù)，瑞利分布是它的一個特例。在較寬擦地角范圍內(nèi)，韋布爾分布比瑞利分布或?qū)?shù)-正態(tài)分布能更好地表征雜波的幅度特性。其概率密度函數(shù)表示為：

式中，r、w分別為韋布爾分布的形狀參數(shù)和尺度參數(shù)。當(dāng)r=2時，韋布爾分布退化為瑞利分布。r越小，PDF曲線越尖銳。韋布爾分布的適用范圍為高分辨率雷達(dá)、低擦地角、高海情條件。

3）對數(shù)-正態(tài)分布

當(dāng)雷達(dá)分辨率較高、擦地角較小時，水平極化的海雜波的幅度通常會表現(xiàn)為對數(shù)-正態(tài)分布[8]，其概率密度函數(shù)表示為：

式中，μ、η2分別為服從正態(tài)分布的lnx的均值和方差，對數(shù)-正態(tài)分布的概率密度曲線由這兩個參數(shù)決定。對數(shù)-正態(tài)分布適用于高分辨率雷達(dá)、低擦地角、高海情條件。

4）復(fù)合K分布

復(fù)合K分布[9-10]模型是由實數(shù)的紋理分量與復(fù)數(shù)的散斑分量的乘積構(gòu)成。調(diào)制生成海雜波的表達(dá)式及其對海表面散射機(jī)理的數(shù)學(xué)描述表述如下。假設(shè)用隨機(jī)變量y服從Gamma分布，表示海雜波的紋理分量，其概率密度函數(shù)表示為：

式中，Γ(v)為伽馬函數(shù)，v、b分別為形狀參數(shù)和尺度參數(shù)。

用隨機(jī)變量x表示海雜波的散斑分量，則x服從瑞利分布，其概率密度函數(shù)表示為：

基于全概率密度函數(shù)，復(fù)合K分布的海雜波的幅度概率密度函數(shù)可以表示為：

3 雜波相關(guān)性模型

1）時間相關(guān)性

通常用功率譜來描述海雜波的時間相關(guān)性[11]，傳統(tǒng)的雜波歸一化功率譜密度函數(shù)有高斯譜、柯西譜及全極譜，其中高斯譜的應(yīng)用最為廣泛，其經(jīng)驗?zāi)Ｐ腿缦拢?/p>

式中，fd為多普勒頻移，也是功率譜的中心頻率，由雷達(dá)與海面的徑向速度決定。s為高斯譜分布的標(biāo)準(zhǔn)差。

2）空間相關(guān)性

海雜波的空間相關(guān)性[12]指徑向相鄰兩雜波距離分辨單元內(nèi)回波信號間的相關(guān)性。海雜波的空間相關(guān)函數(shù)的經(jīng)驗?zāi)Ｐ蚚2]，有：

其中，L為雷達(dá)距離分辨率，m為距離單元個數(shù)。ρ為空間相關(guān)距離，大于該距離的分辨單元不相關(guān)，取值如下。

其中，U為風(fēng)速，θ為雷達(dá)視線與風(fēng)向之間的夾角，g=9.18為重力加速度。

4 仿真實例

基于以上海雜波經(jīng)驗?zāi)Ｐ停疚姆抡嫔闪藭r間相關(guān)瑞利分布序列、時間相關(guān)對數(shù)-正態(tài)分布序列、時間相關(guān)韋布爾分布序列及時間相關(guān)K分布序列[13-16]。仿真結(jié)果如圖1～圖4所示。

圖1 相關(guān)瑞利分布序列概率密度與功率譜

圖2 相關(guān)對數(shù)-正態(tài)分布序列概率密度與功率譜

圖3 相關(guān)韋布爾分布序列概率密度與功率譜

圖4 相關(guān)K分布序列概率密度與功率譜

5 結(jié)束語

相關(guān)海雜波序列的仿真生成對雷達(dá)系統(tǒng)的建模與特性的分析有重要意義。本文首先介紹了3種海雜波幅度均值特性模型及其適用條件，其次介紹了4種海雜波幅度分布經(jīng)驗?zāi)Ｐ图捌鋺?yīng)用場景，最后介紹了海雜波的時間相關(guān)性與空間相關(guān)性模型?；谝陨夏Ｐ筒捎胢atlab語言對海雜波序列進(jìn)行仿真設(shè)計，得到4種不同的相關(guān)海雜波序列仿真結(jié)果圖。