基于圖像特征的雷達(dá)脈內(nèi)調(diào)制類(lèi)型識(shí)別*

高君豐 劉芬芬

（1.中國(guó)人民解放軍92785部隊(duì) 秦皇島 066200）（2.中國(guó)人民解放軍91404部隊(duì) 秦皇島 066000）

1 引言

由于軍事技術(shù)的不斷進(jìn)步和發(fā)展，新體制及復(fù)雜體制雷達(dá)不斷被投入使用，人們接收到的輻射源信號(hào)脈內(nèi)調(diào)制形式復(fù)雜、參數(shù)多變，給識(shí)別帶來(lái)了極大的挑戰(zhàn)。針對(duì)復(fù)雜體制雷達(dá)信號(hào)的脈內(nèi)調(diào)制類(lèi)型識(shí)別已經(jīng)成為電子對(duì)抗領(lǐng)域中亟待解決的關(guān)鍵問(wèn)題。傳統(tǒng)的識(shí)別方法包括時(shí)域自相關(guān)法［1］、相位差分法［2～4］、小波變換［5～6］等，但是這些算法抗噪性能較差，隨著信噪比的降低，識(shí)別性能會(huì)急劇下降。由于不同信號(hào)之間的時(shí)頻分布圖差異性很大，從時(shí)頻圖上可以很容易地判別出信號(hào)的調(diào)制類(lèi)型，所以對(duì)雷達(dá)輻射源信號(hào)進(jìn)行時(shí)頻變換，提取時(shí)頻特征對(duì)信號(hào)分類(lèi)識(shí)別是近年來(lái)的研究熱點(diǎn)［7～8］。

本文利用圖像處理方法對(duì)信號(hào)時(shí)頻圖進(jìn)行預(yù)處理，包括圖像濾波、圖像剪切、圖像二值化等，去掉背景噪聲等冗余信息，截取出圖像中信號(hào)所在部分；然后提取出圖像的中心矩和偽Zernike矩特征，再利用隨機(jī)森林構(gòu)建分類(lèi)器，識(shí)別信號(hào)的調(diào)制樣式。

2 時(shí)頻變換

時(shí)頻分析作為分析和處理非平穩(wěn)信號(hào)的有力工具，反映了信號(hào)能量隨時(shí)間和頻率的分布，可以在時(shí)頻域上更精確地描述信號(hào)［9～10］。最常用的時(shí)頻分布是魏格納－維爾分布（Wigner-Ville，WVD），具有很高的時(shí)頻聚集性，但多分量信號(hào)的WVD變換存在嚴(yán)重的交叉項(xiàng)，造成一些虛假信息的干擾，從而給信號(hào)的時(shí)頻分析帶來(lái)困難。學(xué)者們提出了多種抑制WVD交叉項(xiàng)的方法，有平滑Wigner-Ville分布（SWVD），偽Wigner-Ville分布（PWVD），平滑偽Wigner-Ville分布（SPWVD）。本文所用的SP?WVD采用同時(shí)對(duì)信號(hào)在時(shí)域和頻域加窗平滑的方法來(lái)抑制交叉項(xiàng)。信號(hào)x（t）的SPWVD分布定義為

式中g(shù)（u）和h（τ）分別是時(shí)域和頻域的平滑窗函數(shù)。

本文對(duì)信號(hào)的時(shí)頻變換圖進(jìn)行了歸一化操作，并將信號(hào)能量換算為以dB為單位。即時(shí)頻圖的y軸大小為

圖1給出了信噪比為10dB時(shí)幾種常見(jiàn)脈內(nèi)信號(hào)的時(shí)頻圖，從圖中可以看出，不同信號(hào)的時(shí)頻變換圖差異明顯，肉眼很容易區(qū)分，但要實(shí)現(xiàn)自動(dòng)識(shí)別還要結(jié)合數(shù)字圖像處理技術(shù)進(jìn)一步提取信號(hào)的時(shí)頻特征。

圖1 8類(lèi)圖像的SPWVD時(shí)頻圖

3 圖像預(yù)處理

由于時(shí)頻圖像受噪聲影響較大，并且存在大量的冗余信息，所以在提取特征之前，需要將時(shí)頻圖像進(jìn)行預(yù)處理，以提升圖像的質(zhì)量，保證所提取信號(hào)特征的有效性。

圖2 圖像預(yù)處理過(guò)程

