张梓轩,齐子森,许华,史蕴豪,梁佳

(空军工程大学信息与导航学院,710077,西安)

通信辐射源信号的分类识别作为信号处理领域的研究重点之一,已成为复杂电磁环境下获取电磁情报,实现目标识别的首要步骤,对于构建复杂电磁态势有着关键作用。特定辐射源识别又称射频指纹识别[1],是利用辐射源设备的硬件差异即射频指纹在空间辐射信号中的差异性识别区分辐射源个体的技术。

在辐射源个体识别技术中,核心环节为指纹特征提取和分类器设计,提取的指纹特征即射频指纹,并且这些指纹特征应具有普遍性、可测性、唯一性和稳定性[2-3]。通信辐射源个体识别的方法大致可以分为基于专家经验的识别方法和基于深度学习的识别方法。基于专家经验的识别方法是根据对辐射源畸变的现有认知人工提取指纹特征,并进行分类识别,其中包括基于信号参数统计特征的识别方法、基于信号变换域特征的识别方法和基于发射机机理建模的识别方法。信号参数统计特征包含时域、频域以及高阶谱等。Williams等以GSM和WiMax信号的幅度、相位和频率的均值、方差等数种统计指标为特征,在对手机终端的识别上取得了较好的效果[4]。Reising等以WiFi和WiMax信号为研究对象,提取了信号中导字段的瞬时相位统计特征,对多部同型手机取得了良好的识别效果[5]。信号变换域特征主要包含小波分析、时频分析和经验模态变换。秦嘉利用小波分析的多分辨率优势,从暂态信号中提取了正交小波分析系数特征,达到了识别发射机的目的[6]。Jain等研究了暂态信号的离散Gabor变化和Gabor-Wigner变化,并从中提取幅度统计特征实现了对不同发射机的识别[7]。基于机理建模的识别方法,是从指纹特征的物理本质出发,研究发射机的组成模块和工作原理,通过建立机理模型的方式分析并提取本质特征,主要包括制造容差和漂移容差。任春辉研究并明确定义了通信电台信号个体特征,从发射机硬件实现的角度分析了通信电台个体特征(暂态特征、噪声特征、载频偏差、调制参数的差异和杂散输出的差异等)的产生机理,但未建立指纹特征的机理模型[8]。黄渊凌等系统深入地研究了通信辐射源指纹特征产生机理,阐述了I/Q正交调制器、振荡器、功放、时钟、模拟通路滤波器等发射机模块的指纹机理模型,为设计和评估指纹特征提取算法打下了理论基础,但并没有给出上述多种畸变特征的联合提取方法[9]。基于深度学习的识别方法是通过多层非线性变换对高度复杂的数据进行建模,可以高效地自主逐层提取数据的特征并进行识别。目前,基于深度学习技术的通信辐射源个体识别技术主要分为数据降维后识别、数据二维化后识别和直接从原始数据中识别。王炳程使用全连接神经网络对4种积分双谱组成的特征进行分类,实验结果较SVM、Stacking等分类方法有较大提升[10]。Pan等使用残差网络提取Hilbert-Huang变换后的频谱图特征,然后对辐射源进行识别[11]。Sankhed等设计了一种以信号的I/Q数据作为输入的一维卷积神经网络,使用一维卷积核、最大化池化层和全连接结构设计网络,取得了较好的识别效果[12]。近年来,信息论中熵的相关概念也被引入用于描述辐射源个体差异[13-14]。此外,相空间[15-16]、压缩感知[17-18]以及上述各方法的联合方法[19-20]等也被用作该方向的研究。

