CN106778610B - 一种基于时频图像特征的脉内调制识别方法 - Google Patents

一种基于时频图像特征的脉内调制识别方法 Download PDF

Info

Publication number
CN106778610B
CN106778610B CN201611164892.4A CN201611164892A CN106778610B CN 106778610 B CN106778610 B CN 106778610B CN 201611164892 A CN201611164892 A CN 201611164892A CN 106778610 B CN106778610 B CN 106778610B
Authority
CN
China
Prior art keywords
time
image
frequency
target
representing
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
Application number
CN201611164892.4A
Other languages
English (en)
Other versions
CN106778610A (zh
Inventor
高敬鹏
孔维宇
郜丽鹏
蒋伊琳
赵忠凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201611164892.4A priority Critical patent/CN106778610B/zh
Publication of CN106778610A publication Critical patent/CN106778610A/zh
Application granted granted Critical
Publication of CN106778610B publication Critical patent/CN106778610B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于雷达辐射源信号识别技术领域,具体涉及一种基于时频图像特征的脉内调制识别方法。本发明包括:(1)对雷达信号进行Wigner‑Vill分布时频变换,得到反映信号能量随时间和频率的时频分布图;(2)对步骤(1)得到的时频图像进行图像预处理;(3)对预处理后的图像进行插值打散原轨迹;(4)将步骤(3)中被打散成单个点的轨迹类比成目标跟踪中的目标位置,并引入联合概率数据关联来识别出不同雷达信号的轨迹;(5)运用中心矩提取时频图像的形状特征等。本发明提出了一种将雷达信号通过时频分布变换转换为时频图像,通过数字图像处理将同一时刻到达并且频率互相交叠的情况下,分别识别出雷达信号调制方式的方法。

Description

一种基于时频图像特征的脉内调制识别方法
技术领域
本发明属于雷达辐射源信号识别技术领域,具体涉及一种基于时频图像特征的脉内调制识别方法。
背景技术
雷达辐射源信号识别是电子情报侦察、电子支援侦察和威胁告警系统中的关键处理过程,其识别水平直接关系到雷达对抗设备的技术先进程度。随着现代电子战的激烈对抗,新型复杂体制雷达不断投入使用并逐渐占据主导地位,电磁威胁环境的信号密度高达每秒120万个脉冲以上,雷达工作频率覆盖范围达到了0.1~20GHz,并正在向0.05~140GHz扩展,雷达信号波形在时频等多个域中同时变化,隐身和抗干扰能力大大增强。传统五参数(载频、脉冲到达时间、脉冲幅度、脉冲宽度和脉冲到达方向)识别方法难以适应这样密集、复杂和多变的信号环境,雷达辐射源信号识别面临着前所未有的挑战。近年来关于雷达辐射源信号识别的大量研究成果表明,脉内调制识别将有望在新型复杂体制雷达辐射源信号识别技术与装备方面获得重要突破。现有的雷达辐射源信号脉内调制识别方法主要有小波及小波包变换法、小波脊线法、相像系数法、经验模式分解法、相对无模糊相位重构法、熵特征法、围线积分双谱法、复杂度特征法、分数傅里叶-包络法、分形维数法、模糊函数主脊切面法、瞬时频率派生特征法和模糊函数主脊切面特征法等。现有的方法为使信号之间差异明显,都选择将采样信号做某种变换,对变换后的信号进行特征提取,进而完成对信号的调制方式的识别。这些方法都主要针对少数几种信号,对于其它信号的识别效果还有待于进一步研究,同时这些方法都是只考虑了同一时刻只有单一信号的情况,没有探讨同一时刻多信号频率交叠的情况。
发明内容
本发明的目的是提出一种既能够对不同调制类型的雷达信号、又能够利用时频图像区分同一时刻多信号频率交叠的基于时频图像特征的脉内调制识别方法。
本发明的目的是这样实现的:
(1)对雷达信号进行Wigner-Vill分布时频变换,得到反映信号能量随时间和频率的时频分布图;
对相关函数作滑窗处理,得到时变的局部相关函数公式,
Figure BDA0001182201350000011
当窗函数取时间冲击函数,不加限制,而在时域取瞬时值,
Figure BDA0001182201350000021
对时变局部相关函数作Fourier变换,即可得到WVD,
Figure BDA0001182201350000022
设有n个分量信号,得到多分量信号的WVD,
Figure BDA0001182201350000023
雷达信号经过Wigner-Vill分布时频变换后得到时频图像;
(2)对步骤(1)得到的时频图像进行图像预处理;
对时频图像灰度值进行归一化;设
Figure BDA0001182201350000028
和σ2分别是所有像素点灰度值的均值和方差;
Figure BDA0001182201350000024
Figure BDA0001182201350000025
归一化后像素点的灰度值为:
Figure BDA0001182201350000026
采用自适应维纳滤波器初步去除时频图像的噪声点,对图像进行增强;采用形态学方法对灰度图像进行处理;形态学处理是应用具有一定形态的结构元素对灰度图像进行腐蚀和膨胀的操作,其中使用结构元素b对f的灰度膨胀记为
Figure BDA0001182201350000027
选择半径为3的菱型结构元素对时频图像进行闭运算;
通过检测信号时频分布的起止频率,将没有信号分布的图像区域剪切掉;
采用最近邻域插值法归一化时频图像的纵横比;
将原图定义为一个二维函数F(x,y),(x,y)点的函数值对应像素点的灰度值,f(x,y)为归一化纵横比后的目标图;
Fw、Fh设为原图的宽度和高度;fw、fh设为目标图的宽度和高度;原图坐标(x,y)和目标图坐标(x′,y′)之间公式为:
x=x′(Fw/fw)
y=y′(Fh/fh);
(3)对预处理后的图像进行插值打散原轨迹;
(4)将步骤(3)中被打散成单个点的轨迹类比成目标跟踪中的目标位置,并引入联合概率数据关联来识别出不同雷达信号的轨迹;
引入确认矩阵;确认矩阵定义为:
Ω=(ωjt),j=1,2,...,m(k),t=0,1,...,n
其中ωjt表示量测j是否落入目标t的确认门之内;t=0表示没有目标,此时对应的Ω的列的元素都是1;
设在杂波环境中已有T个目标,则他们的状态方程和测量方程分别表示为:
Figure BDA0001182201350000031
JPDA的状态更新方程为:
Figure BDA0001182201350000032
其中,
Figure BDA0001182201350000033
表示k时刻目标t的状态矢量,
Figure BDA0001182201350000034
表示目标t的状态预测矢量,V′(k)为联合新息:
Figure BDA0001182201350000035
其中,βjt表示目标t与观测j的关联概率,Vj′(k)表示滤波新息;
关联概率表示为:
Figure BDA0001182201350000036
其中,
Figure BDA0001182201350000037
表示在可行事件θj中,观测j是否源于目标t;如果源于目标则等于1,否则为0;在k时刻联合事件θ的条件概率表示为
Figure BDA0001182201350000038
其中,σt为目标检测指示器,如果目标在θ中与观测关联,则σt的值等于1,反之为零;Φ是虚假观测事件的个数,
Figure BDA0001182201350000039
表示目标t的检测概率,V表示航迹有效门体积;
(5)运用中心矩提取时频图像的形状特征;
设大小为M*N的二值图像表示为:f(x,y)∈{1,0},则图像的(p+q)阶原点矩定义为:
Figure BDA00011822013500000310
图像的(p+q)阶中心矩定义为:
Figure BDA0001182201350000041
其中
Figure BDA0001182201350000042
表示水平方向上的质心;
Figure BDA0001182201350000043
表示垂直方向的质心;
其中μ02表示图像在垂直方向上的伸展度;μ20表示图像在水平方向上的伸展度;μ11表示图像的倾斜度;μ03表示图像在垂直方向上的重心偏移度;μ12表示图像垂直伸展的均衡程度;由于时频图的差异性主要体现在垂直方向的频域上,故舍弃μ02,而构造出特征向量;
(6)选取支持向量机分类器的核函数,并根据训练样本的形态特征确定核函数的参数,以训练支持向量机分类器;
其中,支持向量机分类器的核函数为高斯核函数,高斯核函数为:
k(x1,x2)=exp(-γ*|x1-x2|2),
选择任意两类设计1个二分类SVM,总共设计n(n-1)/2个二分类SVM;然后将这个二分类SVM构成一个有向无循环图,该图有n(n-1)/2个非叶子节点和n个叶子节点,每个非叶子节点表示一个二分类SVM分类器,并与下一层的两个节点相连,每个叶子节点表示一个输出;
具体而言,选择上述的高斯核函数作为核函数,核函数参数的选择采用交叉验证的方法;其中,惩罚因子C从2-5成倍增加到25,核函数参数g从2-5成倍增加到25,对于(C,g),选择5折交叉验证,即将训练样本分成5份,挑选1份作为测试样本,剩下4份作为训练样本,重复5次,计算平均正确率,以此衡量当前(C,g)的性能,选择性能相对最好的(C,g)作为最终的参数;最后,将步骤(5)得到的特征向量输入训练好的支持向量机。
本发明的有益效果在于:提出了一种将雷达信号通过时频分布变换转换为时频图像,通过数字图像处理将同一时刻到达并且频率互相交叠的情况下,分别识别出雷达信号调制方式的方法。
附图说明
图1是本发明雷达信号调制识别的流程图。
图2是本发明采用的支持向量机分类器示意图。
具体实施方式
下面结合附图对本发明做进一步描述。
首先将待识别的雷达信号进行Wigner-Vill分布时频变换,得到时频分布图。时频分析反映了信号能量随时间和频率的分布,在时频域精确地描述了信号。通过数字图像处理方法从信号的时频分布中提取出有效的识别特征,并将特征向量输入到分类器完成分类识别任务。
本发明的具体步骤是:
1)对雷达信号进行Wigner-Vill分布时频变换,得到反映信号能量随时间和频率的时频分布图。
2)对时频分布图进行预处理,抑制噪声,将图形分割为单独的点。各个信号时频图像灰度值的动态范围是不一样的,像素点灰度值大的数据对分类识别有着大的影响,为了减少数据间的不平衡性,首先对时频图像灰度值进行归一化,然后使用自适应维纳滤波器去除时频图像的噪声点,对图像进行增强。为了减少后期图像处理的计算量和存储空间,对灰度图像进行二值化,将灰度图转化为二值图,二值化过程中阈值的选取采用一维最大熵法。接着对二值图像进行闭运算(膨胀之后再进行腐蚀运算),使得时频图像中信号分量的轮廓变得光滑,进一步减少噪声。从时频图中可以看出,并非所有区域都分布有信号,我们可以将没有信号分布的图像区域剪切掉,减小冗余信息,更有利于下一步信号特征的提取。最后采用最近邻插值法归一化时频图像的纵横比,使所有信号的时频图像大小都保持一致,并进一步减小数据量。
3)对步骤2预处理后的图像进行插值,将原来连成线的信号时频分布形状打散成单个一系列点迹。
4)通过联合概率数据关联算法对步骤3得到的图像进行轨迹识别,识别出几个轨迹代表可能包含几个不同的雷达信号,并将不同轨迹的信号分离。
5)采用图像处理的算法和工具对时频图像进行处理,提取时频图像的形状特征,本发明采用中心矩来提取形状特征。
6)采用支持向量机分类器识别雷达信号脉内调制方式。支持向量机是一种基于结构风险最小化的分类方法,通过核函数的非线性变换,将待分类样本从样本控件映射到高维特征控件,然后在特征控件中最大化分类间隔,确定最优分类超平面。从每种待识别的雷达信号中抽取一定的数量的样本,对其进行信号特征提取,并输入到对应的二分类SVM分类器中,通过计算,获得该SVM的参数值,直至最终确认整个分类器。将特征向量输入到支持向量机完成基于时频图像特征的脉内调制识别方法。
步骤1:对雷达信号进行Wigner-Vill分布时频变换,得到反映信号能量随时间和频率的时频分布图。
为了体现非平稳信号的局部时变特性,对相关函数作滑窗处理,得到时变的局部相关函数公式(1),
Figure BDA0001182201350000051
当窗函数取时间冲击函数,对不加限制,而在时域取瞬时值,此时如公式(2)
Figure BDA0001182201350000061
对时变局部相关函数作Fourier变换,即可得到WVD,其表达式如公式(3)
Figure BDA0001182201350000062
WVD的时间带宽积达到Heisenberg不确定原理给出的下界,因此具有最理想的时频分辨率。设有n个分量信号,可以得到多分量信号的WVD,如公式(4)
Figure BDA0001182201350000063
雷达信号经过Wigner-Vill分布时频变换后得到时频图像。
步骤2:对步骤1得到的时频图像进行图像预处理。
各个信号的时频图像灰度值的动态范围是不一样的,像素点灰度值大的数据对分类识别有着大的影响,为了减少数据间的不平衡性,需要对时频图像灰度值进行归一化。设
Figure BDA0001182201350000064
和σ2分别是所有像素点灰度值的均值和方差。
Figure BDA0001182201350000065
Figure BDA0001182201350000066
归一化后像素点的灰度值为:
Figure BDA0001182201350000067
接下来采用自适应维纳滤波器初步去除时频图像的噪声点,对图像进行增强。为了使得时频图像的信号分量轮廓变的光滑,降低噪声和交叉项时频点的影响,进一步采用形态学方法对灰度图像进行处理。形态学处理是应用具有一定形态的结构元素对灰度图像进行腐蚀和膨胀的操作,其中使用结构元素b对f的灰度膨胀记为
Figure BDA0001182201350000068
开运算可以去除比结构元素更小的时频独立点和细小的突出部分,而闭运算一般将狭窄的缺口连接起来,填充时频区域之间小的缝隙。考虑到二值化后信号项的时频区域会出现断裂的缝隙,本发明选择半径为3的菱型结构元素对时频图像进行闭运算。
时频图像剪切通过检测信号时频分布的起止频率,将没有信号分布的图像区域剪切掉,减小冗余信息,更有利于下一步信号特征的提取。
为了使所有信号的时频图像大小都保持一致,减小数据量,采用最近邻域插值法归一化时频图像的纵横比。将原图定义为一个二维函数F(x,y),(x,y)点的函数值对应像素点的灰度值,f(x,y)为归一化纵横比后的目标图。
Fw、Fh设为原图的宽度和高度。fw、fh设为目标图的宽度和高度。那么原图坐标(x,y)和目标图坐标(x′,y′)之间公式为:
Figure BDA0001182201350000076
坐标转化过程中,目标图的坐标值可能会出现小数,对此在原图像中选择距离该坐标点最近点灰度值作为目标图中的灰度值。
步骤3:对预处理后的图像进行插值以便打散原轨迹。
图像插值一般都是利用已知邻近像素点的灰度值来产生未知像素点的灰度值,以便由原始图像再生出具有更高分辨率的图像。在此处使用插值是直接将空白的像素点插入图像,以便获得被打散的轨迹。
步骤4:将步骤3中被打散成单个点的轨迹类比成目标跟踪中的目标位置,并引入联合概率数据关联来识别出不同雷达信号的轨迹。
为表示有效回波和各目标跟踪门的复杂关系,引入确认矩阵的概念。确认矩阵定义为:
Ω=(ωjt),j=1,2,...,m(k),t=0,1,...,n (9)
其中ωjt用以表示量测j是否落入目标t的确认门之内。t=0表示没有目标,此时对应的Ω的列的元素都是1,这是因为任意测量都可能源于杂波或虚警。量测落入跟踪门相交区域的情形,对应某些量测可能源于多个目标,联合概率数据关联的目的就是计算每个量测与其可能的各种源目标相关联的概率。
假设在杂波环境中已有T个目标,则他们的状态方程和测量方程分别表示为:
Figure BDA0001182201350000071
JPDA的状态更新方程为:
Figure BDA0001182201350000072
其中,
Figure BDA0001182201350000073
表示k时刻目标t的状态矢量,
Figure BDA0001182201350000074
表示目标t的状态预测矢量,V′(k)为联合新息:
Figure BDA0001182201350000075
其中,βjt表示目标t与观测j的关联概率,Vj′(k)表示滤波新息。
在JPDAF中,由于它考虑了所有目标和航迹的可行联合事件,关联概率可以表示为:
Figure BDA0001182201350000081
其中,
Figure BDA0001182201350000082
表示在可行事件θj中,观测j是否源于目标t。如果源于目标则等于1,否则为0。另外观测还必须服从两个假设:第一,每一个观测有唯一源;第二,对于一个给定的目标,最多只有一个观测以其为源。在k时刻联合事件θ的条件概率可表示为
Figure BDA0001182201350000083
其中,σt为目标检测指示器,如果目标在θ中与观测关联,则σt的值等于1,反之为零。Φ是虚假观测事件的个数,
Figure BDA0001182201350000084
表示目标t的检测概率,V表示航迹有效门体积。
步骤5:运用中心矩提取时频图像的形状特征。
设大小为M*N的二值图像可以表示为:f(x,y)∈{1,0},则图像的(p+q)阶原点矩可以定义为:
Figure BDA0001182201350000085
图像的(p+q)阶中心矩可以定义为:
Figure BDA0001182201350000086
其中
Figure BDA0001182201350000087
表示水平方向上的质心。
Figure BDA0001182201350000088
表示垂直方向的质心。
图像的不同阶数中心矩表征不同的物理意义,其中μ02表示图像在垂直方向上的伸展度;μ20表示图像在水平方向上的伸展度;μ11表示图像的倾斜度;μ03表示图像在垂直方向上的重心偏移度;μ12表示图像垂直伸展的均衡程度。由于时频图的差异性主要体现在垂直方向的频域上,故舍弃μ02,μ02这些描述水平方向时域的特征值,对此可以选用μ02,μ02,μ02,μ02,μ02作为时频图像特征。而构造出特征向量。
步骤6:选取支持向量机分类器的核函数,并根据训练样本的形态特征确定核函数的参数,以训练支持向量机分类器。
其中,支持向量机分类器的核函数为高斯核函数,高斯核函数为:
k(x1,x2)=exp(-γ*|x1-x2|2),其中,采用交叉验证的方法确定所述高斯核函数的参数。
对于n类雷达辐射源信号脉内调制方式识别问题,选择任意两类设计1个二分类SVM,这样,总共可设计n(n-1)/2个二分类SVM;然后将这n(n-1)/2个二分类SVM构成一个有向无循环图,如图2所示,该图有n(n-1)/2个非叶子节点和n个叶子节点,每个非叶子节点表示一个二分类SVM分类器,并与下一层的两个节点相连,每个叶子节点表示一个输出。
具体而言,选择上述的高斯核函数作为核函数,核函数参数的选择采用交叉验证的方法。其中,C(惩罚因子)从2-5~25(成倍增加),g(核函数参数γ)从2-5~25(成倍增加),对于某个(C,g),选择5折交叉验证,即将训练样本分成5份,挑选1份作为测试样本,剩下4份作为训练样本,重复5次,计算平均正确率,以此衡量当前(C,g)的性能,选择性能相对最好的(C,g)作为最终的参数。
最后,将步骤5得到的特征向量输入训练好的支持向量机,得到同一时刻多信号频率交叠的雷达信号脉内调制方式。

