阅读说明：本技术 一种基于sdif与pri变换法结合的雷达信号分选方法 (Radar signal sorting method based on combination of SDIF and PRI transformation method ) 是由 张春杰 刘俞辰 司伟建 邓志安 曲志昱 侯长波 汲清波 杨梦� 于 2019-10-15 设计创作，主要内容包括：本发明属于电子对抗技术领域,具体涉及使搜索效果更快更明显的一种基于SDIF与PRI变换法结合的雷达信号分选方法。本方法包括如下步骤：对雷达信号进行预分组；利用SDIF方法对预分选后的雷达分组建立到达时间级差直方图；依据雷达信号模型时域特点快速搜索提取雷达信号；建立一级差直方图快速分析参差信号；判断是否仍有复杂类型雷达信号残留；查询各模块分选结果。本发明的有益效果在于：预分选与主分选结合,将SDIF与PRI变换法结合作为主分选,SDIF部分对复杂雷达电磁环境中常规信号、参差信号、脉间捷变频以及脉组捷变频信号进行快速而有效的分选,PRI变换法部分分选剩余的抖动信号,各部分算法各司其职,承上启下,组合成有效快速的综合分选算法。(The invention belongs to the technical field of electronic countermeasure, and particularly relates to a radar signal sorting method based on combination of SDIF and PRI transformation methods, which enables a search effect to be faster and more obvious. The method comprises the following steps: pre-grouping radar signals; establishing an arrival time level difference histogram for the pre-sorted radar groups by utilizing an SDIF (Serial digital interface) method; rapidly searching and extracting radar signals according to the time domain characteristics of the radar signal model; establishing a first-level difference histogram to quickly analyze the staggered signals; judging whether complex radar signal residues still exist or not; and inquiring the sorting result of each module. The invention has the beneficial effects that: the SDIF part is used for quickly and effectively sorting conventional signals, staggered signals, pulse-to-pulse frequency agility and pulse group frequency agility signals in a complex radar electromagnetic environment, the PRI conversion method part is used for sorting residual jittering signals, and each part of algorithm plays its own role and is started from the top to the bottom to form an effective and quick comprehensive sorting algorithm.)

一种基于SDIF与PRI变换法结合的雷达信号分选方法

技术领域

本发明属于电子对抗技术领域，具体涉及使搜索效果更快更明显的一种基于SDIF与PRI变换法结合的雷达信号分选方法。

背景技术

随着雷达技术的迅猛发展，战场电磁环境的复杂度进一步增加，战场电磁环境多变，信号种类繁多，雷达信号密度达到每秒百万甚至千万量级，传统预分选技术与单参数分选算法结合愈来愈难以应对复杂的电磁环境，在分选性能与分选速度上亟待改进。传统预分选技术主要是利用雷达信号脉内描述字中的参数，配合容差进行简易分组，易将一部捷变频雷达信号错分为多部雷达信号，并且受设置的PRI分选范围限制，若预设分选范围较大，大PRI信号与小PRI信号易相互干扰；若预设分选范围较小，则会错失有价值的雷达信号。常用的单参数分选算法有扩展关联法、差值直方图法、累积差值直方图法(CDIF)、顺序差值直方图法(SDIF)以及PRI变换法。

扩展关联法是单参数信号分选中最早应用的算法，同时也是最基础的，适用于早期简单的雷达电磁环境；CDIF法是现在普遍应用的一种雷达信号分选算法，基于到达时间差原理，在统计直方图基础上做出较大改进，累积不同级的到达时间差值，设定合理的检测门限，抑制谐波；SDIF是CDIF的改进版本，不再累积到达时间差，计算量相较CDIF减少一半，但检测门限较难设定，且无法分选抖动雷达信号；原始的序列搜索算法时间复杂度高且搜索参差信号困难；原始的参差检验流程复杂；PRI变换法可分选部分抖动信号，抗子谐波能力强，但难以分选参差信号，且运算量巨大，难以实际应用。有学者将改变时间起点与交叠箱的思想引入PRI变换法，使其具备分选大抖动信号能力，但PRI估计精度不佳，易造成错误分选。在实际雷达信号分选系统中，如何综合利用各算法，在提升分选速度的同时，提高分选效果准确性是当前亟待解决的问题。

