雷达辐射源个体识别是通过对不同辐射源信号所携带的无意调制的识别来实现。这些无意调制也称为辐射源的指纹,它是由辐射源自身的硬件缺陷所导致,具有较高的稳定性,且能够唯一地表征雷达辐射源个体[1]。
在当今电子战领域中精确有效的辐射源个体识别是把握战争主动权的关键所在,这也使得它成为当今学术界的热门研究问题。现有的识别方法主要分为2类:1) 基于原始信号域的方法,例如提取信号的瞬时振幅、频率、相位[2-3]以及信号的脉宽、带宽、重复周期等脉冲描述子特征[4-5]进行识别。虽然这类方法计算简单,但抗噪声性能差,并且随着雷达信号调制技术的复杂化,常规参数特征的获取也变得更加困难;2) 基于变换域的方法,将接收到的信号转换为新的信号域进行特征提取,如小波变换[6-7]、Gabor展开、希尔伯特−黄变换[8]、双谱分析[9]等。然而,基于小波的方法很大程度上依赖于对基本小波函数的选择;基于双谱和希尔伯特−黄变换的方法都存在特征维数较高、计算量大的问题。并且这2类方法都需要人工设计提取的特征,大大增加了特征提取的难度。
近些年来,随着深度学习的快速发展,使其在识别领域得到了成功的应用,如利用深度学习进行信号调制样式识别[10]、手写汉字识别[11]、人脸识别等。该类方法能够自动地提取待识别个体的深层次特征,这不仅降低了特征提取的难度、减少了计算量,而且还提高了低信噪比下的识别准确率。因此本文基于信号时频分析并引用深度学习中的变分自编码器,提出了一种新形式的雷达辐射源个体识别算法,该算法将信号域中的辐射源个体识别问题转变到图像域,摒弃了人工提取特征的繁琐,利用变分编码器自动获取图像特征并进行分类识别,仿真结果表明,相较于其他算法该算法降低了特征提取的难度,提高了低信噪比下的识别率。
1 识别系统框架雷达辐射源个体识别系统框架如图1所示。
Download:
|
|
首先将获取的雷达辐射源信号进行预处理:利用时频分析将信号域的一维信号转换为图像域的二维时频图,并对图像进行灰度化、归一化、尺寸调整等操作,以减少图像数据量并突出图像的重要信息。然后利用变分自编码器对预处理后的图像进行特征提取,以获取图像的深层次特征。接着为了降低识别系统的复杂度采用核PCA提取特征中的主成分,最后引用支持向量机作为分类器,将处理后的特征送入支持向量机中识别出不同的辐射源个体。
2 信号预处理 2.1 时频分析雷达辐射源信号作为一种非平稳的信号,其统计量是时变函数,这时只了解信号的时域或频域的全局特性是远远不够的[12]。时频分析反映了信号时域和频域的联合分布信息,而时频图像清晰地描述了信号频率随时间变换的关系,因此对信号进行时频分析是研究非平稳信号的重要方法。
Choi-Williams时频分布由于其良好的时频聚集性、较少的交叉项和时域、频域中较高的识别精度,成为非线性时频分析中常用的方法[13]。考虑到不同辐射源所携带的无意调制可能随机分布在信号的各个频率分量上,因此本文选取Choi-Williams时频分布来获取辐射源信号的二维时频图像, Choi-Williams时频分布为
$\begin{split} {C_x}(t,\omega ) = & \iint\limits_\infty {\sqrt {\frac{\sigma }{{4 {\text{π}} {\tau ^2}}}} \exp \left[ { - \frac{{{{(t - \mu )}^2}}}{{4{\tau ^2}/\sigma }}} \right]x\left( {\mu + \frac{\tau }{2}} \right)} \\ & {x^*}\left( {\mu - \frac{\tau }{2}} \right){{\rm{e}}^{ - j\omega \tau }}{\rm{d}}\mu {\rm{d}}\tau \end{split} $ |
式中:
定义核函数为
$f(\tau ,\mu ) = \exp \left[ { - \frac{{{\mu ^2}}}{{{{4{\tau ^2}} / \sigma }}}} \right]$ |
图2给出了4部辐射源所发射的线性调频信号的时频图像。其中信号的采样频率为80 MHz,初始频率为10 MHz,带宽20 MHz,脉冲宽度5 μs,信噪比0 dB。由图2可以看出不同辐射源的信号差异较小,仅在图像边缘略有不同,而面对这种情况常规的特征提取方法很难奏效,这就需要对图像深层次的特征进行挖掘,这也是本文采用深度学习进行特征提取的主要原因。
Download:
|
|
如果将原始时频图像直接送入变分自编码网络中进行特征提取,图像维度较高,数据间不平衡将会降低图像特征提取效率,因此本文采用图3所示流程对时频图像进行预处理。
Download:
|
|
首先为了降低图片的数据量,本文采用加权平均法将获取的RGB时频图像进行灰度化。
定义灰度级范围为[0, 255],图片中每个像素点的灰度值为
${\rm{Gray}}\left( {i,j} \right) = 0.299B\left( {i,j} \right) + 0.578G\left( {i,j} \right) + 0.114R\left( {i,j} \right)$ |
为了减小图像几何变换的影响加快训练网络的收敛性,本文对灰度图像中的像素点进行归一化处理:
$x = {{\left( {{x_{i,j}} - \bar x} \right)} \Bigg/ {\sqrt {\frac{1}{{n - 1}}\displaystyle \sum\limits_{i = 1}^n {} \displaystyle \sum\limits_{j = 1}^n {({x_{i,j}}} - \bar x{)^2}} }}$ |
式中:
最后为了适应神经网络的输入要求,提高特征提取的速度,需要对归一化后的时频图像进行尺寸调整。本节拟采用局部均值法对归一化图像进行缩小。
设原图像为
$ \begin{array}{l} {{{F}}'}\left( {x,y} \right) = \left[ {\begin{array}{*{20}{c}} {{F_{x - 1,y - 1}}}& \cdots &{{F_{x - 1,y}}} \\ \vdots & & \vdots \\ {{F_{x,y - 1}}}& \cdots &{{F_{x,y}}} \end{array}} \right]\\ \;\;\;\;\;\;{F_{x - 1,y - 1}}=F\left( {i \left( {x - 1} \right) + 1,j \left( {y - 1} \right) + 1} \right)\\ \;\;\;\;\;\;\;\;\;\;\;\;{F_{x - 1,y}}=F\left( {i \left( {x - 1} \right) + 1,j y} \right)\\ \;\;\;\;\;\;\;\;\;\;\;\;{F_{x,y - 1}}=F\left( {i x,j \left( {y - 1} \right) + 1} \right)\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{F_{x,y}} = F\left( {i x,j y} \right) \end{array} $ |
式中:
$g\left( {{x'},{y'}} \right) = {\rm{mean}}\left( {{{{F}}'}\left( {x,y} \right)} \right)$ |
式中
如图4所示,以辐射源1的时频图像为例给出了图像预处理各阶段的处理效果。
Download:
|
|
变分自编码器(variational autoencoder,VAE)是一种将神经网络与贝叶斯概率图相结合的深度生成网络,其具有非常强大的自动学习能力,它可以从一组无标注的数据中学习到该组数据的深层特征,然后用该特征反向表示原始数据,其网络结构如图5所示。
传统自编码器通过编码网络将输入数据转换成一个编码向量,向量的每个维度表示学到的数据的属性,并且每个编码维度均为单值。而VAE以概率的方式描述观察空间,其编码网络的输出是描述潜在空间中每个维度分布的参数[14-16]。
Download:
|
|
假设输入数据为
VAE可被分为识别模型和生成模型2个过程。为了估计出真实后验分布
${\rm{KL}}\left( {{q_\varphi }\left( {{{Z}}/{{X}}} \right)||{Q_\theta }\left( {{{Z}}/{{X}}} \right)} \right) = {E_{{q_\varphi }\left( {{{Z}}/{{X}}} \right)}}\log \frac{{{q_\varphi }\left( {{{Z}}/{{X}}} \right)}}{{{Q_\theta }\left( {{{Z}}/{{X}}} \right)}}$ |
并根据相似程度调整识别模型的网络参数。生成模型则是由
将预处理后的时频图像送入VAE网络中,提取隐层变量
1)在预处理后的时频图像集中,随机选取50%的图像作为训练数据,剩余的50%作为测试数据。
2)训练特征提取模型,将训练集中的图像送到VAE网络中训练网络,得到隐含层的初始化网络参数,并根据网络的损失函数反向调节参数。
3)提取测试特征,将测试集图像送到已经训练好的网络中,提取隐层变量
变分自编码器能够提取到图片中的各类属性特征,但各特征间难免存在冗余,并且随着调制样式的复杂,特征集维度也随之增加。如果不进行特征降维,后续识别过程将会陷入维度灾难[17]。
4.1 核主成分分析核主成分分析(kernel principal component analysis, KPCA)是为了克服传统主成分分析不能有效处理非线性数据的缺点,从而借鉴核函数的思想,采用非线性映射核函数
$\lambda{{ v}} = {{Cv}}$ | (1) |
式中:
${{C}} = \frac{1}{M}\displaystyle \sum\limits_{j = 1}^M {{{\phi}} \left( {{x_j}} \right)} {{\phi}} {\left( {{x_j}} \right)^{\rm{T}}}$ | (2) |
进而有
$\lambda \left( {{{\phi}} \left( {{x_k}} \right){{v}}} \right) = {{\phi}} \left( {{x_k}} \right){{C}}{{{v}}^{}}, {k = 1,2, \cdots ,{{M}}} $ | (3) |
式中
${{v}} = \displaystyle \sum\limits_{i = 1}^M {{a_i}{{\phi}} \left( {{x_i}} \right)} $ | (4) |
将式(2)、(4)代入式(3)得到:
$\begin{split} &\lambda \displaystyle \sum\limits_{i = 1}^M {{a_i}} \left( {{{\phi}} \left( {{x_k}} \right){{\phi}} \left( {{x_i}} \right)} \right) = \\ &\frac{1}{M} \displaystyle \sum\limits_{i = 1}^M {{a_i}}\left( {{{\phi}} \left( {{x_k}} \right) \displaystyle \sum\limits_{j = 1}^M {{{\phi}} \left( {{x_j}} \right)} } \right)\left( {{{\phi}} \left( {{x_j}} \right){{\phi}} \left( {{x_i}} \right)} \right) \end{split}$ |
式中
$M\lambda {{a}} = K{{a}}$ | (5) |
通过式(5)可以求得特征值
鉴于VAE提取的特征中的数据可能为非线性分布,因此本文应用KPCA对其进行主成分分析,以降低数据的维度。具体步骤如下:
1)选定核函数(本文选用高斯径向核函数),并由
2)将
$\eta \left( n \right) =\displaystyle \sum\limits_{m = 1}^g {{\lambda _k}/} \displaystyle \sum\limits_{m = 1}^r {{\lambda _m}} $ |
3)当
${{{V}}'} = {\left[ {{v_1},{v_2}, \cdots ,{v_n}} \right]^{}},\;\;n < r$ |
4)计算降维后的特征
${{Z}} = {\left( {{{{V}}'}} \right)^{\rm{T}}}{{W}}$ |
设置VAE生成模型输入节点为200,则提取的图像特征为200维。以辐射源1的特征为例,表1给出了不同信噪比下利用KPCA进行特征降维的结果。由表1可以看出以
支持向量机(support vector machine,SVM)作为一种常用的分类器通过非线性映射将输入数据映射到高维空间,然后在高维空间中寻找最优分类超平面从而实现数据的分类问题[20]。因此支持向量分类的关键在于最优分类超平面的求解问题。
设样本点集为
${y_i}\left[ {{{w}} {{{x}}_i} + b} \right] \geqslant {1} \;\;\;i = 1,2, \cdots, n$ |
下,求解最小的
$\varGamma = \dfrac{1}{2}||{{w}} |{|^2}$ | (6) |
为了求解上述问题引入拉格朗日乘子
$L\left( {{{w}},b,\alpha } \right) = \frac{1}{2}||{{w}}|{|^2} -\displaystyle \sum\limits_{i = 1}^n {{a_i}\left( {{y_i}\left( {{{{x}}_i}{{w}} + b} \right) - 1} \right){,^{}}{a_i} \geqslant 0} $ |
然后将拉格朗日函数分别关于
$L\left( {{{w}},b,\alpha } \right) = \sum\limits_{i = 1}^n {{a_i} - \frac{1}{2}} \sum\limits_{i,j = 1}^n {{a_i}{a_j}{y_i}{y_j}{{x}}_i^{\rm{T}}{{{x}}_j}} $ |
通过求解
$f\left( {{x}} \right) = {\rm{sgn}} \left[ {\displaystyle \sum\limits_{i = 1}^n {{y_j}{a_i}\left( {{{x}}{{{x}}_i}} \right)} + b} \right]$ |
但在实际情况中分类数据一般不满足线性可分,尤其是对于雷达辐射的特征数据集,因此针对非线性数据引入松弛变量
${y_i}\left[ {{{w}} {{{x}}_i} + b} \right] \geqslant 1 - {\zeta _i},\;\;i = 1,2 ,\cdots ,l$ |
最小间距
$\varGamma = \frac{1}{2}||{{w}} |{|^2} + {{C}}{\displaystyle \sum\limits_{i = 1}^l {{\zeta _i}} _{}}{,_{}}\zeta \geqslant 0$ |
由此计算得到最终的决策函数表示为
$f\left( {{x}} \right) = {\rm{sgn}} \left[ {\displaystyle \sum\limits_{i = 1}^n {{\alpha _i}{y_j}K\left( {{{x}},{{{x}}_i}} \right)} + b} \right]$ |
式中
利用SVM进行分类的步骤如下:
1)将4.2节获得降维后的特征分为训练集和测试集2部分。
2)用训练集特征训练SVM,获得最优的参数。
3)用测试集进行测试,记录识别结果。
6 仿真与分析 6.1 仿真条件为了模拟真实雷达辐射源信号,本文将采用6部信号源进行仿真分析,6部信号源分别为Agilent N5172B、Tektronix AWG70001、Agilent E4438C以及同一型号的Agilent N5172B 2部。并分别以辐射源1~6命名。信号包括常规脉冲信号(continuous wave, CW)、线性调频信号(linear frequency modulation, LFM)、二进制编码压缩信号(binary phase shift keying, BPSK)。其中CW和LFM的载频为500 M、1 G、2 G,带宽为10 M、20 M、30 M。BPSK使用13位巴克码。所采信号的脉宽均为10 μs,占空比为50%。模拟出信噪比为−8 、−4 、0 、4 、8 、12 、16 、20 、24 dB 9种识别环境,并且在每种信噪比、每种信号参数下截取信号脉冲1 000个。在进行识别时信号集中50%的信号用于训练,剩余的50%用于测试。
6.2 实验结果与分析实验1 本文对型号相同但个体不同的辐射源的识别进行仿真。以LFM信号为例,将2部型号相同的辐射源产生的LFM信号送入识别框架中进行识别。不同信噪比下的识别结果如图6所示。
Download:
|
|
仿真结果表明,本文所提算法能够将型号相同但个体不同的辐射源区分开来。当信噪比大于0 dB时识别率均高于93%,并且随着信噪比的升高,识别率也逐渐上升,当信噪比大于12 dB时趋于收敛。从而说明了虽然辐射源型号相同但不同个体间仍存在差异,并进一步说明了本文所提算法对同型号辐射源识别的适用性。
实验2 为了研究信号调制样式对识别结果的影响,将6部信号源的不同调制样式的信号分别送入到识别框架中进行识别。信噪比在−8~ 24 dB范围内的识别效果如图7所示。
Download:
|
|
仿真结果表明,在相同信噪比环境下,调制样式对识别效果影响较小,尤其是当信噪比大于12 dB时,识别率几乎没有差别,均达到99%以上。这说明了雷达辐射源无意调制与信号调制样式无关,也证明了本文算法的普遍适用性。并且对比图6、7可知,当信噪比低于0 dB时个体数量对识别效果有一定影响,但当信噪比高于0 dB时识别效果几乎没有差别,原因在于低信噪比环境下信号时频图像较为模糊,从而影响了VAE对图像特征的提取。
实验3 为了比较本文所提算法与传统算法的抗噪声性能差异,选取信号域中的瞬时相位(instantaneous phase, IP)、瞬时频率(instantaneous frequency, IF),变换域中的双谱特征、小波特征特征与本文算法对比,6部辐射源不同调制信号下的平均识别率如图8所示。
Download:
|
|
由图8可知,当信噪比小于8 dB时,随着信噪比的降低,5种方法的识别率均存在大幅度的下降,但在相同信噪比下本文所提算法的识别率仍要高于其他4种方法。并且当SNR大于12 dB时本文所提算法的识别率逐渐收敛,而基于变换域和基于信号域的方法分别在SNR为20 dB和24 dB时才趋于收敛。由此可证明本文算法抗噪声性能的优越性和较高的系统稳定性。
实验4 将本文算法与文献[7]中基于小波变换的识别方法和文献[9]中基于双谱分析的识别方法进行特征提取复杂度和系统识别复杂度仿真。测试以LFM信号为例,信噪比范围为−4~4 dB,分别提取小波特征、双谱特征和本文时频图像特征送入支持向量机中识别并记录特征提取的耗时,具体的测试环境如表2所示。不同信噪比下的特征提取复杂度仿真结果如表3所示,系统识别复杂度仿真结果如表4所示。
仿真结果表明,相同信噪比下,本文所提算法在特征提取耗时和识别耗时上均为最小,并且随着信噪比的提升耗时也在逐步降低。原因在于本文利用离线学习的VAE网络直接获取图像特征,免去了手动处理特征的繁琐步骤,此外由于采用KPCA对特征进行了降维处理,将高维度的特征数据集转换为低维度的数据集,以此降低了支持向量机的识别分类时间。因此与其他人工特征提取识别算法相比,本文所提算法不仅性能稳定而且实时性高,具有较高的实际应用价值。
7 结论本文将辐射源的个体识别转化为图像的识别,提出了一种基于变分自编码器的雷达辐射源个体识别算法。
1)相较于传统信号域和变换域中提取的特征,本文应用变分自编码器提取的时频图像特征的抗噪声性能更强,对个体更具表征能力。
2)应用KPCA对特征进行降维,将高维度大数据个体特征转化为低维度小样本特征,提高了识别速度,增加了系统识别的实时性。
3)应用变分自编码器自动提取辐射源个体特征,克服了手动特征提取方法的繁琐,不需要对特征进行选择,降低了特征提取的难度,使得识别系统在较低信噪比下仍能取得较高的识别率。
但由于低信噪比下信号时频图像较为模糊,所提取特征分辨力不强从而影响了识别效果,下一步将在图像处理方面进行更深入的研究,提高低信噪比下时频图像的清晰度,进一步提升低信噪比下的识别性能。
[1] | RU Xiaohu, LIU Zheng, HUANG Zhitao, et al. Evaluation of unintentional modulation for pulse compression signals based on spectrum asymmetry[J]. IET radar, sonar and navigation, 2017, 11(4): 656-663. DOI:10.1049/iet-rsn.2016.0248 (0) |
[2] | CONNING M, POTGIETER F. Analysis of measured radar data for specific emitter identification[C]//2010 IEEE International Radar Conference. Washington, USA, 2010: 611–617. (0) |
[3] | RU Xiaohu, LIU Zheng, JIANG Wenli, et al. Recognition performance analysis of instantaneous phase and its transformed features for radar emitter identification[J]. IET radar sonar navigation, 2015, 10(5): 945-952. (0) |
[4] | IGLESIAS V, GRAJAL J, YESTE-OJEDA O, et al. Real-time radar pulse parameter extractor[C]// 2014 IEEE Radar Conference. Cincinnati, Ohio, USA ,2014: 15–22. (0) |
[5] | KAWALEC A, OWCZAREK R. Radar emitter recognition using intrapulse data[C]//15th International Conference on Microwaves, Radar and Wireless Communications. Warsaw, Poland, 2004, 2: 435–438. (0) |
[6] | BERTONCINI C, RUDD K, NOUSAIN B, et al. Wavelet fingerprinting of radio-frequency identification (RFID) tags[J]. IEEE transactions on industrial electronics, 2012, 59(12): 4843-4850. DOI:10.1109/TIE.2011.2179276 (0) |
[7] | GUO Haizhao, ZHANG Xiaonu, YANG Libo, et al. Improved fisher linear discriminant analysis for feature extraction of unintentional modulation on pulse by combining ambiguity function with wavelet transform[C]//2016 International Radar Conference, IET. Guangzhou, China, 2016: 115–123. (0) |
[8] | YUAN Yingjun, HUANG Zhitao, WU Hao, et al. Specific emitter identification based on Hilbert-Huang transform-based time-frequency-energy distribution features[J]. IET communications, 2014, 8(13): 2404-2412. DOI:10.1049/iet-com.2013.0865 (0) |
[9] | KANG Naixin, HE Minghao, HAN Jun Han, et al. Radar emitter fingerprint recognition based on bispectrum and SURF feature[C]//2017 IEEE International Conference on Radar. Guangzhou, China, 2017: 71–76. (0) |
[10] | FANG Cheng, XUE Zhi. Signal classification method based on complete bispectrum and convolutional neural network[J]. Application research of computers, 2018, 35(12): 3766-3769. (0) |
[11] | LIU Lu, YANG Peiliang, SUN Weiwei, et al. Similar handwritten Chinese character recognition based on CNN-SVM[C]//International Conference on Graphics and Signal Processing. Singapore, Singapore, 2017: 81–87. (0) |
[12] | 董爱荣. 小波去噪在移动通信中的应用[D]. 大连: 大连海事大学, 2000. (0) |
[13] | 高敬鹏, 孔维宇, 刘佳琪, 等. 基于时频分析的自适应PCA辐射源调制识别[J]. 应用科技, 2018, 45(5): 33-37. (0) |
[14] | MA Ning, PENG Yu, WANG Shaojun, et al. A weight SAE based hyperspectral image anomaly targets detection[C]//IEEE International Conference on Electronic Measurement and Instruments. Yangzhou, China, 2018: 511–515. (0) |
[15] | CHAN S C, LEE K C. Radar target identification by kernel principal component analysis on RCS[J]. Journal of electromagnetic waves and applications, 2012, 26(1): 64-74. DOI:10.1163/156939312798954900 (0) |
[16] | JIA Jianhua, YANG Ning, ZHANG Chao. Object-oriented feature selection of high spatial resolution images using an improved Relief algorithm[J]. Mathematical and computer modelling, 2013, 58(3–4): 619-626. DOI:10.1016/j.mcm.2011.10.045 (0) |
[17] | LIAO Kuo, FU Jiansheng, YANG Wanlin. Modified reliefF algorithm for radar HRRP target recognition[J]. Journal of electronic measurement and instrument, 2010, 24(9): 831-836. DOI:10.3724/SP.J.1187.2010.00831 (0) |
[18] | ZHANG Yingwei. Enhanced statistical analysis of nonlinear processes using KPCA, KICA and SVM[J]. Chemical engineering science, 2009, 64(5): 801-811. DOI:10.1016/j.ces.2008.10.012 (0) |
[19] | WANG Dongxing, WANG Huibo, BAN Xiaojuan, et al. An adaptive, discrete space oriented wolf pack optimization algorithm for a movable wireless sensor network[J]. Sensors, 2019, 19(19): 4320. DOI:10.3390/s19194320 (0) |
[20] | 龚永罡, 汤世平. 面向大数据的SVM参数寻优方法[J]. 计算机仿真, 2010, 27(9): 204-207. DOI:10.3969/j.issn.1006-9348.2010.09.051 (0) |