何遵文，侯帥，張萬成，張焱

（北京理工大學(xué)信息與電子學(xué)院，北京 100081）

1 引言

特定輻射源識別（SEI,specific emitter identification）是將接收到的脈沖波形與唯一發(fā)射器進(jìn)行關(guān)聯(lián)[1]。5G 使數(shù)據(jù)傳輸速率和可靠性得到了跳躍式的提升，促進(jìn)了物聯(lián)網(wǎng)應(yīng)用的普及和發(fā)展。6G 將進(jìn)一步向萬物互聯(lián)的時代邁進(jìn)，接入物聯(lián)網(wǎng)的無線設(shè)備呈爆炸式增長，設(shè)備的安全性值得關(guān)注。物聯(lián)網(wǎng)中未知的輻射源設(shè)備需要能被準(zhǔn)確地識別身份，而電磁環(huán)境日益復(fù)雜，輻射源數(shù)目龐大，信號密度極高，這為SEI 技術(shù)帶來了更多的挑戰(zhàn)[2]。

用于個體識別的輻射源信號通常分為瞬態(tài)信號和穩(wěn)態(tài)信號[3]。瞬態(tài)信號是發(fā)射體開關(guān)打開時或狀態(tài)改變時的信號，但是瞬態(tài)信號持續(xù)時間較短，在實際通信環(huán)境特別是非合作通信情況下截獲難度較大。相比之下，穩(wěn)態(tài)信號更容易獲取，具備較強(qiáng)的操作性，基于穩(wěn)態(tài)信號的SEI 技術(shù)在實際應(yīng)用中具有更強(qiáng)的適用性[4]。

近年來，深度學(xué)習(xí)在SEI 領(lǐng)域有了廣泛的發(fā)展和應(yīng)用[5]，深度神經(jīng)網(wǎng)絡(luò)能夠通過具有非線性激活函數(shù)的多個隱含層來檢索抽象特征，有利于提取輻射源信號的深層次的細(xì)微特征[6]。因此，一些文獻(xiàn)使用基于深度神經(jīng)網(wǎng)絡(luò)的端到端的SEI 方法，如文獻(xiàn)[7-8]將原始I/Q 信號直接送入神經(jīng)網(wǎng)絡(luò)，完成端到端的SEI 任務(wù)，整體性強(qiáng)。但是這種方法過于依賴神經(jīng)網(wǎng)絡(luò)的設(shè)計，并且對不同的原始數(shù)據(jù)類型都要設(shè)計不同的神經(jīng)網(wǎng)絡(luò)，可擴(kuò)展性不強(qiáng)。另外，常用的卷積神經(jīng)網(wǎng)絡(luò)（CNN，convolutional neural network）更擅長識別二維圖像數(shù)據(jù)，直接處理I/Q數(shù)據(jù)效果會有一定的折扣。

因此，更有效和可靠的方案仍是將智能學(xué)習(xí)與信號處理技術(shù)相結(jié)合，將SEI 任務(wù)分為兩步，一是數(shù)據(jù)的預(yù)處理和特征提取，二是分類器設(shè)計和訓(xùn)練。通過數(shù)據(jù)預(yù)處理和變換以提取穩(wěn)態(tài)信號的細(xì)微特征，細(xì)微特征由器件內(nèi)部的物理特性和器件之間的相互作用產(chǎn)生，整個過程被稱為信號的非有意調(diào)制。信號細(xì)微特征也被稱為射頻指紋，通常表現(xiàn)為不可避免的非平穩(wěn)、非高斯、非線性[9]?，F(xiàn)有的研究主要分為利用信號的時域特征和變換域特征的方法。例如針對信號的時域特征，文獻(xiàn)[10]將瞬時振幅和頻率直接定義為指紋特征，文獻(xiàn)[11]利用經(jīng)驗?zāi)B(tài)分解（EMD,empirical mode decomposition）相關(guān)維數(shù)和小波熵提取信號的細(xì)微特征，文獻(xiàn)[12]基于盒維數(shù)和方差維數(shù)的分形特征進(jìn)行識別，文獻(xiàn)[13]基于實際信號的固有非線性動力學(xué)特性將排列熵提取為信號的射頻指紋。但是，這些利用時域波形的非線性特征的方法容易受到噪聲的影響，識別準(zhǔn)確率波動很大，并且計算量較大，具有一定的局限性?，F(xiàn)在更主流的是基于變換域特征的方法，時頻分析能夠提供信號的時頻聯(lián)合域信息，反映不同信號之間更加細(xì)微的差別[14]。文獻(xiàn)[15]提取動態(tài)小波指紋特征以識別各個射頻發(fā)射器。文獻(xiàn)[16-17]使用短時傅里葉變換時頻圖表示特征，然而，這種方法本質(zhì)上是線性的，對于非線性信號的分析效果并不好。希爾伯特?黃變換（HHT,Hilbert-Huang transform）是一種著名的處理非線性非平穩(wěn)信號的方法[18]，因此，HHT 被廣泛用于SEI 中。文獻(xiàn)[19]將希爾伯特光譜轉(zhuǎn)換為灰度圖像表示特征；文獻(xiàn)[20]使用HHT，并研究了它們在單跳和中繼場景下的各種渠道方案。除此之外，基于高階譜方法也是研究的熱點，高階譜可以保持信號的幅值和相位信息，抑制噪聲對非高斯信號的負(fù)面影響。文獻(xiàn)[21-22]充分利用了雙譜的特點進(jìn)行SEI 任務(wù)。

