发明名称 一种基于CKF的非线性异步多传感器信息融合方法
摘要 本发明公开了一种基于CKF的非线性异步多传感器信息融合方法。分别对各个传感器分别利用CKF估计出各自的状态变量,然后采用细分时间片方法将信息融合中心的时间间隔设定为各传感器间最高精度时间单位,在相应时刻对异步多传感器的估计结果进行判断和融合,得到更加精确的状态变量估计结果。本发明可以增强对异步多传感器信息的利用率,大幅提高多传感器系统中状态变量的估计精度,增强系统的生存能力。
申请公布号 CN103455675A 申请公布日期 2013.12.18
申请号 CN201310397651.4 申请日期 2013.09.04
申请人 哈尔滨工程大学 发明人 高伟;张亚;孙骞;石惠文;周广涛;徐博;鲍桂清;史宏洋;阮双双;赵维珩
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 1.一种基于CKF的非线性异步多传感器信息融合新方法,其特征是:(1)对多传感器系统进行预热,并采集各个传感器的量测信息;(2)建立非线性异步多传感器系统的非线性状态方程和量测方程;(3)采用细分时间片方法,将融合中心的采样间隔设定为各传感器间最高精度时间单位,即取能同时被各传感器采样间隔整除的最大数为融合中心的采样间隔;设任意两传感器的采样时间分别为T<sub>i</sub>,T<sub>j</sub>(i,j=1,2,…,N,且i≠j),N为传感器的总个数;且有:<img file="FDA0000376951360000011.GIF" wi="448" he="129" />则融合中心的采样间隔为:<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;T</mi><mo>=</mo><mi>min</mi><mo>{</mo><msub><mi>&Delta;T</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>T</mi><mi>i</mi></msub><msub><mi>N</mi><mi>i</mi></msub></mfrac><mo>|</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>N</mi><mo>}</mo></mrow></math>]]></maths>(4)根据异步多传感器系统融合中心的采样时刻,在融合中心对各传感器的估计状态进行判断融合:a)k时刻,仅有一个传感器S有量测信息时,利用CKF对该传感器进行状态估计,得到该传感器的状态变量估计值<img file="FDA0000376951360000013.GIF" wi="119" he="64" />和协方差阵P<sub>S</sub>(k),则系统的最优状态估计即为该传感器的状态变量估计值:<maths num="0002"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>S</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>b)k时刻,有m(m∈{2,3,…,N},其中N为传感器的个数)个传感器S<sub>1</sub>,S<sub>2</sub>,…,S<sub>m</sub>有量测信息时,利用CKF对这些传感器分别进行滤波估计,得到各传感器的状态变量估计值<img file="FDA0000376951360000015.GIF" wi="460" he="68" />和协方差阵P<sub>1</sub>(k),P<sub>2</sub>(k),…,P<sub>m</sub>(k);然后对各个传感器的估计值进行加权处理,融合后系统的状态变量估计值为:<maths num="0003"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>D</mi><mn>1</mn></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>D</mi><mn>2</mn></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>D</mi><mi>m</mi></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,D<sub>1</sub>,D<sub>2</sub>,…,D<sub>m</sub>为相应传感器状态变量估计的权值。D<sub>1</sub>,D<sub>2</sub>,…,D<sub>m</sub>可由式(3)~(5)获得:<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>D</mi><mn>1</mn></msub><mo>=</mo><msub><mi>S</mi><mn>1</mn></msub><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>D</mi><mn>2</mn></msub><mo>=</mo><msub><mi>S</mi><mn>2</mn></msub><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>D</mi><mi>m</mi></msub><mo>=</mo><msub><mi>S</mi><mi>m</mi></msub><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>S</mi><mn>1</mn></msub><mo>=</mo><msup><mrow><mo>(</mo><mi>I</mi><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><msubsup><mi>A</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><msubsup><mi>A</mi><mi>m</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><msub><mi>S</mi><mn>1</mn></msub><msub><mi>A</mi><mn>1</mn></msub><msubsup><mi>A</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>m</mi></msub><mo>=</mo><msub><mi>S</mi><mn>1</mn></msub><msub><mi>A</mi><mn>1</mn></msub><msubsup><mi>A</mi><mi>m</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>A</mi><mn>1</mn></msub><mo>=</mo><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>P</mi><mn>1</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>A</mi><mn>2</mn></msub><mo>=</mo><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>P</mi><mn>2</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>A</mi><mi>m</mi></msub><mo>=</mo><msub><mi>P</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>P</mi><mi>m</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,Φ(k+1,k)为异步多传感器系统中k时刻的状态转移矩阵,P<sub>m</sub>(k)为异步多传感器系统中第m个传感器在k时刻的方差阵的估计值;c)k时刻,各传感器均没有量测信息时,便利用前一时刻的估计值进行时间更新;<maths num="0007"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号
您可能感兴趣的专利