CN108121846B - 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 - Google Patents
一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 Download PDFInfo
- Publication number
- CN108121846B CN108121846B CN201611106457.6A CN201611106457A CN108121846B CN 108121846 B CN108121846 B CN 108121846B CN 201611106457 A CN201611106457 A CN 201611106457A CN 108121846 B CN108121846 B CN 108121846B
- Authority
- CN
- China
- Prior art keywords
- clutter
- estimation
- phd
- target
- entropy
- 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 72
- 239000000203 mixture Substances 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims description 21
- 238000009826 distribution Methods 0.000 claims description 16
- 238000002156 mixing Methods 0.000 claims description 12
- 238000012952 Resampling Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 5
- 239000002245 particle Substances 0.000 description 8
- 238000001914 filtration Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000007476 Maximum Likelihood Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 230000008033 biological extinction Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- PXUQTDZNOHRWLI-OXUVVOBNSA-O malvidin 3-O-beta-D-glucoside Chemical compound COC1=C(O)C(OC)=CC(C=2C(=CC=3C(O)=CC(O)=CC=3[O+]=2)O[C@H]2[C@@H]([C@@H](O)[C@H](O)[C@@H](CO)O2)O)=C1 PXUQTDZNOHRWLI-OXUVVOBNSA-O 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公布了一种基于熵惩罚的EM未知杂波估计的PHD多目标跟踪方法,主要解决了概率假设密度多目标跟踪方法中存在的杂波强度未知的问题。所示方法包括:首先采用有限混合模型对未知杂波密度建模;其次分别对混合权重及缺失参数施加熵惩罚因子;然后通过自适应动态系数调节,使得混合模型低权值分量加速消亡,减少了迭代次数,且对初始参数不敏感。该方法在杂波强度未知的环境下,具有跟踪精度高且稳定的优势,提高了PHD滤波器在多目标跟踪中的性能。
Description
技术领域
本发明属于目标跟踪技术领域,特别是未知杂波估计的PHD跟踪方法。
背景技术
多目标跟踪(Multi-Target Tracking,MTT)是目标跟踪的研究热点,其在国防、军事、环境监测、安全监控、智能交通等领域具有广泛的应用价值。多目标跟踪方法主要有基于数据关联和随机有限集两类,前者的主要思想是通过数据关联建立观测信息与目标身份信息的匹配,从而将多目标跟踪问题转化为对单一目标状态的估计。但由于数据关联需要大量的计算,所以这类方法处理起来极为复杂。针对这一问题,Mahler将随机有限集理论(Random Finite Set,RFS)引入到多目标跟踪问题的求解中,提出了基于概率假设密度(Probability Hypothesis Density,PHD)滤波方法。PHD滤波器是RFS理论框架下多目标状态集合的后验概率密度的一阶统计矩,其实现主要分为适用于非线性非高斯情形下的序贯蒙特卡洛PHD(Sequential Monte Carlo PHD,SMC-PHD)滤波器和适用于线性高斯情形下的高斯混合PHD(Gaussian Mixture PHD,GM-PHD)滤波器。这类方法避免了数据关联,从而降低了计算难度,成为了近年来国内外研究热点。
在传统的PHD滤波器中,通常假定杂波的数目服从泊松分布,且在整个区域分布均匀。此种假定在杂波干扰较少的跟踪环境中不会有太大问题,但在复杂环境下,由于杂波个数和分布未知,这种假设使方法的鲁棒性和跟踪精度会严重降低。针对此类问题,目前的研究文献中有很多解决方法,但都存在着各种不足,林晓东等人于2011年提出了改进的基于概率假设密度滤波的多目标跟踪方法,该方法在目标状态中加入标签,除去权重较大标签代表的杂波干扰,此方法能更好地避免杂波干扰,但对新的目标不易发现。Lian F等人于2010年提出了未知杂波强度估计的PHD滤波方法,该方法运用有限混合模型和期望最大化(Expectation Maximization,EM)方法,较全面地解决了多目标跟踪中未知杂波估计的问题,但是杂波模型对初始值要求较敏感,不合理的初始值将影响估计精度。另外,传统EM方法计算量大,收敛速度慢,在复杂杂波背景下进行跟踪难以满足实时性的要求。
发明内容
本发明所解决的技术问题在于提供一种基于熵惩罚EM的未知杂波估计的PHD(Entropy-penalized Expectation Maximization PHD,EPEM-PHD)滤波方法。
实现本发明目的的技术解决方案为:首先,提出了将带有熵惩罚因子的EM方法应用到PHD滤波器中,以加速PHD方法执行效率,提高了方法的实时性。其次,通过动态系数自适应调节使得加速混合模型分量消亡,在减少方法迭代次数的同时,又避免陷入局部最优解。此外,熵惩罚EM方法还具有对初始参数不敏感的优点。通过SMC-PHD滤波器的实现,保证了该方法在非线性情况下也有很好的应用。
本发明与现有技术相比具有以下优点:
1.针对多目标跟踪方法对强实时性的需求,首次提出了将带有熵惩罚因子的EM方法应用到PHD滤波器中,以加速PHD方法执行速度,使得改进后的PHD方法在多目标跟踪中具有实时性强的优势。
2.在熵惩罚EM方法中赋予熵惩罚因子合适的动态调节系数,通过动态系数自适应调节使得加速混合模型分量消亡,在减少方法迭代次数的同时,又避免陷入局部最优解。
3.提出了基于序贯Monte Carlo改进EM的未知杂波估计PHD方法,在完成杂波强度估计的同时,保证了在非线性情况下也能有很好的应用。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为本发明方法的流程图。
图2为无线传感器网络协同跟踪仿真场景图中X坐标真实值与估计值对比。
图3为无线传感器网络协同跟踪仿真场景图中Y坐标真实值与估计值对比。
图4为不同方法跟踪误差对比图。
具体实施方式
依据附图,对本发明的技术方案作具体说明。
所述基于熵惩罚EM的未知杂波估计的PHD滤波方法,包括以下步骤:
步骤1进行基于熵惩罚EM方法的未知杂波估计。
采用随机有限集为杂波强度建模,杂波强度可表示为
k(z)=λc(z) (1)
其中,λ为杂波数目,c(z)为杂波分布。
(1)杂波数目的估计
记H={h1,h2…hk-1},hi表示i时刻的量测个数,i=1,2,…,k-1。设杂波的运动为泊松点过程,即
其中,λ即为平均每帧的杂波个数。其无偏最小方差估计为
(2)杂波分布的估计
采用有限混合模型描述杂波分布
其中,M为分量数目,θm={μm,φm}为分量均值μm与分量方差φm构成的分量m参数集合,θ={θ1,…,θM,π1,…,πM}为混合参数集合,混合权重πm表示混合模型中各分量的比例系数,并且
在有限混合模型的初始步骤中增加一个均匀分布分量,用于拟合来源于背景噪声的稀疏杂波分布,既可实现初始分量对目标区域较为全面的表示,同时又能达到缩减初始分量数目的目的。新的有限混合模型表示如下
由量测集合Z的n个元素相互独立,得到随机有限集合参数θ的极大似然函数为
采用EM方法求参数θ的极大似然估计。令EM方法的缺失参数集合为E={e1,e2,…,en},其中ej为M维向量,emj=0或1,分别表示第j个测量属于或不属于混合模型中的第m分量,且满足
由完整数据集Y={z1,z2,…,zn,E},可得最大似然法求解的相关log似然函数为
由混合模型的各分量均为高斯分布,可得混合模型为
其中,d表示量测维数,δ(z,μ;φ)表示z与μ的马氏距离,emi为缺失参数,μm、φm、πm为混合参数。
步骤2进行基于熵惩罚的改进EM的参数估计
未施加惩罚因子的EM方法通常步骤为
E-Step:求解参数emi,有
M-Step:将似然函数式(7)最大化,以获得混合参数πm、μm和φm的估计值,有
混合权重
分量均值
分量方差
对参数emi与混合权重πm施加惩罚因子,得到新的EM方法。具体为
(1)对混合权重πm施加熵惩罚因子,得到新的M-Step。
混合权重πm是迭代过程中分量删减的依据,πm值越大,则该分量在下一次迭代中越可能被留存下来。而当πm越小,直至小于预先设定的阈值ξ时,该分量就会被删减。取表示所有分量的信息熵。其中,-lnπm表示第m分量中一个数据点所散发的信息量,当时,则认为没有关于πm的信息量。这里通过最大化W,即最小化信息熵实现更好更快地获得πm的信息。在似然函数式(7)中添加惩罚因子以便对EM方法的迭加过程进行调节,可得
其中,α为自适应调节系数,且满足α∈(0,1),有
采用拉格朗日乘子法,将式(5)代入式(15),并对πm求一阶偏导,有
由式(17)、式(13)与式(14)形成新的M-Step。
(2)对缺失参数emi施加熵惩罚因子,得到新的E-Step。
由式(20)形成施加熵惩罚后的新的E-Step。
利用改进后的E-Step和M-Step交替迭代直到收敛,即得到有限混合模型缺失参数emi和混合参数μm、φm和πm的估计结果。
步骤3,得到基于熵惩罚EM的未知杂波估计的PHD滤波方法
根据式(10)和有限混合模型参数估计结果,可得杂波分布c(z)。
PHD算法中,SMC-PHD主要步骤为:
(1)粒子预测
其中,φ(xk,xk-1)=ek|k-1(xk-1)fk|k-1(xk|xk-1)+β(xk|xk-1),为k时刻第i个目标的状态,Jk为出生目标粒子数目。Zk={zk,1,zk,2,...,zk,j,...,zk,n(k)},zk,j表示k时刻第j个量测值,fk|k-1(xk|xk-1)为状态转移密度,ek|k-1(xk)为目标生存概率,β(xk|xk-1)为源自状态xk-1的衍生目标强度,γk(xk)为出生目标强度;
(2)粒子更新
(3)重采样
第1步:对粒子权值求和
将杂波分布c(z)代入式(1)可得估计值k(z),将k(z)代入式(23)和式(24),即可得到EPEM-PHD滤波器。
下面对本发明的方法进行仿真验证:
仿真实验采用序贯Monte Carlo概率假设密度滤波器(SMC-PHD)与所提出的EPEM-PHD比较跟踪性能。仿真实验在AMD FX-7500 2.10GHz处理器、1024MB内存的PC上,使用MATLAB R2012b平台实现。模拟一个[-120,120]×[-100,100]m的x-y二维平面区域,先后出现4个作无规则运动的目标。目标的运动模型为
xk=F·xk-1+G·ωk (22)
其中,xk-1=[xx,vx,xy,vy]T为目标状态,由x,y方向的目标位置及速度[xx,vx]T和[xy,vy]T组成;F为状态转移矩阵,G为噪声输入矩阵,ωk为过程噪声,取高斯白噪声ωk~N(0,Q),F与G分别表示为
目标量测方程Zk为
其中,{V1,k}和{V2,k}为相互独立的零均值高斯白噪声。***采样周期T=1s,仿真时间步长为40。出生目标强度为γk=0.2N(·;x,Q),其正态分布的均值和方差为
目标存活概率Ps(x)取0.95,目标检测概率PD(x)设为1。令平均每帧杂波个数的真实值为λ=50;杂波在设定区域内随机产生。每个目标赋予1000个粒子,随机生成一个随机数作为每帧杂波个数。
实验中所选择的四个目标之中,目标1在第1秒出现,终止于第7秒;目标2在第8秒出现,终止于第25秒;目标3在第12秒出现,终止于第37秒;目标4在第26秒出现,终止于第40秒。实验结果如图2至图4所示。
图2和图3显示了单次Monte Carlo仿真中的X坐标和Y坐标的跟踪效果,从图中可以看出,改进方法能够正确估计多目标状态,且优于PHD方法。跟踪第三个目标在第25秒至第33秒时PHD方法偏离了目标的真实位置,但EPEM-PHD却很好地进行目标跟踪,性能提高的原因在于进行了杂波强度的在线估计。采用Wasserstein距离评价多目标跟踪性能,其表达形式为:
图4显示了通过50次Monte Carlo仿真后,本发明方法和传统PHD多目标跟踪方法的Wasserstein距离。从图中可以看出,本发明方法的Wasserstein距离较小,说明EPEM-PHD方法在多目标跟踪精度上的准确性高于PHD,并且能够在杂波未知的环境下,使得估计结果更接近于真实的分布,具有更好的稳定性和鲁棒性,同时也证明了本发明方法的有效性。
同时,***通过50次Monte Carlo仿真后由CPU平均运算时间得到,PHD方法的平均耗时为7.24秒,而EPEM-PHD平均耗时为7.37秒,耗时增幅为1.80%,是因为对未知杂波进行了估计。增幅并不显著,所以本发明方法保持了较好的实时性。
综上所述,本发明公布的一种基于熵惩罚的EM未知杂波估计的PHD多目标跟踪方法在杂波强度未知的环境下,与传统的PHD滤波方法相比,保持了较好的实时性,具有跟踪精度高且稳定的优势。本发明的基于熵惩罚的EM未知杂波估计的PHD多目标跟踪方法在处理杂波未知时的多目标跟踪方面具有积极的意义。
Claims (3)
1.一种基于熵惩罚的EM未知杂波估计的PHD多目标跟踪方法,包括以下步骤:
第一步,初始化获得k-1时刻的多目标状态集及杂波强度;
第二步,根据杂波强度和PHD预测与更新公式,进行PHD更新,并重采样,得到k时刻多目标状态集;
第三步,根据k时刻的目标量测集与k时刻的多目标状态集对比,得到k时刻的杂波量测集;
第四步,杂波强度估计:由杂波量测个数得到杂波数目估计;根据杂波量测集建立有限混合模型,根据熵惩罚EM方法估计模型参数,从而得到杂波分布估计;然后由杂波数目估计和杂波分布估计得到杂波强度估计, 转至第一步,进入下一时刻;
其中根据杂波量测集建立有限混合模型,根据熵惩罚EM方法估计模型参数,从而得到杂波分布估计具体包括:由有限混合模型描述杂波分布并进行估计,针对目标量测zi,初始化混合模型各项参数;根据施加熵惩罚因子的EM方法估计缺失参数emi和混合参数μm、φm、πm:首先,设变量q=1,由更新混合权重πm,n为目标量测集元素个数,如果πm<ξ,ξ为设定的阈值,则删除混合分量m,同时执行和M为分量数目;然后更新协方差矩阵φm,有
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611106457.6A CN108121846B (zh) | 2016-11-29 | 2016-11-29 | 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611106457.6A CN108121846B (zh) | 2016-11-29 | 2016-11-29 | 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108121846A CN108121846A (zh) | 2018-06-05 |
CN108121846B true CN108121846B (zh) | 2022-01-11 |
Family
ID=62226442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611106457.6A Active CN108121846B (zh) | 2016-11-29 | 2016-11-29 | 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108121846B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110309599B (zh) * | 2019-07-04 | 2022-04-12 | 西北工业大学 | 一种方位测量融合与多目标定位方法 |
CN110909312B (zh) * | 2019-12-18 | 2022-04-22 | 哈尔滨工程大学 | 一种应用于rbmcda跟踪算法的目标消亡判断方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076596B (zh) * | 2013-01-17 | 2015-03-04 | 西安电子科技大学 | 基于先验信息的mimo雷达发射方向图设计方法 |
CN103310115B (zh) * | 2013-06-27 | 2015-12-23 | 西安电子科技大学 | 一种多目标跟踪的杂波估计方法 |
US11443331B2 (en) * | 2014-02-13 | 2022-09-13 | Conduent Business Solutions, Llc | Multi-target tracking for demand management |
CN104794735B (zh) * | 2015-04-02 | 2017-08-25 | 西安电子科技大学 | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 |
CN104750998B (zh) * | 2015-04-09 | 2017-08-25 | 西安电子科技大学 | 基于强度滤波器的被动多传感器目标跟踪方法 |
CN105005983A (zh) * | 2015-04-13 | 2015-10-28 | 西南科技大学 | 一种sar图像背景杂波建模及目标检测方法 |
CN105182291B (zh) * | 2015-08-26 | 2017-08-25 | 西安电子科技大学 | 自适应目标新生强度的phd平滑器的多目标跟踪方法 |
CN105184829B (zh) * | 2015-08-28 | 2017-12-12 | 华中科技大学 | 一种紧密空间目标检测和高精度质心定位方法 |
CN105844217A (zh) * | 2016-03-11 | 2016-08-10 | 南京航空航天大学 | 一种基于量测驱动新生目标强度估计的phd多目标跟踪方法 |
CN106054167B (zh) * | 2016-05-27 | 2018-06-26 | 西安电子科技大学 | 基于强度滤波器的多扩展目标跟踪方法 |
-
2016
- 2016-11-29 CN CN201611106457.6A patent/CN108121846B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108121846A (zh) | 2018-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106487358B (zh) | 一种机动目标转弯跟踪方法 | |
CN110503071B (zh) | 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法 | |
CN105354860B (zh) | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 | |
CN106772353B (zh) | 一种适用于闪烁噪声的多目标跟踪方法及*** | |
CN103902829B (zh) | 传递边缘分布及存在概率的目标跟踪方法与目标跟踪*** | |
Yang et al. | Adaptive probability hypothesis density filter based on variational Bayesian approximation for multi‐target tracking | |
CN104881521B (zh) | 一种标号随机集滤波器分布式融合方法 | |
CN111711432B (zh) | 一种基于ukf和pf混合滤波的目标跟踪算法 | |
Chu et al. | Computation approaches for parameter estimation of Weibull distribution | |
CN108121846B (zh) | 一种基于熵惩罚的em未知杂波估计的phd多目标跟踪方法 | |
CN109687845B (zh) | 一种鲁棒的簇稀疏正则化多任务自适应滤波器网络 | |
CN106054167B (zh) | 基于强度滤波器的多扩展目标跟踪方法 | |
CN111798494B (zh) | 广义相关熵准则下的机动目标鲁棒跟踪方法 | |
CN106772354B (zh) | 基于并行模糊高斯和粒子滤波的目标跟踪方法及装置 | |
CN107797106A (zh) | 一种加速em未知杂波估计的phd多目标跟踪平滑滤波方法 | |
Wang et al. | Analysis of binary data via a centered spatial-temporal autologistic regression model | |
CN104050686B (zh) | 一种密集空间目标跟踪方法 | |
Uslu et al. | Finding optimal value for the shrinkage parameter in ridge regression via particle swarm optimization | |
CN113344970A (zh) | 基于多伯努利的非规则多扩展目标联合跟踪与分类方法 | |
CN106226752B (zh) | 一种扩展海杂波Pareto分布参数估计范围的方法 | |
CN111262556A (zh) | 一种同时估计未知高斯测量噪声统计量的多目标跟踪方法 | |
CN106845016B (zh) | 一种基于事件驱动的量测调度方法 | |
CN111722213A (zh) | 一种机动目标运动参数的纯距离提取方法 | |
Shen et al. | RFS‐based extended target multipath tracking algorithm | |
CN106934124B (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 | ||
GR01 | Patent grant |