发明名称 基于MEMS的小波域多传感器信息融合系统及融合方法
摘要 本发明公开的一种基于MEMS的小波域多传感器信息融合系统及融合方法,包括依次连接的微机电陀螺阵列模块、小波域多尺度分析模块、多尺度多传感器信息融合模块、小波多尺度重构模块及组合微机电惯性测量模块,微机电陀螺阵列模块还通过数据统计性分析模块与多尺度多传感器信息融合模块相连接。通过对各个微机电陀螺的信息进行多尺度小波分解和统计分析,然后在多尺度上进行融合,最后将融合后各尺度上的微机电陀螺信息重构到原始尺度上。使微机电测量设备的精度有效提高,各微机电陀螺仪的协调性得以应用,同时也增强了该系统的可靠性。融合系统性能稳定、工作可靠、体积小、性价比高,可以为各种设备提供精确的角度测量。
申请公布号 CN102023010B 申请公布日期 2012.06.13
申请号 CN201010522825.1 申请日期 2010.10.26
申请人 西安理工大学 发明人 柯熙政;任亚飞
分类号 G01C21/16(2006.01)I;G01C21/18(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 罗笛
主权项 1.一种基于MEMS的小波域多传感器信息融合方法,其特征在于,采用基于MEMS的小波域多传感器信息融合系统,其结构为:包括依次连接的微机电陀螺阵列模块(1)、小波域多尺度分析模块(2)、多尺度多传感器信息融合模块(4)、小波多尺度重构模块(5)及组合微机电惯性测量模块(6),所述的微机电陀螺阵列模块(1)还通过数据统计性分析模块(3)与多尺度多传感器信息融合模块(4)相连接,具体按照以下步骤实施:步骤1:微机电陀螺阵列模块(1)产生各个微机电陀螺仪的信号,传递给小波域多尺度分析模块(2)及数据统计性分析模块(3),一方面,小波域多尺度分析模块(2)对各个微机电陀螺仪的信号进行小波多尺度分析,得到第i+1-N尺度上的细节信号和第N层尺度上的近似信号,传递给多尺度多传感器信息融合模块(4);另一方面,数据统计性分析模块(3)对收集到的各个微机电陀螺仪的信号进行统计分析,计算得到均值和标准差信号,传递给多尺度多传感器信息融合模块(4);所述的小波域多尺度分析模块(2)对各个微机电陀螺仪的信号进行小波多尺度分析,具体按照以下步骤实施:在尺度i+1上,对于输入陀螺仪的信号序列x(i+1,k)∈R<sup>n×1</sup>,(k∈Z),其离散小波变换的分析形势和综合形式分别为:<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>l</mi></munder><mi>h</mi><mrow><mo>(</mo><mn>2</mn><mi>k</mi><mo>-</mo><mi>l</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>d</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>l</mi></munder><mi>g</mi><mrow><mo>(</mo><mn>2</mn><mi>k</mi><mo>-</mo><mi>l</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,h(2k-l)为小波变换中的尺度系数,g(2k-l)为小波变换中的小波系数,选择合适的小波基函数,得到相应的尺度系数和小波系数,由上式得到尺度i上陀螺仪信息中的近似信号x(i,k)和细节信号d(i,k),继续这个过程对尺度i上陀螺仪的近似信号x(i,k)进行离散小波分解,得到尺度i-1上陀螺仪信息中的近似信号x(i-1,k)和细节信号d(i-1,k),重复到最优的分解尺度N上,得到小波离散分解的近似信号x(N,k)和细节信号d(N,k);然后将多尺度分析的结果:第(i+1-N)尺度上的细节信号和第N层尺度上的近似信号送到多尺度多传感器信息融合模块4等待处理;步骤2:多尺度多传感器信息融合模块(4)将收集到的由小波域多尺度分析模块(2)输入的第i+1-N尺度上的细节信号和第N层尺度上的近似信号,以及数据统计性分析模块(3)输入的均值和标准差信号,进行多尺度多传感器信息融合,得到融合后的精确多尺度上的微机电陀螺仪信号,传递给小波多尺度重构模块(5);所述的多尺度多传感器信息融合,具体按照以下步骤实施:假设,传感器的测量模型为:X<sub>i</sub>=X+V<sub>i</sub>,i=1,2,…,m,其中,X表示要估计的传感器的真值,V<sub>i</sub>为均值为0,标准差<img file="FSB00000718518300021.GIF" wi="52" he="59" />为有限值的测量噪声,近似认为传感器测量数据服从高斯分布,m个传感器的标准差分别为<img file="FSB00000718518300022.GIF" wi="298" he="60" />各传感器的加权因子分别为归一化权值W<sub>1</sub>,W<sub>2</sub>,…,W<sub>i</sub>,各个传感器的标准差满足:<maths num="0003"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>W</mi><mi>i</mi></msub><mo>=</mo><mn>1,0</mn><mo>&le;</mo><msub><mi>W</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>融合后的X值和加权因子满足:<maths num="0004"><![CDATA[<math><mrow><mover><mi>X</mi><mo>^</mo></mover><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>W</mi><mi>i</mi></msub><msub><mi>X</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths>则由X<sub>1</sub>,X<sub>2</sub>,…,X<sub>i</sub>彼此独立,时域融合后序列的标准差:<maths num="0005"><![CDATA[<math><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>=</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>(</mo><mi>X</mi><mo>-</mo><mover><mi>X</mi><mo>^</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>=</mo><mi>E</mi><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><msub><mi>W</mi><mi>i</mi></msub><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mover><mi>X</mi><mo>^</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><msub><mi>W</mi><mi>i</mi></msub><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>i</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths>根据多元函数求值理论,求出总均方误差最小时所对应的加权因子为:<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>W</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><msubsup><mi>&sigma;</mi><mi>i</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>,</mo></mrow></math>]]></maths>时域加权融合的最小均方误差为:<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>min</mi><mn>2</mn></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><msubsup><mi>&sigma;</mi><mi>i</mi><mn>2</mn></msubsup></mfrac><mo>,</mo></mrow></math>]]></maths>根据单一尺度上的最优加权信息融合公式,小波域随机信号的重构公式为:<img file="FSB00000718518300034.GIF" wi="1258" he="116" />X<sub>J</sub>是经过J层分解后重构回去的融合数据,<img file="FSB00000718518300035.GIF" wi="72" he="43" />是第J层的尺度函数,ψ<sub>j,k</sub>是从1到J层的小波函数,h<sub>n</sub>(J,k)是第J层上k个近似系数,g<sub>n</sub>(j,k)是1到J层的k个细节系数,经过J层分解后重构回去的融合结果:<img file="FSB00000718518300036.GIF" wi="1171" he="125" />式中,各传感器相应近似信号和细节信号的加权因子分别为Wh<sub>i</sub>和Wg<sub>i</sub>,同理有权值的归一化条件,h<sub>n</sub>,g<sub>n</sub>是独立不相关的,多尺度数据融合重构序列的方差:<maths num="0008"><![CDATA[<math><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>=</mo><msubsup><mi>&sigma;</mi><mi>h</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>g</mi><mn>2</mn></msubsup></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><mo>=</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>V</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Wh</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>J</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>&phi;</mi><mrow><mi>J</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>+</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>D</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Wg</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>g</mi><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>&psi;</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mtext>],</mtext></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><mo>=</mo><mi>E</mi><mo>[</mo><msubsup><mi>X</mi><mi>V</mi><mn>2</mn></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>Wh</mi><mi>i</mi><mn>2</mn></msubsup><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&delta;</mi><mi>i</mi></msub><msub><mi>&phi;</mi><mrow><mi>J</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>]</mo><mo>+</mo><mi>E</mi><mo>[</mo><msubsup><mi>X</mi><mi>D</mi><mn>2</mn></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>Wg</mi><mi>i</mi><mn>2</mn></msubsup><munder><mi>&Sigma;</mi><mi>j</mi></munder><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>j</mi></msup><msub><mi>&delta;</mi><mi>i</mi></msub><msub><mi>&psi;</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>]</mo></mrow></math>]]></maths>关于各加权因子的多元二次函数,结合加权因子的约束条件为:<maths num="0011"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Wh</mi><mi>i</mi></msub><mo>=</mo><mn>1,0</mn><mo>&le;</mo><msub><mi>Wh</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn></mrow></math>]]></maths>和<maths num="0012"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Wg</mi><mi>i</mi></msub><mo>=</mo><mn>1,0</mn><mo>&le;</mo><msub><mi>Wg</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>根据拉格朗日方程组,求取多元函数的极值可得总均方误差最小时所对应的加权因子为:<maths num="0013"><![CDATA[<math><mrow><msubsup><mi>Wh</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><msubsup><mi>&sigma;</mi><mi>hi</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths>和<maths num="0014"><![CDATA[<math><mrow><msubsup><mi>Wg</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><msubsup><mi>&sigma;</mi><mi>gi</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>,</mo></mrow></math>]]></maths>此时所对应的最小均方误差为:<maths num="0015"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>min</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>&sigma;</mi><mrow><mi>h</mi><mi>min</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mrow><mi>g</mi><mi>min</mi></mrow><mn>2</mn></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>J</mi></msup><msubsup><mi>&sigma;</mi><mi>hi</mi><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mn>1</mn><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>j</mi></msup><msubsup><mi>&sigma;</mi><mi>gi</mi><mn>2</mn></msubsup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤3:小波多尺度重构模块(5)对步骤2得到的融合后的精确多尺度上的微机电陀螺仪信号进行小波多尺度重构,得到原始尺度上的微机电陀螺信号,传递给组合微机电惯性测量设备模块(6);所述的小波多尺度重构,具体按照以下步骤实施:<maths num="0016"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>l</mi></munder><mi>h</mi><mrow><mo>(</mo><mn>2</mn><mi>l</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mi>l</mi></munder><mi>g</mi><mrow><mo>(</mo><mn>2</mn><mi>l</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mi>d</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>式中,x(i,l)和d(i,l)分别是尺度i上的近似信号和细节信号,h(2l-k)和g(2l-k)分别是尺度i上的尺度系数和小波系数,x(i+1,k)是经过小波重构得到的尺度i+1上的近似信号,即尺度i+1上陀螺仪的近似信息;步骤4:组合微机电惯性测量设备模块(6)将步骤3得到的原始尺度上的微机电陀螺信号整理成需要的形式输出,完成基于MEMS的小波域多传感器信息融合。
地址 710048 陕西省西安市金花南路5号