申请/专利权人：南京长峰航天电子科技有限公司

申请日：2018-11-30

公开（公告）日：2023-03-24

公开（公告）号：CN109507647B

主分类号：G01S7/36

分类号：G01S7/36

优先权：

专利状态码：有效-授权

法律状态：2023.03.24#授权;2019.04.16#实质审查的生效;2019.03.22#公开

摘要：本发明公开了一种基于PRI直方图的雷达信号分选方法，包括：建立雷达信号脉冲流的PRI直方图；为固定信号和抖动信号设定不同的抖动范围；先进行固定信号提取，再进行抖动信号提取，具体为：将PRI直方图中各PRI交迭箱的箱中值分别与判别门限相比较：若不存在超过判别门限的值，则将PRI直方图中所有PRI交迭箱的箱中值清零；否则，采用超过判别门限的峰值对应的PRI值作为可能的PRI值进行序列搜索；若成功搜索到序列，则从雷达信号脉冲流中剔除已检测到的脉冲，否则，将PRI直方图中所有PRI交迭箱的箱中值清零，进行下一阶次的脉冲提取。本发明能够提高抖动信号分选能力，对复杂环境下的脉冲数据有较好的分选效果。

主权项：1.基于PRI直方图的雷达信号分选方法，其特征在于，所述方法包括：建立雷达信号脉冲流的PRI直方图；分别为待提取的固定信号和抖动信号设定不同的抖动范围；采用下述方法首先进行固定信号提取，然后再进行抖动信号提取：将PRI直方图中各PRI交迭箱的每个箱中值分别与判别门限相比较：若不存在超过判别门限的箱中值，则将PRI直方图中所有PRI交迭箱的箱中值清零；否则，采用超过判别门限的峰值对应的PRI值作为可能的PRI值进行序列搜索；若成功搜索到序列，则从雷达信号脉冲流中剔除已检测到的脉冲，否则，将PRI直方图中所有PRI交迭箱的箱中值清零，进行下一阶次的脉冲提取。

全文数据：基于PRI直方图的雷达信号分选方法技术领域本发明涉及一种基于PRI直方图的雷达信号分选方法，属于电子对抗技术领域。背景技术雷达信号分选指的是从雷达侦察接收机接收到的随机交错的脉冲流中，将同一雷达的脉冲分离出来并选出有用的信号的过程。基于脉冲重复间隔PRI的方法一直是信号分选中研究最多且成果最为丰富的方法之一。PRI搜索法、PRI直方图、PRI变换都是这类方法的典型代表。PRI直方图包括三种典型的方法：TOA差值直方图、积累差直方图CDIF以及时序差直方图SDIF。TOA全差直方图是将每个TOA都与后续的TOA相减，并对每个TOA差值进行累计，得到脉冲序列的全差直方图。该方法主要缺点有:检测判别门限不是最佳判别门限；脉冲丢失严重时虚警较高；只适合定PRI情况，所需计算的差值数较多且谐波问题严重等。CDIF算法是一种基于直方图统计和序列搜索的混合算法，CDIF算法较传统的直方图方法在计算量和抑制谐波方面作了很大的改进，并且由于积累的效果，使得CDIF还具有对干扰脉冲和脉冲丢失不敏感的特点。SDIF是对CDIF提出的一种改进算法，和CDIF相比，两者的主要区别是SDIF针对不同阶的到达时间差直方图不进行累积，相应的检测判别门限也与CDIF不同。SDIF的优点在于相对CDIF减少了计算量，但由于不进行级间积累，使其性能有所下降，该算法最大的不足是不能对重频抖动信号进行处理。在实际系统中，如何确保实时性和复杂环境下分选结果的可靠性，这是非常困难的任务。因此亟须一种实时性好稳定性高的分选方法。发明内容本发明所要解决的技术问题是克服现有技术的缺陷，提供一种基于PRI直方图的雷达信号分选方法，以解决现有技术中导致的上述缺陷或缺陷之一。为达到上述目的，本发明实施例所提供的技术方案为：基于PRI直方图的雷达信号分选方法，包括：建立雷达信号脉冲流的PRI直方图；分别为待提取的固定信号和抖动信号设定不同的抖动范围；采用下述方法首先进行固定信号提取，然后再进行抖动信号提取：将PRI直方图中各PRI交迭箱的每个箱中值分别与判别门限相比较：若不存在超过判别门限的箱中值，则将PRI直方图中所有PRI交迭箱的箱中值清零；否则，采用超过判别门限的峰值对应的PRI值作为可能的PRI值进行序列搜索；若成功搜索到序列，则从雷达信号脉冲流中剔除已检测到的脉冲，否则，将PRI直方图中所有PRI交迭箱的箱中值清零，进行下一阶次的脉冲提取。进一步的，所述固定信号的PRI交迭箱宽度与抖动信号的PRI交迭箱宽度不同。进一步的，所述方法还包括：对所提取的固定信号进行参差分析，以提取可能的参差信号。进一步的，所述信号提取的方法包括如下步骤：步骤1，初始化：令差值阶数C＝1，把待分选PRI区间分为K个PRI交迭箱，确定每个PRI交迭箱的中心位置τl，初始化每个PRI交迭箱的箱中值Dl＝0，l＝1～K，l为PRI交迭箱索引，τl为第l个PRI交迭箱的中心位置；Dl为第l个PRI交迭箱的箱中值；如果此时脉冲流中脉冲总数不大于5个，则结束信号提取；否则，进入下一步；步骤2，令当前采样脉冲序号m＝1，令n＝m+C；步骤3，计算当前差值阶数下第n个采样脉冲到达时间tn与第m个采样脉冲的脉冲到达时间tm的时间差τ，则τ＝tn-tm；如果τ1-ετmin或者τ＞1-ετmax，则转到步骤6；否则，进入下一步；其中，τmin为PRI最大值，τmax为PRI最大值，ε为待提取信号的抖动范围，固定信号和抖动信号提取时抖动范围设置不同；步骤4，计算PRI交迭箱的下限计算PRI交迭箱的上限其中Δτ＝τmax-τminK；步骤5，l＝l1～l2，Dl＝Dl+1；步骤6，m＝m+1，n＝m+C，如果n等于采样脉冲的总个数N，则转到步骤8，否则转到步骤3；步骤7，把直方图中每个PRI交迭箱的箱中值与判别门限比较，如果所有箱中值都没有超过判别门限，则执行C＝C+1，把直方图中所有PRI交迭箱的箱中值清零，转到步骤2，否则，进入下一步；步骤8，对PRI交迭箱中超过判别门限的最大箱中值所对应的附近箱中值执行清零操作，然后进行谐波检测；所述附近箱中值＝最大箱中值×预设比例；步骤9，序列搜索：如果成功搜索到序列，则从脉冲流中剔除已检测信号脉冲，转到步骤1；否则，令C＝C+1，把PRI直方图中所有PRI交迭箱的箱中值清零，转到步骤2。进一步的，所述中心位置τl的计算方法为：第l个交迭箱的宽度为bl＝2ετl。进一步的，所述步骤7判别门限采用下述公式计算确定：Tτ＝xE-Ce-τkN式中：E为脉冲总数；C为差值阶数；k、x均为小于1的正常数；N为交迭箱的划分总数。进一步的，确定序列搜索的PRI值的方法包括：找到过判别门限的最大箱中值；根据预设比例，查找不小于最大箱中值×预设比例的箱中值；计算所查找出的最大箱中值对应交迭箱的中心位置与最小箱中值对应交迭箱的中心位置的平均值，将该平均值作为序列搜索的PRI值。进一步的，所述谐波检测的方法包括：查找PRI交迭箱的最大箱中值对应的PRI值，假如最大箱中值没有过判别门限Tτ，则查找PRI交迭箱中过判别门限的箱中值对应的PRI值；若该过判别门限PRI值是某未过判别门限PRI交迭箱对应的PRI值的谐波，则用未过判别门限PRI交迭箱对应的PRI值进行序列搜索；若过判别门限PRI最小值不是未过判别门限PRI交迭箱对应的PRI值的整数倍，则用过判别门限PRI最小值进行序列搜索；若PRI交迭箱最的箱中值过判别门限，则以过判别门限的第一个箱中值对应的PRI值进行序列搜索。进一步的，所述步骤10的序列搜索方法为：对于固定信号或抖动信号，需搜索连续多个脉冲与PRI值匹配，才认为匹配成功。与现有技术相比，本发明实施例所达到的有益效果是：1实时处理复杂环境下的脉冲信号，分选出正确的信号；2通过改进PRI交迭窗，有效的处理抖动信号，便于工程实现；3通过先提取固定信号，然后在剩余脉冲中提取抖动信号，处理方法简单高效。附图说明图1是本发明实施例提供的固定信号提取流程图。抖动信号提取流程相同，区别在于不需要进行参差分析具体实施方式下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。本发明实施例提供的基于PRI直方图的雷达信号分选方法，首先从雷达信号的脉冲流中提取固定信号，然后从剩余脉冲流中提取抖动信号。分选方法主要包括：PRI直方图建立和序列搜索。PRI直方图建立采用交迭PRI箱的思想，对固定信号和抖动信号取不同的宽度建立不同的PRI交迭箱。序列搜索包括PRI估计和根据PRI进行的序列搜索。如图1所示，是本发明实施例提供的固定信号提取流程图，还包括参差信号提取，具体为，在固定信号提取结束后再进行参差分析，提取出可能的参差信号。抖动信号提取流程与固定信号提取流程基本相同，区别在于不需要进行参差分析。下面对信号提取的方法做进一步详细描述，包括如下步骤：步骤1，初始化：令差值阶数C＝1，把待分选PRI区间分为K个PRI交迭箱，确定每个PRI交迭箱的中心位置τl，初始化每个PRI交迭箱的箱中值Dl＝0，l＝1～K，l为PRI交迭箱索引，τl为第l个PRI交迭箱的中心位置；Dl为第l个PRI交迭箱的箱中值；如果此时脉冲流中脉冲总数不大于5个，则结束信号提取；否则，进入下一步；步骤2，令当前采样脉冲序号m＝1，令n＝m+C；步骤3，计算当前差值阶数下第n个采样脉冲到达时间tn与第m个采样脉冲的脉冲到达时间tm的时间差τ，则τ＝tn-tm；如果τ1-ετmin或者τ＞1-ετmax，则转到步骤6；否则，进入下一步；其中，τmin为PRI最大值，τmax为PRI最大值，ε为待提取信号的抖动范围，固定信号和抖动信号提取时抖动范围设置不同；步骤4，计算PRI交迭箱的下限计算PRI交迭箱的上限其中Δτ＝τmax-τminK；步骤5，l＝l1～l2，Dl＝Dl+1；步骤6，m＝m+1，n＝m+C，如果n等于采样脉冲的总个数N，则转到步骤8，否则转到步骤3；步骤7，把直方图中每个PRI交迭箱的箱中值与判别门限比较，如果所有箱中值都没有超过判别门限，则执行C＝C+1，把直方图中所有PRI交迭箱的箱中值清零，转到步骤2，否则，进入下一步；步骤8，对PRI交迭箱中超过判别门限的最大箱中值所对应的附近箱中值执行清零操作，然后进行谐波检测；所述附近箱中值＝最大箱中值×预设比例；步骤9，序列搜索：如果成功搜索到序列，则从脉冲流中剔除已检测信号脉冲，转到步骤1；否则，令C＝C+1，把PRI直方图中所有PRI交迭箱的箱中值清零，转到步骤2。本实施例中第l个交迭箱的中心位置τl的计算方法为：第l个交迭箱的宽度为bl＝2ετl，对于固定信号，ε可取值为±1％，对于抖动信号，ε可取值为±15％。本实施例中，判别门限计算公式为：Tτ＝xE-Ce-τkN，式中：E为脉冲总数；C为差值阶数；k、x均为小于1的正常数，具体的，k可取0.1，x可取0.2；N为交迭箱的划分总数。本实施例中，确定序列搜索的PRI值的方法包括：找到过判别门限的最大箱中值；根据预设比例，查找不小于最大箱中值×预设比例的箱中值；计算所查找出的最大箱中值对应交迭箱的中心位置与最小箱中值对应交迭箱的中心位置的平均值，将该平均值作为序列搜索的PRI值。预设比例可选取为85％。本实施例中，所述步骤8的谐波检测方法为：查找PRI交迭箱的最大箱中值对应的PRI值，假如最大箱中值没有过判别门限，则查找PRI交迭箱中过判别门限的箱中值对应的PRI值；若该过判别门限PRI值是某未过判别门限PRI交迭箱对应的PRI值的谐波，则用未过判别门限PRI交迭箱对应的PRI值进行序列搜索；若过判别门限PRI最小值不是未过判别门限PRI交迭箱对应的PRI值的整数倍，则用过判别门限PRI最小值进行序列搜索；若PRI交迭箱最的箱中值过判别门限，则以过判别门限的第一个箱中值对应的PRI值进行序列搜索。实施过程最多考虑6倍谐波。本实施例中，所述步骤8的序列搜索方法为：对于固定信号或抖动信号，需搜索连续中间最多允许丢一个脉冲多个脉冲与PRI值匹配，才认为匹配首先产生1个参差、1个固定和2个抖动信号，参差信号重复间隔为PRI1＝1000us210.5us,315.8us,473.6us三参差，固定信号重复间隔为PRI2＝340us，抖动信号重复间隔为PRI3＝400us，PRI4＝600us，抖动范围为±10％；然后进行固定信号提取，步骤1ε取值为±1％，按照步骤1～步骤7的顺序进行PRI直方图建立；PRI直方图建立后，按步骤8和步骤9，进行PRI估计；按步骤10，进行序列搜索，得到输出信号，并且删除相应的脉冲；固定信号提取结束后，进行抖动信号提取，步骤1ε取值为±15％，按照步骤1～步骤7的顺序进行PRI直方图建立；PRI直方图建立后，按步骤8和步骤9，进行PRI估计；按步骤10，进行序列搜索，得到输出信号，并且删除相应的脉冲；经由上述处理，分别得到固定提取信号和抖动提取信号，最后对固定提取的信号进行参差信号提取，满足参差条件的信号为参差信号。上述固定信号提取直方图最大阶数Q取值为8，抖动信号提取直方图最大阶数Q取值为6，剩余脉冲最小个数值固定提取和抖动提取都取6个脉冲。以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。

