一种杂波环境下的认知雷达目标跟踪方法与流程

本发明属于雷达数据处理
技术领域：
，具体涉及一种杂波环境下的认知雷达目标跟踪方法。
背景技术：
：传统的目标跟踪雷达(TraditionalActiveRadar,TAR)仅发射固定波形，在接收端采用自适应信号处理算法提高跟踪性能，而雷达的量测、分辨率等指标很大程度上取决于发射波形参数，当环境产生变化时，仅采用自适应滤波很难达到理想的跟踪效果，因此，自适应技术被应用到发射端，进而出现了自适应雷达(Fore-activeRadar,FAR)，它可以动态地选取发射波形。与此同时，基于知识辅助的雷达信号和数据处理方法也取得了显著的成果与进展，通过使用先验知识可进一步提高雷达的自适应能力。基于这两方面的研究成果，加拿大学者SimonHaykin于2006年首次提出了认知雷达(CognitiveRadar,CR)的概念。认知雷达通过环境的反射回波对环境进行感知，并利用感知的信息和其他先验知识对接收端和发射端进行联合自适应调整，实现对目标有效而稳定的跟踪，从而大幅提高了雷达目标跟踪性能。目前，认知雷达研究的对象主要集中在目标跟踪这一问题上。针对认知雷达目标跟踪这一问题，SimonHaykin等进行了一系列的研究，包括线性和非线性目标跟踪，但其所有的研究都是在假定无杂波这一前提条件下进行的。如果目标处于杂波环境下，目标回波和杂波回波混在一起，现有技术由于无法确定认知雷达感知存储器的输入而无法正常运行。技术实现要素：有鉴于此，本发明针对杂波环境下，目标回波和杂波回波混在一起，现有技术无法确定认知雷达感知存储器的输入而无法正常运行的问题，提供了一种杂波环境下的认知雷达目标跟踪方法。为了解决上述技术问题，本发明公开了一种杂波环境下的认知雷达目标跟踪方法，该方法包括以下步骤：步骤S1：获取当前时刻的所有量测值，并计算落入相关波门内的量测值；步骤S2：计算落入相关波门内的每个量测值来自目标的概率；步骤S3：根据落入相关波门内的每个量测值及其来自目标的概率，计算当前时刻的等效量测值；步骤S4：将等效测量值作为感知存储器的输入，输出得到当前时刻的系统方程参数；步骤S5：根据系统方程参数和等效量测值，进行状态和量测预测，得到状态预测值、量测预测值和状态预测协方差；步骤S6：将获得的状态预测协方差作为执行存储器的输入，输出得到当前时刻的发射波形参数集；步骤S7：对于当前时刻的发射波形参数集，计算每种波形参数下的新息协方差和状态误差协方差；步骤S8：进行最优波形参数选择：从获得的发射波形参数集中选取最佳的发射波形参数，使代价函数最小；步骤S9：根据获得的最优波形参数，计算当前时刻的新息协方差、状态误差协方差以及状态估计值。与现有技术相比，本发明可以获得包括以下技术效果：1、解决了现有的认知雷达目标跟踪技术在杂波环境下无法正常运行的问题。2、提高了杂波环境下的雷达目标跟踪性能。在杂波环境下，本发明方法的目标跟踪性能优于现有的TAR和FAR目标跟踪技术。当然，实施本发明的任一产品并不一定需要同时达到以上所述的所有技术效果。附图说明此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：图1为本发明杂波环境下的认知雷达目标跟踪方法一个实施例的流程图；图2为仿真例中本发明方法与现有的TAR和FAR目标跟踪方法的距离均方根误差曲线；图3为仿真例中本发明方法与现有的TAR和FAR目标跟踪方法的速度均方根误差曲线。具体实施方式以下将配合实施例来详细说明本发明的实施方式，藉此对本发明如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。本发明杂波环境下的认知雷达目标跟踪方法，如图1所示，具体按照以下步骤实施：步骤S1：获取当前时刻的所有量测值，并计算落入相关波门内的量测值。所述的相关波门是指以跟踪目标的量测预测值为中心，用来确定来自目标的量测值可能出现范围的一块区域。相关波门的形状根据实际应用进行设计，在本发明的一个实施例中，优选椭圆(球)相关波门，并可采用以下步骤对计算落入相关波门内的量测值进行详细描述：若当前时刻k等于1，则将当前时刻获取的所有量测值作为落入相关波门内的量测值。若当前时刻k大于等于2，则计算当前时刻获取的每一个量测值与量测预测值之间的相关误差：其中，zk,i为当前时刻获取的第i个量测值，为量测预测值，vi为zk,i与之间的相关误差，N为当前时刻获取的量测值的总数，Sk为新息协方差矩阵，Sk1为Sk的逆矩阵。实施时，和Sk分别由上一时刻计算的步骤S5和S9得到。当上述某个量测值的相关误差小于等于相关门限时，则该量测值落入相关波门内。相关门限的取值决定了量测值落入相关波门内的概率(即门概率)，在本发明的一个实施例中，优选等于16。步骤S2：计算落入相关波门内的每个量测值来自目标的概率：其中，i为落入相关波门内的第i个量测值来自目标的概率，M为落入相关波门内的所有量测值的总数，其中，为相关门限，nz为量测的维数，与相关波门的面(体)积有关，PD为目标检测概率，PG为门概率。在本发明的一个实施例中，设目标检测概率PD0.9，量测维数nz＝2，采用椭圆相关波门，相关门限16，分析可知：PG0.9997。在实施时，若相关波门内没有量测，即M等于0，则跳过步骤S2，直接进入步骤S3。步骤S3：根据落入相关波门内的每个量测值及其来自目标的概率，计算当前时刻的等效量测值。若M大于0，则当前时刻的等效量测为：若M等于0，则当前时刻的等效测量值为：其中，zk为当前时刻的等效量测值，zk,i和i分别为当前时刻落入相关波门内的第i个量测值及其来自目标的概率，为量测预测值。步骤S4：将等效测量值作为感知存储器的输入，输出得到当前时刻的系统方程参数。在本发明的一个实施例中，所述的感知存储器由包含一个输入层、一个隐含层和一个输出层的三层感知神经网络构成，其输入为步骤S3中获得的等效量测zk，输出为与等效量测最匹配的系统方程参数，具体来说，就是状态转移矩阵Fk、量测矩阵Hk和系统噪声协方差矩阵Qk。步骤S5：根据获得的系统方程参数和等效量测值，进行状态和量测预测，得到状态预测值、量测预测值和状态预测协方差，具体如下：其中，和Pk|k1分别为当前时刻的状态预测值、量测预测值和状态预测协方差，xk1为上一时刻的状态估计值，Ek1为上一时刻的状态误差协方差。实施时，若当前时刻k大于等于2，则xk1和Ek1由上一时刻计算的步骤S9得到；若当前时刻k等于1，则xk1和Ek1分别取为步骤S3中获得的等效量测值zk和步骤S4中获得的系统噪声协方差矩阵Qk。步骤S6：将步骤S5中得到的状态预测协方差作为执行存储器的输入，输出得到当前时刻的发射波形参数集。在实施时，所述的执行存储器由一个发射波形库Θ以及包含一个输入层、一个隐含层和一个输出层的三层感知神经网络共同构成。发射波形库Θ中包含多种事先设计好的用于探测环境和目标的发射波形参数，波形库的设计由具体应用环境而定，在本发明的一个实施例中，采用线性调频信号作为发射波形，涉及到的波形参数包括包络持续时间和调频斜率b。三层感知神经网络的输入为步骤S5中获得的状态预测协方差Pk|k1，输出为与之匹配的一个发射波形参数集Θk{θk,1,…,θk,r}，其中，Θk为Θ的子集，r为该子集种包含波形参数的总数量。步骤S7：对于当前时刻的发射波形参数集，计算每种波形参数下的新息协方差和状态误差协方差。在实施时，可采用以下步骤对所述的计算每种波形参数下的新息协方差和状态误差协方差进行详细描述：新息协方差为：状态误差协方差为：其中，Rk(θk,j)和分别为发射波形参数为θk,j时的量测噪声协方差和增益。如步骤S6中所述，在本发明的一个实施例中，采用线性调频信号作为发射波形，则波形参数θk,j由包络持续时间j和调频斜率bj构成，分析可得发射波形参数为θk,j{j,bj}时的量测噪声协方差Rk(θk,j)为：其中，为信噪比，fc为发射信号的中心频率，c为光速。步骤S8：进行最优波形参数选择：从步骤S6中得到的发射波形参数集中选取最佳的发射波形参数，使代价函数最小。在本发明的一个实施例中，所述的代价函数可以选择为状态误差协方差矩阵Ek(θk,j)的迹，则最优波形参数就是让Ek(θk,j)的迹为最小的波形参数：其中，Tr()表示矩阵的迹，为代价函数最小时对应的波形参数，即为最优波形参数。如步骤S7中所述，在本发明的一个实施例中，采用线性调频信号作为发射波形，则最优波形参数为步骤S9：根据步骤S8中得到的最优波形参数，计算当前时刻的新息协方差、状态误差协方差以及状态估计值，具体如下：其中，Sk、Ek和xk分别为当前时刻最优波形参数下的新息协方差、状态误差协方差以及状态估计值，Rk和分别为下的量测噪声协方差和增益。下面通过一个仿真例进一步对本发明方法的效果进行详细说明。仿真例设有一个单目标在二维量测空间中进行线性运动，状态向量为xr,vT，其中，r和v分别表示目标的距离和速度，设目标初始距离为3km，且以200m/s的速度向与雷达相反的方向匀速直线运动，雷达采样间隔为25ms。设杂波密度为0.0001，目标检测概率为PD0.9，采用椭圆相关波门且相关门限16，对应的门概率为PG0.9997。采用线性调频信号作为发射波形，波形参数库为：Θ,b，其中，105:105:104为线性调频信号的包络持续时间，单位为s；b1011:1010:1010∪1010:1010:1011为线性调频信号的调频斜率，单位为Hz/s。分别采用本发明方法以及现有的TAR和FAR方法进行目标跟踪，图2和图3分别为三种方法进行目标跟踪的距离和速度均方根误差曲线，图中结果均为50次蒙特卡洛实验的平均结果。将图2和图3中的跟踪误差曲线分别在时间轴上作统计平均，得到三种方法进行目标跟踪的距离和速度均方根均值如表1所示。表1仿真例中本发明方法与现有的TAR和FAR目标跟踪方法的距离和速度均方根均值均方根误差均值距离(m)速度(m/s)TAR6.90247.2573FAR2.98183.0452本发明0.24052.1813由图2、图3和表1的结果可以看到：在杂波环境下，无论是距离跟踪还是速度跟踪，本发明方法的跟踪误差均明显小于现有的TAR和FAR目标跟踪方法。该仿真例的结果说明了本发明提供的一种杂波环境下的认知雷达目标跟踪方法的有效性，以及该方法相对于现有的TAR和FAR目标跟踪方法具有更优的跟踪性能。上述说明示出并描述了发明的若干优选实施例，但如前所述，应当理解发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围，则都应在发明所附权利要求的保护范围内。当前第1页1 2 3