综上所述可知,现有的基于专家经验和基于深度学习的两类识别方法主要分析和利用射频信号时域、频域以及时频变换域的信息,但应用此类信息的方法存在不同程度的不足。信号参数统计特征的识别方法对非高斯、非平稳信号的分析能力弱[21];信号变化域的识别方法基函数选择困难;机理建模分析的识别方法建模难度大,针对性弱,无法实现多模块畸变联合分析[22]。在通信辐射源空、时、频、能域特征相近甚至相同情况下,时频信息(时域特征、频域特征、时频变换域特征)差异性缩小,利用直接或间接方法分析提取的难度都进一步增大。从原理上分析可知,现有方法识别效果会出现不同程度的恶化,甚至部分方法难以支撑此类通信辐射源个体的识别。

针对上述问题,本文引入电磁波极化域特征,拓展通信辐射源个体识别利用的特征维度,从新特征域分析研究通信辐射源个体识别方法。极化是电磁波电场强度的方向和幅度随时间变化的特性,无线电与其发射天线的极化形式保持一致。不同通信辐射源个体在天线结构和设计参数上的差异导致天线交叉极化鉴别度必然不同,这是天线设计和制造中存在的系统性偏差,可作为通信辐射源个体识别的极化特征来源。在通信辐射源个体处于空间模糊、能量相仿以及频率相近情况下,极化特征作为新维度特征为通信辐射源个体识别提供了新方向。实测分析可知,同一型号刀型天线的交叉极化鉴别度差异微小,体现在射频信号的双极化特征中,极化幅度比最小差值为1/90,极化相位差最小差值为1.5°,容易被通道幅相误差和噪声污染导致混叠。本文通过分析双极化特征的表征方法,构建了用于接收和提取双极化特征的双极化接收系统,并设计了基于自编码器的通道时变幅相误差校正算法,克服了通道时变幅相误差及噪声对双极化特征的影响,实现双极化特征的准确提取,完成了空间模糊、能量相仿及频率相近情况下的通信辐射源个体的分类识别。

目前,国内外公开发表的文献资料中未见直接利用极化域特征进行通信辐射源个体识别的方法,同时利用极化信息进行通信辐射源个体识别的研究体现较少。因此,文中未列举相关论述和研究。本文从原理分析、信号建模、极化接收、特征提取以及实验验证等方面对所提方法的可行性进行了详细论述,最大化呈现采用极化域特征进行通信辐射源个体识别的全流程,并在无法进行实验对比的情况下尽可能验证分析本文方法的有效性。

1 极化特征分析与建模

根据天线设计原理,在天线辐射的电磁波中,主体部分为主极化分量,还有少量与其正交的极化分量,也称交叉或正交极化分量。天线的交叉极化鉴别度是主极化分量的功率电平占全部功率电平的比值。发射系统交叉极化鉴别度越好,系统的性能就越理想,但实际应用中,在器件的固有误差和机械加工误差等因素影响下,天线极化总是存在一定程度的偏差,这一偏差就是通信辐射源极化特征的来源。在通信系统中,双极化天线利用极化分集技术实现了很强的抗干扰性,同时方便实现频率复用,提高信道容量,被学者广为关注研究。极化测量接收机在信号侦收中,通常采用极化-方向图法、幅度-相位法以及多天线法,目前广泛应用于现代接收系统的是幅度-相位法[23]。

1.1 极化特征表征

电磁波的极化是电磁波在空间中传播时任一位置处的电场空间取向随时间变化的性质,通过电场矢量端点随时间变化形成的空间轨迹形状和旋向来描述。根据极化状态差异,可分为完全极化波、部分极化波和非极化波3种。完全极化波是指极化状态随时间按照固定模式变化的电磁波,电场矢量端点随时间随机变化所形成的轨迹为一参数恒定的椭圆,典型信号如单色波等;部分极化波是指极化状态随时间随机变化的电磁波,电场矢量端点描绘的轨迹不再为一恒定的椭圆,而是一条形状和方向都随时间缓慢变化、类似于椭圆的曲线;非极化波是指极化电场矢量分量之间互不相关的电磁波[24]。

假设通信辐射源发射的电磁波为沿z轴正方向传播的单色完全极化电磁波,则电磁波瞬时值的矢量形式[25]为

