CN104794735B - 基于变分贝叶斯期望最大化的扩展目标跟踪方法 - Google Patents
基于变分贝叶斯期望最大化的扩展目标跟踪方法 Download PDFInfo
- Publication number
- CN104794735B CN104794735B CN201510152626.9A CN201510152626A CN104794735B CN 104794735 B CN104794735 B CN 104794735B CN 201510152626 A CN201510152626 A CN 201510152626A CN 104794735 B CN104794735 B CN 104794735B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- msup
- component
- 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 44
- 238000005259 measurement Methods 0.000 claims abstract description 55
- 238000009826 distribution Methods 0.000 claims description 46
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000013138 pruning Methods 0.000 claims description 11
- 238000005192 partition Methods 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 abstract description 6
- 230000007423 decrease Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 11
- 238000004088 simulation Methods 0.000 description 10
- 230000004083 survival effect Effects 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 101001086191 Borrelia burgdorferi Outer surface protein A Proteins 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 238000004235 valence bond calculation Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种基于变分贝叶斯期望最大化VBEM的扩展目标跟踪方法,主要解决传统的扩展目标跟踪领域中,在量测噪声协方差是未知的情况下,目标的跟踪性能将会急剧下降的问题。该方法首先预测了目标状态和量测噪声协方差的联合概率假设密度中高斯逆伽玛分量的相关参数;再对高斯逆伽玛分量参数进行更新;最后通过修建与合并得到扩展目标状态和数目。仿真实验表明,本发明能很好地跟踪未知数目和未知量测噪声协方差下的多扩展目标,且具有较高的跟踪精度,可用于对飞行器和潜艇目标的跟踪。
Description
技术领域
本发明属于信息处理技术领域,特别涉及一种目标跟踪方法,可用于跟踪多扩展目标。
背景技术
在传统的目标跟踪领域里,由于雷达的分辨率有限,因此通常情况下将目标看作是点目标,即一个时刻每个目标只能产生一个量测。近年来,随着雷达探测技术的发展和实际应用的需要,更多的将目标视为扩展目标,即每个目标在每一个时刻都可以产生多个量测。
实际的目标跟踪场景中,目标的数目是无法提前预知的,因此随机集理论的提出极大的满足了目标跟踪理论的需要。在诸多对目标的模型假设中,尤其属扩展目标理论的提出更加贴近目前跟踪理论的需要,并且在实际生活中得到了广泛的应用,近十年来,成为了目标跟踪领域中的研究热点。2003年,Mahler将随机有限集理论应用于多目标跟踪问题,提出了概率假设密度PHD滤波。2005年,Gilholm和Salmond提出一种空间分布服从泊松分布的扩展目标模型。2009年,Mahler推导出了扩展目标PHD滤波,即用每一时刻的量测随机集对目标随机集进行预测、更新,继而可以准确提取目标的运动状态和估计目标的数目。2010年,等给出了扩展目标PHD的高斯混合实现形式。2011年,Orguner等又提出了带势分布的扩展目标PHD(ET-CPHD)滤波,很好的解决了ET-PHD估计目标数目时的缺陷。然而,传统的扩展目标跟踪算法处理的都是量测噪声协方差是已知的情形,实际中当量测噪声协方差未知时,扩展目标的跟踪性能将会急剧下降。
发明内容
本发明的目的在于针对上述问题,提出一种基于变分贝叶斯期望最大化的扩展目标跟踪方法,以提高在量测噪声协方差未知条件下的跟踪性能。
实现本发明的技术关键是:在势概率假设密度滤波框架下,引入变分贝叶斯期望最大化VBEM技术,估计目标状态和量测噪声协方差的联合概率假设密度,实现未知量测噪声协方差下的目标跟踪问题。
VB是一类用于贝叶斯估计和机器学习领域中近似计算复杂积分的方法,本文中VB被用来近似线性高斯***中,扩展目标状态和量测噪声协方差的联合概率假设密度,该算法的主要思想是对扩展目标状态和量测噪声协方差的联合概率假设密度进行参数化近似,并给出其参数化表达形式。在用VB理论对扩展目标状态和量测噪声协方差的联合概率假设密度进行近似的过程中,为了判断所估计的高斯逆伽玛分量相关参数的性能,在VB的基础上又引入了期望最大化EM算法。在期望E步,估计未知参数的期望值,给出当前的参数估计;在最大化M步,重新估计分布参数,以使该似然函数达到最大。M步上得到的参数估计值被用于下一个E步计算,这个过程不断交替进行。这种方法可以广泛的应用于数据有缺损的情况,且具有简单性和稳定性等优点。
本发明利用上述VBEM技术进行扩展目标跟踪的技术步骤包括如下:
(1)在时刻k=0时,初始化扩展目标状态和量测噪声协方差的联合概率假设密度为v0(x,R);
(2)在k≥1时,对扩展目标状态和量测噪声协方差的联合概率假设密度vk-1|k-1(x,R)及用于计算扩展目标数目的势分布Pk-1|k-1(num)进行预测,得到预测的扩展目标联合概率假设密度vk|k-1(x,R)和预测势分布Pk|k-1(num);
(3)对预测的扩展目标状态和量测噪声协方差的联合概率假设密度vk|k-1(x,R)及用于计算扩展目标数目的势分布Pk|k-1(num)进行更新:
(3a)利用变分贝叶斯VB方法对联合概率假设密度vk|k-1(x,R)进行近似,得到用高斯分布的求和形式表示的扩展目标状态的概率假设密度Qx,k|k-1(x)和用逆伽玛分布的求和形式表示的量测噪声协方差的概率假设密度QR,k|k-1(R);
(3b)利用变分贝叶斯期望最大化VBEM方法对扩展目标状态的概率假设密度Qx,k|k-1(x)中的高斯分量和量测噪声协方差的概率假设密度QR,k|k-1(R)中的逆伽玛分量进行迭代更新,得到表示扩展目标运动状态x的高斯分量和表示量测噪声协方差R的逆伽玛分量;
(3c)对步骤(2)预测得到的势分布Pk|k-1(num)进行更新,得到更新后的势分布Pk|k(num);
(4)对更新后的高斯分量和逆伽玛分量进行修剪与合并,并提取合并后的高斯分量和逆伽玛分量的位置和速度作为扩展目标的状态;
(5)对步骤(3)更新得到的势分布Pk|k(num)进行加权平均,得到扩展目标的数目:
(6)重复步骤(2)-(5),继续跟踪扩展目标。
本发明具有以下优点:
第一,本发明引入了变分贝叶斯EM技术,通过估计扩展目标状态和量测噪声协方差的联合概率假设密度,有效估计了各个目标在不同时刻的真实量测噪声,为复杂环境下多扩展目标跟踪场景的分析提供了帮助,保证了CPHD滤波算法可以有效地实现对未知量测噪声协方差环境下扩展目标跟踪。
第二,本发明由于采用先预测高斯逆伽玛分量的相关参数,再对高斯逆伽玛分量参数进行更新,最后通过修剪与合并得到扩展目标状态和数目的跟踪扩展目标的过程,与传统的GM-CPHD滤波算法相比,提高了跟踪精度。
附图说明
图1是本发明的总流程图;
图2是单次实验条件下,用本发明方法跟踪扩展目标的仿真结果图;
图3是在100次实验条件下,用本发明方法与传统GM-CPHD方法估计目标数目的仿真结果对比图;
图4是在100次实验条件下,用本发明方法与传统GM-CPHD方法借助OSPA距离来判断目标跟踪精度的仿真结果对比图。
具体实施方式
下面结合附图对本发明做进一步的描述。
参照图1,本发明的具体实现步骤包括如下:
步骤1,在时刻k=0时,初始化目标状态和量测噪声协方差的联合概率假设密度:
其中,J0表示高斯分量的个数,表示第i个高斯分量的权值,N(·)表示高斯分布,表示第i个高斯分量的均值,表示第i个高斯分量的协方差;IG(·)表示逆伽玛分布,表示第i个逆伽玛分量的协方差,为第i个逆伽玛分量的常量因子,为第i个逆伽玛分量的迭代因子,l表示量测噪声协方差的第l维,d表示量测噪声协方差的维数。
步骤2,在k≥1时,对扩展目标状态和量测噪声协方差的联合概率假设密度vk-1|k-1(x,R)进行预测,得到预测的扩展目标联合概率假设密度vk|k-1(x,R)。
2a)对存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中高斯分量的均值和协方差进行预测,得到预测的存活目标高斯分量的均值和协方差
其中,表示存活目标状态转移矩阵,Qk-1表示存活目标过程噪声协方差,(·)T表示转置;
2b)对存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中逆伽玛分量的常量因子和迭代因子进行预测,得到预测的存活目标逆伽玛分量的常量因子和逆伽玛分量的迭代因子
其中ρl表示遗忘因子,且ρl∈(0,1];
2c)在高斯混合框架下,利用预测的存活目标高斯分量的均值协方差和逆伽玛分量的常量因子迭代因子计算存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R):
其中,PS,k表示扩展目标存活概率,Jk-1表示k-1时刻的高斯分量数,表示k-1时刻第i个高斯分量的权值,N(·)表示高斯分布,表示预测的存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中第i个高斯分量的均值,表示预测的存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中第i个高斯分量的协方差;IG(·)表示逆伽玛分布,d表示量测噪声协方差的维数,表示预测的存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中第i个逆伽玛分量的常量因子,表示预测的存活目标状态和量测噪声协方差的联合概率假设密度vS,k|k-1(x,R)中第i个逆伽玛分量的迭代因子;
2d)对衍生目标状态和量测噪声协方差的联合概率假设密度bk|k-1(x,R)中高斯分量的均值和协方差进行预测,得到预测的衍生目标高斯分量的均值和协方差
其中,i表示k-1时刻高斯逆伽玛分量的第i个,j表示由k-1时刻到k时刻衍生出的高斯逆伽玛分量的第j个,表示衍生目标的状态转移矩阵,表示衍生目标的状态修正量,表示衍生目标过程噪声协方差;
2e)对衍生目标状态和量测噪声协方差的联合概率假设密度bk|k-1(x,R)中逆伽玛分量的常量因子和迭代因子进行预测,得到预测的衍生目标逆伽玛分量的常量因子和迭代因子
2f)在高斯混合框架下,利用预测的衍生目标高斯分量的均值协方差和逆伽玛分量的常量因子迭代因子计算衍生目标状态和量测噪声协方差的联合概率假设密度bk|k-1(x,R):
其中,Jb,k表示k-1时刻到k时刻的衍生目标高斯分量数,表示k时刻第j个衍生目标高斯分量的权值,表示k-1时刻到k时刻由第i个高斯分量衍生得到第j个高斯分量的均值,表示k-1时刻到k时刻由第i个高斯分量衍生得到第j个高斯分量的协方差;表示k-1时刻到k时刻由第i个逆伽玛分量衍生得到第j个逆伽玛分量的常量因子,表示k-1时刻到k时刻由第i个逆伽玛分量衍生得到第j个逆伽玛分量的迭代因子;
2g)计算新生目标状态和量测噪声协方差的联合概率密度γk(x,R):
其中,Jγ,k为k时刻新生目标高斯分量数,为为k时刻新生目标第i个高斯分量的权值,为新生目标第i个高斯分量的状态均值,为新生目标第i个高斯分量的运动状态协方差;为新生目标第i个逆伽玛分量的常量因子,为新生目标第i个逆伽玛分量的迭代因子。
2h)利用步骤2a)到步骤2g)得到的参数,计算扩展目标状态和量测噪声协方差的联合概率假设密度vk|k-1(x,R):
vk|k-1(x,R)=vS,k|k-1(x,R)+bk|k-1(x,R)+γk(x,R)。
步骤3,在k≥1时,对势分布Pk-1|k-1(num)进行预测,得到预测势分布Pk|k-1(num):
其中,Pbirth,k(num-j)表示k-1时刻到k时刻新生num-j个扩展目标的概率,pS,k表示k时刻扩展目标的存活概率,Pk-1|k-1(h)表示k-1时刻有h个扩展目标的概率,!表示阶乘,表示k-1时刻到k时刻有j个扩展目标存活的概率,(1-pS,k)h-j表示k-1时刻到k时刻有h-j个扩展目标消亡的概率。
步骤4.对扩展目标状态和量测噪声协方差的联合概率假设密度vk|k-1(x,R)及预测势分布Pk|k-1(num)进行更新:
4a)利用变分贝叶斯VB方法将联合概率假设密度vk|k-1(x,R)近似为:
式中,Qx,k|k-1(x)为高斯分布的求和形式,
QR,k|k-1(R)为逆伽玛分布的求和形式,
其中,表示第i个高斯分量在第k个时刻的预测权值,i=1,...,Jk,Jk表示第k个时刻扩展目标高斯分量的个数,N(·)表示高斯分布,为第k个时刻预测得到的第i个高斯分量的均值,为第k个时刻预测得到的第i个高斯分量的协方差;IG(·)表示逆伽玛分布,为第k个时刻预测得到的第i个逆伽玛分量的常量因子,为第k个时刻预测得到的第i个逆伽玛分量的迭代因子,l=1,...,d,d表示量测噪声协方差的维数;
4b)利用变分贝叶斯期望最大化VBEM方法对扩展目标状态的概率假设密度Qx,k|k-1(x)中的高斯分量和量测噪声协方差的概率假设密度QR,k|k-1(R)中的逆伽玛分量进行迭代更新:
(4b1)设定逆伽玛分量的常量因子和迭代因子 其中l=1,...,d,d为量测噪声协方差R的维数;
(4b2)根据设定的逆伽玛分量的两个因子,计算量测噪声协方差:
其中n=1,...,N,N为最大迭代次数,diag[...]表示对角化其中的元素;
(4b3)利用量测噪声协方差计算更新因子
其中,表示对量测噪声协方差矩阵进行对角线连接当前单元W的量测数目后的矩阵,blkdiag(·)表示对其中的元素进行对角线连接,|W|表示当前单元W的量测数目;HW表示对观测矩阵Hk垂直连接当前单元W的量测数目后的矩阵,表示k时刻的观测矩阵Hk的转置;表示k-1到k时刻预测的扩展目标高斯分量运动状态协方差,表示矩阵HW的转置;
(4b4)利用更新因子计算增益矩阵
其中[·]-1表示对矩阵求逆;
(4b5)利用增益矩阵计算扩展目标高斯分量运动状态和扩展目标高斯分量运动状态协方差
其中,I表示一个单位矩阵,zW表示某个划分单元W中的所有量测;
(4b6)提取扩展目标高斯分量运动状态的位置信息,将该位置信息用表示;
(4b7)利用该高斯分量位置信息和量测噪声协方差计算量测Yn′由位置处的高斯分量生成的概率γn′i:
其中,Jk表示扩展目标高斯分量数目,Yn′表示当前单元W的第n′个量测,n′=1,...,|W|;N(·)表示高斯分布;πi代表混合系数,Ni表示由位置处的高斯分量产生的有效量测数目,
(4b8)利用量测Yn′由位置处的高斯分量生成的概率γn′i,迭代更新得到扩展目标高斯分量运动状态的位置信息
(4b9)利用扩展目标高斯分量运动状态的位置信息混合系数πi,量测噪声协方差计算最大似然函数L(i)(n):
其中Jk表示扩展目标高斯分量数目;
(4b10)判断|L(i)(n)-L(i)(n-1)|是否小于常量ε=0.01,同时判断当前迭代次数n是否小于最大迭代次数N=100,如果是,则停止迭代;否则,返回步骤(4b2),更新逆伽玛分量迭代因子:
其中,表示对向量中所有元素进行相加,
表示对向量的第j维元素的平方,(·)jj表示取矩阵的对角线元素,表示对迭代因子垂直连接当前单元W的量测数目后的向量,zW表示当前单元的量测;
(4b11)提取扩展目标状态分量扩展目标运动状态协方差迭代因子即,其中扩展目标状态分量中的位置信息用的是步骤(4b8)中迭代更新得到的扩展目标分量运动状态的位置信息
4c)对势分布Pk|k-1(num)进行更新,得到更新后的势分布Pk|k(num)表示如下:
其中,p∠Z表示把量测集Z划分成p个非空子集,W∈p表示第p个非空子集下的某一个单元,Gk|k-1(ρ)表示状态预测概率生成函数,表示状态预测概率生成函数的num阶偏导,GFA(0)表示没有量测时的虚警概率生成函数,ηW表示扩展目标产生量测概率,|p|表示第p个划分下的所有非空单元数目,表示虚警概率生成函数的|W|阶偏导,δnum≥|p|表示当目标数目num大于划分单元|p|时取值为1,否则为0,|Z|=0表示扩展目标没有产生量测,|W|表示各非空单元W中的量测个数,lp,W表示量测划分单元为|p|-1时的虚警常系数,ψp,W表示目标产生量测概率的乘积,ρ表示扩展目标分量没有被检测到的概率:
其中,Π为连乘符号,表示第j个高斯逆伽玛分量在当前所有高斯逆伽玛分量中所占的比重,pk|k-1表示单扩展目标状态转移概率密度函数,pz(z′|·)表示扩展目标量测似然,pFA(z′)表示虚警量测似然,Gz(0|·)表示量测概率生成函数,表示量测概率生成函数的|W|阶偏导,z′∈W表示量测z′属于单元W,PD(·)表示检测概率,W′∈p-W表示p划分下所有单元中除去单元W后剩下的单元,ηW′表示扩展目标虚警量测产生概率,表示状态预测概率生成函数的|p|阶偏导,表示状态预测概率生成函数的|p|-1阶偏导。
步骤5.对更新后的高斯分量和逆伽玛分量进行修剪与合并,其步骤如下:
(5a)设置两个修剪门限T1和T2,一个合并门限U:T1=10-5,T2=120,U=10;设置最大高斯逆伽玛分量数目:Jmax=100;
(5b)计算每个扩展目标分量对应的量测噪声协方差:
(5c)设变量l′=0,对更新后的扩展目标分量进行修剪,得到修剪后的扩展目标分量对应的序号集合I为:
(5d)令l′=l′+1,取表示取最大权值所对应的集合I中元素i′,对修剪后的扩展目标分量中满足合并门限U的分量进行提取,得到适合合并的扩展目标分量对应的序号集合为:
(5e)分别对序号集合中对应的扩展目标分量的权值运动状态常量因子迭代因子协方差进行合并,得到合并后的扩展目标分量的权值运动状态常量因子迭代因子协方差如下:
(5f)将步骤(5c)得到的修剪后的扩展目标分量对应的序号集合I中与步骤(5d)得到的适合合并的扩展目标分量对应的序号集合中相同的元素去除掉,然后判断修剪后的扩展目标分量对应的序号集合I是否为空集,如果不为空集,则返回步骤(5d),否则执行(5g);
(5g)判断变量l′是否大于最大高斯逆伽玛分量数目Jmax,如果l′>Jmax,则将权值对应的高斯逆伽玛分量按从大到小排列,并取前Jmax个权值大于0.5的高斯逆伽玛分量的位置和速度作为扩展目标的状态;如果l′<Jmax,则将所有权值大于0.5对应的高斯逆伽玛分量的位置和速度作为扩展目标的状态。
步骤6.根据步骤4c)更新得到的势分布Pk|k(num)进行加权平均,得到扩展目标的数目:
得到扩展目标的状态和数目以后即可完成对扩展目标的跟踪。
本发明用于跟踪扩展目标的效果可通过以下仿真实验进一步说明:
1.仿真条件
考虑二维平面中做匀速直线运动4个目标,并且在两个位置出现运动轨迹交叉的情况,采样周期为整个观测过程持续40个时刻,目标的运动方程和量测方程分别为:
xk=Fxk-1+Gwk
yk=Hxk+vk
其中,高斯噪声wk的标准差σw=2。扩展目标状态转移矩阵F,输入矩阵G,高斯噪声wk,观测矩阵H设置为:
实际目标跟踪场景中量测噪声vk的协方差是未知的,本次实验中设真实的量测噪声标准差σv=1;
设传统GM-CPHD算法中采用的三个量测噪声标准差分别为σ=0.05,1,50;
设目标产生量测数目服从泊松分布,泊松分布参数β=20,目标产生量测位置服从高斯分布;
设目标存活概率pS,k=0.99,检测概率pD,k=0.98;
设杂波数服从均值为5的泊松分布,杂波在整个观测区域内服从均匀分布。
设置新生目标的状态为新生目标的协方差均为Pγ=diag[10,10,10,10];
设置修剪门限T1=10-5,修剪门限T2=120,合并门限U=10,最大高斯逆伽玛分量数目Jmax=100,遗忘因子ρ=0.9,实验仿真次数为100次。
2.仿真内容与结果
仿真1,采用本发明方法对未知量测噪声协方差下扩展目标运动轨迹的单次实验进行跟踪,结果如图2所示。从图2可以看出,本发明方法可以很好的跟踪扩展目标运动轨迹。
仿真2,采用本发明方法与传统的分别采用不同量测噪声协方差的三个GM-CPHD方法进行100次实验的目标数目估计对比,结果如图3所示。从图3可以看出,本发明方法与传统的给定真实量测噪声协方差下的GM-CPHD方法对目标数目的估计效果一样良好。但是,传统的GM-CPHD方法,当量测噪声协方差未知时,且使用的量测噪声协方差与真实值偏差比较大时,估计的目标数目会发生很大的偏差,同时在k=20,k=35这两个目标运动轨迹出现交叉时刻,目标数目出现漏检;
仿真3,采用本发明方法与传统的分别采用不同量测噪声协方差的三个GM-CPHD方法借助OSPA距离进行100次实验的跟踪精度对比,结果如图4所示。从图4可以看出,本发明方法与传统的给定真实量测噪声协方差下的GM-CPHD算法的跟踪精度一样良好。但是,传统的GM-CPHD方法,当使用的量测噪声协方差与真实的量测噪声协方差偏差比较大时,跟踪精度会变得很差。
实验表明,本发明方法在处理量测噪声协方差是未知条件下的扩展目标跟踪场景时,其跟踪效果优于传统的GM-CPHD扩展目标跟踪方法。
Claims (4)
1.一种基于变分贝叶斯期望最大化的扩展目标跟踪方法,包括如下步骤:
(1)在时刻k=0时,初始化扩展目标状态和量测噪声协方差的联合概率假设密度为v0(x,R);
(2)在k≥1时,对扩展目标状态和量测噪声协方差的联合概率假设密度vk-1|k-1(x,R)及用于计算扩展目标数目的势分布Pk-1|k-1(num)进行预测,得到预测的扩展目标联合概率假设密度vk|k-1(x,R)和预测势分布Pk|k-1(num);
(3)对预测的扩展目标状态和量测噪声协方差的联合概率假设密度vk|k-1(x,R)及用于计算扩展目标数目的势分布Pk|k-1(num)进行更新:
(3a)利用变分贝叶斯VB方法对联合概率假设密度vk|k-1(x,R)进行近似,得到用高斯分布的求和形式表示的扩展目标状态的概率假设密度Qx,k|k-1(x)和用逆伽玛分布的求和形式表示的量测噪声协方差的概率假设密度QR,k|k-1(R);
(3b)利用变分贝叶斯期望最大化VBEM方法对扩展目标状态的概率假设密度Qx,k|k-1(x)中的高斯分量和量测噪声协方差的概率假设密度QR,k|k-1(R)中的逆伽玛分量进行迭代更新,得到表示扩展目标运动状态x的高斯分量和表示量测噪声协方差R的逆伽玛分量:
(3b1)设定逆伽玛分量的常量因子和迭代因子 其中l=1,...,d,d为量测噪声协方差R的维数;
(3b2)根据设定的逆伽玛分量两个因子,计算量测噪声协方差:其中n=1,...,N,N为最大迭代次数,diag[...]表示对角化其中的元素;
(3b3)利用量测噪声协方差计算更新因子
<mrow>
<msubsup>
<mi>S</mi>
<mi>W</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>W</mi>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msup>
<msub>
<mi>H</mi>
<mi>W</mi>
</msub>
<mi>T</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mi>W</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
其中,表示对量测噪声协方差矩阵进行对角线连接当前单元W量测数目后的矩阵,blkdiag(·)表示对其中的元素进行对角线连接,|W|表示当前单元W的量测数目;HW表示对观测矩阵Hk垂直连接当前单元W的量测数目后的矩阵, 表示k时刻的观测矩阵Hk的转置;表示k-1到k时刻预测的扩展目标高斯分量运动状态协方差,表示矩阵HW的转置;
(3b4)利用更新因子计算增益矩阵
<mrow>
<msubsup>
<mi>K</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>H</mi>
<mi>W</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>S</mi>
<mi>W</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
其中[·]-1表示对矩阵求逆;
(3b5)利用增益矩阵计算扩展目标高斯分量运动状态和扩展目标高斯分量运动状态协方差
<mrow>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>K</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>W</mi>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mi>W</mi>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>K</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>H</mi>
<mi>W</mi>
</msub>
<mo>&rsqb;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
其中,I表示一个单位矩阵,zW表示某个划分单元W中的所有量测;
(3b6)提取扩展目标高斯分量运动状态的位置信息,将该位置信息用表示;
(3b7)利用该高斯分量位置信息和量测噪声协方差计算量测Yn′是由位置处的高斯分量生成的概率γn′i:
<mrow>
<msub>
<mi>&gamma;</mi>
<mrow>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&pi;</mi>
<mi>i</mi>
</msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>|</mo>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>R</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</munderover>
<msub>
<mi>&pi;</mi>
<mi>i</mi>
</msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>|</mo>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>R</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,Jk表示扩展目标高斯分量数目,Yn′表示当前单元W的第n′个量测,n′=1,...,|W|;N(·)表示高斯分布;πi代表混合系数,Ni表示由位置处的高斯分量产生的有效量测数目,
(3b8)利用量测Yn′由位置处的高斯分量生成的概率γn′i,迭代更新得到扩展目标高斯分量运动状态的位置信息
<mrow>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
</mrow>
</munderover>
<msub>
<mi>&gamma;</mi>
<mrow>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>Y</mi>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>;</mo>
</mrow>
(3b9)利用扩展目标高斯分量运动状态的位置信息混合系数πi,量测噪声协方差计算最大似然函数L(i)(n):
其中Jk表示扩展目标高斯分量数目;
(3b10)判断当前迭代次数n是否小于最大迭代次数N=100,当前迭代次数n不小于最大迭代次数N=100,则停止迭代,当前迭代次数n小于最大迭代次数N=100,且|L(i)(n)-L(i)(n-1)|小于常量ε=0.01,则停止迭代,否则返回步骤(3b2),更新逆伽玛分量迭代因子:
其中,表示对向量中所有元素进行相加, 表示对向量的第j维元素的平方,(·)jj表示取矩阵的对角线元素,表示对迭代因子垂直连接当前单元W量测数目后的向量,
(3b11)提取扩展目标状态分量扩展目标运动状态协方差迭代因子即,其中扩展目标状态分量中的位置信息用的是步骤(3b8)中迭代更新得到的扩展目标分量运动状态的位置信息
(3c)对步骤(2)预测得到的势分布Pk|k-1(num)进行更新,得到更新后的势分布Pk|k(num);
(4)对更新后的高斯分量和逆伽玛分量进行修剪与合并,并提取合并后的高斯分量和逆伽玛分量的位置和速度作为扩展目标的状态;
(5)对步骤(3)更新得到的势分布Pk|k(num)进行加权平均,得到扩展目标的数目:
(6)重复步骤(2)-(5),继续跟踪扩展目标。
2.根据权利要求1所述的基于变分贝叶斯期望最大化的扩展目标跟踪方法,其中,步骤(3a)所述的利用变分贝叶斯VB方法对联合概率假设密度vk|k-1(x,R)进行近似,按如下公式进行:
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>;</mo>
<msubsup>
<mi>m</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>d</mi>
</munderover>
<mi>I</mi>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&beta;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
其中,Qx,k|k-1(x)为高斯分布的求和形式,表示为
QR,k|k-1(R)为逆伽玛分布的求和形式,表示为
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>J</mi>
<mi>k</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>d</mi>
</munderover>
<mi>I</mi>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&beta;</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
表示第i个高斯分量在第k个时刻的权值,i=1,...,Jk,Jk表示第k个时刻扩展目标高斯分量的个数,N(·)表示高斯分布,为第k个时刻预测得到的第i个高斯分量的均值,为第k个时刻预测得到的第i个高斯分量的协方差;IG(·)表示逆伽玛分布,为第k个时刻预测得到的第i个逆伽玛分量的常量因子,为第k个时刻预测得到的第i个逆伽玛分量的迭代因子,l=1,...,d,d表示量测噪声协方差的维数。
3.根据权利要求1所述的基于变分贝叶斯期望最大化的扩展目标跟踪方法,其中,所述步骤(3c)对势分布Pk|k-1(num)进行更新,得到更新后的势分布Pk|k(num)表示如下:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>&angle;</mo>
<mi>Z</mi>
</mrow>
</msub>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>W</mi>
<mo>&Element;</mo>
<mi>p</mi>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</msub>
<msubsup>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>F</mi>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>W</mi>
</mrow>
<mrow>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
</mrow>
</mfrac>
<mfrac>
<msup>
<mi>&rho;</mi>
<mrow>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>-</mo>
<mrow>
<mo>|</mo>
<mi>P</mi>
<mo>|</mo>
</mrow>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>-</mo>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
<mo>)</mo>
<mo>!</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>&GreaterEqual;</mo>
<mrow>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>F</mi>
<mi>A</mi>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mfrac>
<msup>
<mi>&rho;</mi>
<mrow>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>-</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>-</mo>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>!</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>&GreaterEqual;</mo>
<mrow>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
</mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>&angle;</mo>
<mi>Z</mi>
</mrow>
</msub>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>W</mi>
<mo>&Element;</mo>
<mi>p</mi>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</msub>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mrow>
<mo>|</mo>
<mi>Z</mi>
<mo>|</mo>
</mrow>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msup>
<mi>&rho;</mi>
<mrow>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
</msup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>u</mi>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mrow>
<mo>|</mo>
<mi>Z</mi>
<mo>|</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,p∠Z表示把量测集Z划分成p个非空子集,W∈p表示第p个非空子集下的某一个单元,Gk|k-1(ρ)表示状态预测概率生成函数,表示状态预测概率生成函数的num阶偏导,GFA(0)表示没有量测时的虚警概率生成函数,ηW表示扩展目标产生量测概率,|p|表示第p个划分下的所有非空单元数目,表示虚警概率生成函数的|W|阶偏导,δnum≥|p|表示当目标数目num大于划分单元|p|时取值为1,否则为0,|Z|=0表示扩展目标没有产生量测,|W|表示各非空单元W中的量测个数,lp,W表示量测划分单元为|p|-1时的虚警常系数,ψp,W表示目标产生量测概率的乘积,ρ表示扩展目标分量没有被检测到的概率:
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mo>&Sigma;</mo>
<msubsup>
<mover>
<mi>&omega;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>j</mi>
</msubsup>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>G</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>|</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>&eta;</mi>
<mi>W</mi>
<mo>=</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
<msubsup>
<mi>G</mi>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>|</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
<munder>
<mo>&Pi;</mo>
<mrow>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>&Element;</mo>
<mi>W</mi>
</mrow>
</munder>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>|</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>F</mi>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>F</mi>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msubsup>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>W</mi>
</mrow>
<mrow>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>W</mi>
<mi>F</mi>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msubsup>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>p</mi>
<mo>|</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Pi;</mo>
<mrow>
<msup>
<mi>W</mi>
<mo>&prime;</mo>
</msup>
<mo>&Element;</mo>
<mi>p</mi>
<mo>-</mo>
<mi>W</mi>
</mrow>
</munder>
<msup>
<mi>&eta;W</mi>
<mo>&prime;</mo>
</msup>
</mrow>
其中,Π为连乘符号,表示第j个高斯逆伽玛分量在当前所有高斯逆伽玛分量中所占的比重,pk|k-1表示单扩展目标状态转移概率密度函数,pz(z′|·)表示扩展目标量测似然,pFA(z′)表示虚警量测似然,Gz(0|·)表示量测概率生成函数,表示量测概率生成函数的|W|阶偏导,z′∈W表示量测z′属于单元W,PD(·)表示检测概率,W′∈p-W表示p划分下所有单元中除去单元W后剩下的单元,ηW′表示扩展目标虚警量测产生概率,表示状态预测概率生成函数的|p|阶偏导,表示状态预测概率生成函数的|p|-1阶偏导。
4.根据权利要求1所述的基于变分贝叶斯期望最大化的扩展目标跟踪方法,其中,步骤(4)所述的对更新后的高斯分量和逆伽玛分量进行修剪与合并,按如下步骤进行:
(4a)设置两个修剪门限T1和T2,一个合并门限U:T1=10-5,T2=120,U=10;设置最大高斯逆伽玛分量数目:Jmax=100;
(4b)计算每个扩展目标分量对应的量测噪声协方差:
(4c)设变量l′=0,对更新后的扩展目标分量进行修剪,得到修剪后的扩展目标分量对应的序号集合I为:
(4d)令l′=l′+1,取表示取最大权值所对应的集合I中元素i′,对修剪后的扩展目标分量中满足合并门限U的分量进行提取,得到适合合并的扩展目标分量对应的序号集合为:
(4e)分别对序号集合中对应的扩展目标分量的权值运动状态常量因子迭代因子协方差进行合并,得到合并后的扩展目标分量的权值运动状态常量因子迭代因子协方差如下:
<mrow>
<msubsup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>m</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>&alpha;</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>&alpha;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>&beta;</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>&beta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mover>
<mi>P</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>w</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>m</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>m</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
(4f)将步骤(4c)得到的修剪后的扩展目标分量对应的序号集合I中与步骤(4d)得到的适合合并的扩展目标分量对应的序号集合中相同的元素去除掉,然后判断修剪后的扩展目标分量对应的序号集合I是否为空集,如果不为空集则返回步骤(4d),否则执行(4g);
(4g)判断变量l′是否大于最大高斯逆伽玛分量数目Jmax,如果l′>Jmax,则将权值对应的高斯逆伽玛分量按从大到小排列,并取前Jmax个权值大于0.5的高斯逆伽玛分量的位置和速度作为扩展目标的状态;如果l′<Jmax,则将所有权值大于0.5对应的高斯逆伽玛分量的位置和速度作为扩展目标的状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510152626.9A CN104794735B (zh) | 2015-04-02 | 2015-04-02 | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510152626.9A CN104794735B (zh) | 2015-04-02 | 2015-04-02 | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104794735A CN104794735A (zh) | 2015-07-22 |
CN104794735B true CN104794735B (zh) | 2017-08-25 |
Family
ID=53559514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510152626.9A Expired - Fee Related CN104794735B (zh) | 2015-04-02 | 2015-04-02 | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104794735B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105652250B (zh) * | 2016-01-15 | 2018-01-05 | 西北工业大学 | 一种基于双层期望最大化的机动目标跟踪技术 |
WO2017124299A1 (zh) * | 2016-01-19 | 2017-07-27 | 深圳大学 | 基于序贯贝叶斯滤波的多目标跟踪方法及跟踪*** |
CN105913080B (zh) * | 2016-04-08 | 2019-02-22 | 西安电子科技大学昆山创新研究院 | 基于随机矩阵机动非椭圆扩展目标的联合跟踪与分类方法 |
CN108121846B (zh) * | 2016-11-29 | 2022-01-11 | 南京航空航天大学 | 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 |
CN107797106A (zh) * | 2017-05-08 | 2018-03-13 | 南京航空航天大学 | 一种加速em未知杂波估计的phd多目标跟踪平滑滤波方法 |
CN108519595A (zh) * | 2018-03-20 | 2018-09-11 | 上海交通大学 | 联合多传感器配准与多目标跟踪方法 |
CN108734725B (zh) * | 2018-04-11 | 2020-09-29 | 杭州电子科技大学 | 基于高斯过程的概率数据关联滤波扩展目标跟踪方法 |
CN109284677B (zh) * | 2018-08-16 | 2022-06-03 | 昆明理工大学 | 一种贝叶斯滤波目标跟踪算法 |
CN110909312B (zh) * | 2019-12-18 | 2022-04-22 | 哈尔滨工程大学 | 一种应用于rbmcda跟踪算法的目标消亡判断方法 |
CN112364292B (zh) * | 2020-09-24 | 2024-05-03 | 北京电子工程总体研究所 | 一种基于Ransac的密集目标跟踪方法、装置、设备和介质 |
CN116500575B (zh) * | 2023-05-11 | 2023-12-22 | 兰州理工大学 | 一种基于变分贝叶斯理论的扩展目标跟踪方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708550A (zh) * | 2012-05-17 | 2012-10-03 | 浙江大学 | 一种基于自然图像统计特性的盲去模糊算法 |
CN103235886A (zh) * | 2013-04-25 | 2013-08-07 | 杭州电子科技大学 | 一种基于变分贝叶斯容积强跟踪信息滤波的目标跟踪方法 |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8154600B2 (en) * | 2007-04-20 | 2012-04-10 | Utc Fire & Security Americas Corporation, Inc. | Method and system for distributed multiple target tracking |
-
2015
- 2015-04-02 CN CN201510152626.9A patent/CN104794735B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708550A (zh) * | 2012-05-17 | 2012-10-03 | 浙江大学 | 一种基于自然图像统计特性的盲去模糊算法 |
CN103235886A (zh) * | 2013-04-25 | 2013-08-07 | 杭州电子科技大学 | 一种基于变分贝叶斯容积强跟踪信息滤波的目标跟踪方法 |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
Non-Patent Citations (3)
Title |
---|
An Extended Target CPHD Filter and a Gamma Gaussian Inverse Wishart Implementation;Christian Lundquist et al;《IEEE Journal on Selected Topics in Signal Processing》;20130630;第7卷(第3期);第三部分 * |
Jinlong Yang et al.Adaptive probability hypothesis density filter based on variational Bayesian approximation for multi-target tracking.《IET Rada,Sonar & Navigation》.2013,第7卷(第9期), * |
张俊根 等.高斯混合粒子Cardinalized概率假设密度滤波被动测角多目标跟踪.《控制理论与应用》.2011,第28卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN104794735A (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104794735B (zh) | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 | |
CN104730511B (zh) | 星凸模型下的势概率假设密度多扩展目标跟踪方法 | |
CN101975575B (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
CN103729859B (zh) | 一种基于模糊聚类的概率最近邻域多目标跟踪方法 | |
CN103345577B (zh) | 变分贝叶斯概率假设密度多目标跟踪方法 | |
CN107703480B (zh) | 基于机器学习的混合核函数室内定位方法 | |
CN105354860B (zh) | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 | |
CN104155650A (zh) | 一种基于熵权值法点迹质量评估的目标跟踪方法 | |
CN111127523B (zh) | 基于量测迭代更新的多传感器gmphd自适应融合方法 | |
WO2018098926A1 (zh) | 一种适用于闪烁噪声的多目标跟踪方法及*** | |
CN109195110B (zh) | 基于层次聚类技术和在线极限学习机的室内定位方法 | |
CN108319570B (zh) | 一种异步多传感器空时偏差联合估计与补偿方法及装置 | |
CN103945531B (zh) | 基于信息熵的WLAN室内定位Radio Map更新方法 | |
CN114662535B (zh) | 一种基于变分贝叶斯学习的井下传感器网络目标跟踪方法 | |
Yao et al. | Improved support vector machine regression in multi-step-ahead prediction for rock displacement surrounding a tunnel | |
CN104268567A (zh) | 一种观测数据聚类划分的扩展目标跟踪方法 | |
CN104795063A (zh) | 一种基于声学空间非线性流形结构的声学模型构建方法 | |
CN103885867A (zh) | 一种模拟电路性能的在线评价方法 | |
Li et al. | Message passing based extended objects tracking with measurement rate and extension estimation | |
CN111340853A (zh) | 基于ospa迭代的多传感器gmphd自适应融合方法 | |
CN111262556A (zh) | 一种同时估计未知高斯测量噪声统计量的多目标跟踪方法 | |
Wan et al. | Variational Bayesian learning for robust AR modeling with the presence of sparse impulse noise | |
CN104850856A (zh) | 近邻传播观测聚类的多扩展目标跟踪方法 | |
CN113642785B (zh) | 基于先验信息的空间碎片轨道长期预报方法、***及设备 | |
Yang et al. | TDOA location based on modified Newton method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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: 20170825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |