[发明专利]基于粒子滤波和航迹管理的微弱目标检测前跟踪方法有效
申请号: | 201210106282.4 | 申请日: | 2012-04-02 |
公开(公告)号: | CN102621543A | 公开(公告)日: | 2012-08-01 |
发明(设计)人: | 王国宏;谭顺成;于洪波;关成斌;贾舒宜;曹倩 | 申请(专利权)人: | 中国人民解放军海军航空工程学院 |
主分类号: | G01S13/66 | 分类号: | G01S13/66 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 264001 山*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种基于粒子滤波和航迹管理的微弱目标检测前跟踪方法,属于雷达数据处理领域。基于粒子滤波的微弱目标检测前跟踪方法不能提供目标航迹信息,且在目标信噪比较低时,容易出现漏检的现象。本发明提出的基于粒子滤波和航迹管理的微弱目标检测前跟踪方法立足于解决此类问题。本发明具有结构简单,计算快速,易于硬件实现等优点,同时克服了基于粒子滤波的微弱目标检测前跟踪方法的局限性,通过有效解决漏检问题,进而保证目标航迹的连续性,因此具有较强的工程应用价值和推广前景。 | ||
搜索关键词: | 基于 粒子 滤波 航迹 管理 微弱 目标 检测 跟踪 方法 | ||
【主权项】:
1.基于粒子滤波和航迹管理的微弱目标检测前跟踪方法,其特征包括以下步骤:步骤1:变量初始化K是总仿真时间;T是雷达扫描周期;N是滤波器采用的粒子数;SNR是目标信噪比;Imin和Imax是目标可能的最小和最大强度;vmin和vmax是目标可能的最小和最大速度;γ是目标存在判决门限;η0是目标初始存在概率;L×H表示雷达监测区域的大小;∏k是目标存在状态概率转移矩阵;Fk是目标状态转移矩阵;Qk是过程噪声协方差;步骤2:令k=0,进行粒子初始化,得到初始粒子集其中,包含目标的位置速度和强度信息,是目标存在状态变量,i=1,2,…,N(1)令i=1;(2)目标位置初始化x k l = L × rand ( 1 ) - - - ( 1 ) ]]>y k l = H × rand ( 1 ) - - - ( 2 ) ]]> rand(1)表示按照均匀分布产生一个[0,1]区间上的随机数;(3)目标速度初始化x · k l = v min + ( v max - v min ) × rand ( 1 ) - - - ( 3 ) ]]>y · k l = v min + ( v max - v min ) × rand ( 1 ) - - - ( 4 ) ]]> (4)目标强度初始化I · k l = v min + ( v max - v min ) × rand ( 1 ) - - - ( 5 ) ]]> (5)根据目标初始存在概率η0初始化目标存在状态变量;(6)粒子权重初始化w k l = 1 N - - - ( 6 ) ]]> (7)令i=i+1,若i≤N,转(2);步骤3:令k=k+1,获得k时刻的雷达量测将雷达接收到的信号进行A/D变换,得到k时刻的雷达量测z k { z k ( p , q ) | p = 1,2 , . . . , m ; q = 1,2 , . . . , n } , ]]> 送雷达数据处理计算机;步骤4:对初始粒子集进行预测,得到预测粒子集其中,(1)令i=1;(2)根据前一时刻的目标存在状态变量和目标存在状态概率转移矩阵∏k预测当前时刻目标存在状态变量(3)若且利用方程(1)-方程(5)生成粒子状态(4)若且利用状态转移方程对粒子i的状态进行一步预测x k | k - 1 l = F k x k - 1 l + v k - - - ( 7 ) ]]> 其中,vk为过程噪声,其噪声协方差为Qk;(5)令i=i+1,若i≤N,转(2);步骤5:利用当前时刻的量测对粒子权重进行更新(1)令i=1;(2)确定似然区域的范围C p ( x k | k - 1 l ) = { p 0 - s , . . . , p 0 - 1 , p 0 , p 0 + 1 , . . . , p 0 + s } - - - ( 8 ) ]]>C q ( y k | k - 1 l ) = { q 0 - s , . . . , q 0 - 1 , q 0 , q 0 + 1 , . . . , q 0 + s } - - - ( 9 ) ]]> 其中,p0表示最接近状态矢量元素的整数,q0表示最接近状态矢量元素的整数,s为预设参数;(3)粒子权重更新w k | k - 1 l = Π p ∈ C p ( x k | k - 1 l ) Π q ∈ C q ( x k | k - 1 l ) l ( z k ( p , q ) | x k | k - 1 l ) E k | k - 1 l = 1 1 E k | k - 1 l = 0 - - - ( 10 ) ]]> 其中l ( z k ( p , q ) | x k | k - 1 l ) = exp { - h k ( p , q ) ( x k | k - 1 l ) ( h k ( p , q ) ( x k | k - 1 l ) - 2 z k ( p , q ) ) 2 σ 2 } - - - ( 11 ) ]]>h k ( p , q ) ( x k | k - 1 l ) = Δ x Δ y I k | k - 1 l 2 π Σ 2 exp { - ( pΔ x - x k | k - 1 l ) 2 + ( qΔ y - y k | k - 1 l ) 2 2 Σ 2 } - - - ( 12 ) ]]> 其中,σ2为量测噪声协方差,Δx和Δy表示传感器分辨率,∑表示传感器模糊斑点数;(4)令i=i+1,若i ≤N,转(2);(5)权重归一化w k l = w k | k - 1 l Σ j = 1 N w k | k - 1 j , ]]> i=1,2,…,N (13)步骤6:重采样对权重更新后的粒子集进行系统重采样{ y k l , ( 1 / N ) } l = 1 N = resample { y k | k - 1 l , w k l } l = 1 N - - - ( 14 ) ]]> 得到k时刻的粒子集其中步骤7:目标状态和状态协方差估计(1)计算目标存在概率P ^ k = 1 N Σ l = 1 N E k l - - - ( 15 ) ]]> 若转步骤9;(2)目标状态估计x ^ k = Σ l = 1 N x k l E k l Σ l = 1 N E k l - - - ( 16 ) ]]> (3)状态协方差估计B ^ k = 1 N Σ l = 1 N ( x ^ k - x k l ) ( x ^ k - x k l ) T - - - ( 17 ) ]]> 步骤8:航迹更新(1)将目标状态估计等效为新的量测值,状态协方差估计等效为量测误差;(2)若当前时刻存在目标航迹τk-1或τ′k-1,利用新量测对其进行更新,得到更新的航迹τk,转(5);(3)若当前时刻不存在目标航迹,将新量测作为航迹头,利用3/4逻辑法进行航迹起始;(4)若航迹起始成功,记为τk,设置航迹维持标记flag=0;(5)转步骤10;步骤9:航迹维持和终结(1)若当前时刻不存在目标航迹,转步骤10;(2)若当前时刻存在目标航迹τk-1或τ′k-1,令flag=flag+1,若flag>2,转(4);(3)对τk-1或τ′k-1进行一步预测,得到τ′k,转步骤10;(4)输出航迹τk-3并删除航迹τ′k-1;步骤10:重复步骤3~步骤9,直至雷达关机。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国人民解放军海军航空工程学院,未经中国人民解放军海军航空工程学院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201210106282.4/,转载请声明来源钻瓜专利网。