針對敵我識別系統(tǒng)的個體識別研究

方棉佳，陳 朝

（1.中國人民解放軍93209部隊，北京 100085；2.中國人民解放軍95865部隊，北京 100073）

0 引 言

根據(jù)相關(guān)資料顯示，北約正全力生產(chǎn)、加裝和測試MarkⅫA 敵我識別系統(tǒng)裝備。加裝平臺涵蓋了海軍、空軍、陸軍、海岸警衛(wèi)隊的各類作戰(zhàn)和保障平臺。根據(jù)美國與日本、韓國、中國臺灣之間的相關(guān)軍事協(xié)定，可以推斷日本、韓國、中國臺灣也將具備MarkⅫA 系統(tǒng)能力，以便在局部軍事沖突中，有效獲取美國的軍事情報支援能力。

北約敵我識別系統(tǒng)的威脅越來越大，對敵我識別信號的研究也越來越多。文獻［1］研究了MarkX敵我識別信號的偵察定位技術(shù)，特別是單站定位技術(shù)。而對新型模式5敵我識別信號，普通檢測算法需要高的信噪比，采用同步脈沖相關(guān)檢測算法可在低信噪比下獲得很好的檢測效果［2］。對不同的敵我識別信號模式，需要采用不同的解調(diào)算法。對MarkX 信號，采用脈位解調(diào)技術(shù)。對S模式詢問信號，采用DPSK 解調(diào)技術(shù)［3］。對模式5 敵我識別信號，采用MSK 解調(diào)、Walsh 解擴和RS譯碼技術(shù)［4］。

針對這個問題，本文專門針對敵我識別信號進行個體識別研究。本文通過分析敵我識別信號，發(fā)現(xiàn)其區(qū)別于普通雷達信號的兩個主要特點:一、信號載頻固定不變；二、信號以碼字為單位出現(xiàn)，且脈沖時序格式固定。針對敵我識別信號的這兩個特點，本文設(shè)計的個體識別方法重點包括:一、瞬時頻率可作為有效的個體特征；二、可綜合同一串碼字里多個脈沖的個體識別結(jié)果，提高目標個體識別正確率。

1 敵我識別信號特點分析

北約敵我識別系統(tǒng)技術(shù)標準包括Mark X 和Mark XII，Mark XII 敵我識別系統(tǒng)在Mark X 基礎(chǔ)上增加了模式4 和模式5 兩種模式。Mark 系統(tǒng)技術(shù)體制與國際民航二次雷達一樣，其工作頻點是固定不變的，詢問信號載頻為1 030 MHz，應答信號載頻為1 090 MHz［11］。

Mark X詢問信號一共有6種詢問模式，分別為1，2，3/A，B，C，D 模式。各模式詢問信號均由3 個脈沖P1，P2，P3 組成，根據(jù)不同的模式P1，P3 間隔變化，模式1，2，3/A，B，C，D 對應P1，P3 間隔為3，5，8，17，21和25 μs。Mark X應答信號由12個脈位攜帶信息，對應位置有脈沖代表二進制1，無脈沖代表0。

S模式應答信號由4個前導脈沖和一串數(shù)據(jù)脈沖組成。脈沖寬度為0.5 μs，脈沖時序為0，1，3.5，4.5 μs。數(shù)據(jù)脈沖是脈位調(diào)制的。1 μs 的時間間隔對應于信息的一個bit。前0.5 μs 有脈沖，代表二進制1，后0.5 μs 有脈沖，代表二進制0。其信號時序圖如圖1所示。

圖1 S模式應答信號時序圖

Mark XII 詢問兼容Mark X，并增加了模式4 和模式5。模式4 詢問信號時序如圖2 所示。由4 個前導脈沖和一串數(shù)據(jù)脈沖組成，脈沖寬度為0.5 μs，脈沖間隔為2 μs。

圖2 模式4詢問信號時序圖