對于SEI 任務(wù)的第二步，即分類器的設(shè)計和訓(xùn)練，可以使用傳統(tǒng)機(jī)器學(xué)習(xí)方法和深度學(xué)習(xí)方法。文獻(xiàn)[14]使用支持向量機(jī)（SVM,support vector machine）作為分類器，文獻(xiàn)[13,23-24]在每次實驗的分類過程中，采用K 近鄰判別分類器。深度神經(jīng)網(wǎng)絡(luò)對二維圖像有很好的識別效果，因此一些文獻(xiàn)使用網(wǎng)絡(luò)學(xué)習(xí)能夠體現(xiàn)信號特征的二維圖像來完成SEI 任務(wù)，如文獻(xiàn)[19]構(gòu)造了一個深層的殘差網(wǎng)絡(luò)，用于學(xué)習(xí)希爾伯特光譜圖像中反映的視覺差異以辨別設(shè)備；文獻(xiàn)[24]利用CNN，將I/Q 信號轉(zhuǎn)換成輪廓Stella圖像進(jìn)行分類。

上述方法中，一方面，無論是使用基于時域特征還是變換域特征和高階譜的方法，其大多是使用單一信號處理方法提取其中的一種細(xì)微特征，而實際通信信號復(fù)雜多變，單一的信號特征并不足以全面準(zhǔn)確地表示輻射源信號之間的細(xì)微差別，最終識別準(zhǔn)確率受到限制，具有很大的局限性。另一方面，為了增強(qiáng)識別準(zhǔn)確率，多數(shù)模型使用深度神經(jīng)網(wǎng)絡(luò)作為強(qiáng)分類器；數(shù)據(jù)驅(qū)動的深度學(xué)習(xí)方法需要大量的數(shù)據(jù)，實際通信中受限于時間和人力成本，很難獲得充足的輻射源信號數(shù)據(jù)用于深度學(xué)習(xí)的訓(xùn)練，如果數(shù)據(jù)量過少，較深層和復(fù)雜的神經(jīng)網(wǎng)絡(luò)很容易發(fā)生過擬合，進(jìn)而嚴(yán)重影響最終的識別效果。

基于以上討論，本文提出了一種新型的多特征融合分類方法用于通信輻射源識別。針對信號單一特征的局限性，設(shè)計了基于多通道變換投影的特征提取方法，對原始數(shù)據(jù)集進(jìn)行小波、雙譜、HHT 等多種不同的信號變換，以實現(xiàn)原始信號不同維度細(xì)微特征的提取，豐富信號特征數(shù)據(jù)集的多樣性。將信號變換后的特征矩陣表示為三維張量，為了充分利用深度神經(jīng)網(wǎng)絡(luò)對二維圖像特征提取的優(yōu)勢，選取合適的視角得到二維投影。本文研究了基于多特征融合的雙階段識別分類方法，先利用不同的神經(jīng)網(wǎng)絡(luò)對不同的圖像特征集分別進(jìn)行學(xué)習(xí)和初始分類得到概率表示的特征向量，之后對初始分類的特征向量進(jìn)行融合學(xué)習(xí)，得到最終的分類結(jié)果。

