林 龍,吳 碧

(中國(guó)人民解放軍91388部隊(duì),廣東湛江 524022)

電子偵察是現(xiàn)代電子戰(zhàn)爭(zhēng)的關(guān)鍵部分,而雷達(dá)的工作模式識(shí)別是電子偵察的重要功能。通過雷達(dá)偵察技術(shù)掌握敵方雷達(dá)的信號(hào)參數(shù),依據(jù)雷達(dá)波形的設(shè)計(jì)原理對(duì)敵方雷達(dá)的功能加以推斷,判斷當(dāng)前敵方雷達(dá)的工作模式和威脅程度,為對(duì)敵技術(shù)打擊和硬摧毀提供重要技術(shù)支持,在實(shí)際戰(zhàn)場(chǎng)環(huán)境中起到關(guān)鍵性的作用[1]。

在不同工作模式下,不同平臺(tái)的雷達(dá)所表現(xiàn)出的信號(hào)特征具有多樣性,而且不同型號(hào)雷達(dá)之間存在交集。目前使用的雷達(dá)工作模式識(shí)別技術(shù)大部分存在識(shí)別效率低,更新模型工作量大等不足。傳統(tǒng)基于雷達(dá)數(shù)據(jù)庫(kù)的雷達(dá)工作模式識(shí)別方法對(duì)數(shù)據(jù)庫(kù)的要求較高,它需要大量可靠的經(jīng)驗(yàn)知識(shí),且識(shí)別所需的時(shí)間較長(zhǎng)[2]。文獻(xiàn)[3]通過分析影響相控陣?yán)走_(dá)工作模式的不同特征參數(shù), 選擇合適的相控陣?yán)走_(dá)脈沖描述字中的特征參數(shù)或特征參數(shù)組, 通過模糊聚類方法,來實(shí)現(xiàn)工作模式的識(shí)別,但該方法運(yùn)算量大,收斂慢,實(shí)戰(zhàn)中不利于快速判斷。文獻(xiàn)[4]通過相參脈沖間隔(Coherent Pulse Interval,CPI)特征矩陣,并對(duì)其進(jìn)行關(guān)聯(lián)來識(shí)別雷達(dá)工作模式,但只對(duì)CPI特征矩陣關(guān)聯(lián)進(jìn)行識(shí)別存在一定模糊性, 如對(duì)速度搜索(Velocity Search,VS)搜索模式使用的高重頻(High Pulse Repetition Frequency,HPRF)波形與單目標(biāo)跟蹤模式下使用的單組脈沖重復(fù)間隔波形就難以識(shí)別出工作模式。文獻(xiàn)[5]通過判別敵方機(jī)載雷達(dá)的掃描包絡(luò)進(jìn)行工作模式識(shí)別,但它僅能對(duì)工作狀態(tài)進(jìn)行判斷,無法識(shí)別搜索跟蹤波形。文獻(xiàn)[6]僅能實(shí)現(xiàn)機(jī)載相控陣?yán)走_(dá)工作模式的功能級(jí)仿真識(shí)別,并沒有對(duì)信號(hào)級(jí)仿真進(jìn)行研究。文獻(xiàn)[7] 只能對(duì)四種典型的雷達(dá)空空工作模式進(jìn)行識(shí)別,對(duì)更復(fù)雜的工作模式的判斷能力不足。文獻(xiàn)[4-7]都僅適用于機(jī)載雷達(dá),對(duì)其他平臺(tái)雷達(dá)未做考慮。

本文首先分析和總結(jié)了雷達(dá)搜索和跟蹤的特點(diǎn)。結(jié)合截獲信號(hào)的幅度分布特點(diǎn),利用幅度分布函數(shù)、系統(tǒng)聚類[8-10]、K-S(Kolmogorov-Smirnov)檢驗(yàn)[11-12]等方法對(duì)雷達(dá)的搜索與跟蹤波形進(jìn)行識(shí)別,實(shí)驗(yàn)證明了本方法的有效性。該方法適用于多雷達(dá)平臺(tái),映射匹配更為快捷準(zhǔn)確。

1 信號(hào)分類方法