通過以上介紹可以發(fā)現(xiàn)敵我識別信號區(qū)別于雷達信號的兩個主要特點:一、信號載頻固定不變，詢問采用1 030 MHz，應答采用1 090 MHz；二、信號以脈沖串的形式出現(xiàn)，且脈沖串格式固定，在信號解調(diào)的基礎(chǔ)上，幾乎可以確定一個碼字對應的一組脈沖是來自同一個輻射源的。目前雷達多是跳頻工作，因此單純用瞬時頻率或頻譜作為個體特征并不合適，而用頻率穩(wěn)定度更可靠。目前雷達多采用參差脈沖，結(jié)合跳頻工作模式，很難確定一組脈沖是否來自用一個輻射源，因此只能針對單脈沖進行個體識別。而敵我識別信號頻率固定，頻譜上的無意調(diào)制直接反映在瞬時頻率上，因此可以將瞬時頻率作為個體特征。在信號解調(diào)的基礎(chǔ)上，可以確定一個碼字對應的一組脈沖是來自同一個輻射源的，因此可綜合同一串碼字里多個脈沖的個體識別結(jié)果，提高目標個體識別正確率。

2 針對敵我識別信號的個體識別方法

2.1 特征提取與比較

現(xiàn)有研究常提取的個體指紋特征通常包括包絡(luò)特征、頻率特征和時頻特征。脈沖上升沿可用來表示信號瞬態(tài)特征［12］，瞬時頻率能夠反映信號的脈內(nèi)頻率變化情況，將瞬時頻率與中心載頻作差的結(jié)果可以反映頻率穩(wěn)定度。

本文采用ADS-B 信號進行算法驗證。ADS-B信號是S模式敵我識別信號中的一種類型，在民航飛機上廣泛使用，因此ADS-B 信號容易獲取，為樣本采集提供了極大便利；另一方面，ADS-B 信號所含的地址碼是飛機的唯一身份標識，可以作為標簽對數(shù)據(jù)進行標注，進而可對數(shù)據(jù)樣本進行監(jiān)督學習。因此ADS-B 信號比較適合用來驗證針對敵我識別信號的個體識別算法。

雙向DC/DC變換器按結(jié)構(gòu)類型可分為非隔離型和隔離型兩種。隔離型雙向DC/DC變換器結(jié)構(gòu)復雜、電力元件多、控制復雜導致成本較高，功率密度較低，不適用于船舶的復合儲能裝置中。在非隔離式的Buck型拓撲、Boost型拓撲、Buck-Boost組合型拓撲三種拓撲結(jié)構(gòu)中，Buck-Boost結(jié)構(gòu)具有較小的電感和較低的開關(guān)導通損耗，效率也比另外兩種更高，更符合純電動船舶復合儲能裝置的需求，故選用雙向Buck-Boost型變換器。雙向DC/DC變換器匹配超級電容可以提供穩(wěn)定的電壓輸出，與動力電池并聯(lián)供電時，能夠有效減少環(huán)流問題。

我們提取了15部不同ADS-B 輻射源信號的脈沖上升沿、脈沖下降沿、瞬時頻率、頻率穩(wěn)定度作為個體指紋特征。在使用這些信號進行分類器訓練與識別之前，要進行必要的預處理措施，本文所用預處理措施簡要介紹如下:

缺失值處理:實際采集的敵我識別信號數(shù)據(jù)中，偶爾會出現(xiàn)某個采樣點上I 和Q 數(shù)據(jù)同時為0的情況，當對同時為0 的IQ 數(shù)據(jù)取其幅值（in dB），將出現(xiàn)log（0）的計算，從而導致計算出錯，因此必須將缺失值補全。

脈沖對齊:目前常用的脈沖對齊方法是以曲線擬合后的上升沿第一個過沖點頂部為基準，對脈沖進行對齊［6］。也可以用脈沖上升沿低于脈沖幅度3 dB（或6 dB）點作為基準，對脈沖進行對齊。通過對實際采集的ADS-B 信號進行測試，在數(shù)據(jù)樣本足夠大時，不同的對齊方法對識別結(jié)果影響不大。基于降低算法復雜度的考慮，本文采用后一種脈沖對齊方法。

