本发明涉及一种雷达数据处理方法,特别是涉及一种多目标跟踪滤波方法,适应于密集杂波环境下雷达对多目标的跟踪。
背景技术:
随着现代科技的迅猛发展,空间技术的日新月异,目标监控范围已经发展到陆、海、空、天、电磁五维空间,同时目标所处的电磁环境日益复杂,目标监控的种类和数量也不断增多,给雷达的探测跟踪性能带来了严峻的挑战。由于噪声、杂波以及数据关联和检测等的不确定性,对多目标进行跟踪需要从个数时变的量测中估计出目标个数不确定的各目标状态,这就要求算法能够同时对目标个数和状态进行估计,因此密集杂波环境下的多目标跟踪已成为目标跟踪领域的难点问题之一。如何综合各种量测信息,实现密集杂波环境下雷达对多目标的有效跟踪,对提高雷达的探测跟踪性能具有重要意义。传统的多目标跟踪算法主要有联合概率数据互联(jpda)和多假设跟踪(mht)等,其中基于jpda的多目标跟踪方法由于具有跟踪精度较高的优点,在多目标跟踪领域受到广泛研究。该方法主要通过以下步骤实现:
(1)量测获取;
(2)确认矩阵的拆分;
(3)联合事件概率的计算;
(4)状态估计和协方差的计算。
基于jpda的多目标跟踪方法具有以下三个缺陷:(1)没有利用目标的幅值信息,造成有用量测信息的丢失;(2)要求目标个数固定,受到实际应用的限制;(3)计算量随量测个数的增加而爆炸增长,从而不能适应密集杂波环境。
技术实现要素:
本发明的目的是提出一种结合幅值信息的phd滤波多目标跟踪方法,解决基于jpda的多目标跟踪方法没有充分利用量测信息,要求目标个数固定,以及不适应密集杂波环境等问题。
本发明提出的结合幅值信息的phd滤波多目标跟踪方法的技术方案包括以下步骤:
步骤1:变量初始化
(1)k为总仿真时间,t表示雷达采样间隔;
(2)γ0为目标出现初始数,l0为代表一个目标所需要的粒子数,jk为k时刻搜索新目标的粒子数,lk为k时刻粒子总数;
(3)
(4)γk(x)为目标出现初始分布,κk(z)为杂波分布;
(5)
(6)
步骤2:令k=0,进行滤波器初始化
(1)对任意i∈{1,2,…,l0},根据目标出现初始分布γ0(x)生成粒子
(2)对任意i∈{1,2,…,l0},令
步骤3:令k=k+1,获得k时刻的雷达量测
将雷达接收到的信号进行a/d变换,得到k时刻的雷达扩展量测集
步骤4:粒子集预测
(1)对任意i∈{1,2,…,lk-1},根据k-1时刻粒子状态
并赋予该粒子权重
(2)对任意i∈{lk-1+1,lk-1+2,…,lk-1+jk},根据目标出现初始分布γk(x)采样粒子
(3)将粒子集
步骤5:粒子集更新
(1)对任意i∈{1,2,…,lk-1+jk},根据预测的粒子状态
(2)对任意i∈{1,2,…,lk-1+jk},利用k时刻的扩展量测集
其中
以及
其中rk为量测误差协方差;
(2)计算所有粒子的权重和
并取与
(3)令
(1)对粒子集
(2)对任意
得到第j个目标的状态估计,其中
步骤7:重复步骤3~步骤6,直至雷达关机。
和背景技术相比,本发明提出的结合幅值信息的phd滤波多目标跟踪方法的有益效果说明:
(1)本发明提出的方法有效利用了目标的幅值,解决了jpda方法没有充分利用量测信息的问题,有效提高了多目标的跟踪性能;(2)本发明提出的方法适应于目标个数时变且未知的情况,解决了jpda方法要求目标数固定的问题;(3)本发明提出的方法计算量随量测个数的增加线性增长,解决jpda方法计算量随量测个数的增加爆炸增长的问题,尤其适应于密集杂波环境下的多目标跟踪。
附图说明
附图1是本发明的结合幅值信息的phd滤波多目标跟踪方法的整体流程图,附图中各符号的含义与发明内容部分相应符号的含义相同;
附图2是未结合幅值信息的phd滤波方法估计的多目标状态与真实值以及量测值对比,其中图(a)和图(b)分别为x坐标值对比和y坐标值对比;
附图3是本发明实施例中结合幅值信息的phd滤波多目标跟踪方法估计的多目标状态与真实值以及量测值对比,其中图(a)和图(b)分别为x坐标值对比和y坐标值对比。
具体实施方式
下面结合附图对本发明的结合幅值信息的phd滤波多目标跟踪方法进行详细描述。
不失一般性,假设在任意时刻,目标都在s=[-200,200]×[-200,200]的二维观测区域内运动,且目标可以在该区域随机的出现和消失,总仿真时间为k=50s,采样间隔t=1s;目标初始出现服从泊松模型,其密度函数γk(x)=0.2n(x|x0,qb),n(·|x0,qb)表示均值为x0,协方差为qb的高斯分布,其中x0=[020-2]t和qb=diag([105105]),目标平均信噪比
步骤1:根据以上仿真条件进行变量初始化
(1)总仿真时间k=50s,采样间隔t=1s;
(2)目标出现初始数γ0=0.2,代表一个目标所需要的粒子数l0=800,搜索新目标的粒子数jk=2000;
(3)目标平均信噪比
(4)目标出现初始分布γk(x)和杂波分布κk(z)分别为
γk(x)=0.2n(x|x0,qb)
(5)扩展单目标状态转移矩阵
(6)扩展过程噪声协方差
步骤2:按发明内容部分步骤2所述的方法进行滤波器初始化;
步骤3:按发明内容部分步骤3所述的方法获取当前时刻的量测;
步骤4:按发明内容部分步骤4所述的方法进行粒子集预测;
步骤5:按发明内容部分步骤5所述的方法进行粒子集更新;
步骤6:按发明内容部分步骤6所述的方法进行目标状态估计;
步骤7:循环执行发明内容部分步骤3~步骤6,直至雷达关机。
本实施例中,未结合幅值信息的phd滤波方法几乎不能实现对多目标的跟踪(见附图2),通过有效结合目标的幅值信息,本发明的方法可以实现对多目标的有效跟踪(见附图3);此外,本实施例中平均每帧的噪声点数为μ=1000,雷达虚警概率