发明名称 一种光纤陀螺惯性测量单元的现场标定方法
摘要 一种光纤陀螺惯性测量单元的现场标定方法,本发明涉及一种现场精确标定光纤陀螺惯性测量单元(Fiber Optic Gyro Inertial Measurement Unit,FIMU)误差系数的方法。该方法可在没有精密转台的使用现场,通过FIMU在六个位置上的12次旋转,精确标定出陀螺标度因数、陀螺常值漂移、陀螺安装误差和加速度计常值偏置共15个误差系数。本发明具有精度高且操作简单的特点,不仅提高了FIMU的使用精度,同时也大大提高了标定的效率。
申请公布号 CN100559188C 申请公布日期 2009.11.11
申请号 CN200710064784.4 申请日期 2007.03.26
申请人 北京航空航天大学 发明人 房建成;刘百奇;杨胜;张海鹏;宫晓琳;韩晓英
分类号 G01P21/00(2006.01)I;G01P15/14(2006.01)I 主分类号 G01P21/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 代理人 关 玲;李新华
主权项 1、一种光纤陀螺惯性测量单元的现场标定方法,其特征在于包括以下步骤:(1)将光纤陀螺惯性测量单元放置在一个水平面上,使光纤陀螺惯性测量单元的Z轴向上,启动光纤陀螺惯性测量单元,等待光纤陀螺惯性测量单元预热完毕;(2)将光纤陀螺惯性测量单元在该水平面上逆时针旋转360°,然后保持光纤陀螺惯性测量单元静止,记录旋转及静止过程中光纤陀螺惯性测量单元输出的角速度和加速度;(3)将光纤陀螺惯性测量单元在该水平面上顺时针旋转360°,然后保持光纤陀螺惯性测量单元静止,记录旋转及静止过程中光纤陀螺惯性测量单元输出的角速度和加速度;(4)使光纤陀螺惯性测量单元的Z轴向下放置,重复试验步骤(2)和(3);(5)分别使光纤陀螺惯性测量单元的X轴向上、向下放置,针对每个位置重复试验步骤(2)和(3);(6)分别使光纤陀螺惯性测量单元的Y轴向上、向下放置,针对每个位置重复试验步骤(2)和(3);(7)利用六个位置上的12组数据,根据光纤陀螺惯性测量单元的误差模型,采用光纤陀螺惯性测量单元误差系数计算公式标定出光纤陀螺惯性测量单元的陀螺标度因数、陀螺常值漂移、陀螺安装误差和加速度计常值偏置共四类15个误差系数,所述的光纤陀螺惯性测量单元的误差模型包括角速度通道误差模型和加速度通道误差模型,如下:SU=Mω+DkN=If+B其中,S为光纤陀螺的标度因数,U为光纤陀螺仪的输出,M为光纤陀螺的安装误差,ω为光纤陀螺的输入角速度,D为光纤陀螺的常值漂移,k为加速度计的标度因数,N为加速度计的输出,I为加速度计的安装误差,f为加速度计的输入加速度,B为加速度计的常值偏置;U=[U<sub>x</sub> U<sub>y</sub> U<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,ω=[ω<sub>x</sub> ω<sub>y</sub> ω<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,D=[D<sub>x</sub> D<sub>y</sub> D<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,<maths num="0001"><![CDATA[<math><mrow><mi>S</mi><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>S</mi><mi>x</mi></msub></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd><msub><mi>S</mi><mi>y</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd><msub><mi>S</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>M</mi><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>M</mi><mi>xy</mi></msub></mtd><mtd><msub><mi>M</mi><mi>xz</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>yx</mi></msub></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>M</mi><mi>yz</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>zx</mi></msub></mtd><mtd><msub><mi>M</mi><mi>zy</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>,</mo></mrow></math>]]></maths>N=[N<sub>x</sub> N<sub>y</sub> N<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,f=[f<sub>x</sub> f<sub>y</sub> f<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,B=[B<sub>x</sub> B<sub>y</sub> B<sub>z</sub>]<sup>T</sup><sub>1×3</sub>,<maths num="0003"><![CDATA[<math><mrow><mi>k</mi><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>k</mi><mi>x</mi></msub></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd><msub><mi>k</mi><mi>y</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd><msub><mi>k</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>I</mi><mi>xy</mi></msub></mtd><mtd><msub><mi>I</mi><mi>xz</mi></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>yx</mi></msub></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>I</mi><mi>yz</mi></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>zx</mi></msub></mtd><mtd><msub><mi>I</mi><mi>zy</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>;</mo></mrow></math>]]></maths>所述的光纤陀螺惯性测量单元误差系数计算公式为:<img file="C2007100647840003C5.GIF" wi="1140" he="140" /><maths num="0005"><![CDATA[<math><mrow><msub><mi>D</mi><mi>i</mi></msub><mo>=</mo><msub><mi>S</mi><mi>i</mi></msub><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mrow><mo>(</mo><msup><msub><mi>U</mi><mi>ij</mi></msub><mo>+</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mi>ij</mi></msub><mo>-</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>+</mo></msup><mo>+</mo><msup><msub><mi>U</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo></msup><mo>)</mo></mrow><mi>d&tau;</mi><mo>/</mo><mn>4</mn><mi>t</mi></mrow></math>]]></maths><img file="C2007100647840003C7.GIF" wi="1264" he="115" /><maths num="0006"><![CDATA[<math><mrow><msub><mi>B</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>k</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>N</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>N</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,n,i=x,y,z,n≠i;j表示第j个位置,且i=z时,j=1;i=x时,j=3;i=y时,j=5,上标+和上标-分别表示逆时针旋转和顺时针旋转,t为数据采集时间。
地址 100083北京市海淀区学院路37号