CN109508444B - 区间量测下交互式多模广义标签多伯努利的快速跟踪方法 - Google Patents
区间量测下交互式多模广义标签多伯努利的快速跟踪方法 Download PDFInfo
- Publication number
- CN109508444B CN109508444B CN201811549885.5A CN201811549885A CN109508444B CN 109508444 B CN109508444 B CN 109508444B CN 201811549885 A CN201811549885 A CN 201811549885A CN 109508444 B CN109508444 B CN 109508444B
- Authority
- CN
- China
- Prior art keywords
- target
- particle
- moment
- particles
- label
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005259 measurement Methods 0.000 title claims abstract description 40
- 230000002452 interceptive effect Effects 0.000 title claims abstract description 28
- 239000002245 particle Substances 0.000 claims abstract description 164
- 238000005070 sampling Methods 0.000 claims abstract description 72
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000012546 transfer Methods 0.000 claims abstract description 5
- 238000009826 distribution Methods 0.000 claims description 56
- 230000033001 locomotion Effects 0.000 claims description 16
- 230000004083 survival effect Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 4
- 230000002085 persistent effect Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 37
- 238000001914 filtration Methods 0.000 abstract description 12
- 238000001514 detection method Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000003993 interaction Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
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
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Analysis (AREA)
Abstract
本发明公开一种区间量测下交互式多模广义标签多伯努利的快速跟踪方法,将交互式多模方法与快速算法思想相结合,首先在广义标签多伯努利滤波的框架下,针对目标采样粒子预测阶段,结合区间量测广义似然函数,实现所有粒子对于不同模型的转移预测,随后通过计算模型权概率对粒子进行模型交互,然后通过GLMB滤波更新方程对模型交互后的粒子进行更新。在此基础上结合快速实现方法,将预测与更新相结合,对于每个迭代只需要一个截断过程,降低了算法的计算量,最终解决了机动多目标的检测与跟踪问题。
Description
技术领域
本发明涉及目标跟踪技术领域,具体涉及一种区间量测下交互式多模广义标签多伯努利的快速跟踪方法。
背景技术
非机动目标运动可以用一个固定的模型来描述,但要描述机动目标的运动,可能需要结合具有不同机动特性的运动模型。随着机动目标跟踪技术受到越来越广泛的关注,对机动目标跟踪技术的要求也越来越高。多机动目标跟踪已经成为目标跟踪领域中一个极其困难的问题。
姬红兵在其发表的论文“箱粒子广义标签多伯努利滤波的目标跟踪算法”(西安交通大学学报,2017,51(10):107-112.)中提出了一种箱粒子广义标签多伯努利跟踪算法。该算法利用箱粒子滤波算法近似单目标状态的概率密度,即用一组带权值的均匀分布拟合单目标状态概率密度;最后通过广义标签多伯努利滤波对多目标状态的概率密度进行预测与更新,从更新后的多目标状态概率密度中估计单目标的位置与速度,并且由于单目标的标签互不相同可以实现航迹跟踪。该算法的不足之处在于,无法对强机动目标进行有效跟踪。
Vo等人在其发表的论文“A Generalized Labeled Multi-Bernoulli Filter for Maneuvering Targets”(19th International Conference on Information Fusi on)中通过结合交互式多模概念和标签多伯努利随机有限集(RFS)理论,提出了一种广义标签多伯努利粒子滤波(GLMB)的机动目标跟踪算法,并给出了高斯混合(GM)的实现形式。该算法可以在估计多机动目标状态的同时跟踪不同目标的航迹。该算法的不足之处在于,由于强杂波和多目标环境下的组合***问题,GLMB滤波器实现面临着巨大的计算复杂性。
Vo等人在其发表的论文“An Efficient Implementation of the GeneralizedLabeled Multi-Bernoulli Filter”(IEEE Transactions on Signal Processing,2017,65(8):1975-1987.)中提出了一种基于吉布斯采样的GLMB滤波密度截断算法,将预测和更新集成到一个步骤中,实现了GLMB滤波的高效实现。该算法的不足之处在于,在许多实际应用中,标准量测模型是不够的。虽然传感器检测报告是一个点测量值,但是实际测量会受到未知边界误差分布的影响,使得这种非标准测量需要以区间量测的形式进行。
发明内容
本发明所要解决的是现有多机动目标跟踪方法在区间量测下目标发生较大机动时的目标失跟的问题,提供一种区间量测下交互式多模广义标签多伯努利的快速跟踪方法,其在广义标签多伯努利滤波的基础上结合交互式多模方法和快速算法思想,采用序贯蒙特卡洛(SMC)算法,实现对机动多目标的精准跟踪。
本发明的基本思路是:将交互式多模方法与快速算法思想相结合,首先在广义标签多伯努利滤波的框架下,针对目标采样粒子预测阶段,结合区间量测广义似然函数,实现所有粒子对于不同模型的转移预测,随后通过计算模型权概率对粒子进行模型交互,然后通过GLMB滤波更新方程对模型交互后的粒子进行更新。在此基础上结合快速实现方法,将预测与更新相结合,对于每个迭代只需要一个截断过程,降低了算法的计算量,最终解决了机动多目标的检测与跟踪问题。
为解决上述问题,本发明是通过以下技术方案实现的:
区间量测下交互式多模广义标签多伯努利的快速跟踪方法,包括步骤如下:
步骤1、根据目标运动场景,设定初始时刻目标粒子的状态参数,并用设定的状态参数作为目标粒子的初始分布,采样固定数目的初始分布目标粒子,并将其用标签多伯努利随机集的参数集形式表示,得到初始时刻目标采样粒子标签多伯努利随机集随机集的后验分布;
步骤2、利用前一时刻的目标采样粒子标签多伯努利随机集的后验分布和前一时刻的区间量测数据,用交互式多模方法对目标采样粒子进行预测,得到当前时刻预测的目标采样粒子标签多伯努利随机集的后验分布;
步骤3、利用当前时刻的区间量测数据,计算每个目标采样粒子的广义似然函数,并根据广义似然函数,使用广义标签多伯努利滤波器对当前时刻预测的目标采样粒子标签多伯努利随机集的后验分布进行更新,得到当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布;
步骤4、从步骤4所得到的当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布中,选出其标签权重大于给定阈值的当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布;
步骤5、分别计算步骤5所得到的当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布的势,并找出其中最大势所对应的指标N;从步骤5所得到的当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布中选出标签权重较大的N-1个当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻最终的目标采样粒子标签多伯努利随机集的后验分布;
步骤6、计算当前时刻最终的目标采样粒子标签多伯努利随机集的后验分布的加权和,并将该加权和结果作为当前时刻估计的目标状态;
步骤7、判断所有时刻是否处理完毕:若是,则输出当前时刻估计的目标状态;否则,执行步骤2,处理下一时刻。
与现有技术相比,本发明具有如下特点:
1、将多模型概念和快速算法思想嵌入到广义标签多伯努利算法中,采用粒子滤波方法,不但能精确的跟踪目标的航迹和估计目标的数目,在算法的计算量上也远远小于交互式多模广义标签多伯努利算法,提高了运算性能。
2、结合了交互式多模,针对目标采样粒子预测阶段,实现所有粒子对于不同模型的转移预测,随后通过计算模型权概率对粒子进行模型交互,实现了对机动目标的检测跟踪,克服了现有技术跟踪机动目标存在的失跟问题。
3、结合了快速算法思想,目标量测模型不采用传统传感器量测模型,传感器输出区间量测,区间量测以弥补任何点量测缺陷,全面的考虑到了机动目标的机动性不稳定等因素,快速算法整合预测和更新到一个步骤,对于每个迭代只需要一个截断过程,大大的降低了算法的计算量,有效的实现了区间量测下机动多目标的检测与跟踪。
附图说明
图1为区间量测下交互式多模广义标签多伯努利的快速跟踪方法的流程图。
图2为本实施例的仿真图,其中(a)是目标运动轨迹图,(b)是目标航迹跟踪结果图,(c)是目标数目跟踪效果图,(d)是目标数目跟踪误差图,(e)是目标OSPA距离跟踪误差图,(f)是目标OSPA位置跟踪误差图,(g)是目标OSPA势跟踪误差图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实例,并参照附图,对本发明进一步详细说明。
针对区间量测的不确定性和目标机动性不确定的问题,本发明在广义标签多伯努利滤波器的基础上,使用交互式多模型方法对滤波器中每个目标状态的采样粒子进行预测,再将预测粒子通过引入广义似然函数结合GLMB滤波器更新策略来更新预测粒子。在此基础上结合快速实现方法,将预测与更新相结合,对于每个迭代只需要一个截断过程,降低了算法的计算量,本发明在区间量测环境下,能有效地检测跟踪机动多目标,且对目标状态和数目的估计更加精准。
一种区间量测下交互式多模广义标签多伯努利的快速跟踪方法,其具体步骤如下:
步骤1,初始化目标状态。
根据目标运动场景,设定初始时刻GLMB轨迹表,包括目标粒子(包括持续存活粒子和新生粒子)和状态参数,即目标位置、速度、权重、模型概率、轨迹标签以及轨迹关联历史。初始化GLMB分量假设权重、假设标签、假设势以及势分布。用上述设定的参数作为目标的初始分布,用高斯分布采样固定数目的初始存活粒子及新生粒子,并用标签多伯努利随机集的参数集形式表示。
初始粒子采样过程如下:
令初始时刻k=0,目标初始分布用高斯粒子标签多伯努利随机集的参数集表示形式如下:
其中,H是指标空间,h是指标,表示初始时刻GLMB轨迹表假设标签的假设权重;表示初始时刻第j个目标采样粒子(标签为l)的模型权概率,p为模型个数;表示初始时刻第j个目标采样粒子的状态,x0表示初始时刻目标的横坐标,表示初始时刻目标的水平速度,y0表示初始时刻目标的纵坐标,表示初始时刻目标的垂直速度,T表示转置操作;表示初始时刻第j个目标采样粒子状态对应的状态权重;表示初始时刻的目标采样粒子数目。
步骤2,预测和更新目标状态。
xk=fk-1(xk-1,sk)+vk-1(sk)
其中,fk-1表示k时刻目标的状态转移方程,sk表示k时刻的模型变量,vk-1表示状态噪声。
假定k-1时刻的目标采样粒子标签多伯努利随机集后验分布为:
则持续存活粒子标签多伯努利随机集表示为:
表示从k-1时刻到k时刻持续存活粒子在模型s下的状态预测,表示k时刻新生目标的采样粒子在模型s下的状态预测,表示k时刻持续存活粒子在模型s下的模型权概率,表示新生目标的采样粒子k时刻在模型s下的模型权概率,表示从k-1时刻到k时刻持续存活粒子在模型s下的权重预测,表示k时刻新生目标的采样粒子在模型s下的粒子权重预测。
具体预测方法可由下列步骤来完成。
第2.1步,对k-1时刻更新后的目标采样粒子标签多伯努利随机集后验分布,进行重采样得到k-1时刻存活粒子的采样样本为:
其中N(mi,P),i=1,2,3,4是目标k时刻的新生密度,具体过程表示如下:
其中,Fs,k为模型s对应的状态转移方程,s=1,…,p,p为模型总个数,νk为状态噪声。
模型权概率计算:
第2.4步,计算存活目标和出生目标的预测状态和权重:
第2.5步,合并带标签的权重粒子:
第2.6步,更新目标状态。
假设k时刻预测的目标采样粒子表示为:
则更新的目标采样粒子后验分布为:
其中,表示从k-1时刻到k时刻第j个目标采样粒子状态预测,表示k时刻第j个目标采样粒子状态的更新;表示从k-1时刻到k时刻第j个目标采样粒子的权重的预测,表示k时刻第j个目标采样粒子的权重的更新;表示k时刻预测的目标采样粒子数目,表k时刻的目标采样粒子数目。
利用k时刻目标随机集的广义似然函数,更新k时刻预测完成的目标采样粒子标签多伯努利随机集,得到k时刻目标采样粒子标签多伯努利随机集的后验分布。
具体的更新方法可由下列步骤来完成。
第2.6.1步,利用当前时刻的区间量测数据,计算每个预测目标采样粒子对应的广义似然函数,假设,考虑到在Xk上,检测是独立的,并且杂波与检测无关,多目标似然具体计算如下:
第2.6.2步,根据区间量测似然函数,对k时刻每个预测目标粒子的权重进行更新计算:
步骤3,分量截断。
从当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布中,选出其标签权重大于给定阈值的当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布。
通常情况下,GLMB算法中,前一时刻的假设轨迹的个数到下一时刻会被截断两次,分别是假设轨迹的预测和更新,都需要截断操作。为了提高算法运算效率,本发明算法整合预测和更新步骤,仅需一次截断,大大的提高了算法的运行速度。
假定k-1时刻的目标采样粒子标签多伯努利随机集后验分布为:
则k时刻的目标采样粒子标签多伯努利随机集后验分布为:
其中H是指标空间,h是指标,从k-1时刻分量总数为Hk-1,到k时刻分量总数为Hk,只需要一个截断操作,大大降低了本发明算法的计算复杂度。
步骤4,状态估计。
我们使用边缘多目标估计的次优版本,最大势的后验估计。多目标状态的平均估计取决于估计的势,对更新后的目标粒子计算势分布。
分别计算当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布的势,并找出其中最大势所对应的指标Nk+1;从当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布中选出标签权重较大的Nk个当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻最终的目标采样粒子标签多伯努利随机集的后验分布。
根据当前时刻目标采样粒子标签多伯努利随机集的后验分布,采用加权法计算出目标采样粒子状态的加权和,作为当前时刻真实存在的目标状态。若此时估计出目标后验最大势对应的目标个数为Nk。状态估计具体计算如下:
步骤5,判断所有时刻是否处理完毕,若是,执行步骤6,否则,执行步骤2,处理下一时刻。
步骤6,结束。
本发明在广义标签多伯努利滤波的基础上,使用交互式多模型方法对滤波器中每个目标状态的采样粒子进行预测,再将预测粒子代入到GLMB的算法中进行目标存在概率及分布密度的更新估计。本发明在区间量测环境情况下,能有效地检测跟踪机动多目标,实现对目标状态和数目的估计。
下面结合附图2的仿真图,对本发明的效果做进一步说明。
仿真条件:本发明在Intel(R)Core(TM)i3-2370M [email protected]处理器的电脑上,采用MATLAB R2014a软件完成仿真。
仿真场景设置:为了验证本发明提出的一种区间量测下交互式多模广义标记多伯努利滤波器的快速实现算法能精确的检测和跟踪弱小机动目标,本发明的仿真实验场景为[-2000,2000]×[0,2000]m2二维空间内,整个仿真过程共持续100秒,共五个目标,每个目标的出生死亡时间以及运动状态如表1所述。
表1
目标 | 匀速直线运动 | 左转弯运动 | 右转弯运动 |
1 | 1~20s | 21~65s | 66~100s |
2 | 10~25s | 26~70s | 71~100s |
3 | 15~35s | 36~65s | 66~100s |
4 | 20~35s | 36~57s | 58~80s |
5 | 60~70s | 71~95s | 96~100s |
目标状态方程为:
量测方程为:
其中,表示量测似然函数,(xk,yk)表示目标的位置。其中ωk为零均值且协方差矩阵为的高斯白噪声,其中σr=2.5m和σθ=0.25°。这种传感器提供区间量测,具有区间长度为Δ=[Δr,Δθ]T,其中Δr=50m和Δθ=4°分别是范围和方位角的区间长度。
广义标签多伯努利滤波器初始化时给定五个目标的准确位置,目标初始状态为x1=[1000,-10,1500,-10],x2=[-250,20,1000,3],x3=[-1500,11,250,10],x4=[250,11,750,5],x5=[1000,-50,1500,0]。相关仿真参数设置为:pD,k=0.98和pS,k=0.99,杂波均值λ=10。目标初始模型权概率
为了证明仿真效果,截断阈值Hth=1×10-15,在杂波率为λ=10的情况下,存活及新生粒子数目按照100个粒子进行100次蒙特卡洛实验。计算目标的OSPA(Optimal Sub-Pattern Assignment)距离,OSPA参数设置为c=100m,p=2。图2(a)和图2(b)分别表示杂波率为λ=10的情况下,目标的运动状态和航迹跟踪结果。图2(a)轨迹中的圆圈表示目标初始位置,三角形表示目标最终位置。图2(b)是目标的航迹跟踪结果;由图2(b)可知,本发明的方法能够处理目标机动问题,在目标发生转弯的时刻也能跟踪出目标的准确位置。图2(c)和图2(d)分别是目标势的跟踪效果图和势误差图;由图2(c)和图2(d)可知,在不考虑个别时刻情况下,本发明的方法表现出了相对优势,其对目标势的估计相对更准确一些,估计误差也相对稳定且小于IMM-GLMB算法(区间量测下交互式多模广义标签多伯努利滤波器)和IMM-LMB算法(区间量测下交互式多模标签多伯努利滤波器)。图2(e)、图2(f)和图2(g)分别表示存活及新生粒子数目均为100个的目标OSPA距离跟踪误差图、OSPA位置跟踪误差图和OSPA势跟踪误差图;由图2(e)可知,当目标个数增多的时候,本发明的方法对目标距离跟踪误差相对稳定且小于IMM-GLMB算法和IMM-LMB算法;由图2(f)可知,本发明的方法在估计目标位置时,跟踪性能得到了进一步的提升,减小了目标发生机动时刻的估计偏差,滤波性能较IMM-GLMB算法和IMM-LMB算法两种方法而言更加稳定;由图2(g)可知,本发明的方法在估计目标数目时,跟踪性能明显比IMM-GLMB算法和IMM-LMB算法要好,在目标个数增加时,也能够很好的追踪出目标的个数。
综上所述,从仿真效果图的分析可知,本发明提出的一种区间量测下交互式多模广义标签多伯努利滤波器的快速实现跟踪方法,实现了对区间量测下机动多目标的检测和跟踪。目标跟踪精度高,跟踪性能良好,其性能相对优于交互式多模广义标签多伯努利滤波器、交互式多模标签多伯努利滤波器、交互式多模概率假设密度滤波器和交互式多模势概率假设密度滤波器,本发明的方法在处理机动目标跟踪问题时,在滤波性能上有明显优势,在现实工程应用中更有优势。
需要说明的是,尽管以上本发明所述的实施例是说明性的,但这并非是对本发明的限制,因此本发明并不局限于上述具体实施方式中。在不脱离本发明原理的情况下,凡是本领域技术人员在本发明的启示下获得的其它实施方式,均视为在本发明的保护之内。
Claims (1)
1.区间量测下交互式多模广义标签多伯努利的快速跟踪方法,其特征是,包括步骤如下:
步骤1、根据目标运动场景,设定初始时刻目标粒子的状态参数,并用设定的状态参数作为目标粒子的初始分布,采样固定数目的初始分布目标粒子,并将其用标签多伯努利随机集的参数集形式表示,得到初始时刻目标采样粒子标签多伯努利随机集的后验分布为:
式中,H是指标空间,h是指标,表示初始时刻广义标签多伯努利粒子滤波轨迹表假设标签的假设权重;表示初始时刻第j个目标采样粒子的状态,x0表示初始时刻目标的横坐标,表示初始时刻目标的水平速度,y0表示初始时刻目标的纵坐标,表示初始时刻目标的垂直速度,T表示转置操作;表示初始时刻第j个目标采样粒子的模型权概率,p为模型个数;表示初始时刻第j个目标采样粒子状态对应的状态权重;表示初始时刻的目标采样粒子数目;l为标签;
步骤2、利用前一时刻的目标采样粒子标签多伯努利随机集的后验分布和前一时刻的区间量测数据,用交互式多模方法对目标采样粒子进行预测,得到当前时刻预测的目标采样粒子标签多伯努利随机集的后验分布为:
式中,表示从k-1时刻到k时刻持续存活粒子在模型s下的状态预测,表示k时刻持续存活粒子在模型s下的模型权概率,表示从k-1时刻到k时刻持续存活粒子在模型s下的权重预测,表示持续存活粒子数目;表示k时刻新生目标的采样粒子在模型s下的状态预测,表示新生目标的采样粒子k时刻在模型s下的模型权概率,表示k时刻新生目标的采样粒子在模型s下的粒子权重预测;表示新生粒子数目;l为标签;
具体预测方法可由下列步骤来完成:
步骤2.1、对k-1时刻更新后的目标采样粒子标签多伯努利随机集后验分布,进行重采样得到k-1时刻存活粒子的采样样本;
步骤2.4、计算存活目标和出生目标的预测状态和权重:
步骤2.5、合并带标签的权重粒子:
步骤3、利用当前时刻的区间量测数据,计算每个目标采样粒子的广义似然函数,并根据广义似然函数,使用广义标签多伯努利滤波器对当前时刻预测的目标采样粒子标签多伯努利随机集的后验分布进行更新,得到当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布为:
具体的更新方法可由下列步骤来完成:
步骤3.1、利用当前时刻的区间量测数据,计算每个预测目标采样粒子对应的广义似然函数,多目标似然具体计算如下:
步骤3.2、根据区间量测似然函数,对k时刻每个预测目标粒子的权重进行更新计算:
步骤4、从步骤4所得到的当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布中,选出其标签权重大于给定阈值的当前时刻更新的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布为:
式中,Hk是指标空间H在k时刻分量总数,h是指标,;表示k时刻广义标签多伯努利粒子滤波轨迹表假设标签的假设权重,表示k时刻第j个目标采样粒子的状态,表示k时刻第j个目标采样粒子的模型权概率,表示k时刻第j个目标采样粒子状态对应的状态权重;表示k时刻的目标采样粒子数目;l为标签;
步骤5、分别计算步骤5所得到的当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布的势,并找出其中最大势所对应的指标N;从步骤5所得到的当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布中选出标签权重较大的N-1个当前时刻截断的目标采样粒子标签多伯努利随机集的后验分布,作为当前时刻最终的目标采样粒子标签多伯努利随机集的后验分布;
步骤6、计算当前时刻最终的目标采样粒子标签多伯努利随机集的后验分布的加权和,并将该加权和结果作为当前时刻估计的目标状态;状态估计具体计算如下:
步骤7、判断所有时刻是否处理完毕:若是,则输出当前时刻估计的目标状态;否则,执行步骤2,处理下一时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811549885.5A CN109508444B (zh) | 2018-12-18 | 2018-12-18 | 区间量测下交互式多模广义标签多伯努利的快速跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811549885.5A CN109508444B (zh) | 2018-12-18 | 2018-12-18 | 区间量测下交互式多模广义标签多伯努利的快速跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109508444A CN109508444A (zh) | 2019-03-22 |
CN109508444B true CN109508444B (zh) | 2022-11-04 |
Family
ID=65753547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811549885.5A Active CN109508444B (zh) | 2018-12-18 | 2018-12-18 | 区间量测下交互式多模广义标签多伯努利的快速跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109508444B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11094070B2 (en) * | 2019-04-23 | 2021-08-17 | Jiangnan University | Visual multi-object tracking based on multi-Bernoulli filter with YOLOv3 detection |
CN110084831B (zh) * | 2019-04-23 | 2021-08-24 | 江南大学 | 基于YOLOv3多伯努利视频多目标检测跟踪方法 |
CN110532921B (zh) * | 2019-08-21 | 2022-01-07 | 江南大学 | 基于ssd检测广义标签多伯努利视频多目标跟踪方法 |
CN110967690B (zh) * | 2019-11-12 | 2023-06-16 | 江南大学 | 一种基于多伯努利分布式多传感器多目标跟踪方法 |
CN113409363B (zh) * | 2020-03-17 | 2024-04-05 | 西安电子科技大学 | 一种基于bp-pmbm滤波算法的多目标跟踪方法 |
CN111457916B (zh) * | 2020-03-30 | 2021-05-07 | 中国人民解放军国防科技大学 | 基于拓展标记随机有限集的空间碎片目标跟踪方法和装置 |
CN111504326B (zh) * | 2020-04-30 | 2023-10-27 | 江苏理工学院 | 一种基于t分布的鲁棒glmb多目标跟踪方法 |
CN111812637B (zh) * | 2020-06-02 | 2022-12-02 | 杭州电子科技大学 | 具有类型概率的l-rfs混合目标结构建模与估计方法 |
CN112113572B (zh) * | 2020-09-18 | 2023-04-07 | 桂林电子科技大学 | 一种解决分布式标签融合的多目标跟踪方法 |
CN112305915B (zh) * | 2020-10-28 | 2023-02-07 | 深圳大学 | Tsk迭代回归模型的标签多伯努利多目标跟踪方法及*** |
CN114493885A (zh) * | 2022-03-30 | 2022-05-13 | 支付宝(杭州)信息技术有限公司 | 策略组合的优化方法及装置 |
CN115097437B (zh) * | 2022-06-06 | 2023-06-09 | 哈尔滨工程大学 | 一种基于标签多伯努利检测前跟踪算法的水下目标跟踪轨迹临近交叉解决方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005165791A (ja) * | 2003-12-03 | 2005-06-23 | Fuji Xerox Co Ltd | 対象物の追跡方法及び追跡システム |
CN106199581A (zh) * | 2016-06-30 | 2016-12-07 | 电子科技大学 | 一种随机集理论下的多机动目标跟踪方法 |
CN106408594A (zh) * | 2016-09-28 | 2017-02-15 | 江南大学 | 基于多伯努利特征协方差的视频多目标跟踪方法 |
CN106707272A (zh) * | 2016-07-01 | 2017-05-24 | 电子科技大学 | 一种基于随机集理论的多目标跟踪方法 |
CN107765179A (zh) * | 2017-06-26 | 2018-03-06 | 河海大学 | 一种适用于量测丢失的发电机动态状态估计方法 |
CN111007880A (zh) * | 2019-12-24 | 2020-04-14 | 桂林电子科技大学 | 一种基于汽车雷达的扩展目标跟踪方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10310068B2 (en) * | 2014-12-08 | 2019-06-04 | Northrop Grumman Systems Corporation | Variational track management |
CN104766320B (zh) * | 2015-04-02 | 2017-06-13 | 西安电子科技大学 | 阈值化量测下的多伯努利滤波弱目标检测与跟踪方法 |
CN105354860B (zh) * | 2015-08-26 | 2018-03-06 | 西安电子科技大学 | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 |
CN106910205A (zh) * | 2017-03-03 | 2017-06-30 | 深圳市唯特视科技有限公司 | 一种基于随机有限集滤波器耦合的多目标跟踪方法 |
CN107102295A (zh) * | 2017-04-13 | 2017-08-29 | 杭州电子科技大学 | 基于glmb滤波的多传感器tdoa无源定位方法 |
CN107145612B (zh) * | 2017-06-23 | 2020-11-10 | 北京邮电大学 | 一种主题相关的影响力用户发现和追踪方法 |
CN107677997B (zh) * | 2017-09-28 | 2021-06-29 | 杭州电子科技大学 | 基于GLMB滤波和Gibbs采样的扩展目标跟踪方法 |
CN107730537B (zh) * | 2017-09-29 | 2020-07-07 | 桂林电子科技大学 | 基于箱粒子概率假设密度滤波的弱目标检测与跟踪方法 |
CN107703496B (zh) * | 2017-10-12 | 2021-04-30 | 桂林电子科技大学 | 一种交互式多模伯努利滤波的机动弱目标检测前跟踪方法 |
-
2018
- 2018-12-18 CN CN201811549885.5A patent/CN109508444B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005165791A (ja) * | 2003-12-03 | 2005-06-23 | Fuji Xerox Co Ltd | 対象物の追跡方法及び追跡システム |
CN106199581A (zh) * | 2016-06-30 | 2016-12-07 | 电子科技大学 | 一种随机集理论下的多机动目标跟踪方法 |
CN106707272A (zh) * | 2016-07-01 | 2017-05-24 | 电子科技大学 | 一种基于随机集理论的多目标跟踪方法 |
CN106408594A (zh) * | 2016-09-28 | 2017-02-15 | 江南大学 | 基于多伯努利特征协方差的视频多目标跟踪方法 |
CN107765179A (zh) * | 2017-06-26 | 2018-03-06 | 河海大学 | 一种适用于量测丢失的发电机动态状态估计方法 |
CN111007880A (zh) * | 2019-12-24 | 2020-04-14 | 桂林电子科技大学 | 一种基于汽车雷达的扩展目标跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109508444A (zh) | 2019-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508444B (zh) | 区间量测下交互式多模广义标签多伯努利的快速跟踪方法 | |
CN101975575B (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
CN107402381B (zh) | 一种迭代自适应的多机动目标跟踪方法 | |
CN111693984B (zh) | 一种改进的ekf-ukf动目标跟踪方法 | |
CN111562571B (zh) | 一种未知新生强度的机动多目标跟踪与航迹维持方法 | |
CN107703496A (zh) | 一种交互式多模伯努利滤波的机动弱目标检测前跟踪方法 | |
CN110780290B (zh) | 基于lstm网络的多机动目标跟踪方法 | |
CN109934849B (zh) | 基于轨迹度量学习的在线多目标跟踪方法 | |
CN106054167B (zh) | 基于强度滤波器的多扩展目标跟踪方法 | |
CN106707272B (zh) | 一种基于随机集理论的多目标跟踪方法 | |
CN108871365B (zh) | 一种航向约束下的状态估计方法及*** | |
CN104777465B (zh) | 基于b样条函数任意扩展目标形状及状态估计方法 | |
CN111325776A (zh) | 一种基于变分贝叶斯t分布卡尔曼滤波的phd多目标跟踪方法 | |
CN108320302B (zh) | 基于随机超曲面的CBMeMBer多目标跟踪方法 | |
CN113409363B (zh) | 一种基于bp-pmbm滤波算法的多目标跟踪方法 | |
CN116047495B (zh) | 一种用于三坐标雷达的状态变换融合滤波跟踪方法 | |
CN104849697B (zh) | 一种基于去偏坐标转换的α‑β滤波方法 | |
CN107562837B (zh) | 一种基于道路网的机动目标跟踪算法 | |
CN111262556A (zh) | 一种同时估计未知高斯测量噪声统计量的多目标跟踪方法 | |
CN104931932A (zh) | 一种改进的去偏坐标转换卡尔曼滤波方法 | |
CN111811515B (zh) | 一种基于高斯混合概率假设密度滤波器的多目标航迹提取方法 | |
Lan et al. | Variational nonlinear Kalman filtering with unknown process noise covariance | |
CN115542309A (zh) | 基于imm-vb的非高斯噪声条件下机动目标自适应跟踪方法 | |
CN114910075A (zh) | 基于量测互斥分组的快速标签多伯努利多目标跟踪方法 | |
Nagappa et al. | Fast sequential monte carlo phd smoothing |
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 |