反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)設(shè)計(jì)*

謝奇峰

（空降兵訓(xùn)練基地，廣西 桂林 541003）

反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)設(shè)計(jì)*

謝奇峰

（空降兵訓(xùn)練基地，廣西 桂林 541003）

為了提升某型炮瞄雷達(dá)識(shí)別反輻射導(dǎo)彈的能力，設(shè)計(jì)了一種面向炮瞄雷達(dá)的反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)。系統(tǒng)由雷達(dá)信號(hào)采集終端和信號(hào)處理平臺(tái)兩部分構(gòu)成。信號(hào)采集終端在雷達(dá)同步脈沖的同步下，負(fù)責(zé)采集I路視頻信號(hào)并傳送至信號(hào)處理平臺(tái)，信號(hào)處理平臺(tái)在頻域?qū)崟r(shí)處理I路視頻信號(hào)，判斷是否存在反輻射導(dǎo)彈，并發(fā)出報(bào)警提示。經(jīng)試驗(yàn)驗(yàn)證，該系統(tǒng)可以識(shí)別出高速反輻射導(dǎo)彈。

炮瞄雷達(dá)，反輻射導(dǎo)彈，信號(hào)處理

Abstract：To promoting certain Gun-Pointing radar’s ability to discover ARM，a ARM’s distinguish and alarm system for Gun-Pointing radar has been designed.The system is composed by data acquisition terminal and data processing platform，the former is responsible of collection of video signal and upload with radar’s pulse，the latter is responsible of real-time processing video signal in frequency domain.When find ARM，alarm device sounds the alarm immediately.The system can identify a high-speed ARM target in some trials.

Key words：gun-pointing radar，ARM，information retrieval

0 引言

反輻射導(dǎo)彈在現(xiàn)代戰(zhàn)爭(zhēng)中的作用愈顯重要，嚴(yán)重威脅到地面雷達(dá)的生存，雷達(dá)是否具有抗擊反輻射導(dǎo)彈能力是現(xiàn)代雷達(dá)的重要指標(biāo)之一。然而，一些種類(lèi)炮瞄雷達(dá)的I、II波段接收機(jī)為對(duì)數(shù)接收機(jī)，在跟蹤敵機(jī)的情況下只能監(jiān)視反輻射導(dǎo)彈與載機(jī)的分離過(guò)程，無(wú)法對(duì)發(fā)射后的反輻射導(dǎo)彈進(jìn)行識(shí)別并提供告警提示，地面操縱人員不能有效跟蹤反輻射導(dǎo)彈，這給炮瞄雷達(dá)抗擊反輻射導(dǎo)彈行動(dòng)帶來(lái)很大不利。目前，國(guó)內(nèi)外相關(guān)領(lǐng)域的專(zhuān)家和學(xué)者從新技術(shù)和新方法等方面研究了抗擊反輻射導(dǎo)彈的方法與途徑，有些技術(shù)和方法在預(yù)警雷達(dá)上得到了應(yīng)用。文獻(xiàn)［1-3］研究了采用改進(jìn)的方法，提高預(yù)警雷達(dá)識(shí)別反輻射導(dǎo)彈的能力。文獻(xiàn)［4-5］研究了抗擊反輻射導(dǎo)彈的戰(zhàn)術(shù)措施。關(guān)于炮瞄雷達(dá)識(shí)別發(fā)射后的反輻射導(dǎo)彈的研究還未見(jiàn)諸報(bào)道。針對(duì)炮瞄雷達(dá)在抗擊反輻射導(dǎo)彈性能上的不足，本文緊密結(jié)合炮瞄雷達(dá)特點(diǎn)，提出了在炮瞄雷達(dá)上增加反輻射導(dǎo)彈識(shí)別裝置的解決思路，并構(gòu)建了反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)，為下一步在炮瞄雷達(dá)進(jìn)行適應(yīng)性改造打下基礎(chǔ)。

