法律状态公告日
法律状态信息
法律状态
2019-07-19
专利权的转移 IPC(主分类):G01S13/72 登记生效日:20190701 变更前: 变更后: 申请日:20151207
专利申请权、专利权的转移
2017-12-05
授权
授权
2016-03-02
实质审查的生效 IPC(主分类):G01S13/72 申请日:20151207
实质审查的生效
2016-02-03
公开
公开
技术领域
本发明涉及一种雷达数据处理方法,特别是涉及一种雷达测距模糊情况下的机动多目标跟踪方法,适用于高/中脉冲重复频率(H/MPRF)雷达(如机载脉冲多普勒雷达等)对机动多目标的跟踪。
背景技术
随着电磁干扰等技术的迅猛发展和各类飞行器的大量出现,雷达所处的电磁环境进一步恶化,单位时间内需要处理的量测数据急剧,给传统的雷达目标检测和跟踪带来了严峻的挑战。同时,为了提高自身的生存和突防能力,目标通常会随时进行一定的战术机动,其运动状态也会因不同的作战任务而变,进一步增大了雷达对其进行跟踪的困难。因此,如何实现复杂环境下雷达对机动多目标的有效检测和跟踪,对提高雷达作战效能和生存能力具有重要意义,是目前目标跟踪领域研究的重点和难点问题。目前基于随机有限集的粒子概率假设密度滤波(PPHDF)方法由于具有可以避免多目标数据关联的突出优点,极大的降低多目标跟踪算法的计算复杂度,适应于密集杂波环境下目标数目未知的多目标跟踪场景,在目标跟踪领域得到广泛的关注和研究。该方法主要通过以下步骤实现:
(1)粒子集初始化;
(2)生成用于搜索新目标的搜索粒子集,并对已存在的粒子集进行一步预测得到预测粒子集,再将搜索粒子集和预测粒子集合成新的预测粒子集;
(3)利用新量测对预测粒子集粒子权重进行更新;
(4)对权重更新后的粒子集进行重采样,并得到目标个数和各目标状态估计。
基于PPHDF的多目标跟踪方法存在以下两个缺陷:(1)当目标出现漏检时,重采样会造成粒子多样性的迅速退化,进而造成目标丢失的现象;(2)PD雷达为了无模糊地测量目标速度,通常采用H/MPRF的工作模式,在这种情况下雷达对目标的距离量测是模糊的,PPHDF算法无法直接利用模糊的量测数据对目标进行正确稳定的跟踪。
发明内容
本发明的目的是提出一种同时解距离模糊的改进粒子概率假设密度滤波(IPPHDF)机动多目标跟踪方法,解决一般的PPHDF方法在检测概率较低情况下容易出现目标丢失和无法直接利用模糊量测数据对目标进行跟踪的问题。
本发明提出的同时解距离模糊的IPPHDF机动多目标跟踪方法的技术方案包括以下步骤:步骤1:变量初始化
(1)T是雷达扫描周期,T1,T2,...,TM是脉冲重复周期,R1,R2,...,RM是各PRF对应的最大不模糊距离;
(2)L0为代表1个目标的粒子数,Jk为搜索新目标的粒子数,Sk表示搜索1个消失目标的粒子数,Lk为k时刻滤波器采用的粒子总数;
(3)D0为目标出现的初始分布,γk为平均目标出现概率,PD为目标检测概率;
(4)λk为平均每帧的杂波个数,
(5)Gk为过程噪声分布矩阵,Qk为过程噪声协方差,Rk为量测噪声协方差;
步骤2:初始化粒子集,令k=0和
(1)从初始分布D0中采样粒子
(2)赋予粒子
步骤3:令k=k+1,获得k时刻的雷达模糊量测
(1)设置雷达的工作状态,使雷达依次交替采用某个PRF工作,令
c=mod(k,M)+1
表示雷达在k时刻采用的PRF的索引号,mod(x,y)表示x/y的余数;
(2)将雷达接收到的信号进行A/D变换,得到k时刻的雷达模糊量测集>送雷达数据处理计算机,其中>表示k时刻雷达得到的第q个量测,包含了目标的模糊距离量测
步骤4:已存在粒子集预测,若
(1)根据目标角速度计算目标状态状态转移矩阵
>
(2)利用粒子
>
得到粒子
>
为过程噪声分布矩阵,vk为过程噪声,其噪声协方差为Qk;
(3)赋予粒子
步骤5:生成搜索新目标的粒子集,对任意p∈{Lk-1+Ik+1,…,Lk-1+Ik+Jk}
(1)根据初始分布D0采样“新生”粒子
(2)赋予粒子
步骤6:生成消失目标预测粒子集
(1)若k≤2,令Ik=0,直接转步骤7,否则定义集合
>
其中
(2)若
>
(3)对任意n∈{1,2,…,Ndis,k-1}和任意p∈{Lk-1+(n-1)Sk+1,…,Lk-1+nSk},令>利用
>
得到粒子
步骤7:粒子权重更新
(1)对任意p∈{1,2,…,Lk-1+Ik+Jk},根据量测方程
>
得到预测量测
>
得到模糊预测量测
(2)对任意p∈{1,2,…,Lk-1+Ik+Jk}和任意m∈{1,2,…,Mk},计算新息
>
并计算
>
其中Rk为量测噪声协方差;
(3)对任意m∈{1,2,…,Mk},计算
>
(4)对任意p∈{1,2,…,Lk-1+Ik+Jk},计算粒子权重
>
步骤8:重采样
(1)计算所有粒子的权重和
>
(2)对粒子集
步骤9:目标状态估计
(1)若
(2)若
(3)估计目标的PIN
>
其中Floor(x)表示取小于或等于x最大整数;
步骤10:重复步骤3~步骤9,直至雷达关机。
和背景技术相比,本发明的有益效果说明:
本发明提出的同时解距离模糊的IPPHDF机动多目标跟踪方法,解决一般的PPHDF方法在目标检测概率较低的情况下容易出现目标丢失的问题,且该方法可以直接利用雷达得到的模糊量测数据对机动多目标进行跟踪,在跟踪的同时解决距离量测模糊的问题,扩大了PPHDF算法的适应范围。
附图说明
附图1给出了本发明提出的同时解距离模糊的IPPHDF机动多目标跟踪方法整体流程图;
附图2给出了本发明实施例中三个机动目标“时刻-x位置”和“时刻-y位置”分开表示的真实轨迹示意图,其中:
附图3是给出了本发明实施例中雷达得到的各时刻模糊量测示意图,包含了杂波和雷达对目标的模糊量测;
附图4是本发明实施例中同时解距离模糊的IPPHDF跟踪机动多目标的效果展示,附图中星号“*”表示目标的真实位置,圆圈“o”表示估计的目标状态;
具体实施方式
下面结合附图对本发明提出的同时解距离模糊的IPPHDF方法进行详细描述。
不失一般性,设置一个二维的仿真场景,监测区域S=[-60km,60km]×[-60km,60km],目标可在监测区域内随机出现和消失,平均目标出现概率γk=0.2,目标出现的初始分布D0服从均值为x0和协方差为Qb的正态分布,其中
>
x方向和y方向过程噪声的标准差均为0.01km,转弯率的过程噪声标准差0.35rad/s,目标持续存在的概率为ek|k-1=0.95且与目标状态无关。雷达扫描周期T=1s,可交替采用2个不同的脉冲重复频率交替工作,各脉冲重复周期分别为T1=37μs和T2=48μs,对应的最大不模糊距离分别为R1=5.55km和R2=7.20km,雷达位于点(0km,-10km),检测概率为PD=0.75,平均每帧的杂波个数为λk=4,且在整个监测区域内均匀分布,距离量测误差标准差、多普勒速度量测误差标准差和方位角量测误差标准差分别为0.2km、0.04km/s和0.0087rad/s;代表1个目标的粒子数L0=3000,搜索新目标的粒子数Jk=4000,搜索1个消失目标的粒子数Sk=1000。
其步骤如附图1所示。
(1)根据以上仿真条件进行变量初始化
(2)按发明内容部分步骤2所述的方法进行粒子集初始化;
(3)按发明内容部分步骤3所述的方法获得当前时刻雷达模糊量测集;
(4)按发明内容部分步骤4所述的方法对已存在粒子集进行一步预测;
(5)按发明内容部分步骤5所述的方法生成搜索新目标的粒子集;
(6)按发明内容部分步骤6所述的方法生成消失目标预测粒子集;
(7)按发明内容部分步骤7所述的方法对预测粒子集进行权重更新;
(8)按发明内容部分步骤8所述的方法对权重更新后的粒子集进行重采样;
(9)按发明内容部分步骤9所述的方法得到目标的状态和PIN估计;
(10)循环执行发明内容部分步骤3~步骤9,直至雷达关机。
实施例条件中,且雷达得到的目标距离量测是模糊的(见附图3),本发明提出的同时解距离模糊的IPPHDF方法仍能实现杂波环境下对机动多目标的有效跟踪(见附图4),因此本发明方法克服了一般的PPHDF方法在检测概率较低情况下容易出现目标丢失和无法直接利用模糊量测数据对目标进行跟踪的问题。
机译: 多目标行人跟踪方法,多目标行人跟踪装置和多目标行人跟踪装置
机译: 多目标跟踪系统,多目标跟踪方法和多目标跟踪程序。
机译: 多目标跟踪设备,多目标跟踪方法和多目标跟踪程序