发明名称 一种惯性器件的温度误差补偿方法
摘要 一种惯性器件的温度误差补偿新方法,是指对惯性器件输出电压及供电电压同时进行AD采样,利用两者采样值的比值进行陀螺的角速率信息(或加速度计的加速度信息)解算;再建立角速率信息(或加速度信息)的温度误差模型并进行补偿。该方法可消除惯性器件供电电压及AD转换电路参考电压的温度漂移对器件输出精度的影响,可大大简化惯性器件的温度误差模型,提高器件的温度误差补偿精度。
申请公布号 CN101915588B 申请公布日期 2011.11.09
申请号 CN201010231282.8 申请日期 2010.07.14
申请人 北京航空航天大学 发明人 徐烨烽;张仲毅;李魁;杨国梁;徐韬
分类号 G01C25/00(2006.01)I;G01P21/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉
主权项 1.一种惯性器件的温度误差补偿方法,特征在于其实现步骤如下:第一步,对惯性器件的输出电压及供电电压V<sub>cc</sub>同时进行AD采样,其中,陀螺的输出电压<img file="FSB00000590827300011.GIF" wi="80" he="62" />及加速度计的输出电压<img file="FSB00000590827300012.GIF" wi="80" he="63" />经AD采样后的数字量D<sub>g</sub>及D<sub>a</sub>可表示为:<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>D</mi><mi>g</mi></msub><mo>=</mo><mfrac><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>V</mi><mi>cc</mi></msub><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>V</mi><mi>cc</mi></msub><mi>&omega;</mi></mrow><msub><mi>V</mi><mi>ref</mi></msub></mfrac><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>k</mi><mn>2</mn></msub><mfrac><msub><mi>V</mi><mi>cc</mi></msub><msub><mi>V</mi><mi>ref</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>D</mi><mi>a</mi></msub><mo>=</mo><mfrac><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>V</mi><mi>cc</mi></msub><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mi>V</mi><mi>cc</mi></msub><mi>a</mi></mrow><msub><mi>V</mi><mi>ref</mi></msub></mfrac><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><mi>a</mi><mo>)</mo></mrow><msub><mi>k</mi><mn>2</mn></msub><mfrac><msub><mi>V</mi><mi>cc</mi></msub><msub><mi>V</mi><mi>ref</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>ω表示角速度,a表示加速度,k<sub>1</sub>表示陀螺刻度系数,k<sub>3</sub>表示加速度计的刻度系数;对惯性器件供电电压V<sub>cc</sub>进行AD采样,得到的采样值<img file="FSB00000590827300014.GIF" wi="85" he="66" />为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>=</mo><mfrac><msub><mi>V</mi><mi>cc</mi></msub><msub><mi>V</mi><mi>ref</mi></msub></mfrac><mo>&times;</mo><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>V<sub>ref</sub>为AD采样电路的参考电压,k<sub>2</sub>为AD采样电路的比例系数;第二步,将陀螺输出信号的采样值D<sub>g</sub>、加速度计输出信号的采样值D<sub>a</sub>与供电电压的采样值<img file="FSB00000590827300016.GIF" wi="83" he="67" />相比,求得陀螺的角速度ω及加速度计的加速度a的计算公式为:<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>D</mi><mi>g</mi></msub><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>&omega;</mi></mtd></mtr><mtr><mtd><mfrac><msub><mi>D</mi><mi>a</mi></msub><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><mi>a</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>当外界温度变化时,惯性器件的供电电压由V<sub>cc</sub>变为V<sub>cc</sub>′,AD采样电路的参考电压由V<sub>ref</sub>变为V<sub>ref</sub>′,此时惯性器件的输出电压及供电电压的AD采样值为:<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><msub><mi>D</mi><mi>g</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>k</mi><mn>2</mn></msub><mfrac><msup><msub><mi>V</mi><mi>cc</mi></msub><mo>&prime;</mo></msup><msup><msub><mi>V</mi><mi>ref</mi></msub><mo>&prime;</mo></msup></mfrac></mtd></mtr><mtr><mtd><msup><msub><mi>D</mi><mi>a</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><mi>a</mi><mo>)</mo></mrow><msub><mi>k</mi><mn>2</mn></msub><mfrac><msup><msub><mi>V</mi><mi>cc</mi></msub><mo>&prime;</mo></msup><msup><msub><mi>V</mi><mi>ref</mi></msub><mo>&prime;</mo></msup></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msup><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>&prime;</mo></msup><mo>=</mo><mfrac><msup><msub><mi>V</mi><mi>cc</mi></msub><mo>&prime;</mo></msup><msup><msub><mi>V</mi><mi>ref</mi></msub><mo>&prime;</mo></msup></mfrac><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>将(4)式和(5)式相比得:<maths num="0006"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msup><msub><mi>D</mi><mi>g</mi></msub><mo>&prime;</mo></msup><msup><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>&prime;</mo></msup></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>&omega;</mi></mtd></mtr><mtr><mtd><mfrac><msup><msub><mi>D</mi><mi>a</mi></msub><mo>&prime;</mo></msup><msup><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>&prime;</mo></msup></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><mi>a</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>对比式(3)及式(6)得,<maths num="0007"><![CDATA[<math><mrow><msub><mi>D</mi><mi>g</mi></msub><mo>/</mo><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>=</mo><msup><msub><mi>D</mi><mi>g</mi></msub><mo>&prime;</mo></msup><mo>/</mo><msup><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>&prime;</mo></msup><mo>,</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>D</mi><mi>a</mi></msub><mo>/</mo><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>=</mo><msup><msub><mi>D</mi><mi>a</mi></msub><mo>&prime;</mo></msup><mo>/</mo><msup><msub><mi>D</mi><msub><mi>V</mi><mi>cc</mi></msub></msub><mo>&prime;</mo></msup><mo>,</mo></mrow></math>]]></maths>利用D<sub>g</sub>与<img file="FSB00000590827300026.GIF" wi="95" he="65" />的比值计算得到的ω值不随V<sub>cc</sub>及V<sub>ref</sub>变化;利用D<sub>a</sub>与<img file="FSB00000590827300027.GIF" wi="81" he="65" />的比值计算得到的加速度a值不随V<sub>cc</sub>及V<sub>ref</sub>变化;第三步,根据第二步计算得到的角速度ω中包含陀螺自身的温度误差Δω(T),建立温度误差Δω(T)与温度T之间的模型如式(7)所示,并用该如式(7)的模型对陀螺的输出信号进行温度误差补偿:Δω(T)=fT<sup>2</sup>+bT+cΔT+d             (7)上式(7)中,T为温度,ΔT为温度梯度,f,b,c,d分别为陀螺温度误差模型参数;该陀螺温度误差模型参数可以通过最小二乘拟合的方法得到,具体计算公式为:<maths num="0009"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>f</mi></mtd></mtr><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>c</mi></mtd></mtr><mtr><mtd><mi>d</mi></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>T</mi><mn>1</mn><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mn>2</mn><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mi>n</mi><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mi>n</mi></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mi>n</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo></msup><mo>&times;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>&omega;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>&omega;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>T</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>&Delta;&omega;</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式(8)中,T<sub>i</sub>(i=1,2…n)表示第i时刻温度的采样值,ΔT<sub>i</sub>(i=1,2…n)表示第i时刻温度梯度的采样值,Δω<sub>i</sub>(T)(i=1,2…n)表示第i时刻陀螺温度误差计算值,可由式(6)计算得到,n表示采样点的个数;根据第二步计算得到的加速度a中包含加速度计自身的温度误差Δa(T),建立温度误差Δa(T)与温度T之间的模型如式(9)所示,并用该如式(9)的模型对加速度计的输出信号进行温度误差补偿:Δa(T)=a′T<sup>2</sup>+b′T+c′ΔT+d′     (9)上式(9)中,T为温度,ΔT为温度梯度,a′,b′,c′,d′分别为加速度计温度误差模型参数;该加速度计温度误差模型参数可以通过最小二乘拟合的方法得到,具体计算公式为:<maths num="0010"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>a</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>b</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>c</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>d</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>T</mi><mn>1</mn><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mn>2</mn><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mi>n</mi><mn>2</mn></msubsup></mtd><mtd><msub><mi>T</mi><mi>n</mi></msub></mtd><mtd><mi>&Delta;</mi><msub><mi>T</mi><mi>n</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo></msup><mo>&times;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>a</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>T</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><mi>&Delta;</mi><msub><mi>a</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式(10)中,T<sub>i</sub>(i=1,2…n)表示第i时刻温度的采样值,ΔT<sub>i</sub>(i=1,2…n)表示第i时刻温度梯度的采样值,Δa<sub>i</sub>(T)(i=1,2…n)表示第i时刻加速度计温度误差计算值,可由式(6)计算得到,n表示采样点的个数;第四步,利用第三步所建的温度误差模型,对完成电路漂移补偿后的惯性器件的输出信号进行温度误差补偿。
地址 100191 北京市海淀区学院路37号