发明名称 一种动态旋转调制的陀螺仪寻北方法
摘要 本发明属于惯性测量技术或地球物理实验与仪器领域,具体涉及一种动态旋转调制的陀螺仪寻北方法。目的是在载体存在大姿态角的条件下,且无需调平、无需知道测点纬度的情况下实现快速、精确寻北测量,提高测地作业速度。它使用以角速度Ω匀速转动的转台,转台上安装1个陀螺仪,2个加速度计;首先,利用2个加速度计求解:俯仰角和横滚角;然后,利用陀螺仪求解地理纬度L和航向角;在计算得到初值之后,采用实时递推的最小二乘方法算法提高精度。本发明的优点是无需调平、无需知道测点纬度的情况下实现快速、精确寻北测量,提高测地作业速度。
申请公布号 CN102840856B 申请公布日期 2015.08.26
申请号 CN201110173699.8 申请日期 2011.06.24
申请人 西安测绘研究所 发明人 白云超;田育民;刘思伟;陈炜;马小辉;王成宾;李瑞锋;高扬
分类号 G01C19/38(2006.01)I;G01C19/72(2006.01)I 主分类号 G01C19/38(2006.01)I
代理机构 核工业专利中心 11007 代理人 李烨
主权项 一种动态旋转调制的陀螺仪寻北方法,其特征在于:使用以角速度Ω匀速转动的转台,转台上安装1个陀螺仪,2个加速度计;转台以垂直方向为轴旋转;所述陀螺仪水平安装在支架安装孔内,其敏感轴线与转轴正交,且转台过零位时,陀螺测量坐标系x轴方向指向光电圆光栅尺读数头的零位线;2只加速度计3分别沿双轴陀螺的两个敏感轴的方向,垂直安装在单轴转台支架上方互相垂直的两个安装孔内;设东北天坐标系(OXnYnZn)(n系);旋转台体坐标系(OXbYbZb)(b系)的Xb轴指向光栅零位,Yb轴指向光栅90°方向,Zb轴指向旋转轴;陀螺测量坐标系(OXmYmZm)(m系)的Xm轴指向加速度计1,Ym方向指向加速度计2,Zm轴指向旋转轴;当地地理纬度为L;寻北仪台体的俯仰/横滚/航向角分别为0、γ、φ,求解过程包括如下步骤:首先,利用加速度计求解:俯仰角和横滚角;当惯性测量单元IMU相对台体的瞬时角位置为α=Ω·t+φ<sub>0</sub>时,其中Ω为台体转速,φ<sub>0</sub>为初始相位,读取加速度计测量值为f<sub>mx</sub>、f<sub>my</sub>,变换到台体坐标系后,x轴、y轴上的比力为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>f</mi><mi>bx</mi></msub></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>by</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>f</mi><mi>mx</mi></msub><mi>cos</mi><mi>&alpha;</mi><mo>-</mo><msub><mi>f</mi><mi>my</mi></msub><mi>sin</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>mx</mi></msub><mi>sin</mi><mi>&alpha;</mi><mo>+</mo><msub><mi>f</mi><mi>my</mi></msub><mi>cos</mi><mi>&alpha;</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000683933880000011.GIF" wi="592" he="145" /></maths>则θ=sin<sup>‑1</sup>(f<sub>by</sub>/g)γ=sin<sup>‑1</sup>(‑f<sub>bx</sub>/g/cosθ)然后,利用陀螺仪求解地理纬度L和航向角Φ;所述寻北方法的粗寻北方案采用的就是正交解调法,其数学原理是:ω<sub>g</sub>=‑[(‑cosγsinφ+sinγsinθcosφ)cosL‑sinγcosθsinL]Ω<sub>e</sub>sinα+(cosθcosφcosL+sinθsinL)Ω<sub>e</sub>cosα由于:设姿态角θ、γ均为0时,可以简化为:ω<sub>g</sub>=(sinφcosα+cosφsinα)Ω<sub>e</sub>cosL=Ω<sub>e</sub>cosL·sin(α+φ)令:A=Ω<sub>e</sub>cosL·cosφ,B=Ω<sub>e</sub>cosL·sinφ,则有:ω<sub>g</sub>=(Asinα+Bcosα)分别用sinα、cosα(α=Ω·t+φ<sub>0</sub>)乘以上式,然后两端积分可得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>A</mi><mo>^</mo></mover><mo>=</mo><mfrac><mi>&Omega;</mi><mn>4</mn></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>T</mi></msubsup><msub><mi>&omega;</mi><mi>g</mi></msub><mi>sgn</mi><mo>[</mo><mi>cos</mi><mrow><mo>(</mo><mi>&Omega;t</mi><mo>)</mo></mrow><mo>]</mo><mi>dt</mi></mtd></mtr><mtr><mtd><mover><mi>B</mi><mo>^</mo></mover><mo>=</mo><mfrac><mi>&Omega;</mi><mn>4</mn></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>T</mi></msubsup><msub><mi>&omega;</mi><mi>g</mi></msub><mi>sgn</mi><mo>[</mo><mi>cos</mi><mrow><mo>(</mo><mi>&Omega;t</mi><mo>)</mo></mrow><mo>]</mo><mi>dt</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000683933880000012.GIF" wi="562" he="247" /></maths>于是有:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>&phi;</mi><mo>^</mo></mover><mi>d</mi></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mover><mi>B</mi><mo>^</mo></mover><mo>/</mo><mover><mi>A</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>,</mo><mi>L</mi><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msqrt><mrow><mo>(</mo><msup><mover><mi>A</mi><mo>^</mo></mover><mn>2</mn></msup><mo>+</mo><msup><mover><mi>B</mi><mo>^</mo></mover><mn>2</mn></msup><mo>)</mo></mrow><mo>/</mo><mn>2</mn><msub><mi>&Omega;</mi><mi>e</mi></msub></msqrt><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000683933880000013.GIF" wi="978" he="83" /></maths>即转台旋转一周后即可给出真北方位的概略值。
地址 710054 陕西省西安市雁塔中路1号