E=Exex+Eyey

(1)

式中ex、ey分别为x轴和y轴方向的单位基矢量。Ex、Ey为

Ex=E0xcos(2πf0t-Kz+φx)

(2)

Ey=E0ycos(2πf0t-Kz+φy)

(3)

式中:f0为电磁波频率;K为空间传播常数;E0x和φx分别为x轴正方向上的电场振幅和相位;E0y和φy分别为y轴正方向上的电场振幅和相位。

将式(1)写成复矢量形式

e(t,z)=eej(2πf0t-Kz)

(4)

(5)

式中:e为复电场矢量;e0x和φx分别为x轴正方向上的电场振幅和相位;e0y和φy分别为y轴正方向上的电场振幅和相位。由式(4)、(5)得出,极化信息体现为两个正交方向电场分量的幅度比和相位差。忽略绝对相位信息,按照极化的线极化比方式[26]定义极化比

(6)

式中:tanγ=ey/ex,γ∈ [0,2π);φ=φy-φx,φ∈ [0,2π]。由此对于某一给定的电磁波,其极化状态可由极化比完全描述。

以上给出的(γ,φ)即为电磁波的极化状态相位描述子,可用来表征辐射源个体的极化特征。本文为简化计算量和更直观地表征通信辐射源射频信号的极化特征,直接利用tanγ=ey/ex,γ∈ [0,2π)作为主极化与交叉极化信号的幅度比,即极化幅度比,利用φ=φy-φx,φ∈ [0,2π]作为主极化与交叉极化信号的相位差,即极化相位差。

1.2 极化信号建模

单偶极子可以对入射信号的相应分量产生匹配响应,而任意极化波均可分解为在x、y、z这3个方向上的分量,因此电磁矢量传感器(electromagnetic vector sensor, EMVS)是由相互正交的3个电偶极子和3个磁偶极子构成,可以实现6维电磁场矢量的同时接收。假设由一个EMVS构成的阵元位于坐标中心且3个基本电偶极子方向分别与坐标轴平行,极化相位描述子为(γ,φ)的k个窄带信源sk(k=1,2,…,L)以俯仰角θ和方位角φ沿单位方向矢量u的反方向入射(如图1所示),则u可以表示为

uk=(sinθkcosφk)x+(sinθksinφk)y+(cosθk)z

(7)

在传播方向上的两个正交单位电场分量分别为

uθk=(cosθkcosφk)x+(cosθksinφk)y-(sinθk)z

(8)

uφk=(-sinφk)x+(cosφk)y

(9)

电场和磁场在直角坐标系中的极化-角度导向矢量ak可以表示为

(10)

式中:Ex、Ey、Ez表示3个电场分量;Hx、Hy、Hz表示3个磁场分量;Eθk=sinγkejφ,Eφk=cosγk分别为uθk和uφk方向的极化分量。由于极化特征计算相对值,故阵元方向图对计算无影响,所以式中单位方向图响应设置为1。

假设极化天线的共点接收通道不存在不一致、互耦等因素的影响,且增益均为1,则阵元对入射信号sk的响应为

z(t)=ak·sk(t)+n(t)

(11)

式中:n(t)为高斯白噪声;ak为极化-角度导向矢量。接收到的通信辐射源信号z(t)中包含了极化-角度导向矢量ak,通过计算ak在uθk和uφk两个方向上的极化分量Eθk=sinγkejφ、Eφk=cosγk的比值,即可得到各通信辐射源信号的双极化特征——极化幅度比和极化相位差。由于双极化特征容易被通道幅相误差和噪声污染导致混叠,因此有效利用射频信号中双极化特征的前提是低信噪比条件下,消除通道时变幅相误差精准提取信号的双极化特征。

1.3 极化信号接收与特征提取

