CN108268837B - 基于小波熵和混沌特性的辐射源指纹特征提取方法 - Google Patents
基于小波熵和混沌特性的辐射源指纹特征提取方法 Download PDFInfo
- Publication number
- CN108268837B CN108268837B CN201711494014.3A CN201711494014A CN108268837B CN 108268837 B CN108268837 B CN 108268837B CN 201711494014 A CN201711494014 A CN 201711494014A CN 108268837 B CN108268837 B CN 108268837B
- Authority
- CN
- China
- Prior art keywords
- radiation source
- entropy
- analysis
- calculating
- chaotic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Complex Calculations (AREA)
Abstract
基于小波熵和混沌特性的辐射源指纹特征提取方法,涉及辐射源个体识别领域。包括以下步骤:1)对辐射源信号进行小波熵特征提取,得到特征参数;2)对辐射源信号进行混沌特征分析,提取混沌特征参数,将小波熵特征向量和混沌特征向量结合,得到组合特征向量,输入特征分类器,实现对辐射源个体的识别。克服了目前常用的时频分析方法没有考虑辐射源信号非线性本质的缺陷,既发挥了小波包变换的强时频分辨力特性,提取信号的多尺度局部特征,更通过混沌分析这一非线性分析方法,考量信号整体的非线性情况,从而更精准地反映辐射源信号的特征,使提取出的特征参数具有更强的区分性。
Description
技术领域
本发明涉及辐射源个体识别领域,尤其是涉及一种基于小波熵和混沌特性的辐射源指纹特征提取方法。
背景技术
所谓通信信号的“指纹”,就是指通信辐射源个体的细微特征以信号为载体的表现,通常指的是由于辐射源个体之间的硬件设备差异,而在所发射出的信号上附加了可以区别于其他辐射源的个体特征。在现代化海战中,水下目标识别是先敌发现并有效地对敌进行水声对抗,先敌使用武器攻击,克敌制胜的前提。如果能在复杂的水下环境中有效提取出辐射源的指纹特征,就能将不同辐射源装备区分,通过对个体的分析识别,可以进一步实现通信结构、战略部署的判定,为我军相应的战斗计划部署提供关键依据,具有重要的现实意义。
目前常用的辐射源指纹特征提取方法主要对辐射源信号的时频域进行分析,其缺陷在于这类方法将辐射源信号近似成线性信号来处理,没有考虑信号的非线性本质,因此无法很好地反映其非线性特征。混沌是非线性***一种重要的运动状态,能够很好地表征非线性随机过程内在的特殊规律,混沌分析作为一种非线性的分析方法,能够更好地体现非线性信号的本质特征。
小波变换是一种有效的时频分析方法,能够精细地提取信号在不同尺度下的信息,而小波包变换相比于小波变换具有更强的时频分辨力。信息熵是度量信号不确定度的重要参数,将小波包变换和信息熵分析结合,既能反映信号的多尺度局部特性,又能体现信号整体的复杂程度。
发明内容
本发明的目的在于克服基于时频分析的辐射源指纹特征提取方法没有反映辐射源信号非线性本质的缺陷,提供一种基于小波熵和混沌特性的辐射源指纹特征提取方法。
本发明包括以下步骤:
1)对辐射源信号进行小波熵特征提取,得到特征参数;
2)对辐射源信号进行混沌特征分析,提取混沌特征参数,将小波熵特征向量和混沌特征向量结合,得到组合特征向量,输入特征分类器,实现对辐射源个体的识别。
在步骤1)中,所述对辐射源信号进行小波熵特征提取,得到特征参数的具体方法可为:首先对辐射源信号进行n层小波包分解,得到2n个子频带,计算各子频带的信息熵,得到2n维特征向量。
在步骤2)中,所述对辐射源信号进行混沌特征分析包括四个部分:关联维数分析、Kolmogorov熵分析、Lyapunov指数分析和Hurst指数分析,构成4维混沌特征向量;首先对接收的一维辐射源信号时间序列进行相空间重构,在对重构相空间分析的基础上计算关联维数、Kolmogorov熵和Lyapunov指数,Hurst指数则通过对辐射源信号时间序列采用重标极差方法求得。
本发明克服了目前常用的时频分析方法没有考虑辐射源信号非线性本质的缺陷,既发挥了小波包变换的强时频分辨力特性,提取信号的多尺度局部特征,更通过混沌分析这一非线性分析方法,考量信号整体的非线性情况,从而更精准地反映辐射源信号的特征,使提取出的特征参数具有更强的区分性。
附图说明
图1是本发明提供的基于小波熵和混沌特性的辐射源指纹特征提取方法框图。
图2是本发明提供的小波熵特征提取操作流程图。
图3是本发明实施例的辐射源信号经关联维数分析后的特征曲线图。
图4是本发明实施例的辐射源信号经Kolmogorov熵分析后的特征曲线图。
图5是本发明实施例的辐射源信号经Lyapunov指数分析后的特征曲线图。
图6是本发明实施例的辐射源信号经Hurst指数分析后的特征曲线图。
具体实施方式
为使本发明的目的、方案和优点更加清楚明白,以下实施例将结合附图对本发明作进一步的说明。
图1为根据发明实施方案示出的基于小波熵和混沌特性的辐射源指纹特征提取方法实现框图。对辐射源信号分别提取小波熵特征和混沌特征参数。其中小波熵特征的提取步骤为:对辐射源信号进行n层小波包分解,得到2n个子频带,计算各子频带的信息熵,得到2n维特征向量。混沌特征分析包括关联维数分析、Lyapunov指数分析、Kolmogorov熵分析和Hurst指数分析,提取出4维混沌特征向量。最后将小波熵特征参数和混沌特征参数联合组成的2n+4维特征向量输入分类器进行个体识别。
图2为小波熵特征提取操作流程图。首先对辐射源信号序列X(k)(k=1,2,...,N)进行n层小波包分解,得到2n个子频带重构信号xi(t)(i=0,1,2,...,2n-1,t=1,2,...,m),计算重构信号各离散点能量Ei(t)=|xi(t)|2,进而计算各子频带重构信号的总能量各子频带信息熵通过如下公式计算:
pi(t)为子频带信号各离散点能量占总能量的比例,Si即为所求的第i个子频带的信息熵。通过上述方法求得各子频带的信息熵,即构成2n维小波熵特征向量。
对辐射源信号进行混沌特征分析前,首先对信号进行相空间重构,重构后的高维空间信号包含了丰富的非线性信息。通过选取适当的延迟时间τ和嵌入维数m,将原始辐射源信号序列X(k)(k=1,2,...,N)重构成高维空间序列Xi=[x(i),x(i+1),...,x(i+(m-1)×τ)],i=1,2,...,N-(m-1)×τ。
对辐射源信号进行关联维数分析和Kolmogorov熵分析,通过计算关联积分,绘制关联积分曲线的方法来实现,步骤如下:
(1)对信号进行m维相空间重构后,设定临界距离r,计算相空间中任意两点(Xi,Xj)之间的距离||Xi,Xj||,若小于r,则保留该相点对。重复该步骤,统计距离小于r的相点对的个数,计算其与相点对总数之比,进而得到关联积分函数:
(2)通过以下公式得到关联维数函数D(m):
取不同的嵌入维数m对信号进行相空间重构,在每一种嵌入维数下,取不同的临界距离r,计算相应的关联积分C(r,m),绘制ln[C(r,m)]-lnr曲线。如图3所示,图中曲线从上至下依次表示嵌入维数m从1增加到10的关联积分曲线,随着m的增加,曲线斜率逐渐趋于稳定,此时对应的斜率就等于关联维数D。
(3)采用以下公式计算Kolmogorov熵,记为K:
首先设定嵌入维数m=1,设定初始临界距离rij,按照步骤(1)计算关联积分C(rij,m),不断减小临界距离rij的取值,直到C(rij,m)的值不随着rij减小而变化时,将此时的C(rij,m)值记为C(r,m)。递增m的值,重复上述步骤,计算得到C(r,m+1),并根据式(5)计算K的值,当K的值不再随着m递增而改变,此时的K值即为Kolmogorov熵的值,如图4所示。
最大Lyapunov指数的计算采用小数据量法来实现,步骤如下:
(1)对一维辐射源信号序列X(k)(k=1,2,...N)做FFT变换,计算平均周期P;
di(0)=min(||Xi-Xv||),(i,v=1,2,...,M;|i-v|>P) (6)
(4)对于每一个离散时间步j,计算所有最近邻点对的j步距离di(j)取对数后的平均值y(j):
其中q是非零di(j)的数目,Δt为原始辐射源信号序列的采样间隔。
(5)选取不同的嵌入维数m,重复步骤(2)~(4),绘制不同嵌入维数情况下的y(j)-j曲线。如图5所示为嵌入维数从1到5时的一簇y(j)-j曲线,各条曲线有一段近似平行的部分,对这部分曲线用最小二乘法作出回归直线,该直线的斜率就是最大Lyapunov指数,记为λ1。
Hurst指数的计算采用重标极差法,步骤如下:
(1)将一维辐射源信号序列X(k)(k=1,2,...,N)平均分割成s个长度为l的相邻子区间yi(i=1,2,...,s),N=l·s;第i个子区间的第u项记为yi,u(u=1,2,...,l)。
(2)计算各子区间yi的均值Ei、累积离差Zi、极差Ri和方差Si:
Ri=max{Zi}-min{Zi} (11)
(3)计算所有子区间的Ri/Si的均值,记为(R/S)l,有:
(4)取不同的子区间长度l,得到不同的(R/S)l,绘制lnl-ln(R/S)l曲线,用最小二乘法拟合该曲线斜率,即得到Hurst指数的值H,如图6所示。
以上内容是具体的实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施值局限于这些说明。
Claims (4)
1.基于小波熵和混沌特性的辐射源指纹特征提取方法,其特征在于包括以下步骤:
1)对辐射源信号进行小波熵特征提取,得到特征参数,具体方法为:首先对辐射源信号进行n层小波包分解,得到2n个子频带,计算各子频带的信息熵,得到2n维特征向量;
2)对辐射源信号进行混沌特征分析,提取混沌特征参数,将小波熵特征向量和混沌特征向量结合,得到组合特征向量,输入特征分类器,实现对辐射源个体的识别;所述对辐射源信号进行混沌特征分析包括四个部分:关联维数分析、Kolmogorov熵分析、Lyapunov指数分析和Hurst指数分析,构成4维混沌特征向量;首先对接收的一维辐射源信号时间序列进行相空间重构,在对重构相空间分析的基础上计算关联维数、Kolmogorov熵和Lyapunov指数,Hurst指数则通过对辐射源信号时间序列采用重标极差方法求得;
对辐射源信号进行混沌特征分析,首先对信号进行相空间重构,重构后的高维空间信号包含丰富的非线性信息;通过选取适当的延迟时间τ和嵌入维数m,将原始辐射源信号序列X(k)(k=1,2,...,N)重构成高维空间序列:
Xi=[x(i),x(i+1),...,x(i+(m-1)×τ)],i=1,2,...,N-(m-1)×τ;
对辐射源信号进行关联维数分析和Kolmogorov熵分析,通过计算关联积分,绘制关联积分曲线的方法,步骤如下:
(1)对信号进行m维相空间重构后,设定临界距离r,计算相空间中任意两点(Xi,Xj)之间的距离||Xi,Xj||,若小于r,则保留该相点对,重复该步骤,统计距离小于r的相点对的个数,计算其与相点对总数之比,进而得到关联积分函数:
(2)通过以下公式得到关联维数函数D(m):
取不同的嵌入维数m对信号进行相空间重构,在每一种嵌入维数下,取不同的临界距离r,计算相应的关联积分C(r,m),绘制ln[C(r,m)]-lnr曲线,曲线从上至下依次表示嵌入维数m从1增加到10的关联积分曲线,随着m的增加,曲线斜率逐渐趋于稳定,此时对应的斜率就等于关联维数D;
(3)采用以下公式计算Kolmogorov熵,记为K:
首先设定嵌入维数m=1,设定初始临界距离rij,按照步骤(1)计算关联积分C(rij,m),不断减小临界距离rij的取值,直到C(rij,m)的值不随着rij减小而变化时,将此时的C(rij,m)值记为C(r,m);递增m的值,重复上述步骤,计算得到C(r,m+1),根据计算K的值,当K的值不再随着m递增而改变,此时的K值即为Kolmogorov熵的值。
3.如权利要求1所述基于小波熵和混沌特性的辐射源指纹特征提取方法,其特征在于最大Lyapunov指数的计算采用小数据量法,步骤如下:
(1)对一维辐射源信号序列X(k)(k=1,2,...N)做FFT变换,计算平均周期P;
di(0)=min(||Xi-Xv||),(i,v=1,2,...,M;|i-v|>P)
(4)对于每一个离散时间步j,计算所有最近邻点对的j步距离di(j)取对数后的平均值y(j):
其中q是非零di(j)的数目,Δt为原始辐射源信号序列的采样间隔;
(5)选取不同的嵌入维数m,重复步骤(2)~(4),绘制不同嵌入维数情况下的y(j)-j曲线,嵌入维数从1到5时的一簇y(j)-j曲线,各条曲线有一段近似平行的部分,对这部分曲线用最小二乘法作出回归直线,该直线的斜率就是最大Lyapunov指数,记为λ1。
4.如权利要求1所述基于小波熵和混沌特性的辐射源指纹特征提取方法,其特征在于Hurst指数的计算采用重标极差法,步骤如下:
(1)将一维辐射源信号序列X(k)(k=1,2,...,N)平均分割成s个长度为l的相邻子区间yi(i=1,2,...,s),N=l·s;第i个子区间的第u项记为yi,u(u=1,2,...,l);
(2)计算各子区间yi的均值Ei、累积离差Zi、极差Ri和方差Si:
Ri=max{Zi}-min{Zi}
(3)计算所有子区间的Ri/Si的均值,记为(R/S)l,有:
(4)取不同的子区间长度l,得到不同的(R/S)l,绘制lnl-ln(R/S)l曲线,用最小二乘法拟合该曲线斜率,即得到Hurst指数的值H。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711494014.3A CN108268837B (zh) | 2017-12-31 | 2017-12-31 | 基于小波熵和混沌特性的辐射源指纹特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711494014.3A CN108268837B (zh) | 2017-12-31 | 2017-12-31 | 基于小波熵和混沌特性的辐射源指纹特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108268837A CN108268837A (zh) | 2018-07-10 |
CN108268837B true CN108268837B (zh) | 2020-07-03 |
Family
ID=62772932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711494014.3A Active CN108268837B (zh) | 2017-12-31 | 2017-12-31 | 基于小波熵和混沌特性的辐射源指纹特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108268837B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948731A (zh) * | 2019-03-29 | 2019-06-28 | 成都大学 | 一种通信电台个体识别方法、***、存储介质及终端 |
CN111695444B (zh) * | 2020-05-21 | 2023-06-23 | 哈尔滨工业大学 | 一种基于波原子变换的辐射源个体特征提取方法 |
CN112711026B (zh) * | 2020-12-19 | 2022-09-13 | 西北工业大学 | 基于变微分步长相对变化熵的目标检测与信号频率估计方法 |
CN113255545B (zh) * | 2021-06-03 | 2021-09-21 | 中国人民解放军国防科技大学 | 人工特征和深度特征相结合的通信辐射源个体识别方法 |
CN116257780A (zh) * | 2023-02-22 | 2023-06-13 | 中国人民解放军国防科技大学 | 基于盲信号分离的无监督特征提取和自适应的dbscan聚类方法 |
CN117851808A (zh) * | 2023-12-29 | 2024-04-09 | 中国人民解放军国防科技大学 | 一种基于差分相空间的鲁棒性非线性指纹特征提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650429A (zh) * | 2009-09-04 | 2010-02-17 | 西北工业大学 | 低信噪比水下目标信号的混沌特征参数提取方法 |
CN101732110A (zh) * | 2009-11-19 | 2010-06-16 | 杭州电子科技大学 | 用于手部动作识别的脑电和肌电信号混沌特征融合方法 |
CN103927556A (zh) * | 2014-05-07 | 2014-07-16 | 天津工业大学 | 一种基于小波包和近似熵的心电信号分类方法 |
CN106529478A (zh) * | 2016-11-08 | 2017-03-22 | 中国人民解放军空军工程大学 | 一种联合三维熵特征的雷达辐射源信号识别方法 |
CN106991393A (zh) * | 2017-03-30 | 2017-07-28 | 上海电机学院 | 基于多重分形维数的辐射源指纹特征提取方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010135539A1 (en) * | 2009-05-20 | 2010-11-25 | The Trustees Of The University Of Pennsylvania | Self-adaptive bio-signal and modulation device |
-
2017
- 2017-12-31 CN CN201711494014.3A patent/CN108268837B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650429A (zh) * | 2009-09-04 | 2010-02-17 | 西北工业大学 | 低信噪比水下目标信号的混沌特征参数提取方法 |
CN101732110A (zh) * | 2009-11-19 | 2010-06-16 | 杭州电子科技大学 | 用于手部动作识别的脑电和肌电信号混沌特征融合方法 |
CN103927556A (zh) * | 2014-05-07 | 2014-07-16 | 天津工业大学 | 一种基于小波包和近似熵的心电信号分类方法 |
CN106529478A (zh) * | 2016-11-08 | 2017-03-22 | 中国人民解放军空军工程大学 | 一种联合三维熵特征的雷达辐射源信号识别方法 |
CN106991393A (zh) * | 2017-03-30 | 2017-07-28 | 上海电机学院 | 基于多重分形维数的辐射源指纹特征提取方法 |
Non-Patent Citations (4)
Title |
---|
基于左右手运动想象的脑电信号的分类研究;周宏;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130315(第03期);正文第33-34页第4.3.2节 * |
基于非线性动力学特征指标的癫痫脑电信号分类、定位与预测方法研究;李冬梅;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20170915(第09期);摘要,正文第16页第2.2.1节,第18页第2.2.2节,第21页第2.2.5节,第31页第2.2.2节,第34页第2.2.4节,第44页第2.1,2.2节,第45页第3节 * |
李冬梅.基于非线性动力学特征指标的癫痫脑电信号分类、定位与预测方法研究.《中国优秀硕士学位论文全文数据库 医药卫生科技辑》.2017,(第09期), * |
硫化矿石碓自热过程的非线性及数值仿真研究;潘伟;《中国博士学位论文全文数据库 工程科技I辑》;20121215(第12期);正文第67-69页第4.3节,第70-71页,第74-75页第4.2.2节,附图4-3,附图4-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN108268837A (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108268837B (zh) | 基于小波熵和混沌特性的辐射源指纹特征提取方法 | |
Wong et al. | Clustering learned CNN features from raw I/Q data for emitter identification | |
Ma et al. | Automatic modulation classification under non-Gaussian noise: A deep residual learning approach | |
CN114201988A (zh) | 卫星导航复合干扰信号识别方法和*** | |
CN113259288B (zh) | 基于特征融合与轻量化混合模型的水声调制方式识别方法 | |
CN108549078B (zh) | 一种雷达脉冲信号跨信道合并及检测方法 | |
CN112749633B (zh) | 分离与重构的个体辐射源识别方法 | |
Kong et al. | Waveform recognition in multipath fading using autoencoder and CNN with Fourier synchrosqueezing transform | |
Cheng et al. | Deep learning network based spectrum sensing methods for OFDM systems | |
CN114821251B (zh) | 一种点云上采样网络的确定方法及确定装置 | |
Li et al. | Research on feature importance evaluation of wireless signal recognition based on decision tree algorithm in cognitive computing | |
Dong et al. | Radar signal modulation recognition based on bispectrum features and deep learning | |
Sun et al. | Deep learning based preamble detection and TOA estimation | |
Nguyen et al. | A parallel neural network-based scheme for radar emitter recognition | |
Guven et al. | Classifying LPI radar waveforms with time-frequency transformations using multi-stage CNN system | |
Ding et al. | Radio frequency fingerprint extraction based on singular values and singular vectors of time-frequency spectrum | |
CN110133628A (zh) | 一种短时强干扰下lfm信号的恢复方法 | |
Zheng et al. | Embedded Radio Frequency Fingerprint Recognition Based on A Lightweight Network | |
CN116383719A (zh) | 一种用于lfm雷达的mgf射频指纹识别方法 | |
Akinin et al. | Application of EM algorithm in problems of pattern recognition on satellite images | |
Zhou et al. | Specific emitter identification via feature extraction in Hilbert-Huang transform domain | |
Wei et al. | A multi-resolution channel structure learning estimation method of geometry-based stochastic model with multi-scene | |
Ng et al. | Effect of neighbourhood size selection in SOM-based image feature extraction | |
Guo-Xi | Research on a Convolutional Neural Network Method for Modulation Waveform Classification. | |
CN106650678A (zh) | Gabor小波子带相关结构人脸识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |