一种重频抖动信号和常规信号分离方法与流程

本发明属于雷达数据处理及信号分选技术领域，特别涉及重频抖动信号和常规信号的分离方法；

背景技术：

信号分选是指根据辐射源信号测量参数形成的脉冲描述字(Pulse Descriptor Word，PDW)，通过脉冲串去交错生成辐射源描述字(Emitter Descriptor Word，EDW)；目前，基于脉冲重复间隔(Pulse repetition interval，PRI)的信号分选方法应用最为广泛。

重频抖动信号是复杂辐射源信号的一种重要表现形式，其时域特征可表示如下：t_n＝t_n-1+p(1±ε_n)n＝1,2,...,N-1，其中，p是抖动中心，ε是抖动半径，即抖动信号的PRI是在一定抖动范围内随机选取的值。

抖动信号的分选是一个难题，目前，利用改进的PRI变换法可以有效地分选抖动信号；然而，改进的PRI变换法只在输入信号全是抖动信号的情形下才具有良好的抖动信号分选能力，但是在信号类型较为混杂的情形下，典型如常规信号和抖动信号混合情形，受制于信号间的相互影响，常规信号的存在会降低抖动信号分选准确率；为此需要将常规信号的PDW和抖动信号的PDW进行分离，而实现上述分离的前提和难题，即如何获取输入PDW缓存中包含的潜在的常规信号PRI，从而根据PRI利用成熟的脉冲序列搜索方法提取出所有常规信号的PDW。

本发明的目的在于提供一种重频抖动信号和常规信号分离方法，针对常规信号和抖动信号混合情形，在生成的直方图中确定所有的上阶跃和下阶跃，经过排序和匹配后计算获取潜在的常规信号PRI，为实现常规信号和重频抖动信号的分离提供前提。

技术实现要素：

为实现上述技术目的，本发明采用技术方案包括以下步骤：

步骤1，设信号分选算法关注的PRI范围为[τ_min,τ_max]，将该范围分为K个均匀的小区间，称为PRI小箱，计算各个小箱的中心τ_k；

步骤2，遍历输入PDW缓存中所有脉冲对(t_n,t_m)，计算每个脉冲对的到达时间差，然后进行直方图统计；

步骤3，遍历统计完成的直方图，设置阶跃阈值为T，确定直方图中所有的上阶跃和下阶跃；分别对所有的上阶跃和下阶跃按照直方图从左至右的顺序进行排序，将具有相同序号的上阶跃和下阶跃进行一一配对；基于配对结果计算输入PDW缓存中所包含的潜在的常规信号PRI，其中每个配对结果可计算获得一个常规信号PRI。

本发明的特点和进一步改进在于：

在步骤1中，各个小箱的中心τ_k＝(k-1/2)b+τ_min，其中，k取1到K，b是每个PRI小箱的宽度：b＝(τ_max-τ_min)/K。

步骤2中，对于每个脉冲对(t_n,t_m)，计算其到达时间差|t_n-t_m|进行直方图统计时的小箱起点k_s和终点k_e，起点k_s＝max(k₁,0)，终点k_e＝min(k₂,K)，其中ε是信号分选算法支持分选的抖动信号的最大抖动半径。

步骤3具体子步骤为：

(3.1)遍历统计完成的直方图，设置阶跃阈值为T，确定直方图中所有的上阶跃和下阶跃：设直方图相邻小箱为b_i和b_i+1，i取1到K-1，若b_i+1高度大于b_i，且高度差大于阈值T，则(b_i,b_i+1)称为上阶跃；若b_i高度大于b_i+1，且高度差大于阈值T，(b_i,b_i+1)称为下阶跃；

(3.2)将所有的上阶跃按照直方图从左到右的顺序进行排序，同时将所有的下阶跃也按照直方图从左到右的顺序进行排序，将排序后具有相同序号的上阶跃与下阶跃进行一一配对；

(3.3)每个配对结果对应一个常规信号的PRI，计算方式如下：

设发生上阶跃的两个小箱的小箱中心分别为U₁、U₂，U₁<U₂，则存在某个常规信号，其PRI用p表示，满足：

设发生下阶跃的两个小箱的小箱中心分别为D₁、D₂，D₁<D₂，则存在某个常规信号，其PRI用p表示，满足：

计算得出同时满足上述两个条件的常规信号的PRI，若存在多个合理取值，则确保各个取值为相邻元素差为1的升序序列。

本发明的有益效果为：针对常规信号和抖动信号混合的情形，在生成的直方图中确定所有的上阶跃和下阶跃，经过排序和匹配后计算获取潜在的常规信号PRI，为实现常规信号和重频抖动信号的分离提供良好支持，有利于提高抖动信号分选准确率；

附图说明

图1为本发明的一种重频抖动信号和常规信号分离方法的流程图。

具体实施方式

图1为本发明的重频抖动信号和常规信号分离方法的流程图。结合附图，本发明的具体步骤为：

步骤1，设信号分选算法关注的PRI范围为[τ_min,τ_max]，将该范围分为K个均匀的小区间，称为PRI小箱，计算各个小箱的中心τ_k，其中各个小箱的中心τ_k＝(k-1/2)b+τ_min，k取1到K，b是每个PRI小箱的宽度：b＝(τ_max-τ_min)/K。

步骤2，遍历输入PDW缓存中所有脉冲对(t_n,t_m)，计算每个脉冲对的到达时间差，然后进行直方图统计；

(2.1)对于每个脉冲对(t_n,t_m)，计算其到达时间差|t_n-t_m|进行直方图统计时的小箱起点k_s和终点k_e，起点k_s＝max(k₁,0)，终点k_e＝min(k₂,K)，其中ε是信号分选算法支持分选的抖动信号的最大抖动半径。

步骤3，遍历统计完成的直方图，确定所有的上阶跃和下阶跃；分别对所有的上阶跃和下阶跃按照直方图从左至右的顺序进行排序，将具有相同序号的上阶跃和下阶跃进行一一配对；基于配对结果计算输入PDW缓存中所包含的潜在的常规信号PRI，其中每个配对结果可计算获得一个常规信号PRI：

(3.1)遍历统计完成的直方图，设置阶跃阈值为T，确定直方图中所有的上阶跃和下阶跃：设直方图相邻小箱为b_i和b_i+1，i取1到K-1，若b_i+1高度大于b_i，且高度差大于阈值T，则(b_i,b_i+1)称为上阶跃；若b_i高度大于b_i+1，且高度差大于阈值T，(b_i,b_i+1)称为下阶跃；

(3.2)将所有的上阶跃按照直方图从左到右的顺序进行排序，同时将所有的下阶跃也按照直方图从左到右的顺序进行排序，将排序后具有相同序号的上阶跃与下阶跃进行一一配对；

(3.3)每个配对结果对应一个常规信号的PRI，计算方式如下：

设发生上阶跃的两个小箱的小箱中心分别为U₁、U₂，U₁<U₂，则存在某个常规信号，其PRI用p表示，满足：

设发生下阶跃的两个小箱的小箱中心分别为D₁、D₂，D₁<D₂，则存在某个常规信号，其PRI用p表示，满足：

计算得出同时满足上述两个条件的常规信号的PRI，若存在多个合理取值，则确保各个取值为相邻元素差为1的升序序列。