为了在低信噪比条件下,消除通道时变幅相误差,精准提取信号极化特征,本文借鉴幅度-相位法思路,设计了双极化接收系统用以接收和提取极化特征。该系统以相互正交的两个极化天线构成的极化分集天线接收辐射源信号,通过两个相参通道对信号进行混频、放大、滤波和检波等处理,从而得出两个通道中信号的幅度比和相位差,实现通信辐射源信号的极化特征的接收、提取和输出。系统中的正交双极化天线可以是垂直极化和水平极化,也可以是左旋圆极化和右旋圆极化。系统中的两路信号的混频本振在幅度、相位保持高度一致。同时,射频放大器、中频放大器、滤波器和A/D转换的采样时钟也保持一致。通过调节混频本振的频率,该系统可以调选带宽内的通信辐射源信号,实现宽频带的搜索测频和极化特征提取。该接收系统能够同时处理到达的多个信号,并实时测算出通信辐射源信号的极化特征。双极化接收系统流程如图2所示。

图2 双极化接收系统流程Fig.2 Dual-polarized receiving system workflow

假设通信辐射源发射的信号(平面电磁波)为连续波信号,其表达式为

E=(Excos(ωt-kz+φx))ex+

(Eycos(ωt-kz+φy))ey

(12)

式中:ex和ey分别为主极化方向和交叉方向的单位基矢量;ω=2πf0t;φx和φy分别为主极化和交叉极化分量的初相。该平面波的极化状态可泛化为椭圆极化。

当选择零点距离处的信号时,式(12)可以简化为

E=(Excos(ωt+φx))ex+(Eycos(ωt+φy))ey

(13)

Ex=E0xcos(wt+φx)

(14)

Ey=E0ycos(wt+φy)

(15)

假设本振的稳定频率为f1,相位为φ1,即

Ut=U0cos(2πf1t+φ1)

(16)

主极化和交叉极化通道的信号经过混频和中频滤波后的输出为

Ex1=aE0xcos(2π(f0-f1)t+φx-φ1)

(17)

Ey1=aE0ycos(2π(f0-f1)t+φy-φ1)

(18)

式中a是信号的放大倍数。主极化和交叉极化通道的信号经过正交相位处理,分别输出两路正交的I/Q信号。正交相位处理流程如图3所示。

图3 正交相位处理Fig.3 Orthogonal phase processing

相参信号cos [2π(f0-f1)t]经过90°移相器变为Up=sin [2π(f0-f1)t]。

Ex1、Ey1分别经过正交相位处理,得到两路正交信号

ExI=gaE0xcos(φx-φ)

(19)

ExQ=gaE0xsin(φx-φ)

(20)

EyI=gaE0ycos(φy-φ)

(21)

EyQ=gaE0ysin(φy-φ)

(22)

式中g为混频和低通滤波器的增益。I、Q两路信号经过A/D转换为In、Qn(n=1,2,…,N)。

经过正交双通道线性检波器,得到主极化和交叉极化信号的幅度。

E1=gaE0x

(23)

E2=gaE0y

(24)

正交双通道线性检波流程如图4所示。

图4 正交双通道线性检波Fig.4 Orthogonal dual-channel linear demodulation

经单通道输出的相干检波器,得到主极化和交叉极化信号的相位差

V0=g2a2E0xE0ycos(φy-φx)

(25)

通过计算,得到主极化和交叉极化信号的极化幅度比和极化相位差

(26)

(27)

得到的极化幅度比和极化相位差就是通信辐射源个体的极化特征,从而可以进行分类识别。单通道输出相干检波器处理流程如图5所示。

图5 单通道输出相干检波Fig.5 Single-channel coherent detection

2 通道时变幅相误差校正

通道时变幅相误差会使接收到的极化信息出现较大的偏差,从而影响对信号双极化特征的提取,因此必须对通道时变幅相误差进行校正。通道时变误差校正算法的网络结构基于传统自编码器[27],引入去噪自编码器[28],对输入的扰动数据(含时变幅相误差和噪声)利用编码器进行特征压缩和提取,实现扰动数据的退化处理。按照输入的原始数据(只含噪声的数据)样本分布将扰动数据的部分特征值置零,提取出其中原始数据的特征,利用这些特征进行数据重构,从而实现通道时变幅相误差的校正。