发明内容

本发明的目的在于提供一种基于SDIF与PRI变换法结合的雷达信号分选方法。

本发明的目的是这样实现的：

一种基于SDIF与PRI变换法结合的雷达信号分选方法，包括如下步骤：

步骤1：设置一长一短两段观察时间，1000ms和10ms，利用雷达信号的脉冲描述字如到达方向DOA、载频CF、脉宽PW等参数，建立各参数直方图，对雷达信号进行预分组；

步骤2：根据步骤1预分选结果，利用SDIF方法对预分选后的雷达分组建立到达时间级差直方图，初始级差C为1，当直方图峰值超过门限值时，则判断该到达时间差值为可能的雷达信号PRI；

步骤3：按照步骤2中估计得到的PRI进行序列搜索，依据雷达信号模型时域特点快速搜索提取雷达信号，并且实现一次提取一部完整的参差雷达信号；

步骤4：若搜索到的PRI中存在相同值，则初步认为存在参差信号，将参差信号脉冲串依时域顺序还原，建立一级差直方图快速分析参差信号；

步骤5：判断经过提取的雷达信号电磁空间内是否仍有复杂类型雷达信号残留，如大范围抖动雷达信号，若存在，则采用多级PRI箱结构，如5％，10％，15％，对分组内剩余脉冲序列进行离散形式PRI变换；

步骤6：查询各模块分选结果，完成雷达信号信息输出，如雷达信号类型，PRI等，实现分选。

步骤1具体包括：去除噪声点；利用脉宽、到达角度等参数初步预分选；区分捷变频信号以及对其他类型信号预分选。

步骤2中当超出门限的峰值个数不唯一时，即同时有多个PRI超过门限，先按对应PRI的一半进行脉冲尝试搜索，若均搜索失败，则直接进行下一级差值直方图，再检测是否有PRI值超过门限。

步骤2中门限采用下述公式计算确定：

T_threshold＝α(E-C)e^-λτ

其中λ为单位时间内雷达信号的平均出现概率，E为当前预分组总体雷达脉冲个数，C为正在进行差值级数，(E-C)为脉冲间隔个数，α与脉冲丢失率相关，人为设定，τ为总时间间隔，即脉冲采集时间。

步骤5中PRI交叠箱宽度由基本箱宽度与箱中心抖动值二者之间的较大值决定，基本箱宽度b＝(PRI_max-PRI_min)/K，箱中心抖动值由箱中心位置k_L与箱抖动率ε_i共同决定，其中第L个箱中心位置k_L＝(L-1/2)·(PRI_max-PRI_min)/K+PRI_min，则以箱抖动率ε_i构建的第L个交叠箱宽度为b_iL＝max([b,2ε_ik_iL])。

步骤5中门限设置应考虑观察时间原则、消除子谐波原则以及消除噪声原则，针对第L个位置交叠箱的门限设置应选取三者中的最大值，则第L个交叠箱门限设置

其中α,β,γ为人为设定(如α＝0.34，β＝0.17，γ＝3)，与交叠箱抖动率个数相关，N为当前环境雷达脉冲总数。

本发明的有益效果在于：

1.预分选与主分选结合，将SDIF与PRI变换法结合作为主分选，SDIF部分对复杂雷达电磁环境中常规信号、参差信号、脉间捷变频以及脉组捷变频信号进行快速而有效的分选，PRI变换法部分分选剩余的抖动信号，各部分算法各司其职，承上启下，组合成有效快速的综合分选算法；

2.将统计直方图与多层次分选思想引入传统预分选，利用载频统计直方图实现对捷变频信号的分类，防止增批现象的发生，采用长短观察时间以及两级门限，防止大小范围PRI信号相互干扰；