1 系統(tǒng)的設(shè)計(jì)思路

1.1 系統(tǒng)的原理方案

根據(jù)課題前期研究結(jié)論，反輻射導(dǎo)彈的回波為一串振幅受多普勒頻率調(diào)制的脈沖。輸出包絡(luò)為一正弦信號(hào)，其振幅為為多普勒角頻率，Tr為雷達(dá)脈沖周期，回波是多普勒頻率fd的函數(shù)。要提高炮瞄雷達(dá)對(duì)反輻射導(dǎo)彈的識(shí)別與告警能力，需要考慮來(lái)自I波段接收機(jī)視頻信號(hào)的采集和信號(hào)的處理兩個(gè)問(wèn)題。本文使用DSP設(shè)計(jì)了一個(gè)高速采集電路，采集電路以雷達(dá)上零位脈沖作為采集的外觸發(fā)信號(hào)，在零位脈沖的觸發(fā)下以設(shè)定的速率實(shí)時(shí)采集視頻信號(hào)，并把采集到的數(shù)據(jù)送往信號(hào)處理平臺(tái)。信號(hào)處理平臺(tái)在頻域?qū)崟r(shí)處理采集到的數(shù)據(jù)，分析是否存在反輻射導(dǎo)彈，若檢測(cè)到反輻射導(dǎo)彈立刻向報(bào)警器發(fā)出指令，報(bào)警器工作并發(fā)出警報(bào)聲。反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)方案如圖1所示。

圖1 反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)方案

1.2 雷達(dá)信號(hào)采集思路

為確保炮瞄雷達(dá)的探測(cè)距離，雷達(dá)上設(shè)定的脈沖周期Tr與最大探測(cè)距離Dmax相關(guān)聯(lián)。同時(shí)為消除盲速影響，還設(shè)置了3種脈沖重復(fù)周期，脈沖間隔非均等，故雷達(dá)信號(hào)的采集不能使用一次性方式，必須采用分組的方式，即一個(gè)脈沖周期內(nèi)采集一組數(shù)據(jù)。本文把零位脈沖作為各組采集的觸發(fā)信號(hào)。根據(jù)香農(nóng)采樣定理，綜合考慮脈沖周期、脈沖寬度、目標(biāo)距離、效率等因素，在保證探測(cè)距離的條件下，盡量減少一個(gè)脈沖周期內(nèi)的采集量。各組信號(hào)采集的參數(shù)如式（1）所示。

式中，fs為采樣頻率，N為一個(gè)脈沖寬度被采樣的次數(shù)，為炮瞄雷達(dá)的脈沖寬度，Ts為一個(gè)脈沖周期內(nèi)采集持續(xù)時(shí)間，Nts為一個(gè)脈沖周期內(nèi)的采樣數(shù)量。為了保證精度，N一般不小于2。信號(hào)采集方法如圖2所示。

1.3 信號(hào)檢測(cè)思路

圖2 信號(hào)采集思路

根據(jù)反輻射導(dǎo)彈回波包絡(luò)為正弦信號(hào)的研究結(jié)論，考慮到頻域信號(hào)分析具有明確的含義，獲取信號(hào)特征及其參數(shù)不受信號(hào)周期限制，信號(hào)特征明顯且魯棒性好等特點(diǎn)，反輻射導(dǎo)彈的信號(hào)檢測(cè)在頻域進(jìn)行。利用頻譜分析手段判斷回波的脈沖包絡(luò)是否為正弦形式，若是則認(rèn)為存在反輻射導(dǎo)彈，并進(jìn)一步依據(jù)頻率估測(cè)目標(biāo)的速度。根據(jù)頻域數(shù)據(jù)處理要求，傅里葉變換的采樣點(diǎn)數(shù)N一般取2的整數(shù)次冪。綜合考慮效率和準(zhǔn)確度，本文設(shè)計(jì)了24、25、26點(diǎn)數(shù)據(jù)的離散傅里葉變換算法。根據(jù)傅里葉變換結(jié)果，判斷是否包含反輻射導(dǎo)彈，判斷方法如式（2）所示。