首先將彩色的時(shí)頻圖轉(zhuǎn)換為灰度圖，通過(guò)拉普拉斯銳化［11］，增強(qiáng)圖像的邊緣和灰度跳變的細(xì)節(jié)部分，提高圖像的對(duì)比度，從而使圖像更加清晰。通過(guò)自適應(yīng)維納濾波抑制目標(biāo)信號(hào)和背景中的噪聲，達(dá)到圖像增強(qiáng)的效果［12］。為了減少數(shù)據(jù)間的不均衡性，對(duì)時(shí)頻圖像灰度值進(jìn)行歸一化處理，歸一化后像素點(diǎn)的灰度值為x=(x-)/σ，和σ2分別為所有像素點(diǎn)灰度值的均值和方差。從圖1中可以看出，時(shí)頻圖中大多都是背景，并不是所有的頻率點(diǎn)都存在信號(hào)，沿頻率軸檢測(cè)出信號(hào)時(shí)頻分布的起止頻率，將沒(méi)有信號(hào)分布的圖像區(qū)域剪切掉，能夠減小冗余信息，有利于下一步信號(hào)特征的提取。

經(jīng)過(guò)上述一系列圖像預(yù)處理后，通過(guò)灰度閾值化處理將時(shí)頻圖像轉(zhuǎn)化為二值圖像，二值化使圖像變得簡(jiǎn)單，且數(shù)據(jù)量大為減少，能凸顯出目標(biāo)的輪廓。圖像的二值化處理可以描述為

式中，P(t，f)表示信號(hào)的時(shí)頻圖，二值化過(guò)程中閾值th的選取采用一維最大熵法［13］。

經(jīng)過(guò)預(yù)處理后的時(shí)頻圖像形狀特征更加明顯，如圖3所示。

圖3 預(yù)處理后的時(shí)頻圖

4 特征提取

由圖1和圖3的時(shí)頻圖像可以看出，不同調(diào)制類(lèi)型信號(hào)的幾何形狀具有比較明顯的差異，因此可以提取時(shí)頻圖像的形狀特征用于信號(hào)的分類(lèi)識(shí)別。矩函數(shù)是一種有效的形狀特征描述子，能夠提供大量與圖像有關(guān)的信息，如圖像的形狀、大小和方向等，本文主要選取了矩函數(shù)中的歸一化中心矩以及偽Zernike矩。

4.1 歸一化中心矩

設(shè)大小為M×N的圖像表示為f(x，y)，則圖像的(p+q)階原點(diǎn)矩定義為

(p+q)階中心矩定義為

(p+q)階歸一化中心矩定義為

其中，u00表示圖像面積，本文選用u11、u02、u12、u21、u03這五個(gè)中心矩特征，分別表示圖像的傾斜程度、垂直方向上的伸展度、垂直方向上的均衡程度、水平方向上的均衡程度、垂直方向上的重心偏移度。

4.2 偽Zernike矩

偽Zernike矩是一種正交復(fù)數(shù)矩，對(duì)于連續(xù)函數(shù)f(x，y)，階數(shù)為p，重復(fù)度為q的偽Zernike矩定義為

式中，p為非負(fù)整數(shù)，q為整數(shù)，且滿(mǎn)足|q|≤p。*表示復(fù)數(shù)共軛，Vpq為偽Zernike多項(xiàng)式，表示為

其中，實(shí)值徑向多項(xiàng)式Rpq(ρ)用下式表示：

實(shí)際問(wèn)題中處理的圖像是數(shù)字圖像，偽Zernike矩的離散形式為

式中，I(i，j)是像素點(diǎn) (i，j)處的像素。

為了減少偽Zernike矩的動(dòng)態(tài)范圍，本文對(duì)結(jié)果取對(duì)數(shù)，即pq=ln|Zpq|。并最終采用這7個(gè)偽 Zernike矩特征作為信號(hào)的識(shí)別特征。

綜上，本文共選用12個(gè)形狀特征，作為分類(lèi)識(shí)別的特征。

5 仿真實(shí)驗(yàn)

為了驗(yàn)證本文算法的有效性，利用Matlab仿真了8種雷達(dá)脈內(nèi)調(diào)制信號(hào)，分別是常規(guī)信號(hào)（CW）、線(xiàn)性調(diào)頻信號(hào)（LFM）、非線(xiàn)性調(diào)頻信號(hào)（NLFM）、二相編碼信號(hào)（BPSK）、四相編碼信號(hào)（QPSK）、頻率編碼信號(hào)（FSK）以及兩種復(fù)合信號(hào)（FSK+BPSK復(fù)合、LFM+BPSK復(fù)合）。信號(hào)載頻為400MHz、采樣頻率為1000MHz、脈寬10μs。線(xiàn)性調(diào)頻信號(hào)的調(diào)頻斜率為-8MHz/μs，二相編碼和四相編碼均為32位隨機(jī)碼。非線(xiàn)性調(diào)頻信號(hào)為三階多項(xiàng)式相位信號(hào)，多項(xiàng)式系數(shù)分別為15和-4。頻率編碼是2FSK信號(hào)，頻點(diǎn)值分別是300MHz和330MHz，跳頻序列是［1 0 1 0 1 1］，其中1代表頻點(diǎn)330MHz，0代表300MHz。FSK+BPSK復(fù)合信號(hào)中FSK信號(hào)碼元為300MHz和400MHz，跳頻序列是［1 0 1 0］，其中1代表頻點(diǎn)400MHz，0代表300MHz，BPSK為7位巴克碼［1 1 1 0 0 1 0］。LFM+BPSK復(fù)合信號(hào)中LFM的調(diào)頻率為-8MHz/μs，BPSK為13位巴克碼［1 1 1 1 1 0 0 1 1 0 1 0 1］。

