CN113297969A - 一种雷达波形识别方法及*** - Google Patents
一种雷达波形识别方法及*** Download PDFInfo
- Publication number
- CN113297969A CN113297969A CN202110569180.5A CN202110569180A CN113297969A CN 113297969 A CN113297969 A CN 113297969A CN 202110569180 A CN202110569180 A CN 202110569180A CN 113297969 A CN113297969 A CN 113297969A
- Authority
- CN
- China
- Prior art keywords
- signal
- target component
- time
- neural network
- waveform identification
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013528 artificial neural network Methods 0.000 claims description 31
- 238000009826 distribution Methods 0.000 claims description 25
- 238000004422 calculation algorithm Methods 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000004927 fusion Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 abstract description 17
- 230000007547 defect Effects 0.000 abstract description 5
- 238000012549 training Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 6
- 210000002569 neuron Anatomy 0.000 description 6
- 238000000605 extraction Methods 0.000 description 5
- 230000002776 aggregation Effects 0.000 description 4
- 238000004220 aggregation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 101150056779 rbf-1 gene Proteins 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Signal Processing (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种雷达波形识别方法及***,能够处理类线性调频雷达信号波形识别的问题,首先通过对观测信号进行分数阶傅里叶变换,结合时频分析的识别方法,确定观测信号的WVD中目标分量主脊线条数及其倾斜角度与目标分量宽度的标准差,最后构建一个组合神经网络作为分类器,根据提取的观测信号时频特征参数识别不同雷达信号波形。能过克服传统的时频分析方法存在的弊端,有效提升平均识别率。
Description
技术领域
本发明涉及信号处理技术领域,特别是涉及一种雷达波形识别方法及***。
背景技术
雷达波形识别是根据接收到的雷达信号识别其脉内调制类型的一项技术。在现代战场环境下,高效率地识别敌方雷达波型,可以帮助己方对目标雷达类型、载体、用途、威胁等级具有充分的了解,进而做出下一步的作战决策。雷达波形识别作为电子侦察领域的一项关键技术,在现代信息化战场上发挥着越来越重要的作用。
早期的相关识别算法主要是基于时域,或通过傅里叶变换将其变换到频域进行特征分析,以及高阶统计分析的方法进行雷达波型识别。但是,随着雷达脉内调制方法的日趋复杂,上述传统方法已无法适应具有时变特性的非平稳雷达信号。
发明内容
本发明的目的是提供一种雷达波形识别方法及***,能够克服传统的时频分析法存在的弊端,有效提升平均识别率。
为实现上述目的,本发明提供了如下方案:
一种雷达波形识别方法,所述方法包括:
对观测信号进行分数阶傅里叶变换,提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
本发明还提供了一种雷达波形识别***,所述***包括:
参数计算模块,用于对观测信号进行分数阶傅里叶变换,并提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
波形识别模块,用于根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种雷达波形识别方法及***,考虑到分数阶傅里叶变换对于类线性调频信号具有很好的能量聚集性,利用基于分数阶傅里叶变换结合时频分析的识别方法实现对雷达波形的识别,能够克服传统时频分析方法存在的弊端,有效提升平均识别率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种雷达波形识别方法流程图;
图2为本发明实施例提供的类线性调频信号在分数阶域的能量聚集示意图;
图3为本发明实施例提供的6种类线性调频雷达信号识别准确率示意图;
图4为本发明实施例提供的总体识别准确率示意图;
图5为本发明实施例提供的多个信噪比下信号识别混淆矩阵示意图;
图6为本发明实施例提供的包含多个混合信噪比下的信号识别混淆矩阵示意图;
图7为本发明实施例提供的不同训练样本数目情况下的信号识别率示意图;
图8为本发明实施例提供的一种雷达波形识别***框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
随着时频分析技术的发展,在时域和频域联合分析雷达信号时频特征的方法逐渐成为了主流。现有技术中,存在以Wigner-Ville分布(Wigner-Ville distribution,WVD)作为时频分析的基础,同时为了克服其交叉项的影响,采用Choi-Williams分布(Choi-Williams distribution,CWD)对雷达信号波形进行时频分析,并在此基础上提取奇异值熵和分形维数作为时频特征,以此通过支持向量机(support vector machine,SVM)进行雷达波型的分类识别的方法。实验结果表明,对于8种雷达信号,在信噪比大于1dB时,平均识别率能达到95%以上。还有人采用CWD对雷达信号波形进行时频分析,并通过图像过滤、主成分分析(principal component analysis,PCA)、图像二值化算法、伪Zernike矩计算等方法针对时频图像提取时频特征,然后通过神经网络分类器完成雷达波型的分类识别。实验结果表明,对于8种雷达信号,在信噪比-2dB时,平均识别率能达到94.7%。
现有的识别算法来看,时频分析被广泛应用于雷达波型识别,根据雷达信号的时频分布设计特征提取算法提取时频特征,然后通过设计合理的分类器进行分类识别。但是对于一类特殊的雷达波形——类线性调频雷达信号波形(LFM、Frank码、P1码、P2码、P3码、P4码),其时频分布差异较小,利用传统的时频分析算法进行分类识别,误识别率较高。因此需要设计一种新的识别算法来处理这一类雷达信号。
本发明的目的是提供一种雷达波形识别方法及***,针对现有算法难以处理类线性调频雷达信号波形的问题,考虑到分数阶傅里叶变换对于类线性调频信号具有很好的能量聚集性,提出了一种基于分数阶傅里叶变换结合时频分析的识别算法,能够克服传统时频分析方法存在的弊端,有效提升平均识别率。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示,本实施例中提供了一种雷达波形识别方法,具体包括:
步骤101:对观测信号进行分数阶傅里叶变换预处理,并搜索信号能量在(α,u)平面上的峰值以及对应的最优分数阶变换阶数,根据预处理结果提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
步骤102:根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
本实施例提供的雷达波形识别方法能够有效应对类线性调频信号的波形识别问题。类线性调频雷达信号的时频分布较为相近,采用传统的时频分析方法难以适应其分类任务。类线性调频信号的数学模型表示为:
其中,A表示信号幅度,fc表示载频频率,φi表示调制相位。6种类线性调频信号的不同主要体现在调制相位φi上,其每一类信号的调制相位的数学模型如表1所示:
表1类线性调频雷达信号的调制相位的数学模型
针对类线性信号的特点,本实施例对观测信号进行分数阶傅里叶变换,提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差。
分数阶傅里叶变换定义为:
变换核函数Kp(u,t)定义如下:
分数阶傅里叶变换与Wigner-Ville分布的关系,具体如下:
信号x(t)的Wigner-Ville分布(WVD)定义为:
由分数阶傅里叶变换的时移性质以及逆变换可得:
代入(4)可得:
上式是在(t,ω)坐标系下的表达式,现将其逆时针旋转α角得到新坐标系(u,v),坐标变换关系为:
将式(7)代入式(6)可得:
式(8)表明分数阶傅里叶变换与WVD的关系为:对信号x(t)进行分数阶傅里叶变换得到Xp(u),Xp(u)的WVD是原信号x(t)的WVD的坐标逆时针旋转形式,旋转角度等于分数阶变换角α=pπ/2。
如图2所示,给出了类线性调频信号在分数阶域的能量聚集示意图,图中所示的类线性调频信号的WVD的目标分量为倾斜的主脊线,假设其与时域轴所成角度为β,根据分数阶傅里叶变换与WVD的关系,只要分数阶傅里叶变换的旋转角度α与β正交,信号就会在分数阶域的某个窄带范围内出现能量聚集。
然而,由于在实际应用中无法事先得知信号的WVD的倾斜角β,因此需要以旋转角α为变量,对观测信号进行分数阶傅里叶变换,得到信号在不同分数阶域的能量分布,从而形成信号能量在(α,u)平面上的二维分布。通过对此分布进行二维能量搜索,可得到能量峰值对应的最佳分数阶变换角α0=p0π/2,以及最佳分数阶域频率值u0:
以Frank码为例,其信号能量在(α,u)平面上的二维分布如图4所示。通过信号能量的二维搜索得到能量峰值以及对应的旋转角α0,此时旋转角α0与β正交。对信号进行旋转角为α0的分数阶傅里叶变换,得到对应的最佳分数阶域的能量分布。可得Frank码信号在最佳分数阶域存在两处能量聚集,对应于WVD中存在两条主脊线。因此根据信号的最佳分数阶域的能量聚集情况可判断出信号的WVD中主脊线的条数。同时,由图2可得,只有P2码信号对应的最佳旋转角其余信号对应的最佳旋转角上述所有特征可作为信号分类的依据。
因此,本实施例中在提取所述观测信号时频分布的目标分量主脊线条数参数时,以旋转角α为变量,对观测信号进行分数阶傅里叶变换,得到信号能量在(α,u)平面上的二维分布;
接着选取信号最佳分数阶域中符合预设能量阈值的位置点,本实施例中将预设能量阈值设置为0.6E,即选取信号最佳分数阶域的能量峰值E以及能量值大于0.6E对应的位置点u1,u2,…uk,计算其均值根据均值u0将位置点u1,u2,…uk划分为U1和U2两部分,其中U1包含小于u0的位置点,U2包含大于u0的位置点,分别计算对应均值u′、u″。
由此可将6种类线性调频雷达信号分为三类:LFM、P1码、P4码;Frank码、P3码;P2码。但还是无法区分LFM、P1码、P4码信号,以及Frank码、P3码信号。
由于Frank码、P1码和P2码信号是对LFM信号的频率步进逼近,因此其时频分布有明显的块状结构;P3码和P4码信号是对LFM信号进行采样演化而来,其时频分布更加平滑;LFM信号的时频分布最为平滑,不存在任何其他多余结构。这一特征可通过目标分量宽度的标准差体现。
因为目标分量为与时域轴成一定角度的倾斜直线,直接计算困难。如果旋转时频图像使目标分量与时域轴或频域轴平行,则目标分量的标准偏差的计算将变得更加容易。但是在旋转过程中的插值算法又会带来较大计算量。为解决这一问题,本发明提出一种简便算法,具体如下:
对所述观测信号进行特定变换阶数的分数阶傅里叶变换,使得变换后的信号时频分布的目标分量平行于时域轴,计算目标分量相对于时域轴在每一时域点处的宽度值;
根据所述目标分量相对于时域轴在每一时域点出的宽度值计算目标分量宽度的标准差。
从信号x(t)的WVD逆时针旋转角度θ(θ=pπ/2),得到其阶数为p的分数阶傅里叶变换Xp(u)的WVD。因此在某一特定的变换阶数p下,Xp(u)的WVD中的目标分量平行于时域轴,此时便于计算目标分量的宽度的标准差。因此需要确定特定变换阶数p。
根据图2的假设,如果对信号x(t)进行分数阶变换角度为β(β=pπ/2)的分数阶傅里叶变换,那么得到的Xp(u)的WVD的目标分量平行于时域轴。由于通过二维能量搜索可得到能量峰值对应的最佳分数阶变换角α0,而β与α0正交,因此可得到β:
特定的变换阶数p即可计算得到p=2β/π。
确定了变换阶数p并经过分数阶傅里叶变换得到的Xp(u),其WVD中的目标分量平行于时域轴。此时需要对Xp(u)的WVD进行二值化处理得到二维矩阵B(i,j),i,j=1,2,…,M,其行和为i=1,2,…,M。由于二值化处理后,目标值为1,背景值为0,因此D(i)表征了目标分量相对于时域轴在每一时域点处的宽度值。然后对D(i)进行归一化处理得到:
最后计算目标分量宽度的标准差:
为了得到更精确的分类结果,本实施例中选择组合神经网络作为分类器对雷达波形进行识别。神经网络具有自学习的能力,它可以通过训练分析得到输入数据与输出数据的内在规律,从而预测给定新的输入数据对应的输出值。考虑到上述提取的3类时频特征参数可以区分不同的类线性调频雷达信号,但这些特征参数需要有效地集成以提高整体识别性能。同时为了实现多特征融合并克服过拟合的问题,我们采用组合神经网络作为信号分类器。
组合神经网络包含***输入、单个神经网络(子网络)设计、组合结构与融合规则。输入层与输出层的神经元个数由特征参数的个数和待识别的信号类别决定,因此这里输入、输出的神经元个数分别为3、6。对于单个神经网络的设计,考虑到BP(Backpropagation)神经网络广泛应用于机器学习,并且已经被相关数学理论证明:3层BP神经网络可以任意精度拟合任一种非线性函数,具有将输入非线性映射到输出的强大能力。同时径向基(Radialbasis function,RBF)神经网络是一种可以自适应地确定径向神经元并具有快速收敛速度的局部逼近神经网络。因此这里我们采用不同的BP与RBF神经网络作为子网络,并通过并联组合的结构搭建组合神经网络。组合神经网络的性能很大程度上取决于融合规则。由于经过训练的神经网络的输出可被看成近似后验概率,因此这里我们采用加权后验概率的多数表决算法作为融合规则识别波形:
假设待分类的雷达信号波形为K,子网络的个数为N,输入的特征向量为X,则第n个子网络的第k个神经元输出可表示为:
Onk(X)=P(ck|X)+enk(X) (13)
其中,P(ck|X)表示特征向量X被分类为第k类波形的后验概率,enk(X)表示第n个子网络的第k个神经元输出误差。假设ωnk={ω1k,ω2k,…,ωNk}是第n个子网络的第k个神经元的输出权重,则每个分类子网络分类为同一类型波形的加权输出总和为:
Sk(X)=P(ck|X) (15)
由此,本实施例克服了传统的时频分析方法存在的弊端,有效提升了平均识别率。同时,本实施例提到的分数阶傅里叶变换对类线性调频信号具有很高的能量聚集性,在低信噪比下本发明算法依旧能实现很高的识别率。
同时,本实施例利用分数阶傅里叶变换与WVD的关系,提取信号的时频分布中目标分量主脊线条数以及倾斜角度这两类特征参数,避免了直接处理二维时频分布所带来的较高计算量。同时在计算目标分量宽度的标准差的过程中,通过特定阶数的分数阶傅里叶变换使目标分量平行于时域轴,极大降低了计算难度。因此可较大程度上降低运算复杂度。
为了更好的说明技术方案的效果,本实施例中针对文中研究的6种类线性调频雷达信号,采用美国海军研究生院P.E.Pace等研发的低截获概率雷达产生工具箱(LPIGtoolkit),通过Matlab仿真,每种信号在不同信噪比环境下产生1000个样本,根据本文提出的特征提取算法获取样本信号的特征参数,并送入神经网络进行训练与测试识别。所产生信号的具体参数如下:采样频率为12000Hz,载波频率在3000~3500Hz之间均匀分布。Frank码和P1码信号N取值4~8,P2码信号N取值6、8,P3码和P4码信号脉冲压缩比ρ取值为N×N,其中N∈[4,8],码元速率设为载波频率的0.25。在-6,-4,-2,…12dB的信噪比环境下产生上述信号样本。另外,组合神经网络结构如下:BP1(5-10-6)、BP2(5-12-6)、RBF1(5-10-6)、RBF2(5-12-6)。
实验结果如图3、图4所示,在-2dB信噪比环境下,6种信号的识别率基本达到90%以上(图3),总体识别率达到94.33%(图4)。这说明本文提出的算法在低信噪比环境下能够实现较高识别率,在识别准确性方面性能优良。
图5列举了多个信噪比(-6、-4、-2、0、2dB)下信号识别混淆矩阵,图6示出了包含多个混合信噪比下的信号识别混淆矩阵示意图,由图6可得,P2信号在不同信噪比环境下的识别率基本均能达到100%。这是因为通过基于分数阶傅里叶变换的特征提取算法可以直接将P2信号分选出来,而分数阶傅里叶变换对于类线性调频信号具有很好的能量聚集性,抗噪声干扰性强。这也体现出本文提出的特征提取算法的优越性。同时在低信噪比下,Frank码信号与P3码信号之间容易出现误识别,P1码信号与P4码信号之间容易出现误识别。例如在-6dB信噪比下,Frank码信号有30%的概率误识别为P3码信号,P3码信号有16%的概率误识别为Frank码信号;P1码信号有21%的概率误识别为P4码信号,P4码信号有19%的概率误识别为P1码信号。经分析,此误识别主要因为低信噪比下时频分布受噪声影响较大,非目标分量的成分占比增大,从而影响特征提取与信号识别。
为了进一步研究信噪比对雷达波形识别的影响,图6中示出了包含-6dB、-4dB、-2dB、0dB、2dB的混合信噪比下的信号识别混淆矩阵。在-6dB、-4dB、-2dB、0dB、2dB信噪比下分别产生180个信号,随机组合得到900个信号样本,提取特征参数作为训练集;同理,在相同信噪比下分别产生20个信号,随机组合得到100个信号样本,提取特征参数作为测试集。通过神经网络进行训练与测试。实验结果(图6)表明,信号总体识别率为87.67%,说明本发明提出的信号识别模型在非单一信噪比环境下进行训练与测试,仍能达到较高识别率,能较好适应混合信噪比环境下的信号识别。
为了研究不同训练样本数目对信号识别率的影响,图7中示出了不同训练样本数目情况下的信号识别率情况。其中,每种雷达信号的测试样本数量固定为100个,训练样本的数量由100个增加到900个,每次增加100个,分别在信噪比-6dB、-4dB、-2dB、0dB下进行实验。实验结果(图7)表明,随着训练样本数目的增加,不同信噪比下的信号识别率随之增加。当训练样本数目达到400个时,不同信噪比下的信号识别率增长趋于平稳。因此本实施例提出的算法不需大数据量的训练样本,最少在400个训练样本的情况下即可使分类器达到最佳识别状态。
实施例2
如图8所示,本实施例提供了一种雷达波形识别***,所述***包括:
参数计算模块M1,用于对观测信号进行分数阶傅里叶变换,并提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
波形识别模块M2,用于根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
本说明书中每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种雷达波形识别方法,其特征在于,所述方法包括:
对观测信号进行分数阶傅里叶变换,提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
3.根据权利要求1所述的一种雷达波形识别方法,其特征在于,提取所述观测信号的峰值能量对应的旋转角包括:
以旋转角为变量对所述观测信号进行分数阶傅里叶变换,得到所述观测信号的信号能量在(α,u)平面上的二维分布;
对所述信号能量在(α,u)平面上的二维分布进行搜索,得出峰值能量对应的旋转角。
4.根据权利要求3所述的一种雷达波形识别方法,其特征在于,提取所述观测信号的目标分量主脊线条数参数包括:
选取信号最佳分数阶域的能量分布中符合预设能量阈值的位置点;
计算所述位置点的均值,并根据所述均值将所述位置点以所述位置点的均值为界划分为两部分;
分别计算两部分的均值,并根据所述两部分的均值的比值判断目标分量主脊线条数参数。
5.根据权利要求4所述的一种雷达波形识别方法,其特征在于,所述预设能量阈值为大于0.6倍的峰值能量。
6.根据权利要求3所述的一种雷达波形识别方法,其特征在于,提取所述观测信号的目标分量宽度标准差包括:
对所述观测信号进行特定变换阶数的分数阶傅里叶变换,得到目标分量相对于时域轴在每一时域点处的宽度值;
根据所述目标分量相对于时域轴在每一时域点出的宽度值计算目标分量宽度的标准差。
7.根据权利要求6所述的一种雷达波形识别方法,其特征在于,根据所述观测信号的峰值能量对应的旋转角计算得到所述特定变换阶数。
8.根据权利要求1所述的一种雷达波形识别方法,其特征在于,所述组合神经网络包括BP神经网络与RBF神经网络,所述BP神经网络与所述RBF神经网络并联组合成为所述组合神经网络。
9.根据权利要求1所述的一种雷达波形识别方法,其特征在于,所述组合神经网络的融合规则采用加权后验概率的多数表决算法。
10.根据权利要求1所述的一种雷达波形识别***,其特征在于,所述***包括:
参数计算模块,用于对观测信号进行分数阶傅里叶变换,并提取所述观测信号的时频特征参数;所述时频特征参数包括峰值能量对应的旋转角、目标分量主脊线条数参数以及目标分量宽度标准差;
波形识别模块,用于根据所述时频特征参数,通过组合神经网络分类器对所述观测信号的波形进行识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110569180.5A CN113297969B (zh) | 2021-05-25 | 2021-05-25 | 一种雷达波形识别方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110569180.5A CN113297969B (zh) | 2021-05-25 | 2021-05-25 | 一种雷达波形识别方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113297969A true CN113297969A (zh) | 2021-08-24 |
CN113297969B CN113297969B (zh) | 2022-11-18 |
Family
ID=77324590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110569180.5A Active CN113297969B (zh) | 2021-05-25 | 2021-05-25 | 一种雷达波形识别方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113297969B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115659136A (zh) * | 2022-12-09 | 2023-01-31 | 中铁第一勘察设计院集团有限公司 | 基于神经网络的无线干扰信号波形识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064063A (zh) * | 2011-10-21 | 2013-04-24 | 中国人民解放军海军航空工程学院 | 基于cwd特征的多相编码雷达信号波形自动识别方法 |
CN103746722A (zh) * | 2014-01-02 | 2014-04-23 | 东南大学 | 一种跳频信号跳周期和起跳时间估计方法 |
CN104299611A (zh) * | 2014-09-28 | 2015-01-21 | 江南大学 | 基于时频脊线-Hough变换的汉语声调识别方法 |
CN110187313A (zh) * | 2019-05-31 | 2019-08-30 | 中国人民解放军战略支援部队信息工程大学 | 基于分数阶Fourier变换的雷达信号分选识别方法及装置 |
CN111190157A (zh) * | 2020-01-10 | 2020-05-22 | 中国地质大学(武汉) | 一种ipix雷达回波数据时频分析方法及*** |
CN112560803A (zh) * | 2021-01-22 | 2021-03-26 | 南京航空航天大学 | 基于时频分析与机器学习的雷达信号调制识别方法 |
-
2021
- 2021-05-25 CN CN202110569180.5A patent/CN113297969B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064063A (zh) * | 2011-10-21 | 2013-04-24 | 中国人民解放军海军航空工程学院 | 基于cwd特征的多相编码雷达信号波形自动识别方法 |
CN103746722A (zh) * | 2014-01-02 | 2014-04-23 | 东南大学 | 一种跳频信号跳周期和起跳时间估计方法 |
CN104299611A (zh) * | 2014-09-28 | 2015-01-21 | 江南大学 | 基于时频脊线-Hough变换的汉语声调识别方法 |
CN110187313A (zh) * | 2019-05-31 | 2019-08-30 | 中国人民解放军战略支援部队信息工程大学 | 基于分数阶Fourier变换的雷达信号分选识别方法及装置 |
CN111190157A (zh) * | 2020-01-10 | 2020-05-22 | 中国地质大学(武汉) | 一种ipix雷达回波数据时频分析方法及*** |
CN112560803A (zh) * | 2021-01-22 | 2021-03-26 | 南京航空航天大学 | 基于时频分析与机器学习的雷达信号调制识别方法 |
Non-Patent Citations (3)
Title |
---|
王爽: ""基于时频分析和模糊函数的LPI雷达波形识别算法研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
田润澜: ""低截获概率雷达信号特征分析技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
谢存祥 等: ""基于时频特征提取和残差神经网络的雷达信号识别"", 《***工程与电子技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115659136A (zh) * | 2022-12-09 | 2023-01-31 | 中铁第一勘察设计院集团有限公司 | 基于神经网络的无线干扰信号波形识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113297969B (zh) | 2022-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109254274B (zh) | 一种基于特征融合的雷达辐射源辨识方法 | |
CN108600135B (zh) | 一种信号调制方式的识别方法 | |
Wang et al. | Deep learning-based UAV detection in pulse-Doppler radar | |
CN107301432B (zh) | 基于时频分析的自适应辐射源调制识别方法 | |
Shao et al. | Deep fusion for radar jamming signal classification based on CNN | |
CN105259541B (zh) | 一种多站雷达抗有源欺骗式干扰的方法 | |
CN112560803A (zh) | 基于时频分析与机器学习的雷达信号调制识别方法 | |
Cain et al. | Convolutional neural networks for radar emitter classification | |
Tang et al. | Seamless tracking of apparent point and extended targets using Gaussian process PMHT | |
Nuhoglu et al. | Image segmentation for radar signal deinterleaving using deep learning | |
CN114580461A (zh) | 一种雷达干扰识别网络的训练方法和雷达干扰识别方法 | |
CN113297969B (zh) | 一种雷达波形识别方法及*** | |
Cunxiang et al. | Quasi-LFM radar waveform recognition based on fractional Fourier transform and time-frequency analysis | |
CN114117912A (zh) | 一种数据模型双驱动下的海杂波建模与抑制方法 | |
CN115661869A (zh) | 基于多维特征智能融合的射频信号指纹识别方法和*** | |
Huang et al. | Radar waveform recognition based on multiple autocorrelation images | |
CN104680536A (zh) | 利用改进的非局部均值算法对sar图像变化的检测方法 | |
CN115616503A (zh) | 一种基于卷积神经网络模型的雷达干扰信号类型识别方法 | |
Wang et al. | Detection and parameter estimation of frequency hopping signal based on the deep neural network | |
CN113298138B (zh) | 一种雷达辐射源个体识别方法及*** | |
Delamou et al. | Deep learning-based estimation for multitarget radar detection | |
CN107390164A (zh) | 一种水下分布式多源目标的连续跟踪方法 | |
CN113203990A (zh) | 一种基于相噪特性的干扰源个体识别方法 | |
Tao et al. | Detection and recognition of LPI radar signals using visibility graphs | |
CN113109760B (zh) | 一种基于组稀疏的多线谱联合doa估计和聚类方法及*** |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20221102 Address after: 264001 No. two, 188 Road, Shandong, Yantai Applicant after: NAVAL AVIATION UNIVERSITY OF PLA Address before: 264001 No. two, 188 Road, Shandong, Yantai Applicant before: School of Aeronautical foundation, Naval Aeronautical University PLA |
|
GR01 | Patent grant | ||
GR01 | Patent grant |