Nothing Special   »   [go: up one dir, main page]

CN110673105B - 一种脉冲多普勒雷达解速度模糊的方法 - Google Patents

一种脉冲多普勒雷达解速度模糊的方法 Download PDF

Info

Publication number
CN110673105B
CN110673105B CN201910907887.5A CN201910907887A CN110673105B CN 110673105 B CN110673105 B CN 110673105B CN 201910907887 A CN201910907887 A CN 201910907887A CN 110673105 B CN110673105 B CN 110673105B
Authority
CN
China
Prior art keywords
speed
matrix
value
check
row
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910907887.5A
Other languages
English (en)
Other versions
CN110673105A (zh
Inventor
桂雨洋
徐敏超
陈钟
厉芸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
WUHAN BINHU ELECTRONIC Co.,Ltd.
Original Assignee
Wuhan Binhu Electronic Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Binhu Electronic Co ltd filed Critical Wuhan Binhu Electronic Co ltd
Priority to CN201910907887.5A priority Critical patent/CN110673105B/zh
Publication of CN110673105A publication Critical patent/CN110673105A/zh
Application granted granted Critical
Publication of CN110673105B publication Critical patent/CN110673105B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及雷达信号处理领域,特别是涉及一种脉冲多普勒雷达解速度模糊的方法。本发明包括三个主要步骤,分别为,步骤A:根据雷达参数与目标视在速度,以及第一门限对速度估算区间进行区间划分,筛选出可能的速度矩阵;步骤B:根据第二门限,对可能的速度矩阵进行筛选,获得待判定真伪的速度矩阵;步骤C:对待判定真伪的速度矩阵,进行速度互异性筛选,剔除伪速度,获得按置信度从高到低排列的目标真实速度矩阵。本发明能够在运动目标的回波信息并不理想且不引入距离模糊的情况下,解多普勒速度模糊,从而获得运动目标的真实速度。

Description

一种脉冲多普勒雷达解速度模糊的方法
技术领域
本发明涉及雷达信号处理领域,特别是涉及一种脉冲多普勒雷达解速度模糊的方法。
背景技术
速度模糊问题是脉冲多普勒雷达的固有问题,当雷达发射脉冲的重复频率小于运动目标的最大多普勒频率时,会产生速度模糊;从而使雷达获取不到运动目标的真实速度,导致无法实现雷达的目标测速、目标跟踪等后续功能。
要消除速度模糊,目前普遍的解决方法有两个:一个是提高雷达发射脉冲的重复频率,但过高的脉冲重频又会导致运动目标的距离模糊;另一个方法是使用多个不同的重频同时对运动目标进行探测,并通过多普勒频移公式计算出目标真实速度,但由于真实的目标回波信息的信杂噪比受到外部环境的影响,本身存在误差和畸变,导致此方法在实际工程应用时并不能通过理论公式获得解析解,无法解出运动目标的真实速度。
故需要有一种适用于脉冲多普勒雷达实际使用场景,可以在运动目标的回波信息信杂噪比并不理想的情况下,解多普勒速度模糊的方法。
发明内容
针对现有技术的不足,本发明提供了一种脉冲多普勒雷达解速度模糊的方法,能够在运动目标的回波信息并不理想的情况下,解多普勒速度模糊,从而获得运动目标的真实速度。
为解决上述技术问题,本发明采用的一个技术方案是:一种脉冲多普勒雷达解速度模糊的方法,包括:
步骤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是一种脉冲多普勒雷达解速度模糊的方法具体实施例中步骤A和步骤B流程示意图;
图2是一种脉冲多普勒雷达解速度模糊的方法具体实施例中步骤C流程示意图;
附图中各部件的标记如下:
1、雷达的基本参数和运动目标视在速度;
2、模糊速度、模糊倍数和速度区间分段数等参数;
3、计算可能的速度矩阵V_maybe;
4、计算满足速度估算区间的速度映射矩阵V_map;
5、计算多PRI融合后的速度映射表V_segment_all;
6、搜索速度矩阵V_maybe中满足速度估算区间的有效速度值;
7、进一步搜索速度矩阵V_maybe中满足第二门限deta的可能速度值;
8、计算待判断真伪的速度矩阵V_check;
9、取出拥有最高置信度的速度行;
10、计算速度互异性判别Judge矩阵。
11、创建真伪标志位矩阵;
12、盲速判别;
13、盲速判别后是否进行速度互异性筛选的分支处理;
14、更新V_flag矩阵和V_check矩阵;
15、递归循环,遍历V_check中所有速度行;
16、判断真伪标志位矩阵V_flag;
17、解多普勒速度模糊后得到的真实速度矩阵V_real。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
本发明实施例提供了一种脉冲多普勒雷达解速度模糊的方法,可应用于雷达信号处理算法领域。
请参阅图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矩阵如下:
Figure GDA0003263279230000081
步骤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,记为
Figure GDA0003263279230000091
第二维长度由模糊倍数blur_mul决定,记为
Figure GDA0003263279230000092
第三维长度为5,记为
Figure GDA0003263279230000093
步骤B:根据第二门限,对可能的速度矩阵V_maybe进行筛选,获得待判定真伪的速度矩阵V_check。
步骤4、对可能速度矩阵V_maybe进行重排和映射:
4.1V_maybe的重排:
为便于后续计算,将V_maybe重排为行数为
Figure GDA0003263279230000094
列数为
Figure GDA0003263279230000095
的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矩阵如下:
Figure GDA0003263279230000131
第一行,[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,该结果说明了本方法的正确性和可实施性。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (2)

1.一种脉冲多普勒雷达解速度模糊的方法,其特征在于:包括以下步骤:
步骤A:根据雷达参数与目标视在速度,以及第一门限对速度估算区间进行区间划分,筛选出可能的速度矩阵V_maybe;V_maybe=V_martix+V_blind×blur_mul,其中V_martix为雷达获得的目标视在速度,V_blind为模糊速度,blur_mul为模糊倍数;具体过程为:
步骤1、获得目标视在速度V_martix:
步骤2、计算解多普勒速度模糊所需参数V_blind、blur_mul:
步骤3、计算可能速度矩阵V_maybe:
步骤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中所有列后结束本段递归循环,否则跳转执行步骤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矩阵;创建一个初始化为全0的M行*1列的V_flag矩阵,作为速度真伪标志位矩阵;
步骤11、查询Judge矩阵和V_flag速度真伪标志位矩阵;
当Judge矩阵第m行元素中存在0值,并且V_flag矩阵第m行元素为0时,则继续执行步骤12,否则跳转执行步骤15;
步骤12、进行盲速PRI判别;
根据是否存在盲速PRI的判别结果,执行下列对应处理分支的操作步骤;
步骤13、速度互异性筛选:
步骤13.1、不存在盲速,不进行速度互异性筛选:此时跳转执行步骤15;
步骤13.2、存在盲速,进行速度互异性筛选:
将Judge矩阵第m行中元素为0的列的位置的集合A与V_max_zhi第3列到第7列中元素为0的列的位置集合B作对比:
如果A中有不属于B的元素,继续执行步骤14;否则跳转执行步骤15;
步骤14、更改速度矩阵的置信度和速度真伪标志位矩阵;
步骤15、重复上述步骤11至14,直至遍历所有M行个速度行;
步骤16、遍历完所有M行个速度行后,查询速度真伪标志位矩阵V_flag中是否存在为0;若存在0,则重复步骤9至16;若不存在0,则结束本部分的循环递归操作;
步骤17、结束上述循环递归操作后,获得的V_real矩阵即为解多普勒速度模糊后的目标真实速度矩阵。
2.根据权利要求1的一种脉冲多普勒雷达解速度模糊的方法,其特征在于:步骤A中根据V_blind=frλ/2,其中脉冲重复频率fr=1/PRI,波长λ=C/Fc计算模糊速度V_blind;
根据公式blur_mul=ceil(V_max/V_blind),计算模糊倍数blur_mul;其中,ceil表示向上取整;
根据公式segment_num=2*ceil(V_max/V_seg),划分速度估算区间,计算出速度区间分段数segment_num;其中,ceil表示向上取整,其中V_seg为划分速度区间的间隔值;V_max为速度估算区间最大值的绝对值。
CN201910907887.5A 2019-09-25 2019-09-25 一种脉冲多普勒雷达解速度模糊的方法 Active CN110673105B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910907887.5A CN110673105B (zh) 2019-09-25 2019-09-25 一种脉冲多普勒雷达解速度模糊的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910907887.5A CN110673105B (zh) 2019-09-25 2019-09-25 一种脉冲多普勒雷达解速度模糊的方法

Publications (2)

Publication Number Publication Date
CN110673105A CN110673105A (zh) 2020-01-10
CN110673105B true CN110673105B (zh) 2021-12-10

Family

ID=69078929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910907887.5A Active CN110673105B (zh) 2019-09-25 2019-09-25 一种脉冲多普勒雷达解速度模糊的方法

Country Status (1)

Country Link
CN (1) CN110673105B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111289953B (zh) * 2020-01-14 2022-04-26 北京理工大学 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法
CN111983579B (zh) * 2020-08-14 2021-05-14 成都亘波雷达科技有限公司 一种参差脉冲重复时间消除雷达速度模糊的方法
CN112305528B (zh) * 2020-10-26 2023-07-11 西安电子工程研究所 一种相控阵雷达重复周期优选方法
CN112363144B (zh) * 2020-11-27 2024-07-23 西安空间无线电技术研究所 环扫雷达距离模糊和方位模糊鉴别方法
CN116559853B (zh) * 2023-07-10 2023-09-12 中国电子科技集团公司第十四研究所 一种低重频解速度模糊方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030210175A1 (en) * 2002-02-22 2003-11-13 Bernhard Bickert Method for reducing false alarm rate in radar images
CN103592645A (zh) * 2013-11-22 2014-02-19 中国电子科技集团公司第五十四研究所 一种伪随机码调相连续波雷达速度模糊解算方法
CN107526067A (zh) * 2017-03-07 2017-12-29 中国气象局武汉暴雨研究所 全自动多普勒天气雷达径向速度退模糊算法
JP2018180302A (ja) * 2017-04-13 2018-11-15 キヤノン株式会社 像ブレ補正装置、交換レンズ、カメラアシステム、および像ブレ補正方法
CN109164422A (zh) * 2018-09-26 2019-01-08 西安电子科技大学 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法
CN109642944A (zh) * 2016-07-09 2019-04-16 德克萨斯仪器股份有限公司 用于包括速度模糊度分辨率的mimo雷达中的速度检测的方法和设备
CN109814073A (zh) * 2019-01-21 2019-05-28 西安电子科技大学 一种mtd雷达解模糊测速方法
CN109932700A (zh) * 2019-03-28 2019-06-25 北京润科通用技术有限公司 一种多普勒速度的解模糊方法及装置
CN109975780A (zh) * 2019-04-17 2019-07-05 西安电子工程研究所 基于脉冲多普勒雷达时域回波的直升机型号识别算法
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7764839B2 (en) * 2003-08-14 2010-07-27 Fujifilm Corporation Edge detecting apparatus and method, and image size enlarging and reducing apparatus and method

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030210175A1 (en) * 2002-02-22 2003-11-13 Bernhard Bickert Method for reducing false alarm rate in radar images
CN103592645A (zh) * 2013-11-22 2014-02-19 中国电子科技集团公司第五十四研究所 一种伪随机码调相连续波雷达速度模糊解算方法
CN109642944A (zh) * 2016-07-09 2019-04-16 德克萨斯仪器股份有限公司 用于包括速度模糊度分辨率的mimo雷达中的速度检测的方法和设备
CN107526067A (zh) * 2017-03-07 2017-12-29 中国气象局武汉暴雨研究所 全自动多普勒天气雷达径向速度退模糊算法
JP2018180302A (ja) * 2017-04-13 2018-11-15 キヤノン株式会社 像ブレ補正装置、交換レンズ、カメラアシステム、および像ブレ補正方法
CN109164422A (zh) * 2018-09-26 2019-01-08 西安电子科技大学 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法
CN109814073A (zh) * 2019-01-21 2019-05-28 西安电子科技大学 一种mtd雷达解模糊测速方法
CN109932700A (zh) * 2019-03-28 2019-06-25 北京润科通用技术有限公司 一种多普勒速度的解模糊方法及装置
CN109975780A (zh) * 2019-04-17 2019-07-05 西安电子工程研究所 基于脉冲多普勒雷达时域回波的直升机型号识别算法
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Study on suppressing multipath active interference by sidelobe cancellation;X. Minchao, L. Qing, L. Wenwu, Z. Sanwei and Y. Jingang;《2016 CIE International Conference on Radar (RADAR), 2016》;20171005;1-5 *
一种步进频相参处理检测低信杂比目标的方法;徐敏超;《火控雷达技术》;20181231;1-8 *
脉冲多普勒雷达的速度模糊求解;洪一;《现代电子》;19990131;18-23 *

Also Published As

Publication number Publication date
CN110673105A (zh) 2020-01-10

Similar Documents

Publication Publication Date Title
CN110673105B (zh) 一种脉冲多普勒雷达解速度模糊的方法
CN105277923A (zh) 一种单通道雷达信号分选方法
CN107229033B (zh) 基于高度维分段搜索的多目标到达时间差定位方法
CN103364770B (zh) 基于矩阵填充的雷达目标检测系统及其检测方法
CN105699952B (zh) 海杂波k分布形状参数的双分位点估计方法
CN110146873B (zh) 分布式非相参雷达的目标位置和速度的估计方法
CN112835009B (zh) 一种基于fmcw雷达rai的数据集拓展方法
CN110598738B (zh) 一种基于空间信息的分布式聚类方法
CN111487592B (zh) 一种基于gpu和cpu平台漂移均值聚类的雷达信号分选方法
CN106468771A (zh) 一种低可观测高杂波条件下的多目标检测与跟踪方法
US5396250A (en) Spectral estimation of radar time-of-arrival periodicities
CN108919249A (zh) 一种基于二维局部插值的雷达目标距离联合估计方法
CN107255814A (zh) 一种基于lfmsk波形的雷达目标检测方法
CN103592620A (zh) 长基线时差定位系统中高重频信号解定位模糊方法
CN111999726B (zh) 一种基于毫米波雷达的人员定位方法
JP2010133960A (ja) アンテナ回転からアンテナ回転にかけてのレーダブリップの多次元時間的運動フィルタリングの方法
CN107817501A (zh) 一种可变扫描频率的点云数据处理方法
CN113608193A (zh) 一种基于UNet的雷达多目标距离和速度估计方法
CN105092022B (zh) 基于抹香鲸叫声的水下目标测距单脉冲筛选方法及装置
CN106772357A (zh) 信噪比未知条件下的ai‑phd滤波器
CN111289953B (zh) 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法
US10802109B1 (en) Techniques for building pulse trains of emission sources
CN106324577A (zh) 一种基于标准差椭圆的高分辨雷达检测点聚集方法
CN113390406A (zh) 基于被动多传感器系统的多目标数据关联与定位方法
CN111624575A (zh) 一种未知雷达目标脉冲序列的快速提取方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right

Effective date of registration: 20191231

Address after: No. 103-2, Hong Kong Road, Jiang'an District, Wuhan City, Hubei Province

Applicant after: Xu Huamin

Address before: No. 51, Liufang Road, Donghu New Technology Development Zone, Wuhan City, Hubei Province

Applicant before: Wuhan Binhu Electronic Co., Ltd.

TA01 Transfer of patent application right
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20200609

Address after: 51 No. 430205 Hubei province Wuhan City Liufang high tech Development Zone East Lake Road

Applicant after: WUHAN BINHU ELECTRONIC Co.,Ltd.

Address before: No. 103-2, Hong Kong Road, Jiang'an District, Wuhan City, Hubei Province

Applicant before: Xu Huamin

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant