CN104977001A - 一种应用于个人室内导航***的相对导航方法 - Google Patents
一种应用于个人室内导航***的相对导航方法 Download PDFInfo
- Publication number
- CN104977001A CN104977001A CN201410130948.9A CN201410130948A CN104977001A CN 104977001 A CN104977001 A CN 104977001A CN 201410130948 A CN201410130948 A CN 201410130948A CN 104977001 A CN104977001 A CN 104977001A
- Authority
- CN
- China
- Prior art keywords
- navigation
- accelerometer
- speed
- equation
- error
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明属于导航方法,具体涉及一种应用于个人室内导航***的相对导航方法。它包括:步骤一:初始化,步骤二:惯性相对导航解算,计算速度、位置信息;步骤三:保存必要信息;步骤四:计算投影值积分;步骤五:检测是否静止;步骤六:修正输出;对速度、位置、加速度进行修正,修正后的数据作为本方法的输出,并循环执行步骤二~步骤六,直到导航结束。本发明的效果是:本发明的方法可以准确有效的解算出实验人员行进的轨迹,说明本方法正确、有效。
Description
技术领域
本发明属于导航方法,具体涉及一种应用于个人室内导航***的相对导航方法。
背景技术
由于卫星导航所具有的高精度、低功耗、易携带等独特优势,其在日常生活中的应用日益广泛,然而在室内、地下及工作厂房等密闭空间中由于卫星信号受遮挡卫星用户机不能定位而无法满足用户的导航需求。
个人室内导航***由于其不依赖卫星导航、完全自主的独特优势而受到越来越多的重视,以美国为首的欧美国家对该***的研发起步较早,但尚未出现货架产品,为了填补国内对该技术与产品的研发空白,对个人室内导航技术进行研究。
发明内容
本发明的目的是针对现有技术的缺陷,提供一种应用于个人室内导航***的相对导航方法。
本发明是这样实现的:一种应用于个人室内导航***的相对导航方法,包括下述步骤:
步骤一:初始化
对速度、位置、姿态进行初始化,在第一次计算时使用初始化数据,在后续循环计算中使用外部传感器给出的实时信息;
步骤二:惯性相对导航解算
计算速度、位置信息;
步骤三:保存必要信息
保存加速度计陀螺敏感信息以及加速度计信息在导航坐标系中投影值
步骤四:计算投影值积分
计算加速度计、陀螺敏感信息均方差以及加速度计信息在导航坐标系中投影值积分;
步骤五:检测是否静止
当***同时满足下述三个条件时,判定***静止
1)三轴加速度计敏感信息的均方差均小于0.05m/s^2;
2)三轴陀螺敏感信息的均方差均小于0.001rad/s;
3)一个滤波周期内加速度积分绝对值均小于0.01m/s;
如果判断结果是***处于静止状态,则执行下面步骤六,否则执行步骤二;
步骤六:修正输出
对速度、位置、加速度进行修正,修正后的数据作为本方法的输出,并循环执行步骤二~步骤六,直到导航结束。
如上所述的一种应用于个人室内导航***的相对导航方法,其中,所述的步骤一的初始化包括给定初始化数据,包括
a)初始位置:X=0,Y=0,Z=0;
b)初始速度:Vx=0,Vy=0,Vz=0;
c)初始姿态:俯仰角滚动角航向角Yaw=0.0,其中的asin表示反正弦函数
其中为1秒内加速度的平均值,计算公式为:
其中,和是加速度计敏感参数 的两个输出值,该加速度计敏感参数由外部的加速度计给出。
如上所述的一种应用于个人室内导航***的相对导航方法,其中,所述的步骤二惯性相对导航解算,包括
a)姿态更新
设表征载体姿态的四元数:Q=[q0 q1 q2 q3]T,四元数更新的解析式为:
其中:I为4×4阶单位矩阵; 为k到k+1时刻陀螺敏感到的角增量,由下式获得: 其中Tn=0.005s,即导航解算周期,姿态矩阵的更新计算公式为:
b)速度、位置更新
相对导航速度更新算法为:
相对导航位置更新算法为:
P(K+1)=P(K)+0.5*(V(K+1)+V(K))*Tn (4)
其中为加速度计敏感信息在导航系中投影:V(K)、P(K)为K时刻的速度、位置信息,所述的Tn为导航周期,一般取50毫秒。
如上所述的一种应用于个人室内导航***的相对导航方法,其中,所述的步骤六修正输出中,使用的卡尔曼滤波***方程为
其中,公式(8)中的第一个方程为***误差状态方程,具体形式见步骤(6.1),X(t)为***误差状态量X(t)=[δP,δX,φ,▽,ε]T,F(t)为***状态转移矩阵 W(t)为***噪声矢量,g为当地重力加速度,
公式(8)中的第二个方程为***观测方程,具体形式见步骤(6.2), 为零速观测量,H=[03×3 I3×3 03×3 03×3 03×3],V(t)为观测噪声矢量,
(6.1)***误差方程
a)位置误差方程
b)速度误差方程
c)姿态误差方程
d)加速度计零位误差方程
e)陀螺漂移误差方程
其中,δP=[δPx δPy δPz]T、δV=[δVx δVy δVz]T为位置、速度误差,φ=[φx φy φz]T为姿态误差,▽=[▽x ▽y ▽z]T为机体系加速度计零位误差,ε=[εx εy εz]T为机体系陀螺漂移误差,
(6.2)零速匹配观测方程
在***检测到人体处于静止状态下,利用下式求观测量:
其中,Z为零速匹配卡尔曼滤波器观测量,为零速状态下滤波周期内速度平均值,在此设定卡尔曼滤波周期为50ms,即10个导航周期滤波一次,
观测方程:
本发明的效果是:本发明的方法可以准确有效的解算出实验人员行进的轨迹,说明本方法正确、有效。
附图说明
附图1是工作人员进行实验时行走路线的示意图。
附图2是利用本方法导航得到的轨迹图。
具体实施方式
一种应用于个人室内导航***的相对导航方法,包括下述步骤:
步骤一:初始化
给定初始化数据,包括
d)初始位置:X=0,Y=0,Z=0;
e)初始速度:Vx=0,Vy=0,Vz=0;
f)初始姿态:俯仰角滚动角航向角Yaw=0.0。其中的asin表示反正弦函数
其中为1秒内加速度的平均值,计算公式为:
其中,和是加速度计敏感参数 的两个输出值,该加速度计敏感参数由外部的加速度计给出。初始化之后,在后续计算过程中使用到的速度、姿态、加速度等信息,均由外部相应传感器给出。
步骤二:惯性相对导航解算
c)姿态更新
设表征载体姿态的四元数:Q=[q0 q1 q2 q3]T,四元数更新的解析式为:
其中:I为4×4阶单位矩阵; 为k到k+1时刻陀螺敏感到的角增量,由下式获得: 其中Tn=0.005s,即导航解算周期,姿态矩阵的更新计算公式为:
即用公式(1)计算q0~q3,计算出这个周期的q0~q3,继而得到这个周期的姿态矩阵。下个周期的姿态矩阵再利用这个公式计算。
d)速度、位置更新
相对导航速度更新算法为:
相对导航位置更新算法为:
P(K+1)=P(K)+0.5*(V(K+1)+V(K))*Tn (4)
其中为加速度计敏感信息在导航系中投影:V(K)、P(K)为K时刻的速度、位置信息。所述的Tn为导航周期,一般取50毫秒。
步骤三:保存必要信息
保存加速度计陀螺敏感信息以及加速度计信息在导航坐标系中投影值三类信息。其中加速度计陀螺敏感信息由外部给出,加速度计信息在导航坐标系中投影值由步骤三计算得到。
步骤四:计算投影值积分
计算加速度计、陀螺敏感信息均方差以及加速度计信息在导航坐标系中投影值积分。其中投影值是用原始信息与姿态矩阵的转置矩阵相乘得到,积分是用投影值积分得到的。
步骤五:检测是否静止
由于个人室内导航***采用的是微观性测量单元器件精度较差,短时间内纯惯性导航导航误差发散较快,因此需要通过一定的修正算法抑制误差的快速发散,达到高精度导航的目的,在个人室内导航导航***中,我们充分利用人体行走过程中脚部运动动静交替切换的特点,在脚部处于静止的时间段内,采用零速修正算法对***各项导航误差进行修正。
零速修正的前提是静止状态的检测,零速检测算法需要既能够满足隔离***的运动状态,又能够对震动、扰动等特殊运动状态有一定容忍度,因此检测算法设计如下:
1)每个导航周期实时滑动保存1个滤波周期内(10个导航周期)三轴加速度计陀螺敏感信息以及三轴加速度计信息在导航坐标系中投影
2)每个导航周期实时计算三轴加速度计、陀螺敏感信息的均方差及一个滤波周期内三轴加速度计信息在导航坐标系中投影的积分值。其中,加速度均方差计算公式为:
角速度均方差计算公式为:
三轴加速度计信息在导航坐标系中投影的积分值计算公式为:
如果同时满足以下三个条件,则可以视该滤波周期内***处于静止状态,否则***处于运动状态:
4)三轴加速度计敏感信息的均方差均小于0.05m/s^2;
5)三轴陀螺敏感信息的均方差均小于0.001rad/s;
6)一个滤波周期内加速度积分绝对值均小于0.01m/s;
如果判断结果是***处于静止状态,则执行下面步骤六,否则执行步骤二。
步骤六:修正输出
卡尔曼滤波***方程
其中,公式(8)中的第一个方程为***误差状态方程,具体形式见步骤(6.1),X(t)为***误差状态量X(t)=[δP,δX,φ,▽,ε]T,F(t)为***状态转移矩阵 W(t)为***噪声矢量,g为当地重力加速度。
公式(8)中的第二个方程为***观测方程,具体形式见步骤(6.2), 为零速观测量,H=[03×3 I3×3 03×3 03×3 03×3],V(t)为观测噪声矢量。
(6.1)***误差方程
f)位置误差方程
g)速度误差方程
h)姿态误差方程
i)加速度计零位误差方程
j)陀螺漂移误差方程
其中,δP=[δPx δPy δPz]T、δV=[δVx δVy δVz]T为位置、速度误差,φ=[φx φy φz]T为姿态误差,▽=[▽x ▽y ▽z]T为机体系加速度计零位误差,ε=[εx εy εz]T为机体系陀螺漂移误差,
(6.2)零速匹配观测方程
在***检测到人体处于静止状态下,利用下式求观测量:
其中,Z为零速匹配卡尔曼滤波器观测量,为零速状态下滤波周期内速度平均值,在此设定卡尔曼滤波周期为50ms,即10个导航周期滤波一次。
观测方程:
本领域的技术人员利用本领域的卡尔曼滤波基本知识,结合本申请公开的***方程、***误差方程和零速匹配观测方程可以计算出速度、位置、加速度的误差,然后再利用本领域的通用修正方法,可以得到修正后的速度、位置、加速度。所述修正后的速度、位置、加速度即为本方法的输出值。但是因为导航时持续进行的,因此重复执行步骤二~步骤六,直到导航结束。
Claims (4)
1.一种应用于个人室内导航***的相对导航方法,其特征在于,包括下述步骤:
步骤一:初始化
对速度、位置、姿态进行初始化,在第一次计算时使用初始化数据,在后续循环计算中使用外部传感器给出的实时信息;
步骤二:惯性相对导航解算
计算速度、位置信息;
步骤三:保存必要信息
保存加速度计陀螺敏感信息以及加速度计信息在导航坐标系中投影值
步骤四:计算投影值积分
计算加速度计、陀螺敏感信息均方差以及加速度计信息在导航坐标系中投影值积分;
步骤五:检测是否静止
当***同时满足下述三个条件时,判定***静止
1)三轴加速度计敏感信息的均方差均小于0.05m/s^2;
2)三轴陀螺敏感信息的均方差均小于0.001rad/s;
3)一个滤波周期内加速度积分绝对值均小于0.01m/s;
如果判断结果是***处于静止状态,则执行下面步骤六,否则执行步骤二;
步骤六:修正输出
对速度、位置、加速度进行修正,修正后的数据作为本方法的输出, 并循环执行步骤二~步骤六,直到导航结束。
2.如权利要求1所述的一种应用于个人室内导航***的相对导航方法,其特征在于:所述的步骤一的初始化包括给定初始化数据,包括
a)初始位置:X=0,Y=0,Z=0;
b)初始速度:Vx=0,Vy=0,Vz=0;
c)初始姿态:俯仰角滚动角 航向角Yaw=0.0,其中的asin表示反正弦函数
其中为1秒内加速度的平均值,计算公式为:
其中,和是加速度计敏感参数的两个输出值,该加速度计敏感参数由外部的加速度计给出。
3.如权利要求2所述的一种应用于个人室内导航***的相对导航方法,其特征在于:所述的步骤二惯性相对导航解算,包括
a)姿态更新
设表征载体姿态的四元数:Q=[q0 q1 q2 q3]T,四元数更新的解析式为:
其中:I为4×4阶单位矩阵; 为k到k+1时刻陀螺敏感到的角增量,由下式获得:其中Tn=0.005s,即导航解算周 期,姿态矩阵的更新计算公式为:
b)速度、位置更新
相对导航速度更新算法为:
相对导航位置更新算法为:
P(K+1)=P(K)+0.5*(V(K+1)+V(K))*Tn (4)
其中为加速度计敏感信息在导航系中投影:V(K)、P(K)为K时刻的速度、位置信息,所述的Tn为导航周期,一般取50毫秒。
4.如权利要求3所述的一种应用于个人室内导航***的相对导航方法,其特征在于:所述的步骤六修正输出中,使用的卡尔曼滤波***方程为
其中,公式(8)中的第一个方程为***误差状态方程,具体形式见步骤(6.1),X(t)为***误差状态量X(t)=[δP,δX,φ,▽,ε]T,F(t)为***状态转移矩阵 W(t)为***噪声矢量,g为当地重力加速度,
公式(8)中的第二个方程为***观测方程,具体形式见步骤(6.2), 为零速观测量,H=[03×3 I3×3 03×3 03×3 03×3],V(t)为观测噪声矢量,
(6.1)***误差方程
a)位置误差方程
b)速度误差方程
c)姿态误差方程
d)加速度计零位误差方程
e)陀螺漂移误差方程
其中,δP=[δPx δPy δPz]T、δV=[δVx δVy δVz]T为位置、速度误差,φ=[φx φy φz]T为姿态误差,▽=[▽x ▽y ▽z]T为机体系加速度计零位误差,ε=[εx εy εz]T为机体系陀螺漂移误差,
(6.2)零速匹配观测方程
在***检测到人体处于静止状态下,利用下式求观测量:
其中,Z为零速匹配卡尔曼滤波器观测量,为零速状态下滤波周期内速度平均值,在此设定卡尔曼滤波周期为50ms,即10个导航周期滤波一次,
观测方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130948.9A CN104977001A (zh) | 2014-04-02 | 2014-04-02 | 一种应用于个人室内导航***的相对导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130948.9A CN104977001A (zh) | 2014-04-02 | 2014-04-02 | 一种应用于个人室内导航***的相对导航方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104977001A true CN104977001A (zh) | 2015-10-14 |
Family
ID=54273740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410130948.9A Pending CN104977001A (zh) | 2014-04-02 | 2014-04-02 | 一种应用于个人室内导航***的相对导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104977001A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106289247A (zh) * | 2016-07-26 | 2017-01-04 | 北京长城电子装备有限责任公司 | 基于惯性传感器的室内定位装置 |
CN107345813A (zh) * | 2017-07-07 | 2017-11-14 | 江苏奥斯威尔信息科技有限公司 | 一种基于mt‑pdr与光强信息的室内平面图构建方法 |
CN108088464A (zh) * | 2016-11-22 | 2018-05-29 | 北京自动化控制设备研究所 | 一种闭环修正安装误差的挠曲估计方法 |
CN112798010A (zh) * | 2019-11-13 | 2021-05-14 | 北京三快在线科技有限公司 | 一种视觉惯性里程计vio***的初始化方法、装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445200A (zh) * | 2011-09-30 | 2012-05-09 | 南京理工大学 | 微小型个人组合导航***及其导航定位方法 |
CN102607562A (zh) * | 2012-04-12 | 2012-07-25 | 南京航空航天大学 | 基于载体飞行模态判别的微惯性参数自适应姿态确定方法 |
CN103512575A (zh) * | 2012-06-26 | 2014-01-15 | 北京自动化控制设备研究所 | 一种测绘车用惯导***零速修正方法 |
CN103616030A (zh) * | 2013-11-15 | 2014-03-05 | 哈尔滨工程大学 | 基于捷联惯导解算和零速校正的自主导航***定位方法 |
-
2014
- 2014-04-02 CN CN201410130948.9A patent/CN104977001A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445200A (zh) * | 2011-09-30 | 2012-05-09 | 南京理工大学 | 微小型个人组合导航***及其导航定位方法 |
CN102607562A (zh) * | 2012-04-12 | 2012-07-25 | 南京航空航天大学 | 基于载体飞行模态判别的微惯性参数自适应姿态确定方法 |
CN103512575A (zh) * | 2012-06-26 | 2014-01-15 | 北京自动化控制设备研究所 | 一种测绘车用惯导***零速修正方法 |
CN103616030A (zh) * | 2013-11-15 | 2014-03-05 | 哈尔滨工程大学 | 基于捷联惯导解算和零速校正的自主导航***定位方法 |
Non-Patent Citations (2)
Title |
---|
张金亮等: "基于MEMS惯性技术的鞋式个人导航***", 《惯性技术发展动态发展方向研讨会文集》 * |
殷红: "基于foot-mounted的IMU室内行人航迹推算研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106289247A (zh) * | 2016-07-26 | 2017-01-04 | 北京长城电子装备有限责任公司 | 基于惯性传感器的室内定位装置 |
CN108088464A (zh) * | 2016-11-22 | 2018-05-29 | 北京自动化控制设备研究所 | 一种闭环修正安装误差的挠曲估计方法 |
CN108088464B (zh) * | 2016-11-22 | 2021-07-13 | 北京自动化控制设备研究所 | 一种闭环修正安装误差的挠曲估计方法 |
CN107345813A (zh) * | 2017-07-07 | 2017-11-14 | 江苏奥斯威尔信息科技有限公司 | 一种基于mt‑pdr与光强信息的室内平面图构建方法 |
CN107345813B (zh) * | 2017-07-07 | 2020-08-04 | 江苏奥斯威尔信息科技有限公司 | 一种基于mt-pdr与光强信息的室内平面图构建方法 |
CN112798010A (zh) * | 2019-11-13 | 2021-05-14 | 北京三快在线科技有限公司 | 一种视觉惯性里程计vio***的初始化方法、装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ludwig et al. | Comparison of Euler estimate using extended Kalman filter, Madgwick and Mahony on quadcopter flight data | |
Rotella et al. | State estimation for a humanoid robot | |
Simanek et al. | Evaluation of the EKF-based estimation architectures for data fusion in mobile robots | |
CN110398245B (zh) | 基于脚戴式惯性测量单元的室内行人导航姿态估计方法 | |
CN104764451A (zh) | 一种基于惯性和地磁传感器的目标姿态跟踪方法 | |
CN105865452B (zh) | 一种基于间接卡尔曼滤波的移动平台位姿估计方法 | |
CN101846510A (zh) | 一种基于星敏感器和陀螺的高精度卫星姿态确定方法 | |
CN102116634B (zh) | 一种着陆深空天体探测器的降维自主导航方法 | |
CN103575299A (zh) | 利用外观测信息的双轴旋转惯导***对准及误差修正方法 | |
CN101706287A (zh) | 一种基于数字高通滤波的旋转捷联***现场标定方法 | |
CN104977001A (zh) | 一种应用于个人室内导航***的相对导航方法 | |
Quan et al. | Interlaced optimal-REQUEST and unscented Kalman filtering for attitude determination | |
Gao et al. | Invariant filtering for legged humanoid locomotion on a dynamic rigid surface | |
CN102114918B (zh) | 一种基于多速率敏感器组合定姿的姿控反馈回路 | |
Yang et al. | Quaternion-based Kalman filtering on INS/GPS | |
CN104613964A (zh) | 一种跟踪脚部运动特征的步行者定位方法及*** | |
Bai et al. | A novel plug-and-play factor graph method for asynchronous absolute/relative measurements fusion in multisensor positioning | |
Tang et al. | Exploring the accuracy potential of IMU preintegration in factor graph optimization | |
CN108592902B (zh) | 一种基于多传感器的定位设备及定位方法、***和机械臂 | |
Atchuthan et al. | Odometry based on auto-calibrating inertial measurement unit attached to the feet | |
Lamy-Perbal et al. | An improved pedestrian inertial navigation system for indoor environments | |
Zhou et al. | Vision-based state estimation and trajectory tracking control of car-like mobile robots with wheel skidding and slipping | |
Blachuta et al. | Attitude and heading reference system based on 3D complementary filter | |
Lubbe et al. | State estimation for a hexapod robot | |
Fernando et al. | Observability analysis of position estimation for quadrotors with modified dynamics and range measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151014 |