此外，針對深度學(xué)習(xí)數(shù)據(jù)量不足的問題，本文利用改進(jìn)的輔助分類器（AC,auxiliary classifier）生成對抗網(wǎng)絡(luò)（GAN,generative adversarial network），對不同的圖像特征集進(jìn)行學(xué)習(xí)并擴(kuò)充增強(qiáng)數(shù)據(jù)集以滿足訓(xùn)練需求。基于實測數(shù)據(jù)集對所提方法進(jìn)行驗證，結(jié)果表明該方法相比其他方法具有更高的識別準(zhǔn)確率，在低信噪比下仍有較好的表現(xiàn)。

2 基于多通道變換的特征提取與數(shù)據(jù)增強(qiáng)

SEI 問題本質(zhì)上是一個信號分類問題，將接收機(jī)采集的無線通信信號輸入系統(tǒng)，最終得到這個信號的分類結(jié)果。所提出的SEI 多特征融合分類方法主要包括4 個部分：多特征提取、數(shù)據(jù)擴(kuò)充和增強(qiáng)、初級分類、集成融合，整體流程如圖1 所示。

首先，對信號進(jìn)行多特征提取得到多種不同的特征圖像，并在模型的訓(xùn)練階段使用GAN 對圖像特征集進(jìn)行擴(kuò)充，進(jìn)而訓(xùn)練分類器。在識別階段，直接將特征圖像送入分類器，分類器設(shè)計包括初級分類器和次級分類器兩級結(jié)構(gòu)，采用集成融合學(xué)習(xí)算法將多種初級分類器的識別結(jié)果進(jìn)行再學(xué)習(xí)，通過次級分類器得到最終識別結(jié)果。

考慮系統(tǒng)內(nèi)共有K個輻射源信號發(fā)射器，接收機(jī)接收來自K個發(fā)射機(jī)發(fā)出的理論上相同的信號。接收機(jī)接收到的通信輻射源信號可以表示為

算法1K個輻射源信號的特征提取與數(shù)據(jù)增強(qiáng)

特征提取

輸入對于每個類別的輻射源，輸入時域信號sk(t)的切片樣本，每1 000 個點作為一個信號樣本，每個類別的輻射源信號共n個樣本，共有K個輻射源

輸出對于每個類別的輻射源，輸出3 個二維特征圖像數(shù)據(jù)集，每個數(shù)據(jù)集有n張圖像，共有K個輻射源

1)一共n個樣本

2)對每個輸入樣本進(jìn)行離散小波變換得到小波系數(shù)特征矩陣，并表示為三維圖像形式

3)取投影視角(90,0,0)得到二維小波特征圖像

4)進(jìn)行雙譜變換得到雙譜系數(shù)特征矩陣，并表示為三維圖像形式

5)取投影視角(?45,0,0)得到二維雙譜特征圖像

6)進(jìn)行希爾伯特?黃變換得到HHT 系數(shù)特征矩陣，并表示為三維圖像形式

7)取投影視角(90,0,0)得到二維HHT 特征圖像

8)重復(fù)步驟1)～步驟7)K次，每個類別的輻射源信號都得到3 個二維特征圖像數(shù)據(jù)集，每個數(shù)據(jù)集有n張圖像

數(shù)據(jù)增強(qiáng)

輸入每次輸入一個特征投影圖像數(shù)據(jù)集，包含Kn張圖像，共3 個特征投影圖像數(shù)據(jù)集（小波、雙譜、HHT）

輸出每次輸出一個增強(qiáng)的數(shù)據(jù)集，最終輸出3 個特征圖像數(shù)據(jù)集，每個數(shù)據(jù)集包含K(m+n)張圖像

1)對每個輸入的特征投影圖像數(shù)據(jù)集

2)將圖像輸入AC-GAN 訓(xùn)練至收斂

3)得到m張模擬特征圖像作為補(bǔ)充數(shù)據(jù)集

4)擴(kuò)充數(shù)據(jù)集與原始數(shù)據(jù)集共同構(gòu)成增強(qiáng)數(shù)據(jù)集

5)得到3 個增強(qiáng)圖像數(shù)據(jù)集（小波、雙譜、HHT），每個數(shù)據(jù)集包含K(m+n)張圖像

2.1 基于多通道變換投影的特征提取

