发明名称 基于无线传感器网络的无源多目标检测跟踪方法
摘要 本发明涉及一种基于无线传感器网络的无源多目标检测跟踪方法,其主要技术特点是:包括以下步骤:根据无线传感器网络中不同无线链路接收信号强度,采用扫描圈开窗检测方法构建无源多目标的滑动扫描圈模型,并采用HAC聚类算法进行聚类分析从而提取出检测结果;根据检测结果并采用PHD无源多目标粒子滤波跟踪算法,得到目标数目变化的多目标跟踪结果。本发明设计合理,其检测与跟踪算法具有较高的精度和鲁棒性,能在复杂的多径环境中检测与跟踪多个目标,同时目标检测与跟踪算法的计算复杂度适中,保证检测跟踪系统运行的实时性。
申请公布号 CN103298156B 申请公布日期 2016.01.20
申请号 CN201310233976.9 申请日期 2013.06.13
申请人 北京空间飞行器总体设计部;北京邮电大学 发明人 刘崇华;杨波;张弓;薛建飞;黎杨;刘航;门爱东
分类号 H04W84/18(2009.01)I;H04W64/00(2009.01)I 主分类号 H04W84/18(2009.01)I
代理机构 天津盛理知识产权代理有限公司 12209 代理人 王利文
主权项 一种基于无线传感器网络的无源多目标检测跟踪方法,其特征在于:包括以下步骤:步骤1:根据无线传感器网络中不同无线链路接收信号强度,采用扫描圈开窗检测方法构建无源多目标的滑动扫描圈模型,并采用HAC聚类算法进行聚类分析从而提取出检测结果;步骤2:根据检测结果并采用PHD无源多目标粒子滤波跟踪算法,得到目标数目变化的多目标跟踪结果;所述构建滑动扫描圈模型的方法包括以下步骤:⑴将监控区域中目标的模型近似为圆形扫描圈,使扫描圈能够检测到监控区域的所有位置;⑵统计在每一时刻每个扫描圈中被穿过的链路衰减情况,得到扫描圈中总共有n<sub>l</sub>条被穿过链路、n<sub>a</sub>条衰减链路;⑶通过以下检测似然概率p<sub>c</sub>(i,j)判断扫描圈中是否有目标存在:<img file="FDA0000754187610000011.GIF" wi="267" he="129" />对于第i行,第j列的扫描圈,若检测似然概率p<sub>c</sub>(i,j)大于设定的阈值0.7,则判定为此扫描圈为阳性扫描圈,反之则为阴性扫描圈,其中阳性扫描圈聚集的地方便有可能是目标存在的地方;所述PHD无源多目标粒子滤波跟踪算法包括以下步骤:(1)粒子的初始化以及预测处理;(2)粒子权重的更新处理;(3)目标估算数目的计算以及粒子的重采样处理;所述步骤(1)粒子初始化以及预测处理包括以下步骤:(1)粒子初始化使用高斯AR‑1采样函数和正态采样函数来分别处理上一时刻存活下来的旧粒子以及当前时刻的新生粒子;(2)对于预测阶段,基于k‑1时刻粒子状态<img file="FDA0000754187610000012.GIF" wi="84" he="79" />以及k时刻预测估计的粒子状态<img file="FDA0000754187610000013.GIF" wi="68" he="79" />和检测结果Z<sub>k</sub>,粒子预测权重<img file="FDA0000754187610000014.GIF" wi="112" he="85" />计算过程按以下公式进行:<img file="FDA0000754187610000021.GIF" wi="1104" he="319" />其中,β<sub>k|k‑1</sub>(x,h)=e<sub>k|k‑1</sub>(h)·f<sub>k|k‑1</sub>(x|h),e<sub>k|k‑1</sub>表示目标从k‑1时刻到k时刻目标存活的概率,f<sub>k|k‑1</sub>为高斯AR‑1过程的转移概率密度,<img file="FDA0000754187610000022.GIF" wi="91" he="79" />为k‑1时刻的粒子权重,γ<sub>k</sub>为一个新目标出生的强度函数,其服从泊松点过程,对于M<sub>k‑1</sub>个旧粒子而言,<img file="FDA0000754187610000023.GIF" wi="334" he="78" />为其单个粒子的采样概率,其服从高斯AR‑1分布,而对于J<sub>k</sub>个新生粒子而言,<img file="FDA0000754187610000024.GIF" wi="244" he="78" />为其单个粒子的采样概率,其服从正态分布;所述步骤(2)粒子权重<img file="FDA0000754187610000025.GIF" wi="82" he="78" />的更新,基于k时刻预测估计的粒子状态<img file="FDA0000754187610000026.GIF" wi="74" he="78" />和检测结果Z<sub>k</sub>,按以下公式进行:<img file="FDA0000754187610000027.GIF" wi="1042" he="186" />其中,P<sub>m</sub>为目标没被检测到的概率,c<sub>k</sub>(z)为杂波出现的概率密度函数,<img file="FDA0000754187610000028.GIF" wi="214" he="78" />为每个粒子存在似然度,其表示某一检测目标z由粒子状态<img file="FDA0000754187610000029.GIF" wi="66" he="78" />演变生成的概率,其服从高斯分布,<img file="FDA00007541876100000210.GIF" wi="116" he="85" />为粒子预测权重,W<sub>k</sub>(z)为权重更新因子,其定义如下:<img file="FDA00007541876100000211.GIF" wi="966" he="97" />所述步骤(3)目标估算数目的计算以及粒子的重采样处理方法包括以下步骤:(1)当前时刻所有粒子的权重更新完毕之后,将其求和并向上取整得到估算目标数值<img file="FDA00007541876100000212.GIF" wi="114" he="88" />(2)粒子通过分层重采样机制以粒子权重大小为优先条件将现有的M<sub>k‑1</sub>+J<sub>k</sub>个粒子重采样获得M<sub>k</sub>个粒子,其M<sub>k</sub>的数值是通过<img file="FDA00007541876100000213.GIF" wi="254" he="88" />原则计算取得,其中d<sub>c</sub>为每个目标固定分配的粒子数。
地址 100094 北京市海淀区友谊路104号