2.1 极化通道误差建模

理想情况下,构成阵元的EMVS的3个电偶极子和3个磁偶极子的取向是完全相同的。但是,在实际应用中,天线平台的高强度振动会导致天线发生不规则形变,致使阵元所对应的通道并不完全平行,出现各通道的不一致性,从而产生通道时变幅相误差。

(28)

(29)

式中:ρt为时变的幅度误差;ejφt为时变相位误差。

由于通道时变幅相误差是天线载体以一定规律振动造成的,故可参考天线振动模型[29]进行设置

C(x,t)=Z0(x)+α(t)+q1(t)Z1(t)

(30)

式中:α(t)=ω0cos(ωt)为坐标原点处引入的强制振动;Z0(x)为低频的弯曲模式;q1为振动模式系数;Z1(x)为高频的第一振动模式。本文为简化模型,只考虑引入的强制振动,因此式(30)可化简为

C(t)=α(t)=ω0cos(ωt)

(31)

式中C(t)即为本文在天线搭载平台上加载的振动,也就是接收信号中的通道时变幅相误差。

因此,在共点接收通道存在不一致且增益不同为1的情况下,含有时变幅相误差的信号数据表示为

(32)

2.2 校正算法框架

本文算法包含两个阶段。第一阶段为网络训练阶段,主要包括3个模块:数据处理模块、参数学习模块和数据校正模块。数据处理模块的输入为从极化天线采集到的扰动数据和原始数据,数据处理模块进行维数变换、功率归一化、切片和加噪等操作,并划分训练集、测试集和验证集。参数学习模块主要由编码器和解码器组成,编码器对输入的扰动数据进行特征提取和特征压缩,从而得到数据的深层特征表示,解码器利用提取出的高度抽象的特征进行数据重构。数据校正模块通过最小化重构数据与原始数据的差异,实现网络参数的优化更新,如图6所示。第二阶段为测试阶段,即实际应用阶段,主要由数据处理模块和数据校正模块组成。该阶段利用训练好的网络对数据处理模块得到的扰动数据集进行深层特征提取和数据重构,实现通道时变幅相误差校正。在数据集准备过程中,扰动数据为极化天线在振动条件下采集到的数据,原始数据为极化天线在静止状态(无振动)下采集到的数据。

2.3 校正网络结构

为实现通道时变幅相误差校正,同时确保数据校正精度,设计了8层网络结构:输入层、3层卷积层、3层全连接层和输出层,如图7所示。数据处理模块的输入层将采集到的数据(扰动数据和原始数据)进行维度变换,由M×N维转换为M×2×N×1维。考虑到数据为时序数据,将网络参数学习模块的编码过程设计为3层卷积神经网络,以便更好地提取时序信息。通过二维卷积,以same模式对M×2×N×1维扰动数据进行卷积操作,保证充分提取特征的前提下,卷积后的样本尺度不变。卷积核尺寸都为1×4,特征通道数分别为8、4、2,分别提取出M×8维、M×4和M×2维特征,之后采用ReLU函数进行分别激活。ReLU激活函数采用具有分段线性的整流线性单元[30]促进了梯度的反向传播,降低了激活函数的计算量,而且其部分激活特性相当于对网络施加了稀疏正则化,一定程度上能够提高网络的鲁棒性和泛化能力。网络参数学习模块的解码过程采用3层全连接神经网络来重构数据,并通过反向传播算法最小化重构数据与原始数据的差异,通过输入的扰动数据和原始数据的双驱动,使网络学习到原始数据的深层规律,达到对通道时变幅相误差的实时校正。为有效利用编码部分提取到的潜在特征并方便数据重构,将3层全连接网络的节点数分别设计为350、250、200,前两层的激活函数为ReLu,最后一层的激活函数为tanh。经过3层全连接层输出的重构数据尺度为M×2×N×1维,再经过数据尺度变换得到M×N维重构数据。

