发明内容
针对现有技术的不足,本发明提供了一种脉冲多普勒雷达解速度模糊的方法,能够在运动目标的回波信息并不理想的情况下,解多普勒速度模糊,从而获得运动目标的真实速度。
为解决上述技术问题,本发明采用的一个技术方案是:一种脉冲多普勒雷达解速度模糊的方法,包括:
步骤A:根据雷达参数与目标视在速度,以及第一门限对速度估算区间进行区间划分,筛选出可能的速度矩阵V_maybe。
步骤1、获得目标视在速度V_martix:
整理雷达各个脉冲重复周期PRI获得的过检测门限的视在速度,形成速度矩阵,记为V_martix。
多重频多载频联合的有益效果是:充分利用雷达的固有设计参数,通过多重频多载频对同一目标的多次检测,获取该运动目标的速度集中性;且所使用的多个重频、多个载频的数目、数值均可以按实际雷达自身的载频参数进行设定和调整,具有泛用性和灵活性。
步骤2、计算各类解多普勒速度模糊所需参数:
计算模糊速度V_blind、模糊倍数blur_mul和速度分段数segment_num。
第一门限即为划分速度区间的间隔值V_seg,该值不能太小,一般V_seg取值需要大于雷达的测速精度;且不能太大,取值过大会使过多的真伪速度落在同一个速度区间,在后续筛选中存在遗漏真速度的风险,从而失去了第一门限的意义。
第一门限的有益效果是:在不遗漏真速度的情况下,对真伪速度进行初步的筛选和划分;通过速度区间的划分,使真伪速度尽可能的落在不同的速度区间,便于后续第二门限的进一步筛选,减少了第二门限判别时的计算量。
步骤3、计算可能速度矩阵V_maybe:
计算出可能速度矩阵V_maybe,V_maybe为3维矩阵。
步骤B:根据第二门限,对可能的速度矩阵V_maybe进行筛选,获得待判定真伪的速度矩阵V_check。
步骤4、对可能速度矩阵V_maybe进行重排和映射:
将V_maybe重排为2维矩阵V_reshape。
将V_reshape中满足速度估算区间[-V_max,V_max]要求的速度值映射到速度估算区间;将映射后的速度矩阵记为V_map。
并用V_segment储存V_map中每个元素位于速度分段区间的位置。
步骤5、多PRI速度映射表融合:
将V_segment中所有行求和,得到V_segment_all;找出V_segment_all中≥3的元素位置储存为pos_maybe矩阵。
以下步骤为循环递归操作,为了方便描述,特此定义参数n:初始值为1,取值为1到N的自然数,每执行一次步骤7后自加1。
步骤6、筛选满足速度估算区间的速度矩阵:
查找V_segment的第pos_maybe[n]列中元素值>0的行数,记为pri_pos。
将V_reshape的第pri_pos列中满足速度估算区间的可能速度值映射到速度估算区间;将映射后的速度矩阵记为V_yinshe。
找出V_yinshe中与pos_maybe[n]数值最接近的列数位置,将该位置的速度值存入V_lose_check矩阵;并将该速度在原速度矩阵V_martix中的位置,存入V_lose_check_pos矩阵。
步骤7、筛选V_lose_check中满足第二门限的速度矩阵:
将V_lose_check中满足速度估算区间的速度组成一个V_match矩阵;统计V_match中元素的个数记为V_zhi。
将V_match矩阵所有元素中的最大值和最小值作差,第二门限deta取值为2;如果该差值大于第二门限deta,跳转执行步骤6;如果该差值小于第二门限deta,跳转执行步骤8。
第二门限的有益效果是:第二门限通过分析第一门限筛选后的速度的集中性,对可能的真实速度进行筛选,进一步剔除假速度。
步骤8、计算待判定真伪的速度矩阵V_check:
将可能的速度V_match平均值和其置信度V_zhi存入V_check矩阵的第k行(k初始化为1,每执行一次步骤8后加1)。
直到遍历pos_maybe中所有列后(即n=N)结束本段递归循环,否则跳转执行步骤6。
将递归循环完成后获得的V_check的行数记为M行,列数为7列。
步骤C:对待判定真伪的速度矩阵V_check,进行速度互异性筛选,剔除伪速度,获得按置信度从高到低排列的目标真实速度矩阵V_real。
以下步骤为循环递归操作,为了方便描述,特此定义参数m:初始值为1,取值为1到M的自然数,每次执行步骤15后自加1,每次执行步骤16后初始化为1。
步骤9、取出置信度最高的速度行:
取出V_check矩阵中置信度最高的行,记为V_max_zhi,并存入V_real矩阵第k行(k初始化为1,每执行一次步骤9后加1)。
步骤10、计算速度互异性判别矩阵Judge矩阵:
将V_max_zhi和V_check矩阵,从矩阵的第3列起作差,得到Judge矩阵。
创建速度真伪标志位矩阵。
步骤11、查询Judge矩阵和V_flag速度真伪标志位矩阵:
当Judge矩阵第m行元素中存在0值,并且V_flag矩阵第m行元素为0时,则继续执行步骤12,否则跳转执行步骤15。
步骤12、进行盲速PRI判别:
如果V_max_zhi的第3列到第7列中存在0值,存在0,执行步骤13.2;否则执行步骤13.1。
根据是否存在盲速PRI的判别结果,执行下列对应处理分支的操作步骤。
步骤13、速度互异性筛选:
步骤13.1、不存在盲速,不进行速度互异性筛选:
跳转执行步骤15。
步骤13.2、存在盲速,进行速度互异性筛选:
将Judge矩阵第m行中元素为0的列的位置的集合A与V_max_zhi第3列到第7列中元素为0的列的位置集合B作对比:
如果A中有不属于B的元素,继续执行步骤14;否则跳转执行步骤15。
速度互异性筛选的有益效果是:运用一个真实目标速度能且只能被一个PRI检测出来的原理,进一步剔除假速度。
步骤14、更改速度矩阵的置信度和速度真伪标志位矩阵:
更改V_check和V_flag中第m行的值。
步骤15、重复上述步骤11至14,直至遍历所有M行个速度行。
步骤16、遍历完M行后,查询速度真伪标志位V_flag中是否存在0。若存在,则重复步骤9至16;若不存在,则结束本部分的循环递归操作。
步骤17、结束上述循环递归操作后,获得的V_real矩阵即为解多普勒速度模糊后的目标真实速度矩阵。
本发明的有益效果是:本发明对可能的速度矩阵进行步骤A、步骤B、步骤C三次速度筛选,具有较高的置信度和抗信噪比波动性;与传统的提高重频解决多普勒速度模糊方法相比,避免了提高重频有可能导致距离模糊的风险;与单纯的使用多普勒频移公式进行理论计算的方法相比,极大的消除了目标回波信号信噪比不理想时对解速度模糊结果准确性的影响;通过对速度聚集性的分析与判别,实现了脉冲多普勒雷达的解速度模糊。本发明中,重频、载频的数目和数值都可以按实际雷达自身的载频参数进行设定和调整,使该方法能够适用于各体制的脉冲多普勒雷达系统。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
本发明实施例提供了一种脉冲多普勒雷达解速度模糊的方法,可应用于雷达信号处理算法领域。
请参阅图1和图2,一种脉冲多普勒雷达解速度模糊的方法包括如下步骤:
本实施例中使用的雷达脉冲重复周期PRI个数为5个,记为PRI1、PRI2、PRI3、PRI4、PRI5,分别为461us、481us、500us、528us、549us;对应的载频记为Fc1、Fc2、Fc3、Fc4、Fc5,分别为5550MHz、5550MHz、5550MHz、5550MHz、5550MHz;并假设共有3个需要解多普勒速度模糊的运动目标,其真实速度记为V_real1、V_real2和V_real3,分别为550m/s、-100m/s、212m/s。
解多普勒速度的速度估算区间为[-V_max,V_max]=[-1000,1000],速度区间间隔值为V_seg=2;即解出来的目标真实速度范围区间为[-1000,1000],精度为2。
下面将本方法实施步骤分为3个主要步骤,详细描述本方法的实施流程:
步骤A:根据雷达参数与目标视在速度,以及第一门限对速度估算区间进行区间划分,筛选出可能的速度矩阵V_maybe。
步骤1、获得目标视在速度V_martix:
整理雷达各个脉冲重复周期PRI获得的过检测门限的视在速度,形成一个3行*5列的速度矩阵,记为V_martix;其中每个PRI获得的速度均排列成一列,当某个PRI获得的目标视在速度个数不足3个时,在矩阵中将该位置用-1填充表示,便于后续的处理运算。
本实施例中使用的V_martix矩阵如下:
步骤2、计算各类解多普勒速度模糊所需参数:
步骤2.1、计算模糊速度V_blind:
根据多普勒频率公式V_blind=frλ/2,其中脉冲重复频率fr=1/PRI,波长λ=C/Fc,计算出模糊速度V_blind。
步骤2.2、计算模糊倍数blur_mul:
根据公式blur_mul=ceil(V_max/V_blind),计算出速度估算区间内的模糊倍数blur_mul;其中,ceil表示向上取整。
步骤2.3、计算速度分段数segment_num:
划分速度估算区间,根据公式segment_num=2*ceil(V_max/V_seg),计算出速度区间分段数segment_num;其中,ceil表示向上取整。
第一门限即为划分速度区间的间隔值V_seg,该值不能太小,一般V_seg取值需要大于雷达的测速精度;且不能太大,取值过大会使过多的真伪速度落在同一个速度区间,在后续筛选中存在遗漏真速度的风险,从而失去了第一门限的意义。
步骤3、计算可能速度矩阵V_maybe:
根据公式V_maybe=V_martix+V_blind×blur_mul计算出可能速度矩阵,V_martix中为-1填充的速度不参与计算,为便于后续的速度筛选,V_maybe中该速度对应的位置暂时先用一个任意大值填充,如1000。
V_maybe是一个三维矩阵,第一维是同一PRI下速度;第二维是每个目标真实速度对应的解模糊后的可能速度;第三维是PRI。在本实施例中V_maybe第一维长度为3,记为
第二维长度由模糊倍数blur_mul决定,记为
第三维长度为5,记为
步骤B:根据第二门限,对可能的速度矩阵V_maybe进行筛选,获得待判定真伪的速度矩阵V_check。
步骤4、对可能速度矩阵V_maybe进行重排和映射:
4.1V_maybe的重排:
为便于后续计算,将V_maybe重排为行数为
列数为
的2维矩阵,记为V_reshape;列代表速度,行代表PRI。
4.2V_reshape的映射:
找出V_reshape中满足速度估算区间[-V_max,V_max]的速度,根据速度区间分段数segment_num和速度区间间隔值V_seg,将满足估算区间要求的速度值映射到速度估算区间;将映射后的速度矩阵记为V_map。
并建立一个行数PRI个数*列数速度区间分段数segment_num的速度区间映射表,记为V_segment;将V_map中每个元素位于速度分段区间的位置在该表中通过置‘1’标记出来,代表该PRI在该速度区间可能检测出真实的运动目标。
步骤5、多PRI速度映射表融合:
将V_segment中所有行求和,得到一个行数1*列数segment_num的矩阵,记为V_segment_all;该矩阵代表所有PRI在各个速度区间测得的可能目标速度的个数。
运用3/5准则,找出V_segment_all中≥3的元素位置,位置值储存在一个1行*N列的矩阵pos_maybe中;代表该目标速度至少被5个PRI中的任意3个所检测出来,才认为该速度有可能是真实的目标速度,才对该速度进行后续的处理。
以下步骤为循环递归操作,为了方便描述,特此定义参数n:初始值为1,取值为1到N的自然数,每次执行步骤7后自加1。
步骤6、筛选满足速度估算区间的速度矩阵:
以pos_maybe中第n列储存的速度位置作为列数序号,查找V_segment的该列中元素值>0的行数,即得到检测出该速度的PRI序号值,记为pri_pos。
返回到V_reshape矩阵中,找出V_reshape中pri_pos所代表的PRI序号值该列中满足速度估算区间[-V_max,V_max]的可能速度值,并将该速度值映射到速度估算区间;将映射后的速度矩阵记为V_yinshe。
以pos_maybe中第n列储存的速度位置与V_yinshe中的每个元素进行对比,找出与之数值最接近的列数位置,将该位置的速度值存入V_lose_check矩阵;并找到该速度位置在原速度矩阵V_martix中的位置,即该速度对应原速度矩阵第几个PRI的第几个速度,并将其结果存入V_lose_check_pos矩阵,行数表示该速度是由第几个PRI检测出;V_lose_check中的速度值和V_lose_check_pos中的位置值一一对应。
步骤7、筛选V_lose_check中满足第二门限的速度矩阵:
将V_lose_check矩阵所有元素中满足速度估算区间[-V_max,V_max]的速度组成一个V_match矩阵;统计V_match中元素的个数,此个数为该速度为可能真实速度的置信度,记为V_zhi。
同时将V_match矩阵所有元素中的最大值和最小值作差,第二门限deta取值为2;如果该差值大于第二门限deta,则说明该组速度不具有聚集性,其中存在伪速度,并跳转执行步骤6(查找通过步骤5生成的V_segment_all矩阵和通过步骤4.1生成的V_reshape矩阵,获取有效速度值V_lose_check矩阵);如果该差值小于第二门限deta,则说明该组速度具有聚集性,故跳转执行步骤8。
第二门限用于分析经过第一门限筛选后的运动目标可能速度的集中性,对可能的真实速度进行筛选,更好的剔除伪速度。
步骤8、计算待判定真伪的速度矩阵V_check:
创建一个列数为7的V_check矩阵,将V_match中所有元素求平均值的结果存入V_check的第k行第1列中(k初始化为1,每执行一次步骤8后加1);将V_zhi存入V_check的第k行第2列中;将V_lose_check_pos的第一列的5个数分别存入V_check的第k行第3至7列。
直到遍历pos_maybe中所有列后(即n=N)结束本段递归循环,否则跳转执行步骤6。
上述递归循环完成后获得的V_check的行数记为M行,列数为7列;下面将对V_check矩阵进行速度互异性筛选,并剔除其中的伪速度。
步骤C:对待判定真伪的速度矩阵V_check,进行速度互异性筛选,剔除伪速度,获得按置信度从高到低排列的目标真实速度矩阵V_real。
以下步骤为循环递归操作,为了方便描述,特此定义参数m:初始值为1,取值为1到M的自然数,每次执行步骤15后自加1,每次执行步骤16后初始化为1。
步骤9、取出置信度最高的速度行:
取出V_check矩阵中置信度最高的速度行,当存在相同置信度的速度行时,可按照信号检测幅度排序,采纳信号幅度最大组合的置信度;记为V_max_zhi,并将该行存入V_real矩阵第k行(k初始化为1,每执行一次步骤9后加1)。
步骤10、计算速度互异性判别矩阵Judge矩阵:
将V_max_zhi和V_check矩阵,从矩阵的第3列起,每行各个元素分别作差,得到一个M行*5列的Judge矩阵。
同时创建一个初始化为全0的M行*1列的V_flag矩阵,作为速度真伪标志位矩阵。
步骤11、查询Judge矩阵和V_flag速度真伪标志位矩阵:
当Judge矩阵第m行元素中存在0值,并且V_flag矩阵第m行元素为0时,则继续执行步骤12,否则跳转执行步骤15。
步骤12、进行盲速PRI判别:
查询置信度最高的速度行V_max_zhi,如果该行的第3列到第7列的元素中存在0值,则说明5个PRI中存在盲速PRI的情况;反之则不存在盲速PRI的情况。
根据是否存在盲速PRI的判别结果,执行下列对应处理分支的操作步骤。
步骤13、速度互异性筛选:
所谓速度互异性,就是一个真实的目标速度不可能被同一PRI检测出两次,若被同一PRI检测出两次,那么这两个被检测出的目标速度中至少有一个速度是伪速度。
步骤13.1、不存在盲速,不进行速度互异性筛选:
不存在盲速PRI,则不存在一个速度被同一PRI检测出两次的情况,不需要进行速度互异性筛选;此时跳转执行步骤15。
步骤13.2、存在盲速,进行速度互异性筛选:
将Judge矩阵第m行中元素为0的列的位置的集合A与V_max_zhi第3列到第7列中元素为0的列的位置集合B作对比:
如果A中有不属于B的元素,则说明第m行速度不满足速度互异性条件,该行对应的速度为伪速度,需要被剔除,故继续执行步骤14;
如果A中的元素均属于B,则说明第m行速度满足速度互异性条件,该速度将在后续的循环递归操作中进行解析和判决,故跳转执行步骤15。
步骤14、更改速度矩阵的置信度和速度真伪标志位矩阵:
将V_check中第m行的置信度更改为0,并且将V_flag中第m行的数值改为-1。代表该速度行对应的速度是伪速度。
步骤15、重复上述步骤11至14,直至遍历所有M行个速度行。
步骤16、遍历完所有M行个速度行后,查询速度真伪标志位矩阵V_flag,看矩阵中是否存在为0的值。若存在0,则重复步骤9至16;若不存在0,则结束本部分的循环递归操作。
步骤17、结束上述循环递归操作后,获得的V_real矩阵即为解多普勒速度模糊后的目标真实速度矩阵,每行对应一个解出来的目标真实速度;其中第1列为速度值,第2列为置信度,第3列至第7为各个PRI检测到的速度值序号。
本实施例中,经过多重频多载频联合解速度模糊获得的V_real矩阵如下:
第一行,[550,5,1,1,1,1,1];第二行,[212,4,3,3,0,3,3];第三行,[-100,3,2,2,2,0,0]。可以看到,3个运动目标的真实速度550、212、-100均已经被解析出来;三个真实速度置信度分别为5,4,3,该结果说明了本方法的正确性和可实施性。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。