?

時(shí)頻圖像二維EMD分解在LPI信號(hào)識(shí)別中的應(yīng)用*1

張帆1，刁鳴1，楊承志1,2

(1.哈爾濱工程大學(xué) 信息與通信工程學(xué)院，黑龍江 哈爾濱150001;

2. 空軍航空大學(xué) 信息對(duì)抗系，吉林 長春130022)

摘要：針對(duì)LPI信號(hào)分類識(shí)別問題中，時(shí)頻圖像受噪聲干擾嚴(yán)重的問題，提出了一種基于二維快速經(jīng)驗(yàn)?zāi)Ｊ椒纸?FBEMD)的圖像降噪算法，并利用該算法實(shí)現(xiàn)對(duì)LPI信號(hào)的分類。首先利用時(shí)頻分析方法，獲得待分類信號(hào)的時(shí)頻分布圖像；使用二維EMD分解算法對(duì)圖像降噪；截取包含時(shí)頻信息的圖像部分，通過主分量分析法提取特征矢量；最后采用RBF神經(jīng)網(wǎng)絡(luò)完成信號(hào)的分類識(shí)別任務(wù)。對(duì)常見的LPI雷達(dá)信號(hào)進(jìn)行仿真，結(jié)果表明較低信噪比情況下，該方法仍能獲得較好的分類結(jié)果。當(dāng)信噪比為-2 dB時(shí)，采用二維EMD降噪算法，平均正確識(shí)別率能夠達(dá)到93%。

關(guān)鍵詞：LPI雷達(dá)信號(hào)；時(shí)頻分布；二維EMD分解；主分量分析

0引言

近年來，雷達(dá)技術(shù)不斷發(fā)展，信號(hào)頻率范圍更大，信號(hào)調(diào)制特征更復(fù)雜，傳統(tǒng)雷達(dá)信號(hào)的五大參數(shù)難以有效描述雷達(dá)信號(hào)，信號(hào)識(shí)別要求難以滿足。同時(shí)在復(fù)雜的戰(zhàn)場(chǎng)電磁環(huán)境下，截獲的信號(hào)必然存在大量噪聲干擾，較低的信噪比使得雷達(dá)信號(hào)的識(shí)別更加困難。

近年來，利用時(shí)頻分析技術(shù)與數(shù)字圖像處理技術(shù)相結(jié)合，實(shí)現(xiàn)對(duì)LPI信號(hào)的識(shí)別和分類的方法，引起了大量學(xué)者的關(guān)注。文獻(xiàn)[1-2]使用了一種將時(shí)頻圖像二值化，利用圖像處理形態(tài)學(xué)中的腐蝕和膨脹降噪，完成分類的算法；文獻(xiàn)[3]提出一種利用二維高斯低通濾波器處理時(shí)頻圖像，完成分類的算法；文獻(xiàn)[4]利用自適應(yīng)維納濾波對(duì)圖像處理，二值化后提取局部二值模式紋理特征完成分類。還有大量研究，重點(diǎn)關(guān)注對(duì)時(shí)頻圖像特征的提取，如利用偽Zernik矩陣[2]、時(shí)頻圖像代數(shù)特征[5]、二維主分量分析[6]等。這些算法中，對(duì)LPI(low probability of intercept)雷達(dá)信號(hào)時(shí)頻圖像的降噪處理研究不深，降噪同時(shí)無法有效保留圖像的高頻細(xì)節(jié)信息，對(duì)分類正確率有很大的影響。

2003年J. C Nunes.[7]將基于經(jīng)驗(yàn)的模式分解算法(empirical mode decomposition，EMD)推廣到二維領(lǐng)域，提出了二維EMD算法。近幾年來，二維EMD分解被大量學(xué)者深入研究，廣泛應(yīng)用于圖像去噪[8]、圖像融合[9]等領(lǐng)域。

針對(duì)上述問題，本文提出一種利用時(shí)頻圖像二維EMD分解和重構(gòu)算法的LPI信號(hào)分類方法。

1基于時(shí)頻圖像處理的LPI信號(hào)分類方法

本文首先通過二維EMD分解和重構(gòu)對(duì)獲取的LPI雷達(dá)信號(hào)的時(shí)頻圖像進(jìn)行預(yù)處理；然后利用奇異值分解與主分量分析相結(jié)合的方法提取時(shí)頻圖像特征矢量， 根據(jù)特征矢量完成對(duì)LPI雷達(dá)信號(hào)的調(diào)

制識(shí)別。具體流程如圖1所示。

1.1時(shí)頻分析

時(shí)頻分析能把時(shí)域信號(hào)映射到二維的時(shí)頻平面上，能夠同時(shí)從時(shí)域和頻域描述信號(hào)，對(duì)信號(hào)描述更加精確。Wigner-Ville分布(WVD)是一種具有邊緣特性、時(shí)移頻移不變性、實(shí)值性等特性的時(shí)頻分析方法，而且對(duì)線性調(diào)頻信號(hào)有著良好的能量集中性。因此，文中選擇了加窗的偽Wigner-Ville分布(PWVD)進(jìn)行分析,其定義為

(1)

式中：s(t)為受噪聲干擾的信號(hào)；h(t)為窗函數(shù)。本文采用了7種LPI雷達(dá)信號(hào)，包括BPSK信號(hào)、COSTAS信號(hào)、FRANK信號(hào)、FMCW信號(hào)、P1信號(hào)、P2信號(hào)、PT1信號(hào)， 其PWVD時(shí)頻圖像如圖2所示。從時(shí)頻圖像中，可以發(fā)現(xiàn)有明顯的交叉項(xiàng)，但本文采用該種時(shí)頻分析方法，主要用于完成對(duì)不同降噪算法的對(duì)比，故本文中不對(duì)交叉項(xiàng)的抑制進(jìn)行詳細(xì)研究。

1.2基于二維快速EMD算法的時(shí)頻圖像預(yù)處理算法

根據(jù)對(duì)LPI雷達(dá)信號(hào)的PWVD結(jié)果的分析，可以發(fā)現(xiàn)時(shí)頻圖像中，不僅包含與信號(hào)有關(guān)的時(shí)域和頻域特征，也包含大量的高頻噪聲、無信號(hào)區(qū)域等冗余信息。噪聲干擾使分類能力降低；冗余信息使得時(shí)頻圖像數(shù)據(jù)尺寸過大，直接進(jìn)行分類大大增加計(jì)算量。因此，降低噪聲干擾、去除冗余信息是對(duì)LPI雷達(dá)信號(hào)正確、迅速分類必須完成的任務(wù)。

1.2.1二維快速EMD算法

完成EMD分解算法，需要首先找到最佳的局部均值。傳統(tǒng)EMD分解確定信號(hào)的極大極小值點(diǎn)，通過插值算法構(gòu)成上下包絡(luò)線，最后利用上下包絡(luò)線求均值?？焖俣S經(jīng)驗(yàn)?zāi)Ｊ椒纸馑惴?，基于矩陣變換，減少了求極值的計(jì)算，提高了EMD算法效率。FBEMD(two-dimensional fast EMD)的核心思想是利用均值濾波矩陣組求矩陣鄰域均值[10]：

圖1　基于圖像處理的LPI信號(hào)分類流程Fig.1　LPI signal classification process based on image processing

(2)

當(dāng)x=1時(shí)，P,Q為式(2)所示的矩陣，令k=r=1/3,P▽Q可得矩陣每一點(diǎn)的三鄰域均值。利用該三鄰域均值表示傳統(tǒng)二維EMD算法的均值包絡(luò)曲面，可以穩(wěn)定迅速地實(shí)現(xiàn)二維EMD算法，在速度上比傳統(tǒng)算法快2～3個(gè)數(shù)量級(jí)。

2.2.2基于二維快速EMD算法的時(shí)頻圖像降噪算法

高斯白噪聲經(jīng)過EMD獲得的每個(gè)IMF仍高度

近似服從高斯分布；對(duì)噪聲圖像進(jìn)行二維EMD分解，獲得多個(gè)IMF，第1個(gè)IMF主要由噪聲信號(hào)和非常少量的高頻圖像細(xì)節(jié)信息組成[11]。利用這一特點(diǎn)，可以對(duì)原始時(shí)頻圖像進(jìn)行自適應(yīng)的降噪重構(gòu)。將待分解的時(shí)頻圖像表示為

(3)

式中：IMFi∈Rm×n是第i次分解的模態(tài)函數(shù)分量，Rl∈Rm×n是l次分解后的剩余量。FBEMD降噪算法可描述為：

(1) 初始化：i=1,IMFi=A,R1=A；

(2) 計(jì)算均值圖像Hi=P×IMFi×QT，從而求得IMFi=Ri-Hi,Ri+1=Hi；

