CN114510676B - 一种评估连续时间航迹跟踪算法性能的度量方法 - Google Patents
一种评估连续时间航迹跟踪算法性能的度量方法 Download PDFInfo
- Publication number
- CN114510676B CN114510676B CN202210143954.2A CN202210143954A CN114510676B CN 114510676 B CN114510676 B CN 114510676B CN 202210143954 A CN202210143954 A CN 202210143954A CN 114510676 B CN114510676 B CN 114510676B
- Authority
- CN
- China
- Prior art keywords
- track
- time
- real
- measurement
- consistent
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 26
- 238000000691 measurement method Methods 0.000 title abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 35
- 238000012423 maintenance Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 19
- 230000001915 proofreading effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 101001086191 Borrelia burgdorferi Outer surface protein A Proteins 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
- Debugging And Monitoring (AREA)
Abstract
本发明公开了一种评估连续时间航迹跟踪算法性能的度量方法,通过一种“距离”度量基于时间轨迹函数(Trajectory Function of Time,T‑FOT)模型的真实航迹与量测航迹之间的差异。该方法旨在实时、快速、合理地体现出不同跟踪算法生成的量测航迹之间的优劣之分,从而选择出最优的跟踪算法实现对目标的实时跟踪。根据航迹维持时间的差异,将量测航迹划分为一致维持航迹与不一致偏差航迹两部分,在符合数学与运动学规律的前提下,通过对一致维持航迹部分的积分运算,叠加不一致偏差航迹部分的惩罚项,得到各量测航迹与真实航迹之间的差异值。通过上述方式,本发明能够对连续时间航迹跟踪算法实现实时的性能评估,且设计简便,计算复杂度小,计算效率高。
Description
技术领域
本发明属于目标跟踪技术领域,具体涉及一种航迹跟踪算法性能的度量方法。
背景技术
航迹估计属于目标跟踪研究的核心问题,目前国内外主流算法实现的是离散时间航迹估计,因其易受机动目标、未知输入等因素影响,一个替代方案——连续时间航迹估计已成为相关领域学者的研究热点。评估目标跟踪算法性能是否良好的关键在于度量估计航迹和真实航迹之间的“距离”,广受认可的综合度量方法包括最优子模式分配(optimalsub-pattern assignment,OSPA)、OSPA2、广义最优子模式分配(generalized optimalsubpattern assignment,GOSPA)等,但以上方法均只适用于对离散时间航迹的度量,在连续时间航迹估计的研究趋势下,一种能够度量连续航迹跟踪性能的方法至关重要。
文献“A Complete Optimal Subpattern Assignment(COSPA)Metric,Fusion2020 in July 6-9,2020”提出了一种称为完整最优子模式分配的度量方法。该方法是基于OSPA度量的改进方法,通过对距离误差和基数误差的阈值进行单独控制,达到了较全面地评估目标跟踪算法整体性能的效果。此外,选择小于或等于每个不可分配点的惩罚的截止值,能够克服OSPA度量对两个有限点集存在一空集的情况不敏感的局限性,同时保留其用于评估多目标滤波算法性能的所有优点。然而,文献所述方法是基于两个有限点集的度量方法,仅在评估离散时间航迹领域有一定效用,并不适用于实现连续时间航迹跟踪算法性能的评估。倘若量测航迹为平滑分段多项式拟合函数,该方法则毫无效用,因此在该领域的相关评估算法有待进一步扩展。
发明内容
为了克服现有技术的不足,本发明提供了一种评估连续时间航迹跟踪算法性能的度量方法,通过一种“距离”度量基于时间轨迹函数(Trajectory Function of Time,T-FOT)模型的真实航迹与量测航迹之间的差异。该方法旨在实时、快速、合理地体现出不同跟踪算法生成的量测航迹之间的优劣之分,从而选择出最优的跟踪算法实现对目标的实时跟踪。根据航迹维持时间的差异,将量测航迹划分为一致维持航迹与不一致偏差航迹两部分,在符合数学与运动学规律的前提下,通过对一致维持航迹部分的积分运算,叠加不一致偏差航迹部分的惩罚项,得到各量测航迹与真实航迹之间的差异值。通过上述方式,本发明能够对连续时间航迹跟踪算法实现实时的性能评估,且设计简便,计算复杂度小,计算效率高。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤1:将真实航迹与量测航迹建模为T-FOT形式:
xt=F(t;λ)+et
式中,xt表示t时刻目标的状态,λ表示连续时间函数曲线系数,F(t;λ)表示量测航迹连续时间函数曲线,et表示F(t;λ)与真实航迹连续时间函数曲线f(t)的误差;
步骤2:航迹维持时间校对;
对比量测航迹与真实航迹的起始时间和终止时间,将量测航迹与真实航迹时间序列中重合时域部分Tconsistent的航迹定义为一致维持航迹Φconsistent,重合时域Tconsistent以外的航迹定义为不一致偏差航迹Φerror;
步骤3:对于一致维持航迹Φconsistent,计算量测航迹与真实航迹之间的一致维持航迹位置误差E:
步骤4:对于不一致偏差航迹Φerror,定义平均航迹误差与不一致偏差航迹位置误差C分别为:
其中,Terror表示非重合时域;
步骤5:将一致维持航迹位置误差E与不一致偏差航迹位置误差C叠加,得到最终的连续时间航迹度量值:D=E+C,作为评估连续时间航迹跟踪算法性能的指标。
本发明的有益效果如下:
本发明提出一种评估连续时间航迹跟踪算法性能的度量方法,真实航迹与不同跟踪算法生成的量测航迹均基于T-FOT模型建立为连续时间航迹拟合曲线,主要任务是衡量真实航迹与量测航迹之间的差异,并以一种“距离”的方式表达出来,能够实时、快速、合理地确定出与真实航迹最为相似的量测航迹。此方法主要实现了连续时间航迹跟踪算法性能的评估,更加贴合目前相关领域的发展趋势。
附图说明
图1是本发明方法的流程图。
图2是本发明方法的原理图。
图3是本发明实施例的仿真场景。
图4是本发明实施例的仿真结果。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图1所示,一种评估连续时间航迹跟踪算法性能的度量方法,包括如下步骤:
步骤1:将真实航迹与量测航迹建模为T-FOT形式:
xt=F(t;λ)+et
式中,xt表示t时刻目标的状态,λ表示连续时间函数曲线系数,F(t;λ)表示量测航迹连续时间函数曲线,et表示F(t;λ)与真实航迹连续时间函数曲线f(t)的误差;
步骤2:航迹维持时间校对;
对比量测航迹与真实航迹的起始时间和终止时间,将量测航迹与真实航迹时间序列中重合时域部分Tconsistent的航迹定义为一致维持航迹Φconsistent,重合时域Tconsistent以外的航迹定义为不一致偏差航迹Φerror;
步骤3:对于一致维持航迹Φconsistent,计算量测航迹与真实航迹之间的一致维持航迹位置误差E:
步骤4:对于不一致偏差航迹Φerror,定义平均航迹误差与不一致偏差航迹位置误差C分别为:
步骤5:将一致维持航迹位置误差E与不一致偏差航迹位置误差C叠加,得到最终的连续时间航迹度量值:D=E+C,作为评估连续时间航迹跟踪算法性能的指标。
如图2所示,实线为真实航迹,虚线为量测航迹,可观察到两航迹时域不重合。需要对斜线阴影部分(重合时域航迹)计算一致维持航迹位置误差E,对网格阴影部分(非重合时域航迹)计算得到不一致偏差航迹位置误差C,两误差值相加即可得到最终的连续时间航迹度量值D。
具体实施例:
1、设置仿真场景。
参照图3,仿真场景设置为一维坐标系下,目标随时间发生位置移动,真实轨迹是以时间t为自变量的一维曲线。不同算法对真实轨迹进行估计,拟合出三条不同的量测航迹,本实施例主要是为了衡量真实航迹与量测航迹之间的距离误差,因此所有航迹曲线均为预设好的拟合函数以便说明本发明的核心步骤。真实航迹与三条量测航迹曲线函数如下所示:
xtrue=0.2t2+0.3t+1,t∈(-1,4) (1)
x1=0.2t2+0.1t+4,t∈(-1,4.5) (2)
x2=0.3t2+0.1t,t∈(-2,3) (3)
x3=0.3t2+2,t∈(-1.5,5) (4)
其中,xtrue为真实航迹拟合函数,x1、x2、x3分别为量测航迹1、量测航迹2、量测航迹3的拟合函数。
2、根据航迹时域一致的区间计算一致维持航迹Φconsistent的位置误差。
三条量测航迹起始时间和终止时间均不相同,需要找出与真实航迹时间对应的航迹区间,由于仿真场景以t为自变量,可通过简单的时间对比得到一致时间航迹区间。通过对一致航迹区间进行定积分,计算出量测航迹与真实航迹的一致航迹位置误差,其计算公式如下:
式中,Tconsistent表示量测航迹与真实航迹的重合时域。
3、计算不一致偏差航迹Φerror的位置误差。
各量测轨迹的不一致时间航迹区间不同,本实施例通过对一致时间航迹区间位置误差进行关于时间的均值计算,所得结果作为各量测航迹的惩罚项叠加到距离误差中,能够做到合理且均匀,三条量测航迹的惩罚项分别为:
式中,Terrori为各量测航迹与真实航迹时间不一致的航迹区间,为各量测航迹的平均航迹误差。
惩罚叠加后各量测航迹误差值分别为:
D1=E1+C1 (11)
D2=E2+C2 (12)
D3=E3+C3 (13)
4、对比总位置误差值确定最佳量测航迹。
量测航迹位置误差值越接近0,其与真实航迹的相似度越高。量测航迹1、2、3的距离误差值如图4所示,可明显看出量测航迹3的分数最接近于0,因此该量测航迹与真实航迹最为相似。
Claims (1)
1.一种评估连续时间航迹跟踪算法性能的度量方法,其特征在于,包括如下步骤:
步骤1:将真实航迹与量测航迹建模为T-FOT形式:
xt=F(t;λ)+et
式中,xt表示t时刻目标的状态,λ表示连续时间函数曲线系数,F(t;λ)表示量测航迹连续时间函数曲线,et表示F(t;λ)与真实航迹连续时间函数曲线f(t)的误差;
步骤2:航迹维持时间校对;
对比量测航迹与真实航迹的起始时间和终止时间,将量测航迹与真实航迹时间序列中重合时域部分Tconsistent的航迹定义为一致维持航迹Φconsistent,重合时域Tconsistent以外的航迹定义为不一致偏差航迹Φerror;
步骤3:对于一致维持航迹Φconsistent,计算量测航迹与真实航迹之间的一致维持航迹位置误差E:
步骤4:对于不一致偏差航迹Φerror,定义平均航迹误差与不一致偏差航迹位置误差C分别为:
其中,Terror表示非重合时域;
步骤5:将一致维持航迹位置误差E与不一致偏差航迹位置误差C叠加,得到最终的连续时间航迹度量值:D=E+C,作为评估连续时间航迹跟踪算法性能的指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210143954.2A CN114510676B (zh) | 2022-02-17 | 2022-02-17 | 一种评估连续时间航迹跟踪算法性能的度量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210143954.2A CN114510676B (zh) | 2022-02-17 | 2022-02-17 | 一种评估连续时间航迹跟踪算法性能的度量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114510676A CN114510676A (zh) | 2022-05-17 |
CN114510676B true CN114510676B (zh) | 2024-07-02 |
Family
ID=81551751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210143954.2A Active CN114510676B (zh) | 2022-02-17 | 2022-02-17 | 一种评估连续时间航迹跟踪算法性能的度量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114510676B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116910577B (zh) * | 2023-05-09 | 2024-02-13 | 中国人民解放军91404部队第340所 | 航空兵蓝军模拟战术运用相似度评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977022A (zh) * | 2014-04-04 | 2015-10-14 | 西北工业大学 | 多目标跟踪系统性能评估仿真方法 |
CN110990756A (zh) * | 2019-12-03 | 2020-04-10 | 中国电子科技集团公司第二十八研究所 | 基于可拓型逐层分析的空中目标航迹质量评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3179355B2 (ja) * | 1996-11-07 | 2001-06-25 | 三菱電機株式会社 | 多目標追尾方法及び多目標追尾装置 |
JP5107282B2 (ja) * | 2009-02-27 | 2012-12-26 | 株式会社東芝 | 目標追尾装置 |
CN111896946B (zh) * | 2020-06-29 | 2023-08-18 | 西北工业大学 | 一种基于航迹拟合的连续时间目标跟踪方法 |
CN113008222B (zh) * | 2021-02-20 | 2023-03-31 | 西北工业大学 | 一种基于连续时间轨迹函数的航迹约束目标跟踪方法 |
-
2022
- 2022-02-17 CN CN202210143954.2A patent/CN114510676B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977022A (zh) * | 2014-04-04 | 2015-10-14 | 西北工业大学 | 多目标跟踪系统性能评估仿真方法 |
CN110990756A (zh) * | 2019-12-03 | 2020-04-10 | 中国电子科技集团公司第二十八研究所 | 基于可拓型逐层分析的空中目标航迹质量评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114510676A (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114510676B (zh) | 一种评估连续时间航迹跟踪算法性能的度量方法 | |
CN107133406B (zh) | 一种电力系统静态电压稳定域边界的快速搜索方法 | |
CN113312413B (zh) | 基于电力大数据的企业碳排放异常监测方法 | |
CN110286332A (zh) | 一种基于多新息理论的电动汽车动力电池soc估计方法 | |
CN112584311B (zh) | 一种基于wknn融合的室内三维空间指纹定位方法 | |
CN106958643B (zh) | 齿轮轮齿接触性能分析方法 | |
CN111427007B (zh) | 基于不完全量测下集员滤波的矿井人员安全状态估计方法 | |
CN114186476A (zh) | 基于历史工况的锅炉燃烧在线优化方法 | |
CN111553060B (zh) | 基于类噪声扰动的电网惯量常态化连续估计方法 | |
CN110231620A (zh) | 一种噪声相关系统跟踪滤波方法 | |
CN116995653A (zh) | 一种基于ddtw距离的低压台区相户拓扑识别方法 | |
Yang et al. | Research on indoor location algorithm based on WIFI | |
CN111144705B (zh) | 一种基于带时标信息采集的全网同断面数据处理方法 | |
CN105761489A (zh) | 一种交通流预测的三次指数平滑最优方法 | |
CN117380424A (zh) | 一种数字化喷涂优化方法及系统 | |
CN115603845B (zh) | 一种基于一致性策略的最优分布式滤波器协同目标跟踪方法 | |
CN111178528A (zh) | 一种应用于无线充电系统的精英遗传算法改进方法 | |
CN115859626A (zh) | 针对周期运动目标的自适应无迹卡尔曼滤波器设计方法 | |
CN115034040A (zh) | 一种多情景重现期暴雨雨型的演示方法 | |
CN114091341A (zh) | 一种基于igwo-bp的变压器故障诊断方法 | |
CN112084455B (zh) | 一种s参数的最小二乘插值方法 | |
CN109408983B (zh) | 一种基于高程的路基压实自动分层方法 | |
Lu et al. | An approximation method combined with conformal transformation for the fast measurement of blade section curves | |
Djuric et al. | Two-dimensional simulation of liquid metal spray deposition onto a complex surface | |
CN111145537B (zh) | 一种出行生成量预测方法及系统 |
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 |