发明名称 采用奇异摄动方法进行圆锥误差与划船误差补偿的捷联惯导系统
摘要 本发明公开了一种采用奇异摄动方法进行圆锥误差与划船误差补偿的捷联惯导系统,该捷联惯导系统是在惯性测量单元与姿态计算单元之间嵌入圆锥误差补偿奇异摄动模型;在惯性测量单元与导航计算单元之间嵌入划船误差补偿奇异摄动模型。本发明的补偿是在一个采样周期T内载体运动的变化是一种摄动过程,利用惯性元件在起、始采样点的实际输出值建立强迫奇异摄动模型。根据强迫奇异摄动模型对圆锥误差、划船误差进行补偿,获得采样周期内姿态更新旋转矢量以及速度更新量,通过调节可变参数逼近准确的姿态更新量和速度更新量,从而达到提高捷联惯性导航计算单元算法精度的要求。
申请公布号 CN101294811B 申请公布日期 2010.06.09
申请号 CN200810113368.3 申请日期 2008.05.29
申请人 北京航空航天大学 发明人 富立;王玲玲;刘文丽
分类号 G01C21/16(2006.01)I;G01C21/18(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 周长琪
主权项 1.一种采用奇异摄动方法进行圆锥误差与划船误差补偿的捷联惯导系统,其特征在于:该捷联惯导系统是在惯性测量单元与姿态计算单元之间嵌入圆锥误差补偿奇异摄动模型;在惯性测量单元与导航计算单元之间嵌入划船误差补偿奇异摄动模型;惯性组件用于输出载体的角速度<img file="F2008101133683C00011.GIF" wi="73" he="49" />载体的加速度<img file="F2008101133683C00012.GIF" wi="61" he="62" />所述角速度<img file="F2008101133683C00013.GIF" wi="71" he="51" />加速度<img file="F2008101133683C00014.GIF" wi="35" he="63" />经信号处理单元进行滤波处理后输出滤波后角速度<img file="F2008101133683C00015.GIF" wi="68" he="50" />滤波后加速度<img file="F2008101133683C00016.GIF" wi="59" he="62" />滤波后角速度<img file="F2008101133683C00017.GIF" wi="41" he="50" />经圆锥误差补偿奇异摄动模型补偿后输出补偿后的角速度<img file="F2008101133683C00018.GIF" wi="33" he="37" />给姿态计算单元;姿态计算单元对接收的补偿后的角速度<img file="F2008101133683C00019.GIF" wi="33" he="38" />进行解析获得该角速度<img file="F2008101133683C000110.GIF" wi="34" he="38" />下的位姿信息并输出给导航计算单元;滤波后加速度<img file="F2008101133683C000111.GIF" wi="28" he="63" />经划船误差补偿奇异摄动模型补偿后输出补偿后的加速度<img file="F2008101133683C000112.GIF" wi="24" he="51" />给导航计算单元;导航计算单元对接收的位姿信息、补偿后的加速度<img file="F2008101133683C000113.GIF" wi="24" he="51" />进行导航解析,从而获得导航信息输出;所述的圆锥误差补偿奇异摄动模型A为:<maths num="0001"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&epsiv;</mi><mn>1</mn></msub><mover><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mover><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo></mover><mo>+</mo><msub><mover><mi>&beta;</mi><mo>&RightArrow;</mo></mover><mn>1</mn></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mi>n</mi></msub><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><mo>,</mo><msub><mover><mi>w</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>w</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中ε<sub>1</sub>表示第一奇异摄动小参数,其取值为[0,1];<img file="F2008101133683C000115.GIF" wi="33" he="50" />表示在6~2000Hz高动态环境下快变状态的角加加速度矢量;<img file="F2008101133683C000116.GIF" wi="33" he="49" />表示在6~2000Hz高动态环境下慢变状态的角加速度矢量;<img file="F2008101133683C000117.GIF" wi="39" he="63" />表示在一个采样周期T内载体的角振动频率f<sub>w</sub>有关参数矢量;即<img file="F2008101133683C000118.GIF" wi="350" he="101" />N表示与角振动频率f<sub>w</sub>相关的一个常量,该常量为浮点数;T表示一个采样周期;<img file="F2008101133683C000119.GIF" wi="47" he="50" />表示陀螺在当前采样点n的实际输出的角速率w<sub>1</sub>;<img file="F2008101133683C000120.GIF" wi="75" he="52" />表示陀螺在下一时刻采样点n+1的实际输出的角速率w<sub>2</sub>;所述的划船误差补偿奇异摄动模型B为:<maths num="0002"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&epsiv;</mi><mn>2</mn></msub><mover><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mover><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo></mover><mo>+</mo><msub><mover><mi>&beta;</mi><mo>&RightArrow;</mo></mover><mn>2</mn></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mi>n</mi></msub><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub><mo>,</mo><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>f</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中ε<sub>2</sub>表示第二奇异摄动小参数,其取值为[0,1];<img file="F2008101133683C000122.GIF" wi="24" he="62" />表示在6~2000Hz高动态环境下快变状态的载体加加速度矢量;<img file="F2008101133683C000123.GIF" wi="24" he="61" />表示在6~2000Hz高动态环境下慢变状态的载体加速度矢量;<img file="F2008101133683C000124.GIF" wi="47" he="67" />表示在一个采样周期T内载体的线振动频率f<sub>a</sub>有关参数矢量;即<img file="F2008101133683C00021.GIF" wi="337" he="104" />M表示与线振动频率f<sub>a</sub>相关的一个常量,该常量为浮点数;T表示一个采样周期;<img file="F2008101133683C00022.GIF" wi="35" he="63" />表示加速度计在当前采样点n的实际输出的加速度f<sub>1</sub>;<img file="F2008101133683C00023.GIF" wi="62" he="63" />表示加速度计在下一时刻采样点n+1的实际输出的加速度f<sub>2</sub>。
地址 100083 北京市海淀区学院路37号