單一的信號特征很難準(zhǔn)確表示輻射源信號之間的差別，而過多的特征融合可能會影響模型性能。因此選取小波變換、雙譜變換、HHT 這3 種有代表性并被證明有效的信號處理方法提取3 種不同的信號變換域特征，小波變換是時間（空間）頻率的局部化分析，更能關(guān)注信號的區(qū)域性信息[25]；雙譜變換可以更多地關(guān)注信號的相位信息，同時可以用來描述信號的不對稱程度和非線性程度[21]；HHT理論上最適合處理非線性非平穩(wěn)信號，且更能關(guān)注一些突變信號的細(xì)微特征，提供有關(guān)信號的非線性和非平穩(wěn)特性的大量信息[18]。以上3 種特征可以從不同維度利用各自的優(yōu)勢，豐富信號特征數(shù)據(jù)集的多樣性。二維圖像更有利于深度學(xué)習(xí)和分類，因此，本文設(shè)計一種變換投影的方法對信號變換后的三維特征圖像選取合適的視角進(jìn)行投影，進(jìn)而將其表示為二維特征圖像。

2.1.1 基于小波的特征提取

小波變換（WT,wavelet transform）通過對時間（空間）頻率進(jìn)行局部化分析[25]，能夠凸顯信號區(qū)域性的細(xì)微特征，在個體識別領(lǐng)域有著獨特的優(yōu)勢。Morlet 小波為復(fù)值小波，在時域和頻域都有很好的聚集性，本文采集的實測數(shù)據(jù)均為復(fù)信號，因此選取Morlet 小波作為基小波函數(shù)用于時頻分析。

利用尺度參數(shù)scale 來衡量小波的頻率f，兩者之間的轉(zhuǎn)換關(guān)系為

其中，F(xiàn)s為信號的采樣頻率，F(xiàn)c為小波的中心頻率。為使小波尺度圖的頻率范圍為(0,Fs/2)，尺度范圍應(yīng)為(2Fc,+∞)。實際中，尺度應(yīng)為等差數(shù)列形式，即scale=2Fc,…,2Fctotalscale，其中totalscale是預(yù)先設(shè)定好對信號進(jìn)行小波變換時所用尺度序列的長度。完成參數(shù)設(shè)定后，小波變換的過程可以表示為

其中，sk(t)為原始時間信號；為小波基函數(shù)，具有平移和伸縮特性；scale 為尺度參數(shù)，控制小波函數(shù)的伸縮；τ為平移參數(shù)，控制小波函數(shù)的平移。尺度對應(yīng)于頻率，平移量τ對應(yīng)于時間。用表示得到的頻率?時間?小波系數(shù)矩陣，三維圖像如圖2(a)所示。

2.1.2 基于雙譜的特征提取

雙譜的物理定義可以近似為信號偏態(tài)在頻域上的分解，用來描述信號的不對稱性和非線性程度，測量信號偏離高斯分布的程度。在實際通信過程中，輻射源個體的一些細(xì)微特征更多以非平穩(wěn)、非線性、非高斯的形式附加在信號中，例如由于發(fā)射系統(tǒng)中使用的電路和設(shè)備的不同而產(chǎn)生的不同的發(fā)射端噪聲，可以作為區(qū)分輻射源信號的指紋特征。接收系統(tǒng)中大量散射體和熱噪聲形成的雜波呈高斯分布，通過雙譜變換，可以有效地抑制一定數(shù)量的高斯噪聲，使發(fā)射機(jī)噪聲充分反映在雙譜信息中。雙譜定義為

2.1.3 基于HHT 的特征提取

HHT 相比傅里葉變換和小波變換這類傳統(tǒng)方法，徹底擺脫了線性和平穩(wěn)性束縛，更適用于分析非線性非平穩(wěn)信號[20]。為了使瞬時頻率具有明顯的物理意義，HHT 首先需要進(jìn)行經(jīng)驗?zāi)B(tài)分解，也通常被稱為一個篩選過程，依據(jù)信號特點自適應(yīng)地把任意一個復(fù)雜信號分解為一系列本征模態(tài)函數(shù)（IMF,intrinsic mode function）[18]。這個過程可以表示為

原始輸入信號sk(t)被表示成一系列IMF 分量ci(t)和一個殘余項rn(t)的組合。之后對每個IMF 分量進(jìn)行希爾伯特變換來提供時頻分布

構(gòu)造解析函數(shù)

其中，ai(t)、?i(t)和ωi(t)分別表示瞬時振幅、相位和頻率，計算式分別為

進(jìn)而希爾伯特譜可以表示為

2.1.4 信號特征表示

圖2 信號變換生成的特征圖像

