发明名称 一种MEMS惯性测量单元冷启动时温度误差补偿方法
摘要 本发明属于惯性器件的误差补偿领域,具体涉及一种MEMS惯性测量单元冷启动时温度误差补偿方法。本发明包括(1)选取温度误差补偿模型;温补数据采集;温补数据预处理;温补误差模型参数寻优;温度补偿。本发明针对MEMS惯性器件中的陀螺和加速度计提出了一种改进的温度误差补偿模型和相应的、简单易行的工作流程。经过温度补偿的MEMS惯性测量单元,在冷启动时,将受自身温升和外界温度变化的影响变小,这样极大地缩短了MEMS惯性测量单元的准备时间,提高了将其作为导航系统的快速性。
申请公布号 CN103940427A 申请公布日期 2014.07.23
申请号 CN201410085446.9 申请日期 2014.03.11
申请人 哈尔滨工程大学 发明人 何昆鹏;韩继韬;曾建辉;周雪梅;于文浩;王兴彬;单飞;王晨阳;胡守雷;王志强
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 一种MEMS惯性测量单元冷启动时温度误差补偿方法,其特征在于:(1)选取温度误差补偿模型:MEMS加速度计的静态数学模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>A</mi><mi>ind</mi></msub><mo>=</mo><mfrac><mi>E</mi><mi>K</mi></mfrac><mo>=</mo><mi>B</mi><mo>+</mo><msub><mi>a</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000475138300000011.GIF" wi="356" he="132" /></maths>式中:A<sub>ind</sub>—仪表指示的加速度值;E—加速度计的输出;K—标度因数;B—偏值;a<sub>i</sub>—输入加速度;加速度计模型系数偏值B、标度因数K与温度T的数学模型:B(T)=S<sub>B0</sub>+S<sub>B1</sub>TK(T)=S<sub>K0</sub>+S<sub>K1</sub>T式中:S<sub>B0</sub>—加速度计偏值的零阶温度模型系数;S<sub>B1</sub>—加速度计偏值的一阶温度模型系数;S<sub>K0</sub>—加速度计标度因数的零阶温度模型系数;S<sub>K1</sub>—加速度计标度因数的一阶温度模型系数;得到不同温度时的加速度计输出E(T):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>K</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>K</mi><mn>1</mn></mrow></msub><mi>T</mi></mrow></mfrac><mo>=</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mi>T</mi><mo>+</mo><msub><mi>a</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000475138300000012.GIF" wi="576" he="139" /></maths>在T为室温T0时:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mi>B</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mi>i</mi></msub><mo>&DoubleLeftRightArrow;</mo><mfrac><mi>E</mi><mi>K</mi></mfrac><mo>=</mo><mi>B</mi><mo>+</mo><msub><mi>a</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000475138300000013.GIF" wi="691" he="139" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>K</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>K</mi><mn>1</mn></mrow></msub><mi>T</mi></mrow></mfrac><mo>-</mo><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>E</mi><mi>K</mi></mfrac><mo>-</mo><mi>B</mi></mrow>]]></math><img file="FDA0000475138300000014.GIF" wi="786" he="140" /></maths>得温度补偿公式:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>E</mi><mo>=</mo><mo>[</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>K</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>K</mi><mn>1</mn></mrow></msub><mi>T</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>B</mi><mo>-</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mi>T</mi><mo>)</mo></mrow><mo>&times;</mo><mi>K</mi><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mfrac><msub><mi>S</mi><mrow><mi>K</mi><mn>0</mn></mrow></msub><mi>K</mi></mfrac><mo>+</mo><mfrac><msub><mi>S</mi><mrow><mi>K</mi><mn>1</mn></mrow></msub><mi>K</mi></mfrac><mi>T</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>BK</mi><mo>-</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mi>K</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mi>K</mi><mo>&times;</mo><mi>T</mi></mtd></mtr><mtr><mtd><mo>&DoubleLeftRightArrow;</mo><mi>E</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>a</mi><mo>+</mo><mi>bT</mi></mrow></mfrac><mo>&times;</mo><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>+</mo><mi>c</mi><mo>+</mo><mi>dT</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000475138300000015.GIF" wi="1602" he="342" /></maths>针对惯性器件温度补偿模型:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mi>a</mi><mo>+</mo><mi>b</mi><mo>&CenterDot;</mo><mi>T</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>c</mi><mo>+</mo><mi>d</mi><mo>&CenterDot;</mo><mi>T</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000475138300000016.GIF" wi="514" he="138" /></maths>式中:E表示室温下惯性器件的输出值;T表示惯性器件当前温度;E(T)表示T温度下惯性器件输出值;a+b·T表示惯性器件标度因数补偿,a表示其零次项系数,b表示其一次项系数;(c+d·T)表示惯性器件零偏补偿,c表示其零次项系数,d表示其一次项系数;ΔT为前n分钟,n=1,2,3……,惯性器件温度的变化值,建立改进后的惯性器件温度补偿模型如下:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mi>a</mi><mo>+</mo><mi>b</mi><mo>&CenterDot;</mo><mi>T</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>c</mi><mo>+</mo><mi>d</mi><mo>&CenterDot;</mo><mi>T</mi><mo>)</mo></mrow><mo>+</mo><mi>e</mi><mo>&CenterDot;</mo><mi>&Delta;T</mi></mrow>]]></math><img file="FDA0000475138300000021.GIF" wi="655" he="140" /></maths>式中:e为温度变化率系数;(2)温补数据采集:MEMS惯性测量单元的x轴指天,开启系统电源,以1Hz的频率,采集惯性器件和温度传感器输出数据;y轴指天,采集惯性器件和温度传感器输出数据;z轴指天,采集惯性器件和温度传感器输出数据;(3)温补数据预处理:对上一步采集得到惯性器件和温度传感器输出数据进行预处理,得到当前惯性器件温度T,T温度下惯性器件输出值E(T)、及前n分钟,n=1,2,3……,惯性器件温度的变化值ΔT;(4)温补误差模型参数寻优:获得加速度计不同温度下的输出数据E(T),以及对应的温度数据T和温度变化数据ΔT;(5)温度补偿:获得加速度计温度误差模型参数a,b,c,d,e,当惯性测量单元正常工作后,可得到加速度计输出数据E(T)和对应的温度数据T,温度变化ΔT,最终得到加速度计温度补偿后的输出值。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室