典型的雷達(dá)偵察系統(tǒng)所獲取的是直接來自于雷達(dá)的發(fā)射脈沖信號(hào)的極化及調(diào)制信息,這些參數(shù)一般組合到一起,成為脈沖描述字(PDW)。由于環(huán)境的復(fù)雜性,雷達(dá)偵察系統(tǒng)所獲取的信息通常為是較長(zhǎng)的數(shù)據(jù)串,同時(shí)數(shù)據(jù)“輻射源”可能較多,信號(hào)復(fù)雜。為有效提取所需信號(hào),這里先采用系統(tǒng)聚類法對(duì)數(shù)據(jù)進(jìn)行初步的分類。

作為聚類分析中應(yīng)用最為廣泛的一種方法,系統(tǒng)聚類法的基本原理可簡(jiǎn)單理解為:首先，將一定數(shù)量的樣本各自看成一類,然后，根據(jù)樣本的相似程度,將相似程度最高的兩類進(jìn)行合并,然后，考慮合并后的類與其他類之間的相似程度,再進(jìn)行合并。重復(fù)這一過程,直至將所有的樣本合并為一類[13]。這里使用的是歐氏距離。其定義為:假設(shè)有兩個(gè)n維樣本x1=(x11,x12,…x1n)和x2=(x21,x22,…x2n),則它們的歐氏距離為

(1)

結(jié)合偵察場(chǎng)景數(shù)據(jù),聚類分析具體的操作步驟如下:

1)選擇某一工作模式下的波束掃描時(shí)間內(nèi)變化較小的參數(shù)作為聚類分析的特征量(如載頻、脈寬、帶寬等)。

2)根據(jù)每個(gè)樣本的特征量所組成的歐氏距離,選擇合理的門限,對(duì)場(chǎng)景數(shù)據(jù)進(jìn)行初步聚類。

3)初步聚類完成后,把聚類后的信號(hào)根據(jù)脈沖到達(dá)時(shí)間(TOA)的間隔大小再次進(jìn)行分割,分割后的信號(hào)作為單個(gè)信號(hào)組。

4)把分割好的信號(hào)組按照幅度變化情況進(jìn)行合并，將在連續(xù)時(shí)間內(nèi)的幅度變化不大的信號(hào)組合并到一起。在發(fā)射信號(hào)時(shí),由于抗干擾或者解模糊的需要,會(huì)發(fā)射組間頻率捷變或者組間重頻參差信號(hào),把這些相關(guān)的信號(hào)組歸納到一起有利于觀察信號(hào)的幅度、脈沖重復(fù)周期(Pluse Repetition Time,PRT)及載頻的變化情況,能夠更好地確定雷達(dá)信號(hào)對(duì)應(yīng)的雷達(dá)功能。

從原始數(shù)據(jù)中是無法直接獲取PRT這一參數(shù)的,但原始數(shù)據(jù)中會(huì)有TOA參數(shù)。只要把信號(hào)合理分組,在每一組信號(hào)內(nèi)的相鄰信號(hào)的TOA之差即可作為一個(gè)PRT,而且,可用PRT驗(yàn)證信號(hào)組合前的信號(hào)分類是否正確。

2 雷達(dá)波形識(shí)別

搜索和跟蹤是雷達(dá)基本功能,針對(duì)不同的應(yīng)用場(chǎng)景、不同的體制雷達(dá),搜索與跟蹤所采用的雷達(dá)波形參數(shù)可能范圍存在較大的重疊區(qū)域,對(duì)于多功能雷達(dá)尤為明顯,因此，通常上述功能僅從雷達(dá)參數(shù)上是無法識(shí)別的。本文將結(jié)合雷達(dá)的資源調(diào)度特點(diǎn)完成所偵察對(duì)象的搜索與跟蹤功能識(shí)別分類。

2.1 雷達(dá)基本波形分類

雷達(dá)基本波形根據(jù)重復(fù)頻率可以分為低重頻(Low Pulse Repetition Frequency,LPRF)波形、中重頻(Medium Pulse Repetition Frequency,MPRF)波形和HPRF波形;根據(jù)功能模式可以分為中重頻邊搜索邊測(cè)距(Medium Range While Search,MRWS)模式和高重頻邊搜索邊測(cè)距(High Range While Search,HRWS)模式。常見雷達(dá)的基本波形[12]如下:預(yù)警LPRF、預(yù)警MRWS、預(yù)警HRWS、預(yù)警VS、火控搜索LPRF、火控搜索MRWS(參差)、火控搜索HRWS(線性調(diào)頻)、火控跟蹤MPRF、火控跟蹤HPRF或VS、火控搜索HPRF(參差)[14]。

