改进搜索策略的MM‑PPHDF机动多目标跟踪方法技术

本发明专利技术公开了一种改进搜索策略的MM‑PPHDF机动多目标跟踪方法，属于雷达数据处理领域。基于MM‑PPHDF的机动多目标跟踪方法存在一个明显的缺陷，即当目标出现漏检时，重采样会造成粒子多样性的迅速退化，进而造成目标丢失的现象，因此该算法难以适应目标检测概率较低时的机动多目标跟踪。本发明专利技术提出的改进搜索策略的MM‑PPHDF机动多目标跟踪方法立足于解决此类问题。本发明专利技术具有结构简单，计算快速，对强非线性非高斯系统具有良好的适应性等优点，实施例的实验仿真结果也表明本发明专利技术可以有效克服基于MM‑PPHDF方法应用的局限性，因此本发明专利技术公开的改进搜索策略的MM‑PPHDF机动多目标跟踪方法具有较强的工程应用价值和推广前景。

【技术实现步骤摘要】
改进搜索策略的MM-PPHDF机动多目标跟踪方法
本专利技术涉及一种雷达数据处理方法，特别是涉及一种低检测概率情况下雷达对机动多目标的跟踪方法。
技术介绍
电磁环境的日益复杂和各类飞行器的大量出现，给雷达的检测跟踪性能带来了严峻的挑战。同时，为了提高生存和突防能力，目标通常会随时进行一定的机动，其运动状态也会因不同的作战任务而变，进一步增大了雷达对其进行跟踪的困难。因此，实现复杂电磁环境下雷达对机动多目标的有效跟踪，对提高雷达作战效能和生存能力具有重要意义。目前基于粒子概率假设密度滤波(PPHDF)的多目标跟踪方法通过将量测和目标状态建模为随机集，可以在避免目标和量测之间的关联问题的同时对目标个数和目标状态进行估计，极大的降低了多目标跟踪算法的复杂性和计算量，因此在多目标跟踪领域得到广泛的关注，而PPHDF通过借鉴多模型粒子滤波(MMPF)的思想，形成了一种基于多模型粒子概率假设密度滤波(MM-PPHDF)的机动多目标跟踪方法。该方法主要通过以下步骤实现：(1)将雷达接收机输出的目标回波信号进行A/D变换，送雷达数据处理计算机执行以下步骤；(2)对目标运动模式进行预测，并利用预测的目标运动模式和相应的目标状态转移方程对已存在的粒子集进行预测；(3)生成用于搜索新目标的搜索粒子集，将预测粒子集和搜索粒子集合成新的预测粒子集；(4)利用新量测对预测粒子集粒子权重进行更新；(5)对权重更新后的粒子集进行重采样，并得到目标个数和各目标状态估计。基于MM-PPHDF的机动多目标跟踪方法存在一个比较明显的缺陷，即当目标出现漏检时，重采样会造成粒子多样性的迅速退化，进而造成目标丢失的现象。因此，该算法难以适应目标检测概率较低时的机动多目标跟踪。
技术实现思路
本专利技术的目的是提出一种改进搜索策略的MM-PPHDF机动多目标跟踪方法，解决一般的MM-PPHDF方法在检测概率较低情况下容易出现目标丢失的问题。本专利技术提出的改进搜索策略的MM-PPHDF机动多目标跟踪方法的技术方案包括以下步骤：步骤1：变量初始化K是雷达关机时刻；T是雷达扫描周期；L0为代表1个目标的粒子数；D0为目标出现的初始分布；Jk为搜索新目标的粒子数；Sk表示搜索1个消失目标的粒子数；Lk为k时刻滤波器采用的粒子总数；γk为平均目标出现概率；PD为目标检测概率；λk为平均每帧的杂波个数；为k时刻估计的目标个数；是目标初始运动模式概率；∏m是目标运动模式概率转移矩阵；为状态转移矩阵；Gk为过程噪声分布矩阵；Qk为过程噪声协方差；Rk为量测噪声协方差；步骤2：初始化粒子集，令k＝0和对任意p∈{1,2,…,L0}(1)根据目标初始运动模式概率初始化目标运动模式变量(2)从初始分布D0中采样粒子其中表示粒子代表的目标状态，包含了目标的位置和速度信息；(3)赋予粒子权重步骤3：令k＝k+1，获得k时刻的雷达量测将雷达接收到的信号进行A/D变换，得到k时刻的雷达量测集送雷达数据处理计算机，其中表示k时刻雷达得到的第q个量测，包含了目标的距离量测多普勒速度量测以及方位量测等信息，而Mk则表示k时刻雷达得到的量测个数；步骤4：已存在粒子集预测，若直接转步骤5，否则对任意p∈{1,2,…,Lk-1}(1)根据上一时刻的目标运动模式变量和目标运动模式概率转移矩阵∏m预测当前时刻的目标运动模式变量(2)根据选择目标状态状态转移矩阵，其中对应于目标匀速直线运动对应于目标坐标转动，其中分别对应于顺时针坐标转动和逆时针坐标转动，αm>0是一个典型的机动加速度；(3)利用粒子的状态进行一步预测得到粒子其中为过程噪声分布矩阵，vk为过程噪声，其噪声协方差为Qk；(4)赋予粒子权重步骤5：生成搜索新目标的粒子集，对任意p∈{Lk-1+Ik+1,…,Lk-1+Ik+Jk}(1)根据目标初始运动模式概率初始化目标运动模式变量(2)根据初始分布D0采样“新生”粒子(3)赋予粒子权重步骤6：生成消失目标预测粒子集(1)若k≤2，令Ik＝0，直接转步骤7，否则定义集合其中，表示k-2时刻存在而k-1时刻消失的第n个目标的状态，Ndis,k-1表示消失的目标数；(2)若令Ik＝0，直接转步骤7，否则令Ik＝Ndis,k-1Sk，对任意根据估计的目标运动模式和相应的状态转移方程进行一步预测(3)对任意n∈{1,2,…,Ndis,k-1}和任意p∈{Lk-1+(n-1)Sk+1,…,Lk-1+nSk}，令利用进行一步预测得到粒子并赋予该粒子权重步骤7：更新对任意p∈{1,2,…,Lk-1+Ik+Jk}，利用量测集Zk对粒子权重进行更新其中gk(z|x)为量测似然函数；步骤8：重采样(1)计算所有粒子的权重和(2)对粒子集进行重采样，得到新的粒子集其中Round(x)表示取与x最接近的整数；步骤9：目标状态和运动模式估计(1)若直接转步骤10，否则采用聚类分析的方法将粒子集划分为个类其中Lk,n表示第n个类包含的粒子数，满足(2)估计目标状态(3)估计目标的运动模式其中函数表示集合中元素m的个数；步骤10：重复步骤3～步骤8，直至雷达关机。和
技术介绍
相比，本专利技术的有益效果说明：本专利技术提出的改进搜索策略的MM-PPHDF多目标跟踪方法，可以解决一般的MM-PPHDF方法在目标检测概率较低的情况下容易出现目标丢失的问题，提高了MM-PPHDF算法的适应范围。附图说明附图1是本专利技术提出的改进搜索策略的MM-PPHDF机动多目标跟踪方法的整体流程图；附图2是本专利技术实施例中三个机动目标的真实轨迹示意图；附图3是本专利技术实施例中改进搜索策略的MM-PPHDF方法跟踪机动多目标的效果展示，附图中星号“*”表示目标的真实位置，圆圈“o”表示估计的目标状态；附图4是本专利技术实施例中一般的MM-PPHDF方法跟踪机动多目标的效果展示，附图中星号“*”表示目标的真实位置，圆圈“o”表示估计的目标状态。具体实施方式下面结合附图对本专利技术提出的改进搜索策略的MM-PPHDF方法进行详细描述。不失一般性，设置一个二维的仿真场景，监测区域S＝[-60km,60km]×[-60km,60km]，总仿真时间K＝50s。假设目标都服从以下3种运动模型中的一种：匀速直线运动模式，顺时针坐标转动模式和逆时针坐标转动模式，目标初始运动模式概率可在3种运动模式之间自由转换，各种运动模型之间的转换可描述为3态的1阶Markov链，模式概率转移矩阵机动加速度αm＝1.08。目标可在监测区域内随机出现和消失，平均目标出现概率γk＝0.2，目标出现的初始分布D0服从均值为x0和协方差为Qb的正态分布，这里取x0＝[30km0.2km/s30km-0.1km/s]T和Qb＝diag([1km0.5km/s1km0.5km/s])，x方向和y方向过程噪声的标准差均为0.01km，目标持续存在的概率为ek|k-1＝0.95且与目标状态无关。雷达位于点(0km,-10km)，检测概率为PD＝0.75，平均每帧的杂波个数为λk＝4，且在整个监测区域内均匀分布，距离量测误差标准差、多普勒速度量测误差标准差和方位角量测误差标准差分别为0.2km、0.04km/s和0.0087rad，雷达采样间隔T＝1s；代表1个目标的粒子数L0＝3000，搜索新目标的粒子数本文档来自技高网...