(3) 重復(fù)步驟(2)操作，直至獲得合適的IMF分量；

(4) 將最后一次的剩余量賦值給Rl.

(5) 重構(gòu)時(shí)頻圖像矩陣為

(4)

圖3表示SNR=-6 dB情況下，F(xiàn)MCW信號(hào)的時(shí)頻圖像，其中圖3a)為原始PWVD時(shí)頻圖像；圖3b)為降噪重構(gòu)后的PWVD時(shí)頻圖；圖3c)為圖像二維EMD分量IMF1。經(jīng)過二維EMD降噪算法，時(shí)頻圖像噪聲得到有效抑制，時(shí)頻集聚性更強(qiáng)，有利于圖像特征的提取。

圖2　不同LPI雷達(dá)信號(hào)PWVD時(shí)頻分布圖像Fig.2　PWVD distribution images of different LPI radar signals

圖3　SNR=-6 dB時(shí)FMCW信號(hào)的時(shí)頻圖像Fig.3　PWVD distribution images of FMCW signals when SNR= -6 dB

1.3時(shí)頻圖像特征提取

1.3.1圖像裁剪

時(shí)頻圖像在經(jīng)過二維EMD分解與重構(gòu)后，噪聲得到有效抑制，但是仍然存在大量冗余信息，首先需要對(duì)圖像進(jìn)一步提取和裁剪。圖像裁剪步驟如下：

(1) 計(jì)算邊際頻率分布，邊際頻率分布將信號(hào)的瞬時(shí)能量表示為頻率的函數(shù)[3]

(5)

即對(duì)時(shí)頻圖像各頻率的時(shí)間值相加，儲(chǔ)存為列矢量。

(2) 歸一化邊際頻率分布，利用直方圖計(jì)算獲得閾值，確定截取圖像的起始和截止頻率。

(3) 根據(jù)起始頻率和截止頻率，截取圖像，然后將圖像調(diào)整為固定大小50×200，圖4為SNR=-6 dB情況下，F(xiàn)MCW信號(hào)截取后的灰度圖像。

(4) 將圖像轉(zhuǎn)化為10 000×1的列矢量，用于后續(xù)的特征提取。

圖4　SNR=-6 dB時(shí)FMCW信號(hào)截取的灰度圖像Fig.4　Grayscale image of FMCW signal whenSNR=-6 dB

1.3.2圖像特征提取

圖像裁剪后，時(shí)頻圖像的冗余信息已經(jīng)得到很大程度的消除，但是10 000×1的列矢量仍屬于高維數(shù)據(jù)矢量，不適用于實(shí)際工程。因此，進(jìn)一步通過主分量分析(principal component analysis,PCA)和奇異值分解相結(jié)合的方式對(duì)數(shù)據(jù)降維[12]，具體過程如下：

(1) 利用一系列訓(xùn)練信號(hào)，經(jīng)過圖像預(yù)處理后構(gòu)建訓(xùn)練矩陣X,X的維度為N×P。

(2) 對(duì)訓(xùn)練矩陣X進(jìn)行奇異值分解，得X=UWVH，則U=XVW-1，其中U為N×N酉矩陣，V為P×P酉矩陣，W為非負(fù)奇異值構(gòu)成的N×P矩陣。

(3) 構(gòu)建投影矩陣A，選取特征矩陣U中低于閾值的特征值置0，閾值選取最大特征值的0.001倍，建立一個(gè)非零特征值降序排列的投影矩陣A。

(4) 最后利用公式Y(jié)=AHX，獲得訓(xùn)練矩陣和測(cè)試信號(hào)的低維空間投影，完成圖像特征的提取。

2仿真實(shí)驗(yàn)與結(jié)果分析

徑向基函數(shù)(radial basis function,RBF)神經(jīng)網(wǎng)絡(luò)，由于有著結(jié)構(gòu)簡單、收斂性好、速度快等優(yōu)點(diǎn)，目前得到了廣泛的使用[13]。本文中，RBF神經(jīng)網(wǎng)絡(luò)的輸入為測(cè)試信號(hào)經(jīng)過圖像預(yù)處理、特征提取后的低維數(shù)據(jù)，輸出為本實(shí)驗(yàn)采用的7種常用的LPI雷達(dá)信號(hào)。

所有信號(hào)載頻由900 Hz至1 400 Hz均勻分布,采樣頻率為10 kHz。BPSK信號(hào)采用13位Barker碼；FMCW調(diào)制帶寬為250 Hz和500 Hz，調(diào)制周期為20 ms；Frank信號(hào)、P1信號(hào)、P2信號(hào)，設(shè)置脈沖壓縮比為64，編碼周期為1；COSTAS跳頻信號(hào)，跳頻序列設(shè)為{3,2,6,4,5,1}×200 Hz，{2,4,8,5,10,9,7,3,6,1}×200 Hz；多時(shí)編碼信號(hào)PT1，相位狀態(tài)數(shù)設(shè)為2，步進(jìn)段數(shù)為4。從2 dB到10 dB，每隔2 dB生成10個(gè)各種參數(shù)的信號(hào)，最終生成共280個(gè)訓(xùn)練信號(hào)。利用本文提出的LPI雷達(dá)信號(hào)分類識(shí)別算法，對(duì)這280個(gè)訓(xùn)練信號(hào)進(jìn)行處理，利用提取的特征參數(shù)訓(xùn)練并保存網(wǎng)絡(luò)。

文中進(jìn)行了2組實(shí)驗(yàn)，一組為固定信噪比情況下，測(cè)試不同信號(hào)，研究不同信號(hào)測(cè)試的識(shí)別混淆問題；另一組為固定測(cè)試信號(hào)種類，改變信噪比，研究不同信噪比情況下信號(hào)的識(shí)別正確率。

第1組實(shí)驗(yàn)，在0 dB信噪比條件下，每種信號(hào)按以上信號(hào)參數(shù)產(chǎn)生100個(gè)測(cè)試信號(hào)，利用本文方法進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果如表1所示，當(dāng)信號(hào)信噪比為0 dB情況下，識(shí)別效果良好，平均正確識(shí)別率在99%以上。

第2組實(shí)驗(yàn)，每種LPI雷達(dá)信號(hào)，信噪比從-10 dB到10 dB每隔2 dB產(chǎn)生100個(gè)測(cè)試信號(hào)，對(duì)測(cè)試信號(hào)進(jìn)行分類識(shí)別實(shí)驗(yàn)。

圖5為第2組實(shí)驗(yàn)的實(shí)驗(yàn)結(jié)果，在信噪比SNR>-2 dB時(shí)，對(duì)于本文采用的7種LPI雷達(dá)信號(hào)都能有效識(shí)別，平均正確識(shí)別率達(dá)到了93%以上，圖像相似度較大的Frank多相碼信號(hào)和P1多相碼信號(hào)也可以有效區(qū)分；當(dāng)SNR<-6 dB時(shí)，由于PWVD時(shí)頻分析方法受交叉項(xiàng)影響嚴(yán)重，對(duì)于部分LPI雷達(dá)信號(hào)已經(jīng)無法有效分析，正確識(shí)別率較低。

圖5　同一LPI信號(hào)不同信噪比情況下的識(shí)別率Fig.5　Average correct recognition rate of the　　　 same LPI radar signal in different SNR

本文對(duì)幾種LPI雷達(dá)信號(hào)時(shí)頻圖像降噪方進(jìn)行了對(duì)比。其中文獻(xiàn)[1-2]采用二值化后腐蝕膨脹的算法降噪，文獻(xiàn)[3]采用二維高斯低通濾波器降噪。實(shí)驗(yàn)過程中，3種降噪算法采用相同的仿真條件，即使用本文采用的信號(hào)，時(shí)頻分析方法和特征提取方法，所得實(shí)驗(yàn)結(jié)果的不同主要由降噪算法的不同所致。圖6表示3種方法在這種實(shí)驗(yàn)條件下的平均正確識(shí)別率。

在信噪比低于10 dB的情況下，本文方法明顯優(yōu)于其他2種方法。文獻(xiàn)[1-2]采用的二值化圖像，利用圖像腐蝕膨脹降噪，無法獲取圖像的一些細(xì)節(jié)信息，如圖2中的Frank編碼和P1編碼的時(shí)頻圖像二值化后，為斜率近似的直線，難以辨識(shí)；另一方面二值化圖像采用主分量分析法， 無法獲取最優(yōu)主元，也是識(shí)別率較低的原因。文獻(xiàn)[3]采用二維高斯低通濾波器，一方面頻率矩陣值ω1，ω2設(shè)置憑經(jīng)驗(yàn)設(shè)置，難以獲得適應(yīng)各種信噪比條件的最優(yōu)值；另一方面，采用低通濾波器，在去除高頻噪聲的同時(shí)，也使圖像喪失了相當(dāng)一部分高頻細(xì)節(jié)信息，降低信號(hào)辨識(shí)度。