权利要求：1.基于PRI直方图的雷达信号分选方法，其特征在于，所述方法包括：建立雷达信号脉冲流的PRI直方图；分别为待提取的固定信号和抖动信号设定不同的抖动范围；采用下述方法首先进行固定信号提取，然后再进行抖动信号提取：将PRI直方图中各PRI交迭箱的每个箱中值分别与判别门限相比较：若不存在超过判别门限的箱中值，则将PRI直方图中所有PRI交迭箱的箱中值清零；否则，采用超过判别门限的峰值对应的PRI值作为可能的PRI值进行序列搜索；若成功搜索到序列，则从雷达信号脉冲流中剔除已检测到的脉冲，否则，将PRI直方图中所有PRI交迭箱的箱中值清零，进行下一阶次的脉冲提取。2.根据权利要求1所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述固定信号的PRI交迭箱宽度与抖动信号的PRI交迭箱宽度不同。3.根据权利要求1所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述方法还包括：对所提取的固定信号进行参差分析，以提取可能的参差信号。4.根据权利要求1所述的基于PRI直方图的雷达信号分选方法，其特征是在于，所述信号提取的方法包括如下步骤：步骤1，初始化：令差值阶数C＝1，把待分选PRI区间分为K个PRI交迭箱，确定每个PRI交迭箱的中心位置τl，初始化每个PRI交迭箱的箱中值Dl＝0，l＝1～K，l为PRI交迭箱索引，τl为第l个PRI交迭箱的中心位置；Dl为第l个PRI交迭箱的箱中值；如果此时脉冲流中脉冲总数不大于5个，则结束信号提取；否则，进入下一步；步骤2，令当前采样脉冲序号m＝1，令n＝m+C；步骤3，计算当前差值阶数下第n个采样脉冲到达时间tn与第m个采样脉冲的脉冲到达时间tm的时间差τ，则τ＝tn-tm；如果τ1-ετmax，则转到步骤6；否则，进入下一步；其中，τmin为PRI最大值，τmax为PRI最大值，ε为待提取信号的抖动范围，固定信号和抖动信号提取时抖动范围设置不同；步骤4，计算PRI交迭箱的下限计算PRI交迭箱的上限其中Δτ＝τmax-τminK；步骤5，l＝l1～l2，Dl＝Dl+1；步骤6，m＝m+1，n＝m+C，如果n等于采样脉冲的总个数N，则转到步骤8，否则转到步骤3；步骤7，把直方图中每个PRI交迭箱的箱中值与判别门限比较，如果所有箱中值都没有超过判别门限，则执行C＝C+1，把直方图中所有PRI交迭箱的箱中值清零，转到步骤2，否则，进入下一步；步骤8，对PRI交迭箱中超过判别门限的最大箱中值所对应的附近箱中值执行清零操作，然后进行谐波检测；所述附近箱中值＝最大箱中值×预设比例；步骤9，序列搜索：如果成功搜索到序列，则从脉冲流中剔除已检测信号脉冲，转到步骤1；否则，令C＝C+1，把PRI直方图中所有PRI交迭箱的箱中值清零，转到步骤2。5.根据权利要求4所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述中心位置τl的计算方法为：第l个交迭箱的宽度为bl＝2ετl。6.根据权利要求4所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述步骤7判别门限采用下述公式计算确定：Tτ＝xE-Ce-τkN式中：E为脉冲总数；C为差值阶数；k、x均为小于1的正常数；N为交迭箱的划分总数。7.根据权利要求4所述的基于PRI直方图的雷达信号分选方法，其特征在于，确定序列搜索的PRI值的方法包括：找到过判别门限的最大箱中值；根据预设比例，查找不小于最大箱中值×预设比例的箱中值；计算所查找出的最大箱中值对应交迭箱的中心位置与最小箱中值对应交迭箱的中心位置的平均值，将该平均值作为序列搜索的PRI值。8.根据权利要求4所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述谐波检测的方法包括：查找PRI交迭箱的最大箱中值对应的PRI值，假如最大箱中值没有过判别门限Tτ，则查找PRI交迭箱中过判别门限的箱中值对应的PRI值；若该过判别门限PRI值是某未过判别门限PRI交迭箱对应的PRI值的谐波，则用未过判别门限PRI交迭箱对应的PRI值进行序列搜索；若过判别门限PRI最小值不是未过判别门限PRI交迭箱对应的PRI值的整数倍，则用过判别门限PRI最小值进行序列搜索；若PRI交迭箱最的箱中值过判别门限，则以过判别门限的第一个箱中值对应的PRI值进行序列搜索。9.根据权利要求1所述的基于PRI直方图的雷达信号分选方法，其特征在于，所述步骤10的序列搜索方法为：对于固定信号或抖动信号，需搜索连续多个脉冲与PRI值匹配，才认为匹配成功。

百度查询： 南京长峰航天电子科技有限公司 基于PRI直方图的雷达信号分选方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD