韓壯志+宋春吉+李欽+陳燕

摘 要： 針對連發(fā)火炮的復(fù)雜多目標環(huán)境，傳統(tǒng)的單元平均恒虛警處理方法不能滿足信號檢測的要求。為了充分利用目標的運動信息，通過距離?多普勒二維恒虛警方法估計背景噪聲，采用修正的單元平均恒虛警處理剔除目標之間的相互干擾，并結(jié)合雙門限檢測對連發(fā)火炮多目標回波進一步處理，引入動態(tài)閾值增強檢測的魯棒性。工程實踐表明，該方法提高了回波信號的信噪比和信雜比，有效地實現(xiàn)了連發(fā)火炮多目標的檢測，為下一步目標參數(shù)提取提供基礎(chǔ)。

關(guān)鍵詞： 多目標； 距離?多普勒二維； 修正單元平均恒虛警； 雙門限檢測

中圖分類號： TN955+.1?34 文獻標識碼： A 文章編號： 1004?373X（2017）15?0006?04

Abstract： For the complex multi?target environment of continuous?emission artillery， the CA?CFAR processing method can′t meet the requirement of signal detection. In order to make full use of the moving information of the target， the distance?Doppler two?dimensional CFAR method is used to estimate the background noise. The modified CA?CFAR processing is adopted to eliminate the mutual interference among the targets. The dual?threshold detection is combined to further process the multi?target echo of continuous?emission artillery. The dynamic threshold is introduced to enhance the detection robustness. The engineering practice result shows that the method can improve the signal?to?noise ratio and signal?to?clutter ratio of the echo signal， realize the continuous?emission artillery multi?target detection effectively， and provide the basis for the target parameter extraction of next step.

Keywords： multi?target； distance?Doppler two?dimension； modified CA?CFAR； dual?threshold detection

0 引 言

恒虛警處理是雷達信號檢測的關(guān)鍵。雷達信號檢測是在各種噪聲和雜波干擾的環(huán)境中進行的[1]，雜波的分布特性通常是未知的，選擇恰當?shù)暮闾摼幚矸椒軌蚝芎玫毓烙嬰s波功率，實現(xiàn)目標檢測。經(jīng)典的恒虛警處理方法有均值類恒虛警處理（ML?CFAR）和有序統(tǒng)計恒虛警處理（OS?CFAR）[2?4]，但是針對復(fù)雜的多目標環(huán)境都存在不同程度的局限性。

由于彈丸相對雷達徑向變速運動，所以目標在時頻圖上的信息是一條斜線。本文考慮運用距離?多普勒二維信息對背景噪聲進行估計，逐頻率點進行時間維上的恒虛警。為了避免多目標之間的相互影響，對傳統(tǒng)單元平均恒虛警進行修正減少目標干擾[5]。結(jié)合雙門限檢測算法，極大地降低了回波信號的虛警率。工程實踐表明，該方法可以實現(xiàn)連發(fā)火炮多目標的檢測。

1 恒虛警處理在連發(fā)火炮多目標檢測中的改進

測速雷達發(fā)射偽隨機碼調(diào)相連續(xù)波信號進行火炮彈丸探測，設(shè)定雷達與火炮位置的關(guān)系如圖1所示。

本文利用傳統(tǒng)的檢測后跟蹤步驟實現(xiàn)彈丸的跟蹤濾波。首先，獲取單個距離門回波數(shù)據(jù)并進行雜波與噪聲抑制；其次，經(jīng)過干擾抑制后的數(shù)據(jù)進行目標檢測，一般為固定門限或恒虛警檢測；最后，單距離門數(shù)據(jù)檢測后的準彈跡點進行距離門之間相互關(guān)聯(lián)，形成多距離門航跡。適當?shù)暮闾摼幚矸椒ㄊ菍γ恳粋€距離門內(nèi)數(shù)據(jù)進行目標檢測的關(guān)鍵。恒虛警處理方法通過對參考單元的功率電平按照一定準則得到待測單元雜波功率電平的近似估計，根據(jù)規(guī)定的虛警率計算檢測門限。不斷改進的恒虛警處理方法的核心在于根據(jù)實際情況調(diào)整對參考單元功率的計算方法[6?9]。

1.1 傳統(tǒng)的臨近單元平均恒虛警

傳統(tǒng)的臨近單元平均恒虛警處理方法是通過對參考單元的功率電平求和取平均來估計待測單元的雜波功率電平，算法如圖2所示。

待測單元相鄰兩個保護電平用來防止目標能量泄露到參考單元中影響檢測效果。為待測單元雜波功率電平的估計，通過對個參考單元的處理得到。為門限系數(shù)，它和的乘積作為參考門限電平。當待測單元的電平值超過門限值，認為有目標；反之，認為無目標。

可以構(gòu)建單元平均恒虛警處理的檢測門限為：

式中為與虛警概率有關(guān)的門限乘積系數(shù)，其表達式如下：

式中：表示虛警概率；表示所有參考單元的數(shù)目。采用單元平均恒虛警處理方法需要滿足兩個條件：

① 臨近參考單元所含雜波的統(tǒng)計特性與待檢單元的一致；

② 參考單元中不包含任何目標，僅存在干擾噪聲[10]。

1.2 修正單元平均恒虛警

連發(fā)火炮是復(fù)雜的多目標環(huán)境，由于干擾和其他目標的存在，檢測性能有所損失，傳統(tǒng)的臨近單元平均恒虛警處理方法不能滿足信號檢測的要求，因而本文采用一種修正的臨近單元平均恒虛警處理方法。為了滿足條件①，保持雜波干擾在距離維的均勻性，參考距離單元的數(shù)量不宜太多，本文中前導(dǎo)參考窗和滯后參考窗都包含4個距離單元；為了滿足條件②，并結(jié)合實驗需要，本文首先對前導(dǎo)和滯后參考單元的功率電平值進行平均，剔除電平值超出平均值倍的距離單元（假定序號為值視檢測效果調(diào)整）的影響，因此，修正的單元平均恒虛警處理表達式為：

式中為剔除包含強目標參考單元后的參考單元數(shù)量。

2 雙門限檢測

高射頻彈丸多目標進行恒虛警處理主要有兩個目的：消除噪聲帶來的虛警；消除雜波以及由于偽碼調(diào)相回波信號相關(guān)旁瓣形成的虛警。

由于在一定積累時間內(nèi)，目標相對于雷達的徑向速度會發(fā)生變化，使得目標在距離?多普勒譜中的表現(xiàn)形式不是一個孤立的點，而是在距離維和多普勒維上都有一定程度的擴展，表現(xiàn)為在RD譜中占據(jù)多個單元格的面目標。目標一維CFAR檢測算法只在檢測單元所在的距離向或多普勒向選取參考單元，不能充分地利用目標在距離?多普勒譜上的信息，而目標二維CFAR檢測能更好地對噪聲背景進行估計[11]。因此，針對彈丸目標的運動性，在進行目標CFAR檢測時，本文采用目標二維CFAR檢測算法。將一維時域回波數(shù)據(jù)轉(zhuǎn)化為二維的時頻數(shù)據(jù)矩陣，其數(shù)值大小對應(yīng)信號幅度，如圖3所示。橫軸對應(yīng)多普勒即速度單元，縱軸對應(yīng)距離即時間單元。

在二維區(qū)域上進行CFAR檢測，其檢測假設(shè)為：

式中：和分別表示多普勒單元和距離單元；表示對應(yīng)的目標在多普勒維和距離維的回波信號；為背景噪聲或雜波；和是對應(yīng)距離?多普勒頻率圖的幅度。

將二維數(shù)據(jù)矩陣映射到檢測空間用以進行后續(xù)的恒虛警處理以及目標檢測。在修正單元平均恒虛警的基礎(chǔ)上，采用檢測雙門限。門限一為低門限，用來消除噪聲功率帶來的虛警變化，算法基本流程如圖4所示。

由于每個距離門的數(shù)據(jù)是一個數(shù)據(jù)矩陣（行列），橫軸對應(yīng)頻率點數(shù)縱軸對應(yīng)時間點數(shù)對每一列的所有數(shù)據(jù)進行求和，然后求平均值。為了加強檢測的魯棒性，本文引入動態(tài)閾值得到該列數(shù)據(jù)的門限，表示如下：

門限計算完畢后，將該列中每一個數(shù)據(jù)與門限進行比較，大于等于該門限的數(shù)據(jù)保留下來，如果小于門限則將該數(shù)據(jù)剔除，用0進行代替。該列比較完畢后，對數(shù)據(jù)矩陣下一個頻率點進行同樣處理，遍歷所有頻率點。

由于信道的不一致性以及目標特性的復(fù)雜性，造成目標之間干擾嚴重，采用單一門限處理后的結(jié)果有時不太理想，因此本文增加了恒虛警門限二。恒虛警門限二的基礎(chǔ)是一組彈丸回波信號幅度的抖動保持在一定的范圍內(nèi)。在門限一處理的基礎(chǔ)上，根據(jù)目標的動態(tài)變化范圍首先找到該距離門內(nèi)各個單元數(shù)據(jù)的最大值，然后用最大值與門限二增量的差值作為門限2。高于門限的數(shù)據(jù)不變，低于門限的數(shù)據(jù)清零。算法基本流程如圖5所示。

