发明名称 一种基于MEMS惯性器件的室内环境精确实时定位方法
摘要 一种基于MEMS惯性器件的室内环境精确实时定位方法,是一种采用MEMS惯性器件通过建立回归模型估算步长获取人员精确定位信息的室内惯性定位算法。本发明在航向角的确定中,利用梯度下降算法补偿姿态误差,增加了行人航向的精度,减小了计算量;本发明还提出利用“三条件判断法”与合加速度融合来判断行人步态的变化,从而提高算法的适应性。该发明不受外界环境干扰,成本低,简单实用,精度高,可应用于室内、隧道、矿井下等特殊场所。
申请公布号 CN105628027A 申请公布日期 2016.06.01
申请号 CN201610094555.6 申请日期 2016.02.19
申请人 中国矿业大学 发明人 李世银;张楠;沈胜强;程瑞良;杨瑞鑫
分类号 G01C21/18(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 重庆百润洪知识产权代理有限公司 50219 代理人 刘立春
主权项 一种基于MEMS惯性器件的室内环境精确实时定位方法,其特征在于:具体步骤如下:步骤A、计算相对位移量,利用公式S=L*N,其中,L为步长,N为步数,具体步骤如下:步骤A1、步长L用行走频率WF和加速度方差AV相关的线性方程表示:<img file="FDA0000926603340000011.GIF" wi="490" he="70" />其中,α、β和γ是步长计算的相关参数,可通过实验训练获得;实验采集行人不同步态下的相关数据,描点划线最后得到一条拟合曲线,即是所需线性方程;步骤A2、加速度的波动可以从一定程度上表征步数N的变化,通过加速度的输出值可检测出行走的步数;航向角可以通过角速度积分获取;步骤B、判断行人步态的变化:同一个线性回归模型并不适用于不同的行走模式,将会带来较大的定位误差,因此在计算步长的同时,应当判断行人的行走模式,识别不同人的步态变化;实验数据表明,静止、行走、跑步三种行走模式下,合加速度方差的数值差异较为明显;通过设定加速度方差的阈值可用来判断行人的行走模式;此外,“零速度检测”阶段可采用基于合加速度、合加速度方差和角速度值的多条件判断方法,精度较高;步骤C、梯度下降算法实现姿态融合,提高航向精确度:包括如下具体步骤:步骤C1、利用加速度计表示四元数误差矩阵:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mi>g</mi></msub><mrow><mo>(</mo><mover><mi>q</mi><mo>^</mo></mover><mo>,</mo><mover><mi>a</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>-</mo><msub><mi>q</mi><mn>0</mn></msub><msub><mi>q</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>a</mi><mi>x</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>0</mn></msub><msub><mi>q</mi><mn>1</mn></msub><mo>+</mo><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>a</mi><mi>y</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>2</mn><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>-</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>a</mi><mi>z</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000926603340000012.GIF" wi="637" he="326" /></maths>步骤C2、对四元数误差矩阵求导,获得雅克比矩阵<img file="FDA0000926603340000013.GIF" wi="159" he="71" />根据公式<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>&dtri;</mo><msub><mi>f</mi><mi>g</mi></msub><mrow><mo>(</mo><mover><mi>q</mi><mo>^</mo></mover><mo>,</mo><mover><mi>a</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><msubsup><mi>J</mi><mi>g</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mover><mi>q</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>*</mo><msub><mi>f</mi><mi>g</mi></msub><mrow><mo>(</mo><mover><mi>q</mi><mo>^</mo></mover><mo>,</mo><mover><mi>a</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000926603340000014.GIF" wi="557" he="85" /></maths>得梯度;步骤C3、求出经过梯度下降法计算的目标姿态:<img file="FDA0000926603340000021.GIF" wi="573" he="174" />其中,μ<sub>t</sub>代表步长,其和实际物体运动的角速度和采样时间相关,在低速运动下,应满足下式<img file="FDA0000926603340000022.GIF" wi="318" he="125" />其中,α为动态变量,β为四元数微分方程求解的姿态算法的收敛速度;由此的最终的姿态融合公式为:<img file="FDA0000926603340000023.GIF" wi="771" he="174" />以上为陀螺仪和加速度计的梯度下降法姿态融合算法;同理也可得陀螺仪和地磁计的姿态融合公式。
地址 221116 江苏省徐州市铜山区大学路1号