式中，X（k）為N點(diǎn)離散傅里葉變換值，kmax為傅里葉變換的最大值所對(duì)應(yīng)的k值。

由于數(shù)據(jù)采集的分段性，進(jìn)行傅里葉變換之前需要對(duì)采集到的N組數(shù)據(jù)重新組合，形成Nts（一個(gè)脈沖周期內(nèi)的采樣數(shù)量）組具有N個(gè)采集點(diǎn)的數(shù)據(jù)組合。具有N個(gè)采集點(diǎn)的數(shù)據(jù)組合作為傅里葉變換的輸入數(shù)據(jù)。數(shù)據(jù)組合方式由式（3）中的縱列確定，如第 1 組數(shù)據(jù)為（x11、x21、…、xN1）。

信號(hào)的檢測(cè)從第1組數(shù)據(jù)的傅里葉變換開(kāi)始，直到檢測(cè)出反輻射導(dǎo)彈信號(hào)，所以最多需要進(jìn)行Nts次傅里葉變換。發(fā)現(xiàn)反輻射導(dǎo)彈信號(hào)后，還需要計(jì)算出導(dǎo)彈的距離。導(dǎo)彈距離D由式（4）確定。

式中，n為檢測(cè)到反輻射導(dǎo)彈信號(hào)的數(shù)據(jù)組合序號(hào)。為了提高檢測(cè)的可靠性，數(shù)據(jù)采集模塊在一個(gè)脈沖寬度內(nèi)必須完成2次以上采集。在第n組數(shù)據(jù)中檢測(cè)到反輻射導(dǎo)彈信號(hào)后，再次對(duì)第n+1組數(shù)據(jù)進(jìn)行傅里葉變換以確認(rèn)檢測(cè)結(jié)果。

2 具體設(shè)計(jì)方案

根據(jù)式（1）得到的fs，設(shè)計(jì)了采集速度最高為40 MSPS的采集電路，以保證對(duì)每個(gè)脈沖至少能采集2次以上。采集電路對(duì)Ⅰ波段接收機(jī)的信號(hào)進(jìn)行采集，采集到規(guī)定的數(shù)據(jù)后將之送往信號(hào)處理平臺(tái)。信號(hào)處理平臺(tái)對(duì)信號(hào)進(jìn)行分析與處理，若發(fā)現(xiàn)反輻射導(dǎo)彈，發(fā)送報(bào)警指令，報(bào)警裝置接收到指令后發(fā)出警報(bào)。

2.1 信號(hào)采集終端設(shè)計(jì)

信號(hào)采集終端實(shí)際上是一個(gè)信號(hào)采集和聲光報(bào)警集成裝置，進(jìn)行信號(hào)的采集，并根據(jù)報(bào)警指令實(shí)現(xiàn)聲光報(bào)警。信號(hào)采集終端主要包括邏輯控制器、模數(shù)轉(zhuǎn)換器、存儲(chǔ)器、數(shù)據(jù)傳輸通道、報(bào)警器等單元。信號(hào)采集終端連接示意如圖3所示。A/D轉(zhuǎn)換電路選用TLC5540高速模數(shù)轉(zhuǎn)換芯片，采用說(shuō)明書(shū)上的典型電路方案［6-7］，在炮瞄雷達(dá)零位脈沖的作用下，芯片啟動(dòng)持續(xù)Ts時(shí)長(zhǎng)的采樣。

圖3 信號(hào)采集部分硬件連接示意圖

2.2 信號(hào)處理平臺(tái)設(shè)計(jì)

信號(hào)處理平臺(tái)具備信號(hào)采集的控制和信號(hào)處理等功能。信號(hào)處理平臺(tái)設(shè)計(jì)的重點(diǎn)是功能函數(shù)的實(shí)現(xiàn)。功能函數(shù)主要包括定時(shí)響應(yīng)函數(shù)OnTimer、波形繪制函數(shù)PaintWave、頻譜繪制函數(shù)PaintPinPu、傅里葉變換函數(shù)FFT、延遲函數(shù)WSDelay和報(bào)警函數(shù)BaoJing。為提高代碼移植性，屏蔽復(fù)雜的協(xié)議細(xì)節(jié)，把底層交互代碼封裝成dll和lib文件，通過(guò)編程接口訪(fǎng)問(wèn)信號(hào)采集部分的資源?？紤]到交互需求，信號(hào)處理程序采用對(duì)話(huà)框類(lèi)型。主交互界面在資源模板中完成設(shè)計(jì)，如圖4所示。

圖4 信號(hào)處理平臺(tái)主界面

2.3 相關(guān)算法的設(shè)計(jì)

2.3.1 信號(hào)采集算法

當(dāng)信號(hào)處理平臺(tái)發(fā)出“打開(kāi)設(shè)備”命令后，A/D轉(zhuǎn)換電路進(jìn)入工作狀態(tài)，檢測(cè)是否存在零位脈沖。當(dāng)檢測(cè)到零位脈沖后，開(kāi)始以fs采樣頻率采集I路信號(hào)，采集Nts個(gè)采樣點(diǎn)后暫停采集，等待下一個(gè)零位脈沖。當(dāng)下一個(gè)零位脈沖產(chǎn)生之后，繼續(xù)以fs采樣頻率采集Nts個(gè)采樣點(diǎn)，等待下一個(gè)零位脈沖。當(dāng)A/D轉(zhuǎn)換電路采集到32組數(shù)據(jù)后，一次性把這些數(shù)據(jù)上傳至信號(hào)處理平臺(tái)。數(shù)據(jù)采集過(guò)程如圖5所示。

圖5 數(shù)據(jù)采集流程

數(shù)據(jù)的上傳有兩種方式，一種是每采集完一組數(shù)據(jù)就上傳給信號(hào)處理平臺(tái)，另一種是采集完32組數(shù)據(jù)后一次性上傳。第1種方式的實(shí)時(shí)性比較好，可以不間斷地對(duì)數(shù)據(jù)進(jìn)行采集、處理和顯示，但對(duì)數(shù)據(jù)的傳輸速率和信號(hào)處理平臺(tái)的處理速度有很高的要求，要求傳輸速率至少為采集速率的3倍以上，信號(hào)處理平臺(tái)進(jìn)行一次FFT變換的時(shí)間需小于520 μs。第2種方案的穩(wěn)定性比較好，無(wú)數(shù)據(jù)丟失，對(duì)硬件要求相對(duì)較低，但存在短暫的延遲。經(jīng)計(jì)算，第2種方式中進(jìn)行一次FFT變換的時(shí)間可放寬至520×32 μs，延遲時(shí)間在0.017 s之內(nèi)，3倍音速換算成距離約為17 m，能滿(mǎn)足要求。本文采用第2種數(shù)據(jù)傳輸方式，以降低對(duì)數(shù)據(jù)傳輸速率和數(shù)據(jù)處理速度的要求，提高系統(tǒng)穩(wěn)定性。N取32點(diǎn)基于以下兩點(diǎn)考慮，一是減少一次FFT變換的計(jì)算量；二是提高檢測(cè)的靈敏度，N點(diǎn)從32到258，距離延遲將從17 m增加到136 m。

2.3.2 目標(biāo)檢測(cè)算法

反輻射導(dǎo)彈檢測(cè)算法采取時(shí)域抽取基2FFT算法。設(shè)采樣后時(shí)域中的離散信號(hào)為x（n），其中n=0、1、…、N-1。將序列x（n）按n為奇、偶數(shù)分為x1（n）、x2（n）兩組序列，用兩個(gè)N/2點(diǎn)傅里葉變換來(lái)完成一個(gè)N點(diǎn)傅里葉變換的計(jì)算。設(shè)序列x（n）的長(zhǎng)度為N，且滿(mǎn)足：N=2M，M 為自然數(shù)。按 n的奇偶把 x（n）分解為兩個(gè)N/2點(diǎn)的子序列，如式（5）所示：

