一种基于频谱的雷达信号调制识别方法与流程

1.本发明涉及雷达信号调制识别的技术领域，尤其涉及一种基于频谱的雷达信号调制识别方法。


背景技术：

2.雷达调制类型识别的目的是识别截获雷达信号的调制类型，为判断雷达辐射源设备功能、特性，协成波形干扰或抗干扰任务。随着雷达技术的快速发展，雷达使用数量和种类越来越多，雷达信号调制方式越来越复杂，电磁环境也变得更加恶劣。在这样的电子环境下，繁多的信号脉内调制形式、较大的信号噪声比(snr)范围，使得通过单个特征仍然无法提取到信号的本质特征，不仅让雷达信号的调制方式越来越多样化,而且使雷达信号正常工作的信噪比(snr)门限也越来越低,造成了一些只针对少数几种信号、抗噪性差的传统雷达脉内调制识别算法失效。因此准确地识别雷达信号调制方式是一项重要而艰巨的任务。
3.随着现代电子和雷达技术的迅猛发展，雷达信号波形的自动识别被广泛的应用到了认知雷达、雷达发射机识别以及威胁检测和分析等系统中，成为了电子对抗水平的重要衡量标志。因此，在电子战中对雷达信号进行高效准确的识别变得越来越重要。
4.目前，申请号为202210061559.x的中国发明专利公开了一种高效灵活的雷达频谱显示方法，能够根据频谱显示区像素宽度建立了相应的链表数据结构，并对采样点合理的抽象使雷达频谱显示。但是在军事领域雷达信号对雷达信号特征的提取时，待提取的信号频率是完全已知的，不适用于实际场景，会导致某些特征非常重要的雷达信号未被提取出来，识别的准确度低，而且由于军事领域雷达信号训练识别样本少，无法对雷达信号的调制方式进行精确的识别，识别的准确度骤降。


技术实现要素：

5.本发明解决的技术问题是：在军事领域雷达信号对雷达信号特征的提取时，待提取的信号频率是完全已知的，不适用于实际场景，会导致某些特征非常重要的雷达信号未被提取出来，而且由于军事领域雷达信号训练识别样本少，无法对雷达信号的调制方式进行精确的识别，识别的准确度骤降。
6.为解决上述技术问题，本发明提供如下技术方案：一种基于频谱的雷达信号调制识别方法，包括获取时域信号；对所述时域信号进行间隔取样，获取取样信息；计算所述取样信息之间的关联度；基于所述关联度对雷达信号调制识别。
7.作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：获取时域信号包括：对雷达信号进行分解，雷达信号包括回波信号和发射信号，雷达间隔发射雷达信号，按照雷达发射信号的时间对回波信号和发射信号进行采集。
8.作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：将雷达的回波信号与发射信号输入混频器进行混频，得到一段中频信号，将所述中频信号作为时域信号进行间隔取样。
9.作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：对所述时域信号进行间隔取样，获取取样信息包括：
10.其中，表示采样间隔，x(n)表示一组有限个信号采样值，x(k)表示一组有限个频谱采样值，nk表示截取长度，e表示自然对数的底数，i是虚数单位，n点表示频率区间[0,2n]上的间隔采样。
[0011]
作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：通过改变频率区间的距离，分别对不同的频率区间进行间隔采样，获取不同频率区间的取样信息。
[0012]
作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：计算取样信息之间的关联度包括：建立关联向量度模型，其计算表达式为：
[0013]
其中，c(q,r)为关联积分，q为奇异测度，r为关联积分的观测尺度，x={x1,x2,...x
n-1
},表示新构造的相空间序列，h(x)表示阶跃函数，x表示参考点的值；计算取样信息中每个测试点到参考点的距离，给定r的值，当相对距离小于r的值时，称其为关联向量对，然后对关联的向量对进行计数，得到的相对的统计距离。
[0014]
作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：由于关联积分和多重分形谱是近似对数的关系，对r的值进行计算的计算表达式为：
[0015]
其中，dq表示多重分形维数。
[0016]
作为本发明所述的基于频谱的雷达信号调制识别方法的一种优选方案，其中：基于所述关联度对雷达信号调制识别包括：基于关联向量度模型计算频谱各个波段雷达信号的关联度区间;根据关联积分所在的关联度区间对雷达信号进行识别。
[0017]
本发明的有益效果：能够根据雷达信号的时频特征，对相邻两个采样点的信号特征的关联度进行计算，对时频特征信息进行深入挖掘，在样本量较少的情况下，将雷达信号时频特征提取出来，提高对雷达信号调制识别能力。通过将雷达发射脉冲信号与回波信号一起输入一个混频器进行混频，得到一段中频信号，然后对取样信息中的关联度进行计算，根据关联积分所在的关联度区间对雷达信号进行识别，以识别不同的目标物信号，有利于
提高识别准确率。
附图说明
[0018]
图1为本发明一个实施例提供的一种基于频谱的雷达信号调制识别方法的基本流程示意图。
[0019]
图2为本发明一个实施例提供的一种基于频谱的雷达信号调制识别方法雷达信号4ghz时的输出频谱。
[0020]
图3为本发明一个实施例提供的一种基于频谱的雷达信号调制识别方法雷达信号2.1ghz时的输出频谱。
[0021]
图4为本发明一个实施例提供的一种基于频谱的雷达信号调制识别方法的频响特性示意图。
实施方式
[0022]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合说明书附图对本发明的具体实施方式做详细的说明，显然所描述的实施例是本发明的一部分实施例，而不是全部实施例。
实施例1
[0023]
参照图1，为本发明的一个实施例，提供了一种基于频谱的雷达信号调制识别方法，包括：s1：获取时域信号包括：对雷达信号进行分解，雷达信号包括回波信号和发射信号，雷达间隔发射雷达信号，按照雷达发射信号的时间对回波信号和发射信号进行采集。
[0024]
将雷达的回波信号与发射信号输入混频器进行混频，得到一段中频信号，将所述中频信号作为时域信号进行间隔取样。
[0025]
雷达信号识别通过分析不同调制波形信号脉内特征，得到目标物反射的特有属性，实现雷达调制类型识别。雷达在工作时间段内，其发射的雷达信号，可以看做一个个具有时间间隔的小工作时间段组成的，而在小工作时间段内雷达发射脉冲信号，雷达信号的频率在小工作时间段内具有周期性变化，通过将小工作时间段内雷达发射脉冲信号与回波信号一起输入一个混频器进行混频，得到一段中频信号，中频信号包含了目标物的信息，对中频信号后续再经过滤波、放大、模数转换和测频等处理后就可以获得目标物的信息，根据中频信号的时频特征信息凸显不同调制类型信号的局部差异，从而表征目标信号，所以将所述中频信号作为时域信号进行间隔取样。
[0026]
对所述时域信号进行间隔取样，获取取样信息包括：
[0027]
其中，表示采样间隔，x(n)表示一组有限个信号采样值，x(k)表示一组有限个
频谱采样值，nk表示截取长度，e表示自然对数的底数，i是虚数单位，n点表示频率区间[0,2n]上的间隔采样。
[0028]
通过改变频率区间的距离，分别对不同的频率区间进行间隔采样，获取不同频率区间的取样信息。
[0029]
可以对同一个小工作时间段对应的中频信号进行间隔采样，其采样间隔小于脉冲时长，或者对连续多个一个小工作时间段对应的中频信号进行间隔采样，其采用间隔大于脉冲信号发射间隔。
[0030]
能够根据雷达信号的时频特征，对相邻两个采样点的信号特征的关联度进行计算，对时频特征信息进行深入挖掘，在样本量较少的情况下，将雷达信号时频特征提取出来，提高对雷达信号调制识别能力。
[0031]
s3：计算所述取样信息之间的关联度包括：建立关联度模型，其计算表达式为：
[0032]
其中，c(q,r)为关联积分，q主要区分不同的度量空间测度，r为关联积分的观测尺度，r也可以看作是时序序列的时延，r也可以作为观察阈值，x={x1,x2,...x
n-1
},表示新构造的相空间序列，q在关联积分中表示奇异测度，在相空间中就代表的是嵌入的维数，h(x)表示阶跃函数，x表示参考点的值；结合多重分形的定义可知，关联积分和多重分形谱是近似对数的关系。
[0033]
计算取样信息中每个测试点到参考点的距离，给定r的值，当相对距离小于r的值时，称其为关联向量对，然后对关联的向量对进行计数，得到的相对的统计距离。
[0034]
关联积分原来常用于时间序列的计算之中，核心的算法思想是选取一个用于计算相对距离的参考点，然后计算序列中每个测试点到参考点的距离。然后对于给定的阈值r，当相对距离小于阈值r的值时，称其为关联向量对。然后对关联向量对进行计数，得到的相对的统计距离。这个距离就是我们所求的关联积分。
[0035]
若r取值过大会导致观测的尺度过大，很难对不同信号的特征进行区分。若r的取值过小则观测尺度过小，函数值受到环境噪声影响较大并且计算复杂度增大。因此阈值r的选取将影响多重分形维数的计算。在r的有效的阈值区间之内，关联向量对才有意义。
[0036]
由于关联积分和多重分形谱是近似对数的关系，对r的值进行计算的计算表达式为：
[0037]
其中，dq表示多重分形维数，通过不断调整r值，可以不断得到关联函数c与多重分形维数d之间的函数关系式。
[0038]
s4：基于所述关联度对雷达信号调制识别包括：
基于所述关联度对雷达信号调制识别包括：基于关联向量度模型计算频谱各个波段雷达信号的关联度区间;根据关联积分所在的关联度区间对雷达信号进行识别。
[0039]
由于雷达回波信号不仅包括发射信号的反射信号，还有杂波信号和干扰信号，导致回波信号频谱复杂，需要分析不同调制波形信号脉内特征，得到目标物反射信号的特有属性，实现雷达调制类型识别。可以预先对需要识别的目标物的雷达信号的关联度区间进行计算，确定其关联度区间；在雷达工作中，通过将雷达发射脉冲信号与回波信号一起输入一个混频器进行混频，得到一段中频信号，然后对取样信息中的关联度进行计算，根据关联积分所在的关联度区间对雷达信号进行识别，以识别不同的目标物信号，有利于提高识别准确率。
实施例2
[0040]
参照图2至图4，为本发明另一个实施例，该实施例不同于第一个实施例的是，提供了一种基于频谱的雷达信号调制识别方法的实验验证，为对本方法中采用的技术效果加以验证说明，本实施例采用传统技术方案与本发明方法进行对比测试，以科学论证的手段对比试验结果，以验证本方法所具有的真实效果。
[0041]
经实验测试采集雷达信号在4ghz时的输出频谱如图2所示，雷达信号在2.1ghz时的输出频谱如图3所示。
[0042]
在4ghz时，输出1.0ghz到1.4ghz的信号在500mhz带宽内杂散抑制度均大于55dbc，但是在输出频率为1.3ghz和1.4ghz时，主频对带外1200mhz-400mhz和1200mhz-500mhz处的抑制度只有约52dbc，而且在2.1ghz之外的杂散抑制只有5dbc左右，所以dds输出还需要滤波器对远端杂散进行抑制，低通滤波器xecjlfcn3216-1500对2.1g处抑制度优于20dbc，带通滤波器bpf004，对800mhz和2.1ghz的抑制度分别有20dbc和35dbc，其频响特性如图4所示。
[0043]
鉴相器选用成都振芯的gm4704b，通过adisimpll软件以鉴相频率为2.5mhz，本振锁相源的相位噪声约为-98dbc/hz@1khz、-100dbc/hz@10khz、-98dbc/hz@100khz，鉴相泄露杂散抑制度优于75dbc，跳频时间小于30us。经过二倍频至x波段后，相位噪声恶化6db，鉴相泄露杂散恶化6db，故最终输出x波段信号的相位噪声约为-92dbc/hz@1khz、-94dbc/hz@10khz、-92dbc/hz@100khz，鉴相泄露杂散抑制度优于69dbc。
[0044]
gm4912c脉冲响应时间由两部分时间累加而成，第一个是fpga对脉冲调制信号的响应延时约5个fpga工作周期（4.ghz的1/16作为fpga脉冲产生基准时钟）约35ns，第二个是gm4912c在线性调频模式下的脉冲响应延时为392个sysclk（4.ghz周期约0.208ns）时钟周期约82ns，这两个时间相加得到脉冲响应延时约为117ns。
[0045]
对雷达信号进行分解需要抑制其本振信号，需要抑制三阶交调信号即（lo+2if），根据混频器特性，中频信号输入功率为-8dbm时，其发射输出功率约(-8-9)=-17dbm，本振驱动功率约13dbm，发射端口泄露的本振信号功率约(13-40)=-27dbm，三阶交调信号(2if+lo)功率约-35dbm(实测值），所以混频后滤波器对低端本振泄露信号的抑制度要大于(55-(27-17))=45dbc、对高端三阶交调信号的抑制度要大于(55-(35-17))=37dbc才能满足55dbc的指标要求。
[0046]
因本振带宽较宽，为800mhz，方案中混频后采用了分段滤波，选用中电十三所的mems滤波器，采用两级滤波器级联的方式可以满足指标要求。混频并经一级mems滤波后，发射信号功率为-22dbm，第一级放大器增益为24db，输出p-1为9dbm，再经mems滤波器后第二级放大器增益为16db，输出p-1为11dbm，此时发射信号功率约为12dbm，然后经过10db定向耦合器，耦合端输出rfx_c，直通端再经过两级数控衰减器bw178d后输出。
[0047]
应当认识到，本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现，其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而，若需要，该程序可以以汇编或机器语言实现。在任何情况下，该语言可以是编译或解释的语言。此外，为此目的该程序能够在编程的专用集成电路上运行。
[0048]
应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

技术特征：
1.一种基于频谱的雷达信号调制识别方法，其特征在于，包括：获取时域信号；对所述时域信号进行间隔取样，获取取样信息；计算所述取样信息之间的关联度；基于所述关联度对雷达信号调制识别。2.如权利要求1所述的基于频谱的雷达信号调制识别方法，其特征在于：获取时域信号包括：对雷达信号进行分解，雷达信号包括回波信号和发射信号，雷达间隔发射雷达信号，按照雷达发射信号的时间对回波信号和发射信号进行采集。3.如权利要求2所述的基于频谱的雷达信号调制识别方法，其特征在于：将雷达的回波信号与发射信号输入混频器进行混频，得到一段中频信号，将所述中频信号作为时域信号进行间隔取样。4.如权利要求3所述的基于频谱的雷达信号调制识别方法，其特征在于：对所述时域信号进行间隔取样，获取取样信息包括：其中，表示采样间隔，x(n)表示一组有限个信号采样值，x(k)表示一组有限个频谱采样值，nk表示截取长度，e表示自然对数的底数，n点表示频率区间[0,2n]上的间隔采样。5.如权利要求4所述的基于频谱的雷达信号调制识别方法，其特征在于：通过改变频率区间的距离，分别对不同的频率区间进行间隔采样，获取不同频率区间的取样信息。6.如权利要求5所述的基于频谱的雷达信号调制识别方法，其特征在于：计算取样信息之间的关联度包括：建立关联向量度模型，其计算表达式为：其中，c(q,r)为关联积分，q为奇异测度，r为关联积分的观测尺度，x={x1,x2,...x
n-1
},表示新构造的相空间序列，h(x)表示阶跃函数；计算取样信息中每个测试点到参考点的距离，给定r的值，当相对距离小于r的值时，称其为关联向量对，然后对关联的向量对进行计数，得到的相对的统计距离。7.如权利要求6所述的基于频谱的雷达信号调制识别方法，其特征在于：r值计算表达式为：其中，d
q
表示多重分形维数。
8.如权利要求7所述的基于频谱的雷达信号调制识别方法，其特征在于：基于所述关联度对雷达信号调制识别包括：基于关联向量度模型计算频谱各个波段雷达信号的关联度区间；根据关联积分所在的关联度区间对雷达信号进行识别。

技术总结
本发明公开了一种基于频谱的雷达信号调制识别方法，涉及雷达信号调制识别的技术领域，包括获取时域信号；对所述时域信号进行间隔取样，获取取样信息；计算所述取样信息之间的关联度；基于所述关联度对雷达信号调制识别。获取时域信号包括：对雷达信号进行分解，雷达信号包括回波信号和发射信号，雷达间隔发射雷达信号，按照雷达发射信号的时间对回波信号和发射信号进行采集。能够根据雷达信号的时频特征，对相邻两个采样点的信号特征的关联度进行计算，对时频特征信息进行深入挖掘，在样本量较少的情况下，将雷达信号时频特征提取出来，提高对雷达信号调制识别能力。提高对雷达信号调制识别能力。提高对雷达信号调制识别能力。


技术研发人员：徐杰 王智 梁海兴 王俊 马德明
受保护的技术使用者：南京誉葆科技股份有限公司
技术研发日：2023.07.13
技术公布日：2023/8/14