发明名称 一种基于随机集理论的分布式融合方法
摘要 该发明公开了一种随机集理论下分布式融合方法,该方法属于多传感器融合领域,它特别涉及了随机集理论下多目标跟踪和分布式多传感器融合技术领域。通过两步近似推导出闭合解形式的多目标伯努利分布式融合的表达式,为实现多目标伯努利分布式融合提供了先驱条件;近似步骤8多目标伯努利的近似可以实现多传感器网络中传感器个数大于2的分布式融合和带反馈工作模式的分布式融合下。本发明的优点在于提供了多目标伯努利分布式融合闭合解形式表达式的先驱条件,并可以在传感器个数大于2的多传感器网络中和带反馈工作模式的分布式融合下实现分布式融合。
申请公布号 CN104680002A 申请公布日期 2015.06.03
申请号 CN201510070611.8 申请日期 2015.02.10
申请人 电子科技大学 发明人 易伟;王佰录;李溯琪;杨亚;崔国龙;孔令讲;杨晓波
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 一种随机集理论下分布式融合方法,它包括以下步骤:步骤1、选定多传感器融合准则:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&pi;</mi><mi>&omega;</mi></msub><mrow><mo>(</mo><msup><mi>X</mi><mi>k</mi></msup><mo>|</mo><msubsup><mi>Z</mi><mn>1</mn><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>Z</mi><mn>2</mn><mi>k</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&pi;</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>k</mi></msup><mo>|</mo><msubsup><mi>Z</mi><mn>1</mn><mi>k</mi></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>1</mn></msub></msup><msub><mi>&pi;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>k</mi></msup><mo>|</mo><msubsup><mi>Z</mi><mn>2</mn><mi>k</mi></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>2</mn></msub></msup></mrow><mrow><msub><mrow><mo>&Integral;</mo><mi>&pi;</mi></mrow><mn>1</mn></msub><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>k</mi></msup><mo>|</mo><msubsup><mi>Z</mi><mn>1</mn><mi>k</mi></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>1</mn></msub></msup><msub><mi>&pi;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>k</mi></msup><mo>|</mo><msubsup><mi>Z</mi><mn>2</mn><mi>k</mi></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>2</mn></msub></msup><mi>&delta;</mi><msup><mi>X</mi><mi>k</mi></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000670183170000011.GIF" wi="1028" he="170" /></maths>其中,<img file="FDA0000670183170000012.GIF" wi="246" he="84" />表示第s(s=1,2)的个传感器k时刻的后验概率分布;<img file="FDA0000670183170000013.GIF" wi="254" he="82" />表示融合后的后验概率密度分布;X表示目标状态集合X={x<sub>1</sub>,…,x<sub>n</sub>},x<sub>n</sub>表示第n个目标的状态;Z表示传感器的量测集合;ω<sub>s</sub>表示该融合准则的参数,满足0≤ω<sub>s</sub>≤1,ω<sub>1</sub>+ω<sub>2</sub>=1,这个参数决定了其相应后验合分布在融合时的权重,δ表示集合变量的微分符号;步骤2、各传感器接收回波信号,并采用多目标伯努利滤波器进行本地滤波,由此各传感器得到的本地后验概率密度分布都为多目标伯努利分布;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&pi;</mi><mi>s</mi></msub><mrow><mo>(</mo><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>}</mo><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><msub><mi>i</mi><mi>s</mi></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>s</mi></msub></munderover><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msub><mi>i</mi><mi>s</mi></msub></msubsup><munder><mi>&Sigma;</mi><mrow><mn>1</mn><mo>&le;</mo><msubsup><mi>i</mi><mi>s</mi><mn>1</mn></msubsup><mo>&NotEqual;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&NotEqual;</mo><msubsup><mi>i</mi><mi>s</mi><mi>n</mi></msubsup><mo>&le;</mo><msub><mi>M</mi><mi>s</mi></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mfrac><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup><mrow><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup></mrow></mfrac><mo>,</mo><mi>s</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>S</mi></mrow>]]></math><img file="FDA0000670183170000014.GIF" wi="1253" he="173" /></maths>其中,M<sub>s</sub>为第s个传感器伯努利分量个数;<img file="FDA0000670183170000015.GIF" wi="60" he="77" />表示第s个传感器第i伯努利分量的存在概率,<img file="FDA0000670183170000016.GIF" wi="139" he="78" />为其相应的概率密度函数;S为传感器个数。步骤3、定义<img file="FDA0000670183170000017.GIF" wi="216" he="113" />为分数阶指数次幂和<img file="FDA0000670183170000018.GIF" wi="172" he="90" />称为实数的分数阶指数次幂的求和;设目标间的状态是分离的,将各传感器的多目标伯努利分布的分数阶指数次幂的形式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&pi;</mi><mi>s</mi></msub><msup><mrow><mo>(</mo><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>}</mo><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>s</mi></msub></msup><mo>=</mo><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><msub><mi>i</mi><mi>s</mi></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>s</mi></msub></munderover><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msub><mi>i</mi><mi>s</mi></msub></msubsup><munder><mi>&Sigma;</mi><mrow><mn>1</mn><mo>&le;</mo><msubsup><mi>i</mi><mi>s</mi><mn>1</mn></msubsup><mo>&NotEqual;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&NotEqual;</mo><msubsup><mi>i</mi><mi>s</mi><mi>n</mi></msubsup><mo>&le;</mo><msub><mi>M</mi><mi>s</mi></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mi>s</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mfrac><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup><mrow><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>s</mi></msub></msup></mrow>]]></math><img file="FDA0000670183170000019.GIF" wi="1135" he="205" /></maths>化简得到实数的分数阶指数次幂的求和的形式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&pi;</mi><mi>s</mi></msub><msup><mrow><mo>(</mo><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>}</mo><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>s</mi></msub></msup><mo>&ap;</mo><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><msub><mi>i</mi><mi>s</mi></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>s</mi></msub></munderover><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msub><mi>i</mi><mi>s</mi></msub></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>s</mi></msub></msup><munder><mi>&Sigma;</mi><mrow><mn>1</mn><mo>&le;</mo><msubsup><mi>i</mi><mi>s</mi><mn>1</mn></msubsup><mo>&NotEqual;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&NotEqual;</mo><msubsup><mi>i</mi><mi>s</mi><mi>n</mi></msubsup><mo>&le;</mo><msub><mi>M</mi><mi>s</mi></msub></mrow></munder><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mi>s</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mfrac><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup><mrow><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mi>s</mi><msubsup><mi>i</mi><mi>s</mi><mi>j</mi></msubsup></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>s</mi></msub></msup></mrow>]]></math><img file="FDA00006701831700000110.GIF" wi="1228" he="188" /></maths>步骤4、获得多目标伯努利分布的广义协方差交叉信息融合闭合表达式的闭合表达式;4.1、将广义协方差交叉信息融合表达式的分母项定义为一常数K;4.2、将步骤2得到的传感器1和传感器2的多目标伯努利分布式带入广义协方差交叉信息融合表达式的分子项,获得广义协方差交叉信息融合表达式分子项的闭合表达式;<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>&pi;</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>}</mo><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>1</mn></msub></msup><msub><mi>&pi;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>}</mo><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>2</mn></msub></msup></mtd></mtr><mtr><mtd><mo>&ap;</mo><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><msub><mi>i</mi><mn>1</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mn>1</mn><msub><mi>i</mi><mn>1</mn></msub></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>1</mn></msub></msup><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><msub><mi>i</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mn>2</mn><msub><mi>i</mi><mn>2</mn></msub></msubsup><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>2</mn></msub></msup><munder><mi>&Sigma;</mi><mrow><mn>1</mn><mo>&le;</mo><msubsup><mi>i</mi><mn>1</mn><mn>1</mn></msubsup><mo>&NotEqual;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&NotEqual;</mo><msubsup><mi>i</mi><mn>1</mn><mi>n</mi></msubsup><mo>&le;</mo><msub><mi>M</mi><mn>1</mn></msub></mrow></munder><munder><mi>&Sigma;</mi><mrow><mn>1</mn><mo>&le;</mo><msubsup><mi>i</mi><mn>2</mn><mn>1</mn></msubsup><mo>&NotEqual;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&NotEqual;</mo><msubsup><mi>i</mi><mn>2</mn><mi>n</mi></msubsup><mo>&le;</mo><msub><mi>M</mi><mn>2</mn></msub></mrow></munder><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mn>1</mn><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mfrac><msubsup><mi>r</mi><mn>1</mn><msubsup><mi>i</mi><mn>1</mn><mi>j</mi></msubsup></msubsup><mrow><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mn>1</mn><msubsup><mi>i</mi><mn>1</mn><mi>j</mi></msubsup></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>1</mn></msub></msup><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mn>2</mn><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mfrac><msubsup><mi>r</mi><mn>2</mn><msubsup><mi>i</mi><mn>2</mn><mi>j</mi></msubsup></msubsup><mrow><mn>1</mn><mo>-</mo><msubsup><mi>r</mi><mn>2</mn><msubsup><mi>i</mi><mn>2</mn><mi>j</mi></msubsup></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>&omega;</mi><mn>2</mn></msub></msup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000670183170000021.GIF" wi="1754" he="271" /></maths>步骤5、建立传感器1到传感器2的假设航迹映射关系集合;步骤6、利用步骤5建立的假设航迹映射集合,将步骤4得到的闭合表达式转换成非标号版本的δ‑广义多目标伯努利分布;步骤7、通过对非标号版本的δ‑广义伯努利分布进行集合积分,利用积分为1的特点,获得常数项K的闭合形式表达式,并将其代入步骤6的非标号版本的δ‑广义伯努利分布;步骤8、利用多目标伯努利分布与δ‑广义伯努利分布的一阶统计特性相匹配特点,将步骤6得到的δ‑广义伯努利分布近似为多目标伯努利分布,该分布为传感器1和传感器2的融合多目标伯努利分布;步骤9、采用与步骤4~步骤8相同的方法将传感器1和传感器2的融合多目标伯努利分布与传感器3的多目标伯努利分布进行融合;按照该方法进一步融合后序所有传感器的多目标伯努利分布;步骤10、利用序列蒙特卡洛方法实现步骤9得到的多目标伯努利分布分布式融合算法。通过上面的步骤,就可以得到基于广义协方差交叉信息融合准则的多目标伯努利分布式融合的闭合表达式,并可以实现多目标伯努利分布分布式融合。
地址 611731 四川省成都市高新区(西区)西源大道2006号