3.利用序列差直方图过程中的级差直方图对正在抽取的信号脉冲进行初步参差信号判定，在遍历所有脉冲到达时间是否存在预测到达时间范围内时，结合雷达信号序列脉冲模型时域特点，大大减少遍历时间复杂度，使其具有对多重频参差雷达信号快速搜索的能力，提高搜索速度与搜索效果；

4.将统计直方图思想引入参差检验，通过还原原始重频参差雷达信号，利用统计直方图快速得到子周期信息，低复杂度，时间复杂度与空间复杂度均为O(n)，提升参差检验实时性；

5.首次提出多级PRI箱结构应用于PRI变换法，使抖动信号尽可能进入与其抖动率匹配的PRI箱中，弥补了单一PRI箱估计精度与估计范围不能兼得的缺点，在保证分选大范围抖动信号的同时，提高PRI估计精度。

附图说明

图1是本发明的整体分选算法流程图；

图2是本发明的预分选算法流程图；

图3是本发明的SDIF部分算法流程图；

图4是本发明实施例中的PRI变换法部分流程图；

图5是复杂电磁环境参数及分选结果图。

具体实施方式

下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

本发明实施例提供的基于SDIF与PRI变换法结合的雷达信号分选方法，整体流程如图1所示，首先进行预分选，预分选流程如图2所示，大大降低信号密度，减轻后续主分选压力，再利用SDIF算法在预分选结果中估计提取出普通重频信号，参差信号以及捷变频信号，并进行参差检验，SDIF部分算法流程如图3所示，如果预分组内仍存在抖动信号，则利用PRI变换法对其分选，PRI变换法部分流程如图4所示。

本发明步骤如下：

步骤1：设置一长一短两段观察时间(如1000ms和10ms)，利用雷达信号的脉冲描述字如到达方向DOA、载频CF、脉宽PW等参数，建立各参数直方图，对雷达信号进行预分组，实现捷变频信号与其余类型雷达信号的区分、以及各雷达信号的预分选；

步骤1-1：去除噪声点。噪声点的载频与脉宽等雷达信号参数均是随机分布的，按容差分别建立载频统计直方图与脉宽统计直方图，根据两直方图纵轴高度，设置门限分别清除孤立噪声点与干扰点，如载频统计直方图容差为10MHz，脉宽统计直方图容差为0.5μs，将高度低于10的脉冲点作为噪声点，在针对短观察时间预分选内将高度低于150的点作为干扰点，在针对长观察时间预分选内将高度高于2000的点作为干扰点，从脉冲描述字中删除；

步骤1-2：利用脉宽、到达角度等参数初步预分选。根据去噪后的脉宽、到达角度统计直方图，结合容差，依据脉宽、到达角度等参数将相似信号脉冲点分组，初步稀释信号。容差的设置应尽可能囊括信号脉宽、到达角度的抖动范围，以防止初步预分选造成错分，如设置脉宽容差为3μs，到达角度容差为0.5°；

步骤1-3：区分捷变频信号以及对其他类型信号预分选。在脉宽、到达角度预分选结果中利用载频统计直方图，设置门限将捷变频信号区分开来，将低于门限的信号判断为捷变频信号，将高于门限的信号依据载频容差划分为多个不同载频的其他雷达信号分组，如设置载频容差为100MHz，依据观察时间设置门限为：

T_threshold＝βT/τ β∈(0,1)

其中T为观察时间，τ为当前分选PRI范围中值，预先设置长短观察时间内的两种PRI范围，如短观察时间内PRI范围为0～200μs，长观察时间内PRI范围为200～2000μs，β为人为设定参数。

步骤2：根据步骤1预分选结果，利用SDIF方法对预分选后的雷达分组建立到达时间级差直方图(初始级差C为1，最高级差为8)，当直方图峰值超过门限值时，则判断该到达时间差值为可能的雷达信号PRI，当超出门限的峰值个数不唯一时，即同时有多个PRI超过门限，先按对应PRI的一半进行脉冲尝试搜索，即按此PRI搜索，是否可以搜索到5个脉冲，若均搜索失败，则直接进行下一级差值直方图，再检测是否有PRI值超过门限。

门限采用下述公式计算确定：

T_threshold＝α(E-C)e^-λτ

