CN105445731B - 一种基于高精度测速的雷达全自动跟踪方法 - Google Patents
一种基于高精度测速的雷达全自动跟踪方法 Download PDFInfo
- Publication number
- CN105445731B CN105445731B CN201510796261.3A CN201510796261A CN105445731B CN 105445731 B CN105445731 B CN 105445731B CN 201510796261 A CN201510796261 A CN 201510796261A CN 105445731 B CN105445731 B CN 105445731B
- Authority
- CN
- China
- Prior art keywords
- mark
- msub
- flight path
- mrow
- radar
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 230000000977 initiatory effect Effects 0.000 claims abstract description 5
- 230000001419 dependent effect Effects 0.000 claims description 4
- 230000004304 visual acuity Effects 0.000 claims description 4
- 241001269238 Data Species 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 abstract description 2
- 230000003993 interaction Effects 0.000 description 1
- 238000009333 weeding Methods 0.000 description 1
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
- G01S13/00—Systems 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/66—Radar-tracking systems; Analogous 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
本发明公开了一种基于高精度测速的雷达全自动跟踪方法,利用大于1024阶MTD滤波器积累目标,确保速度分辨力和测速精度;利用CFAR技术完成恒虚警检测,过门限的目标输出,称为点迹;计算点迹的多普勒速度;根据点迹的多普勒速度、航迹径向速度、点迹位置和航迹位置完成点迹、航迹的高关联度相关;剩余点迹根据点迹多普勒速度和点迹位置完成点迹之间的高关联度相关;完成航迹滤波;完成航迹起始;完成杂波(虚警)剔除。本发明使大量的虚警、杂波无法与真实目标相关,因此解决了密集杂波背景下无法进行全自动跟踪的行业难题。该方法适用于可以进行长时间相参积累的雷达,具有算法复杂度低,容易实现的优点。
Description
技术领域
本发明属于雷达数据处理技术领域,具体涉及一种解决对空情报雷达无法实现全自动跟踪问题的基于高精度测速的雷达全自动跟踪方法。
背景技术
雷达系统中,雷达回波经滤波后进行恒虚警检测,对于情报雷达通常要求检测概率为0.5,虚警概率为10-6,但是在复杂电磁环境下往往会有很多的虚假目标,给操作员的工作带来压力。
雷达数据处理主要包括点迹处理和航迹两个阶段,在算法上主要是相关算法和滤波跟踪算法。雷达数据处理的性能是两个算法互相作用的结果,但主要问题出现在相关算法上。在密集杂波背景下,当目标出现机动或者目标暂消以后,在目标预测位置的周围往往会出现虚假点迹,从而造成误相关,这就是当前我国地面情报雷达无法做到全自动跟踪的根本原因。
发明内容
为了解决上述技术问题,本发明利用高精度测速的方法,提出了一种在满足雷达虚警概率和检测概率要求的前提下实现雷达的自动跟踪的方法。
本发明所采用的技术方案是:一种基于高精度测速的雷达全自动跟踪方法,其特征在于,包括以下步骤:
步骤1:利用N阶MTD滤波器积累雷达IQ视频数据,确保速度分辨力和测速精度;其中雷达的工作波长为λ,雷达的重复频率为fr,多普勒频率为fd,频率分辨力为根据公式v=fd×λ/2,得雷达的多普勒速度分辨力为
步骤2:利用CFAR技术完成恒虚警检测,过门限的目标输出,称为点迹;
步骤3:计算点迹的多普勒速度;
步骤4:根据点迹的多普勒速度、航迹径向速度、点迹位置和航迹位置完成点迹、航迹的高关联度相关;
步骤5:剩余点迹根据点迹多普勒速度和点迹位置完成点迹之间的高关联度相关;
步骤6:根据卡尔曼滤波器对相关的航迹进行航迹滤波和预测;
步骤7:根据卡尔曼滤波器对相关的点迹进行航迹起始;
步骤8:对于未相关的点迹进行剔除。
作为优选,步骤2是对MTD后的距离-多普勒矩阵求模后进行CFAR计算,完成恒虚警检测;其中目标处于通道号为n的多普勒通道中,对应的多普勒频率为距离为Rd。
作为优选,步骤3中根据公式v=fd×λ/2求出点迹的多普勒速度vd。
作为优选,步骤4中航迹与点迹的相关规则为:
其中:a为目标的加速度,vd为点迹径向速度,vt为航迹径向速度,T为点迹和航迹更新的时间差,Rt为航迹滤波距离,Rd为点迹距离,dR1为根据点迹和航迹的速度推算出航迹的运动距离,dR2为根据点迹和航迹的距离计算出得航迹运动距离,δ为相关门限,δR为雷达测距精度,vf为雷达速度分辨力,Pf为相关门限因子;
当|dR2-dR1|≤δ时,点迹和航迹相关,否则点迹和航迹不相关。
本发明突破了传统的雷达数据处理的思维定势,创新性的引入了点迹的多普勒速度作为相关的重要因素。该方法使大量的虚警、杂波无法与真实目标相关,因此解决了密集杂波背景下无法进行全自动跟踪的行业难题。该方法适用于可以进行长时间相参积累的雷达,具有算法复杂度低,容易实现的优点。
附图说明
附图1:为本发明实施例的基本流程图;
附图2:为本发明实施例在密集杂波背景下的全自动跟踪实测效果图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
本发明使用的设备包括:雷达信号处理机、数据处理机,还包括显示器、鼠标、键盘等人机交互设备。信号处理机完成雷达的MTD、CFAR和点迹测速,数据处理机完成点迹-航迹相关、点迹-点迹相关、航迹滤波、航迹起始和杂波剔除工作。
设雷达重复周期为1000Hz,雷达工作频率为200MHz(λ为1.5米),此时速度分辨力为,vf=0.37m/s,测速精度为0.3m/s,测距精度为20m。基于高精度测速的全自动跟踪算法如下:
为说明问题,假设空中有2个目标被雷达掌握,目标1的径向速度为200m/s,已经形成稳定航迹;目标2的径向速度为100m/s,尚未形成航迹;同时存在若干虚警。
请见图1,本发明提供的一种基于高精度测速的雷达全自动跟踪方法,包括以下步骤:
步骤1:雷达IQ视频数据经2048阶MTD后,目标得到积累,并在特定的滤波器内得到响应。根据公式fd=2v/λ,目标1的多普勒频率fd应为266.7Hz,在滤波器的546通道响应最大;目标2的多普勒频率fd应为133.3Hz,在滤波器的273通道响应最大。
步骤2:对MTD后的距离-多普勒矩阵求模后进行CFAR计算,将过门限的视频数据输出,其中目标1对应的多普勒通道为546,距离为Rd1;目标2对应的多普勒通道为273,距离为Rd2。
步骤3:根据点迹所处的多普勒通道,根据公式v=fd×λ/2求出点迹的多普勒速度,目标1、目标2对应的多普勒速度分别为266.6m/s和133.3m/s。设目标1所在的位置Rd1=15000米,目标2所在的位置Rd2=50000米。
步骤4:设航迹的速度vt=250m/s,所在的位置为Rt=15250;航迹与点迹的相关规则为:
计算a=-16.6m/s2,dR1=258.3m,dR2=250m,δ=0.37×1+20=20.37,|dR2-dR1|=8.3<δ。在此例中目标1的点迹与目标1的航迹相关;目标1的航迹与目标2的点迹以及其它杂波不会相关。
步骤5:按照第4步的公式进行点迹之间的相关,可以得到目标2本次点迹和上一个周期的点迹相关,目标2本次点迹和其它杂波点不相关;其它杂波点之间因很难满足相关的要求,因此不会相关,即便本周期因极小概率相关了,下一个周期再次相关的概率将会更为微小。
步骤6:根据卡尔曼滤波器对目标1进行航迹滤波和预测。
步骤7:根据卡尔曼滤波器对目标2进行航迹起始。
步骤8:对于未相关的点迹进行剔除。
请见图2,是本发明在密集杂波背景下的全自动跟踪实测效果图,在此时刻存在182个点迹,实际空中目标只有18个,因为对点迹的测速精度达到了5m/s,因此在使用本专利的算法进行全自动跟踪时,没有出现航迹与点迹误相关的情况;没有出现航迹误起始的情况;也没有出现真实目标没有起始航迹的情况。
本发明基于对点迹进行高精度测速的方法,在点迹和航迹相关时考虑距离和速度的二维相关因素,能在航迹暂消和机动时以极大概率避免虚警的误相关,实现雷达全自动跟踪。从而提高雷达的自动化和智能化水平,有利于实现雷达的少人、无人值守。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (4)
1.一种基于高精度测速的雷达全自动跟踪方法,其特征在于,包括以下步骤:
步骤1:利用N阶MTD滤波器积累雷达IQ视频数据,确保速度分辨力和测速精度;其中雷达的工作波长为λ,雷达的重复频率为fr,多普勒频率为fd,频率分辨力为根据公式v=fd×λ/2,得雷达的多普勒速度分辨力为
步骤2:利用CFAR技术完成恒虚警检测,过门限的目标输出,称为点迹;
步骤3:计算点迹的多普勒速度;
步骤4:根据点迹的多普勒速度、航迹径向速度、点迹位置和航迹位置完成点迹、航迹的高关联度相关;
步骤5:剩余点迹根据点迹多普勒速度和点迹位置完成点迹之间的高关联度相关;
步骤6:根据卡尔曼滤波器对相关的航迹进行航迹滤波和预测;
步骤7:根据卡尔曼滤波器对相关的点迹进行航迹起始;
步骤8:对于未相关的点迹进行剔除。
2.根据权利要求1所述基于高精度测速的雷达全自动跟踪方法,其特征在于:步骤2是对MTD后的距离-多普勒矩阵求模后进行CFAR计算,完成恒虚警检测;其中目标处于通道号为n的多普勒通道中,对应的多普勒频率为点迹距离为Rd。
3.根据权利要求1所述基于高精度测速的雷达全自动跟踪方法,其特征在于:步骤3中根据公式v=fd×λ/2求出点迹的多普勒速度vd。
4.根据权利要求1所述基于高精度测速的雷达全自动跟踪方法,其特征在于:步骤4中航迹与点迹的相关规则为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>d</mi>
</msub>
</mrow>
<mi>T</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>dR</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>d</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>T</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&CenterDot;</mo>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>dR</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>f</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>&delta;</mi>
<mi>R</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>|</mo>
<msub>
<mi>dR</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>dR</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>&le;</mo>
<mi>&delta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中:a为目标的加速度,vd为点迹的多普勒速度,vt为航迹径向速度,T为点迹和航迹更新的时间差,Rt为航迹滤波距离,Rd为点迹距离,dR1为根据点迹和航迹的速度推算出航迹的运动距离,dR2为根据点迹和航迹的距离计算出的航迹运动距离,δ为相关门限,δR为雷达测距精度,vf为雷达的多普勒速度分辨力,Pf为相关门限因子;
当|dR2-dR1|≤δ时,点迹和航迹相关,否则点迹和航迹不相关。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510796261.3A CN105445731B (zh) | 2015-11-18 | 2015-11-18 | 一种基于高精度测速的雷达全自动跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510796261.3A CN105445731B (zh) | 2015-11-18 | 2015-11-18 | 一种基于高精度测速的雷达全自动跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105445731A CN105445731A (zh) | 2016-03-30 |
CN105445731B true CN105445731B (zh) | 2017-11-28 |
Family
ID=55556149
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510796261.3A Expired - Fee Related CN105445731B (zh) | 2015-11-18 | 2015-11-18 | 一种基于高精度测速的雷达全自动跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105445731B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108020838B (zh) * | 2016-11-02 | 2021-08-20 | 惠州市德赛西威汽车电子股份有限公司 | 一种毫米波雷达信号在自适应巡航中的处理方法 |
US20190107615A1 (en) * | 2017-10-05 | 2019-04-11 | GM Global Technology Operations LLC | Method of tracking an object |
CN110031834A (zh) * | 2018-01-12 | 2019-07-19 | 西安艾索信息技术有限公司 | 一种改进的多目标雷达航迹处理方法 |
CN108646235B (zh) * | 2018-05-03 | 2022-03-15 | 西安电子工程研究所 | 基于点迹位置聚合识别空间散布固定杂波的方法 |
CN112639524B (zh) * | 2020-04-30 | 2022-05-17 | 华为技术有限公司 | 一种目标检测方法及装置 |
CN111751815B (zh) * | 2020-06-10 | 2024-04-09 | 安徽瞭望科技有限公司 | 一种基于距离微分的虚假航迹剔除方法 |
CN114035165A (zh) * | 2021-11-29 | 2022-02-11 | 南京国睿防务系统有限公司 | 一种炮位雷达的杂波抑制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000376A (zh) * | 2007-01-08 | 2007-07-18 | 清华大学 | 双基地合成孔径雷达的双门限恒虚警运动目标检测方法 |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN103472440A (zh) * | 2013-08-12 | 2013-12-25 | 武汉滨湖电子有限责任公司 | 基于点迹质量判决和航迹质量判决的全自动数据处理方法 |
-
2015
- 2015-11-18 CN CN201510796261.3A patent/CN105445731B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000376A (zh) * | 2007-01-08 | 2007-07-18 | 清华大学 | 双基地合成孔径雷达的双门限恒虚警运动目标检测方法 |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN103472440A (zh) * | 2013-08-12 | 2013-12-25 | 武汉滨湖电子有限责任公司 | 基于点迹质量判决和航迹质量判决的全自动数据处理方法 |
Non-Patent Citations (3)
Title |
---|
CORRELATION OF TARGET REPORTS FROM MTD PROCESSOR AND DATA ANALYSIS;Yang Wen-ling et al.;《Radar 97 (Conf. Publ. No. 449)》;19971016;535-539 * |
WIND TURBINE CLUTTER DETECTION IN SCANNING ATC RADAR;Weikun He et al.;《Integrated Communication, Navigation, and Surveillance Conference (ICNS), 2015》;20150423;U4-1-U4-8 * |
舰载雷达数据处理的杂波抑制方法;徐晓群 等;《现代雷达》;20040831;第26卷(第8期);36-38 * |
Also Published As
Publication number | Publication date |
---|---|
CN105445731A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105445731B (zh) | 一种基于高精度测速的雷达全自动跟踪方法 | |
Xu et al. | Research on computer vision-based for UAV autonomous landing on a ship | |
Moyer et al. | A multi-dimensional Hough transform-based track-before-detect technique for detecting weak targets in strong clutter backgrounds | |
CN104215951B (zh) | 一种在海杂波背景下慢速小目标的检测系统及其检测方法 | |
CN105510896B (zh) | 一种集中式多雷达数据处理的加权最近邻域数据关联方法 | |
CN102621542B (zh) | 基于多模粒子滤波和数据关联的机动微弱目标检测前跟踪方法 | |
CN103336947B (zh) | 基于显著性和结构性的红外运动小目标识别方法 | |
CN104133211B (zh) | 一种多普勒频率变换雷达目标分类识别方法 | |
CN105301584B (zh) | 同时解距离模糊的ipphdf机动多目标跟踪方法 | |
CN104714225B (zh) | 一种基于广义似然比的动态规划检测前跟踪方法 | |
CN103345587B (zh) | Ads‑b监视数据与雷达航迹的直觉模糊关联方法、装置 | |
CN105093215A (zh) | 基于多普勒信息的雷达对低空慢速小目标的跟踪方法 | |
CN109100696B (zh) | 基于点状杂波图的慢速运动目标消除方法 | |
CN109086648A (zh) | 一种融合目标检测和特征匹配的目标跟踪方法 | |
CN101840579A (zh) | 一种利用视频分割与粒子滤波实现多目标跟踪方法 | |
CN116794647A (zh) | 基于轨迹跟踪的导线和塔架分类 | |
CN104459686B (zh) | 一种基于Hough变换与蚁群相似度的目标检测与跟踪方法 | |
CN103985127A (zh) | 一种密集恒星背景的弱小目标检测方法及装置 | |
Li et al. | The application of Yolov4 and a new pedestrian clustering algorithm to implement social distance monitoring during the COVID-19 pandemic | |
CN108414992B (zh) | 一种基于相位信息杂波图的目标检测方法 | |
CN117331071A (zh) | 一种基于毫米波雷达与视觉多模态融合的目标检测方法 | |
Wang et al. | A survey of hand gesture recognition based on FMCW radar | |
Niu et al. | UAV detection based on improved YOLOv4 object detection model | |
CN114114242A (zh) | 一种基于点选航迹及航迹外推的雷达目标追赶跟踪方法 | |
CN105353352B (zh) | 改进搜索策略的mm‑pphdf机动多目标跟踪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: No. 718 Luoyu Road, Hongshan District, Wuhan City, Hubei Province 430074 China Shipbuilding Heavy Industry 709 Institute Patentee after: China Shipbuilding Lingjiu Electric (Wuhan) Co.,Ltd. Address before: No. 718 Luoyu Road, Hongshan District, Wuhan City, Hubei Province 430074 China Shipbuilding Heavy Industry 709 Institute Patentee before: CSIC (WUHAN) LINCOM ELECTRONICS Co.,Ltd. |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171128 |