區(qū)分基本波形的依據(jù)是其對(duì)應(yīng)的平臺(tái)及工作模式。一般情況下,基本波形都有對(duì)應(yīng)的參數(shù)設(shè)置范圍。從截獲的雷達(dá)信息可以一定程度上判斷出其所屬的基本波形種類,如機(jī)載火控HPRF波形與機(jī)載多功能速度搜索波形的參數(shù)值就很接近,故很難從截獲的波形參數(shù)值上判斷波形類別,特別是對(duì)跟蹤及搜索狀態(tài)的判斷。

2.2 雷達(dá)波形幅度特性

作為實(shí)際應(yīng)用中常用的兩種工作模式,搜索與跟蹤在波形參數(shù)設(shè)計(jì)時(shí)通常會(huì)使用有一定區(qū)別的系統(tǒng)參數(shù)。從資源調(diào)度的角度出發(fā),兩者還有一個(gè)重要區(qū)別為:數(shù)據(jù)率的不同[4]。通常情況下，數(shù)據(jù)率指單位時(shí)間內(nèi)回波對(duì)某波位的重訪次數(shù)或者某波位的重訪時(shí)間間隔。為了方便推導(dǎo),本文對(duì)數(shù)據(jù)率進(jìn)行擴(kuò)展定義。下面首先給出數(shù)據(jù)率的廣義定義:一定時(shí)間頻率內(nèi)數(shù)據(jù)出現(xiàn)的頻率,這里稱之為廣義數(shù)據(jù)率,如不加以說明,后文提到的均為廣義數(shù)據(jù)率[15]。

從雷達(dá)偵察的角度對(duì)搜索、跟蹤功能下的截取回波特點(diǎn)進(jìn)行分析,主要特點(diǎn)可總結(jié)為:

1)搜索:由于目標(biāo)位置的未知,雷達(dá)波束掃描范圍較大,偵察雷達(dá)所截獲的波形幅度變化范圍較大。同時(shí)其幅值分布滿足正態(tài)分布,通常采用單樣本K-S檢驗(yàn)來檢測(cè)。

2)跟蹤:特定目標(biāo)位置已知,雷達(dá)波束掃描集中于小范圍或者特定目標(biāo)的鎖定,而非大范圍的掃描,因此偵察雷達(dá)所截獲的雷達(dá)波形幅值會(huì)集中于小范圍內(nèi),甚至某幾個(gè)點(diǎn)的幅值(點(diǎn)的個(gè)數(shù)與跟蹤目標(biāo)的個(gè)數(shù)相關(guān))。

2.3 K-S檢驗(yàn)

單樣本的K-S檢驗(yàn)是用來檢驗(yàn)一個(gè)數(shù)據(jù)的觀測(cè)經(jīng)驗(yàn)分布是否是已知的理論分布[16]。檢驗(yàn)方法為:以樣本數(shù)據(jù)的累計(jì)頻數(shù)分布與特定理論分布比較,若兩者間的差距很小,則推論該樣本取自某特定分布族[17]。

首先,假設(shè)檢驗(yàn)問題為:H0:樣本所來自的總體分布服從某特定一維連續(xù)分布F;H1:樣本所來自的總體分布不服從某特定一維連續(xù)分布F,F為理論分布的分布函數(shù)[18]。

構(gòu)造檢驗(yàn)統(tǒng)計(jì)量

(2)

H0為真,則Z分布收斂于Kolmogorov分布,即當(dāng)樣本取自一維連續(xù)分布F時(shí),

(3)

注:當(dāng)F為連續(xù)分布時(shí),隨機(jī)變量K的分布不依賴于F。

2.4 雷達(dá)波形識(shí)別流程

在判斷雷達(dá)搜索跟蹤波形時(shí),先計(jì)算出每個(gè)幅度的平均頻率,然后判斷場(chǎng)景中包含的基本波形是否只有一種。若該場(chǎng)景中只有一種基本波形,波形識(shí)別流程如圖1所示;若該場(chǎng)景中包含多種基本波形,波形識(shí)別流程如圖2所示。

圖1 單一基本波形識(shí)別流程

3 仿真驗(yàn)證

假設(shè)仿真場(chǎng)景中,雷達(dá)處于搜索加跟蹤狀態(tài),且由于反偵察、增大作用距離等需要,在發(fā)射的時(shí)候會(huì)采用參差重頻、頻率捷變等技術(shù)。在參數(shù)選擇上,PRT為:7、9、10、11、12、13、14、15、17(單位：us);載頻捷變采用載頻為:9.49、9.53、9.55、9.61、9.64、9.66、9.68、9.70、9.74(單位：GHz)。

圖2 多種基本波形識(shí)別流程

3.1 信號(hào)分類驗(yàn)證仿真

對(duì)信號(hào)組進(jìn)行隨機(jī)抽取,觀察分類后的信號(hào)組的幅度、PRT、載頻的變化情況,如圖3所示。從圖3可以看出,在某一組信號(hào)內(nèi)幅度基本不變,PRT與載頻都有參差。圖中仿真結(jié)果驗(yàn)證了基于聚類的信號(hào)分類方法的合理性。

3.2 基本波形識(shí)別仿真

對(duì)3.1中已分類的雷達(dá)波形做初步的信號(hào)仿真,進(jìn)行參數(shù)識(shí)別。對(duì)于分類后的每組信號(hào)我們稱為子波,根據(jù)該子波的參數(shù)特性,我們對(duì)其進(jìn)行初步劃分判斷,其依據(jù)主要是常見雷達(dá)子波功能的參數(shù)范圍,如載頻、PRT等。

圖3 分類后的若干信號(hào)組的參數(shù)變化情況

圖4 子波序號(hào)識(shí)別結(jié)果

仿真中場(chǎng)景的子波識(shí)別結(jié)果如圖4所示,子波序號(hào)在9號(hào)和10號(hào)兩個(gè)基本波形中,其中大部分子波屬于9號(hào)波形。根據(jù)前文基本波形分類,有必要對(duì)兩個(gè)波形進(jìn)行驗(yàn)證,而且對(duì)于9號(hào)波形而言,還需要對(duì)其究竟是搜索還是跟蹤狀態(tài)進(jìn)行進(jìn)一步判定。

3.3 搜索跟蹤波形識(shí)別判斷仿真

對(duì)于雷達(dá)搜索與跟蹤波形的判斷,本文主要結(jié)合波形幅度統(tǒng)計(jì)分布及相應(yīng)的概率分布進(jìn)行仿真判斷。首先對(duì)9、10號(hào)子波的幅度進(jìn)行統(tǒng)計(jì)分析如圖5、6。

圖5 9號(hào)波形幅度統(tǒng)計(jì)圖

圖6 10號(hào)波形幅度統(tǒng)計(jì)圖

對(duì)圖5和圖6進(jìn)行分析,兩種波形的幅度統(tǒng)計(jì)結(jié)果如表1所示。

表1 幅度統(tǒng)計(jì)結(jié)果

然后,根據(jù)兩種波形的幅度統(tǒng)計(jì)，計(jì)算它們的正態(tài)分布概率,如圖7、圖8所示。

圖7 9號(hào)波形正態(tài)分布概率圖

圖8 10號(hào)波形正態(tài)分布概率圖

從圖7可以看出,9號(hào)波形的幅度并不符合正態(tài)分布,為跟蹤波形;從圖8可以看出,10號(hào)波形的幅度更符合正態(tài)分布,為搜索波形。

根據(jù)上述判斷結(jié)果,通過幅度統(tǒng)計(jì)分布及概率分布雙重標(biāo)準(zhǔn)的方式,可推導(dǎo)出9號(hào)波形處于單目標(biāo)跟蹤狀態(tài),10號(hào)波形處于搜索狀態(tài),即雷達(dá)處于搜索加跟蹤的工作狀態(tài),識(shí)別結(jié)果正確。

4 結(jié)束語

本文根據(jù)雷達(dá)信號(hào)的功能特點(diǎn),首先利用了系統(tǒng)聚類法等數(shù)學(xué)工具,結(jié)合實(shí)際情況,如時(shí)序的連續(xù)性、參數(shù)的變化情況等因素對(duì)雷達(dá)信號(hào)進(jìn)行了有效分類,特別是通過分析雷達(dá)在不同功能下的波形幅度分布情況,結(jié)合K-S檢驗(yàn)方法對(duì)雷達(dá)信號(hào)分類后的子波進(jìn)行功能識(shí)別判斷,給出了雷達(dá)工作模式識(shí)別的詳細(xì)流程。仿真實(shí)驗(yàn)證明了該方法的正確性,在雷達(dá)偵察領(lǐng)域具備良好的應(yīng)用價(jià)值。