其中λ为单位时间内雷达信号的平均出现概率，E为当前预分组总体雷达脉冲个数，C为正在进行差值级数，(E-C)为脉冲间隔个数，α与脉冲丢失率相关，人为设定，τ为总时间间隔，即脉冲采集时间。

步骤3：按照步骤2中估计得到的PRI进行序列搜索，依据雷达信号模型时域特点快速搜索提取雷达信号，并且实现一次提取一部完整的参差雷达信号；

步骤3-1：依据步骤2估计PRI对应的级差直方图，在直方图中搜索符合估计PRI的脉冲个数；

步骤3-2：确定脉冲描述字内所估计PRI的起点脉冲位置；

步骤3-3：进行序列搜索，在搜索过程中，将本次搜索到的脉冲位置存储，由于下次预测脉冲位置必在本次搜索到的脉冲位置之后，在下次搜索时，直接从所存储位置开始搜索；

步骤3-4：在本次搜索结束之后，将本次搜索提取到的脉冲个数与级差直方图中寻找到的符合PRI个数比较，依据二者倍数关系初步判断是否为参差信号，若为参差信号，则在起点脉冲的下一个脉冲再次进行搜索。

步骤4：若搜索到的PRI中存在相同值，则初步认为存在参差信号，将参差信号脉冲串依时域顺序还原，建立一级差直方图快速分析参差信号；

步骤4-1：依据步骤2得到的PRI估计值，在一定容差内，根据相同的PRI值个数，初步判断参差数，如容差为1μs；

步骤4-2：利用序列搜索提取到的脉冲流，根据到达时间先后，还原出原始的参差雷达信号；

步骤4-3：对原始参差雷达信号依次进行一级差，并进行直方图统计；

步骤4-4：依据参差数，若参差数为3，依次选取频数较高的3个到达时间差值作为子周期。

步骤5：判断经过提取的雷达信号电磁空间内是否仍有复杂类型雷达信号残留，如大范围抖动雷达信号，若存在，则采用多级PRI箱结构(如5％，10％，15％)对分组内剩余脉冲序列进行离散形式PRI变换；

步骤5-1：判断当前雷达信号分组电磁空间经过SDIF分选提取后，空间内脉冲数是否超过门限(如设置为20)，若超过则进行PRI变换法分选，否则跳过步骤5；

步骤5-2：分别创建K个由多个不同抖动率构建的交叠PRI箱，如多个抖动率依次为ε₁、ε₂、ε₃(ε₁<ε₂<ε₃)，如K取300，分别创建300个5％、10％以及15％的PRI交叠箱，初始化PRI变换数组、标记脉冲流数组、时间起点数组以及标记时间起点数组，数组长度均为K；

步骤5-3：令n＝2，m＝1，τ＝t_n-t_m，判断到达时间差τ是否在预设的PRI范围内，若不在，则m＝m+1再进行判断，直至m＝n-1，若以上到达时间差均不在预设PRI范围内，令n＝n+1继续进行PRI范围判断；

步骤5-4：若到达时间差τ在预设的PRI范围内，则从抖动率为ε₁的第一个交叠箱开始，遍历各箱；

步骤5-5：首先判断到达时间差τ是否在第L个以ε₁抖动率构建的交叠箱范围内，若在其中，则将标记脉冲流数组中的第L个位置C_L加一，并判断标记时间起点数组的第L个位置是否为0，若不为0，则不做处理；若为0，则将该位置标记为时间起点O_k，在标记时间起点数组相应位置置1；

步骤5-6：计算相位值，利用公式相位η₀＝(t_n-O_k)/τ求得相位，其中当(v＝η₀+0.4999…)向下取整为1或大于2且|(η₀/v)-1|<ζ₀时(ζ₀为人为设定的正值，如0.03)，令当前t_n作为时间起点。

步骤5-7：进行PRI变换，对PRI变换数组D_k进行更新，D_k＝D_k+exp(2πiη₀)；

