CN110231603B - 一种基于gmti的快速解算目标速度的方法 - Google Patents
一种基于gmti的快速解算目标速度的方法 Download PDFInfo
- Publication number
- CN110231603B CN110231603B CN201910570010.1A CN201910570010A CN110231603B CN 110231603 B CN110231603 B CN 110231603B CN 201910570010 A CN201910570010 A CN 201910570010A CN 110231603 B CN110231603 B CN 110231603B
- Authority
- CN
- China
- Prior art keywords
- speed
- prf
- target
- gmti
- target speed
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请属于机载火控雷达领域,特别涉及一种基于GMTI的快速解算目标速度的方法,包括如下步骤:检测出目标的模糊频率分别为F1、F2以及主杂波频率分别为FZ1、FZ2;分别计算目标相对主杂波的多普勒频移F1‑FZ1+M1f1和F2‑FZ2+M2f2;分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的速度Vt 1和Vt 2;判断速度Vt 1和Vt 2的差值满足解模糊门限阈值Vth;如果满足,则速度Vt 1和Vt 2为目标速度;不满足,则重新确定步骤二中的模糊重数,并重复步骤三到步骤四。本申请的基于GMTI的快速解算目标速度的方法,在解算模糊重数时,根据地面动目标的速度分布特点,以主杂波所在频率为起始进行小范围遍历,能够快速解算出地面动目标的速度信息,软硬件资源需求少。
Description
技术领域
本申请属于机载火控雷达领域,特别涉及一种基于GMTI的快速解算目标速度的方法。
背景技术
地面运动目标显示(Ground Moving Target Indicator,GMTI)是现代机载雷达空地工作模式的一种。地面动目标的周围经常存在着各种背景,例如山、建筑物、树木、鸟群等,这些背景所产生的回波称为杂波。杂波和运动目标在雷达显示器上同时显示时,由于杂波的功率太强而难以观测到目标。如果目标处在杂波背景内,弱的目标淹没在强杂波中,发现目标十分困难;即使目标不在杂波背景内,要在成片的杂波中很快地分辨出运动目标也是十分不容易的。此外,载机本身的运动会使杂波谱大幅度平移和加宽,而且载机快速机动也会导致杂波谱特性快速变化。
一种抑制地面杂波有效的方法是动目标显示(Moving Target Indicator,MTI)。MTI是指利用杂波抑制滤波器来抑制各种杂波,提高雷达信号的信杂比,以利于运动目标检测的技术。其实质是利用动目标回波的多普勒频移来区分动目标和固定目标。常用的MTI滤波器,比如两脉冲对消器和三脉冲对消器等都可以归为正向传输的延迟线滤波器,由n节延时线、n+1个权重系数形成器和求和器构成。当杂波频谱是分布在零频率和脉冲重复频率所有谐波上的谱线时(机载雷达需要考虑地面相对载机的多普勒频率),MTI滤波器可以对主杂波进行抑制。
地面运动目标显示(GMTI)采用具有足够低的脉冲重复频率波形。考虑到低脉冲重复频率波形具有测速模糊的缺点,在进行地面动目标速度求解的过程中,工程师尝试了很多不同的方法,其中主流的有两类:
第一类是通过距离微分解算目标经向速度。这种解算方法得到的速度精度取决于多次距离测量值的精度,误差比较大。
第二类是通过延拓和相关处理的策略解模糊重数、进而得到目标经向速度的方法。大范围地遍历解算模糊重数会导致过多时间和资源的浪费。
因此,急需找到一种基于GMTI的快速解算目标速度的方法。
发明内容
为了解决上述技术问题至少之一,本申请提供了一种基于GMTI的快速解算目标速度的方法。
本申请公开了一种基于GMTI的快速解算目标速度的方法,包括如下步骤:
步骤一、分别对低脉冲重复频率波形PRF1和PRF2回波进行MTI滤波、脉冲压缩后进行快速傅里叶变换,检测出目标的模糊频率分别为F1、F2以及主杂波频率分别为FZ1、FZ2,其中,PRF1对应的多普勒频率为f1,PRF2对应的多普勒频率为f2;
步骤二、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的多普勒频移F1-FZ1+M1f1和F2-FZ2+M2f2,其中,M1和M2分别为低脉冲重复频率波形PRF1和PRF2的模糊重数;
步骤三、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的速度Vt 1和Vt 2;
步骤四、判断步骤三中速度Vt 1和Vt 2的差值满足解模糊门限阈值Vth;如果满足,则速度Vt 1和Vt 2为目标速度;如果不满足,则重新确定步骤二中的模糊重数,并重复步骤三到步骤四。
根据本申请的至少一个实施方式,在所述步骤一中,f1和f2设置为相差在预定数值内的互质数。
根据本申请的至少一个实施方式,在所述步骤三中,分别通过如下公式(3)计算速度Vt 1:
根据本申请的至少一个实施方式,在所述步骤三中,分别通过如下公式(4)计算速度Vt 2:
根据本申请的至少一个实施方式,在所述步骤四中,是通过如下关系式(7)来判断速度Vt 1和Vt 2的差值是否满足解模糊门限阈值Vth:
根据本申请的至少一个实施方式,所述预定数值为10Hz。
本申请至少存在以下有益技术效果:
本申请的基于GMTI的快速解算目标速度的方法,在解算模糊重数时,根据地面动目标的速度分布特点,以主杂波所在频率为起始进行小范围遍历,其方法简单、易于工程实现,能够快速解算出地面动目标的速度信息,软硬件资源需求少。
附图说明
图1是本申请基于GMTI的快速解算目标速度的方法中GMTI原理框图;
图2是当杂波和目标同时存在时的雷达回波的典型功率谱密度示意图;
图3是MTI滤波器的频率响应示意图;
图4是MTI滤波器的输出示意图;
图5是三重PRF解速度模糊(以零频率为起始遍历)示意图;
图6是两重PRF解GMTI速度模糊(以主杂波中心频率为起始遍历)示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
GMTI流程如图1所示。首先对多通道的雷达原始回波数据进行通道校准和运动补偿等前期处理。接着,进行MTI杂波对消。其中MTI滤波器滤除杂波的过程如图2-图4所示,在脉冲雷达中MTI滤波器就是利用杂波与运动目标的多普勒频率的差异,使得滤波器的频率响应在杂波谱的位置形成“凹口”,以抑制杂波,而让动目标回波通过后的损失尽量小。地面雷达的杂波谱通常集中于直流(多普勒频率fd=0)和雷达重复频率fr的整数倍处,但是在机载雷达上还需要考虑地面相对载机的多普勒频率fp。最后,进行地面动目标检测,解算目标角度、距离和速度等信息。
在地面动目标速度信息解算的过程中,考虑到低重复频率脉冲波形具有测距不模糊、测速模糊的特点,工程师主要运用了两类方法。
第一类是通过距离微分解算目标经向速度。这种解算方法得到的速度精度取决于多次距离测量值的精度,误差比较大。
第二类是通过延拓和相关处理的策略解模糊重数、进而得到目标经向速度的方法。解多普勒速度模糊的方法同解算距离模糊类似,即延拓和相关处理。如图5所示,速度延拓就是对检测到的模糊径向速度增加一组带符号的整数倍的PRF速度(第一盲速),即
M∈[-J … 0 ... K] (2);
式中,代表第一盲速(PRF对应的相对速度);Fc是经过质心处理的目标回波多普勒滤波器序号;NFFT是多普勒滤波器组中的滤波器个数,即FFT的点数;[-J…0...K]代表覆盖感兴趣目标最大正负多普勒速度范围的多普勒模糊数。低脉冲重复频率波形回波的速度模糊重数求解需要进行大范围的遍历,时间和资源浪费严重。
鉴于以上两种常用的解地面动目标速度的方法都存在缺陷,本发明提出了一种基于GMTI的快速解算目标速度的方法,包括如下步骤:
步骤一、分别对低脉冲重复频率波形PRF1和PRF2回波进行MTI滤波、脉冲压缩后进行快速傅里叶变换,检测出CAFAR目标的模糊频率分别为F1、F2以及主杂波频率分别为FZ1、FZ2,其中,PRF1对应的多普勒频率为f1,PRF2对应的多普勒频率为f2;其中,f1和f2设置为相差在预定数值内的互质数,预定数值优选为10Hz。
步骤二、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的多普勒频移F1-FZ1+M1f1和F2-FZ2+M2f2,其中,M1和M2分别为低脉冲重复频率波形PRF1和PRF2的模糊重数。
步骤三、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的速度Vt 1和Vt 2;
具体地,分别通过如下公式(3)、(4)计算速度Vt 1和Vt 2:
步骤四、判断步骤三中速度Vt 1和Vt 2的差值满足解模糊门限阈值Vth;如果满足,则速度Vt 1和Vt 2为目标速度;如果不满足,则重新确定步骤二中的模糊重数,并重复步骤三到步骤四。
具体地,在步骤四中,是通过如下关系式(7)来判断速度Vt 1和Vt 2的差值是否满足解模糊门限阈值Vth:
进一步,如图6所示,地面动目标相对主杂波的速度,为目标绝对速度在载机和目标之间连线的投影,一般速度比较小,即有:
max(Vt 1,Vt 2)<Vmax (5);
其中Vmax为地面目标最大的绝对速度。
根据多普勒效应:
一般公路上车辆速度在40m/s(144km/h),即Vmax=40m/s以内,在求解下式时M1和M2在{-3,-2,-1,0,1,2,3}中的整数遍历,可以找到地面动目标速度的准确解:
其中Vth为解算速度误差阈值。
最终,遍历至M1=0、M2=0时,求解满足|Vt 1-Vt 2|≤Vth。
综上所述,本申请的基于GMTI的快速解算目标速度的方法,在解算模糊重数时,根据地面动目标的速度分布特点,以主杂波所在频率为起始进行小范围遍历,其方法简单、易于工程实现,能够快速解算出地面动目标的速度信息,软硬件资源需求少。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种基于GMTI的快速解算目标速度的方法,其特征在于,包括如下步骤:
步骤一、分别对低脉冲重复频率波形PRF1和PRF2回波进行MTI滤波、脉冲压缩后进行快速傅里叶变换,检测出目标的模糊频率分别为F1、F2以及主杂波频率分别为FZ1、FZ2,其中,PRF1对应的多普勒频率为f1,PRF2对应的多普勒频率为f2;
步骤二、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的多普勒频移F1-FZ1+M1f1和F2-FZ2+M2f2,其中,M1和M2分别为低脉冲重复频率波形PRF1和PRF2的模糊重数;
步骤三、分别计算低脉冲重复频率波形PRF1和PRF2条件下目标相对主杂波的速度Vt 1和Vt 2;
步骤四、判断步骤三中速度Vt 1和Vt 2的差值满足解模糊门限阈值Vth;如果满足,则速度Vt 1和Vt 2为目标速度;如果不满足,则重新确定步骤二中的模糊重数,并重复步骤三到步骤四。
2.根据权利要求1所述的基于GMTI的快速解算目标速度的方法,其特征在于,在所述步骤一中,f1和f2设置为相差在预定数值内的互质数。
6.根据权利要求2所述的基于GMTI的快速解算目标速度的方法,其特征在于,所述预定数值为10Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910570010.1A CN110231603B (zh) | 2019-06-27 | 2019-06-27 | 一种基于gmti的快速解算目标速度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910570010.1A CN110231603B (zh) | 2019-06-27 | 2019-06-27 | 一种基于gmti的快速解算目标速度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110231603A CN110231603A (zh) | 2019-09-13 |
CN110231603B true CN110231603B (zh) | 2023-01-06 |
Family
ID=67857500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910570010.1A Active CN110231603B (zh) | 2019-06-27 | 2019-06-27 | 一种基于gmti的快速解算目标速度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110231603B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111308465B (zh) * | 2020-03-26 | 2022-07-12 | 中国气象局气象探测中心 | 一种基于多普勒速度方位滤波识别空中生物的方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4275396A (en) * | 1979-10-12 | 1981-06-23 | Jacomini Omar J | Helicopter rotating blade detection system |
US4529985A (en) * | 1982-08-10 | 1985-07-16 | The United States Of America As Represented By The Secretary Of Transportation | MTI Velocity discriminator |
CN101907702A (zh) * | 2010-06-24 | 2010-12-08 | 西安电子科技大学 | 用于mimo雷达的两维多脉冲对消器 |
CN101975939B (zh) * | 2010-09-28 | 2013-07-17 | 北京航空航天大学 | 一种基于压缩感知理论的多普勒解模糊处理方法 |
CN104239683B (zh) * | 2014-07-23 | 2017-04-19 | 西安空间无线电技术研究所 | 一种基于改变信号调频率的解距离—速度模糊方法 |
CN106569212B (zh) * | 2016-11-09 | 2019-06-18 | 西安空间无线电技术研究所 | 一种多通道sar-gmti距离模糊杂波抑制方法 |
CN107271997B (zh) * | 2017-06-05 | 2020-05-01 | 西北工业大学 | 机载多通道cssar地面运动目标运动参数估计方法 |
CN109814073B (zh) * | 2019-01-21 | 2022-11-11 | 西安电子科技大学 | 一种mtd雷达解模糊测速方法 |
-
2019
- 2019-06-27 CN CN201910570010.1A patent/CN110231603B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110231603A (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111913158B (zh) | 复杂杂波背景下检测低慢小目标的雷达信号处理方法 | |
US4536764A (en) | Method of counting multiple targets in the post detection processing of a radar | |
CN108508439B (zh) | 双机载sar对目标协同成像立体定位的方法 | |
WO2022068097A1 (zh) | 一种改进调频连续波雷达目标检测的降噪方法 | |
CN108089166B (zh) | 一种基于毫米波测云雷达的自适应频域检测方法 | |
CN110988874B (zh) | 一种复杂运动目标的isar成像方法 | |
CN111624574A (zh) | 弱目标探测的目标检测方法、系统、存储介质和装置 | |
CN109061640B (zh) | 一种用于顺轨干涉sar海流反演的方位模糊抑制方法 | |
CN109116326B (zh) | 一种基于中位数估计的自适应雷达海杂波抑制方法 | |
CN110231603B (zh) | 一种基于gmti的快速解算目标速度的方法 | |
CN110146890B (zh) | 一种时频域单通道sar慢速目标检测方法 | |
CN108919250B (zh) | 一种基于多谱线精确插值的低小慢动目标处理方法 | |
CN110133646B (zh) | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 | |
Kovregin et al. | A Unified Method for Observation of an Air Object with a Complex Spectrum in Radar with Quasi-Continuous Radiation | |
CN113221062B (zh) | 一种小型无人机载BiSAR系统的高频运动误差补偿方法 | |
CN108196238B (zh) | 高斯背景下基于自适应匹配滤波的杂波图检测方法 | |
CN111796288B (zh) | 一种基于杂波频谱补偿技术的三坐标雷达动目标处理方法 | |
CN108020835B (zh) | 一种星载合成孔径雷达sar地面运动目标指示gmti的强杂波抑制方法及装置 | |
CN111123269B (zh) | 用于无人机避障雷达的地面杂波抑制方法、模块及装置 | |
US20240069152A1 (en) | Techniques for mitigating interference in radar signals | |
CN109143235A (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
KR101844151B1 (ko) | 고도 클러터를 고려한 이동 플랫폼 레이더의 이동 타겟 탐지 장치 및 방법 | |
CN116449326A (zh) | 宽带多目标平动参数估计与补偿方法 | |
CN115267715A (zh) | 一种强运动杂波环境下弱小目标检测方法 | |
CN114488145A (zh) | 一种基于极化分类的图像域dpca虚警抑制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |