发明名称 基于粒子群优化的传感器网络对机动目标的协同跟踪方法
摘要 本发明公开了一种基于粒子群优化的传感器网络对机动目标的协同跟踪方法。其实现步骤是:(1)根据目标当前的位置,采用粒子群优化算法对传感器网络中的传感器位置进行重新部署;(2)利用重新部署后的传感器对目标的预测位置进行预探测,获得关于目标位置的期望的传感器观测值集合;(3)将期望传感器观测值集合中观测值随机分成N<sub>C</sub>个期望传感器观测值组合;(4)采用二值粒子群优化算法求取上述N<sub>C</sub>个期望传感器观测值组合中的最优传感器观测值组合;(5)用最优观测值组合所对目标进行真实跟踪探测,并将跟踪结果进行融合。本发明能在非高斯非线性环境下保证目标跟踪精度,降低网络能耗,可用于传感器网络对机动目标的协同跟踪。
申请公布号 CN102967857B 申请公布日期 2014.04.16
申请号 CN201210501347.5 申请日期 2012.11.28
申请人 西安电子科技大学 发明人 刘峥;刘钦;陈熠;谢荣;刘韵佛
分类号 G01S13/66(2006.01)I 主分类号 G01S13/66(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于粒子群优化的传感器网络对机动目标的协同跟踪方法,包括如下步骤:1)根据目标当前的位置,采用粒子群优化算法对传感器网络中的动态传感器的位置进行重新部署;2)利用重新部署后的动态传感器和静态传感器对目标的预测位置进行预探测,获得关于目标位置的期望传感器观测值集合<img file="FDA0000442838300000011.GIF" wi="178" he="94" />其中,k为观测时间,j为传感器编号,N为动态传感器总数,M为静态传感器总数;3)将期望传感器观测值集合<img file="FDA0000442838300000012.GIF" wi="160" he="91" />中的N+M个观测值进行随机组合,形成N<sub>C</sub>个期望传感器观测值组合,其中,<img file="FDA0000442838300000013.GIF" wi="353" he="87" />q表示第q个传感器观测值,<img file="FDA0000442838300000014.GIF" wi="105" he="63" />是从N+M个期望观测值中选择的第q个观测值的排列数;4)采用二值粒子群优化算法求取上述N<sub>C</sub>个期望传感器观测值组合中的最优传感器观测值组合:(4a)在N<sub>C</sub>个期望传感器观测值组合中随机选择N<sub>P</sub>个传感器观测值组合,作为N<sub>P</sub>个粒子群位置向量集<img file="FDA0000442838300000015.GIF" wi="196" he="86" />其中Y<sub>i</sub>(t)为粒子群中第i个粒子在第t代的位置向量,<img file="FDA0000442838300000016.GIF" wi="553" he="83" /><img file="FDA0000442838300000017.GIF" wi="107" he="74" />为粒子i在第t代的位置向量Y<sub>i</sub>(t)的第l个分量,<img file="FDA0000442838300000018.GIF" wi="107" he="74" />代表一个传感器选择向量,<img file="FDA0000442838300000019.GIF" wi="161" he="73" />代表该位置上对应的传感器观测值被选择,<img file="FDA00004428383000000110.GIF" wi="167" he="73" />代表该位置上对应的传感器观测值未被选择;(4b)求取N<sub>P</sub>个粒子群位置向量集<img file="FDA00004428383000000111.GIF" wi="164" he="81" />的适应度函数值F:<maths num="0001"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>E</mi><mn>0</mn></msub><mo>-</mo><msub><mi>D</mi><mi>&alpha;</mi></msub></mtd><mtd><msub><mi>D</mi><mi>&alpha;</mi></msub><mo>></mo><msub><mi>D</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>S</mi></msub><mo>-</mo><msub><mi>D</mi><mn>0</mn></msub></mtd><mtd><msub><mi>D</mi><mi>&alpha;</mi></msub><mo>&le;</mo><msub><mi>D</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,E<sub>0</sub>为网络能耗门限,D<sub>α</sub>为当前选择的粒子所代表的期望观测值组合的Rényi信息增量,D<sub>0</sub>为Rényi信息增量门限,E<sub>S</sub>为当前选择的粒子所代表的期望观测值组合的Rényi信息增量;(4c)令上述N<sub>P</sub>个适应度函数值中最小值所对应的粒子为全局最优粒子,利用全局最优粒子通过如下公式对当前代粒子进行变异,获取下一代新粒子;<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>y</mi><mi>i</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>exchange</mi><mrow><mo>(</mo><msubsup><mi>y</mi><mi>i</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd><mi>rand</mi><mo>&lt;</mo><mn>2</mn><mo>&times;</mo><mo>|</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mrow></mfrac><mo>-</mo><mn>0.5</mn><mo>|</mo></mtd></mtr><mtr><mtd><msubsup><mi>y</mi><mi>i</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>rand</mi><mo>&GreaterEqual;</mo><mn>2</mn><mo>&times;</mo><mo>|</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>l</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mrow></mfrac><mo>-</mo><mn>0.5</mn><mo>|</mo></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中<img file="FDA0000442838300000022.GIF" wi="161" he="72" />为第t+1代的位置向量Y<sub>i</sub>(t+1)的第l个分量,rand表示[0,1]内均匀分布的随机数,exchange(·)表示0,1取反运算,<img file="FDA0000442838300000023.GIF" wi="107" he="74" />为粒子i在第t代的位置向量Y<sub>i</sub>(t)的第l个分量<img file="FDA0000442838300000024.GIF" wi="136" he="73" /><img file="FDA0000442838300000025.GIF" wi="154" he="73" />为粒子i在第t+1代的速度向量;(4d)返回步骤(4b),直至到达迭代次数的最大值,获得全局最优传感器观测值组合;5)用最优观测值组合所对应的传感器组合对目标进行真实跟踪探测,并将跟踪结果进行融合,得到k+1时刻目标状态估计的融合值x<sub>k+1/k+1</sub>:<maths num="0003"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>{</mo><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>x</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>x</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msubsup><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>x</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mrow><mo>{</mo><msubsup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>}</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mtd></mtr></mtable></mfenced></math>]]></maths>其中,P<sub>k+1/k+1</sub>为k+1时刻的误差协方差矩阵,x<sub>k+1/k</sub>为系统在k时刻对k+1时刻目标状态的预测值,P<sub>k+1/k</sub>为系统在k时刻对k+1时刻目标状态的误差协方差矩阵,N<sub>j</sub>为最优观测值组合中所包含的传感器个数,x<sub>j,k+1/k+1</sub>为第j个传感器的估计值,P<sub>j,k+1/k+1</sub>为第j个传感器的误差协方差矩阵,[]<sup>-1</sup>为矩阵求逆。
地址 710071 陕西省西安市太白南路2号