王 旭 蔡興雨 朱思橋 朱永杰

(西安電子工程研究所 西安 710100)

基于點(diǎn)跡特征的氣象雜波抑制技術(shù)

王 旭 蔡興雨 朱思橋 朱永杰

(西安電子工程研究所 西安 710100)

對(duì)空情報(bào)雷達(dá)在復(fù)雜環(huán)境中經(jīng)常遭受各種雜波的影響，主要包括地雜波與氣象雜波等等。地雜波主要會(huì)對(duì)較低空域的目標(biāo)造成掩蓋、產(chǎn)生虛警，一般采用合適的MTI、MTD結(jié)合雜波圖等雜波抑制手段即可以較有效地抑制；而氣象雜波主要會(huì)影響較高空域目標(biāo)的探測(cè)，由于通常處于運(yùn)動(dòng)狀態(tài)、雜波譜較寬等因素，較難抑制。針對(duì)此問(wèn)題，提出一種基于點(diǎn)跡特征的氣象雜波抑制技術(shù)，首先對(duì)已檢測(cè)點(diǎn)跡進(jìn)行幅度歸一預(yù)處理，然后利用實(shí)測(cè)數(shù)據(jù)分析出氣象雜波與目標(biāo)點(diǎn)跡在幅度起伏和速度上具有較大差異，根據(jù)幅度均值、方差和速度等特征優(yōu)化合適的分類器參數(shù)，對(duì)滿足氣象特性的航跡進(jìn)行剔除。結(jié)果驗(yàn)證了所提方法的有效性。

氣象雜波；點(diǎn)跡特征；航跡；歸一化幅度

Abstract: Air intelligence radar is often effected by various clutter in complex environment, including ground clutter and meteorological clutter, and the former will mainly cover target at lower altitude, and generate false alarm, which can be suppressed effectively by using appropriate clutter suppressing means like MTI and MTD combining clutter map; and the latter will effect detecting target at higher altitude which is difficult to be suppressed due to moving status of clutter and wider clutter spectrum. Aiming at this issue, a meteorological clutter suppression technology based on plot characteristic is presented. Firstly, the detected plot is performed amplitude normalization processing, and then using measured data analyzes bigger difference between meteorological clutter and target plot in amplitude fluctuation and speed. An appropriate classifier parameter is optimized according to amplitude average value, variance and speed. The track which met meteorological clutter characteristics is eliminated. The result verifies effectiveness of the presented method.

Keywords：meteorological clutter; plot characteristic; track; normalized amplitude

0 前言

雷達(dá)回波中通常包含有目標(biāo)的反射回波、噪聲以及自然環(huán)境中不需要的反射回波分量。如果目標(biāo)處于雜波背景中，弱目標(biāo)會(huì)淹沒(méi)于強(qiáng)雜波中，造成漏警；強(qiáng)雜波也會(huì)造成虛警，從而產(chǎn)生虛假航跡以及影響真實(shí)目標(biāo)航跡的建立；此外過(guò)多的雜波會(huì)引起終端的過(guò)載，嚴(yán)重影響雷達(dá)的工作性能。因此，雜波抑制是雷達(dá)相關(guān)領(lǐng)域一個(gè)很重要的問(wèn)題[1-6]。

不同應(yīng)用領(lǐng)域的雷達(dá)系統(tǒng)對(duì)雜波的定義可能不同，對(duì)于空情雷達(dá)而言，來(lái)自地物、海面、云雨以及敵人釋放的箔條等均為雜波，可以分為固定雜波與運(yùn)動(dòng)雜波[1-3]。固定雜波一般位于低空區(qū)域，主要會(huì)對(duì)較低空域的目標(biāo)造成掩蓋、產(chǎn)生虛警，一般采用合適的MTI、MTD結(jié)合雜波圖等雜波抑制手段即可以較有效地抑制，對(duì)于中高空目標(biāo)可以采用空域?yàn)V波進(jìn)一步抑制；運(yùn)動(dòng)雜波中的氣象對(duì)雷達(dá)的影響主要表現(xiàn)在兩個(gè)方面，一是電磁波穿過(guò)氣象區(qū)域產(chǎn)生的衰減，二是氣象顆粒對(duì)電磁波散射產(chǎn)生的回波[7]。由于氣象雜波通常處于運(yùn)動(dòng)狀態(tài)而且雜波譜較寬，常規(guī)MTI和MTD技術(shù)無(wú)法有效對(duì)其抑制，目前多采用自適應(yīng)MTI技術(shù)[1-5，8]，即首先估計(jì)出雜波的多普勒頻移，再移動(dòng)MTI濾波器的凹口使之對(duì)準(zhǔn)雜波的多普勒頻移，或者將運(yùn)動(dòng)雜波頻譜補(bǔ)償至固定雜波譜位置上[1]。但當(dāng)雜波多普勒頻移估計(jì)不準(zhǔn)時(shí)，可能抑制效果不理想，此外，抑制該多普勒頻移雜波的同時(shí)，也會(huì)同時(shí)抑制該多普勒頻移附近的目標(biāo)。而由于氣象雜波的運(yùn)動(dòng)，剩余的氣象虛警會(huì)逐漸建立較穩(wěn)定的航跡，而且與目標(biāo)航跡比較相似，以致雷達(dá)操作人員造成誤判，影響后續(xù)的處理。

