CN103729637A - 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 - Google Patents
基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 Download PDFInfo
- Publication number
- CN103729637A CN103729637A CN201310753640.5A CN201310753640A CN103729637A CN 103729637 A CN103729637 A CN 103729637A CN 201310753640 A CN201310753640 A CN 201310753640A CN 103729637 A CN103729637 A CN 103729637A
- Authority
- CN
- China
- Prior art keywords
- formula
- target
- gauss
- represent
- item
- 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.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
技术领域
本发明属于目标跟踪方法技术领域,涉及一种扩展目标概率假设密度滤波方法,具体涉及一种基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法。
背景技术
目标跟踪就是结合传感器接收到的量测,对目标当前状态进行估计的过程。传统的多目标跟踪方法,如:联合数据关联、多假设跟踪等都是基于数据关联的,随着目标数目的增加,这些跟踪方法会存在“组合爆炸”、计算量呈指数级增长的问题;而近年来提出的基于随机有限集的概率假设密度(Probability Hypothesis Density,PHD)滤波方法,较为有效的避免了数据关联这个难题,能够直接对复杂环境中目标数未知且随时间变化的多目标的状态和个数进行联合估计,逐渐成为多目标跟踪领域的研究热点。但由于PHD滤波的递推公式中含有多个积分,所以它在一般情况下难以获得解析解。
随后,Mahler和Vo等学者分别给出了适用于非线性非高斯条件的序贯蒙特卡罗PHD(Sequential Monte Carlo PHD,SMC-PHD,又称为particle-PHD)滤波和适用于线性高斯条件的高斯混合PHD(Gaussian MixturePHD,GM-PHD)滤波“Ba-Ngu Vo,Wing-Kin Ma.The Gaussian mixtureprobability hypothesis density filter,IEEE Transactions on Signal Processing,2006,54(11):4091-4104.”。在此基础上,Daniel C和Clark D等别对SMC-PHD滤波方法和GM-PHD滤波方法的收敛性进行了证明。Mahler通过将PHD滤波进行扩展,得到了势概率假设密度滤波(cardinalized PHD,CPHD)滤波,该滤波方法可以获得更多的关于多目标状态集合的势分布的信息,能够得到比PHD滤波更加精确的估计,但其时间复杂度较大。类似于PHD滤波,在线性高斯系统下的近似实现方法称为GM-CPHD,在非线性非高斯系统下CPHD的近似实现方法称为SMC-CPHD。
在多目标跟踪系统中,一般都假设目标为点目标,即在每一时刻每个目标至多产生一个量测。但是,当传感器与目标之间的距离较小或传感器分辨率较高时,单个目标(例如一架飞机)的不同发射点(例如机头、机翼和机尾等不同部位)可同时产生多个观测,此时称该目标为一个扩展目标。在假设扩展目标数目服从泊松分布的条件下,Mahler于2009年在文献“GranstromK,Lundquist C,Orguner U.A Gaussian mixture PHD filter for extended targettracking.Proceedings of the13th International Conference on Infermation Fusion,2010,1-8.”中讨论了扩展目标的量测模型,并给出了扩展目标PHD(Extended-target PHD,EPHD)形式解法。紧接着,Granstrom等人给出了EPHD的高斯混合实现方法,称之为高斯混合扩展目标概率假设密度滤波,简称GM-EPHD滤波方法。该方法能够较好地解决线性高斯条件下多个扩展目标的跟踪问题。此外,Granstrom等还提出了一种简单的划分量测集的方法。随后,连峰等在“连峰,韩崇昭,刘伟峰,元向辉。高斯混合扩展目标概率假设密度滤波器的收敛性分析.自动化学报.2012(08):1343-1352.”一文中对扩展目标的GM-PHD滤波方法的收敛性进行了分析。在放松了量测率服从泊松分布条件的假设下,Orguner U和Lundquist C等给出了扩展目标的GM-CPHD滤波方法。与CPHD方法一样,该方法虽然能提高扩展目标的跟踪精度,但其过程较为复杂,且时间复杂度增大。Feldman等学者采用随机矩阵的方法为扩展目标建立椭圆形轮廓,然后将其看作状态变量的一部分进行估计和跟踪。随后,Wieneke等将上述方法推广到多个扩展目标的跟踪情形,使用期望最大化方法对扩展目标的椭圆形状和运动情况进行递推估计。
已有的扩展目标GM-PHD滤波方法只适用于线性高斯系统,扩展目标的扩展卡尔曼概率假设密度(Extended Kalman EPHD,EK-EPHD)滤波方法也只在系统为弱非线性时,能取得较好的滤波效果,而对于强非线性系统,其误差性能退化。针对该问题,本发明方法在扩展目标GM-PHD滤波算法的基础上,结合容积卡尔曼滤波(Cubature Kalman Filter,CKF)算法,给出了一种适用于非线性系统下的扩展目标跟踪方法,该方法能够解决滤波系统中非线性函数雅可比矩阵不存在或难以求解的跟踪问题,为非线性系统下扩展目标的跟踪提供了一种新的实现途径。
发明内容
本发明的目的在于提供一种基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法,解决EK-EPHD处理不了的非线性函数雅克比矩阵不存在或难以求解的跟踪问题。
本发明所采用的技术方案是,基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法,具体按照以下步骤实施:
步骤1、预先设定k-1时刻后验强度的高斯混合形式,得到第i个高斯项的均值和协方差;
步骤2、对步骤1得到的高斯混合形式中第i个高斯项的权值、均值和协方差进行预测;
本发明的特点还在于,
步骤1中的高斯混合形式为:
步骤2具体按照以下步骤实施:
步骤2.3、根据状态方程计算第i个高斯项的第l个预测容积点(l=1,2,...,m),具体按照以下算法实施:
步骤2.4、计算得到存活目标的进一步预测均值和方差,具体按照以下算法实施:
式(7)中,Qk-1为过程噪声方差;
步骤2.5、对各高斯项的权值进行进一步预测,具体按照以下算法实施:
对于新生目标的高斯项,其权值预测为:
对于衍生目标的高斯项,其权值预测为:
对于存活目标的高斯项,其权值预测为:
式(10)中,pS表示目标存活概率;
根据上述的算法,进一步预测强度的高斯混合形式为:
式(11)中,Jk|k-1表示一步预测的高斯项数。
步骤3具体按照以下步骤实施:
步骤3.1、量测更新分为两部分:一部分为未检测到的目标强度,另一部分为检测到的目标强度;
若更新未检测到目标强度,则按照以下算法实施:
若未检测到的目标后验强度可表示为:
若检测到的目标强度,具体按照以下算法实施:
则新的容积点为:
式(15)中,标准容积点ξl标号l=1,2,…,m,而m等于状态向量维数的2倍;
步骤3.2、利用量测函数计算出传递的容积点,具体按照以下算法实施:
式(16)中,h(·)为目标的量测函数;
式(17)中,absW表示当前时刻的每个量测划分里的元胞W所包含的量测的个数,其中,当前时刻的每个量测划分里的元胞W为已知量,repeat(A,m,n)函数表示创建一个m×n矩阵,矩阵中的每个值都为A;
步骤3.4、估计量测预测,具体按照以下算法实施:
步骤3.5、更新积分点的状态和协方差,具体按照以下算法实施:
式(25)中,blkdiag(.)表示块对角矩阵,|Wk|表示当前元胞Wk的所包含的量测值的个数;
步骤3.6、更新相应的权值:
式(26)中,wp是划分的权重,且满足:
式(27)中,p和p'均表示对量测集Z的划分,W和W'分别表示p和p'中的一个子集元胞;
Г(i)=e-γ(i)(γ(i))|W| (28);
式(28)中,|W|表示集合的势;
式(29)中,δ|W|,1表示若|W|=1,则δ|W|,1=1,否则,δ|W|,1=0;
量测更新后的后验强度为如下的高斯混合形式:
本发明的有益效果在于:
(1)本发明的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法针对非线性系统下的扩展目标跟踪,不需要计算非线性函数的雅克比矩阵,滤波精度与EK-EPHD方法几乎相同;
(2)本发明的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法能够解决EK-EPHD处理不了的非线性函数雅克比矩阵不存在或难以求解的跟踪问题,该方法可以用于非线性系统中的扩展目标跟踪。
附图说明
图1是单个目标的真实轨迹及分别采用本发明的CK-EPHD滤波方法及EK-EPHD滤波方法的滤波估计对比图;
图2是在整个仿真时间里的目标量测值图谱;
图3是采用传统EK-GMPHD滤波方法估计出来的目标个数图;
图4是在整个时间段内目标的真实个数及分别采用本发明的CK-EPHD滤波方法和传统的EK-EPHD滤波方法估计得到的结果对比图;
图5是采用本发明的CK-EPHD滤波方法和传统的EK-EPHD滤波方法对应的OSPA距离对比图;
图6是多个扩展目标的真实轨迹与采用本发明的CK-EPHD滤波方法、传统的EK-EPHD滤波方法滤波后得到的估计对比图;
图7是在整个仿真时间里的多个扩展目标量测值图谱;
图8是采用传统的EK-GMPHD滤波方法得到的多个扩展目标估计个数图;
图9是目标的真实个数及分别采用本发明的CK-EPHD滤波方法和传统的EK-EPHD滤波方法得到的结果对比图;
图10是在整个仿真时间里本发明的CK-EPHD滤波方法和传统的EK-EPHD滤波方法对应的OSPA距离对比图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
扩展目标与一般意义下的点目标不同,点目标假设每个时刻至多产生一个量测,而扩展目标则假设每一时刻可能会产生多个量测。
通常,在未特别说明的情况下,扩展目标的状态是指它质心的运动状态,具体包括扩展目标质心的位置、速度等;
相应地,k时刻得到的量测集为
量测集由真实扩展目标产生的量测和杂波构成,且二者不可区分,在多个扩展目标的跟踪系统中,扩展目标状态的动态模型为:
其中,νk是零均值高斯白噪声,协方差为Qk,f(.)表示状态转移函数。
量测模型可写为:
其中,wk是零均值高斯白噪声,协方差为Rk,h(.)表示量测函数。
Granstrom等已证明,若扩展目标PHD的预测能够写为高斯混合形式,则量测更新后的PHD也为高斯混合形式,进而各个时刻的预测和更新PHD也都服从高斯混合形式,这一递推实现的过程即为GM-EPHD滤波。这与传统的点目标GM-PHD滤波方法不同,由于扩展目标在每一时刻都会有多个量测,所以在滤波估计前必须对量测集进行划分,然后再进行一步预测和量测更新。
本发明的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法,具体按照以下步骤实施:
步骤1、预先设定k-1时刻后验强度的高斯混合形式为:
步骤2、对步骤1得到的高斯混合形式中第i个高斯项的权值、均值和协方差进行预测:
步骤2.3、根据状态方程计算第i个高斯项的第l个预测容积点(l=1,2,...,m),具体按照以下算法实施:
步骤2.4、计算得到存活目标的进一步预测均值和方差,具体按照以下算法实施:
式(7)中,Qk-1为过程噪声方差;
步骤2.5、对各高斯项的权值进行进一步预测,具体按照以下算法实施:
对于新生目标的高斯项,其权值预测为:
对于衍生目标的高斯项,其权值预测为:
对于存活目标的高斯项,其权值预测为:
式(10)中,pS表示目标存活概率;
根据上述的算法,进一步预测强度的高斯混合形式为:
式(11)中,Jk|k-1表示一步预测的高斯项数;
步骤3.1、量测更新分为两部分:一部分为未检测到的目标强度,另一部分为检测到的目标强度;
若更新未检测到目标强度,则按照以下算法实施:
若未检测到的目标后验强度可表示为:
若检测到的目标强度,具体按照以下算法实施:
则新的容积点为:
式(15)中,标准容积点ξl标号l=1,2,…,m,而m等于状态向量维数的2倍;
步骤3.2、利用量测函数计算出传递的容积点,具体按照以下算法实施:
式(16)中,h(·)为目标的量测函数;
式(17)中,absW表示当前时刻的每个量测划分里的元胞W所包含的量测的个数,其中,当前时刻的每个量测划分里的元胞W为已知量,repeat(A,m,n)函数表示创建一个m×n矩阵,矩阵中的每个值都为A;
步骤3.4、估计量测预测,具体按照以下算法实施:
步骤3.5、更新积分点的状态和协方差,具体按照以下算法实施:
式(25)中,blkdiag(.)表示块对角矩阵,|Wk|表示当前元胞Wk的所包含的量测值的个数;
步骤3.6、更新相应的权值:
式(26)中,wp是划分的权重,且满足:
式(27)中,p和p'均表示对量测集Z的划分,W和W'分别表示p和p'中的一个子集元胞,采用右上标“'”来对同一公式中分子分母出现的具有相同含义的符号进行区分;
Г(i)=e-γ(i)(γ(i))|W| (28);
式(28)中,|W|表示集合的势;
式(29)中,δ|W|,1表示若|W|=1,则δ|W|,1=1,否则,δ|W|,1=0;
量测更新后的后验强度为如下的高斯混合形式:
本发明的方法结合了容积卡尔曼滤波与GM-EPHD算法,是一种新的适用于非线性系统下的CK-EPHD滤波方法。
由于高斯项数目会随着跟踪步数的增加而急速增加,计算量不断增大。为了降低计算量,可以通过控制每一步的高斯项数目解决该问题。合理控制高斯项数目的典型方法如文献“Ba-Ngu Vo,Wing-Kin Ma.The Gaussianmixture probability hypothesis density filter,IEEE Transactions on SignalProcessing,2006,54(11):4091-4104.”中所述的合并与裁剪方法。
以下通过实施例验证本发明的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法的有效性和实用性:
采用对比实验的形式,在二维空间中分别给出单个扩展目标和多个扩展目标的跟踪实验,具体的仿真条件详见每个实验的描述:
实施例1:单个扩展目标的跟踪实验
预先设定目标的状态转移方程与量测方程分别为:
xk=Fkxk-1+Gkνk (33);
当目标做匀速运动,即采用匀速模型时,
当目标做常数转弯运动,即采用匀速转弯模型时,
状态噪声转移矩阵为:
观测噪声νk和量测噪声wk都是零均值高斯白噪声,其协方差分别为Qk=diag([0.5,0.5]),量测噪声标准差σθ=5×(π/180),σr=8,采样周期T=1s,转弯速率Ω=(π/80)rad/s,量测向量zk包括偏转角和径向距离两个分量。
设定新生目标随机集的PHD为:
γk(x)=0.1N(x;mγ,Pγ);
上式中,mγ=[250 6 250 -15]T,Pγ=diag([100 25 100 25]);
杂波均匀分布于观测区域,数目服从均值为λk=10的泊松分布,扩展目标的量测值个数服从均值为γ(i)=3的泊松分布。
目标的存活概率和检测概率分别为pS=0.99和pD=0.90,合并门限Uprun=4,修剪门限Tprun=1e-5,势误差与状态误差的调节因子cospa=70,最优子模型分配距离取2阶距离pospa=2,最大高斯项数为Jmax=100,整个仿真时间为100s,监视区域为[-π/2,π/2]rad×[0,1500]m。
如图1所示,图中显示了目标的真实轨迹与两种方法的滤波估计结果,由图1可知,使用本发明的CK-EPHD滤波方法不仅能达到EK-EPHD的滤波效果,而且能够去除掉EK-EPHD所不能达到的杂波,例如,在位置(1200,-300)左右有三个杂波点,用EK-EPHD方法不能过滤掉,但本发明的CK-EPHD滤波方法却能达到;如图2所示,是在整个监控区域的目标量测集。
如图3所示,给出了用传统EK-GMPHD滤波方法估计出来的目标个数,从图3中我们可看出传统的滤波方法完全不能准确地估计目标个数,其平均估计目标数约为10,而真实目标数为1。
如图4所示,图4中显示了在整个时间段内目标的真实个数及本发明的CK-EPHD滤波和传统EK-EPHD滤波的估计结果,由图中可知,本发明的CK-EPHD的估计结果与传统EK-EPHD的基本一致,在54时刻这两种滤波方法的估计结果都出现了误差,但本发明的CK-EPHD滤波方法在下一时刻就能够正确的估计出结果,而EK-EPHD依旧估计错误。
如图5所示,给出本发明的CK-EPHD滤波方法和传统的EK-EPHD滤波方法对应的OSPA距离,可以从图中看到,两种方法的OSPA距离也基本一致,其中在55时刻左右的CK-EPHD的OSPA距离要比EK-EPHD的低一点;并且由仿真实验可知,在整个仿真步骤里,CK-EPHD方法的平均OSPA距离要比EK-EPHD的略小些。
实施例2
多个扩展目标的跟踪
目标的状态转移方程和量测方程与实施例1相同,参数设置也与其相同;
设定在整个仿真区域有3个目标,目标1在k=1时刻出现,消亡于k=100时刻;目标2在k=11时刻出现并消亡于k=100时刻;目标3在k=66时刻出现并消亡于k=100时刻;目标1和目标2均做匀速运动,目标3做转弯运动;
设新生目标随机集的PHD为:
上式中, Pγ=diag([100 25 100 25]T);
杂波均匀分布于整个观测区域,数目服从均值为λk=5的泊松分布,量测个数服从γ(i)=5的泊松分布,其他各仿真参数与单个扩展目标实验相同。
由图6可知,本发明的CK-EPHD滤波方法能够达到与传统EK-EPHD滤波方法相同的效果。
如图7所示,给出了在监控区域内的量测集。
如图8所示,图中显示了用传统EK-GMPHD滤波方法得到的多个扩展目标估计个数,从图中可以看到,它的滤波效果很差,其估计目标数始终比真实目标数大很多,即过估计。
如图9所示,图中给出了目标的真实个数及本发明的CK-EPHD滤波方法和EK-EPHD滤波方法的估计结果,由图中可知,本发明的CK-EPHD滤波方法得到的估计目标数与传统的EK-EPHD滤波方法的滤波结果几乎完全一致,只在23时刻CK-EPHD估计出了准确的数目,而传统的EK-EPHD滤波方法没有。
如图10所示,图中给出了两种方法的OSPA距离,由图10可知,两种方法的OSPA距离几乎一样,只在少有的几个时刻本发明的CK-EPHD滤波方法的OSPA距离比EK-EPHD的稍低。由仿真实验可知,在整个仿真步骤里,本发明的CK-EPHD滤波方法的平均OSPA距离比传统的EK-EPHD滤波方法的略小些。
将以上两个实施例进行对比,得出以下结果;
不论目标的轨迹是直线运动还是转弯运动,本发明的CK-EPHD方法都能够达到传统EK-EPHD方法的滤波效果,并且其平均OSPA距离都比EK-EPHD的要低一点儿。
Claims (4)
3.根据权利要求1所述的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法,其特征在于,所述步骤2具体按照以下步骤实施:
式中,ξl为第l个标准容积点,[1]l表示第l列,[1]与状态向量的维数有关,当状态向量维数为2时,[1]表示下列点集: 标准容积点的标号l=1,2,…,m,而m等于状态向量维数的2倍;
步骤2.3、根据状态方程计算第i个高斯项的第l个预测容积点(l=1,2,...,m),具体按照以下算法实施:
步骤2.4、计算得到存活目标的进一步预测均值和方差,具体按照以下算法实施:
式(7)中,Qk-1为过程噪声方差;
步骤2.5、对各高斯项的权值进行进一步预测,具体按照以下算法实施:
对于新生目标的高斯项,其权值预测为:
对于衍生目标的高斯项,其权值预测为:
对于存活目标的高斯项,其权值预测为:
式(10)中,pS表示目标存活概率;
根据上述的算法,进一步预测强度的高斯混合形式为:
式(11)中,Jk|k-1表示一步预测的高斯项数。
4.根据权利要求1所述的基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法,其特征在于,所述步骤3具体按照以下步骤实施:
步骤3.1、量测更新分为两部分:一部分为未检测到的目标强度,另一部分为检测到的目标强度;
若更新未检测到目标强度,则按照以下算法实施:
若未检测到的目标后验强度可表示为:
若检测到的目标强度,具体按照以下算法实施:
则新的容积点为:
式(15)中,标准容积点ξl标号l=1,2,…,m,而m等于状态向量维数的2倍;
步骤3.2、利用量测函数计算出传递的容积点,具体按照以下算法实施:
式(16)中,h(·)为目标的量测函数;
式(17)中,absW表示当前时刻的每个量测划分里的元胞W所包含的量测的个数,其中,当前时刻的每个量测划分里的元胞W为已知量,repeat(A,m,n)函数表示创建一个m×n矩阵,矩阵中的每个值都为A;
步骤3.4、估计量测预测,具体按照以下算法实施:
步骤3.5、更新积分点的状态和协方差,具体按照以下算法实施:
式(20)中,[z1,…,z|W|]T为相应元胞W中的真实量测值,为滤波增益,具体按照以下算法实施:
式(25)中,blkdiag(.)表示块对角矩阵,|Wk|表示当前元胞Wk的所包含的量测值的个数;
步骤3.6、更新相应的权值:
式(26)中,wp是划分的权重,且满足:
式(27)中,p和p'均表示对量测集Z的划分,W和W'分别表示p和p'中的一个子集元胞;
Г(i)=e-γ(i)(γ(i))|W| (28);
式(28)中,|W|表示集合的势;
式(29)中,δ|W|,1表示若|W|=1,则δ|W|,1=1,否则,δ|W|,1=0;
量测更新后的后验强度为如下的高斯混合形式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310753640.5A CN103729637B (zh) | 2013-12-31 | 2013-12-31 | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310753640.5A CN103729637B (zh) | 2013-12-31 | 2013-12-31 | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103729637A true CN103729637A (zh) | 2014-04-16 |
CN103729637B CN103729637B (zh) | 2017-01-11 |
Family
ID=50453702
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310753640.5A Expired - Fee Related CN103729637B (zh) | 2013-12-31 | 2013-12-31 | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103729637B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN104778358A (zh) * | 2015-04-09 | 2015-07-15 | 西安工程大学 | 多传感器存在监测区域部分重叠的扩展目标跟踪方法 |
CN104849702A (zh) * | 2015-04-30 | 2015-08-19 | 中国民航大学 | 利用ads-b数据的gm-ephd滤波雷达系统误差联合估计方法 |
CN105321379A (zh) * | 2014-07-31 | 2016-02-10 | 霍尼韦尔国际公司 | 基于传感器轨迹id在phd滤波器中合并强度 |
CN106372646A (zh) * | 2016-08-30 | 2017-02-01 | 上海交通大学 | 基于srck‑gmcphd滤波的多目标跟踪方法 |
CN106443661A (zh) * | 2016-09-08 | 2017-02-22 | 河南科技大学 | 基于无迹卡尔曼滤波的机动扩展目标跟踪方法 |
CN106526585A (zh) * | 2016-10-26 | 2017-03-22 | 中国人民解放军空军工程大学 | 基于高斯粒子势概率假设密度滤波的目标检测前跟踪方法 |
CN106556818A (zh) * | 2016-11-18 | 2017-04-05 | 辽宁工业大学 | 一种用于单目标跟踪的低计算复杂度贝努利滤波器 |
CN107703504A (zh) * | 2017-10-12 | 2018-02-16 | 中国电子科技集团公司第二十九研究所 | 一种基于随机集的多点定位目标跟踪方法 |
CN108363054A (zh) * | 2018-02-07 | 2018-08-03 | 电子科技大学 | 用于单频网络和多路径传播的被动雷达多目标跟踪方法 |
CN108519595A (zh) * | 2018-03-20 | 2018-09-11 | 上海交通大学 | 联合多传感器配准与多目标跟踪方法 |
CN108717702A (zh) * | 2018-04-24 | 2018-10-30 | 西安工程大学 | 基于分段rts的概率假设密度滤波平滑方法 |
WO2019047455A1 (zh) * | 2017-09-08 | 2019-03-14 | 深圳大学 | 一种适用于闪烁噪声的多机动目标跟踪方法及系统 |
CN110488273A (zh) * | 2019-08-30 | 2019-11-22 | 成都纳雷科技有限公司 | 一种基于雷达的车辆跟踪检测方法及装置 |
US10605607B2 (en) | 2014-07-31 | 2020-03-31 | Honeywell International Inc. | Two step pruning in a PHD filter |
US11175142B2 (en) | 2014-07-31 | 2021-11-16 | Honeywell International Inc. | Updating intensities in a PHD filter based on a sensor track ID |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065037A (zh) * | 2012-11-13 | 2013-04-24 | 杭州电子科技大学 | 非线性系统基于分散式容积信息滤波的目标跟踪方法 |
CN103217175A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 一种自适应容积卡尔曼滤波方法 |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
-
2013
- 2013-12-31 CN CN201310753640.5A patent/CN103729637B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065037A (zh) * | 2012-11-13 | 2013-04-24 | 杭州电子科技大学 | 非线性系统基于分散式容积信息滤波的目标跟踪方法 |
CN103217175A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 一种自适应容积卡尔曼滤波方法 |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
Non-Patent Citations (4)
Title |
---|
王品等: "一种非线性GM-PHD滤波新方法", 《电子学报》 * |
连峰等: "高斯混合扩展目标概率假设密度滤波器的收敛性分析", 《自动化学报》 * |
陈里铭: "基于随机有限集的多目标跟踪方法研究", 《中国博士学位论文全文数据库信息科技辑》 * |
陈里铭等: "基于数值积分卡尔曼-概率假设密度滤波的多说话人跟踪方法", 《信号处理》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105321379A (zh) * | 2014-07-31 | 2016-02-10 | 霍尼韦尔国际公司 | 基于传感器轨迹id在phd滤波器中合并强度 |
US10605607B2 (en) | 2014-07-31 | 2020-03-31 | Honeywell International Inc. | Two step pruning in a PHD filter |
US11175142B2 (en) | 2014-07-31 | 2021-11-16 | Honeywell International Inc. | Updating intensities in a PHD filter based on a sensor track ID |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN104778358A (zh) * | 2015-04-09 | 2015-07-15 | 西安工程大学 | 多传感器存在监测区域部分重叠的扩展目标跟踪方法 |
CN104778358B (zh) * | 2015-04-09 | 2017-11-10 | 西安工程大学 | 多传感器存在监测区域部分重叠的扩展目标跟踪方法 |
CN104849702A (zh) * | 2015-04-30 | 2015-08-19 | 中国民航大学 | 利用ads-b数据的gm-ephd滤波雷达系统误差联合估计方法 |
CN104849702B (zh) * | 2015-04-30 | 2017-10-27 | 中国民航大学 | 利用ads‑b数据的gm‑ephd滤波雷达系统误差联合估计方法 |
CN106372646A (zh) * | 2016-08-30 | 2017-02-01 | 上海交通大学 | 基于srck‑gmcphd滤波的多目标跟踪方法 |
CN106372646B (zh) * | 2016-08-30 | 2020-07-14 | 上海交通大学 | 基于srck-gmcphd滤波的多目标跟踪方法 |
CN106443661A (zh) * | 2016-09-08 | 2017-02-22 | 河南科技大学 | 基于无迹卡尔曼滤波的机动扩展目标跟踪方法 |
CN106443661B (zh) * | 2016-09-08 | 2019-07-19 | 河南科技大学 | 基于无迹卡尔曼滤波的机动扩展目标跟踪方法 |
CN106526585B (zh) * | 2016-10-26 | 2019-01-25 | 中国人民解放军空军工程大学 | 基于高斯粒子势概率假设密度滤波的目标检测前跟踪方法 |
CN106526585A (zh) * | 2016-10-26 | 2017-03-22 | 中国人民解放军空军工程大学 | 基于高斯粒子势概率假设密度滤波的目标检测前跟踪方法 |
CN106556818B (zh) * | 2016-11-18 | 2018-11-23 | 辽宁工业大学 | 一种用于单目标跟踪的低计算复杂度贝努利滤波器 |
CN106556818A (zh) * | 2016-11-18 | 2017-04-05 | 辽宁工业大学 | 一种用于单目标跟踪的低计算复杂度贝努利滤波器 |
WO2019047455A1 (zh) * | 2017-09-08 | 2019-03-14 | 深圳大学 | 一种适用于闪烁噪声的多机动目标跟踪方法及系统 |
CN107703504A (zh) * | 2017-10-12 | 2018-02-16 | 中国电子科技集团公司第二十九研究所 | 一种基于随机集的多点定位目标跟踪方法 |
CN108363054A (zh) * | 2018-02-07 | 2018-08-03 | 电子科技大学 | 用于单频网络和多路径传播的被动雷达多目标跟踪方法 |
CN108363054B (zh) * | 2018-02-07 | 2021-07-02 | 电子科技大学 | 用于单频网络和多路径传播的被动雷达多目标跟踪方法 |
CN108519595A (zh) * | 2018-03-20 | 2018-09-11 | 上海交通大学 | 联合多传感器配准与多目标跟踪方法 |
CN108717702B (zh) * | 2018-04-24 | 2021-08-31 | 西安工程大学 | 基于分段rts的概率假设密度滤波平滑方法 |
CN108717702A (zh) * | 2018-04-24 | 2018-10-30 | 西安工程大学 | 基于分段rts的概率假设密度滤波平滑方法 |
CN110488273A (zh) * | 2019-08-30 | 2019-11-22 | 成都纳雷科技有限公司 | 一种基于雷达的车辆跟踪检测方法及装置 |
CN110488273B (zh) * | 2019-08-30 | 2022-08-05 | 成都纳雷科技有限公司 | 一种基于雷达的车辆跟踪检测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103729637B (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103729637A (zh) | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 | |
CN101770024B (zh) | 多目标跟踪方法 | |
CN104237879B (zh) | 一种雷达系统中的多目标跟踪方法 | |
CN103472445B (zh) | 一种针对多目标场景的检测跟踪一体化方法 | |
CN103345577B (zh) | 变分贝叶斯概率假设密度多目标跟踪方法 | |
CN101975575B (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
Yang et al. | Comparison of unscented and extended Kalman filters with application in vehicle navigation | |
CN110503071A (zh) | 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法 | |
CN105182291A (zh) | 自适应目标新生强度的phd平滑器的多目标跟踪方法 | |
CN103759732B (zh) | 一种角度信息辅助的集中式多传感器多假设跟踪方法 | |
CN104155650A (zh) | 一种基于熵权值法点迹质量评估的目标跟踪方法 | |
CN111722214B (zh) | 雷达多目标跟踪phd实现方法 | |
CN105761276B (zh) | 基于迭代ransac自适应新生目标强度估计的gm-phd多目标跟踪方法 | |
CN109886305A (zh) | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 | |
CN103017771B (zh) | 一种静止传感器平台的多目标联合分配与跟踪方法 | |
CN106407677A (zh) | 一种测量数据丢失情况下的多目标跟踪方法 | |
CN104713560A (zh) | 基于期望最大化的多源测距传感器空间配准方法 | |
CN104833949A (zh) | 一种基于改进距离参数化的多无人机协同无源定位方法 | |
CN104268567A (zh) | 一种观测数据聚类划分的扩展目标跟踪方法 | |
CN107797106A (zh) | 一种加速em未知杂波估计的phd多目标跟踪平滑滤波方法 | |
CN106054167A (zh) | 基于强度滤波器的多扩展目标跟踪方法 | |
CN115204212A (zh) | 一种基于stm-pmbm滤波算法的多目标跟踪方法 | |
CN107102293A (zh) | 基于滑窗累积密度估计的未知杂波无源协同定位方法 | |
CN111711432B (zh) | 一种基于ukf和pf混合滤波的目标跟踪算法 | |
CN109509207A (zh) | 一种对点目标和扩展目标进行无缝跟踪的方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170111 Termination date: 20181231 |