用N/2點(diǎn)X1（k）和X2（k）表示N點(diǎn)傅里葉變換X（k），如式（6）所示：

式中，WNk=e-j2πrk/N，為旋轉(zhuǎn)因子。由于X1（k）和X2（k）均以N/2為周期，且，X（k）又可表示為：

這樣將N點(diǎn)傅里葉變換分解為兩個(gè)N/2點(diǎn)的傅里葉變換。對(duì)每個(gè)N/2點(diǎn)傅里葉變換繼續(xù)進(jìn)行同樣的分解，直到只有最基本的2點(diǎn)傅里葉變換為止，就得到了時(shí)域抽取基2 FFT算法?；? FFT算法的計(jì)算量使乘法計(jì)算由N2次降為Nlog2N/2次，加法次數(shù)由N（N-1）次降為Nlog2N次。經(jīng)分析，該算法的第1級(jí)有1個(gè)組，往后逐級(jí)二分直到有N/2個(gè)組，每組所乘的旋轉(zhuǎn)因子是同樣的，并且每組的每個(gè)蝶形輸入數(shù)據(jù)之間的間隔是逐級(jí)減小到2。根據(jù)這個(gè)規(guī)律，基2算法的實(shí)現(xiàn)如圖6所示。

圖6 基2FFT算法流程

3 測(cè)試結(jié)果

3.1 測(cè)試條件

反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)的測(cè)試條件如下：首先疊加雷達(dá)信號(hào)噪聲和目標(biāo)回波，對(duì)疊加的模擬信號(hào)進(jìn)行采集，轉(zhuǎn)化為數(shù)字信號(hào)并上傳至信號(hào)處理平臺(tái)，觀(guān)察分析時(shí)頻域數(shù)據(jù)。以任意波信號(hào)發(fā)生器（AWG）產(chǎn)生反輻射導(dǎo)彈回波信號(hào)，把雷達(dá)上Ⅰ波段接收機(jī)送來(lái)的I路視頻信號(hào)作為噪聲源，用同相加法器疊加回波信號(hào)和噪聲這兩路信號(hào)，把疊加之后的信號(hào)作為反輻射導(dǎo)彈回波信號(hào)，見(jiàn)圖1。

采用LM358系列中的LM358AM運(yùn)算放大器構(gòu)成同相加法器。I路視頻噪聲信號(hào)和反輻射導(dǎo)彈模擬信號(hào)分別經(jīng)200 kΩ電阻后接入同相輸入端，反相輸入端經(jīng)200 kΩ電阻后接地。由于同相加法器中的4個(gè)電阻的阻值均為200 kΩ，輸出端的信號(hào)為同相輸入端的信號(hào)之和，從而實(shí)現(xiàn)信號(hào)疊加。

3.2 測(cè)試結(jié)果分析

圖7 反輻射導(dǎo)彈信號(hào)波形

圖8 反輻射導(dǎo)彈信號(hào)幅頻

圖7為采集到的反輻射導(dǎo)彈回波信號(hào)與未開(kāi)高壓時(shí)雷達(dá)信號(hào)疊加在一起的時(shí)域圖，在圖中除了存在較明顯的定周期負(fù)脈沖零位信號(hào)，其他信號(hào)雜亂無(wú)章。圖8～圖10分別為經(jīng)過(guò)FFT變換的幅頻、功率譜和對(duì)數(shù)功率譜，在這些圖中能分辨出明顯的反輻射導(dǎo)彈信號(hào)和信號(hào)頻率。對(duì)比可以看出：在時(shí)域無(wú)法判斷是否存在反輻射導(dǎo)彈；經(jīng)過(guò)32點(diǎn)FFT變換到頻域后，無(wú)反輻射導(dǎo)彈時(shí)不會(huì)形成尖峰，告警裝置不報(bào)警，存在反輻射導(dǎo)彈時(shí)出現(xiàn)了明顯的尖峰，超過(guò)了式（3）中設(shè)定的門(mén)限值，報(bào)警裝置告警，實(shí)現(xiàn)對(duì)來(lái)襲導(dǎo)彈的識(shí)別與告警。