原始信號經(jīng)過3 種變換得到了三維特征矩陣，為了更好地表示特征并利于后續(xù)進(jìn)行深度學(xué)習(xí)，將三維特征矩陣取投影生成二維特征圖像。同一類別的不同輻射源信號樣本是基于同一信號的不同時間切片生成的，而不同時間切片的信號差異可能會影響同類信號樣本的相似度。因此對于小波和希爾伯特時頻特征矩陣，取投影視角(90,0,0)表示投影角度為90°將三維圖像投影在YZ平面上，淡化時間維度的影響，更加突出頻率和系數(shù)的關(guān)系，這樣更能保證同類信號樣本之間的相似度。對于雙譜三維特征矩陣，雙譜在2 個頻率主軸形成的平面上其第一、三象限包含的內(nèi)容較多，根據(jù)雙譜的對稱性，投影視角(?45,0,0)包含三維雙譜圖像的輪廓、紋理特征最多，且完全對稱，是對2 個頻率主軸形成的平面上第一、三象限內(nèi)容的完全體現(xiàn)，因此選用投影視角(?45,0,0)表征的雙譜特征最完備[27]。

最終得到3 種特征的二維投影圖像Iik，分別如圖2(b)、圖2(d)和圖2(f)所示。

2.2 基于AC-GAN 的數(shù)據(jù)增強(qiáng)擴(kuò)充

后續(xù)用于識別的初級分類器使用深度殘差網(wǎng)絡(luò)，網(wǎng)絡(luò)層數(shù)和參數(shù)較多，因此在模型的訓(xùn)練階段需要大量的數(shù)據(jù)提供給網(wǎng)絡(luò)學(xué)習(xí)，受限于信號接收采集的時間長度等因素，原始數(shù)據(jù)量難以滿足訓(xùn)練需求。

因此，所提方法使用AC-GAN 來進(jìn)行數(shù)據(jù)增強(qiáng)，通過這種方式能夠?qū)崿F(xiàn)基于有限樣本的數(shù)據(jù)擴(kuò)充，生成和實測數(shù)據(jù)具有相似分布與特征的大量訓(xùn)練數(shù)據(jù)，從而為后級分類識別提供支撐。增強(qiáng)后的數(shù)據(jù)集可表示為Di′。傳統(tǒng)GAN 只能生成單一類別的圖像，AC-GAN 將類別標(biāo)簽作為條件信息和噪聲一起放入生成器后生成樣本，判別器在判斷輸入樣本是否為真的同時，利用輔助分類器來判斷輸入樣本的所屬類別，因此AC-GAN 可以生成不同類別的圖像。使用3 個AC-GAN 對3 種投影圖像數(shù)據(jù)集進(jìn)行擴(kuò)充。AC-GAN 的結(jié)構(gòu)和工作流程如圖3 所示。

圖3 AC-GAN 的結(jié)構(gòu)和工作流程

生成器有2 個輸入，一個是投影圖像所屬的輻射源類別標(biāo)簽k，k=0,…,K?1；另一個是隨機(jī)噪聲數(shù)據(jù)，用于生成模擬的投影圖像Ifake。網(wǎng)絡(luò)可以自動學(xué)習(xí)3種不同投影圖像的內(nèi)在特征，進(jìn)而生成對應(yīng)的模擬投影圖像。判別器分別給出圖像（真假）的概率分布P(S|I)以及分類標(biāo)簽的概率分布P(K|I)。

AC-GAN 的目標(biāo)函數(shù)包含2 個部分。第一部分是面向輸入的投影圖像真實與否的代價函數(shù)LS

第二部分是對輸入的投影圖像分類準(zhǔn)確性的代價函數(shù)LK

在優(yōu)化過程中，希望判別器能夠盡可能地區(qū)分原始真實的投影圖像和模擬的投影圖像，并且能有效分類，即令LS+LK盡可能地最大。另外希望生成器生成的模擬投影圖像被盡可能認(rèn)為是真實的，即令LK?LS盡可能地最大。迭代數(shù)次后，判別器對生成器生成的模擬投影圖像識別準(zhǔn)確率趨于穩(wěn)定，模擬的投影圖像已經(jīng)可以滿足數(shù)據(jù)要求，此時生成器生成的模擬投影圖像構(gòu)成補(bǔ)充數(shù)據(jù)集，與原始數(shù)據(jù)集一起構(gòu)成擴(kuò)充后的增強(qiáng)數(shù)據(jù)集。