步骤5-8：判断到达时间差τ是否进入第L个以ε₂抖动率构建的交叠箱范围内，再进行步骤5-5、步骤5-6以及步骤5-7，之后判断其是否进入第L个以ε₃抖动率构建的交叠箱范围内，以此类推，当遍历L位置的多个交叠箱后，L＝L+1，当所有交叠箱遍历之后，m＝m+1；

步骤5-9：若PRI变换数组中存在D_k超过门限值，在超过门限的D_k部分中寻找连续区域内的峰值，峰值代表的交叠箱中心即为PRI变换法估计PRI结果。

PRI交叠箱宽度由基本箱宽度与箱中心抖动值二者之间的较大值决定，基本箱宽度b＝(PRI_max-PRI_min)/K，箱中心抖动值由箱中心位置k_L与箱抖动率ε_i共同决定，其中第L个箱中心位置k_L＝(L-1/2)·(PRI_max-PRI_min)/K+PRI_min，则以箱抖动率ε_i构建的第L个交叠箱宽度为b_iL＝max([b,2ε_ik_iL])。

门限设置应考虑观察时间原则、消除子谐波原则以及消除噪声原则，针对第L个位置交叠箱的门限设置应选取三者中的最大值，则第L个交叠箱门限设置

其中α,β,γ为人为设定(如α＝0.34，β＝0.17，γ＝3)，与交叠箱抖动率个数相关，N为当前环境雷达脉冲总数。

步骤6：查询各模块分选结果，完成雷达信号信息输出，如雷达信号类型，PRI等，实现分选。

步骤6-1：将各模块输出信息汇总，如针对捷变频雷达信号分组的SDIF结果、常规雷达信号SDIF结果、参差雷达信号参差分析结果以及PRI变换法分选结果等；

步骤6-2：将汇总的雷达信号结果输出，如雷达信号类型，PRI等。

使用本方法在PC(i5 M460+6GB RAM，Win7+VS2010)上对复杂雷达环境进行测试，加入十部雷达信号与500个噪声点，分选5671个脉冲点仅耗时30ms，所有信号均成功分选，其中大抖动信号参数估计误差小于0.5％，如图5。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。

综上，本发明公开了一种基于SDIF与PRI变换法结合的雷达信号分选方法，该方法将SDIF与PRI变换法结合，对预分选、序列搜索算法、参差检验流程以及PRI变换法进行改进，提升分选算法性能，提高分选速度。

本发明公开了一种基于SDIF与PRI变换法结合的雷达信号分选方法，能够在复杂信号环境中快速完成信号分选，分选效果良好。该方法特点包括：预分选与主分选结合，将SDIF与PRI变换法结合作为主分选，SDIF部分对复杂雷达电磁环境中常规信号、参差信号以及捷变频信号进行快速而有效的分选，PRI变换法部分分选剩余的抖动信号；提出针对捷变频信号的载频预分选设计，解决预分选时PRI预设分选范围限制问题；优化SDIF中的序列搜索算法；提高参差检验流程速度；在PRI变换法中设计多级PRI箱结构，使改进的PRI变换法在保持分选大抖动信号能力的同时，提高PRI估计精度。可证明该方法分选算法性能良好，可快速完成分选，易于工程实现。

本发明的目的在于提供一种基于SDIF与PRI变换法结合的雷达信号分选方法，提出针对捷变频信号的载频预分选设计，解决预分选时结果受PRI预设分选范围限制问题，优化SDIF中的序列搜索算法，使搜索效果更快，更明显，提高参差检验流程速度，在修正PRI变换法中设计多级PRI箱结构，使改进的修正PRI法在保持分选大抖动信号能力的同时，提高PRI估计精度。

本发明的目的是这样实现的：

步骤1：设置一长一短两段观察时间(如1000ms和10ms)，利用雷达信号的脉冲描述字如到达方向DOA、载频CF、脉宽PW等参数，建立各参数直方图，对雷达信号进行预分组，实现捷变频信号与其余类型雷达信号的区分、以及各雷达信号的预分选；

步骤2：根据步骤1预分选结果，利用SDIF方法对预分选后的雷达分组建立到达时间级差直方图(初始级差C为1)，当直方图峰值超过门限值时，则判断该到达时间差值为可能的雷达信号PRI；