標準化處理:不同輻射源個體的脈沖信號幅度不一樣，同一個輻射源個體不同時刻的脈沖幅度也不一樣，因此需要對脈沖包絡(luò)作標準化處理，消除具體幅值的影響，只關(guān)注包絡(luò)的變化特征。本文對脈沖包絡(luò)作0-1 標準化處理，使結(jié)果落到［0，1］區(qū)間，轉(zhuǎn)換函數(shù)如下:

式中，max為最大樣本值，min為最小樣本值。

圖3 顯示了標準化后的脈沖上升沿和下降沿包絡(luò)。選取了3 部不同的ADS-B 輻射源，用3 種不同顏色表示，每個輻射源選取100個脈沖的數(shù)據(jù)進行展示。從圖中可以看出不同輻射源上升沿和下降沿包絡(luò)在局部上是有些差異的。

圖3 標準化上升沿和下降沿包絡(luò)特征

圖4 顯示了3 部不同ADS-B 輻射源的瞬時頻率特征。可以看出，不同ADS-B 輻射源瞬時頻率特征差異比較明顯。

圖4 瞬時頻率特征

圖5 顯示了3 部不同ADS-B 輻射源的頻率穩(wěn)定度特征。可以看出，相較于瞬時頻率特征，不同ADS-B輻射源的頻率穩(wěn)定度特征差異明顯減小。

圖5 頻率穩(wěn)定度特征

2.2 多脈沖綜合識別

提取了輻射源個體信號特征后，利用不同輻射源個體間信號特征的差異，結(jié)合分類器對不同個體進行分類識別。近年來隨著機器學習成為研究熱點，多種分類器被用于個體識別領(lǐng)域，包括但不限于支持向量機（SVM）、隨機森林、神經(jīng)網(wǎng)絡(luò)、深度學習等算法［13］。本文并不打算研究不同分類器的個體識別效果，而是注重于對多個脈沖的識別結(jié)果進行融合，給出目標個體識別結(jié)果的綜合判決。

本文選取人工神經(jīng)網(wǎng)絡(luò)作為個體識別的分類器。人工神經(jīng)網(wǎng)絡(luò)由輸入層、隱藏層和輸出層組成。神經(jīng)網(wǎng)絡(luò)的預測輸出表達式為

式中:k表示第k個測試數(shù)據(jù)；j∈(1,2,…,M)，M為輸出層神經(jīng)元個數(shù)，代表有M個不同的類；βj為第j個神經(jīng)元的輸入，θj為第j個神經(jīng)元的閾值。f（）為激活函數(shù)，本文采用Sigmoid 函數(shù)作為激活函數(shù)，由于其取值范圍在0～1之間，因此可以將理解為第k個數(shù)據(jù)被判別為第j類的概率。

取最大的pkj對應的j作為第k個數(shù)據(jù)的識別結(jié)果，記為yk，即假設(shè)一串敵我識別信號碼字包含K個脈沖，每個脈沖都有一個個體識別結(jié)果yk∈(1,2,…,M)，其中M為輻射源個數(shù)，k∈(1,2,…,K)。對這K個脈沖的識別結(jié)果進行綜合處理，有兩種處理方法。

第一種，直接對最終識別結(jié)果進行綜合處理，在K個識別結(jié)果中，找到出現(xiàn)最多次的那一個結(jié)果作為對該輻射源的最終個體識別結(jié)果。

第二種，先對識別概率作綜合，再得到識別結(jié)果。記pj為該輻射源為第j類的概率，則有

然后取最大的pj對應的j作為該輻射源最終的個體識別結(jié)果y=arg max(pj)。

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

3.1 信號樣本準備與分類器訓練

