基于多域融合的空间辐射源细微特征识别.pdf

中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4 1 2 6-1 3 6C h i n e s e S p a c e S c i e n c e a n d T e c h n o l o g yI S S N1 0 0 0-7 5 8 X C N1 1-1 8 5 9/Vh t t p:z g k j.c a s t.c nD O I:1 0.1 6 7 0 8/j.c n k i.1 0 0 0-7 5 8 X.2 0 2 3.0 0 6 1基于多域融合的空间辐射源细微特征识别王晓晗1,2,闫毅1,*,范亚楠1,李雪1,牟娇11.中国科学院 国家空间科学中心,北京1 0 0 1 9 02.中国科学院大学,北京1 0 0 0 4 9摘 要:辐射源识别主要关注辐射源的个体差异,通过信号处理手段,获取辐射源信号上与设备硬件相关的特征参数,从而掌握辐射源设备的型号、工作状态等信息。由于通信信号愈发复杂,单一域特征识别无法全面反映信号的细微差异,直接对信号处理结果进行识别存在大量数据冗余,导致识别效果不佳。为提高空间辐射源的识别效果,提出了一种基于多域特征融合的辐射源识别方法,通过对信号进行H i l b e r t-H u a n g变换和高阶谱分析,并提取变换结果的均值、方差、峰度、偏度和信息熵,将其拼接、融合形成空间辐射源信号的多域特征,利用S VM分类器进行分类,实现多域融合的辐射源细微特征识别。仿真结果表明,使用多域融合方法进行辐射源识别,在2 0d B的信噪比下可达到9 5.4 4%的正确率,与传统基于边际谱信息熵的分类识别方法对比,正确率提升9.4 1%。对比基于HHT边际谱、边际谱信息熵、双谱投影、双谱矩形积分、四阶累积量切片谱的分类识别方法,本方法的识别效果均有提升。关键词:空间通信;特征融合;空间辐射源识别;H i l b e r t-H u a n g变换;高阶谱分析中图分类号:T N 9 7 5 文献标识码:A收稿日期:2 0 2 2-0 2-2 8;修回日期:2 0 2 2-0 4-2 5;录用日期:2 0 2 2-0 8-1 0;网络出版时间:2 0 2 2-1 1-0 8 1 5:0 6基金项目:中国科学院战略高技术创新基金项目(GQ R C-1 9-1 4)*通信作者.E-m a i l:y a n y i n s s c.a c.c n引用格式:王晓晗,闫毅,范亚楠,等.基于多域融合的空间辐射源细微特征识别J.中国空间科学技术,2 0 2 3,4 3(4):1 2 6-1 3 6.WAN G X H,YAN Y,F AN Y N,e t a l.S p a c e e m i t t e r f i n e f e a t u r e i d e n t i f i c a t i o n b a s e d o n m u l t i-d o m a i n f u s i o nJ.C h i n e s e S p a c e S c i e n c e a n d T e c h n o l o g y,2 0 2 3,4 3(4):1 2 6-1 3 6(i n C h i n e s e).S p a c e e m i t t e r f i n e f e a t u r e i d e n t i f i c a t i o n b a s e d o n m u l t i-d o m a i n f u s i o nWA N G X i a o h a n1,2,Y A N Y i1,*,F A N Y a n a n1,L I X u e1,MO U J i a o11.N a t i o n a l S p a c e S c i e n c e C e n t e r,C A S,B e i j i n g 1 0 0 1 9 0,C h i n a2.U n i v e r s i t y o f C h i n e s e A c a d e m y o f S c i e n c e s,B e i j i n g 1 0 0 0 4 9,C h i n aA b s t r a c t:Em i t t e r i d e n t i f i c a t i o n m a i n l y f o c u s e s o n i n d i v i d u a l d i f f e r e n c e s o f e m i t t e r s,a n d o b t a i n s c h a r a c t e r i s t i c p a r a m e t e r s r e l a t e d t o e q u i p m e n t h a r d w a r e o n e m i t t e r s i g n a l s t h r o u g h s i g n a l p r o c e s s i n g m e a n s,s o a s t o m a s t e r t h e m o d e l a n d w o r k i n g s t a t u s o f e m i t t e r e q u i p m e n t.A s c o mm u n i c a t i o n s i g n a l s b e c o m e m o r e a n d m o r e c o m p l e x,f e a t u r e i d e n t i f i c a t i o n i n a s i n g l e d o m a i n c a n n o t f u l l y r e f l e c t t h e f i n e d i f f e r e n c e s o f s i g n a l s,a n d t h e r e i s a l a r g e a m o u n t o f d a t a r e d u n d a n c y i n d i r e c t i d e n t i f i c a t i o n o f s i g n a l p r o c e s s i n g r e s u l t s,r e s u l t i n g i n p o o r i d e n t i f i c a t i o n e f f e c t.I n o r d e r t o i m p r o v e t h e i d e n t i f i c a t i o n e f f e c t o f s p a c e e m i t t e r,a n e m i t t e r i d e n t i f i c a t i o n m e t h o d w a s p r o p o s e d b a s e d o n t h e m u l t i-d o m a i n f e a t u r e f u s i o n.B y p e r f r o m i n g H i l b e r t-H u a n g t r a n s f o r m a n d h i g h e r o r d e r s p e c t r u m a n a l y s i s o f t h e s i g n a l,t h e m e a n,v a r i a n c e,k u r t o s i s,s k e w n e s s a n d i n f o r m a t i o n e n t r o p y o f t h e t e a n s f o r m e d r e s u l t s a r e e x t e a c t e d,S VM c l a s s i f i e r i s u s e d t o c l a s s i f y e m i t t e r s a n d r e a l i z e t h e f i n e f e a t u r e i d e n t i f i c a t i o n o f m u l t i-d o m a i n f u s i o n.T h e s i m u l a t i o n r e s u l t s s h o w t h a t t h e a c c u r a c y o f m u l t i-d o m a i n f e a t h e r f u s i o n m e t h o d c a n r e a c h 9 5.4 4%u n d e r t h e S N R o f 2 0d B,a n d t h a t t h e a c c u r a c y i s i m p r o v e d b y 9.4 1%c o m p a r e d w i t h t h e t r a d i t i o n a l c l a s s i f i c a t i o n m e t h o d b a s e d o n m a r g i n a l s p e c t r u m i n f o r m a t i o n e n t r o p y.C o m p a r e d w i t h t h e c l a s s i f i c a t i o n a n d i d e n t i f i c a t i o n m e t h o d s b a s e d o n HHT m a r g i n a l s p e c t r u m,m a r g i n a l 王晓晗,等:基于多域融合的空间辐射源细微特征识别1 2 7 s p e c t r u m i n f o r m a t i o n e n t r o p y,b i s p e c t r a l p r o j e c t i o n,b i s p e c t r a l r e c t a n g u l a r i n t e g r a l a n d f o u r t h-o r d e r c u m u l a n t s l i c e s p e c t r u m,t h e i d e n t i f i c a t i o n e f f e c t o f t h e p r o p o s e d m e t h o d i s i m p r o v e d.K e y w o r d s:s p a c e c o mm u n i c a t i o n;f e a t u r e f u s i o n;s p a c e e m i t t e r i d e n t i f i c a t i o n;H i l b e r t-H u a n g t r a n s f o r m;h i g h e r o r d e r s p e c t r a1 引言辐射源个体识别在电子对抗和通信安全等方面起到至关重要的作用,是非合作通信侦察中的一项关键技术1。辐射源识别主要根据所截获射频信号的差异判断所属辐射源2,根据辐射源工作状态掌握其载体的情况,从而有针对性地对其进行监视和管理3-4。卫星通信不受地域地形的限制,具有通信距离远、全球覆盖、安全可靠的优势,研究空间辐射源识别技术,对通信卫星的有效的侦察、监测和管理具有重要意义。由于辐射源设备的个体差异,在信号上不可避免地存在着不影响信息传递的细微特征差异,这些可检测、可重现的细微差别被称为信号细微特征或辐射源个体细微特征5。射频设备中电子元件的容差是形成细微特征的原因,元件在生产时实际值与标称值存在差异,在工作过程中温度、湿度、老化也会引起参数的变化。目前认为,空间辐射源细微特征的主要成因是卫星转发器中功率放大器的非线性效应6。通信卫星与地面站的距离遥远,为了抵消路径损耗,保证功率放大器较高的工作效率,卫星转发器中的功率放大器通常工作在饱和状态或邻近饱和状态,导致输出的信号产生严重的非线性失真。由放大器产生的非线性产物主要有谐波、交调、互调和杂散等7-9。因此,由于元件容差的存在,系统的非线性特性不同,对信号的影响也随之不同,这给依据信号中提取细微特征,进而识别出相应的辐射源个体提供了可行性。目前,国内外针对辐射源识别问题,主要针对暂态信号与稳态信号两类信号进行研究。暂态信号是指辐射源在开、关机和模式转换时由于工作状态不稳定所产生的能体现辐射源特征的信号,因其持续时间短且受噪声干扰影响大,不易于进行辐射源识别。稳态信号是指辐射源稳定工作时所辐射的信号,存在于整段信号的信息中,易于获取并且能够更全面描述发射机信号的特征,因此是辐射源识别的研究重点。近年来,众多学者针对稳态信号的提取提出了多种处理方法,主要包括贝叶斯估计1 0、信号包 络1 1、二 维 时 频 特 征1 2、W i g n e r-V i l l e分布1 3、小波变换1 4、固有时间尺度分解1 5、变分模态分解1 6、H i l b e r t-H u a n g变换1 7-1 8、双谱分析1 92 1等理论。文献2 2 利用雷达信号的累积量作为个体识别的特征,使用KNN分类器对3部同型号机载雷达的实测数据进行分类识别,取得良好的分类效果。文献2 3 使用离散傅里叶变换的功率谱构建了谱域特征,结果表明谱域特征比时域特征性能好,与小波域特征性能接近,并略有提高。文献2 4 使用G a b o r变换系数构建了 基 于G a b o r变 换 的 特 征,结 果 表 明 基 于G a b o r变换的特征比时域和谱域特征更为有效。直接提取信号变换的结果进行辐射源识别,存在大量的无效或冗余数据。无效数据对识别效果贡献不足,甚至会产生负面影响,冗余数据会增大计算量,影响辐射源识别效率。这些方法大多使用单一信号处理手段提取特征,难以全面准确地表示辐射源信号之间的细微差别,存在一定的局限性。随着信号体制和电磁环境的复杂化,单一域的特征难以完整体现信号的特性,将这些特征进行融合,建立更高维度的特征向量集合,可以全面表达信号的特征,获得更优的识别效果。为解决现有分类方法2 5-2 7数据冗余和特征单一对分类效果的负面影响,本文综合现有方法与成果,提出一种多域特征融合(m u l t i-d o m a i n f e a t u r e f u s i o n,MF F)的空间辐射源识别方法,通过对信号进行H i l b e r t-H u a n g变换和高阶谱分析,然后提取HHT边际谱和双谱的特征并将其拼接融合为信号的多域特征,采用S VM分类器进行分类识别,进行仿真实验,验证方法的性能。1 2 8 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.42 辐射源细微特征识别2.1 辐射源细微特征分析处理流程通信设备细微特征提取与识别过程如图1所示,辐射源设备发送射频信号,通信接收机收到信号后,送入预处理机中,预处理的主要作用是对信号进行去噪处理,以削弱噪声的影响,加强有用信息;预处理后的数据将送入主处理机进行信号处理和特征提取,生成特征数据集送入分类识别模块;在分类识别模块,将已标注的特征数据集送入分类器中进行训练,再将待识别信号特征数据集送入训练后的分类器中,得到分类识别结果。图1 通信设备细微特征提取与识别流程F i g.1 I d e n t i f i c a t i o n p r o c e s s o f e m i t t e r f i n e f e a t u r e s2.2 功放模型及细微特征功率放大器是通信发射机中的重要器件,为保证较高的工作效率,卫星转发器内部的高功率放大器通常工作在饱和状态或邻近饱和状态,导致输出的信号存在严重的非线性失真,是辐射源细微特征的主要来源。一个双音复信号x(t)可以表示为:xt =2i=1aie x p j2 fit+i 式中:ai为信号的幅度;fi为信号的频率;i为信号的相位。假设功率放大器输入信号为x(t),输出信号y(t)可表示为:yt =2i=1A ai e x p j2 fit+i +ai 式中:A ai 为功率放大器由于非线性特性引起的幅度AM-AM失真;ai 为相位AM-PM失真。卫星通信系统中,常采用行波管放大器作为转发器内部的高功率放大器。S a l e h模型是针对行波管提出的失真模型2 8,对幅度和相位失真均有描述:A r =1r1+1r2,r =2r21+2r2 图2给出了功率放大器的AM-AM、AM-P M特性曲线,体现了功率放大器的非线性特性。图2 功率放大器的AM-AM和AM-PM曲线F i g.2 AM-AM a n d AM-PM c u r v e s o f p o w e r a m p l i f i e r信号经过非线性功率放大器,会产生多种非线性产物,这些产物中蕴含着信号的细微特征。由功放产生的细微特征主要包含谐波、交调、互调、杂散等成分,影响信号的幅度、频率和相位。对于每一台辐射源设备,由于其内部的电子元件存在差异,导致信号的细微特性不同,从而形成辐射源个体识别的依据。3 特征设计与分类识别方法3.1 H i l b e r t-H u a n g变换H i l b e r t-H u a n g变 换 利 用 经 验 模 态 分 解(e m p i r i c a l m o d e d e c o m p o s i t i o n,EMD)对原始信号进行处理,分解原始信号得到若干个本征模态函数(i n t r i n s i c m o d e f u n c t i o n,I MF),再对得到的I MF进行H i l b e r t变换得到瞬时频率,从而获得信号的H i l b e r t时频谱2 9。在经验模态分解过程中,信号s(t)被分解为I MF分量ci(t)与残留分量rn(t)之和,即:s(t)=ni=1cit +rnt 图3为信号波形及其I MF分量,其中s i g n a l为原始信号s(t),之后依次为EMD分解之后得王晓晗,等:基于多域融合的空间辐射源细微特征识别1 2 9 到I MF分量ci(t),分别称作I MF 1I MF 6,最后为残留分量rn(t),每一个I MF分量代表了原始信号中存在的一种模态分量。I MF 1频率最高,特征时间尺度最小。随着分解层数的增加,I MF分量的特征时间尺度增大,频率降低,多分量信号中的不同频率成分被分离出来。图3 信号及其I MF分量F i g.3 S i g n a l a n d I MF c o m p o n e n t sEMD分解后的各个I MF分量经过H i l b e r t变换可以求得瞬时频率。原始信号可以表示为:st =R eni=1Ait ejit =R eni=1Ait ejit dt 反映了信号幅值与时间、频率的关系,信号幅值可以表 示 为 时 间t和 瞬 时 频 率的 函 数H,t ,从而获得信号总的HHT时频谱:H,t =ni=1Hi,t 图4所示为仿真信号的HHT时频谱图,由于EMD能够实现对信号成分的精细分析,H(,t)频率分辨率很高,信号的不同频率分量以及低频噪声分量清晰地显示出来,能较好反映信号的细微特征。HHT边际谱由HHT时频谱沿时间轴进行积分得到:J =T0H,t dt图4 信号的HHT时频谱F i g.4 H i l b e r t-H u a n g s p e c t r u m o f s i g n a l 图5为 仿 真 信 号 的HHT边 际 谱 图,是H i l b e r t时频谱在时间轴上的积分,它反映了每一个频率点上的幅值分布,代表信号的幅值在整个频率段上随频率变化的情况,反映了概率意义上幅值在整个数据跨度上的积累幅值。与傅里叶分析不同,傅里叶变换的频谱中某一频率处能量的存在代表该频率在整个时间轴上的存在,而HHT边际谱数值的大小表示该频率出现的概率高低,该频率出现 的具体时刻由 图 4中 的HHT时频谱表示。图5 信号的HHT边际谱F i g.5 H i l b e r t-H u a n g m a r g i n a l s p e c t r u m o f s i g n a l3.2 高阶谱分析高阶谱是二维功率谱的推广,是高阶累积量的傅里叶变换,可以提供二阶统计量所没有的相位信息,适用于对信号细微特征的研究。定义一个随机序列为:x n ,x n+1 ,x n+k-1 其高阶累积量ck x1,2,k-1 满足:1 3 0 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.41=-k-1=-ck x1,2,k-1 则k阶谱是k阶累积量的k-1维离散时间傅里叶变换,即:Sk x1,2,k-1 =+1=-+k-1=-ck x1,2,k-1 e-j 11+k-1k-1 因此,三阶谱又称为双谱,四阶谱又称为三谱,分别定义如下:Bx1,2 =+1=-+2=-c3x1,2 e-j 11+22 Tx1,2,3 =+1=-+2=-+3=-c4x1,2,3 e-j 11+22+33 双谱是阶数最低的高阶谱,在特征提取、故障检测方面有着较为广泛的应用,形成了较为成熟的理论和估计方法。图6是非线性信号在S N R=0 d B条件下的双谱图。图7所示的双谱图是用等高线来定性刻画的,完整地保留了信号的双谱特征,最大限度地保留了信号原有的信息,方便观察不同辐射源信号高阶谱的差异。图6 信号的双谱F i g.6 B i s p e c t r u m o f s i g n a l3.3 MF F辐射源识别方法单一域特征往往无法全面体现信号复杂丰富的特征信息,导致识别效果不佳。虽然单一域特征在分类识别中存在局限性,但将多种特征进行拼接融合,可以发挥特征在不同域的作用。多域融合特征设计的目的是尽量获得足够多的特征,构建新的多域特征向量空间,以更全面描述辐射源信号的细微特征,并通过多域特征表征不同信号的差异性,将辐射源识别的问题转换到高维特征空间中的判决区域划分,以实现更优的识别效果。图7 通过F F T估计的双谱F i g.7 B i s p e c t r u m e s t i m a t e d v i a t h e F F T m e t h o d统计量是样本的函数,用于描述或者构建样本之间的关系,从样本数据中提取信息并推断总体信息。信息熵可以表征样本数据能量的分布情况。有效的辐射源细微特征应当具有较强的判别性以及较低的冗余性。根据上述原则,本文设计了辐射源多域特征集合,如图8所示。对于一个随机变量X,数学期望为D(X)。其方差定义为:D X =EX-E X 2 峰度定义为:K u r tX =E X-E X 4E X-E X 2 2-3偏度定义为:S k e wX =E X-E X 3E X-E X 32 使用EMD算法分解信号,并进行H i l b e r t变换,求得HHT边际谱J 。基于HHT边际谱的特征设计包括以下统计特征。王晓晗,等:基于多域融合的空间辐射源细微特征识别1 3 1 图8 辐射源多域特征集合F i g.8 M u l t i-d o m a i n f e a t u r e s e t o f e m i t t e r1)HHT边际谱均值AVGJ为:AVGJ=E J 2)HHT边际谱方差VA RJ为:VA RJ=D J 3)HHT边际谱信息熵H s eJ为:H s eJ=-fJ l nJ 4)HHT边际谱峰度K u r tJ为:K u r tJ=K u r tJ 5)HHT边际谱偏度S k e wJ为:S k e wJ=S k e wJ 定义双谱矩阵为:Bx1,2 =ijbi j双谱投影为:P=ibi j基于双谱的特征设计包括以下统计特征。1)双谱均值AVGB为:AVGB=E Bx1,2 2)双谱方差VA RJ为:VA RB=D Bx1,2 3)双谱信息熵H s eB为:H s e=-ijbi jl nbi j 4)双谱投影峰度K u r tP为:K u r tP=K u r tP 5)双谱投影偏度S k e wP为:S k e wP=S k e wP 在特征融合的基础上,本文设计了基于M F F的辐射源识别方法,流程如图9所示。原始信号集合包含已标注信号和未知信号,将其送入非线性元件并加以高斯白噪声,得到有噪声的失真信号集合;对失真信号集合分别进行H i l b e r t-H u a n g变换和高阶谱分析,得到信号的HHT边际谱和双谱;然后进行特征提取和特征融合,提取HHT边际谱的均值、方差、峰度、偏度和信息熵作为信号HHT时频域特征,提取双谱均值、双谱方差、双谱投影峰度、双谱投影偏度和双谱信息熵作为信号的双谱域特征,再将其拼接融合得到多域特征集合;最后将已标注的多域特征送入S VM分类器中进行训练,将待识别的多域特征送入训练好的分类器中,得到分类识别结果。图9 基于MF F的辐射源识别方法流程F i g.9 MF F-b a s e d e m i t t e r i d e n t i f i c a t i o n m e t h o d4 试验设置与结果分析4.1 试验设置为了评价多域融合的辐射源识别方法,本次1 3 2 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4仿真试验对发射端的非线性失真进行了仿真,功率放大器采用S a l e h模型,模拟了信号的谐波、交调、互调、杂散等非线性产物特征,将本文提出的方法与传统基于HHT边际谱、边际谱信息熵、双谱投影、双谱矩形积分、四阶累积量切片谱的分类识别方法进行了比较。实际的通信卫星工作在L、C、S、K u及K a频段,在试验设置中,使用不失代表性的简单模型,参考I n m a r s a t和T h u r a y a卫星通信系统,卫星转发器带宽3 6MH z,共划分为1 2 0个子带,多路转发器相邻通道的频率间隔为3 0 0k H z,1 0MH z、1 0.3MH z对应原始信号的两个射频频点,发送到卫星转发器的两个通道中,然后经过变频、放大到下行1 6 2 6MH z、1 6 2 6.3MH z两个频点,在卫星地面站接收到高频信号,将其下变频至1 0MH z、1 0.3MH z的 中 频 频 段,在 中 频 接 收 机 选 取8 0MH z的频率对信号进行采样,选取每段信号的长度为N=4 0 0。设置辐射源个数为4个,辐射源输入相同的信号并加以高斯白噪声,信噪比的设置S N R0,2 0d B,其中,S N R仿真步长为1 d B。在仿真信号中,每个信噪比下选取5 0段构成有标签的训练数据集,并在同样的信噪比设置下分别选取5 0段构成测试数据集。然后进行预处理,对训练集和测试集信号进行小波阈值去噪,对小波系数进行阈值处理。因为小波变换后有用信号的能量主要集中在少数绝对值较大的小波系数上和特定的频率范围内,而噪声的能量却分布于整个小波域内。小波分解后的有用信号的系数幅值会大于噪声的系数幅值,采用阈值处理的方法可实现信号去噪。4.2 试验结果分析为验证试验结果的有效性,进行如下试验,将本文提出的多域融合识别方法与传统方法进行对比,将训练集输入到S VM分类器中,训练不同的分类器,将测试集输入分类器中进行分类,计算正确率,表1为各种方法分类器的输入。边际谱信息熵描述信号边际谱的不确定性;四阶累积量切片谱是四阶谱在一维的投影;双谱投影和双谱矩形路径积分通过不同的积分路径将二维双谱通过积分转化为一维数据进行处理;传统方法常把上述特征作为衡量信号细微特征的重要参数。表1 试验设置T a b l e 1 E x p e r i m e n t s e t t i n g sC l a s s i f i c a t i o n m e t h o dC l a s s i f i e r i n p u tM e t h o d AHHT m a r g i n a l s p e c t r u mM e t h o d BM a r g i n a l s p e c t r a l e n t r o p yM e t h o d CA x i a l i n t e g r a l b i s p e c t r u mM e t h o d DS q u a r e i n t e g r a l b i s p e c t r u mM e t h o d EF o u r t h-o r d e r c u m u l a n t s l i c e s p e c t r u mM e t h o d o f t h i s a r t i c l eM u l t i-d o m a i n f e a t u r e f u s i o n S VM算法最初是为二值分类问题设计的,当处理多类问题时,就需要构造合适的多类分类器。本文采用一对一 S VM,在任意两类样本之间设计一个S VM分类器,k个类别的样本就需要设计k(k-1)/2个S VM分类器。当对测试数据进行分类时,对所有的分类器分别进行测试,并保存结果,最后得票最多的类别即为该样本的类别。对去噪后的信号进行HHT变换,得到信号的HHT边际谱,图1 0是在S N R=2 0 d B下4台发射机信号的HHT边际谱。图1 0 S N R=2 0 d B信号的HHT边际谱F i g.1 0 H i l b e r t-H u a n g m a r g i n a l s p e c t r u m o f s i g n a l w h e n S N R=2 0 d B提取HHT边谱的均值、方差、峰度、偏度、信息熵特征,为方便显示,固定三维图像的X轴代表HHT边际谱的均值,Y轴代表方差,Z轴分别代表HHT边际谱的偏度、峰度和信息熵,绘制特征空间分布散点图,如图1 11 3所示。王晓晗,等:基于多域融合的空间辐射源细微特征识别1 3 3 图1 1 HHT边谱均值-方差-偏度F i g.1 1 M e a n v a r i a n c e a n d s k e w n e s s o f H i l b e r t-H u a n g m a r g i n a l s p e c t r u m图1 2 HHT边谱均值-方差-峰度F i g.1 2 M e a n v a r i a n c e a n d k u r t o s i s o f H i l b e r t-H u a n g m a r g i n a l s p e c t r u m图1 3 HHT边谱均值-方差-熵F i g.1 3 M e a n v a r i a n c e a n d e n t r o p y o f H i l b e r t-H u a n g m a r g i n a l s p e c t r u m然后计算信号的双谱,图1 4为S N R=2 0 d B下4台发射机信号的双谱图像。提取双谱均值、双谱方差、双谱投影峰度、双谱投影偏度和双谱熵,为方便显示,固定三维图像图1 4 S N R=2 0 d B 信号的双谱图像F i g.1 4 B i s p e c t r u m o f s i g n a l w h e n S N R=2 0 d B的X轴代表双谱的均值,Y轴代表双谱方差,Z轴分别代表双谱投影峰度、双谱投影偏度和双谱熵,绘制特征空间分布散点图,如图1 5 1 7所示。图1 5 双谱均值-方差-偏度F i g.1 5 M e a n v a r i a n c e a n d s k e w n e s s o f B i s p e c t r u m图1 6 双谱均值-方差-峰度F i g.1 6 M e a n v a r i a n c e a n d k u r t o s i s o f B i s p e c t r u m1 3 4 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4图1 7 双谱均值-方差-熵F i g.1 7 M e a n v a r i a n c e a n d e n t r o p y o f B i s p e c t r u m可以看出,不同发射机信号的HHT边际谱和双谱图像具有细微差异,这些细微区别是辐射源发射机硬件的不同引起的,是发射机固有特性在信号上的体现。因此HHT边际谱和双谱可以用于提取个体信号的特征,满足辐射源个体识别的基本要求。在HHT边际谱和双谱的特征分布散点图中,以不同颜色表示不同辐射源个体的特征。除个别异常值外,特征总体呈现明显的同类聚集性和类间分离性,具有良好区分度,能够反映不同的辐射源个体的差异。图1 8为4台发射机的四阶累积量切片谱。四阶累积量是一个有3个参数变量的四维函数,三谱即四阶谱,是四阶累积量的三维离散时间傅里叶变换,但计算量大且耗时长,在有限数据情况下相对于相关运算其估计方差较大。四阶累积量切片谱实质上是四阶谱在一维的投影,四阶累积量切片谱比四阶谱明显简化了计算的复杂度,但保留了其抑制高斯噪声的固有优点,利于实际应用中的图形显示。对基于HHT边际谱、边际谱信息熵、双谱投影、双谱矩形积分、四阶累积量切片谱的分类识别方法及本节所提方法的分类结果进行对比分析,具体的分类结果随S N R的变化特性曲线如图1 9所示。由图 1 9可得,对比基于HHT边际谱、边际谱信息熵、双谱投影、双谱矩形积分、四阶累积量切片谱的分类识别方法,本文提出的MF F方法有着更好的性能。MF F方法自S N R=9d B开始,分类正确率开始显著提高,自S N R=1 3d B图1 8 四阶累积量切片谱F i g.1 8 F o u r t h-o r d e r c u m u l a n t s l i c e s p e c t r u m图1 9 不同方法分类效果曲线F i g.1 9 C l a s s i f i c a t i o n e f f e c t c u r v e s o f d i f f e r e n t m e t h o d s开始,分类正确率的升高趋于平缓,S N R1 6 d B时,可以达到9 0%以上的分类正确率。在S N R=2 0 d B时,基于HHT边际谱方法的分类正确率为6 5.4 5%,基于边际谱信息熵方法的分类正确率为8 6.0 3%,基于双谱投影方法的分类正确率为7 9.9 8%,基于双谱矩形路径积分方法的分类正确率为7 0.1 5%,基于四阶累积量切片谱方法的分类正确率为8 3.7 7%,本文方法的正确率达到9 5.4 4%,同对比组中效果最好的基于边际谱信息熵方法相比,正确率有9.4 1%的提升。从分类效果曲线图中可看出,分类正确率主要受到特征提取和信噪比高低的影响。特征是决定辐射源个体识别效果的关键,特征提取的目的就是寻找表征显著类别差异的信息,提取出最能表征原始数据本质信息的少量综合特征。直接提王晓晗,等:基于多域融合的空间辐射源细微特征识别1 3 5 取信号变换的结果进行辐射源识别,存在大量的无效或冗余数据,并且单一域的特征难以完整体现信号的特性,这些都对辐射源识别带来困难。辐射源的特征十分细微,容易淹没在噪声中。噪声对识别正确率有着很大影响,不同方法的分类正确率都随着信噪比的增大而提高,在1 0 d B信噪比附近,各方法的正确率开始显著提高。5 结论本文分析了空间辐射源个体细微特征产生的机理,建立了发射端非线性模型,提取了辐射源信号的多域特征,实现基于MF F方法的辐射源设备细微特征的识别。本文主要完成的研究工作如下:1)分析了细微特征的产生机理,对谐波、交叉调制、相互调制、杂散输出等非线性产物做