发明名称 基于加速度计和陀螺仪的载体转速与转角补偿测量方法
摘要 本发明是一种基于加速度计和陀螺仪的载体转速与转角补偿测量方法,属于惯性传感器与测量技术领域。本发明方法同时采用三轴加速度计和单轴陀螺仪对载体转动的转角进行估计,充分利用其在测量频域上的互补性,得到准确的载体转角估计结果。同时,通过该结果计算得到的载体转速,相比于直接利用陀螺仪测量得到的载体转速,可以保持与载体转角持续的一致性。本方法的传感器可以安装于载体旋转平面上的任意位置,不依赖于载体确定的旋转中心和旋转半径,因此不需要对载体进行物理尺寸、旋转中心、旋转半径的测量,并可对旋转中心和旋转半径变化的运动载体进行测量,这极大地简化了测量的复杂性,提高了该方法的适用范围。
申请公布号 CN103913171A 申请公布日期 2014.07.09
申请号 CN201410175383.6 申请日期 2014.04.28
申请人 北京理工大学 发明人 王卫江;高巍;荣政
分类号 G01C21/16(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 基于加速度计和陀螺仪的载体转速与转角补偿测量方法,其特征在于:具体包括如下步骤:第一步:将传感器和相关测量电路安装于载体任意位置上,使得加速度计x、y轴平面与载体转动平面平行即可,确定载体的初始转角并记为θ<sub>est</sub>(0);所述传感器包括三轴加速度计、单轴陀螺仪;所述相关测量电路包括信号采集电路和处理电路,采集电路与传感器连接,处理电路与采集电路连接;第二步:确定传感器数据采样周期T和互补滤波器数字滤波方法中的互补滤波器时间常数τ,令载体开始转动,测量电路以采样周期T读取三轴加速度计和单轴陀螺仪的测量值,当前时刻n的测量值包括加速度计x轴方向上的加速度值a<sub>x</sub>(n)、y轴方向上的加速度值a<sub>y</sub>(n)、陀螺仪角速度测量值rate(n),采用移动平均方法对采集到的加速度计和陀螺仪测量值进行平稳去噪;第三步:利用加速度计测量值a<sub>x</sub>(n)和a<sub>y</sub>(n)计算加速度计所测载体转角θ<sub>acc</sub>(n),利用陀螺仪测量值rate(n)计算陀螺仪所测载体转角θ<sub>gyro</sub>(n);具体方法为:所述的利用加速度计测量值计算载体转角是利用重力加速度在x轴和y轴上分量的三角函数关系求得,当载体在非水平面内做旋转运动时,重力加速度在加速度计的x轴和y轴上产生分量,在不考虑其他外力情况下,利用该x轴和y轴上的加速度分量根据三角关系求得:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mrow><mi>tan</mi><mi>&theta;</mi></mrow><mi>acc</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000497669130000011.GIF" wi="413" he="144" /></maths>继而得到:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>acc</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>arctan</mi><mfrac><mrow><msub><mi>a</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000497669130000012.GIF" wi="468" he="148" /></maths>所述的利用陀螺仪测量值计算载体转角是利用转动角速度对时间积分求得;利用上一时刻转角估计值以及当前时刻角速度测量值,并将积分运算离散化为求和运算以后可得:θ<sub>gyro</sub>(n)=θ<sub>est</sub>(n‑1)+rate(n)×T;第四步:根据传感器数据采样周期T和互补滤波器时间常数τ计算滤波器权重系数W,采用互补滤波器数字滤波方法计算载体转角估计值θ<sub>est</sub>(n);互补滤波器数字滤波方法的具体过程为:根据传感器数据采样周期和互补滤波器时间常数τ计算滤波器权重系数W,由该滤波器时间常数<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&tau;</mi><mo>=</mo><mfrac><mrow><mi>W</mi><mo>&times;</mo><mi>T</mi></mrow><mrow><mn>1</mn><mo>-</mo><mi>W</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000497669130000013.GIF" wi="218" he="126" /></maths>得到:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>W</mi><mo>=</mo><mfrac><mi>&tau;</mi><mrow><mi>&tau;</mi><mo>+</mo><mi>T</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000497669130000021.GIF" wi="243" he="116" /></maths>再利用互补滤波器将加速度计测量计算所得转角θ<sub>acc</sub>(n)和陀螺仪测量计算所得转角θ<sub>gyro</sub>(n)权重相加,得到载体转角估计值:θ<sub>est</sub>(n)=W×θ<sub>gyro</sub>(n)+(1‑W)×θ<sub>acc</sub>(n);第五步:利用当前时刻载体转角估计值θ<sub>est</sub>(n)、上一时刻载体转角估计值θ<sub>est</sub>(n‑1)和采样周期T计算载体转动角速度ω(n);所述的载体转动角速度ω(n)的求解是利用载体转角的估计值对时间微分得到,将微分运算离散化后得到:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>est</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&theta;</mi><mi>est</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mi>T</mi></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0000497669130000022.GIF" wi="549" he="127" /></maths>
地址 100081 北京市海淀区中关村南大街5号