圖6　不同時(shí)頻圖像降噪方法對(duì)比Fig.6　Contrast of different time-frequency　　　image de-noising methods

輸入信號(hào)BPSKFrankCOSTASFMCWP1P2PT1BPSK9811Frank994COSTAS100FMCW100P195P2100PT12100

3結(jié)束語

本文針對(duì)LPI雷達(dá)信號(hào)時(shí)頻圖像受噪聲干擾嚴(yán)重的問題，使用一種基于二維EMD快速分解的圖像降噪方法，有效地降低噪聲干擾，使不同LPI雷達(dá)信號(hào)的時(shí)頻圖像辨識(shí)度更佳。處理后的時(shí)頻圖像，經(jīng)過圖像裁剪、特征提取，利用神經(jīng)網(wǎng)絡(luò)完成了LPI雷達(dá)信號(hào)分類識(shí)別。在仿真條件相同的條件下，對(duì)比了幾種LPI時(shí)頻圖像降噪的算法，證明本文方法能夠更好地降低時(shí)頻圖像噪聲，提高LPI雷達(dá)信號(hào)分類識(shí)別的正確率。但當(dāng)信噪比低于-4 dB時(shí)，識(shí)別效果不佳，因?yàn)椴糠中盘?hào)在較低信噪比下，PWVD時(shí)頻分析方法受交叉項(xiàng)影響嚴(yán)重，無法獲得有效的時(shí)頻圖像，有必要進(jìn)一步研究對(duì)交叉項(xiàng)的抑制。

參考文獻(xiàn)：

[1]熊坤來, 羅景青, 吳世龍. 基于時(shí)頻圖像和神經(jīng)網(wǎng)絡(luò)的LPI雷達(dá)信號(hào)調(diào)制識(shí)別[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2012, 31(5): 230-233.

XIONG Kun-lai, LUO Jing-qing, WU Shi-long. Modulation Identification of LPI Radar Signals Based on Time Frequency Image and Neural Network[J]. Journal of Projectiles，Rockets，Missiles and Guidance, 2011, 31(5): 230-233.

[2]白航, 趙擁軍, 胡德秀, 等. 基于Choi-Williams時(shí)頻圖像特征的雷達(dá)輻射源識(shí)別[J]. 數(shù)據(jù)采集與處理, 2012, 27(4): 480-485.

BAI Hang, ZHAO Yong-jun, HU De-xiu, et al. Radar Emitter Recognition Based on Image Feature of Choi-Williams Time-Frequency Distribution[J]. Journal of Data Acquisition and Processing, 2012, 27(4): 480-485.

[3]ZILBERMAN E R. Autonomous Time-Frequency Cropping and Feature-Extraction Algorithms for Classification of LPI Radar Modulations[R]. Naval Postgraduate School Monterey CA, 2006.

[4]白航, 趙擁軍, 胡德秀, 等. 時(shí)頻圖像局部二值模式特征在雷達(dá)信號(hào)分類識(shí)別中的應(yīng)用[J]. 宇航學(xué)報(bào), 2013, 34(1): 139-146.

BAI Hang, ZHAO Yong-jun, HU De-xiu, et al. Radar Signal Recognition Based on the Local Binary Pattern Feature of Time-Frequency Image[J]. Journal of Astronautics, 2013, 34(1): 139-146.

[5]郝士琦, 周建國, 冷蛟鋒, 等. 基于時(shí)頻圖像代數(shù)特征的雷達(dá)信號(hào)脈內(nèi)調(diào)制方式分類研究[J]. 電子對(duì)抗, 2012(4):22-26.

HAO Shi-qi, ZHOU Jian-guo, LENG Jiao-feng, et al. Research on Classifications of Intra-Pulse Modulation Based on Algebra Characteristics of TFD Images [J]. Electronic Warfare, 2012(4):22-26.

[6]冀貞海, 朱偉強(qiáng), 趙力, 等. 二維主分量分析的脈內(nèi)調(diào)制識(shí)別算法研究[J]. 電光與控制, 2009, 16(11): 33-37.

JI Zhen-hai, ZHU Wei-qiang, ZHAO Li, et al. Intra-Pulse Modulation Recognition Method Based on Two-Dimensional Principal Component Analysis[J]. Electronics Optics and Control, 2009, 16(11): 33-37.

[7]NUNES J C, BOUAOUNE Y, DELECHELLE E, et al. Image Analysis by Bidimensional Empirical Mode Decomposition [J]. Image and Vision Computing, 2003, 21(12):1019-1026.

[8]戴桂平. 基于二維EMD和小波閾值的掌紋圖像去噪[J]. 計(jì)量學(xué)報(bào), 2011, 32(4): 368-372.

DAI Gui-ping. Palm Print Image De-noising Based on 2-D EMD and Wavelet Thresholding[J]. Acta Metrologica Sinica, 2011, 32(4): 368-372.

[9]趙春暉, 任龍濤, 萬建. 一種基于二維EMD的圖像融合方法[J]. 應(yīng)用科技, 2009, 36(9): 15-19.

ZHAO Chun-hui, REN Long-tao, WAN Jian. A Method of Image Fusion Based on the Two-Dimensional Empirical Mode Decomposition[J]. Applied Science and Techology, 2009, 36(9): 15-19.

[10]郭珈, 王孝通, 徐曉剛, 等. 基于均值濾波矩陣組的二維信號(hào)快速經(jīng)驗(yàn)?zāi)Ｊ椒纸夥椒╗J]. 電子學(xué)報(bào), 2013, 40(10): 1980-1983.