將信噪比從-8dB變化到12dB，每一類(lèi)信號(hào)隔2dB生成300個(gè)樣本，共生成26400個(gè)樣本。其中，每一類(lèi)信號(hào)每個(gè)信噪比下200個(gè)樣本作為訓(xùn)練集，共計(jì)17600個(gè)樣本；每一類(lèi)信號(hào)每個(gè)信噪比下剩余的100個(gè)樣本作為測(cè)試集，共計(jì)8800個(gè)樣本。經(jīng)過(guò)圖像預(yù)處理，提取特征后，送入RF分類(lèi)器，得到這8種信號(hào)的混淆矩陣如表1所示。

表1 混淆矩陣

每一種信號(hào)共有1100個(gè)樣本，從表1可以看出，CW信號(hào)和相位編碼信號(hào)容易相互識(shí)別錯(cuò)誤，這是因?yàn)樵谛旁氡容^低時(shí)，三種信號(hào)的時(shí)頻特征比較相像。這從第二節(jié)的時(shí)頻圖也可以看出，當(dāng)信噪比較低時(shí)，受噪聲的影響，二相編碼和四相編碼時(shí)頻圖上的跳變被湮沒(méi)，容易識(shí)別錯(cuò)誤；而常規(guī)信號(hào)由噪聲引起錯(cuò)誤的跳變，容易誤識(shí)別為相位編碼信號(hào)。除了這三種信號(hào)之外的其余信號(hào)的時(shí)頻特征比較明顯且不相似，識(shí)別結(jié)果均比較好。

為了更加清楚地說(shuō)明本文算法的有效性，圖4給出了每種信號(hào)在不同信噪下的正確識(shí)別率。

圖4 正確識(shí)別率

從圖4可以看出，在信噪比低至-6dB時(shí)，LFM、NLFM、FSK、FSK+BPSK、LFM+BPSK信號(hào)仍能達(dá)到98%以上的識(shí)別正確率；NS信號(hào)在信噪比大于-2dB時(shí)，識(shí)別正確率達(dá)到90%以上；BPSK、QPSK信號(hào)在信噪比不小于-2dB時(shí)，識(shí)別正確率大于85%；且隨著信噪比的增加，所有信號(hào)的識(shí)別正確率大致也是上升的，BPSK和QPSK由于本身時(shí)頻特征比較相像，且跳變特征沒(méi)有很顯著，導(dǎo)致有輕微的波動(dòng)，但大體上識(shí)別正確率曲線(xiàn)還是隨著信噪比上升的。

為了測(cè)試隨機(jī)森林分類(lèi)器的性能，本文分別利用隨機(jī)森林、支持向量機(jī)（SVM）、K近鄰（KNN）三種分類(lèi)方法對(duì)時(shí)頻特征不明顯的相位編碼信號(hào)進(jìn)行分類(lèi)識(shí)別。信號(hào)樣本同上，在不同信噪比下的識(shí)別率如圖5所示，可以看出，隨機(jī)森林的識(shí)別準(zhǔn)確率最高，明顯優(yōu)于SVM和KNN。

圖5 不同分類(lèi)器的分類(lèi)結(jié)果

6 結(jié)語(yǔ)

本文針對(duì)雷達(dá)中復(fù)雜調(diào)制類(lèi)型的識(shí)別問(wèn)題，運(yùn)用數(shù)字圖像處理技術(shù)和形狀特征理論，提出了一種對(duì)時(shí)頻圖像預(yù)處理后提取圖像矩特征，再利用隨機(jī)森林分類(lèi)器識(shí)別的方法。對(duì)8種常見(jiàn)信號(hào)的仿真結(jié)果表明，本文算法在低信噪比下也能獲得很高的識(shí)別準(zhǔn)確率，具有較強(qiáng)的抗噪能力。并且對(duì)于目前比較復(fù)雜的復(fù)合調(diào)制類(lèi)型，識(shí)別率也很高，具有一定的工程應(yīng)用價(jià)值。