强杂波环境下相控阵雷达的GM-PHD目标跟踪方法与流程

本发明属于雷达数据处理领域，具体涉及一种强杂波环境下相控阵雷达的gm-phd目标跟踪方法。
背景技术：
：雷达多目标跟踪一直以来都是一个重要的理论和实践问题。近年来，随着统计理论的不断发展，mahler提出了基于概率假设密度(phd)滤波器的雷达目标跟踪方法。相比与联合概率数据互联算法(jpda)和多假设跟踪算法(mht)等传统的多目标跟踪方法而言，该方法有效避免了数据关联带来的困难和麻烦，受到了国内外学者的广泛关注和研究。传统phd滤波在进行目标跟踪时，一般在固定位置进行新生目标的起航，然而在实际应用时通常难以得到目标起始位置的先验信息，因此很难正确起始航迹。为改善这一问题，有学者提出了基于量测的phd新生目标强度计算方法，在量测处设立新生目标以解决该问题。该方法虽可解决无先验信息条件下phd滤波中的目标起始问题，但在强杂波环境下会起始过多的虚假航迹，导致phd滤波的效率降低。此外,对于相控阵雷达体制下的目标跟踪而言，随着雷达功能的切换，在tws过程中加入tas等，会随之带来雷达扫描周期不固定的问题，对航迹的预测和更新带来影响，因此要将phd滤波方法应用于相控阵雷达，必需加以有效措施来解决该问题。技术实现要素：本发明的目的在于提供一种强杂波环境下相控阵雷达的gm-phd目标跟踪方法，解决相控阵雷达变扫描周期情况下的目标跟踪问题，同时提高强杂波环境下phd滤波的效率，减少虚假航迹的产生。实现本发明目的的技术解决方案为：一种强杂波环境下相控阵雷达的gm-phd目标跟踪方法，包括以下步骤：步骤1、对量测数据进行预处理，标记杂波区域；步骤2、然后根据时间戳外推航迹，利用筛选得到的量测数据更新航迹数据；步骤3、对剩余量测数据与临时航迹关联处理，计算新生目标的强度。与现有技术相比，本发明的显著优点为：(1)本发明通过提取每帧数据的时间戳作为时间基础进行航迹的预测，而不是采用固定的处理周期，能有效适应相控阵雷达变扫描周期情况下的目标跟踪；(2)本发明通过对量测的统计处理对可能来源于杂波的量测进行自动识别和标记，并且在计算phd滤波的新生目标强度时，通过设置起航条件和约束条件来对可能产生新目标的量测进行选择和确认，有效抑制了虚假目标的产生，提高了phd滤波的计算效率。附图说明图1是本发明强杂波环境下相控阵雷达的gm-phd目标跟踪方法流程图。图2是采用本发明方法对仿真数据的目标轨迹跟踪结果示意图。图3是采用本发明方法对仿真数据的目标数估计结果示意图。图4是采用本发明方法对实测数据的目标轨迹跟踪结果示意图。具体实施方式本发明提供一种杂波环境下相控阵雷达的gm-phd目标跟踪方法，该方法的基本思路是不根据预设的扫描周期进行航迹的预测，而是在接收每一帧数据后提取当前时间戳来为航迹预测提供时间基础，因此可适用于相控阵雷达扫描周期变化的特点；此外，在预处理阶段，本发明通过对量测的统计处理对可能来源于杂波的量测进行自动识别和标记，并且在phd滤波中通过设置约束条件来对产生新生目标的潜在量测进行选择，以抑制杂波、减少虚假航迹。如图1所示，本发明的目标跟踪方法包括以下步骤：步骤1、对量测数据进行预处理，标记杂波区域；接收雷达信号处理器输入的目标量测位置信息，提取当前数据的时间戳tk，并对量测点进行预处理，标记和剔除可能来源于杂波的量测。将雷达接收的量测数据划分网格，采用滑窗法实时统计和计算最近5-10个周期的以下4个统计量：1)每帧数据中的平均量测数s1；2)有量测点落入的网格数目s2；3)当前量测数据中落入每个网格点的量测数s3；4)根据s1和s2计算的到网格平均量测点数s4＝s1/s2；比较s3和s4，将s3大于s4的网格判别为杂波网格，并对该部分网格域置上标识；再将杂波网格的s3值与设定的杂波判别门限lup和llow比较，当s3≤llow时该网格判定为弱杂波区，当lup＞s3＞llow时该网格判定为中杂波区，当s3≥lup该网格判定为强杂波区；对弱杂波区采用3点累积航迹起始条件，对于中杂波区再增加1-3点的航迹起始所需的积累点数，而对落入强杂波区域的量测直接剔除；lup和llow的取值参考下式进行设置：max(s3)＞lup＞min(s3)＞llow＞s4假设剔除强杂波后，第k帧数据的量测集合为其中mk为量测数，中各分量分别为量测的距离、方位角和多普勒。步骤2:根据时间戳外推航迹，利用筛选得到的量测数据更新航迹数据；(1)目标状态预测假设在处理完第k-1帧数据后，可靠航迹的高斯混合phd(gm-phd)为其中，和分别为目标随机集强度第i个高斯成分的权值、均值和方差，n(·；m,p)为均值为m，方差为p的高斯分布，nk-1为k-1时刻的可靠航迹数。同时，可靠航迹的标签集为可靠航迹的时间戳为tk-1，即航迹上最新点的时间为tk-1。将所有可靠航迹都外推至tk时刻，预测phd可表示为其中，dk|k-1为tk时刻的预测概率假设密度函数，z1:k-1为至tk-1时刻的量测序列，ps为目标存活概率，fk|k-1为状态转移矩阵，fk|k-1中的时间间隔为δt＝tk-tk-1，qk|k-1为系统噪声协方差。在预测后，各高斯项的标签不发生变化，因此可靠航迹的标签集也保持不变，即预测后各可靠航迹的时间戳变更为tk。(2)量测筛选。根据所设定的量测筛选条件，对接收到的量测进行筛选，以确定可用于更新可靠航迹的量测。首先通过计算量测与航迹间的欧式距离进行粗筛选，然后通过计算量测与航迹的新息加权范数进一步精确筛选。(2a)粗筛选将量测zk中的位置信息由极坐标转换为直角坐标下的量测值计算量测与可靠航迹预测值间的欧式距离其中当δd(i,j)大于预设门限时，将该量测筛除，其余量测进入下一步精筛选。(2b)精筛选可靠航迹与量测间的新息加权范数为其中σr、σθ和σd分别为距离测量误差、方位角测量误差和多普勒测量误差。当小于预设门限时，认为该量测可能源于可靠航迹，可用于下一步可靠航迹的状态更新。而通过筛选排除的量测则与临时航迹做关联处理，以计算phd滤波中的新生目标信息。(3)可靠航迹的目标状态更新。假设筛选得到的量测集合为利用筛选得到的量测集合对目标状态进行更新其中其中，pd为检测概率，κ(z)为杂波强度，hk(·)为系统观测方程，rk为测量噪声协方差矩阵，z为量测集合zc,k中的元素，i为单位矩阵。更新后的目标标签集由预测部分和量测更新部分组成，预测部分标签集与tk|k-1相同，而一个目标状态与nc,k个量测作用会产生nc,k个高斯项，即每个目标标签产生nc,k副本。更新后的目标标签集为其中为t(i)的nc,k个副本标签。对更新后的高斯项进行枝剪和合并处理，对处理后仍具有相同标签的分量，保持权值大的目标的标签不变，对权值小的分量重新分配新的标签。步骤3：对剩余量测数据与临时航迹关联处理，计算新生目标的强度。根据上一步设定的量测筛选条件得到的与当前可靠航迹都不相关联的量测，与临时航迹做关联处理，满足条件的临时航迹转为可靠航迹，并计算phd滤波中的新生目标强度。量测与临时航迹的关联按以下多假设航迹起始的思路进行：(a)对于长度为1的临时航迹(即航迹头)，用速度法建立相关波门，对落入初始相关波门内的所有量测均建立临时航迹，即尝试起始多条航迹；(b)对于长度大于1的临时航迹进行外推，以外推点为中心建立相关波门，对落入相关波门内的所有量测均与临时航迹互联，即尝试起始多条航迹；为提高关联的可靠性、减少误关联，在步骤(a)和(b)中还可测设置一定的约束条件。考虑目标的运动特征，本发明在量测和临时航迹互联时还采用如下3个约束条件：最大速度约束：最小速度约束：最大方向夹角约束：其中，rk为目标的位置向量，δt为时间间隔。(c)若持续一段时间没有量测与某条临时航迹关联，则撤销此临时航迹；(d)未与临时航迹产生关联的量测点作为航迹头，加入临时航迹文件中；(e)持续上述步骤，直至某条临时航迹累计关联上一定数量的量测点(具体数量可根据杂波情况调节)，则将该临时航迹转为可靠航迹，并以航迹的起始状态作为phd滤波中新生目标的信息：其中，nk,b是k时刻满足航迹起始条件的临时航迹数；为新生目标的权值，可根据杂波强度设定，在强杂波区域可减小该值，而在弱杂波区域可适当提高；和为根据航迹上的历史点计算得到的目标初始状态和初始协方差。同时，将新生目标加入目标标签集中。下面结合实例对本发明的效果做进一步说明。实施例首先采用本发明方法对仿真数据进行处理，仿真数据是在强杂波环境下对多目标进行跟踪，设置6个目标，跟踪时间为0-60s，目标初始状态和出现时刻如表1所示。表1目标出现时刻和初始状态出现时刻(s)初始位置(m)速度(m/s)0(400,-300)(-10,5)0(400,-600)(-7,-4)20(-800,-200)(12,7)20(-200,800)(15,-10)40(100,50)(-20,-15)40(-200,500)(15,-5)从图2中可以看到，采用本发明中的方法可在杂波较强的情况下，对多目标进行较好的跟踪。图3为跟踪过程中对目标数估计的结果。可以看到本发明方法对目标数目进行了准确的估计，有效抑制了强杂波可能带来的虚假航迹。由于引入了前后量测点的信息来进行新生目标强度的计算，因此在图3中可以看到本发明对新生目标的响应会出现一定的滞后。实际中，可根据杂波强度调节航迹起始的点数和约束条件中的阈值，以提升对新生目标的响应速度。如图4所示，采用本发明对某相控阵雷达的实测数据进行处理，场景中包含一架进行机动转弯飞行的无人机目标和若干地面车辆，且回波数据中包含了大量杂波数据。相控阵雷达首先工作于tws跟踪模式，在发现无人机目标后，加入tas跟踪，而对其余目标保持tws跟踪。从处理结果可以看到，本发明方法可在密集的杂波环境下对多目标进行稳定跟踪，且跟踪轨迹连续，并有效抑制虚假目标的产生，同时也能够在雷达添加tas目标，扫描周期变化的情况下，对各目标良好的跟踪。当前第1页12