3 多特征融合雙階段分類識別

通過多特征提取與增強(qiáng)步驟，基于小波變換、雙譜變換和HHT 分別獲取了原始輻射源信號的某個方面的特征，并將信號特征轉(zhuǎn)換為了二維圖像?；贏C-GAN 進(jìn)一步實現(xiàn)了數(shù)據(jù)集的增強(qiáng)?；谠鰪?qiáng)后的特征數(shù)據(jù)集，在分類識別階段采用雙階段分類識別方法，融合多種特征識別結(jié)果以提升SEI 的準(zhǔn)確率。

算法2模型訓(xùn)練和分類識別

模型訓(xùn)練

輸入每次輸入一個增強(qiáng)圖像數(shù)據(jù)集，一共3個增強(qiáng)圖像數(shù)據(jù)集

輸出每次輸出對每張圖像的預(yù)測結(jié)果，用于訓(xùn)練反饋

1)將小波特征數(shù)據(jù)集輸入殘差網(wǎng)絡(luò)1 進(jìn)行訓(xùn)練，輸出初始分類結(jié)果Pcwt，并保存模型參數(shù)

2)將雙譜特征數(shù)據(jù)集輸入殘差網(wǎng)絡(luò)2 進(jìn)行訓(xùn)練，輸出初始分類結(jié)果Pbis，并保存模型參數(shù)

3)將HHT 特征數(shù)據(jù)集輸入殘差網(wǎng)絡(luò)3 進(jìn)行訓(xùn)練，輸出初始分類結(jié)果Phht，并保存模型參數(shù)

4)將Pcwt、Pbis和Phht輸入次級分類器進(jìn)行訓(xùn)練，保存模型參數(shù)

分類識別

輸入每次輸入待識別信號的一張?zhí)卣鲌D像，一共3 張?zhí)卣鲌D像

輸出每次輸出對待識別信號的最終分類預(yù)測結(jié)果Pfinal

1)將待識別信號的小波特征圖像輸入殘差網(wǎng)絡(luò)1 輸出初始預(yù)測結(jié)果Pcwt

2)將待識別信號的雙譜特征圖像輸入殘差網(wǎng)絡(luò)2 輸出初始預(yù)測結(jié)果Pbis

3)將待識別信號的HHT 特征圖像輸入殘差網(wǎng)絡(luò)3 輸出初始預(yù)測結(jié)果Phht

4)將Pcwt、Pbis和Phht輸入次級分類器進(jìn)行最終識別，輸出最終分類預(yù)測結(jié)果Pfinal

3.1 基于殘差網(wǎng)絡(luò)的初級分類

把原始信號特征轉(zhuǎn)化為二維圖像表征，分類器實現(xiàn)方面選用了當(dāng)前在圖像分類領(lǐng)域較成熟的深度學(xué)習(xí)模型ResNet18，殘差網(wǎng)絡(luò)ResNet因其高效和易優(yōu)化而被廣泛利用在深度學(xué)習(xí)的研究中[28]，網(wǎng)絡(luò)的參數(shù)經(jīng)過公開數(shù)據(jù)集ImageNet 預(yù)訓(xùn)練，這樣不僅可以進(jìn)一步避免過擬合等問題，還充分利用了網(wǎng)絡(luò)中已有的參數(shù)，使網(wǎng)絡(luò)具有更好的適應(yīng)性。更改網(wǎng)絡(luò)的全連接層，將輸出類別更改為K個，輸入增強(qiáng)數(shù)據(jù)集D′j訓(xùn)練網(wǎng)絡(luò)全連接層的權(quán)重，完成訓(xùn)練保存完整模型。在用網(wǎng)絡(luò)進(jìn)行SEI識別時，將圖像Iik輸入網(wǎng)絡(luò)，網(wǎng)絡(luò)的輸出是Softmax 層輸出的各個類別概率，可以表示為Pi=[p0,p1,…,pk]。3 個網(wǎng)絡(luò)輸出的預(yù)測分別為Pcwt、Pbis和Phht。

3.2 基于集成融合的次級分類

次級分類器需要完成最終的分類，使用集成學(xué)習(xí)中的Stacking 方法，將3 個神經(jīng)網(wǎng)絡(luò)分類器的初級分類結(jié)果進(jìn)行融合，即使某一個初級分類器預(yù)測錯誤，也能通過另外2 個分類器將錯誤糾正回來，最大化整體模型的分類能力。為了降低模型整體過擬合的風(fēng)險，次級學(xué)習(xí)器選擇結(jié)構(gòu)相對簡單的線性回歸模型（LRM,linear regression model），降低模型的整體復(fù)雜度，減少計算開支。

