发明名称 基于无线传感器网络接收信号强度的无源目标追踪方法
摘要 本发明涉及一种基于无线传感器网络接收信号强度的无源目标追踪方法,其技术特点是:包括以下步骤:无线传感器网络中的节点进行数据广播,所有节点接听并记录数据;节点计算接收信号强度值与静态网络数据的差值的绝对值,通过无源目标追踪模型进行目标位置似然度计算;用粒子滤波算法和参数估计算法并结合目标位置似然度,进行目标的实时追踪。本发明充分利用无线传感器网络所发出的射频信号具有一定的穿透性的特点,少量障碍物不会阻挡节点之间的相互通信,且无线信号的传播不受光照影响,可广泛应用于夜晚或是或是仓库、地下室等缺乏光照的室内场景。
申请公布号 CN102938876B 申请公布日期 2015.01.07
申请号 CN201210520764.4 申请日期 2012.12.07
申请人 李云鹏;陈曦 发明人 李云鹏;陈曦;马克科茨
分类号 H04W4/02(2009.01)I;H04W4/04(2009.01)I;H04W64/00(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W4/02(2009.01)I
代理机构 天津盛理知识产权代理有限公司 12209 代理人 王利文
主权项 一种基于无线传感器网络接收信号强度的无源目标追踪方法,其特征在于:包括以下步骤: 步骤1:无线传感器网络中的节点进行数据广播,所有节点接听并记录数据; 步骤2:节点计算接收信号强度值与静态网络数据的差值的绝对值,通过无源目标追踪模型进行目标位置似然度计算; 步骤3:用粒子滤波算法和参数估计算法并结合目标位置似然度,进行目标的实时追踪;该步骤包括以下处理过程: ⑴初始化:赋予参数θ<sub>0</sub>到任意初始值,将计数器b设为0,从状态的先验概率分布p(x<sub>0</sub>)抽样出状态粒子<img file="FDA0000570167690000011.GIF" wi="54" he="71" />;⑵设置计数器k=1,2…,N,N为粒子数; ⑶执行状态更新: <img file="FDA0000570167690000012.GIF" wi="759" he="93" />f为系统模型,v<sub>k</sub>为第k步的系统噪声,h为测量模型,s<sub>k</sub>为第k步的测量噪声;⑷粒子滤波,计算状态估计粒子<img file="FDA0000570167690000013.GIF" wi="84" he="84" />和权重<img file="FDA0000570167690000014.GIF" wi="449" he="86" />其中θ<sub>b‑1</sub>为第b‑1步的参数值,<img file="FDA0000570167690000015.GIF" wi="88" he="65" />()为第b‑1步的状态的后验概率公式,Z<sub>k</sub>为第k步的实际观测值,归一化<img file="FDA0000570167690000016.GIF" wi="310" he="211" />再采样<img file="FDA0000570167690000017.GIF" wi="240" he="78" />以获得新的均匀权重<img file="FDA0000570167690000018.GIF" wi="212" he="128" />从<img file="FDA0000570167690000019.GIF" wi="300" he="78" />中采样出<img file="FDA00005701676900000110.GIF" wi="126" he="78" />更新后权重为<img file="FDA00005701676900000111.GIF" wi="450" he="148" />归一化<img file="FDA00005701676900000112.GIF" wi="334" he="211" />⑸执行线上期望最大化步骤(online EM算法): <img file="FDA0000570167690000021.GIF" wi="873" he="140" />式中:Ψ为充分统计量,Ω为充分统计量的集合,充分统计量权重<img file="FDA0000570167690000026.GIF" wi="89" he="55" />满足<img file="FDA0000570167690000022.GIF" wi="463" he="168" /><img file="FDA0000570167690000023.GIF" wi="176" he="132" />和<img file="FDA0000570167690000024.GIF" wi="264" he="136" />⑹设置变量count=k mod L,L为线上期望最大化步骤的周期长度;如果count=0,则执行最大化处理:<img file="FDA0000570167690000025.GIF" wi="274" he="85" />计数器b加1;⑺循环执行步骤⑵至⑹直至所有粒子均处理结束。 
地址 610066 四川省成都市锦江区喜树街780号锦江城市花园一期8栋一单元2105