步骤3：按照步骤2中估计得到的PRI进行序列搜索，依据雷达信号模型时域特点快速搜索提取雷达信号，并且实现一次提取一部完整的参差雷达信号；

步骤4：若搜索到的PRI中存在相同值，则初步认为存在参差信号，将参差信号脉冲串依时域顺序还原，建立一级差直方图快速分析参差信号；

步骤5：判断经过提取的雷达信号电磁空间内是否仍有复杂类型雷达信号残留，如大范围抖动雷达信号，若存在，则采用多级PRI箱结构(如5％，10％，15％)对分组内剩余脉冲序列进行离散形式PRI变换；

步骤6：查询各模块分选结果，完成雷达信号信息输出，如雷达信号类型，PRI等，实现分选。

本发明还包括这样一些结构特征：

1.步骤1具体包括：

步骤1-1：去除噪声点。噪声点的载频与脉宽等雷达信号参数均是随机分布的，按容差分别建立载频统计直方图与脉宽统计直方图，根据两直方图纵轴高度，设置门限分别清除孤立噪声点与干扰点，如载频统计直方图容差为10MHz，脉宽统计直方图容差为0.5μs，将高度低于10的脉冲点作为噪声点，在针对短观察时间预分选内将高度低于150的点作为干扰点，在针对长观察时间预分选内将高度高于2000的点作为干扰点，从脉冲描述字中删除；

步骤1-2：利用脉宽、到达角度等参数初步预分选。根据去噪后的脉宽、到达角度统计直方图，结合容差，依据脉宽、到达角度等参数将相似信号脉冲点分组，初步稀释信号。容差的设置应尽可能囊括信号脉宽、到达角度的抖动范围，以防止初步预分选造成错分，如设置脉宽容差为3μs，到达角度容差为0.5°；

步骤1-3：区分捷变频信号以及对其他类型信号预分选。在脉宽、到达角度预分选结果中利用载频统计直方图，设置门限将捷变频信号区分开来，将低于门限的信号判断为捷变频信号，将高于门限的信号依据载频容差划分为多个不同载频的其他雷达信号分组，如设置载频容差为100MHz，依据观察时间设置门限为：

T_threshold＝βT/τ β∈(0,1)

其中T为观察时间，τ为当前分选PRI范围中值，预先设置长短观察时间内的两种PRI范围，如短观察时间内PRI范围为0～200μs，长观察时间内PRI范围为200～2000μs，β为人为设定参数。

2.步骤2中当超出门限的峰值个数不唯一时，即同时有多个PRI超过门限，先按对应PRI的一半进行脉冲尝试搜索，若均搜索失败，则直接进行下一级差值直方图，再检测是否有PRI值超过门限。

3.步骤2中门限采用下述公式计算确定：

T_threshold＝α(E-C)e^-λτ

其中λ为单位时间内雷达信号的平均出现概率，E为当前预分组总体雷达脉冲个数，C为正在进行差值级数，(E-C)为脉冲间隔个数，α与脉冲丢失率相关，人为设定，τ为总时间间隔，即脉冲采集时间。

4.步骤3具体包括：

步骤3-1：依据步骤2估计PRI对应的级差直方图，在直方图中搜索符合估计PRI的脉冲个数；

步骤3-2：确定脉冲描述字内所估计PRI的起点脉冲位置；

步骤3-3：进行序列搜索，在搜索过程中，将本次搜索到的脉冲位置存储，由于下次预测脉冲位置必在本次搜索到的脉冲位置之后，在下次搜索时，直接从所存储位置开始搜索；

步骤3-4：在本次搜索结束之后，将本次搜索提取到的脉冲个数与级差直方图中寻找到的符合PRI个数比较，依据二者倍数关系初步判断是否为参差信号，若为参差信号，则在起点脉冲的下一个脉冲再次进行搜索。

5.步骤4具体包括：

步骤4-1：依据步骤2得到的PRI估计值，在一定容差内，根据相同的PRI值个数，初步判断参差数，如容差为1μs；