3 個初級分類器基于完整的訓(xùn)練集進(jìn)行訓(xùn)練，為了防止過擬合，每一個初級分類器需要通過5 折交叉檢驗的方式進(jìn)行訓(xùn)練。將訓(xùn)練集樣本隨機(jī)分為5 份，每折訓(xùn)練取其中4 份對一個初級分類器進(jìn)行訓(xùn)練，取剩下的一份數(shù)據(jù)進(jìn)行預(yù)測分類，保留分類結(jié)果。同時也對測試集所有樣本進(jìn)行分類，直到5 份訓(xùn)練集樣本都被分類，得到該殘差網(wǎng)絡(luò)對所有訓(xùn)練集樣本的分類結(jié)果，此時測試集樣本被分類了5 次，最終對分類結(jié)果取平均值。訓(xùn)練集和測試集的使用流程如圖4 所示。

訓(xùn)練集和測試集中每個樣本的分類結(jié)果都表示為Pi=[p0,p1,…,pk]。將3 個初級分類器對原訓(xùn)練集的分類結(jié)果整合在一起，用作次級分類器的訓(xùn)練樣本，而測試集的分類結(jié)果用作新的測試樣本對次級分類器進(jìn)行測試。對同一個樣本，3 個初級分類器有3 個分類結(jié)果，次級分類器分配這3 個結(jié)果的權(quán)重以進(jìn)行取舍，得到最終分類結(jié)果為Pfinal。

4 實驗測試和結(jié)果分析

4.1 實驗環(huán)境配置和數(shù)據(jù)準(zhǔn)備

為了驗證所提SEI 框架和方法的可行性，并使應(yīng)用場景更符合實際，本文選用了5 個型號為Tenda U6的USB 無線網(wǎng)卡作為輻射源發(fā)射器。將無線網(wǎng)卡安裝在電腦主機(jī)作為信號發(fā)射端，工作頻率為2.462 GHz，帶寬為40 MHz，最大發(fā)射功率為20 dBmW，均使用IEEE 802.11b 協(xié)議。接收機(jī)使用NI-PXI8106機(jī)柜和NI 矢量信號分析儀連接天線，固定于接收端15 cm 處采集無線網(wǎng)卡信號，采集的信號為I/Q 雙路復(fù)信號，信號采樣率為100 MSPS（million sample per second）。為了對比模型在不同信噪比（SNR,signal-to-noise ratio）下的性能，對原始信號加入高斯白噪聲模擬噪聲環(huán)境，得到5～25 dB 每隔5 dB 的不同SNR 條件下的無線網(wǎng)卡信號采樣數(shù)據(jù)。

將每1 000 個點的信號采樣數(shù)據(jù)作為一個輻射源信號樣本，最終識別準(zhǔn)確率定義為用于測試的正確分類的樣本數(shù)量在整體測試樣本數(shù)量中的占比。

實驗在MATLAB 環(huán)境進(jìn)行信號處理工作，深度學(xué)習(xí)環(huán)境為Python3.6+Pytorch1.0。

4.2 實驗結(jié)果分析

4.2.1 多層融合分類方法的整體表現(xiàn)

首先對整個多層融合框架的性能表現(xiàn)進(jìn)行測試。選取一種單一特征提取方法[19]、2 種多特征提取方法[4,23]與本文方法進(jìn)行對比。文獻(xiàn)[19]作為對比方法1 僅采用HHT 方法進(jìn)行識別。與本文方法類似，文獻(xiàn)[4]同樣使用多特征提取思想，提取信號的能量分布信息、分形特征和雙譜特征等多種雜散特征，并利用SVM 進(jìn)行分類，以此作為對比方法2。文獻(xiàn)[23]利用PCA 算法對樣本進(jìn)行降維，實現(xiàn)特征提取，然后利用K 近鄰、隨機(jī)森林和神經(jīng)網(wǎng)絡(luò)構(gòu)建混合分類模型，同樣使用多分類器融合的設(shè)計，以此作為對比方法3。結(jié)果如表1 和圖5 所示。

圖4 訓(xùn)練集和測試集的使用流程

