发明名称 一种基于熵权值法点迹质量评估的目标跟踪方法
摘要 该发明涉及一种基于点迹质量评估的最近邻域目标跟踪方法,属于雷达目标跟踪技术领域,特别涉及一种基于熵权值的多评估指标信息融合的目标跟踪方法。本发明提出了依据点迹的多种特征综合评估改进最近邻域关联度方法。在数据关联时,对落入波门内的点迹综合考虑位置、幅度、多普勒、工作频率四种量测指标,参考预测值或先验信息得到各指标下的评估值,运用熵权值法确定权值,得到综合评估值,利用质量评估中最优量测代替传统方法中离预测点距离最小的量测作为航迹更新。从而提高传感器信息利用率,提高在复杂环境下微弱目标跟踪性能。
申请公布号 CN104155650A 申请公布日期 2014.11.19
申请号 CN201410401795.7 申请日期 2014.08.14
申请人 电子科技大学 发明人 孔令讲;李雯雯;吴健;刘羽锐;易伟;崔国龙;李溯琪;杨建宇;杨晓波
分类号 G01S13/66(2006.01)I 主分类号 G01S13/66(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 一种基于熵权值法点迹质量评估的目标跟踪方法,该方法包括:步骤1、初始化目标状态X(k)及目标协方差矩阵P(k),k=0;目标的状态转移方程为:X(k+1)=FX(k)+ω(k)其中,F为状态转移矩阵,ω(k)为零均值的高斯过程噪声,其协方差矩阵为Q,量测模型为:Z(k)=HX(k)+v(k)Z(k)是目标量测,H为观测矩阵,v(k)为零均值的高斯观测噪声,其协方差矩阵为R;步骤2、根据k‑1时刻目标t的状态X(k‑1)及其协方差矩阵P(k‑1),得到k时刻目标t的状态预测值:X(k|k‑1)=FX(k‑1)量测预测值:Z(k|k‑1)=HX(k|k‑1)协方差一步预测值P(k|k‑1):P(k|k‑1)=FP(k‑1)F<sup>T</sup>+Q;步骤3、从雷达接收机中读取第k帧位置量测集合<img file="FDA0000554202900000011.GIF" wi="355" he="85" />z<sub>i</sub>(k)为k时刻的第i个量测,同时,并获取第k帧量测幅度集合<img file="FDA0000554202900000012.GIF" wi="362" he="85" />多普勒集合<img file="FDA0000554202900000013.GIF" wi="381" he="85" />工作频率集合<img file="FDA0000554202900000014.GIF" wi="373" he="85" />且每个位置量测值对应三个知识辅助量测指标a(k)、d(k)、f(k),计算量测值z(k)是否满足下式,若满足则作为候选回波,并统计满足要求的回波个数,[z(k)‑Z(k|k‑1)]<sup>T</sup>S<sup>‑1</sup>(k)×[z(k)‑Z(k|k‑1)]≤γ其中S(k)为k时刻的新息协方差S(k)=HP(k|k‑1)H<sup>T</sup>+R,γ为预先设定的波门大小;步骤4、根据步骤3统计得到落入波门内的候选目标个数,来更新目标航迹,若候选回波个数为零,即没有量测落入相关波门内,用步骤2的预测值更新航迹,转入步骤8.1;若落入相关波门内的量测值只有一个,则该量测值可直接用于航迹更新,转入步骤8.2;若有一个以上的回波落入被跟踪目标的相关波门内,则对波门内点迹进行质量评估,选取最优点,转入步骤5;步骤5、计算各候选回波的空间位置评估值、幅度评估值、多普勒评估值、工作频率评估值,得到各评估指标下的评估向量,再将评估向量组合为一个评估矩阵A;步骤6、计算出各项评估指标的熵权法权值,组成权值向量;步骤6.1设x<sub>ij</sub>为点迹质量矩阵A中的元素,i=1,2,3,...,m<sub>k</sub>,m<sub>k</sub>为落入波门内的有效量测总个数,j=1,2,3,4分别代表四项量测指标,然后进行数据平移x<sub>ij</sub>=x<sub>ij</sub>+1,因利用熵值求权数时要用对数进行计算,进行数据平移可避免取对数时无意义;步骤6.2计算第j项指标下第i个因子的值在此项指标中所占的比重<img file="FDA0000554202900000021.GIF" wi="255" he="199" />步骤6.3计算第j项指标的熵值<img file="FDA0000554202900000022.GIF" wi="465" he="137" />步骤6.4计算第j项指标的差异系数:<img file="FDA0000554202900000023.GIF" wi="334" he="214" />步骤6.5计算权值<img file="FDA0000554202900000024.GIF" wi="473" he="202" />得到四项评估指标权值向量;步骤7、将步骤5得到的评估矩阵与步骤6得到的权值向量相乘,得到最终质量评估值V=Aw,找到评估值V中最大值对应的位置i<sub>op</sub>,得到波门内最优量测z<sub>op</sub>(k);步骤8、计算目标t的滤波增益矩阵K(k)=P(k|k‑1)H(k)S(k)<sup>‑1</sup>,得到状态更新表达式x(k|k),进行航迹更新;步骤8.1若候选回波个数为零,即没有量测落入相关波门内,利用预测点更新航迹;x(k|k)=x(k|k‑1)P(k|k)=P(k|k‑1)步骤8.2若落入相关波门内的量测值只有一个,则该量测值可直接用于航迹更新;x(k|k)=x(k|k‑1)+K(z)[z(k)‑HX(k|k‑1)]P(k|k)=P(k|k‑1)‑K(k)HP(k|k‑1)步骤8.3若波门内有多个点,则选取经过点迹质量评估后用最优点更新航迹x(k|k)=x(k|k‑1)+K(z)[z<sub>op</sub>(k)‑Z(k|k‑1)]P(k|k)=P(k|k‑1)‑K(k)HP(k|k‑1)。
地址 611731 四川省成都市高新区(西区)西源大道2006号