一种基于目标模态估计的临近空间高超声速目标跟踪方法技术

本发明专利技术提供一种基于目标模态估计的临近空间高超声速目标跟踪方法，利用交互多模型跟踪算法跟踪目标，并实时动态估计目标运动模态，通过统计目标特性判断目标运动模态；最后，根据所估计的目标运动模态转入相应的单模匹配跟踪，避免了多模型间的竞争，解决了现有临近空间高超声速跟踪算法存在计算复杂，模型竞争大导致跟踪精度低等问题。本发明专利技术具有计算量小和跟踪精度高，有效的提升了跟踪系统的整体性能。

【技术实现步骤摘要】
一种基于目标模态估计的临近空间高超声速目标跟踪方法
本专利技术属于雷达信号处理技术，特别涉及对临近空间高超声速典型运动目标的跟踪技术。
技术介绍
临近空间指距离地面20～100Km高的空域，临近空间高超声速飞行器速度可达4～20马赫。目前，临近空间高超声速目标跟踪是跟踪领域的一个热点。临近空间机动目标运动的复杂性和运行环境的多变性，造成难以对此类目标建立准确跟踪模型，隐身技术的发展更加大了对临近空间高超声速目标跟踪的难度。目前，基于交互式多模型(IMM)结构跟踪算法，是一种公认最有效临近空间高超声速目标的跟踪算法。在文献“ResearchofMethodforTrackingHighSpeedandHighlyManeuveringTarget,InternationalConferenceonITSTelecommunicationsProceedings,1236-1239,2006”中，提出基于交互多模结构的跟踪算法，该算法利用多个不同机动模型交互对临近空间高超声速目标进行跟踪，相较于单模跟踪具有较广的覆盖范围，有较大的机动适应性；但该算法模型间相互竞争导致跟踪精度较差，且计算量较大。
技术实现思路
本专利技术所要解决的技术问题是，提供一种跟踪精度较高，且计算量较小的临近空间高超声速目标跟踪方法。本专利技术为解决上述技术问题采用的技术方案是，一种基于目标模态估计的临近空间高超声速目标跟踪方法，包括以下步骤：步骤1、将雷达采集到的量测数据进行目标航迹起始；步骤2、利用交互多模跟踪IMM算法估计匀速运动模型、匀加速运动模型和拐弯运动模型所对应的下一帧目标状态、状态协方差矩阵以及相应的模型转移概率；所述目标状态包括目标位置、速度与加速度；当连续统计L帧下匀速运动模型、匀加速运动模型和拐弯运动模型对应的转移概率、目标速度和加速度后能够确定目标当前所处的运动模态，则进入步骤3，否则返回步骤2；若在L帧时间内，匀速模型转移概率ucv一直保持最大，速度变化量Δvk≤τv，加速度变化量Δak≤τa，则判断目标为匀速运动阶段；若在L帧时间内，匀加速模型转移概率uca一直保持最大，速度变化量Δvk≥τv，加速度变化量Δak≤τa，则判断目标为匀加速运动阶段；若在L帧时间内，匀加速模型转移概率uct一直保持最大，速度变化量Δvk≥τv，加速度Δak≥τa，则判断目标处于跳跃运动阶段；其中，τv表示速度波动门限，τa表示加速度波动门限；步骤3、单模匹配跟踪：3-1确定目标运动模态后，转入单模匹配跟踪；初始化单模参数、目标初始状态；3-2利用单模匹配跟踪算法估计下一帧目标状态；3-3计算当前时刻k的归一化残差平方εv为：vx为当前量测残差，Sx表示当前残差协方差矩阵，下标x表示与目标真实运动相匹配的运动模型，即x∈{CV,CA,CT}，CV表示匀速运动，CA表示匀加速运动，CT表示跳跃运动，(·)T表示矩阵转置；当归一化残差平方εv小于等于门限τ时，认为目标运动状态不变，返回步骤3-2，当归一化残差平方εv大于门限τ时，认为目标运动状态发生改变，返回到步骤2。本专利技术利用交互多模型跟踪算法跟踪目标，并实时动态估计目标运动模态，通过统计目标特性判断目标运动模态；最后，根据所估计的目标运动模态转入相应的单模匹配跟踪，避免了多模型间的竞争，解决了现有临近空间高超声速跟踪算法存在计算复杂，模型竞争大导致跟踪精度低等问题。本专利技术的有益效果是，具有计算量小和跟踪精度高，有效的提升了跟踪系统的整体性能。附图说明图1基于交互多模型和目标运动模态动态估计的跟踪算法流程图。图2新跟踪算法和交互多模跟踪算法跟踪误差曲线图。图3新跟踪算法和交互多模跟踪算法在不同检测概率下跟踪成功概率曲线。图4新跟踪算法和交互多模跟踪算法在不同采样率下跟踪成功概率曲线。具体实施方式采用计算机仿真的方法进行验证本专利技术效果，通过MATLAB-R2010b，实施步骤如图1所示：步骤1、雷达数据转换被跟踪目标进入地基雷达扫描区域，雷达所采集到数据为量测数据(r,θ,)和目标的多普勒信息(Doppler)fd，r表示目标距离、θ表示目标方位、表示目标俯仰角。利用修正量测无偏转换法MUCMKF将极坐标系量测数据(r,θ,)转换到直角坐标系下，得到转换后的量测数据其中，分别表示转换后量测在直角坐标系下目标位置；其中，补偿因子分别表示方位角和俯仰角度的均方差，θm,εm分别表示方位角和俯仰角。步骤2、目标航迹起始2.1储存3帧转换后的量测数据，第k时刻数据Z(k)＝{z1(k),z2(k),z3(k),...zm1(k)}，第k+1时刻得到的数据集为Z(k+1)＝{z1(k+1),z2(k+1),z3(k+1),...zm2(k+1)}，第k+2时刻得到的数据集为Z(k+2)＝{z1(k+2),z2(k+2),z3(k+2),...zm3(k+2)}，遍历关联连续三帧数据，利用高度约束hmin<h<hmax和速度约束vmin<v<vmax和Doppler信息fd剔除部分杂波点迹。其中，zi(k)表示第k时刻第i个量测，m1,m2,m3分别表示3个时刻的量测个数，hmin,hmax分别表示最低和最高飞行高度，vmin,vmax为最小和最大飞行速度。2.2将余下的连续3帧数据投影到两个平面(x-z和y-z平面)，再将两平面数据变换到参数空间，ρ＝xcosθ+ysinθ，其中，(θ,ρ)为参数空间中的坐标，(x,y)为直角坐标系下观测数据；再利用修正霍夫Hough变换算分别对两平面量测数据进行航迹起始，对比关联两平面起始的航迹，进而确定目标航迹起始。为了尽可能抑制虚假航迹，这里优选使用修正霍夫Hough变换。本领域其它的航迹起始方法对临近空间目标航迹起始会导致虚假航迹增多，不利于目标跟踪，比如直接航迹起始、修正逻辑法。步骤3、利用交互多模跟踪IMM算法跟踪目标3.1目标航迹起始后，首先初始化目标状态各模型转移概率u＝[ucvucauct]和各模型状态协方差矩阵Pcv,Pca,Pct。假设前3个时刻的量测值分别为Z(1),Z(2),Z(3)，且Z(i)＝[x(i),y(i),z(i)]T，观察噪声协方差矩阵R；其中，R通过量测转换方程直求得；系统状态初始化：其中，分别表示x轴、y轴和z轴上的位置、速度和加速度；[ucv,uca,uct],[Pcv,Pca,Pc本文档来自技高网...

【技术保护点】
一种基于目标模态估计的临近空间高超声速目标跟踪方法，其特征在于，包括以下步骤：步骤1、将雷达采集到的量测数据进行目标航迹起始；步骤2、利用交互多模跟踪IMM算法估计匀速运动模型、匀加速运动模型和拐弯运动模型所对应的下一帧的模型转移概率以及目标状态；所述目标状态包括目标位置、速度与加速度；当连续统计L帧下匀速运动模型、匀加速运动模型和拐弯运动模型对应的转移概率、目标速度和加速度后能够确定目标当前所处的运动模态，则进入步骤3，否则返回步骤2；若在L帧时间内，匀速模型转移概率ucv一直保持最大，速度变化量△vk≤τv，加速度变化量△ak≤τa，则判断目标为匀速运动阶段；若在L帧时间内，匀加速模型转移概率uca一直保持最大，速度变化量△vk≥τv，加速度变化量△ak≤τa，则判断目标为匀速运动阶段；若在L帧时间内，拐弯模型转移概率uct一直保持最大，速度变化量△vk≥τv，加速度△ak≥τa，则判断目标处于跳跃运动阶段；其中，τv表示速度波动门限，τa表示加速度波动门限；步骤3、单模匹配跟踪：3‑1确定目标运动模态后，转入单模匹配跟踪；初始化单模参数、目标初始状态；3‑2利用单模匹配跟踪算法预测下一帧目标状态；3‑3计算当前时刻k的归一化残差平方εv为：vx为当前量测残差，Sx表示当前残差协方差矩阵，下标x表示与目标真实运动相匹配的运动模型，即x∈{CV,CA,CT}，CV表示匀速运动，CA表示匀加速运动，CT表示跳跃运动，(·)T表示矩阵转置；当归一化残差平方εv小于等于门限τ时，认为目标运动状态不变，返回步骤3‑2，当归一化残差平方εv大于门限τ时，认为目标运动状态发生改变，返回到步骤2。...

【技术特征摘要】
1.一种基于目标模态估计的临近空间高超声速目标跟踪方法，其特征在于，包括以下步骤：步骤1、将雷达采集到的量测数据进行目标航迹起始；步骤2、利用交互多模跟踪IMM算法估计匀速运动模型、匀加速运动模型和拐弯运动模型所对应的下一帧的模型转移概率以及目标状态；所述目标状态包括目标位置、速度与加速度；当连续统计L帧下匀速运动模型、匀加速运动模型和拐弯运动模型对应的转移概率、目标速度和加速度后能够确定目标当前所处的运动模态，则进入步骤3，否则返回步骤2；若在L帧时间内，匀速模型转移概率ucv一直保持最大，速度变化量△vk≤τv，加速度变化量△ak≤τa，则判断目标为匀速运动阶段；若在L帧时间内，匀加速模型转移概率uca一直保持最大，速度变化量△vk≥τv，加速度变化量△ak≤τa，则判断目标为匀加速运动阶段；若在L帧时间内，拐弯模型转移概率uct一直保持最大，速度变化量△vk≥τv，加速度△ak≥τa，则判断目标处于跳跃运动阶段；其中，τv表示速度波动门限，τa表示加速度波动门限；步骤3、单模匹配跟踪：3-1确定目标运动模态后，转入单模匹配跟踪；初始化单模参数、目标初始状态；3-2利用单模匹配跟踪算法预测下一帧目标状态；3-3计算当前时刻k的归一化残差平方εv...

【专利技术属性】
技术研发人员：易伟，董天发，苟清松，郝凯利，崔国龙，孔令讲，杨建宇，李小龙，夏玫，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51