針對(duì)現(xiàn)有方法無(wú)法有效抑制氣象雜波的問(wèn)題，本文提出一種基于點(diǎn)跡特征的氣象雜波抑制技術(shù)，首先對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行分析，統(tǒng)計(jì)出氣象雜波與目標(biāo)點(diǎn)跡的幅度起伏特性和速度均值，根據(jù)訓(xùn)練樣本幅度均值、方差和速度優(yōu)化合適的分類器參數(shù)，利用推導(dǎo)出的均值和方差更新公式實(shí)時(shí)計(jì)算已建航目標(biāo)的幅度均值、幅度方差和速度，對(duì)判別為氣象的航跡進(jìn)行剔除。結(jié)果驗(yàn)證了所提方法的有效性。

1 氣象雜波

考慮雷達(dá)波束照射某個(gè)方向，波束內(nèi)的云等氣象現(xiàn)象會(huì)產(chǎn)生后向散射。在云中，后向散射體通常指云層中的很多顆粒狀的水滴、冰晶。每個(gè)顆粒即可以看成一個(gè)散射體，雷達(dá)回波信號(hào)即是每個(gè)散射體回波的矢量和。不同類型的云所含有的水滴與冰晶密度不同，因此散射特性有所區(qū)別，其有效截面積可以表示為[2，7]

σc=σ0V

(1)

由于風(fēng)速的影響，雷達(dá)雜波單元內(nèi)各散射體具有一定的速度，因此雜波的功率譜也具有一定的分布，通常對(duì)于氣象雜波的譜特性一般符合高斯分布模型[1，2，7]，即

(2)

其中Pc為雜波功率，σf為雜波譜寬的標(biāo)準(zhǔn)偏差，fdc表示雜波的多普勒頻移，氣象雜波的多普勒譜如圖1所示。

圖1 固定雜波與氣象雜波譜特性示意圖

`從圖中可以看出氣象雜波的中心頻率不在零頻，因此采用常規(guī)MTI或MTD無(wú)法有效進(jìn)行抑制。由于氣象雜波的運(yùn)動(dòng)，剩余的氣象虛警會(huì)逐漸建立較穩(wěn)定的航跡，如圖2所示，而且與目標(biāo)航跡比較相似，以致雷達(dá)操作人員造成誤判，影響后續(xù)的處理。目前較有效的方式是采用自適應(yīng)MTI，首先估計(jì)出雜波的多普勒頻移，再對(duì)MTI濾波器系數(shù)乘以一組頻移系數(shù)，將凹口對(duì)準(zhǔn)雜波的中心頻率，或者將運(yùn)動(dòng)雜波頻譜補(bǔ)償至固定雜波譜位置上。由于多普勒頻率僅僅為點(diǎn)跡中的一種特征，當(dāng)雜波多普勒頻移估計(jì)不準(zhǔn)時(shí)，有可能抑制效果不理想，此外，抑制該多普勒頻移雜波的同時(shí)，也會(huì)同時(shí)抑制該多普勒頻移附近的目標(biāo)。因此考慮采用點(diǎn)跡的多種特征對(duì)氣象雜波進(jìn)行抑制。

圖2 氣象雜波形成的虛假航跡

2 氣象雜波抑制

雷達(dá)檢測(cè)出的點(diǎn)跡通常具有幅度、速度、距離、方位等信息。而不同的目標(biāo)、雜波往往分布在不同的空間位置上，具有不同的幅度大小、起伏特性，不同的速度大小，或者在空間上占據(jù)不同的大小。如果聯(lián)合點(diǎn)跡的這些特征或者部分特征，相信可以提高氣象雜波抑制性能。下面采用點(diǎn)跡的幅度均值、方差與速度特征來(lái)抑制氣象雜波航跡。

2.1 雷達(dá)方程與預(yù)處理

設(shè)雷達(dá)的發(fā)射功率為Pt，雷達(dá)發(fā)射天線增益為Gt，接收天線增益為Gr，信號(hào)波長(zhǎng)為λ，目標(biāo)或雜波散射截面積為σ，其距離雷達(dá)的距離為R，則雷達(dá)接收端的回波功率可以表示為[1]

(3)

由于不同的雷達(dá)截面積、不同的氣象和不同的距離將引起不同的回波強(qiáng)度，雷達(dá)通常接收幅度變化很大的回波，超過(guò)固定增益接收機(jī)的動(dòng)態(tài)范圍[1，3]。而增加接收機(jī)動(dòng)態(tài)范圍，會(huì)大為增加接收機(jī)的設(shè)計(jì)難度。為了壓縮接收機(jī)的動(dòng)態(tài)范圍，并防止近程強(qiáng)雜波使得接收機(jī)過(guò)載飽和，通常采用近程增益控制(STC)，STC使雷達(dá)接收機(jī)的靈敏度隨時(shí)間(距離)變化，在近距離衰減的很大，在遠(yuǎn)距離時(shí)保持原來(lái)的增益和靈敏度。因此經(jīng)過(guò)信號(hào)處理后檢測(cè)出的第i個(gè)點(diǎn)跡的幅度可以表示為

(4)

其中Gp表示接收端信號(hào)處理的得益，F(xiàn)stc(R)表示STC控制規(guī)律，其為距離R的函數(shù)?？梢钥闯鳇c(diǎn)跡幅度Ai為距離Ri、目標(biāo)散射截面積σi以及STC控制規(guī)律Fstc(Ri)的函數(shù)。對(duì)于檢測(cè)出的目標(biāo)點(diǎn)跡，距離Ri已知，因此為了消除距離對(duì)點(diǎn)跡幅度的影響，將所有檢測(cè)出的點(diǎn)跡歸一至相同的距離R0，即

(5)

2.2 點(diǎn)跡特征分析

從實(shí)測(cè)數(shù)據(jù)中選擇已建航的清潔區(qū)域的目標(biāo)點(diǎn)跡作為目標(biāo)的訓(xùn)練樣本，從無(wú)目標(biāo)影響時(shí)所錄的實(shí)測(cè)數(shù)據(jù)中選擇合適的氣象雜波點(diǎn)跡作為氣象雜波訓(xùn)練樣本。按照式(5)分別對(duì)其航跡中的點(diǎn)跡幅度做歸一化預(yù)處理，得到相應(yīng)的歸一化幅度，然后統(tǒng)計(jì)其幅度均值和方差。

圖2-圖4所示為從多組樣本中包含的點(diǎn)跡信息，圖5為多組樣本中統(tǒng)計(jì)出的各特征，從圖中可以看出氣象雜波的幅度較弱、起伏小，而且速度也比目標(biāo)小。其余樣本也具有類似結(jié)果，因此考慮對(duì)于已建航的目標(biāo)，選擇幅度均值、幅度標(biāo)準(zhǔn)差及速度均值作為特征進(jìn)行分類。

2.3 均值與方差更新公式

由于雷達(dá)波束掃描一次，目標(biāo)航跡和虛假航跡會(huì)不斷有新點(diǎn)跡出現(xiàn)，從而引起幅度均值與方差的變化，這要求能夠?qū)崟r(shí)地計(jì)算均值與方差。因此下面推導(dǎo)了均值與方差的更新公式，即用上一時(shí)刻的均值、方差以及當(dāng)前的幅度來(lái)計(jì)算當(dāng)前時(shí)刻的均值與方差。

圖3 氣象航跡的點(diǎn)跡特征

設(shè)第n個(gè)時(shí)刻的樣本表示為xn，第n時(shí)刻得到的均值為mn，方差為vn。則均值的更新公式可以直接表示為：

(6)

第n時(shí)刻的方差可以表示為：

從而有

因此方差更新公式為

(7)

2.4 處理方式

根據(jù)以上分析，基于點(diǎn)跡特征的氣象雜波抑制技術(shù)可總結(jié)如下：

a)對(duì)已檢測(cè)點(diǎn)跡的幅度按照式(5)進(jìn)行歸一化預(yù)處理；

b)數(shù)處按照式(6)、式(7)對(duì)建航的目標(biāo)統(tǒng)計(jì)其幅度均值、方差和速度均值；

c)根據(jù)幅度均值、方差和速度均值，按照一定的準(zhǔn)則對(duì)航跡進(jìn)行分類；

圖4 目標(biāo)的點(diǎn)跡特征

d)終端根據(jù)數(shù)處的統(tǒng)計(jì)結(jié)果及顯示使能按鈕，對(duì)可能的氣象雜波的航跡進(jìn)行區(qū)別顯示或不顯示。

此外，對(duì)于剛建航的目標(biāo)(點(diǎn)跡較少的情況)，判為氣象雜波的可疑航跡，可以標(biāo)記為可疑。

2.4.1 線性分類器

通過(guò)圖6可以看出兩類數(shù)據(jù)在特征空間下具有線性可分性，因此可以采用如下的線性分類器[9]

(8)

其中αi為松弛項(xiàng)，xi∈RN×1表示第i個(gè)訓(xùn)練樣本，N表示特征的維數(shù)，li∈{1，-1}，L表示樣本數(shù)。由于線性分類器不具有拒判能力，因此結(jié)合門限濾波的方式(對(duì)線性分類器分類結(jié)果中判為氣象的樣本進(jìn)行判決，低于門限的判為氣象，予以濾除)，同時(shí)滿足條件1與條件2，判為氣象航跡

condition1：wTx+b≤0

(9)

condition2：xn≤Thn，n=1，…，N

根據(jù)訓(xùn)練樣本，優(yōu)化得到如下參數(shù)

w=[-0.124 -0.065 -0.341]T，b=30.704，Th=[88.932 8.360 60.442]T。

2.4.2 一類分類器

考慮到一類分類器對(duì)樣本數(shù)據(jù)具有較好的描述，可以減小庫(kù)外樣本被錯(cuò)分的可能，因此設(shè)計(jì)了基于最小化超球體積[10]的一類分類器

圖5 點(diǎn)跡特征統(tǒng)計(jì)

圖6 兩類樣本在特征空間中的可分性

(10)

其中xi表示第一類訓(xùn)練樣本，xj表示第二類訓(xùn)練樣本，x0表示第一類樣本的均值，αi為松弛項(xiàng)。對(duì)滿足如下條件的樣本，判為氣象航跡。

(x-x0)TA(x-x0)≤1

(11)

根據(jù)訓(xùn)練樣本，優(yōu)化得到如下參數(shù)

x0=[69.727 4.431 28.988]T，

Th=[69.727 4.431 28.988]T，

2.4.3 門限判決方式

通過(guò)圖5與圖6可以看出，氣象航跡具有較小的幅度均值和方差，且速度較小，因此考慮采用設(shè)置門限的方式，對(duì)滿足如下條件的航跡判為氣象航跡。

x≤Th

根據(jù)訓(xùn)練樣本得到的門限如下

Th=[86.391 7.956 58.715]T。

圖7所示為采用三種方式得到的分類結(jié)果，可以看出三種分類方式能夠有效的對(duì)氣象航跡與目標(biāo)航跡進(jìn)行區(qū)分。線性分類器計(jì)算量較小，但可能由于訓(xùn)練數(shù)據(jù)不足等原因，造成錯(cuò)分；一類分類器計(jì)算量與其它兩類相比稍大，但拒判性能較好；門限判決方式的計(jì)算量小，可能存在錯(cuò)分，應(yīng)盡量將門限選小。

3 實(shí)測(cè)數(shù)據(jù)處理結(jié)果

按照第2.2節(jié)所說(shuō)明的方式選擇實(shí)測(cè)數(shù)據(jù)作為訓(xùn)練樣本。圖8所示為訓(xùn)練數(shù)據(jù)中的一組，從圖8.a中可以看出，氣象雜波航跡在空間上許多區(qū)域均有分布，航跡方向較一致，這是由于氣象雜波的運(yùn)動(dòng)受風(fēng)向、風(fēng)速影響的原因。與圖8.b中目標(biāo)航跡相比，非常相似，因此嚴(yán)重影響了雷達(dá)的正常工作性能以及雷達(dá)操作人員的決策。

圖7 分類結(jié)果

圖8 作為訓(xùn)練用的實(shí)測(cè)數(shù)據(jù)

按照第2.4節(jié)的處理方式，根據(jù)訓(xùn)練樣本的幅度均值、標(biāo)準(zhǔn)差和速度均值等特征，設(shè)定合適的參數(shù)，并應(yīng)用于受氣象雜波污染的實(shí)測(cè)數(shù)據(jù)。圖9所示為氣象雜波航跡剔除前后的對(duì)比圖(采用線性分類器，其他兩種方式效果類似)，對(duì)比圖9.a和圖9.b，可以看出，目標(biāo)飛行區(qū)域中的航跡得到較好的保留，而氣象雜波得到了較明顯的抑制。

圖9 氣象雜波抑制前后的對(duì)比結(jié)果

4 結(jié)束語(yǔ)

目前針對(duì)氣象雜波的抑制方法多基于雜波的多普勒特性，但當(dāng)雜波多普勒頻移估計(jì)不準(zhǔn)時(shí)，雜波抑制性能會(huì)嚴(yán)重下降，而且虛警建立的氣象航跡與目標(biāo)航跡相似，嚴(yán)重影響雷達(dá)的正常工作性能。針對(duì)氣象雜波抑制問(wèn)題，本文結(jié)合雜波速度均值與幅度起伏特性，提出一種基于點(diǎn)跡特征的氣象雜波抑制技術(shù)，首先利用實(shí)測(cè)數(shù)據(jù)分析出氣象雜波與目標(biāo)點(diǎn)跡在幅度起伏特性和速度上具有較大差異，根據(jù)訓(xùn)練樣本中幅度均值、方差和速度等特征，得到合適分類器參數(shù)，對(duì)滿足氣象特性的航跡進(jìn)行剔除。本文方法對(duì)于點(diǎn)跡特征的選擇，只考慮了幅度的均值、方差和速度，而氣象雜波在空間上通常具有一定的擴(kuò)展型，因此結(jié)合點(diǎn)跡的空間分布特性可能可以更好地提高抑制效果，這有待進(jìn)一步研究。

[1] 丁鷺飛，耿富錄，陳建春.雷達(dá)原理[M].北京：電子工業(yè)出版社，2014.

[2] Skolnik Merrill I，左群生，徐國(guó)良等譯.雷達(dá)系統(tǒng)導(dǎo)論(第3版)[M].北京：電子工業(yè)出版社，2014.

[3] Skolnik Merrill I.雷達(dá)手冊(cè)(第3版)[M]. 北京：電子工業(yè)出版社，2010.

[4] Candan Cagatay，Yilmaz Ozgur A. Efficient Methods of Clutter Suppression for Coexisting Land and Weather Clutter Systems[J]. IEEE Transactions on Aerospace and Electronic Systems. 2009，45(4)：1641-1650.

[5] Short Robert D. An Adaptive MTI for Weather Clutter Suppression[J]. IEEE Transactions on Aerospace and Electronic Systems.1982，18(5)：552-562.

[6] 葛鳳翔，孟華東，彭應(yīng)寧. 雜波譜中心和譜寬估計(jì)方法[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版). 2002，42(7)：941-944

[7] 楊尚軍，雷達(dá)模擬器氣象雜波的建模與仿真[D]. 國(guó)防科技大學(xué)碩士學(xué)位論文. 2006：33-39.

[8] 伊銳，唐緯. S波段兩坐標(biāo)雷達(dá)反氣象的一種方法[J]. 現(xiàn)代雷達(dá). 2001，23(增)：34-36.

[9] 李彥兵，基于微多普勒效應(yīng)的運(yùn)動(dòng)車輛分類研究[D]. 西安電子科技大學(xué)博士學(xué)位論文. 2012：57-62.

[10] Boyd S.，Vandenberghe L. Convex Optimization[M]. Cambridge，UK：Cambridge Univ. Press，2004：221-229，422-432.

TechnologyofMeteorologicalClutterSuppressionBasedonPlotCharacteristic

Wang Xu， Cai Xingyu， Zhu Siqiao， Zhu Yongjie

(XI’an Electronic Engineering Research Institute， Xi’an 710100)

TN95

A

1008-8652(2017)02-006-09

2017-04-25

王 旭(1987-)，男，博士研究生。研究方向?yàn)槔走_(dá)系統(tǒng)、信號(hào)處理、發(fā)射波形設(shè)計(jì)。