图6 极化通道时变误差校正模型Fig.6 Time-varying error correction model for polarized channel

图7 时变误差校正网络结构Fig.7 Time-varying error correction network structure

2.4 损失函数定义及优化

(33)

(34)

3 分类判决器设计

为直观验证极化特征对于辐射源个体识别的有效性,本文设计了基于阈值规则的硬判决分类器。分类判决分为两步:第一步为极化幅度比或极化相位差特征阈值判决,利用提取出的极化幅度比和极化相位差两个极化特征,根据对双极化天线的极化特征实测数据,设置极化幅度比和极化相位差两个判决阈值,计算特征集中样本间的单一极化特征的欧式距离,经过与阈值比对进行硬判决,得到判决结果;第二步为双极化特征联合阈值判决,将极化幅度比和极化相位差两个极化特征联合,计算极化特征集中样本间的欧式距离,通过与计算得出的联合阈值进行比对判决,最终得出判决结果。具体的判决步骤如下。

4 通信辐射源个体分类识别流程

采用极化特征的通信辐射源个体分类识别主要分为3部分:极化特征分析与建模、通道时变幅相误差校正以及辐射源分类判决。第一部分是极化特征分析与建模,首先经过极化原理分析,提出辐射源天线交叉极化鉴别度的差异性,然后通过双极化特征表征推导出辐射源极化特征相关量,再由极化信号建模分析出极化信号在接收机前端的信号形式,最后通过双极化接收系统实现极化信号接收与双极化特征提取。第二部分是通道时变幅相误差校正,设计了基于自编码器的时变幅相误差校正算法,利用深度神经网络对含有时变幅相误差扰动及噪声的数据进行校正,重构受扰动和噪声影响的微小极化特征。第3部分是辐射源分类判决,利用设计的基于阈值规则的硬判决器,对提取的极化幅度比和极化相位差两个极化特征进行硬判决,通过极化幅度比或极化相位差特征判决和极化幅度比-极化相位差联合特征判决相结合,将极化特征的差异性充分利用,从而实现通信辐射源个体的分类识别。

本文提出的采用极化特征的通信辐射源个体识别方法的流程图如图8所示,具体流程如下。

第1步极化特征分析与建模。

(1)通过对通信辐射源天线的交叉极化鉴别度差异性分析,推导其发射电磁波的极化特征表征。

(2)利用极化特征表征涉及的相关参数,进行极化信号建模。

(3)根据极化信号建模,设计双极化接收系统进行极化特征提取,具体流程如下。

①通过极化分集天线实现主极化和交叉极化两个方向的极化信号接收。

②采用可调本振对宽带内的通信辐射源信号进行搜索测频并实现混频。

③采用相参振荡器进行主极化和交叉极化信号的正交相位处理。

④经过A/D转换实现数字采样。

⑤利用正交双通道线性检波器进行包络检波,提取主极化和交叉极化信号的幅度值。

⑥利用单通道输出的相干检波器进行相位检波,提取主极化和交叉极化信号的相位值。

⑦通过计算得出辐射源信号的极化幅度比和极化相位差两个极化特征。

第2步通道时变幅相误差校正。

(1)通过分析通道时变幅相误差,建立通道时变幅相误差模型。

(2)利用原始数据(有噪声不含时变幅相误差的信号数据)以及扰动数据(同时含有时变幅相误差和噪声的信号数据)对校正网络进行预训练。

(3)保存训练好的网络结构和参数。

(4)利用训练好的校正网络对扰动数据(含有时变幅相误差和噪声的信号数据)进行校正。

第3步辐射源分类判决。

(1)根据在暗室对某刀型天线的交叉极化鉴别度的测量结果,得出主极化和交叉极化的最小极化幅度比阈值和最小极化相位差阈值,并计算得出联合阈值。