可以看出，本文方法的最終識別率最高，在識別準(zhǔn)確率和抗噪表現(xiàn)上都優(yōu)于其他方法。對比方法1使用單一特征提取算法HHT，整體表現(xiàn)不如融合模型；對比方法2 雖然同樣使用多特征提取，時域方法在特征的選擇上具有局限性，且其分類器使用SVM 的分類效果弱于深度神經(jīng)網(wǎng)絡(luò)對圖像的分類效果；對比方法3 雖然使用了多分類器集成，但使用PCA 方法降維并未充分提取表示信號差別的細(xì)微特征，整體表現(xiàn)不佳。本文方法使用基于變換域的多特征提取，并使用深度神經(jīng)網(wǎng)絡(luò)集成識別特征圖像的方法，具有更好的性能表現(xiàn)。

表1 實測不同方法的輻射源識別準(zhǔn)確率

圖5 不同方法在不同信噪比下的識別效果對比

4.2.2 使用AC-GAN 擴(kuò)充數(shù)據(jù)的影響

比較使用AC-GAN 進(jìn)行不同比例的數(shù)據(jù)擴(kuò)充的最終識別準(zhǔn)確率，按擴(kuò)充數(shù)據(jù)占原始數(shù)據(jù)的比例為0～90%進(jìn)行對比試驗，結(jié)果如圖6 所示。從圖6可以看出，當(dāng)數(shù)據(jù)擴(kuò)充比例為50%時，模型的識別精度不再隨著數(shù)據(jù)擴(kuò)充比例的增加而明顯提高，而是逐漸趨于穩(wěn)定。使用AC-GAN 擴(kuò)充數(shù)據(jù)集的模型在抗噪表現(xiàn)和識別準(zhǔn)確率都表現(xiàn)更好，特別是在高信噪比的情況下，模型整體識別效果有顯著提升，說明使用AC-GAN 對數(shù)據(jù)集進(jìn)行增強(qiáng)，能夠提升整體識別準(zhǔn)確率，從而改善數(shù)據(jù)量不充足條件下模型整體識別率不高的問題。

圖6 模型使用AC-GAN 進(jìn)行不同比例的數(shù)據(jù)擴(kuò)充的最終識別準(zhǔn)確率

4.2.3 衰落信道對信號識別的影響

本節(jié)使用室外典型場景多徑衰落信道模型對輻射源信號進(jìn)行了處理，選用瑞利衰落信道模型，將輻射源產(chǎn)生的信號通過多徑傳播模型以模擬其在室外場景下的傳播特性[29]。由圖7 可以看出，衰落信道環(huán)境對輻射源識別效果產(chǎn)生了負(fù)面影響，但模型仍可較好地完成分類識別任務(wù)，即使是在瑞利衰落信道這種相對惡劣的信道環(huán)境下，本文方法在20 dB 信噪比條件下的識別準(zhǔn)確率也有80%以上，相比其他對比方法，本文方法有明顯的優(yōu)勢。通過上述結(jié)果表明，本文方法能夠適用于復(fù)雜無線信道環(huán)境的應(yīng)用。

5 結(jié)束語

本文提出了一種基于多特征提取、集成融合學(xué)習(xí)和生成對抗網(wǎng)絡(luò)的通信輻射源融合分類識別方法。通過對原始信號進(jìn)行小波變換、雙譜變換、HHT 等信號變換處理，得到3 種可反映信號細(xì)微特征的系數(shù)矩陣，對系數(shù)矩陣的三維圖像進(jìn)行二維投影以得到特征圖像數(shù)據(jù)集，并采用AC-GAN對圖像數(shù)據(jù)集進(jìn)行增強(qiáng)。此外，本文設(shè)計了一種多特征融合雙階段分類識別方法，利用多個殘差網(wǎng)絡(luò)對不同圖像特征集分別進(jìn)行訓(xùn)練學(xué)習(xí)，得到初始分類結(jié)果，并通過融合得到最終分類結(jié)果。結(jié)果表明，相對于其他相關(guān)的SEI 方法，本文方法在不同信噪比下的平均識別率較高，即使在低信噪比條件下也有較好的表現(xiàn)。通過引入基于AC-GAN 的數(shù)據(jù)增強(qiáng)機(jī)制，能夠提升整體識別準(zhǔn)確率，改善數(shù)據(jù)量不足而影響識別效果的問題。

圖7 在衰落信道環(huán)境下輻射源信號的識別效果對比