一种在雷达脉冲严重丢失情况下计算真实重频值的方法与流程

本发明涉及雷达信号分选
技术领域：
，尤其是一种在雷达脉冲严重丢失情况下计算真实重频值的方法。
背景技术：
：雷达信号分选是电子侦收设备信号处理的重要环节之一，只有对信号进行正确分选后，才能对雷达信号参数进行有效分析和提取。目前许多分类方法用于雷达信号分选，其中直方图方法(包括累积差值直方图和序列差值直方图)属于雷达信号分选的常用方法之一。在复杂电磁环境中电子侦收设备难以接收雷达发射的所有脉冲，不可避免出现丢失脉冲情况，这会造成直方图方法出现重频的分频现象，进而影响后续分选结果的准确性。为了消除重频分频不利影响，一种减少重频分频的最大公约数的现有技术被提出，该技术鲁棒性较好，适合于工程应用。但是，消除重频分频不利影响的现有技术只适用于分频数互质的情况，实际情况中存在分频数不互质的情况，例如真实重频码值为500，可能会出现4分频、6分频和7分频，分别产生2000、3000、3500重频分频码值，现有方法能从4分频和7分频中得到真实重频500，然而从4分频和6分频中会得到错误的重频值。技术实现要素：本发明所要解决的技术问题是：针对上述消除重频分频不利影响的现有技术只适用于分频数互质的情况，而不适用于分频数不互质的情况的问题，提供一种在雷达脉冲严重丢失情况下计算真实重频值的方法。本发明采用的技术方案如下：一种在雷达脉冲严重丢失情况下计算真实重频值的方法，包括如下步骤：步骤1，输入重频分频值tpri1和tpri2，占空比的最小值dutymin和最大值dutymax，以及脉宽pw；步骤2，根据重频分频值tpri1和tpri2，占空比的最小值dutymin和最大值dutymax，以及脉宽pw；分别计算预设分频数k1和k2的最小值和最大值；计算公式如下：kimin＝tprii×dutymin÷pw,i＝1,2kimax＝tprii×dutymax÷pw,i＝1,2步骤3，循环计算round(tprii/ki),ki∈[kimin,kimax]，结果存入缓存数组pi[]，其中i＝1，2；步骤4，比对p1[]和p2[]，找出满足round(tpri1/k1)≈round(tpri2/k2)的第一个值作为真实重频值。其中，当不存在数据测量误差时，则设定当|round(tpri1/k1)-round(tpri2/k2)|＝0时，round(tpri1/k1)≈round(tpri2/k2)。其中，当存在数据测量误差时，则设定当|round(tpri1/k1)-round(tpri2/k2)|≤1时，round(tpri1/k1)≈round(tpri2/k2)。综上所述，由于采用了上述技术方案，本发明的有益效果是：本发明的在雷达脉冲严重丢失情况下计算真实重频值的方法能够在分频数不互质的情况下计算真实重频值，并且鲁棒性较好，适用于工程应用。附图说明为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。图1为本发明在雷达脉冲严重丢失情况下计算真实重频值的方法的流程框图。具体实施方式为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明，即所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。实施例1设真实重复间隔tpri＝500，由于脉冲丢失的原因，得到虚假的重复间隔测量值tpri1＝2000和tpri2＝3000，其中tpri1的分频数是n1＝4，tpri2对应的分频数是n2＝6，两个分频数不是互质的。已知该雷达占空比的最小值dutymin＝6，最大值dutymax＝15，则占空比范围[dutymin，dutymax]＝[6，15](单位％)，脉宽pw＝30。注：tpri、tpri1、tpri2、pw在工程中被量化的整数值。由于本实施例是不存在数据测量误差的情况，因此，在本实施例中的现有技术和本发明的在雷达脉冲严重丢失情况下计算真实重频值的方法中均设定：当|round(tpri1/k1)-round(tpri2/k2)|＝0时，表示round(tpri1/k1)≈round(tpri2/k2)。根据现有技术，预设分频数最大值为16，得到如表1所示的计算结果。表1，现有技术的计算结果：分频数n12345678910111213141516round(tpri1/k1)20001000667500400333286250222200182167154143133125round(tpri2/k2)300015001000750600500429375333300273250231214200188根据表1，得到计算结果为1000、500、333、250、200时都满足round(tpri1/k1)≈round(tpri2/k2)，选择满足条件的第一个数作为计算结果，即，选择1000作为现有技术的计算结果。然而真实重频的正确值是500。因此当前测试表明当分频数不是互质时，现有技术会产生不正确结果。由此，本发明提供一种在雷达脉冲严重丢失情况下计算真实重频值的方法，包括：(1)输入重频分频值tpri1＝2000和tpri2＝3000，占空比的最小值dutymin＝6，最大值dutymax＝15(单位％)，以及脉宽pw＝30；(2)根据重频分频值tpri1＝2000和tpri2＝3000，占空比的最小值dutymin＝6，最大值dutymax＝15(单位％)，以及脉宽pw＝30分别计算预设分频数k1和k2的最小值和最大值；计算公式如下：kimin＝tprii×dutymin÷pw,i＝1,2kimax＝tprii×dutymax÷pw,i＝1,2计算得到tpri1对应的预设分频数k1的最小值k1min＝4，预设分频数k1的最大值k1max＝10，tpri2对应的预设分频数k2的最小值k2min＝6，预设分频数k2的最大值k2max＝15。(3)循环计算round(tprii/ki),ki∈[kimin,kimax]，结果存入缓存数组pi[]，其中i＝1，2；得到如表2所示的计算结果：表2，本发明的计算结果：分频数n123456789101112131415round(tpri1/k1)------500400333286250222200----------round(tpri2/k2)----------500429375333300273250231214200(4)比对p1[]和p2[]，找出满足round(tpri1/k1)≈round(tpri2/k2)的第一个值作为真实重频值。根据表2，得到计算结果为500、333、250、200时都满足round(tpri1/k1)≈round(tpri2/k2)，选择满足条件的第一个数作为计算结果，最后将500作为计算结果。对照真实重频值的正确值500，得出当前测试数据表明本发明可以得到正确结果。实施例2为了测试算法的鲁棒性和进一步测试算法的准确性，本实施例考虑数据测量误差的影响。由于数据测量误差，因此重复间隔测量值tpri1和tpri2会在一定范围内抖动。设抖动误差为0.2％，tpri1的抖动范围是[1996，2004]，tpri2的抖动范围是[2994，3006]，其余参数与实施例1一致。由于本实施例是存在数据测量误差的情况，因此，在本实施例中的现有技术和本发明的在雷达脉冲严重丢失情况下计算真实重频值的方法中均设定：当|round(tpri1/k1)-round(tpri2/k2)|≤1时，表示round(tpri1/k1)≈round(tpri2/k2)。本发明以及现有技术得到的计算结果分别如表3和表4所示。表3，存在误差情况下本发明的计算结果：表4，存在误差情况下现有技术的计算结果：为了表示方便，表3的计算结果用表示，表4的计算结果用表示，其中m＝1996，1997，…，2004，n＝2994，2995，…，3006，tpri＝500是真实重频值。用下式去衡量算法计算的准确性。计算得到accu1＝2.6％，accu2＝9404.2％。accul值越小表明算法的准确性越高和鲁棒性越强。因为accu1远小于accu2，所以该测试表明本文算法的准确性和鲁棒性都优于现有技术。需要说明的是，在本发明的在雷达脉冲严重丢失情况下计算真实重频值的方法中的计算结果判断条件为当|round(tpri1/k1)-round(tpri2/k2)|≤δ时，表示round(tpri1/k1)≈round(tpri2/k2)，其中，δ在不存在数据测量误差时为0，在存在数据测量误差时，可以根据实际需求进行设定，本发明的实施例2将其优选为1。以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。当前第1页12