【技术实现步骤摘要】
本专利技术属于信息处理
,涉及一种多扩展目标跟踪方法,具体涉及一种基于箱粒子CPHD的多扩展目标跟踪方法,可用于车辆编队和人群等系统中的跟踪。
技术介绍
随着现代先进传感器如相控阵雷达和逆合成孔径雷达等探测设备分辨率的不断提高,每个运动目标的回波信号可能会同时占据多个距离分辨单元,即每个目标可能产生多个量测,称这样的目标为扩展目标。扩展目标跟踪不仅可以提供目标的运动状态,包括位置、速度、加速度,还可以提供目标的扩展形态,包括目标大小、形状等多方面信息,因此引起了国内外学者的广泛关注。在现实跟踪场景中,存在杂波和虚警,往往无法预知扩展目标的个数,对数目未知且变化的多扩展目标跟踪已成为目标跟踪领域具有挑战的热门课题。多扩展目标跟踪方法可分为基于数据关联的多扩展目标跟踪方法和基于随机集的多扩展目标跟踪方法两大类。在基于数据关联的多扩展目标跟踪方法中,由于多扩展目标会产生大量量测,每个量测与每个扩展目标进行关联假设,会带来巨大的数据计算量,导致跟踪实时性较差。基于随机集的多扩展目标跟踪方法,将目标状态集合和量测集合均建模成随机集,跟踪过程中避免了复杂的数据关联,具有良好的实时性。Mahler结合扩展目标空间分布模型和随机有限集理论提出的扩展目标PHD滤波算法有效地避免了复杂的数据关联,对多扩展目标进行有效跟踪;随后,Orguner等又提出了多扩展目标的CPHD滤波算法。箱粒子(Box-Particle)滤波是由Abdallah等提出的一种基于区间分析的蒙特卡罗滤波算法,能够有效处理量测的三种不确定性。在实际应用中,箱粒子滤波只需要用几十个的箱粒子代替即可达 ...
【技术保护点】
一种基于箱粒子CPHD的多扩展目标跟踪方法,包括如下步骤:(1)将k时刻的多扩展目标量测集合Zk划分成多个划分单元并利用包含函数[f]分别对多个划分单元进行包含,得到k时刻的量测箱粒子集其中k表示滤波时刻,mk表示k时刻的划分单元个数,Wj表示第j个划分单元内量测分布的空间范围,表示k时刻的量测箱粒子状态;(2)假设k‑1时刻的存活箱粒子集为初始势分布为pk‑1(n),将k‑1时刻的新生箱粒子集并入存活箱粒子集得到k‑1时刻的初始箱粒子集其中和分别表示k‑1时刻的存活箱粒子状态、新生箱粒子状态和初始箱粒子状态,和分别表示k‑1时刻的存活箱粒子强度、新生箱粒子强度和初始箱粒子强度,Nk‑1,s、Nk‑1,b和Nk‑1分别表示k‑1时刻的存活箱粒子数目、新生箱粒子数目和初始箱粒子数目,n表示目标个数;(3)对得到的k‑1时刻的初始箱粒子集进行预测,得到k时刻的预测箱粒子集同时对假设的k‑1时刻的初始势分布pk‑1(n)进行势分布预测,得到k时刻的预测势分布pk|k‑1(n),其中,和分别表示k时刻的预测箱粒子状态和预测箱粒子强度;(4)对k时刻的预测箱粒子集和预测势分布pk|k‑1(n)进 ...
【技术特征摘要】
1.一种基于箱粒子CPHD的多扩展目标跟踪方法,包括如下步骤:(1)将k时刻的多扩展目标量测集合Zk划分成多个划分单元并利用包含函数[f]分别对多个划分单元进行包含,得到k时刻的量测箱粒子集其中k表示滤波时刻,mk表示k时刻的划分单元个数,Wj表示第j个划分单元内量测分布的空间范围,表示k时刻的量测箱粒子状态;(2)假设k-1时刻的存活箱粒子集为初始势分布为pk-1(n),将k-1时刻的新生箱粒子集并入存活箱粒子集得到k-1时刻的初始箱粒子集其中和分别表示k-1时刻的存活箱粒子状态、新生箱粒子状态和初始箱粒子状态,和分别表示k-1时刻的存活箱粒子强度、新生箱粒子强度和初始箱粒子强度,Nk-1,s、Nk-1,b和Nk-1分别表示k-1时刻的存活箱粒子数目、新生箱粒子数目和初始箱粒子数目,n表示目标个数;(3)对得到的k-1时刻的初始箱粒子集进行预测,得到k时刻的预测箱粒子集同时对假设的k-1时刻的初始势分布pk-1(n)进行势分布预测,得到k时刻的预测势分布pk|k-1(n),其中,和分别表示k时刻的预测箱粒子状态和预测箱粒子强度;(4)对k时刻的预测箱粒子集和预测势分布pk|k-1(n)进行量测更新,实现步骤如下:(4a)根据得到的k时刻的量测箱粒子集和k时刻的预测箱粒子集求得似然函数(4b)利用求得的似然函数对得到的k时刻的预测箱粒子集进行量测更新,得到k时刻的更新箱粒子集(4c)利用求得的似然函数对得到的k时刻的预测势分布pk|k-1(n)进行量测更新,得到k时刻的势分布pk(n);(5)对得到的k时刻的量测箱粒子集进行搜索,并判断是否搜索到对预测箱粒子强度贡献最大的量测箱[zk],若是,执行步骤(6),否则,将得到的k时刻的更新箱粒子集输出;(6)利用搜索到的量测箱[zk],对得到的k时刻的更新箱粒子集进行收缩,得到k时刻的收缩箱粒子集其中表示收缩箱粒子状态;(7)对得到的k时刻的势分布pk(n)进行估计,得到目标数目估计(8)利用随机子划分方法,对得到的k时刻的收缩箱粒子集进行重采样,得到k时刻的重采样箱粒子集为其中和分别表示k时刻的重采样箱粒子状态和重采样箱粒子强度,N表示重采样次数;(9)对得到的k时刻的重采样箱粒子状态集进行点化,得到k时刻的箱中心粒子状态集对该箱中心粒子状态集进行k-means聚类,得到目标的状态估计其中表示k时刻的箱中心粒子状态;(10)判断k时刻是否为跟踪最后时刻,若是,则跟踪结束;否则保存k时刻的重采样箱粒子集并作为下个时刻的存活箱粒子集,令k=k+1,执行步骤(2)。2.根据权利要求1中所述的基于箱粒子CPHD的多扩展目标跟踪方法,其特征在于:步骤(2)中所述的k-1时刻新生箱粒子集是通过对上一时刻的量测箱粒子周围采样产生的,其状态服从高斯分布权值设置为其中,Pb表示箱粒子新生概率,Σ表示新生箱粒子状态分布协方差。3.根据权利要求1中所述的基于箱粒子CPHD的多扩展目标跟踪方法,其特征在于:步骤(3)中所述的对得到的k-1时刻的初始箱粒子集进行预测,是通过系统预测模型和强度预测方程实现的,其中,fk|k-1表示状态转移函数,PS表示箱粒子存活概率。4.根据权利要求1中所述的基于箱粒子CPHD的多扩展目标跟踪方法,其特征在于:步骤(3)中所述的对假设的k-1时刻的初始势分布pk-1(n)进行势分布预测,是通过马尔科夫转移矩阵和势分布预测方程实现的,其中,pbirth(n)表示从k-1时刻到k时刻有n个目标新生的概率,d表示目标消失概率。5.根据权利要求1中所述的基于箱粒子CPHD的多扩展目标跟踪方法,其特征在于:步骤(4a)中所述的似然函数采用以下公式求得: L k ( [ zp k ( j ) ] | [ x k | k - 1 ( i ) ] ) = | [ h c p ] ( [ x k | k - 1 ...
【专利技术属性】
技术研发人员:宋骊平,宋志龙,姬红兵,苗雨,程轩,
申请(专利权)人:西安电子科技大学,
类型:发明
国别省市:陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。