步骤4-2：利用序列搜索提取到的脉冲流，根据到达时间先后，还原出原始的参差雷达信号；

步骤4-3：对原始参差雷达信号依次进行一级差，并进行直方图统计；

步骤4-4：依据参差数，依次选取频数较高的几个到达时间差值作为子周期。

6.步骤5具体包括：

步骤5-1：判断当前雷达信号分组电磁空间经过SDIF分选提取后，空间内脉冲数是否超过门限(如设置为20)，若超过则进行PRI变换法分选，否则跳过步骤5；

步骤5-2：分别创建K个由多个不同抖动率构建的交叠PRI箱，如多个抖动率依次为ε₁、ε₂、ε₃(ε₁<ε₂<ε₃)，如K取300，分别创建300个5％、10％以及15％的PRI交叠箱，初始化PRI变换数组、标记脉冲流数组、时间起点数组以及标记时间起点数组，数组长度均为K；

步骤5-3：令n＝2，m＝1，τ＝t_n-t_m，判断到达时间差τ是否在预设的PRI范围内，若不在，则m＝m+1再进行判断，直至m＝n-1，若以上到达时间差均不在预设PRI范围内，令n＝n+1继续进行PRI范围判断；

步骤5-4：若到达时间差τ在预设的PRI范围内，则从抖动率为ε₁的第一个交叠箱开始，遍历各箱；

步骤5-5：首先判断到达时间差τ是否在第L个以ε₁抖动率构建的交叠箱范围内，若在其中，则将标记脉冲流数组中的第L个位置C_L加一，并判断标记时间起点数组的第L个位置是否为0，若不为0，则不做处理；若为0，则将该位置标记为时间起点O_k，在标记时间起点数组相应位置置1；

步骤5-6：计算相位值，利用公式相位η₀＝(t_n-O_k)/τ求得相位，其中当(v＝η₀+0.4999…)向下取整为1或大于2且|(η₀/v)-1|<ζ₀时(ζ₀为人为设定的正值，如0.03)，令当前t_n作为时间起点。

步骤5-7：进行PRI变换，对PRI变换数组D_k进行更新，D_k＝D_k+exp(2πiη₀)；

步骤5-8：判断到达时间差τ是否进入第L个以ε₂抖动率构建的交叠箱范围内，再进行步骤5-5、步骤5-6以及步骤5-7，之后判断其是否进入第L个以ε₃抖动率构建的交叠箱范围内，以此类推，当遍历L位置的多个交叠箱后，L＝L+1，当所有交叠箱遍历之后，m＝m+1；

步骤5-9：若PRI变换数组中存在D_k超过门限值，在超过门限的D_k部分中寻找连续区域内的峰值，峰值代表的交叠箱中心即为PRI变换法估计PRI结果。

7.步骤5中PRI交叠箱宽度由基本箱宽度与箱中心抖动值二者之间的较大值决定，基本箱宽度b＝(PRI_max-PRI_min)/K，箱中心抖动值由箱中心位置k_L与箱抖动率ε_i共同决定，其中第L个箱中心位置k_L＝(L-1/2)·(PRI_max-PRI_min)/K+PRI_min，则以箱抖动率ε_i构建的第L个交叠箱宽度为b_iL＝max([b,2ε_ik_iL])。

8.根据权利要求7中所述的方法，其特征在于：步骤5中门限设置应考虑观察时间原则、消除子谐波原则以及消除噪声原则，针对第L个位置交叠箱的门限设置应选取三者中的最大值，则第L个交叠箱门限设置

其中α,β,γ为人为设定(如α＝0.34，β＝0.17，γ＝3)，与交叠箱抖动率个数相关，N为当前环境雷达脉冲总数。

9.步骤6具体包括：

步骤6-1：将各模块输出信息汇总，如针对捷变频雷达信号分组的SDIF结果、常规雷达信号SDIF结果、参差雷达信号参差分析结果以及PRI变换法分选结果等；

步骤6-2：将汇总的雷达信号结果输出，如雷达信号类型，PRI等。