发明名称 挠性陀螺仪最优八位置标定方法
摘要 本发明公开了一种挠性陀螺仪最优八位置标定方法,是将挠性陀螺仪安装在三轴位置速率转台上,挠性陀螺仪与数据采集设备相连,数据采集设备与计算机相连;本发明在最优八位置的标定中公开了特定的方位,对获得的最优八位置漂移系数与挠性陀螺静态误差补偿模型G<sub>0</sub>进行的测量值补偿有效地提高了挠性陀螺仪的输出。在惯导测试中心挠性陀螺测试过程中分别采用传统八位置方法和最优八位置方法得到的漂移系数,利用两种方法得到的漂移系数分别客观的对空间另外三位置的陀螺输出补偿后的评价结果,由陀螺测量值剩余平方和可见,利用挠性陀螺仪最优八位置试验设计方法求解的漂移系数进行补偿后的结果较传统八位置方法提高了4~8倍。
申请公布号 CN100554884C 申请公布日期 2009.10.28
申请号 CN200810101156.3 申请日期 2008.02.28
申请人 北京航空航天大学 发明人 富立;王玲玲;刘文丽
分类号 G01C25/00(2006.01)I;G01C19/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 北京永创新实专利事务所 代理人 周长琪
主权项 1、一种挠性陀螺仪最优八位置标定方法,是在将挠性陀螺仪安装在三轴位置速率转台上,挠性陀螺仪与数据采集设备相连,数据采集设备与计算机相连;所述计算机内安装有位置测量软件;其特征在于有下列标定执行步骤:第一步:标定最优八位置方位第一方位:挠性陀螺仪的X测量轴指向“天”,挠性陀螺仪的Y测量轴指向“西”,挠性陀螺仪的Z自转轴指向“北”;第二方位:挠性陀螺仪的X测量轴指向“地”,挠性陀螺仪的Y测量轴指向“北”,挠性陀螺仪的Z自转轴指向“东”;第三方位:挠性陀螺仪的X测量轴指向“北”,挠性陀螺仪的Y测量轴指向“天”,挠性陀螺仪的Z自转轴指向“东”;第四方位:挠性陀螺仪的X测量轴指向“北”,挠性陀螺仪的Y测量轴指向“地”,挠性陀螺仪的Z自转轴指向“西”;第五方位:挠性陀螺仪的X测量轴指向“东”,挠性陀螺仪的Y测量轴指向“南”,挠性陀螺仪的Z自转轴指向“地”;第六方位:挠性陀螺仪的X测量轴指向“南”,挠性陀螺仪的Y测量轴指向“东”,挠性陀螺仪的Z自转轴指向“天”;第七方位:挠性陀螺仪的X测量轴指向“西”,挠性陀螺仪的Y测量轴指向“南”,挠性陀螺仪的Z自转轴指向“天”;第八方位:挠性陀螺仪的X测量轴指向“南”,挠性陀螺仪的Y测量轴指向“西”,挠性陀螺仪的Z自转轴指向“地”;第二步:获取漂移系数(A)对传统八位置下的数据进行挠性陀螺静态误差模型G<sub>1</sub>的基于最小二乘法解析获得传统八位置漂移系数;(B)对最优八位置下的数据进行挠性陀螺静态误差模型G<sub>1</sub>的基于最小二乘法解析获得最优八位置漂移系数;所述挠性陀螺静态误差模型<maths num="0001"><![CDATA[<math><mrow><msub><mi>G</mi><mn>1</mn></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>i</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>U</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>U</mi><mn>1</mn></msub></mtd><mtd><msub><mi>U</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mn>1</mn></msub></mtd><mtd><msub><mi>V</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&omega;</mi><mi>Y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>X</mi></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>U</mi><mn>3</mn></msub></mtd><mtd><msub><mi>U</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mn>3</mn></msub></mtd><mtd><msub><mi>V</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mi>X</mi></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mi>Y</mi></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>U</mi><mn>5</mn></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mn>5</mn></msub></mtd></mtr></mtable></mfenced><msub><mi>a</mi><mi>Z</mi></msub><mo>,</mo></mrow></math>]]></maths>其中,<maths num="0002"><![CDATA[<math><mrow><msub><mi>U</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mi>cos</mi><mrow><mo>(</mo><mi>&epsiv;</mi><mo>+</mo><mi>&xi;</mi><mo>)</mo></mrow></mrow><mrow><msub><mrow><mo>(</mo><mi>SF</mi><mo>)</mo></mrow><mi>Y</mi></msub><mi>cos</mi><mi>&xi;</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>U</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mi>sin</mi><mrow><mo>(</mo><mi>&epsiv;</mi><mo>+</mo><mi>&xi;</mi><mo>)</mo></mrow></mrow><mrow><msub><mrow><mo>(</mo><mi>SF</mi><mo>)</mo></mrow><mi>Y</mi></msub><mi>cos</mi><mi>&xi;</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>V</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mo>-</mo><mi>sin</mi><mi>&epsiv;</mi></mrow><mrow><msub><mrow><mo>(</mo><mi>SF</mi><mo>)</mo></mrow><mi>X</mi></msub><mi>cos</mi><mi>&xi;</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>V</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mi>cos</mi><mi>&epsiv;</mi></mrow><mrow><msub><mrow><mo>(</mo><mi>SF</mi><mo>)</mo></mrow><mi>X</mi></msub><mi>cos</mi><mi>&xi;</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>U<sub>0</sub>=U<sub>1</sub>×D(X)<sub>F</sub>+U<sub>2</sub>×D(Y)<sub>F</sub>,V<sub>0</sub>=V<sub>1</sub>×D(X)<sub>F</sub>+V<sub>2</sub>×D(Y)<sub>F</sub>,U<sub>3</sub>=U<sub>1</sub>×D(X)<sub>X</sub>+U<sub>2</sub>×D(Y)<sub>X</sub>,U<sub>4</sub>=U<sub>1</sub>×D(X)<sub>Y</sub>+U<sub>2</sub>×D(Y)<sub>Y</sub>,V<sub>3</sub>=V<sub>1</sub>×D(X)<sub>X</sub>+V<sub>2</sub>×D(Y)<sub>X</sub>,V<sub>4</sub>=V<sub>1</sub>×D(X)<sub>Y</sub>+V<sub>2</sub>×D(Y)<sub>Y</sub>,U<sub>5</sub>=U<sub>1</sub>×D(X)<sub>Z</sub>+U<sub>2</sub>×D(Y)<sub>Z</sub>,V<sub>5</sub>=V<sub>1</sub>×D(X)<sub>Z</sub>+V<sub>2</sub>×D(Y)<sub>Z</sub>,式中:i<sub>x</sub>表示挠性陀螺仪X测量轴的力矩器电流所对应的脉冲数,i<sub>y</sub>表示挠性陀螺仪Y测量轴的力矩器电流所对应的脉冲数,ω<sub>X</sub>表示地球自转角速度在挠性陀螺仪X测量轴上的分量,ω<sub>Y</sub>表示地球自转角速度在挠性陀螺仪Y测量轴上的分量,a<sub>X</sub>表示挠性陀螺仪X测量轴上的加速度分量,a<sub>Y</sub>表示挠性陀螺仪Y测量轴上的加速度分量,a<sub>Z</sub>表示挠性陀螺仪Z自转轴上的加速度分量,(SF)<sub>X</sub>表示挠性陀螺仪X测量轴的力矩器刻度系数,(SF)<sub>Y</sub>表示挠性陀螺仪Y测量轴的力矩器刻度系数,ε表示挠性陀螺仪的力矩器X轴与挠性陀螺仪的壳体X轴之间的夹角,ξ表示挠性陀螺仪的力矩器Y轴与挠性陀螺仪的壳体Y轴之间的夹角;第三步:对最优八位置方位的测量值进行补偿利用挠性陀螺静态误差补偿模型G<sub>0</sub>与最优八位置漂移系数对挠性陀螺仪输出测量值进行补偿获得补偿后的测量值;所述挠性陀螺静态误差补偿模型为<maths num="0006"><![CDATA[<math><mrow><msub><mi>G</mi><mn>0</mn></msub><mo>=</mo><mfenced open='{' close='}'><mtable><mtr><mtd><mi>D</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>F</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>X</mi></msub><msub><mi>a</mi><mi>X</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>Y</mi></msub><msub><mi>a</mi><mi>Y</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>Z</mi></msub><msub><mi>a</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mi>D</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mi>F</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mi>X</mi></msub><msub><mi>a</mi><mi>X</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mi>Y</mi></msub><msub><mi>a</mi><mi>Y</mi></msub><mo>+</mo><mi>D</mi><msub><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mi>Z</mi></msub><msub><mi>a</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>式中,D(X)表示挠性陀螺仪X测量轴的漂移量,D(Y)表示挠性陀螺仪Y测量轴的漂移量,D(X)<sub>F</sub>表示挠性陀螺仪绕X测量轴与加速度无关的漂移系数,D(Y)<sub>F</sub>表示挠性陀螺仪绕Y测量轴与加速度无关的漂移系数,D(X)<sub>X</sub>表示X测量轴中挠性陀螺仪绕X测量轴与加速度一次方有关的漂移系数,D(X)<sub>Y</sub>表示X测量轴中挠性陀螺仪绕Y测量轴与加速度一次方有关的漂移系数,D(X)<sub>Z</sub>表示X测量轴中挠性陀螺仪绕Z自转轴与加速度一次方有关的漂移系数,D(Y)<sub>X</sub>表示Y测量轴中挠性陀螺仪绕X测量轴与加速度一次方有关的漂移系数,D(Y)<sub>Y</sub>表示Y测量轴中挠性陀螺仪绕Y测量轴与加速度一次方有关的漂移系数,D(Y)<sub>Z</sub>表示Y测量轴中挠性陀螺仪绕Z自转轴与加速度一次方有关的漂移系数,a<sub>X</sub>表示挠性陀螺仪X测量轴上的加速度分量,a<sub>Y</sub>表示挠性陀螺仪Y测量轴上的加速度分量,a<sub>Z</sub>表示挠性陀螺仪Z自转轴上的加速度分量。
地址 100083北京市海淀区学院路37号