发明名称 水下多目标跟踪方法
摘要 本发明公开了一种水下多目标跟踪方法,无线传感器网络定位系统包括多个传感器节点、簇头节点以及控制系统,系统定位采用一种改进重采样非扩维无味粒子滤波算法实现在基于传感器网络配置下的目标定位,该算法充分利用当前时刻方位测量的信息来估计粒子后验概率密度函数,克服了PF算法中粒子选择具有盲目性的缺点,在很少粒子数目的情况下就能达到比PF算法近似且更好的定位精度,并且不需要扩展向量中状态维数,这样大大降低了计算复杂度,提高了计算时间,有利于系统的实时实现。
申请公布号 CN103645487B 申请公布日期 2017.04.05
申请号 CN201310655173.2 申请日期 2013.12.06
申请人 江苏科技大学 发明人 王彪;冯科磊;朱志宇
分类号 G01S19/42(2010.01)I 主分类号 G01S19/42(2010.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 楼高潮
主权项 一种水下多目标跟踪方法,基于无线传感器网络定位系统,所述无线传感器网络定位系统包括多个传感器节点(1)、簇头节点(2)以及控制系统(3),其特征在于,水下多目标跟踪方法包括以下步骤:1)无线传感器网络定位系统中的传感器节点(1)自我编号,标示唯一身份,所述传感器节点(1)以自带的GPS装置完成时间同步与自身定位;2)依据传感器节点(1)通信距离对网络进行簇划分,并且选择簇头节点(2),完成网络初始化;3)所述传感器节点(1)对所在簇内进行监测,依据采集信号判断是否有目标出现,当出现目标后,利用阵列信号处理技术对目标的方位进行测量;4)传感器节点(1)将测量到的目标相对于传感器节点(1)的方位角及自身所在位置经纬度信息通过无线方式发送给簇头节点(2);5)根据k‑1时刻的估计值<img file="FDA0001121613420000011.GIF" wi="174" he="70" />预测k时刻的测量值<img file="FDA0001121613420000012.GIF" wi="253" he="70" />其中<img file="FDA0001121613420000013.GIF" wi="1302" he="78" />x<sub>p</sub>(k‑1),x<sub>v</sub>(k‑1)分别表示目标在x轴方向的位置和速度,y<sub>p</sub>(k‑1),y<sub>v</sub>(k‑1)分别表示目标在y轴方向的位置和速度,z<sub>p</sub>(k‑1),z<sub>v</sub>(k‑1)分别表示目标在z轴方向的位置和速度;<img file="FDA0001121613420000014.GIF" wi="213" he="70" />表示k时刻传感器节点(1)测量到的方位角,包含两个角度即:水平方位角<img file="FDA0001121613420000015.GIF" wi="341" he="127" />以及俯仰角<img file="FDA0001121613420000016.GIF" wi="734" he="143" />其中,x(k),y(k),z(k)表示目标在k时刻的位置坐标,x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>表示1号传感器节点所部署的位置坐标;6)拆分得到与目标t所有可能的互联集合L,其中t表示k时刻的目标,t=1,2,…,T;T表示目标总数;7)对于所有可能的互联集合L,for i=1到N<sub>s</sub>,N<sub>s</sub>表示系统中传感器的个数;①根据预测<img file="FDA0001121613420000017.GIF" wi="214" he="71" />确定目标t确认区域,产生确认矩阵Ω=[w<sub>jt</sub>],j=1,2,…,m<sub>k</sub>;t=0,1,…,T,其中m<sub>k</sub>表示k时刻落入跟踪波门中的有效回波量测数目,w<sub>jt</sub>表示量测j落入目标t的确认门内,当量测j落入目标t的跟踪门时w<sub>jt</sub>=1,否则w<sub>jt</sub>=0;确认矩阵用来表示有效回波和各目标跟踪门的关系;②把确认矩阵拆成互联矩阵;③根据式<img file="FDA0001121613420000021.GIF" wi="1499" he="136" />计算互联概率<img file="FDA0001121613420000022.GIF" wi="163" he="71" />其中,Pr()表示括号内所指事件发生的概率;θ<sub>i</sub>(k)是第i个联合互联事件,i=1,2,…,n<sub>k</sub>,表示m<sub>k</sub>个量测与各个目标相匹配的一种可能;<img file="FDA0001121613420000023.GIF" wi="118" he="71" />具体表示了此联合事件中θ<sub>i</sub>(k)量测j来源于目标t;θ<sub>jt</sub>(k)指第j个量测与第t个目标互联的事件;<img file="FDA0001121613420000024.GIF" wi="62" he="80" />为目标检测指示,代表了联合互联事件θ<sub>i</sub>(k)中目标是否被检测到;<img file="FDA0001121613420000025.GIF" wi="123" he="63" />表示给定测量序列Z<sup>k</sup>的条件下事件<img file="FDA0001121613420000026.GIF" wi="110" he="63" />的概率,<img file="FDA0001121613420000027.GIF" wi="142" he="71" />是第i个传感器的第l<sub>i</sub>个量测与目标t互联的单个传感器联合概率数据互联事件的概率;8)按照改进重采样非扩维UPF算法根据L计算<img file="FDA0001121613420000028.GIF" wi="162" he="78" />所述<img file="FDA0001121613420000029.GIF" wi="131" he="79" />是k时刻第t个目标的改进重采样非扩维UPF的状态估计;所述改进重采样非扩维UPF算法包括:①采用非扩维UKF产生粒子建议分布;②将需要复制的粒子和被抛弃粒子进行组合,生成新的粒子,使新粒子同时包含被复制粒子和被抛弃粒子的信息,保证了粒子包含足够的信息来表示目标状态在状态空间的分布;9)根据式<img file="FDA00011216134200000210.GIF" wi="910" he="137" />计算<img file="FDA00011216134200000211.GIF" wi="155" he="71" />得到k时刻目标状态估计,式中<img file="FDA00011216134200000212.GIF" wi="131" he="71" />是k时刻第t个目标的改进重采样非扩维UPF的状态估计,<img file="FDA00011216134200000213.GIF" wi="141" he="70" />是第i个传感器的第l<sub>i</sub>个量测与目标t互联的单个传感器联合概率数据互联事件的概率,<img file="FDA00011216134200000214.GIF" wi="127" he="71" />指由多传感器的联合事件关联概率<img file="FDA00011216134200000215.GIF" wi="147" he="62" />和对k时刻第t个目标的改进重采样非扩维UPF的状态估计,得出的最终多传感器系统中第t个目标的状态估计;10)跳转到步骤5)进行下一时刻的循环。
地址 212003 江苏省镇江市梦溪路2号