发明名称 一种空中微小抛撒物运动参数测量方法及系统
摘要 一种用于空中微小抛撒物运动的测量方法及系统,其步骤为:(1)将由三轴微陀螺、三轴微加速度计和三轴微磁强计构成的微测量单元安装于抛撒物中,测量并记录抛撒物从抛撒时刻到落地静止后一段时间内的传感器数据;(2)以抛撒物落地静止后的位置作为相对的初始位置,对抛撒物落地静止后一段时间的传感器数据进行计算,得到抛撒物落地静止后的姿态和三轴微陀螺的零偏;(3)以抛撒物飞行过程中记录的三轴微陀螺和三轴微加速度计数据,用逆时间惯性导航法得到抛撒物在整个飞行过程中的实时运动参数。测量系统包括微测量单元、非易失存储器以及嵌入式处理器。本发明具有结构简单紧凑、原理简单、测量精度高、可靠性好、抗干扰能力强等优点。
申请公布号 CN101644579B 申请公布日期 2011.08.03
申请号 CN200910044123.4 申请日期 2009.08.17
申请人 中国人民解放军国防科学技术大学 发明人 李涛;齐振恒;吴美平;吴文启;曹聚亮;罗兵;练军想
分类号 G01C21/18(2006.01)I;G01C21/08(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 湖南兆弘专利事务所 43008 代理人 赵洪;周长清
主权项 1.一种用于空中微小抛撒物运动的测量方法,其特征在于步骤为:(1)、将由三轴微陀螺、三轴微加速度计和三轴微磁强计构成的微测量单元安装于抛撒物中,测量并记录抛撒物从抛撒时刻到落地静止后一段时间内的传感器数据;读取传感器数据写入非易失存储器中;(2)、以抛撒物落地静止后的位置作为相对的初始位置,采用三轴微磁强计和三轴微加速度计构成静态数字罗盘,对抛撒物落地静止后一段时间的传感器数据进行计算,得到抛撒物落地静止后的姿态,并计算三轴微陀螺的零偏;即,读取非易失存储器中的传感器数据,依据读取的非易失存储器中的传感器数据得到抛撒物的落地静止时刻;依据落地静止时刻将传感器数据分成飞行段和落地静止段;根据落地静止段的陀螺数据计算出陀螺的零偏并对陀螺数据进行修正;计算抛撒物落地静止后的位置、速度且根据静止的三轴微加速度计和三轴微磁强计算静止后的姿态;(3)、在已知落地静止后运动参数的基础上,根据抛撒物飞行段的三轴微陀螺和三轴微加速度计数据,采用逆时间惯性导航法得到抛撒物在整个飞行过程中的各种实时运动参数;所述逆时间惯性导航法,是在得到静止后运动参数的的基础上,根据飞行阶段的三轴微陀螺和三轴微加速度计的数据,进行逆时间惯性导航计算,得到抛撒物的运动参数;其具体步骤为:a.确定实时姿态:采用式(12)的四元数更新进行计算<img file="FSB00000417968900011.GIF" wi="1259" he="91" />其中,q<sub>k</sub>表示T<sub>k</sub>时刻的姿态四元数,<img file="FSB00000417968900012.GIF" wi="67" he="55" />为由三轴微陀螺测量值计算得到的从T<sub>k</sub>到T<sub>k+1</sub>时刻的旋转四元数,<sup>*</sup>表示四元数的共轭,o表示四元数相乘;b.确定实时速度:采用式(18)计算<maths num="0001"><![CDATA[<math><mrow><msub><mi>V</mi><mi>k</mi></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><mrow><mo>(</mo><msubsup><mi>C</mi><mi>b</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>A</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><msup><mi>g</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mi>&Delta;t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中V<sub>k</sub>表示T<sub>k</sub>时刻的速度,<img file="FSB00000417968900014.GIF" wi="203" he="54" />为T<sub>k+1</sub>时刻的方向余弦矩阵,A<sub>f</sub>(T<sub>k+1</sub>)为T<sub>k+1</sub>时刻的三轴微加速度计的测量值,g<sup>n</sup>(k+1)为T<sub>k+1</sub>时刻的重力矢量,Δt为采样步长;c.确定实时位置:采用式(21)计算P<sub>k</sub>=P<sub>k+1</sub>-V<sub>k+1</sub>*Δt    (21)其中,P(T<sub>k</sub>)表示T<sub>k</sub>时刻的位置向量。
地址 410073 湖南省长沙市砚瓦池正街47号中国人民解放军国防科学技术大学机电工程与自动化学院313室