(2)利用样本集中提取到的极化幅度比和极化相位差特征,分别计算一维空间中样本间的欧式距离。

(3)利用样本集中提取到的极化幅度比-极化相位差联合特征,计算二维空间中样本间的欧式距离。

(4)将计算得到的欧式距离分别与对应阈值进行比对判决,根据判决得出样本判决结果,实现辐射源分类。

图8 采用极化特征的通信辐射源分类识别流程Fig.8 Communication emitter classification and recognition process using polarization features

5 实验分析

为充分验证所提方法的有效性,对中间校正网络和整体识别效果分别进行了实验验证。实验基于仿真数据和实采数据在计算机上进行。仿真数据利用Matlab软件生成。参照实际双极化天线结构和通信流程,模拟生成主极化和交叉极化两个通道的二进制相移键控(BPSK)调制信号,得到两组相互正交的I/Q信号,构建实验数据集。实采数据通过微波暗室采集得到,共有两类数据。一类为原始数据,即在暗室中极化天线静止状态下采集的高信噪比数据。另一类为扰动数据,即在暗室中极化天线振动时采集到的数据。

5.1 通道时变幅相误差校正网络训练验证

为验证在实际应用环境中采集到的数据可以用作校正网络的有效训练,模拟理想状态和实际状态中的数据采集环境,分别使用无噪声原始数据和20 dB信噪比原始数据分别对网络进行训练,得到理想网络与实际网络。

图9为实际网络训练收敛曲线。通过图9中的(a)和(b)说明实际网络在20 dB信噪比的数据下可以快速收敛,在主极化和交叉极化两个通道上都可以达到很好的训练效果,证明了网络在不同通道上的适应性较好,可以实现对双极化通道时变幅相误差的校正。

(a)主极化通道

(b)交叉极化通道

按照表1中的参数采集10 dB信噪比下的扰动数据,在以上两个完成训练的网络中分别进行数据重构,结果如图10和图11所示。

表1 模拟数据集参数

实际网络和理想网络在10 dB信噪比下扰动数据重构结果分别如图10和11所示。两图对比说明,在实际环境中采集20 dB原始数据训练的网络可以基本达到与理想无噪声原始数据训练的网络同等的效果,证明本文校正网络具有消噪能力,在实际环境中可以用作数据重构和消噪处理。

(a)实部

(b)虚部

(a)实部

(b)虚部

5.2 通道时变幅相误差校正效果验证

为了验证校正网络在实际情况下对极化通道的校正效果,在表1参数前提下,将接收端增益分别设置为-5、0、5、10、15、20 dB采集数据,利用5.1小节中20 dB数据预训练好的网络逐一进行数据重构,100次蒙特卡罗实验的结果如图12～13所示。

图12是主极化通道在10 dB信噪比下对扰动数据的重构结果。从图12中原始数据、重构数据和扰动数据的对比可以得出,主极化通道中在信噪比10 dB下的扰动数据经过校正网络后,可以很好地拟合原始数据,同时网络对数据进行了一定程度的噪声抑制,达到了比原始数据更好的效果。

(a)实部

(b)虚部

从图13中原始数据、重构数据和扰动数据对比同样可以得出,交叉极化通道中网络能够很好地压制噪声并对扰动数据进行重构,证明了网络对时变幅相误差的校正能力。

(a)实部

(b)虚部

表2是不同信噪比下主极化和交叉极化通道的归一化实部和虚部均方误差,表2中的数据进一步展示了校正网络的性能。

表2 不同信噪比下各通道实部和虚部均方误差

可以看出,-5～20 dB的范围内,网络能够对主极化和交叉极化两个通道进行不同程度的数据重构和噪声抑制。5 dB以上时:主极化通道校正数据的实部(Ⅰ路信号)均方误差能够控制0.3%以内,虚部(Q路信号)均方误差能够基本保持在60%左右;交叉极化通道校正数据的实部(Ⅰ路信号)均方误差能够控制0.5%以内,虚部(Q路信号)均方误差能够基本保持在6%左右。证明了网络具有较好的鲁棒性,适用的信噪比范围较大。

