一种反射点不确定场景下多路径融合目标检测算法制造技术

本发明专利技术公开了一种反射点不确定场景下多路径融合目标检测算法，属于雷达技术领域。本发明专利技术以迭代的方式同时估计目标和反射点的状态，用一批数据使用滑动窗口方式用于航迹初始化。每个迭代由两个步骤组成：预测和更新。在预测阶段，通过上个时刻最后一次迭代的估计状态来预测反射点的状态，在此基础上，通过多遍网格搜索全局最优解来估计目标的状态，在更新阶段，利用更新后的目标状态来重新估计反射点状态。重复上述步骤，迭代过程基于给定的最大迭代次数或判断是否过门限，如果迭代次数达到最大迭代次数或过门限则停止迭代，输出目标位置和反射点的位置。本发明专利技术提高了目标状态估计精度，对低空小目标进行有效的检测。

【技术实现步骤摘要】
一种反射点不确定场景下多路径融合目标检测算法
本专利技术属于雷达
，具体涉及一种反射点不确定场景下多路径融合目标检测算法。
技术介绍
低小慢目标指的是飞行高度较低、体积相对较小、飞行速度慢的小型飞行器，典型的低小慢目标包括无人机、航模、滑翔伞等。最近几年随着我国民用无人机的发展，无人机性能达到国际一流水平，普通民众也能消费得起民用无人机。与此同时也带来了一些安全上的隐患，比如在民航机场经常会出现未被允许飞行的小型无人机，这对于机场的正常运行造成严重的影响，甚至会影响飞机的起飞与降落。另外，在一些活动场所中经常出现一些无人机，同样严重影响到了活动的安全保障。因此，如何有效管理和控制无人机，成为了当前雷达监测系统亟待解决的问题。为了提高在低信噪比、高杂波场景下目标的检测与跟踪性能，在检测系统中应尽可能挖掘传感器接收量测中的目标信息。在一些应用场景中，由于电磁波传播存在多径效应，传感器在接收回波中同一帧中能产生多个目标量测。比如在城市环境中，高楼大厦随处可见，对于监测系统来说目标很容易被遮挡，存在视觉盲区，再者电磁波信号在传播的时候不仅可以通过目标直接反射回到接收机，还会通过建筑物反射回到接收机，这种情况下目标会有多条传播路径，如果不对其进行处理的话会影响目标的正常检测。换个角度来说，通过其他路径反射回来的目标信号也是目标信息的一部分，如果能将这一部分目标信息利用起来，势必会增加目标检测性能。由于目标在不停的运动，反射点的位置一般也会随着目标的运动而运动，通常情况下只知道建筑物的大概位置，但是并不知道反射点在建筑物的具体位置。现有的基于多路径目标检测算法假定反射点已知，量测模式确定，如果在反射点未知场景下使用传统的多路径ML-PMHT(最大似然概率多假设跟踪)目标检测算法，会使得目标估计误差增大，而且容易产生虚假航迹，失跟率高。现有的多路径ML-PMHT目标检测算法，并没有将反射点不确定性建模到算法中，由于反射点初始位置不确定，导致量测误差增加，因此通过搜索优化得到的目标状态距离真实目标状态差距较大，并且容易将由杂波形成的对数似然比(LLR)误认为是目标形成的对数似然比(LLR)，从而造成虚假航迹。
技术实现思路
本专利技术的目的是克服ML-PMHT算法在反射点不确定多径环境中应用的不足，提供一种反射点不确定场景下高效多路径信息融合的最大似然概率多假设跟踪(ML-PMHT)的目标检测算法，具体来说是一种基于城市复杂环境的信号多路径特点，以及反射点不确定的传播特点，对低空小目标(比如小型飞行器)进行检测的方法。本专利技术所述方法通过对目标、多路径观测函数显式建模，多并对多路径观测信息进行融合处理，从而对低信噪比、高杂波、反射点不确定场景下的低空小目标的检测。本专利技术所述方法的思路是，以迭代的方式同时估计目标和反射点的状态。用一批数据使用滑动窗口方式用于航迹初始化。在每次迭代中，使用相同的测量集。每个迭代由两个步骤组成：预测和更新。在预测阶段，通过上个时刻最后一次迭代的估计状态来预测反射点的状态，在此基础上，通过多遍网格搜索全局最优解来估计目标的状态，由于LLR的计算使用的是估计的反射点而不是使用真实反射点位置来计算的，因此LLR的全局最优并不能得到真实的目标状态。在更新阶段，利用更新后的目标状态来重新估计反射点状态。重复上述步骤，该迭代过程基于给定的最大迭代次数或判断是否过门限，如果迭代次数达到最大迭代次数或过门限则停止迭代，输出目标位置和反射点的位置。从而实现城市低信噪比、高杂波、反射点不确定场景下低空小目标的检测。本专利技术所提出的技术问题是这样解决的：一种反射点不确定场景下多路径融合目标检测算法，包括以下步骤：步骤1.初始化多路径ML-PMHT算法环境参数：初始化观测环境各项参数：监控空间V，反射点数L，反射点初始状态表示均值为协方差为的高斯分布，1≤l≤L，每个反射点对应一个路径，量测通过传播路径l来源于目标n的先验概率为πn，l，量测来源于杂波的先验概率π00，当前迭代次数Iter＝1，最大迭代次数为Imax；目标初始化场景下参数向量分别表示为：其中，k为批处理的采样时刻数，1≤k≤Nw，Nw为批处理长度，X为Nw个时刻的目标状态参数，Xf为Nw个时刻反射点状态参数，Z为Nw个时刻的量测集合；在k时刻目标的状态向量表示为：Xk＝{x1，k，...，xn，k，...，xN，k}n＝1，...，N在k时刻反射点状态向量表示为：在k时刻收到的量测集合表示为：其中，xn，k为k时刻第n个目标的状态，N为目标的个数，表示k时刻第l个反射点的状态，L为反射点的个数，zk，j表示传感器k时刻收到的第j个量测，mk表示传感器k时刻收到的量测个数；对目标的运动过程和对不同传播路径的观测方程分别表示为：xn，k＝Fn，kxn，k-1+vn，k其中，Fn，k为目标运动矩阵，xn，0为第n个目标运动的初始状态，hl为xn，k经路径l的观测函数；zn，k，l为xn，k通过路径l产生的目标量测；vn，k和ωn，l，k分别为目标过程噪声和观测噪声，都是均值为零的高斯白噪声，协方差矩阵分别为Qn，k和Rn，l，k；步骤2.预测反射点状态与方差：其中，表示k时刻第l个反射点预测状态，Fl，k表示反射点的运动矩阵，vl，k表示反射点过程噪声，是零均值的高斯白噪声，协方差矩阵为表示k时刻反射点l的预测状态协方差，表示反射点l的初始预测状态协方差；是测量模型对反射点状态的雅可比矩阵，其中表示关于的偏导，上标T表示转置；步骤3.根据预测的反射点状态，搜索所有可能的目标状态：步骤3-1.令N′＝1，目标状态集合X′为空向量，步骤3-2.搜索单目标LLR全局最大值；其中，pl[zk，j|xn′，k]表示通过路径l以第n′个目标为中心的高斯概率；其中，表示高斯概率密度函数，高斯变量zk，j的均值为协方差为Rn′，l，k，且有：步骤3-3.判断单目标LLR全局最大值是否大于检测门限，若是，确认为新目标，即将加入到向量X′中，转至步骤3-4；否则，结束检测过程，转至步骤4；步骤3-4.从量测集合中剔除与已确认为新目标关联的量测，剔除依据为目标后验关联概率最高的量测将被剔除；ML-PMHT后验关联概率wj，n′，l，k计算公式为：步骤3-5.令N′＝N′+1，转至步骤3-2；此时搜索状态空间为4N′维，由于有(N′-1)个目标已经确认，在新的搜索过程无需搜索已经搜出的4(N′-1)维，所以新的搜索阶段只需要对未知的4维状态进行搜寻；步骤4.更新反射点状态：步骤4-1.计算k时刻的量测zk，j通过第l条路径来源第n个目标的后验概率：步骤4-2.计算每个目标在每条路径下的综合量测及综合协方差...

