多重搜索粒子概率假设密度滤波的多目标跟踪方法
【技术领域】
[0001] 本发明涉及一种雷达数据处理方法，特别是涉及一种低检测概率情况下雷达对多 目标的跟踪方法。
【背景技术】
[0002] 粒子概率假设密度滤波（Particle probability hypothesis density filter，PPHDF)是一种密集杂波环境下对多目标进行跟踪的有效方法。PPHDF通过将量测 和目标状态建模为随机集，可以非常方便的从个数时变的量测中估计出个数时变且未知的 目标状态，即可同时对目标个数和目标状态进行估计；同时，PPHDF可以避免目标和量测之 间的关联问题，极大的降低了多目标跟踪算法的复杂性和计算量。因此，PPHDF在多目标跟 踪领域得到了广泛的关注和研究。该方法主要通过以下步骤实现：
[0003] (1)初始化，得到初始粒子集；
[0004] (2)对已存在的粒子集进行一步预测得到预测粒子集，并生成用于搜索新目标的 搜索粒子集，将预测粒子集和搜索粒子集合成新的预测粒子集；
[0005] (3)利用新量测对预测粒子集粒子权重进行更新；
[0006] (4)对权重更新后的粒子集进行重采样，并得到目标个数和各目标状态估计。
[0007] 基于PPHDF的多目标跟踪方法存在一个明显的缺陷，即当目标出现漏检时，重采 样会造成粒子多样性的迅速退化，进而造成目标丢失的现象，因此该算法难以适应目标检 测概率较低时的多目标跟踪。
【发明内容】
[0008] 本发明的目的是提出一种多重搜索粒子概率假设密度滤波（MS-PPHDF)的多目标 跟踪方法，解决一般的PPHDF方法在目标检测概率较低的情况下容易出现目标丢失的问 题。
[0009] 本发明提出的MS-PPHDF方法的技术方案包括以下步骤：
[0010] 步骤1 :变量初始化
[0011] K是雷达关机时刻；
[0012] T是雷达扫描周期；
[0013] L。为代表1个目标的粒子数；
[0014] D。为目标出现的初始分布；
[0015] Jk为搜索新目标的粒子数；
[0016] &表示搜索1个消失目标的粒子数；
[0017] LkS k时刻滤波器采用的粒子总数；
[0018] "为平均目标出现概率；
[0019] PD为目标检测概率；
[0020] λ k为平均每帧的杂波个数；
[0021] Fk为状态转移矩阵；
[0022] GkS过程噪声分布矩阵；
[0023] Qk为过程噪声协方差；
[0024] Rk为量测噪声协方差；
[0025] 步骤2 :令k = 0,初始化粒子集
[0026] 对任意p e {1，2,…，L。}，从初始分布D。中米样粒子4，.并赋予该粒子权重
_.得到初始粒子集
|其中，
表示粒子代表的目标状 态，包含了目标的位置(.<，<)和速度(却,.<)信息；
[0027] 步骤3 :令k = k+Ι，获得k时刻的雷达量测
[0028] 将雷达接收到的信号进行A/D变换，得到k时刻的雷达量测集
，送雷达数据处理计算机，其中
'表示k时刻雷达得 到的第q个量测，包含了目标的距离量测^?、多普勒速度量测4以及方位量测繆等信息， 而Mk则表示k时刻雷达得到的量测个数；
[0029] 步骤4 :预测
[0030] (1)若k彡2,令Ik= 0,直接转（4)，否则定义集合
[0032] 其中，表示k-2时刻存在而k-Ι时刻消失的第η个目标的状态，Ndls, k i表示 消失的目标数；
[0033] (2)若 Xdlsp =0，令 Ik= 0,直接转（4)，否则令 I k= N dls,k Λ，对任意
，根据状态转移方程
[0035] 进行一步预测，其中
[0037] (3)对任意 n
e {1，2,…，Ndis.k和任意 p e {Lk i+Oi-DSk+l，…，Lk i+nSk}，根据
[0039] 采样粒子，并赋予该粒子权重^，^ = )，其中
[0041] 为过程噪声分布矩阵，vk为过程噪声，其噪声协方差为Qk;
[0042] (4)对任意p e {Lk i+Ik+l，…，Lk片以以，根据初始分布D。采样"新生"粒子4-,， 并赋予该粒子权重
[0043] (5)对任意 p e {1，2,…，Lk J，根据
[0045] 采样粒子·^4，并赋予该粒子权重
[0046] 步骤5 :更新
[0047] 对任意p e {1，2,…，Lk flk+JJ，利用量测集Zk对粒子权重进行更新
[0049] 其中
[0052] gk (z | X)为量测似然函数；
[0053] 步骤6 :重采样
[0054] (1)计算所有粒子的权重和
[0056] (2)对粒子集
进行重采样，得到新的粒子集
为k时刻估计的目标数， Round (X)表示 取与X最接近的整数；步骤7 :目标状态估计
[0057] 若/^ = 0，:直接转步骤8,否则采用K-均值聚类分析的方法将粒子集{<1·^划分 为&个类
，Lkin表示第η个类包含的粒子数，满 足
：.，.则估计的目标状态为
[0059] 步骤8 :重复步骤3~步骤8,直至雷达关机。
[0060] 和【背景技术】相比，本发明的有益效果说明：
[0061 ] 本发明提出的MS-PPHDF多目标跟踪方法，可以解决一般的PPHDF方法在目标检测 概率较低的情况下容易出现目标丢失的问题，提高了 PPHDF算法的适应范围。
【附图说明】
[0062] 附图1是本发明提出的MS-PPHDF多目标跟踪方法的整体流程图；
[0063] 附图2是本发明实施例中MS-PPHDF方法跟踪多目标的效果展示，附图中星号 表示目标的真实位置，圆圈" Ο "表示估计的目标状态；
[0064] 附图3是本发明实施例中一般的PPHDF方法跟踪多目标的效果展示，附图中星号 " "表示目标的真实位置，圆圈" 〇 "表示估计的目标状态。
【具体实施方式】
[0065] 下面结合附图对本发明提出的MS-PPHDF方法进行详细描述。
[0066] 不失一般性，设置一个二维的仿真场景，监测区域S =
[_60km，60km] X [_60km，60km]，总仿真时间K = 50s。假设目标可在监测区域内随机出现 和消失，平均目标出现概率γ k= 〇. 2,目标出现的初始分布D。服从均值为X。和协方差为Q b的正态分布，这里取x〇= [30km 0.2km/s 0.5km/s lkm 0. 5km/s])，x方向和y方向过程噪声的标准差均为0. 01km，目标持续存在的概率为ek|k i = 0. 95且与目标状态无关。雷达位于点（0km，-10km)，检测概率为PD= 0. 75,平均每帧的杂 波个数为Ak= 4,且在整个监测区域内均匀分布，距离量测误差标准差、多普勒速度量测误 差标准差和方位角量测误差标准差分别为0. 2km、0. 04km/s和0. 0087rad，雷达采样间隔T =Is ;代表1个目标的粒子数LQ= 3000,搜索新目标的粒子数J k= 4000,搜索1个消失目 标的粒子数sk= 1000。
[0067] 其步骤如附图1所示。
[0068] (1)根据以上仿真条件进行变量初始化
[0069] 由以上仿真条件可知
[0071] 雷达量测方程为
[0073] 其中Wk为相互独立的零均值高斯白噪声，其协方差为
[0075] (2)按
【发明内容】
部分步骤2所述的方法进行粒子集初始化，得到初始粒子集；
[0076] (3)按
【发明内容】
部分步骤3所述的方法获得当前时刻雷达量测；
[0077] (4)按
【发明内容】
部分步骤4所述的方法对粒子集进行一步预测，得到预测的粒子 集；
[0078] (5)按
【发明内容】
部分步骤5所述的方法对预测的粒子集的权重进行更新；
[0079] (6)按
【发明内容】
部分步骤6所述的方法对权重更新后的粒子集进行重采样；
[0080] (7)按
【发明内容】
部分步骤7所述的方法得到目标的状态估计；
[0081] (8)循环执行
【发明内容】
部分步骤3~步骤8,直至雷达关机。
[0082] 实施例条件中，在目标检测概率PD= 0. 75的情况下，本发明提出的MS-PPHDF方法 仍能实现杂波环境下对多目标的有效跟踪（见附图2)，而一般的PPHDF方法在目标出现漏 检时出现了目标的丢失的情况（见附图3)，因此本发明方法可以有效的改善一般PPHDF方 法在目标检测概率较低情况下的多目标跟踪性能（附图2和附图3对比）。
【主权项】
1.多重搜索粒子概率假设密度滤波的多目标跟踪方法，其特征包括以下步骤： 步骤1 :变量初始化 K是雷达关机时刻； T是雷达扫描周期； L。为代表1个目标的粒子数； D。为目标出现的初始分布； Jk为搜索新目标的粒子数； &表示搜索1个消失目标的粒子数； LkS k时刻滤波器采用的粒子总数； YkS平均目标出现概率； Pd为目标检测概率； Xk为平均每帧的杂波个数； Fk为状态转移矩阵； GkS过程噪声分布矩阵； Qk为过程噪声协方差； Rk为量测噪声协方差； 步骤2 :令k = 0,初始化粒子集 对任意P e {1，2,…，L。}，从初始分布D。中米样粒子碟，并赋予该粒子权重得到初始粒子集其中表示粒子代表的目标状 态，包含了目标的位置(<，<)和速度(冗，咒）信息； 步骤3 :令k = k+Ι，获得k时刻的雷达量测 将雷达接收到的信号进行A/D变换，得到k时刻的雷达量测_送雷达数据处理计算机，其中表示k时刻雷达得到的第q个量测，包含了目 标的距离量测<、多普勒速度量测爾以及方位量测頌等信息，而Mk则表示k时刻雷达得到 的量测个数； 步骤4 :预测 (1)若k < 2,令Ik= 0,直接转（4)，否则定义集合其中，：表示k-2时刻存在而k-Ι时刻消失的第η个目标的状态，N dlsik i表示消失 的目标数； ⑵若直接转（4)，否则令I k= Ndls,k ^，对任意根据状态转移方程进行一步预测，其中(3) 对任意 n e {1，2, · · ·，Ndis,k J 和任意 p e {Lk !+(n-DSk+l，···，Lk i+nSk}，根据采样粒子，并赋予该粒子权重其中为过程噪声分布矩阵，vk为过程噪声，其噪声协方差为Q k; (4) 对任意根据初始分布D。采样"新生"粒子并 赋予该粒子权重(5) 对任意p e {1，2,…，Lk J，根据采样粒子并赋予该粒子权重步骤5 :更新 对任意P e {1，2,…，Lk Mk+Jj，利用量测集Zk对粒子权重进行更新gk (z I X)为量测似然函数； 步骤6 :重采样 (1) 计算所有粒子的权重和(2) 对粒子集：:进行重采样，得到新的粒子集其 中为k时刻估计的目标数，表示取与X最接近的整 数；步骤7:目标状态估计 若< =O，直接转步骤8,否则采用K-均值聚类分析的方法将粒子集划分为& 个类表示第η个类包含的粒子数，满足则估计的目标状态为步骤8 :重复步骤3~步骤8,直至雷达关机。
【专利摘要】本发明公开了一种多重搜索粒子概率假设密度滤波的多目标跟踪方法，属于雷达数据处理领域。基于粒子概率假设密度滤波的多目标跟踪方法存在一个明显的缺陷，即当目标出现漏检时，重采样会造成粒子多样性的迅速退化，进而造成目标丢失的现象，因此该算法难以适应目标检测概率较低时的多目标跟踪。本发明提出的多重搜索粒子概率假设密度滤波即立足于解决此类问题。本发明具有结构简单，计算快速，易于硬件实现，同时克服了基于一般的粒子概率假设密度滤波方法应用的局限性，对非线性非高斯系统具有较强的适应性，因此具有较强的工程应用价值和推广前景。
【IPC分类】G01S13/66, G01S7/02, G01S13/56
【公开号】CN105353353
【申请号】CN201510791334
【发明人】谭顺成, 王国宏, 吴巍, 于洪波
【申请人】中国人民解放军海军航空工程学院
【公开日】2016年2月24日
【申请日】2015年11月17日