发明名称 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置
摘要 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置,属于惯性导航技术领域。该装置由惯性测量单元、旋转平台、控制单元和计算单元构成;惯性测量单元通过集成的陀螺仪和加速度计分别测量角速率和加速度;计算单元通过角速率和加速度分别计算惯性测量单元的倾角数值;旋转平台控制惯性测量单元旋转以改变姿态;通过旋转前后不同姿态下的倾角数值和陀螺漂移的数学关系实现陀螺漂移的估计,并补偿漂移以提高姿态定位精度。本发明的特点在于:无需磁场计或GPS等外部传感器的辅助,仅通过控制惯性测量单元旋转即可估计陀螺仪的漂移,减少长时间内姿态角尤其是航向角的积累误差。
申请公布号 CN103712622A 申请公布日期 2014.04.09
申请号 CN201310750475.8 申请日期 2013.12.31
申请人 清华大学 发明人 吴哲明;孙振国;张文增;陈强
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种基于惯性测量单元旋转的陀螺漂移估计补偿方法,其特征在于该方法包括以下步骤:1)通过惯性测量单元所集成的陀螺仪和加速度计测量角速率和加速度;通过陀螺仪测量的角速率计算时间间隔Δt内惯性测量单元的倾角的变化数值 <mrow> <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> </msub> <mo>:</mo> </mrow> <mrow> <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> </msub> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>sin</mi> <mi></mi> <mi>&phi;</mi> <mi>tan</mi> <mi>&theta;</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi></mi> <mi>&phi;</mi> <mi>tan</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <mi>&phi;</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi></mi> <mi>&phi;</mi> <mi>sec</mi> <mi>&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>z</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>&Delta;t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>式中 <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&omega;</mi> <mi>z</mi> </msub> </mtd> </mtr> </mtable> </mfenced>是陀螺仪测量得到的角速率,φ为横滚角,θ为俯仰角, <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> </msub>为通过角速率计算的倾角变化值;2)通过加速度计测量的加速度利用公式(2)计算惯性测量单元的倾角,并获得时间间隔Δt内倾角变化 <mrow> <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>a</mi> </msub> <mo>:</mo> </mrow> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>a</mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>A</mi> <mi>mx</mi> </msub> <mo>,</mo> <mo>-</mo> <msub> <mi>A</mi> <mi>mz</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>arcsin</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>A</mi> <mi>my</mi> </msub> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>A</mi> <mi>m</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>3)式中是Amx、Amy、Amz分别为三轴加速度计测量的加速度,Am=[Amx  Amy  Amz], <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>a</mi> </msub>为通过加速度计算的倾角变化值;4)计算2)获得的倾角变化 <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> </msub>和3)获得的倾角变化 <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>a</mi> </msub>的差值r: <mrow> <mi>r</mi> <mo>=</mo> <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> </msub> <mo>-</mo> <msub> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>&Delta;&phi;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mi>a</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>5)确定横滚角为φ、俯仰角为θ时倾角变化的差值r和陀螺漂移b的数学关系式:r=Vb   (4)式中, <mrow> <mi>V</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>sin</mi> <mi></mi> <mi>&phi;</mi> <mi>tan</mi> <mi>&theta;</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi></mi> <mi>&phi;</mi> <mi>tan</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mi>&phi;</mi> <mi>sec</mi> <mi>&theta;</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>cos</mi> <mi></mi> <mi>&phi;</mi> <mi>sec</mi> <mi>&theta;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>6)控制惯性测量单元旋转,旋转轴方向的单位向量为u,旋转角度为Δα,旋转后的新姿态下的横滚角为φ′、俯仰角为θ′;7)重复步骤1)~5),确定横滚角为φ′、俯仰角为θ′时倾角变化的差值r’和陀螺漂移b的数学关系式:r′=V′b   (5)式中, <mrow> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>sin</mi> <mi></mi> <msup> <mi>&phi;</mi> <mo>&prime;</mo> </msup> <mi>tan</mi> <msup> <mi>&theta;</mi> <mo>&prime;</mo> </msup> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi></mi> <msup> <mi>&phi;</mi> <mo>&prime;</mo> </msup> <mi>tan</mi> <msup> <mi>&theta;</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <msup> <mi>&phi;</mi> <mo>&prime;</mo> </msup> <mi>sec</mi> <msup> <mi>&theta;</mi> <mo>&prime;</mo> </msup> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>cos</mi> <mi></mi> <msup> <mi>&phi;</mi> <mo>&prime;</mo> </msup> <mi>sec</mi> <msup> <mi>&theta;</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>8)联立旋转前后不同姿态下的数学关系式,通过最小二乘法计算陀螺仪的漂移并补偿陀螺仪角速率,输出校正后的姿态测量结果。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室