【技术保护点】
改进搜索策略的MM‑PPHDF机动多目标跟踪方法，其特征包括以下步骤：步骤1：变量初始化K是雷达关机时刻；T是雷达扫描周期；L0为代表1个目标的粒子数；D0为目标出现的初始分布；Jk为搜索新目标的粒子数；Sk表示搜索1个消失目标的粒子数；Lk为k时刻滤波器采用的粒子总数；γk为平均目标出现概率；PD为目标检测概率；λk为平均每帧的杂波个数；为k时刻估计的目标个数；是目标初始运动模式概率；∏m是目标运动模式概率转移矩阵；为状态转移矩阵；Gk为过程噪声分布矩阵；Qk为过程噪声协方差；Rk为量测噪声协方差；步骤2：初始化粒子集，令k＝0和对任意p∈{1,2,…,L0}(1)根据目标初始运动模式概率初始化目标运动模式变量(2)从初始分布D0中采样粒子其中表示粒子代表的目标状态，包含了目标的位置和速度信息；(3)赋予粒子权重步骤3：令k＝k+1，获得k时刻的雷达量测将雷达接收到的信号进行A/D变换，得到k时刻的雷达量测集送雷达数据处理计算机，其中表示k时刻雷达得到的第q个量测，包含了目标的距离量测多普勒速度量测以及方位量测等信息，而Mk则表示k时刻雷达得到的量测个数；步骤4：已存在粒子集预测，若直接转步骤5，否则对任意p∈{1,2,…,Lk‑1}(1)根据上一时刻的目标运动模式变量和目标运动模式概率转移矩阵∏m预测当前时刻的目标运动模式变量(2)根据选择目标状态状态转移矩阵，其中Fk(1)=1T00001000001T00001000001]]>对应于目标匀速直线运动Fk(j)=1sin(Ωk(j)T)Ωk(j)0-1-cos(Ωk(j)T)Ωk(j)0cos(Ωk(j)T)0-sin(Ωk(j)T)01-cos(Ωk(j)T)Ωk(j)1sin(Ωk(j)T)Ωk(j)0sin(Ωk(j)T)0cos(Ωk(j)T),j=2,3]]>对应于目标坐标转动，其中Ωk(2)=-αm(x·k-1p)2+(y·k-1p)2,Ωk(3)=αm(x·k-1p)2+(y·k-1p)2]]>分别对应于顺时针坐标转动和逆时针坐标转动，αm>0是一个典型的机动加速度；(3)利用粒子的状态进行一步预测xk,k-1p=Fk-1(mk|k-1p)xk-1p+Gkvk]]>得到粒子其中Gk=T2/20T00T2/20T]]>为过程噪声分布矩阵，vk为过程噪声，其噪声协方差为Qk；(4)赋予粒子权重wk|k-1p=wk-1p;]]>步骤5：生成搜索新目标的粒子集，对任意p∈{Lk‑1+Ik+1,…,Lk‑1+Ik+Jk}(1)根据目标初始运动模式概率初始化目标运动模式变量(2)根据初始分布D0采样“新生”粒子(3)赋予粒子权重wk|k-1p=γk/(Jk+Ik);]]>步骤6：生成消失目标预测粒子集(1)若k≤2，令Ik＝0，直接转步骤7，否则定义集合Xdis,k-2=Δ{xdis,k-2n|n=1,2,...,Ndis,k-1}]]>其中，表示k‑2时刻存在而k‑1时刻消失的第n个目标的状态，Ndis,k‑1表示消失的目标数；(2)若令Ik＝0，直接转步骤7，否则令Ik＝Ndisk‑1Sk，对任意根据估计的目标运动模式和相应的状态转移方程进行一步预测xpre,k-1n=Fk-1(m^k-2,n)xdis,k-2n]]>(3)对任意n∈{1,2,…,Ndis,k‑1}和任意p∈{Lk‑1+(n‑1)Sk+1,…,Lk‑1+nSk}，令mk|k-1p=m^k-2,n,]]>利用进行一步预测xk|k-1p=Fk(mk|k-1p)xpre,k-1n+Gkvk]]>得到粒子并赋予该粒子权重步骤7：更新对任意p∈{1,2,…,Lk‑1+Ik+Jk}，利用量测集Zk对粒子权重进行更新wk|kp=[1-PD+Σz∈Zkψk,z(xk|k-1p)λk+Ck(z)]wk|k-1p]]>其中ψk,z(xk|k-1p)=PDgk(z|xk|k-1p)]]>Ck(z)=Σp=1Lk-1+Ik+Jkψk,z(xk|k-1p)wk|k-1p]]>gk(z|x)为量测似然函数；步骤8：重采样(1)计算所有粒子的权重和N^k|k=Σp=1Lk-1+Ik+Jkwk|kp]]>(2)对粒子集进行重采样，得到新的粒子集{xkp,mkp,wkp/N^k|k}p=1Lk,]]>其中N^k=Round(N...

【技术特征摘要】
1.改进搜索策略的MM-PPHDF机动多目标跟踪方法，其特征在于，包括以下步骤：步骤1，变量初始化：K是雷达关机时刻；T是雷达扫描周期；L0为代表1个目标的粒子数；D0为目标出现的初始分布；Jk为搜索新目标的粒子数；Sk表示搜索1个消失目标的粒子数；Ik为搜索消失目标的粒子数，初始值为零；Lk为k时刻滤波器采用的粒子总数；γk为平均目标出现概率；PD为目标检测概率；λk为平均每帧的杂波个数；为k时刻估计的目标个数；是目标初始运动模式概率；∏m是目标运动模式概率转移矩阵；为状态转移矩阵；Gk为过程噪声分布矩阵；Qk为过程噪声协方差；Rk为量测噪声协方差；步骤2，初始化粒子集，令k＝0和对任意p∈{1,2,…,L0}：(1)根据目标初始运动模式概率初始化目标运动模式变量(2)从初始分布D0中采样粒子其中表示粒子代表的目标状态，包含了目标的位置和速度信息；(3)赋予粒子权重步骤3，令k＝k+1，获得k时刻的雷达量测：将雷达接收到的信号进行A/D变换，得到k时刻的雷达量测集送雷达数据处理计算机，其中表示k时刻雷达得到的第q个量测，包含了目标的距离量测多普勒速度量测以及方位量测信息，而Mk则表示k时刻雷达得到的量测个数；步骤4，已存在粒子集预测，若直接转步骤5，否则对任意p∈{1,2,…,Lk-1}：(1)根据上一时刻的目标运动模式变量和目标运动模式概率转移矩阵Πm预测当前时刻的目标运动模式变量(2)根据选择目标状态状态转移矩阵，其中对应于目标匀速直线运动对应于目标坐标转动，其中分别对应于顺时针坐标转动和逆时针坐标转动，αm＞0是一个典型的机动加速度；(3)利用粒子的状态进行一步预测得到粒子其中...

【专利技术属性】
技术研发人员：谭顺成，王国宏，吴巍，于洪波，
申请(专利权)人：中国人民解放军海军航空工程学院，
类型：发明
国别省市：山东;37