本文利用敵我識別信號偵察設(shè)備采集了15架民航飛機的ADS-B 信號數(shù)據(jù)，每架飛機的脈沖樣本數(shù)為1 000～3 000不等。同一架飛機的ADS-B信號解調(diào)出的地址碼相同，不同飛機的地址碼不同，可以用地址碼給信號樣本打上標簽。同一個碼字一串脈沖的基帶IQ樣本數(shù)據(jù)存成一個文件，如圖6所示為一串ADS-B 信號的脈沖包絡(luò)，可以看出其信號格式與圖1所示相同。

圖6 ADS-B信號脈沖包絡(luò)

依次對IQ 數(shù)據(jù)進行缺失值處理、脈沖檢測、對齊、特征提取、標準化處理，得到脈沖包絡(luò)、瞬時頻率、頻率穩(wěn)定度等個體指紋特征數(shù)據(jù)，作為輸入樣本訓練神經(jīng)網(wǎng)絡(luò)。本文采用Matlab 自帶的神經(jīng)網(wǎng)絡(luò)模塊，構(gòu)建了一個包含輸入層、15 個隱含層和輸出層的神經(jīng)網(wǎng)絡(luò)，其中輸出層包含15 種輸出結(jié)果，分別表示15個不同的輻射源個體。特征數(shù)據(jù)按70∶15∶15 的比例隨機分為訓練集、驗證集和測試集，將訓練集和驗證集輸入到神經(jīng)網(wǎng)絡(luò)對其進行訓練，用測試集對訓練好的神經(jīng)網(wǎng)絡(luò)分類器進行性能測試。

3.2 不同特征識別效果比較

用提取的脈沖上升沿和下降沿包絡(luò)、脈沖瞬時頻率、頻率穩(wěn)定度，以及綜合這三種特征分別訓練神經(jīng)網(wǎng)絡(luò)分類器，并對其分類效果進行測試，試驗結(jié)果如圖7 所示。需要指出的是這里的測試結(jié)果是單脈沖的識別結(jié)果，并未對多脈沖結(jié)果作綜合。從測試結(jié)果可以看出瞬時頻率特征的識別效果要優(yōu)于脈沖包絡(luò)和頻率穩(wěn)定度特征，而綜合多種特征的個體識別效果表現(xiàn)最優(yōu)。另外可以看出輻射源個數(shù)越多，個體識別準確率越低，這是顯而易見的。

圖7 不同指紋特征的個體識別準確率

3.3 多脈沖綜合識別測試

敵我識別各種模式信號中，一串碼字對應的脈沖數(shù)最少的是MarkX 詢問和模式4應答模式，包含3 個脈沖；脈沖數(shù)最多的是S 模式應答（含ADSB），脈沖數(shù)不固定但最多可達112個脈沖。對不同脈沖個數(shù)的綜合識別準確率試驗結(jié)果如圖8所示。試驗條件為:輻射源個數(shù)為15 個，結(jié)合脈沖包絡(luò)、瞬時頻率和頻率穩(wěn)定度三種特征作為輸入特征。由圖8可以看出，多脈沖綜合識別正確率顯著高于單脈沖識別正確率，綜合50 個脈沖的識別正確率高達98.2%；而且僅僅綜合3 個脈沖的識別正確率就從81%（單脈沖）大幅提升至91.6%。同時，我們對3.2 節(jié)所描述的兩種綜合方法進行了比較，識別概率綜合方法要優(yōu)于識別結(jié)果綜合方法，但差距不大。

圖8 多脈沖綜合識別結(jié)果

4 結(jié)束語

本文針對敵我識別信號頻率和脈沖時序固定的特點，提出一種針對敵我識別系統(tǒng)的個體識別方法。通過分析與測試，發(fā)現(xiàn)瞬時頻率比較適合作為敵我識別系統(tǒng)個體指紋特征。同時，結(jié)合敵我識別信號解調(diào)，對同一個碼字包含的多個脈沖進行綜合識別，綜合多個脈沖的個體識別正確率顯著高于單脈沖個體識別正確率。后續(xù)工作將對具有脈內(nèi)調(diào)制的敵我識別信號，包括S模式詢問和模式5信號，開展個體識別研究。