圖9 反輻射導(dǎo)彈信號(hào)功率譜

圖10 反輻射導(dǎo)彈信號(hào)對(duì)數(shù)功率頻譜

4 結(jié)論

根據(jù)數(shù)據(jù)采集和信號(hào)處理相關(guān)理論，設(shè)計(jì)了一種反輻射導(dǎo)彈識(shí)別與告警系統(tǒng)。該系統(tǒng)具有從炮瞄雷達(dá)上采集Ⅰ波段接收機(jī)的I路視頻信號(hào)，對(duì)采集到的信號(hào)進(jìn)行實(shí)時(shí)分析與處理，檢測(cè)出反輻射導(dǎo)彈回波信號(hào)，向報(bào)警裝置發(fā)出指令，并根據(jù)信號(hào)特征計(jì)算目標(biāo)距離和速度。該裝置已經(jīng)在某型炮瞄雷達(dá)上進(jìn)行多次試驗(yàn)。從試驗(yàn)數(shù)據(jù)來(lái)看，該裝置能檢測(cè)出反輻射導(dǎo)彈信號(hào)，并通過(guò)報(bào)警裝置發(fā)出警報(bào)。下一步的研究工作是：采用集成化設(shè)計(jì)方法提高系統(tǒng)的可靠性，改進(jìn)信號(hào)處理算法，減少虛警率。

［1］劉曉明，王建東，王旭東.基于參數(shù)化時(shí)頻分析的反輻射導(dǎo)彈檢測(cè)［J］.計(jì)算機(jī)應(yīng)用，2010，30（11）：3108-3114.

［2］方前學(xué)，王永良，王首勇.高斯混合分布雜波下ARM的Wald 檢測(cè)方法［J］.現(xiàn)代雷達(dá)，2010，32（2）：61-65.

［3］陳小龍，關(guān)鍵，黃勇，等.分?jǐn)?shù)階Fourier變換在動(dòng)目標(biāo)檢測(cè)和識(shí)別中的應(yīng)用：回顧和展望［J］.信號(hào)處理，2013，29（1）：85-93.

［4］羅金亮，袁可亮.防空兵火控雷達(dá)抗反輻射導(dǎo)彈研究［J］.火控雷達(dá)技術(shù)，2011，40（2）：26-29.

［5］羅金亮.防空兵火控雷達(dá)抗反輻射導(dǎo)彈對(duì)策研究［J］.飛航導(dǎo)彈，2011，40（2）：26-29.

［6］柴曉飛.基于并行交替式數(shù)據(jù)采集系統(tǒng)的虛擬示波器研究［D］.天津：天津工業(yè)大學(xué)，2007.

［7］熊菡.手持嵌入式設(shè)備虛擬示波器系統(tǒng)設(shè)計(jì)［J］.工業(yè)儀表與自動(dòng)化裝置，2013，43（2）：76-80.

Design of ARM’s Distinguish and Alarm System

XIE Qi-feng
（Airborne Training Base，Guilin 541003，China）

TJ761

A

10.3969/j.issn.1002-0640.2017.09.037

1002-0640（2017）09-0166-05

2016-08-06

2016-09-08

空軍武器裝備資助項(xiàng)目（KJ2012277）

謝奇峰（1978- ），男，湖南湘鄉(xiāng)人，博士。研究方向：武器系統(tǒng)的控制、仿真。