5.3 采用极化特征的通信辐射源个体识别效果验证

为验证本文方法的整体效果,通过模拟数据和实采数据分别对极化特征在通信辐射源个体识别的效果进行了验证。对同一厂家、同一批次、同一型号的刀型天线进行交叉极化鉴别度实测,得出双极化天线的主极化和交叉极化参数最小差异值。极化幅度比最小差异值归一化为1/90,极化相位差最小差异值为1.5°。模拟设置两个极化特征不同的辐射源个体,极化幅度比分别为1∶10和1∶9,极化相位差分别为2°和3.5°,其余参数完全相同。在表1参数前提下,将接收端的增益分别设置为-5、0、5、10、15、20 dB采集数据,对采集到的400 帧(40 000采样点)数据进行100次蒙特卡罗实验,得到辐射源分类识别结果,如图14～16所示。个体识别准确是指个体类别分类正确,识别准确率是指个体样本识别正确的数量占个体总样本数的比例。

(a)极化幅度比特征

(b)极化相位差特征

在硬判决规则下分别采用极化幅度比和极化相位差特征,对不同信噪比下通信辐射源个体原始数据和校正数据进行识别,识别准确率如图14所示。可以看出:随着信噪比的提升,本文方法对辐射源的识别准确率逐步提高。在10 dB时,利用极化幅度比和极化相位差特征的识别准确率分别可达98%和75%以上;在20 dB时,利用幅度比极化特征的个体识别准确率达到100%。重构数据较原始数据识别效果更好,证明了本文方法的有效性和可靠性。

10 dB和20 dB信噪比下,采用极化幅度比、极化相位差以及极化幅度比和极化相位差时,个体极化特征分布情况分别如图15和图16所示。

由图15和图16可以得出,单纯的极化幅度比特征相较于极化相位差特征在个体区分度上效果更好。在信噪比较低时,极化幅度比特征仍然可以清晰地区分出两类目标,而极化相位差特征分类界限不明显。在信噪比为10 dB以上时,综合使用极化幅度比和极化相位差两个特征可以较好地实现通信辐射源个体的分类识别。

(a)极化幅度比特征

(c)极化幅度比-极化相位差特征

(a)极化幅度比特征

(b)极化相位差特征

(c)极化幅度比-极化相位差特征

除此之外,在暗室理想条件下,对同一厂家、同一批次、同一型号的10副某刀型天线进行数据实采,并使用本文方法进行分类识别。结果表明,本文方法对10个天线个体的识别准确率均达到99%以上。在实采数据上叠加高斯白噪声,模拟现实应用场景中接收机在信噪比10 dB下采集到的辐射源信号。分类识别实验结果表明,本文方法对10个天线个体的识别率均达到了95%以上。

6 结 论

本文针对通信辐射源个体在空、时、频、能域特征下难以准确分类识别的问题,提出了一种采用极化特征的通信辐射源个体识别方法。推导了支撑通信辐射源个体识别的双极化特征表征方法,构建了双极化接收系统,实现了极化信号的接收和双极化特征的提取,并通过建模分析极化通道的时变幅相误差,设计了基于自编码器的通道时变幅相误差校正算法,克服了双极化特征对通道噪声与时变幅相误差敏感的问题,有效实现了采用极化特征的通信辐射源个体识别。在仿真实验与暗室实测条件下,验证了方法的有效性。通过分析可知,低信噪比条件和变调制样式情况下,不同辐射源个体的极化特征混叠明显。下一步将考虑以下3个方向的研究:①低信噪比条件下有效提取并运用极化特征进行辐射源个体的分类识别;②面对不同调制样式场景下,采用极化特征对通信辐射源个体进行有效识别;③跳频通信个体在不同组网方式下利用极化特征实现跳频网台的有效分选。