Claims (1)

1.一种基于时频图像特征的脉内调制识别方法,其特征在于,包括如下步骤:
(1)对雷达信号s(t)进行Wigner-Vill分布时频变换,其中t为时间,得到反映信号能量随时间和频率的时频分布图;
对相关函数作滑窗处理,得到时变的局部相关函数公式,
Figure FDA0002326066790000011
当窗函数取时间冲击函数δ(t),不加限制,而在时域取瞬时值,
Figure FDA0002326066790000012
对时变局部相关函数作Fourier变换,即可得到WVD,
Figure FDA0002326066790000013
设有a个单分量信号sε,ε=1,2,…,a,得到多分量信号的WVD,
Figure FDA0002326066790000014
雷达信号经过Wigner-Vill分布时频变换后得到时频图像;
(2)对步骤(1)得到的时频图像进行图像预处理;
对时频图像灰度值进行归一化;设
Figure FDA0002326066790000015
和σ2分别是所有像素点灰度值的均值和方差;
Figure FDA0002326066790000016
Figure FDA0002326066790000017
归一化后像素点的灰度值为:
Figure FDA0002326066790000018
采用自适应维纳滤波器初步去除时频图像的噪声点,对图像进行增强;采用形态学方法对灰度图像进行处理;形态学处理是应用具有一定形态的结构元素对灰度图像进行腐蚀和膨胀的操作,其中使用结构元素αb对输入元素F的灰度膨胀记为
Figure FDA0002326066790000019
选择半径为3的菱型结构元素对时频图像进行闭运算;
通过检测信号时频分布的起止频率,将没有信号分布的图像区域剪切掉;
采用最近邻域插值法归一化时频图像的纵横比;
将原图定义为一个二维函数F(x,y),(x,y)点的函数值对应像素点的灰度值,f(x′,y′)为归一化纵横比后的目标图;
Fw、Fh设为原图的宽度和高度;fw、fh设为目标图的宽度和高度;原图坐标(x,y)和目标图坐标(x′,y′)之间公式为:
Figure FDA0002326066790000021
(3)对预处理后的图像进行插值打散原轨迹;
(4)将步骤(3)中被打散成单个点的轨迹类比成目标跟踪中的目标位置,并引入联合概率数据关联来识别出不同雷达信号的轨迹;
引入确认矩阵;确认矩阵定义为:
Ω=(ωjb),j=1,2,…,mk;b=0,1,…,B
其中ωjb表示观测j是否落入目标b的确认门之内;b=0表示没有目标,此时对应的Ω的列的元素都是1;设在杂波环境中已有B个目标,则他们的状态方程和测量方程分别表示为:
Figure FDA0002326066790000022
JPDA的状态更新方程为:
Figure FDA0002326066790000023
其中,
Figure FDA0002326066790000024
表示k时刻目标b的状态矢量,
Figure FDA0002326066790000025
表示目标b的状态预测矢量,V′(k)为联合新息:
Figure FDA0002326066790000026
其中,βjb表示目标b与观测j的关联概率,V′j(k)表示滤波新息;
关联概率表示为:
Figure FDA0002326066790000027
其中,
Figure FDA0002326066790000028
表示在可行事件θj中,观测j是否源于目标b;如果源于目标则等于1,否则为0;在k时刻联合事件θ的条件概率表示为
Figure FDA0002326066790000031
其中,σb为目标检测指示器,如果目标在θ中与观测关联,则σb的值等于1,反之为零;Φ是虚假观测事件的个数,
Figure FDA0002326066790000032
表示目标b的检测概率,V表示航迹有效门体积;
(5)运用中心矩提取时频图像的形状特征;
设大小为M*N的二值目标图像表示为:f(x,y)∈{1,0},则图像的(p+q)阶原点矩mpq定义为:
Figure FDA0002326066790000033
图像的(p+q)阶中心矩μpq定义为:
Figure FDA0002326066790000034
其中p=1且q=0时原点矩为m10,p=0且q=0时原点矩为m00,则
Figure FDA0002326066790000035
表示水平方向上的质心;p=0且q=1时原点矩为m01,则
Figure FDA0002326066790000036
表示垂直方向的质心;
其中p=0且q=2时μ02表示图像在垂直方向上的伸展度;p=2且q=0时μ20表示图像在水平方向上的伸展度;p=1且q=1时μ11表示图像的倾斜度;p=0且q=3时μ03表示图像在垂直方向上的重心偏移度;p=1且q=2时μ12表示图像垂直伸展的均衡程度;由于时频图的差异性主要体现在垂直方向的频域上,故舍弃μ02,而构造出特征向量;
(6)选取支持向量机分类器的核函数,并根据训练样本的形态特征确定核函数的参数,以训练支持向量机分类器;
其中,支持向量机分类器的核函数为高斯核函数,高斯核函数为:h(x1,x2)=exp(-γ*|x1-x2|2)选择任意两类设计1个二分类SVM,总共设计n(n-1)/2个二分类SVM;然后将这个二分类SVM构成一个有向无循环图,有向无循环图有n(n-1)/2个非叶子节点和n个叶子节点,每个非叶子节点表示一个二分类SVM分类器,并与下一层的两个节点相连,每个叶子节点表示一个输出;
具体而言,选择上述的高斯核函数作为核函数,核函数参数的选择采用交叉验证的方法;其中,惩罚因子C从2-5成倍增加到25,核函数参数g从2-5成倍增加到25,对于(C,g),选择5折交叉验证,即将训练样本分成5份,挑选1份作为测试样本,剩下4份作为训练样本,重复5次,计算平均正确率,以此衡量当前(C,g)的性能,选择性能相对最好的(C,g)作为最终的参数;最后,将步骤(5)得到的特征向量输入训练好的支持向量机。
CN201611164892.4A 2016-12-16 2016-12-16 一种基于时频图像特征的脉内调制识别方法 Active CN106778610B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611164892.4A CN106778610B (zh) 2016-12-16 2016-12-16 一种基于时频图像特征的脉内调制识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611164892.4A CN106778610B (zh) 2016-12-16 2016-12-16 一种基于时频图像特征的脉内调制识别方法