門限二的計算公式為：

3 實驗結(jié)果與分析

本實驗采用實測3連發(fā)步槍信號模擬連發(fā)火炮信號，實驗結(jié)果如圖6～圖8所示。

圖6是某距離門短時傅里葉變換后得到的時間速度圖。明暗代表能量的大小。3條明亮的斜線代表3發(fā)彈經(jīng)過該距離門時的回波信號。具有較強的噪聲，彈跡表現(xiàn)不是一個孤立的點，而是在距離維和多普勒維上都有一定程度的擴展，在時頻圖中形成面目標，導(dǎo)致無法準確提取此時彈丸的速度和時間信息。明亮的豎線是由偽隨機編碼的碼重復(fù)周期帶來的干擾，能量超過彈丸回波，嚴重影響后續(xù)的數(shù)據(jù)處理。

圖7是按照修正的單元平均恒虛警處理后的結(jié)果，從圖7中可以看出背景噪聲大部分被濾除，可以清晰地看到3條明亮的斜線，但仍然不是孤立的點跡，一部分噪聲不能濾除，很多亮點散布。值得注意的是，明亮的豎線由于能量恒定，通過臨近單元能量門限的劃定可以將其濾除。

圖8是檢測雙門限的處理結(jié)果，此時門限一修正量dB，門限二修正量與圖7相比，線目標已經(jīng)基本變成點目標，可以較為準確地提取出3發(fā)彈丸在此距離門內(nèi)的時間和速度信息。結(jié)合目標檢測算法，最終提取出此時3發(fā)彈丸的運動參數(shù)如表1所示。

4 結(jié) 論

為了充分利用彈丸目標的運動信息，本文在距離?多普勒二維進行目標檢測處理。在傳統(tǒng)單元平均恒虛警處理的基礎(chǔ)上，針對連發(fā)火炮多目標環(huán)境提出修正的單元平均恒虛警處理技術(shù)，剔除不滿足功率條件的參考單元后進行檢測。以恒虛警處理為基礎(chǔ)，結(jié)合雙門限實現(xiàn)多目標的檢測，為下一步的目標參數(shù)提取提供基礎(chǔ)。實際工程應(yīng)用中，剔除單元能量的選定、門限動態(tài)閾值的設(shè)定可以根據(jù)實際情況做出相應(yīng)調(diào)整。

參考文獻

[1] 黃坤，劉詠，宋石玉.恒虛警處理中基于噪聲采樣的檢測門限設(shè)定[J].電訊技術(shù)，2010，50（8）：103?105.

[2] 秦玉華，郝程鵬.基于自動刪除平均和單元平均的恒虛警檢測器[J].探測與控制學(xué)報，2011，33（1）：15?17.

[3] 陳建軍，黃孟俊，邱偉.海雜波下的雙門限恒虛警目標檢測新方法[J].電子學(xué)報，2011，39（9）：2136?2140.

[4] 程曉莉，王曉暉，梁車平.幾種恒虛警處理方法及性能比較[J].電子元器件應(yīng)用，2010，12（3）：59?60.

[5] 喬宏章，張軍.泄漏電纜周界監(jiān)視技術(shù)研究[J].無線電工程，2013，43（3）：43?46.

[6] 許江湖，張明敏，王平波.基于自動刪除算法的最小選擇恒虛警檢測器[J].武漢理工大學(xué)學(xué)報（交通科學(xué)與工程版），2006，30（6）：1065?1068.

[7] 簡濤，蘇峰，何友，等.復(fù)合高斯雜波下距離擴展目標的自適應(yīng)檢測[J].電子學(xué)報，2012，40（5）：991?994.

[8] 杜鵬飛，張祥軍.單元平均恒虛警率檢測中的一個新結(jié)論[J].現(xiàn)代雷達，2007，29（2）：60?62.

[9] 桂任舟.利用二維恒虛警進行非均勻噪聲背景下的目標檢測[J].武漢大學(xué)學(xué)報，2012，37（3）：355?357.

[10] RICHARDS M A.雷達信號處理基礎(chǔ)[M].北京：電子工業(yè)出版社，2012.

[11] 梁建.高頻地波雷達目標二維CFAR檢測及軟件實現(xiàn)[D].青島：中國海洋大學(xué)，2014.