【技术保护点】
1.一种反射点不确定场景下多路径融合目标检测算法，其特征在于，包括以下步骤：/n步骤1.初始化多路径ML-PMHT算法环境参数：/n初始化观测环境各项参数：监控空间V，反射点数L，反射点初始状态

【技术特征摘要】
1.一种反射点不确定场景下多路径融合目标检测算法，其特征在于，包括以下步骤：
步骤1.初始化多路径ML-PMHT算法环境参数：
初始化观测环境各项参数：监控空间V，反射点数L，反射点初始状态表示均值为协方差为的高斯分布，1≤l≤L，每个反射点对应一个路径，量测通过传播路径l来源于目标n的先验概率为πn,l，量测来源于杂波的先验概率π00，当前迭代次数Iter＝1，最大迭代次数为Imax；
目标初始化场景下参数向量分别表示为：









其中，k为批处理的采样时刻数，1≤k≤Nw，Nw为批处理长度，X为Nw个时刻的目标状态参数，Xf为Nw个时刻反射点状态参数，Z为Nw个时刻的量测集合；
在k时刻目标的状态向量表示为：
Xk＝{x1,k,...,xn,k,...,xN,k}n＝1,...,N
在k时刻反射点状态向量表示为：



在k时刻收到的量测集合表示为：



其中，xn,k为k时刻第n个目标的状态，N为目标的个数，表示k时刻第l个反射点的状态，L为反射点的个数，zk,j表示传感器k时刻收到的第j个量测，mk表示传感器k时刻收到的量测个数；
对目标的运动过程和对不同传播路径的观测方程分别表示为：
xn,k＝Fn,kxn,k-1+νn,k



其中，Fn,k为目标运动矩阵，xn,0为第n个目标运动的初始状态，hl为xn,k经路径l的观测函数；zn,k,l为xn,k通过路径l产生的目标量测；νn,k和ωn,l,k分别为目标过程噪声和观测噪声，都是均值为零的高斯白噪声，协方差矩阵分别为Qn,k和Rn,l,k；
步骤2.预测反射点状态与方差：






其中，表示k时刻第l个反射点预测状态，Fl,k表示反射点的运动矩阵，vl,k表示反射点过程噪声，是零均值的高斯白噪声，协方差矩阵为表示k时刻反射点l的预测状态协方差，表示反射点l的初始预测状态协方差；是测量模型对反射点状态的雅可比矩阵，其中表示关于的偏导，上标T表示转置；
步骤3.根据预测的反射点状态，搜索所有可能的目标状态：
步骤3-1.令N'＝1，目标状态集合X'为空向量，
步骤3-2.搜索单目标LLR全局最大值；






其中，pl[zk,j|xn',k]表示通过路径l以第n'个目标为中心的高斯概率；<...

【专利技术属性】
技术研发人员：唐续，光昌国，王代维，程建，杨阳，董平，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51