Publications (2)

Publication Number Publication Date
CN106778610A CN106778610A (zh) 2017-05-31
CN106778610B true CN106778610B (zh) 2020-04-07

Family

ID=58892886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611164892.4A Active CN106778610B (zh) 2016-12-16 2016-12-16 一种基于时频图像特征的脉内调制识别方法

Country Status (1)

Country Link
CN (1) CN106778610B (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301432B (zh) * 2017-07-11 2020-09-25 哈尔滨工程大学 基于时频分析的自适应辐射源调制识别方法
CN108280394A (zh) * 2017-12-15 2018-07-13 西安电子科技大学 基于时频分布奇异值分解的雷达信号识别方法
CN110163040B (zh) * 2018-03-28 2022-12-30 信阳农林学院 非高斯杂波中雷达辐射源信号识别技术
CN108805039A (zh) * 2018-04-17 2018-11-13 哈尔滨工程大学 联合熵与预训练cnn提取时频图像特征的调制识别方法
CN108764331A (zh) * 2018-05-25 2018-11-06 哈尔滨工程大学 基于分数阶傅里叶变换的联合分类器多信号调制识别方法
CN109086749A (zh) * 2018-09-14 2018-12-25 汪新群 一种民用无人机系统被动式探测方法
CN109544479B (zh) * 2018-11-27 2021-06-15 中国科学院深圳先进技术研究院 一种时频图像的增强方法及相关设备
CN109684995A (zh) * 2018-12-22 2019-04-26 中国人民解放军战略支援部队信息工程大学 基于深度残差网络的特定辐射源识别方法及装置
CN110147812A (zh) * 2019-04-04 2019-08-20 中国人民解放军战略支援部队信息工程大学 基于扩张残差网络的雷达辐射源识别方法及装置
CN110244271B (zh) * 2019-05-17 2022-06-10 中国人民解放军战略支援部队信息工程大学 基于多重同步压缩变换的雷达辐射源分选识别方法及装置
CN110109063A (zh) * 2019-05-20 2019-08-09 吉林大学 一种基于深度信念网络的辐射源重频调制类型识别方法
CN111474524A (zh) * 2020-04-22 2020-07-31 安徽华可智能科技有限公司 一种雷达干扰装备干扰效果监测与决策支持系统
CN111427018A (zh) * 2020-04-22 2020-07-17 安徽华可智能科技有限公司 一种雷达干扰装备干扰效果评估方法
CN111679264B (zh) * 2020-06-19 2023-09-15 中国人民解放军空军工程大学 一种间歇采样转发干扰抑制方法、装置和设备
CN111898531B (zh) * 2020-07-30 2024-03-26 中国电子科技集团公司第三十六研究所 卫星通信信号识别方法、装置及电子设备
CN112115821B (zh) * 2020-09-04 2022-03-11 西北工业大学 一种基于小波近似系数熵的多信号智能调制模式识别方法
CN112712047B (zh) * 2021-01-08 2022-09-16 自然资源部第一海洋研究所 一种基于图像处理的海洋哺乳动物回声定位信号检测方法
CN112906591A (zh) * 2021-03-02 2021-06-04 中国人民解放军海军航空大学航空作战勤务学院 基于多级跳线残差网络的雷达辐射源识别方法
CN114036982B (zh) * 2021-11-04 2024-07-02 西安建筑科技大学 雷达辐射源个体识别方法、系统、设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6897802B1 (en) * 2003-11-10 2005-05-24 Raytheon Company Fusion of shape and multiscale features for unknown target rejection
CN101587186A (zh) * 2008-05-22 2009-11-25 赵力 一种雷达脉内调制信号的特征提取方法
CN102279390A (zh) * 2011-05-06 2011-12-14 西南交通大学 一种低信噪比雷达辐射源信号脉内调制识别方法
CN102645649A (zh) * 2012-05-14 2012-08-22 重庆大学 基于雷达目标距离像时频特征提取的雷达目标识别方法
CN103064063A (zh) * 2011-10-21 2013-04-24 中国人民解放军海军航空工程学院 基于cwd特征的多相编码雷达信号波形自动识别方法
CN104200214A (zh) * 2014-08-22 2014-12-10 中国人民解放军国防科学技术大学 使用Walsh变换的功率谱来描述图像局部特征区域的方法
CN105354592A (zh) * 2015-10-22 2016-02-24 上海无线电设备研究所 基于分类的最优时频分布设计与目标识别方法
CN105426842A (zh) * 2015-11-19 2016-03-23 浙江大学 基于支持向量机的表面肌电信号多类手部动作识别方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6897802B1 (en) * 2003-11-10 2005-05-24 Raytheon Company Fusion of shape and multiscale features for unknown target rejection
CN101587186A (zh) * 2008-05-22 2009-11-25 赵力 一种雷达脉内调制信号的特征提取方法
CN102279390A (zh) * 2011-05-06 2011-12-14 西南交通大学 一种低信噪比雷达辐射源信号脉内调制识别方法
CN103064063A (zh) * 2011-10-21 2013-04-24 中国人民解放军海军航空工程学院 基于cwd特征的多相编码雷达信号波形自动识别方法
CN102645649A (zh) * 2012-05-14 2012-08-22 重庆大学 基于雷达目标距离像时频特征提取的雷达目标识别方法
CN104200214A (zh) * 2014-08-22 2014-12-10 中国人民解放军国防科学技术大学 使用Walsh变换的功率谱来描述图像局部特征区域的方法
CN105354592A (zh) * 2015-10-22 2016-02-24 上海无线电设备研究所 基于分类的最优时频分布设计与目标识别方法
CN105426842A (zh) * 2015-11-19 2016-03-23 浙江大学 基于支持向量机的表面肌电信号多类手部动作识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Automatic Radar Waveform Recognition;Jarmo Lunden 等;《IEEE journal of selected topics in signal processing》;20070630;第1卷(第1期);第124-136页 *
Autonomous Time-Frequency Morphological Feature Extraction Algorithm for LPI Radar Modulation Classification;E.R.Zilberman 等;《2006 International Conference on Image Processing》;20070220;第2321-2324页 *
基于图像特征的雷达信号脉内调制识别方法;龚文斌 等;《电光与控制》;20080430;第15卷(第4期);第45-49页 *
基于视频分析的雷达信号调制方法识别技术研究;刘康宁;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150215;第2015年卷(第02期);第I136-945页 *

Also Published As

Publication number Publication date
CN106778610A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106778610B (zh) 一种基于时频图像特征的脉内调制识别方法
Wang et al. Deep learning-based UAV detection in pulse-Doppler radar
Shao et al. Deep fusion for radar jamming signal classification based on CNN
EP0888560B1 (en) Improved method of moment estimation and feature extraction for devices which measure spectra as a function of range or time
Malof et al. A large-scale multi-institutional evaluation of advanced discrimination algorithms for buried threat detection in ground penetrating radar
Lim et al. Detection and localization of people inside vehicle using impulse radio ultra-wideband radar sensor
CN107688170A (zh) 一种基于随机森林的雷达目标航迹起始方法
CN112130142B (zh) 一种复杂运动目标微多普勒特征提取方法及系统
CN103886606B (zh) 基于联合广义伽玛分布参数的sar图像分割方法
JP6773606B2 (ja) レーダ装置及びそのレーダ信号処理方法
Qiao et al. Human activity classification based on micro-Doppler signatures separation
Nuhoglu et al. Image segmentation for radar signal deinterleaving using deep learning
KR101234192B1 (ko) 바이스태틱 레이더의 송수신 각도 추출 장치
Xie et al. Fast ship detection from optical satellite images based on ship distribution probability analysis
Zhang et al. Target detection for HFSWR based on an S 3 D algorithm
Reichman et al. Improvements to the Histogram of Oriented Gradient (HOG) prescreener for buried threat detection in ground penetrating radar data
CN113960587B (zh) 基于类别信息反馈的毫米波雷达多目标跟踪方法
Sinha et al. Estimation of Doppler profile using multiparameter cost function method
Mostafa et al. Segmentation by classification for through-the-wall radar imaging using polarization signatures
Yan et al. Floating small target detection in sea clutter using mean spectral radius
Xia et al. Target detection in sea clutter via contrastive learning
Aftab et al. A Gaussian process convolution particle filter for multiple extended objects tracking with non-regular shapes
CN107729903A (zh) 基于区域概率统计和显著性分析的sar图像目标检测方法
Pelliza et al. Optimal Canny’s parameters regressions for coastal line detection in satellite-based SAR images
CN109271902B (zh) 复杂背景下基于时域经验模态分解的红外弱小目标检测方法

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