GUO Jia, WANG Xiao-tong, XU Xiao-gang, et al. Fast 2D Empirical Mode Decomposition Based on Mean Value Matrix Filter Sets[J]. Acta Electronica Sinica, 2012, 40(10): 1980-1983.

[11]郭聳, 顧國昌, 李常有, 等. 利用 EMD 的自適應(yīng)圖像去噪[J]. 計(jì)算機(jī)工程與應(yīng)用, 2013, 49(8): 12-16.

GUO Song, GU Guo-chang, LI Chang-you, et al. Adaptive Image Denoising Based on EMD[J].Computer Engineering and Applications, 2013, 49(8):12-16.

[12]FARGUES M P. Investigation of Feature Dimension Reduction Schemes for Classification Applications[R].Naval Postgraduate School Monterey CA Dept of Electrical and Computer Engineering, 2001.

[13]孫志強(qiáng)，葛哲學(xué)．神經(jīng)網(wǎng)絡(luò)理論與MATLAB 7實(shí)現(xiàn)[M]．北京：電子工業(yè)出版社，2005：79-81．

SUN Zhi-qiang, GE Zhe-xue. Neural Network Theory and the Realization of MATLAB7[M]. Beijing: Electronic Industry Press, 2005：79-81.

Radar Signal Classification Based on Two-Dimensional Empirical Mode Decomposition of Time-Frequency Images

ZHANG Fan1, DIAO Ming1, YANG Cheng-zhi1，2

(1. Harbin Engineering University, College of Communication and Information Technology,Heilongjiang Harbin 150001, China 2. Aviation University of Air Force, Department of Communication Counter,Jilin Changchun 130022, China)

Abstract:To classify low probability of intercept (LPI) radar signals correctly, a novel method based on two-dimensional fast empirical mode decomposition (FBEMD) algorithm is proposed. The method is used to reduce the time-frequency image noise. Firstly, time-frequency images are obtained by using PWVD distribution. Intercept image part with signal information and extract feature vector by principal component analysis (PCA). Finally, the radial basis function(RBF) neural network is used to automatically classify LPI radar signals. Simulation results show that the method is effective at low signal noise ratio(SNR). By using FBEMD noise reduction algorithm, the average correct recognition rate can reach 93% when SNR is -2 dB.

Key words:low probability of intercept(LPI) radar signals; time-frequency distribution; two-dimensional empirical mode decomposition(EMD); principal component analysis

中圖分類號(hào)：TN957.52；TP391.9

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1009-086X(2015)-05-0172-06

doi:10.3969/j.issn.1009-086x.2015.05.028

通信地址：150001黑龍江省哈爾濱市南崗區(qū)哈爾濱工程大學(xué)信息與通信工程系E-mail：zhangfan8908@gmail.com

作者簡介：張帆(1989-)，男，山東日照人。碩士生，研究方向?yàn)槔走_(dá)信號(hào)處理。

*收稿日期：2014-06-05；修回日期：2014-08-28