发明名称 一种基于变系数α‑β滤波器的目标跟踪方法与跟踪系统
摘要 本发明适用于多传感器信息融合领域,提供了一种基于变系数α‑β滤波器的目标跟踪方法与跟踪系统。所述方法步骤如下:首先根据前一时刻的目标矩及狄拉克项,预测前一时刻已经存在的目标在当前时刻的目标矩及狄拉克项,并为当前时刻新生的目标指定相应的目标矩及狄拉克项,再结合当前时刻的位置测量并利用变系数α‑β滤波器确定更新的狄拉克项,再对更新的狄拉克项进行裁减与合并,得到当前时刻的目标矩及狄拉克项,当前时刻的目标矩及狄拉克项作为下一时刻递归的输入,最后提取权重系数大于第一阈值的狄拉克项作为当前时刻的输出,所输出狄拉克项中的目标状态值为当前时刻一个目标的状态值。
申请公布号 CN104063615B 申请公布日期 2017.02.15
申请号 CN201410315478.3 申请日期 2014.07.03
申请人 深圳大学 发明人 刘宗香;谢维信;李良群;李丽娟
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 深圳市恒申知识产权事务所(普通合伙) 44312 代理人 陈健
主权项 一种基于变系数α‑β滤波器的目标跟踪方法,其特征在于,包括以下步骤:步骤S1:根据前一时刻目标矩<img file="FDA0001113220610000011.GIF" wi="590" he="135" />及狄拉克项<img file="FDA0001113220610000012.GIF" wi="358" he="79" />预测前一时刻已经存在的目标在当前时刻的狄拉克项<img file="FDA0001113220610000013.GIF" wi="414" he="86" />其中,i=1,2,…,J<sub>k‑1</sub>,δ表示狄拉克分布,x表示目标的状态,<img file="FDA0001113220610000014.GIF" wi="84" he="63" />表示前一时刻第i个狄拉克项的权重系数,<img file="FDA0001113220610000015.GIF" wi="94" he="63" />表示前一时刻第i个狄拉克项的目标状态值,J<sub>k‑1</sub>表示前一时刻狄拉克项的数目,<img file="FDA0001113220610000016.GIF" wi="65" he="62" />表示第i个狄拉克项对应的新生狄拉克项产生的时间,<img file="FDA0001113220610000017.GIF" wi="101" he="71" />表示第i个预测的狄拉克项的权重系数,<img file="FDA0001113220610000018.GIF" wi="118" he="71" />表示第i个预测的狄拉克项的目标状态值,<img file="FDA0001113220610000019.GIF" wi="86" he="71" />表示第i个预测的狄拉克项对应的新生狄拉克项产生的时间,并且<img file="FDA00011132206100000110.GIF" wi="994" he="87" />其中,F<sub>k‑1</sub>为状态转移矩阵,p<sub>sk</sub>为目标幸存的概率;预测的目标矩p<sub>k|k‑1</sub>(x)为<img file="FDA00011132206100000111.GIF" wi="678" he="150" />以及,指定当前时刻新生目标的狄拉克项为<img file="FDA00011132206100000112.GIF" wi="478" he="87" />i=1,2,…,J<sub>γk</sub>,新生目标的目标矩为<img file="FDA00011132206100000113.GIF" wi="678" he="151" />其中,<img file="FDA00011132206100000114.GIF" wi="118" he="70" />表示新生的狄拉克项的权重系数,<img file="FDA00011132206100000115.GIF" wi="131" he="67" />表示新生的狄拉克项的目标状态值,<img file="FDA00011132206100000116.GIF" wi="102" he="70" />表示新生狄拉克项产生的时间,J<sub>γk</sub>为新生狄拉克项的个数,i为索引号,并且<img file="FDA00011132206100000117.GIF" wi="206" he="85" />t表示当前时刻的时间;步骤S2:根据预测的前一时刻已经存在的目标在当前时刻的狄拉克项<img file="FDA00011132206100000118.GIF" wi="389" he="69" />及目标矩p<sub>k|k‑1</sub>(x),以及指定的当前时刻新生目标的狄拉克项<img file="FDA0001113220610000021.GIF" wi="449" he="71" />及目标矩p<sub>γk</sub>(x),以及当前时刻的位置测量,利用变系数α‑β滤波器确定出当前时刻更新的狄拉克项<img file="FDA0001113220610000022.GIF" wi="331" he="87" />i=1,2,…,J<sub>k‑1</sub>+J<sub>γk</sub>,j=1,2,…,n<sub>k</sub>+1,其中,<img file="FDA0001113220610000023.GIF" wi="78" he="71" />为更新的狄拉克项的权重系数、<img file="FDA0001113220610000024.GIF" wi="86" he="71" />为更新的狄拉克项的目标状态值,<img file="FDA0001113220610000025.GIF" wi="68" he="71" />为更新的狄拉克项产生的时间,n<sub>k</sub>为当前时刻位置测量的个数,所述当前时刻的位置测量包含由当前时刻目标产生的位置测量和由当前时刻杂波产生的位置测量,i和j是索引号,并且当j≤n<sub>k</sub>时,<img file="FDA0001113220610000026.GIF" wi="1998" he="235" /><img file="FDA0001113220610000027.GIF" wi="230" he="88" />其中,t表示当前时刻的时间,K<sub>i</sub>为增益矩阵,z<sub>j</sub>表示当前时刻n<sub>k</sub>个位置测量中的第j个位置测量,H<sub>k</sub>为观测矩阵,R<sub>j</sub>为观测噪声的方差矩阵,c为放大因子,R<sub>0</sub>是为探测新目标而设置的协方差,p<sub>Dk</sub>为目标的检测概率,λ<sub>C</sub>为当前时刻观测空间中杂波的密度,当j=n<sub>k</sub>+1时,<img file="FDA0001113220610000028.GIF" wi="582" he="86" /><img file="FDA0001113220610000029.GIF" wi="437" he="85" />K<sub>i</sub>取变系数α‑β滤波器的增益矩阵,并且<img file="FDA00011132206100000210.GIF" wi="459" he="295" />其中,T为当前时刻与前一时刻的时间差,α<sub>k</sub>和β<sub>k</sub>是两个时变的系数,并且<img file="FDA00011132206100000211.GIF" wi="474" he="116" /><img file="FDA00011132206100000212.GIF" wi="477" he="95" />步骤S3:将当前时刻更新的狄拉克项<img file="FDA00011132206100000213.GIF" wi="333" he="85" />i=1,2,…,J<sub>k‑1</sub>+J<sub>γk</sub>,j=1,2,…,n<sub>k</sub>+1表述为<img file="FDA00011132206100000214.GIF" wi="327" he="86" />q=1,…,J<sub>k|k</sub>;其中,q为索引号,取值为1至J<sub>k|k</sub>,J<sub>k|k</sub>为更新的狄拉克项的个数,<img file="FDA00011132206100000215.GIF" wi="77" he="63" />为第q个狄拉克项的权重,<img file="FDA00011132206100000216.GIF" wi="81" he="67" />为第q个狄拉克项的目标状态值,<img file="FDA0001113220610000031.GIF" wi="54" he="63" />为第q个狄拉克项对应的新生狄拉克项产生的时间,并且J<sub>k|k</sub>=(n<sub>k</sub>+1)(J<sub>k‑1</sub>+J<sub>γk</sub>),q=(J<sub>k‑1</sub>+J<sub>γk</sub>)(j‑1)+i,<img file="FDA0001113220610000032.GIF" wi="806" he="87" />i=1,2,…,J<sub>k‑1</sub>+J<sub>γk</sub>,j=1,2,…,n<sub>k</sub>+1;以及,从当前时刻更新的狄拉克项<img file="FDA0001113220610000033.GIF" wi="323" he="79" />q=1,…,J<sub>k|k</sub>中裁减掉权重系数<img file="FDA0001113220610000034.GIF" wi="75" he="70" />小于第二阈值的狄拉克项;以及,从裁减后余下的狄拉克项中将距离d<sub>ij</sub>小于第三阈值的狄拉克项合并成一个;其中,合并距离d<sub>ij</sub>为<img file="FDA0001113220610000035.GIF" wi="712" he="86" />多个狄拉克项的合并方法为:<img file="FDA0001113220610000036.GIF" wi="734" he="143" />其中,L为合并狄拉克项上标形成的集合,上标T和b分别表示矩阵的转置和合并后狄拉克项的索引号;合并狄拉克项的时间<img file="FDA0001113220610000037.GIF" wi="40" he="62" />取合并前狄拉克项中权重系数最大狄拉克项的时间;以及,将裁减与合并后的狄拉克项<img file="FDA0001113220610000038.GIF" wi="259" he="86" />i=1,…,J<sub>k</sub>作为当前时刻的狄拉克项,并将裁减与合并后的狄拉克项的加权和作为当前时刻的目标矩<img file="FDA0001113220610000039.GIF" wi="533" he="141" />其中,J<sub>k</sub>为当前时刻狄拉克项的个数;步骤S4:根据当前时刻的目标矩及狄拉克项,提取权重系数大于第一阈值的狄拉克项作为当前时刻的输出,并将所输出狄拉克项中的目标状态值为当前时刻一个目标的状态值。
地址 518000 广东省深圳市南山区南海大道3688号