CN106405510B - 一种基于伪滑窗l判决准则的航迹删除方法 - Google Patents
一种基于伪滑窗l判决准则的航迹删除方法 Download PDFInfo
- Publication number
- CN106405510B CN106405510B CN201610801724.5A CN201610801724A CN106405510B CN 106405510 B CN106405510 B CN 106405510B CN 201610801724 A CN201610801724 A CN 201610801724A CN 106405510 B CN106405510 B CN 106405510B
- Authority
- CN
- China
- Prior art keywords
- track
- sliding window
- pseudo
- tracking
- association
- 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
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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- 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/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
-
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
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
本发明提供一种基于伪滑窗的航迹删除方法,考虑了在窄波束边扫描变跟踪雷达跟踪系统中,不同目标被观测到的时间间隔不同的情况,跟踪的同时根据雷达波束扫描规律实时判定目标是否被雷达回波照射,雷达实时跟踪情况将滑窗区分为有效和伪滑窗,对航迹进行管理的时候剔除掉伪滑窗,以忽略雷达照射不到目标所在区域时其对航迹删除决策的影响,若航迹在连续L个有效的滑窗内没关联到量测,则认为该航迹为虚假航迹并删除。由于航迹删除过程中,将划窗划分为有效和伪滑窗,只关注有效滑窗,这种方法称为伪滑窗L准则。本发明在跟踪的同时,根据雷达波束的实时监控区域,对滑窗动态调整,并据此对航迹做出合理的判决。
Description
技术领域
本发明属于雷达信号处理技术领域,特别设计边扫描变跟踪雷达系统多目标跟踪技术领域。
背景技术
边扫描变跟踪雷达是一种多目标跟踪雷达,在跟踪目标的同时分出一部分波束对空域进行搜索。由于一部雷达系统的时间能量资源有限,要实现对更多批次目标的搜索和跟踪,需要尽可能提高雷达系统的时间能量资源利用率。为此,出现了各种系统资源自适应调度算法。其中,针对多目标跟踪问题提出的自适应调整数据率的目标跟踪算法一般采用的策略为:在满足跟踪精度要求情况下,根据目标的运动状态自适应选择合适的数据更新率,尽可能减少对每个目标的累计照射时间,以消耗最少的系统时间能量资源。一般,当目标机动性较大时,采用高采样率进行扫描;目标机动性较小时采用低采样率。尽管交互式多模型滤波器可以用于机动性目标跟踪问题,但当目标的机动性很大时,交互式多模型算法很难反应目标的机动性,因此仍然需要较高的采样率。存在的问题是:当抽样频率很高时,由于跟踪波束每次只能照射到部分区域,需要连续多个扫描波束才能扫完整个监控区域。于是,目标可能连续多帧检测不到。
在杂波干扰背景下,跟踪过程中可能会将杂波误认为一个目标并建立一条航迹,这就需要在跟踪过程中判断哪些是虚假航迹并予以删除,以提高跟踪性能。一个可行的方法就是在航迹更新过程中根据某种准则删除掉一部分航迹,而这个准则就是所谓的航迹删除准则。在传统跟踪算法中,一种常用的航迹删除技术是基于滑窗L判决准则航迹删除算法:当一条航迹连续几帧都没有关联到点迹,而只能靠外推预测位置进行更新时,那么这条航迹很可能就是一条虚假航迹。但由于在实际的边扫描边跟踪系统中,雷达需要多次扫描才能覆盖整个区域,目标在这多个扫描间隔内可能只能被观测到一次。即使是已经被锁定的目标,也可能由于目标数量过大而与很多目标共享一个雷达波束。因此,目标很可能在连续多帧内都检测不到,那么根据传统的航迹删除方法,真实航迹就可能会被误认为是虚假航迹而删除。
发明内容
本发明所要解决的技术问题是,提供一种能够跟踪多个采样间隔随时间变化的目标的航迹确认从而进行航迹删除的方法。
本发明为解决上述技术问题所采用的技术方案是,一种基于伪滑窗L判决准则的航迹删除方法,包括以下步骤:
步骤1、初始化参数,包括数据关联门限γ,航迹删除阈值L;
步骤2、当观测的帧数k达到航迹判决基数时,对第k帧量测zk与所有航迹进行数据关联;
步骤3、根据关联结果更新航迹的状态:若航迹关联到量测,则基于卡尔曼滤波算法更新目标状态;否则,不对目标状态进行更新;
步骤4、根据关联结果更新航迹对应的滑窗:
若第i条航迹关联到量测,则更新k时刻的第i条航迹关联标志Mt(k,i)=1,i=1,2,…,m,m为航迹条数;
若第i条航迹关联未到量测,且目标的预测位置不在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=0;1、0均表示有效滑窗;
若第i条航迹关联未到量测,且目标的预测位置在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=inf;inf表示伪滑窗;
步骤5、查找1到k时刻第i条航迹的航迹关联标志Mt(:,i)提取有效滑窗向量M';
步骤6、对有效滑窗向量M'根据L判决准则进行航迹删除。
本发明考虑了在窄波束边扫描变跟踪雷达跟踪系统中,不同目标被观测到的时间间隔不同的情况,跟踪的同时根据雷达波束扫描规律实时判定目标是否被雷达回波照射,雷达实时跟踪情况将滑窗区分为有效和伪滑窗,对航迹进行管理的时候剔除掉伪滑窗,以忽略雷达照射不到目标所在区域时其对航迹删除决策的影响,若航迹在连续L个有效的滑窗内没关联到量测,则认为该航迹为虚假航迹并删除。由于航迹删除过程中,将划窗划分为有效和伪滑窗,只关注有效滑窗,这种方法称为伪滑窗L准则。本发明在跟踪的同时,根据雷达波束的实时监控区域,对滑窗动态调整,并据此对航迹做出合理的判决。
本发明的有益效果是,有效地解决了在窄波束边扫描变跟踪雷达跟踪系统中传统航迹删除方法可能将真实航迹被误认为虚假航迹的问题,简单易行,可以在边扫描变跟踪雷达体系中及微弱目标跟踪中实现虚假航迹的删除,降低了航迹误删了概率。
附图说明
图1为本发明的流程框图;
图2为雷达检测后得到的点迹和真实航迹图;
图3为本发明与传统基于滑窗L判决准则航迹删除技术的成功跟踪航迹条数对比图;
图4为本发明与传统基于滑窗L判决准则航迹删除技术的跟踪精度对比图。
具体实施方式
为了方便描述本发明的内容,首先对以下术语进行解释:
1:L判决准则
现有用于判断航迹是否为虚假航迹的逻辑判决准则,当目标连续L帧内没有被观测到,则认为该航迹为虚假航迹。
2:伪滑窗
滑窗在本发明中指用于标记n个假设事件的向量,若第k(k=1,2,……,n)个事件为真,则第k个滑窗标记为1,否则,则标记为0。本发明中增加了一种中间状态inf,中间状态的滑窗认为是无效的,进行相关处理时会删除无效滑窗,因此称为伪滑窗。有效滑窗指标记为1或0的滑窗。
3:恒虚警检测
在实际干扰环境下提供可预知的检测和虚警,从而可以准确地设定检测门限的一种现有技术,具有这种特性的检测器称为恒虚警检测器。
4:数据关联
在跟踪目标时,除了附加噪声的不确定性外,观测还具有杂波造成的来源不确定性。对量测数据和估计航迹进行配对的过程称为数据关联,数据关联为一种常用的数据处理方法。
本发明提出了一种基于伪滑窗L判决准则的航迹删除技术,步骤如下:
步骤1、初始化参数,包括数据关联门限γ,航迹删除阈值L,假设已经成功起始航迹,时间帧k=3;
步骤2、第k帧量测与所有航迹进行数据关联,n为检测到的量测个数,如果满足下式,则认为关联成功:
v'(k)S-1(k)v(k)<γ2
其中,’表示转置,v(k)表示量测zk与预测位置的差值,分别为k时刻的状态预测,H表示状态转移矩阵
S(k)表示v(k)的协方差矩阵;这里使用量测与预测位置的马氏距离小于固定门限值γ则认为量测和航迹来自同一目标的关联方式,也可以使用现有其他数据关联方式进行量测与航迹的关联方法;
步骤3、根据关联结果更新航迹的状态:若航迹关联到量测,则基于卡尔曼滤波算法更新目标状态;否则,不对目标状态进行更新;
步骤4、根据关联结果更新航迹对应的滑窗:
若第i条航迹关联到量测,则更新k时刻的第i条航迹关联标志Mt(k,i)=1,i=1,2,…,m,m为航迹条数;
若第i条航迹关联未到量测,且目标的预测位置不在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=0;
若第i条航迹关联未到量测,且目标的预测位置在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=inf,inf为无穷大;
步骤5、查找1到k时刻第i条航迹的航迹关联标志Mt(:,i)中0或1元素的索引或者删除Mt(:,i)中inf的索引,根据提取出的索引或者删除后剩下的索引,就能提取出第i条航迹对应的有效滑窗向量M';
步骤6、基于伪滑窗L判决准则进行航迹删除,本发明的基于伪滑窗L判决是指在删除了伪滑窗的前提下只对有效滑窗进行的L判决:
若有效滑窗向量M'长度大于或等于航迹删除阈值L,则对第i条航迹继续维持;航迹删除阈值L通常取值为4到6;
若有效滑窗向量M'长度小于航迹删除阈值L,则删除第i条航迹;
步骤7、判断当前观测帧数k,是否满足k=K,如是,观测结束,否则令k=k+1,返回步骤2进行下一帧的观测,K为总的观测帧数。这里,观测时刻与观测帧数的起点、终点与计数方式相同,均用变量k表示。
本发明主要采用计算机仿真的方法进行验证,所有步骤、结论都在MATLAB-R2013b上验证正确。如图2所示,共有2个目标在一个二维平面上做匀速直线运动,目标被扫描到的时间间隔是不固定的。图3为100次蒙特卡洛试验下雷达跟踪到的目标个数,结果表明:传统航迹删除方法下,目标仅被跟踪到不足0.1条,而本发明达到1.6条,大大提高了成功跟踪概率。图4为跟踪精度对比,由于成功跟踪概率低,传统跟踪算法的均方根误差(RMSE)维持在最大跟踪误差附近;而本发明的均方根误差显然大大降低。其中,本发明的RMSE曲线起伏的原因是雷达连续多帧观测不到目标导致误差积累,检测到目标时,状态估计得到修正,于是误差骤然降低。
通过本发明具体实施方式可以看出,本发明可以很好的实现对抽样间隔随时间变化的目标的跟踪。同时,本发明还可以解决集中式跟踪系统多传感器异步通信时相同的问题。可以应用于水下被动定位声呐跟踪等领域。
Claims (1)
1.一种基于伪滑窗L判决准则的航迹删除方法,其特征在于,包括以下步骤:
步骤1、初始化参数,包括数据关联门限γ,航迹删除阈值L;
步骤2、当观测的帧数k达到航迹判决基数时,对第k帧量测zk与所有航迹进行数据关联;
步骤3、根据关联结果更新航迹的状态:若航迹关联到量测,则基于卡尔曼滤波算法更新目标状态;否则,不对目标状态进行更新;
步骤4、根据关联结果更新航迹对应的滑窗:
若第i条航迹关联到量测,则更新k时刻的第i条航迹关联标志Mt(k,i)=1,i=1,2,…,m,m为航迹条数;
若第i条航迹关联未到量测,且目标的预测位置不在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=0;Mt(k,i)为1或0时,均表示有效滑窗;
若第i条航迹关联未到量测,且目标的预测位置在当前雷达波束照射区域,则更新k时刻的第i条航迹关联标志Mt(k,i)=inf;inf表示无穷大,Mt(k,i)=inf时,表示伪滑窗;
步骤5、查找1到k时刻第i条航迹的航迹关联标志Mt(:,i)提取有效滑窗向量M';
步骤6、对有效滑窗向量M'根据L判决准则进行航迹删除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610801724.5A CN106405510B (zh) | 2016-09-05 | 2016-09-05 | 一种基于伪滑窗l判决准则的航迹删除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610801724.5A CN106405510B (zh) | 2016-09-05 | 2016-09-05 | 一种基于伪滑窗l判决准则的航迹删除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106405510A CN106405510A (zh) | 2017-02-15 |
CN106405510B true CN106405510B (zh) | 2019-01-11 |
Family
ID=57999458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610801724.5A Expired - Fee Related CN106405510B (zh) | 2016-09-05 | 2016-09-05 | 一种基于伪滑窗l判决准则的航迹删除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106405510B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291534B (zh) * | 2016-09-21 | 2018-07-27 | 电子科技大学 | 一种改进的航迹确认方法 |
CN109521419B (zh) * | 2017-09-20 | 2021-07-20 | 比亚迪股份有限公司 | 基于车辆雷达的目标跟踪方法和装置 |
CN108072863B (zh) * | 2017-12-15 | 2020-12-15 | 北京航空航天大学 | 一种星载sar发射时序编码方法及装置 |
CN110031834A (zh) * | 2018-01-12 | 2019-07-19 | 西安艾索信息技术有限公司 | 一种改进的多目标雷达航迹处理方法 |
CN110456341B (zh) * | 2019-09-11 | 2021-09-28 | 安徽隼波科技有限公司 | 一种基于双滑窗的雷达目标航迹质量管理方法 |
CN110888137B (zh) * | 2019-12-06 | 2022-07-15 | 西北工业大学 | 基于滑窗式双门限检测的角度跟踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614817A (zh) * | 2009-06-24 | 2009-12-30 | 北京航空航天大学 | 一种基于地面动目标指示雷达系统的多目标跟踪方法 |
CN101944234A (zh) * | 2010-07-23 | 2011-01-12 | 中国科学院研究生院 | 特征迹驱动的多目标跟踪方法及装置 |
CN103885057A (zh) * | 2014-03-20 | 2014-06-25 | 西安电子科技大学 | 自适应变滑窗多目标跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4865958B2 (ja) * | 2001-05-23 | 2012-02-01 | 株式会社トクホン | 鎮痛抗炎症局所作用型の貼付剤 |
-
2016
- 2016-09-05 CN CN201610801724.5A patent/CN106405510B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614817A (zh) * | 2009-06-24 | 2009-12-30 | 北京航空航天大学 | 一种基于地面动目标指示雷达系统的多目标跟踪方法 |
CN101944234A (zh) * | 2010-07-23 | 2011-01-12 | 中国科学院研究生院 | 特征迹驱动的多目标跟踪方法及装置 |
CN103885057A (zh) * | 2014-03-20 | 2014-06-25 | 西安电子科技大学 | 自适应变滑窗多目标跟踪方法 |
Non-Patent Citations (2)
Title |
---|
Architectures and algorithms for track association and fusion;Chee-Yee Chong et al.;《IEEE Aerospace and Electronic Systems Magazine》;20000131;第15卷(第1期);全文 * |
无源多基地雷达中一种综合航迹处理方法;吕小永等;《电子学报》;20160430;第44卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106405510A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405510B (zh) | 一种基于伪滑窗l判决准则的航迹删除方法 | |
CN108490410B (zh) | 一种两坐标雷达对海目标联合检测跟踪方法 | |
CN110927712B (zh) | 跟踪方法和装置 | |
CN106199584B (zh) | 一种基于量测存储的航迹起始方法 | |
CN109655822A (zh) | 一种改进的航迹起始方法 | |
CN106291534B (zh) | 一种改进的航迹确认方法 | |
CN109946671A (zh) | 一种基于双门限判决的水下机动弱目标检测跟踪方法 | |
CN107436434B (zh) | 基于双向多普勒估计的航迹起始方法 | |
CN106204641B (zh) | 一种基于多假设检验的动态规划检测前跟踪方法 | |
CN105866769B (zh) | 一种并行运算的多目标检测前跟踪方法 | |
CN110133612A (zh) | 一种基于跟踪反馈的扩展目标检测方法 | |
CN104569963A (zh) | 用于超宽带穿墙雷达运动目标一维检测与跟踪的方法 | |
CN107102293B (zh) | 基于滑窗累积密度估计的未知杂波无源协同定位方法 | |
Solonar et al. | General construction principles and performance features of trajectory processing by data from one radar data source | |
CN113126086A (zh) | 一种基于状态预测积累的生命探测雷达弱目标检测方法 | |
Bi et al. | Improved multi-target radar TBD algorithm | |
JP5623253B2 (ja) | 多目標追尾装置 | |
Feng et al. | Cardinality balanced multi-target multi-Bernoulli filter for target tracking with amplitude information | |
Chandrakanth et al. | Target tracking in blind range of radars with deep learning | |
Yong et al. | Radar track-before-detect algorithm of multitarget based on the dynamic programming | |
Xu et al. | Efficient centralized track initiation method for multistatic radar | |
Jiang | Multiple-model Bernoulli filters—Part I: A Gaussian mixture implementation | |
Zheng et al. | A fast visual tracking method via spatio-temporal context learning for unmanned rotorcrafts fixed-pointed landing | |
Zhang et al. | Improved Gaussian Mixture PHD for close multi-target tracking | |
Zhang et al. | Track Initiation Method Based on Deep Learning and Logic Method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190111 Termination date: 20210905 |
|
CF01 | Termination of patent right due to non-payment of annual fee |