CN111256727A - 一种基于Augmented EKF的提高里程计精度的方法 - Google Patents
一种基于Augmented EKF的提高里程计精度的方法 Download PDFInfo
- Publication number
- CN111256727A CN111256727A CN202010101440.1A CN202010101440A CN111256727A CN 111256727 A CN111256727 A CN 111256727A CN 202010101440 A CN202010101440 A CN 202010101440A CN 111256727 A CN111256727 A CN 111256727A
- Authority
- CN
- China
- Prior art keywords
- state quantity
- agv
- odometer
- observation
- matrix
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000003190 augmentative effect Effects 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 230000007704 transition Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C22/00—Measuring distance traversed on the ground by vehicles, persons, animals or other moving solid bodies, e.g. using odometers, using pedometers
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
- Measurement Of Distances Traversed On The Ground (AREA)
Abstract
本发明提出了一种基于Augmented EKF的提高里程计精度的方法,包括如下步骤:S1、将移动机器人AGV的初始状态量Xk‑1以及初始状态量Xk‑1的协方差矩阵Pk‑1输入至预测模型中,获得预测状态量和预测状态量的协方差矩阵S2、将预测状态量和预测状态量的协方差矩阵作为更新过程的输入,输入至更新模型中,获得更新状态量S3、将获得的更新状态量导入AGV运动学模型,里程计根据动态的AGV物理模型参数进行航位推算。通过将AGV的物理模型参数拟合到卡尔曼滤波中,再通过扩展卡尔曼滤波融合全局定位的信息,计算出AGV物理模型参数,确保里程计在进行航位推算时,能够动态校准AGV物理模型参数,减少噪声的干扰,提高了里程计航位推算的精确度。
Description
技术领域
本发明涉及动态系统估算技术领域,更具体地说是涉及一种基于Augmented EKF的提高里程计精度的方法。
背景技术
目前,里程计作为对移动机器人AGV进行航位推算的主要传感器,具有稳定性高、成本低的特点;但是,AGV中里程计的精度主要取决于对AGV物理模型的估算,而在AGV硬件制作工艺和环境温度导致的硬件形变等因素的限制下,其物理模型实际上是动态变化的,传统的基于里程计的AGV航位推算精度受限于AGV错误的物理模型,导致航位推算结果随着时间的累计误差不断增大,极大的降低了对AGV定位的精度。
卡尔曼滤波是一种高效的自回归滤波,能够在存在诸多不确定性情况的组合信息中估计动态系统的状态,非常适合不断变化的系统,且内存占用较小、速度快,是实时问题和嵌入式系统的理想选择。
因此,如何将AGV的物理模型参数拟合到卡尔曼滤波中,根据AGV的全局定位动态估算物理模型的参数值,以期提高里程计的精度,是本领域技术人员亟待解决的问题。
发明内容
有鉴于此,本发明提供了一种基于Augmented EKF的提高里程计精度的方法,通过将AGV的物理模型参数拟合到卡尔曼滤波中,再通过扩展卡尔曼滤波融合全局定位的信息,计算出AGV物理模型参数,以期提高里程计精度。
为了实现上述目的,本发明采用如下技术方案:
一种基于Augmented EKF的提高里程计精度的方法,包括如下步骤:
具体的,设定初始状态量Xk-1=[x y θ rL rR d];其中,x,y为AGV的初始位置坐标值,θ为AGV的偏转角度值,rL和rR分别表示AGV两个左轮和右轮的半径,d表示AGV左轮和右轮之间的距离。
具体的,预测阶段的卡尔曼滤波的动态方程为其中,k为时间点,vk为在k时间点获得的里程计线速度,wk为在k时间点获得的里程计的转向速度;预测状态量的不确定性用协方差更新方程表示,其中,为预测状态量的协方差,A为状态转换Jacobian矩阵,Q为过程噪声矩为阵,在AGV状态转移的过程中,除了状态具有不确定性外,移动也具有不确定性,将移动的不确定性定义为过程噪声矩阵Q,该噪声源于外部干扰和系统误差,例如AGV左轮或者右轮打滑等,总体的不确定度即过程不确定度之和。
具体的,所述初始状态量Xk-1=[x y θ rL rR d]在卡尔曼滤波中的动态方程为:
具体的,在每次循环中的Jacobian矩阵A为:
其中,
wL为左轮的角速度,wR为右轮角速度,和具体的,AGV的角度与偏转角度值θ直接通过全局观测获得,当输入观测结果时,执行更新过程;观测结果和状态量之间的转换公式为:其中,H为状态量和观测结果的转换关系,z为观测状态,u为观测噪声。状态观测的物理意义为并非所有的状态都可以观测到,传感器观测状态量存在局限,例如,状态量中包含距离和速度,如果传感器观测距离,那么速度就无法观测,另外,在观测中还存在尺度的变化,因此观测结果和状态量之间需要进行转换。
具体的,所述状态量和观测结果的转换关系H=[I3 O3];其中,I3为单位矩阵,O3为零矩阵。
具体的,将所述状态转换Jacobian矩阵A和状态量和观测结果的转换关系H带入更新模型中,动态获得rL、rR和d的值。
本发明具有以下有益效果:
经由上述的技术方案可知,基于现有技术,本发明提出了一种基于Augmented EKF的提高里程计精度的方法,确保里程计在进行航位推算时,能够动态校准AGV物理模型参数,减少噪声的干扰,提高了里程计航位推算的精确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为基于Augmented EKF的提高里程计精度的方法流程图;
图2为雷达与地图完全匹配的初始图;
图3(a)(b)校准前后只基于里程计的3次360度逆时针旋转的航位推算的对比图,图3(a)为未经过校准的只基于里程计的3次360度逆时针旋转的航位推算图,图3(b)为经过校准的只基于里程计的3次360度逆时针旋转的航位推算图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了如下方法:
一种基于Augmented EKF的提高里程计精度的方法,该方法分为预算阶段和更新阶段;
当有里程计数据输入时执行预测阶段,预测阶段过程如下:
首先设定初始状态量Xk-1=[x y θ rL rR d];其中,x,y为AGV的初始位置坐标值,θ为AGV的偏转角度值,rL和rR分别表示AGV两个左轮和右轮的半径,d表示AGV左轮和右轮之间的距离。
具体的,初始状态量Xk-1=[x y θ rL rR d]在卡尔曼滤波中的动态方程为:
卡尔曼滤波中所有状态均为高斯分布,分布的中心值就是状态值估计值,不确定性即为协方差矩阵,因此,预测状态量的不确定性用协方差更新方程表示,其中,为预测状态量的协方差,A为状态转换Jacobian矩阵,Q为过程噪声矩阵,该噪声源于外部干扰和系统误差,例如AGV轮子打滑,由于预测过程符合全概率法则,因此,总体的不确定度为过程不确定度的叠加。
在每次循环中的Jacobian矩阵A为:
其中,
当有全局观测值输入时执行更新步骤,更新阶段步骤如下:
AGV的位置(x,y)与偏转角度值θ直接通过全局观测获得,全局观测的位置和偏转角度值可以从自适应蒙特卡洛粒子滤波(AMCL)中导出,也可以从其他全局定位方法,例如超宽带UWB,或者基于路标的定位系统中导出。当输入观测结果时,执行更新过程;观测结果和状态量之间的转换公式为:其中,H为状态量和观测结果的转换关系,z为观测状态,u为观测噪声。具体的,状态观测的物理意义即不能观测到AGV的全部状态,例如状态量中包含距离和速度两个量,如果传感器观测距离量,则无法观测速度量,另外还存在观测尺度上的变化,因此观测结果和状态量之间需要转换,为该过程的转换公式。
其中,状态量和观测结果的转换关系H=[I3 O3];其中,I3为单位矩阵,O3为零矩阵。
S3、将获得的AGV物理模型参数rL、rR和d的值导入AGV运动学模型,里程计根据动态的AGV物理模型参数进行航位推算。
如图2所示,为雷达与地图完全匹配的初始图,白色点运为雷达读数。
如图3(a)(b)所示,为校准前后的只基于里程计的3次360度逆时针旋转的航位推算图,图3(a)为未经过校准的只基于里程计的3次360度逆时针旋转的航位推算图,图3(b)为经过校准的只基于里程计的3次360度逆时针旋转的航位推算图。很显然校准前如图3(a)所示,里程计对AGV的航位推算与初始图2有较大偏差,而经卡尔曼滤波校准后,如图3(b)所示,里程计对AGV的航位推算的校准较精确。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
2.根据权利要求1所述的一种基于Augmented EKF的提高里程计精度的方法,其特征在于,设定初始状态量Xk-1=[x y θ rL rR d];其中,x,y为AGV的初始位置坐标值,θ为AGV的偏转角度值,rL和rR分别表示AGV两个左轮和右轮的半径,d表示AGV左轮和右轮之间的距离。
7.根据权利要求6所述的一种基于Augmented EKF的提高里程计精度的方法,其特征在于,所述状态量和观测结果的转换关系H=[I3 O3];其中,I3为单位矩阵,O3为零矩阵。
9.根据权利要求8所述的一种基于Augmented EKF的提高里程计精度的方法,其特征在于,将所述状态转换Jacobian矩阵A和状态量和观测结果的转换关系H带入更新模型中,动态获得rL、rR和d的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010101440.1A CN111256727B (zh) | 2020-02-19 | 2020-02-19 | 一种基于Augmented EKF的提高里程计精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010101440.1A CN111256727B (zh) | 2020-02-19 | 2020-02-19 | 一种基于Augmented EKF的提高里程计精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111256727A true CN111256727A (zh) | 2020-06-09 |
CN111256727B CN111256727B (zh) | 2022-03-08 |
Family
ID=70952933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010101440.1A Active CN111256727B (zh) | 2020-02-19 | 2020-02-19 | 一种基于Augmented EKF的提高里程计精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111256727B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112082546A (zh) * | 2020-08-20 | 2020-12-15 | 北京自动化控制设备研究所 | 光纤管道测量装置数据后处理方法 |
CN112683268A (zh) * | 2020-12-08 | 2021-04-20 | 中国铁建重工集团股份有限公司 | 一种基于扩展卡尔曼滤波的巷道实时定位导航方法及系统 |
CN114462545A (zh) * | 2022-02-15 | 2022-05-10 | 集美大学 | 一种基于语义slam的地图构建方法及装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102556075A (zh) * | 2011-12-15 | 2012-07-11 | 东南大学 | 一种基于改进扩展卡尔曼滤波的车辆运行状态估计方法 |
CN104908548A (zh) * | 2014-03-11 | 2015-09-16 | 丰田自动车株式会社 | 车辆状态推定装置、车辆控制装置及车辆状态推定方法 |
DE102016011366A1 (de) * | 2016-09-21 | 2017-04-06 | Daimler Ag | Verfahren zur Positionsbestimmung eines Fahrzeugs |
CN108955688A (zh) * | 2018-07-12 | 2018-12-07 | 苏州大学 | 双轮差速移动机器人定位方法及系统 |
CN109596128A (zh) * | 2019-01-14 | 2019-04-09 | 哈尔滨工程大学 | 一种基于多水听器提高多auv协同定位性能的方法 |
CN109974697A (zh) * | 2019-03-21 | 2019-07-05 | 中国船舶重工集团公司第七0七研究所 | 一种基于惯性系统的高精度测绘方法 |
CN109993113A (zh) * | 2019-03-29 | 2019-07-09 | 东北大学 | 一种基于rgb-d和imu信息融合的位姿估计方法 |
CN110068334A (zh) * | 2019-04-23 | 2019-07-30 | 科罗玛特自动化科技(苏州)有限公司 | 一种磁导航agv的高精度定位方法 |
CN110602647A (zh) * | 2019-09-11 | 2019-12-20 | 江南大学 | 基于扩展卡尔曼滤波和粒子滤波的室内融合定位方法 |
-
2020
- 2020-02-19 CN CN202010101440.1A patent/CN111256727B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102556075A (zh) * | 2011-12-15 | 2012-07-11 | 东南大学 | 一种基于改进扩展卡尔曼滤波的车辆运行状态估计方法 |
CN104908548A (zh) * | 2014-03-11 | 2015-09-16 | 丰田自动车株式会社 | 车辆状态推定装置、车辆控制装置及车辆状态推定方法 |
DE102016011366A1 (de) * | 2016-09-21 | 2017-04-06 | Daimler Ag | Verfahren zur Positionsbestimmung eines Fahrzeugs |
CN108955688A (zh) * | 2018-07-12 | 2018-12-07 | 苏州大学 | 双轮差速移动机器人定位方法及系统 |
CN109596128A (zh) * | 2019-01-14 | 2019-04-09 | 哈尔滨工程大学 | 一种基于多水听器提高多auv协同定位性能的方法 |
CN109974697A (zh) * | 2019-03-21 | 2019-07-05 | 中国船舶重工集团公司第七0七研究所 | 一种基于惯性系统的高精度测绘方法 |
CN109993113A (zh) * | 2019-03-29 | 2019-07-09 | 东北大学 | 一种基于rgb-d和imu信息融合的位姿估计方法 |
CN110068334A (zh) * | 2019-04-23 | 2019-07-30 | 科罗玛特自动化科技(苏州)有限公司 | 一种磁导航agv的高精度定位方法 |
CN110602647A (zh) * | 2019-09-11 | 2019-12-20 | 江南大学 | 基于扩展卡尔曼滤波和粒子滤波的室内融合定位方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112082546A (zh) * | 2020-08-20 | 2020-12-15 | 北京自动化控制设备研究所 | 光纤管道测量装置数据后处理方法 |
CN112082546B (zh) * | 2020-08-20 | 2023-01-10 | 北京自动化控制设备研究所 | 光纤管道测量装置数据后处理方法 |
CN112683268A (zh) * | 2020-12-08 | 2021-04-20 | 中国铁建重工集团股份有限公司 | 一种基于扩展卡尔曼滤波的巷道实时定位导航方法及系统 |
CN114462545A (zh) * | 2022-02-15 | 2022-05-10 | 集美大学 | 一种基于语义slam的地图构建方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111256727B (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111256727B (zh) | 一种基于Augmented EKF的提高里程计精度的方法 | |
CN102620886B (zh) | 两步在轨辨识组合航天器转动惯量估计方法 | |
CN109459019B (zh) | 一种基于级联自适应鲁棒联邦滤波的车载导航计算方法 | |
CN112013836B (zh) | 一种基于改进自适应卡尔曼滤波的航姿参考系统算法 | |
CN109885883B (zh) | 一种基于gk聚类算法模型预测的无人车横向运动的控制方法 | |
CN101819041A (zh) | 自进化anfis与ukf结合的gps/mems-ins组合定位误差动态预测方法 | |
CN110515381B (zh) | 用于定位机器人的多传感器融合算法 | |
CN103776449B (zh) | 一种提高鲁棒性的动基座初始对准方法 | |
CN110174907B (zh) | 一种基于自适应卡尔曼滤波的人体目标跟随方法 | |
CN113280808A (zh) | 一种移动机器人定位精度提升方法及系统 | |
Qian et al. | IMM-UKF based land-vehicle navigation with low-cost GPS/INS | |
Guo et al. | A stochastic model-based fusion algorithm for enhanced localization of land vehicles | |
Lin et al. | Combined adaptive robust Kalman filter algorithm | |
CN110243363B (zh) | 一种基于低成本imu与rfid技术结合的agv实时定位方法 | |
CN113306573B (zh) | 一种自动驾驶车辆学习型路径跟踪预测控制方法 | |
CN108680162B (zh) | 一种基于渐进无迹卡尔曼滤波的人体目标跟踪方法 | |
CN106871905B (zh) | 一种非理想条件下高斯滤波替代框架组合导航方法 | |
CN117246359A (zh) | 车辆轨迹的预测方法、装置、设备、存储介质及产品 | |
CN110672127A (zh) | 阵列式mems磁传感器实时标定方法 | |
Krantz et al. | Non-uniform dead-reckoning position estimate updates | |
CN103616027B (zh) | 一种基于改进msd的重力匹配方法 | |
CN105549003A (zh) | 一种汽车雷达目标跟踪方法 | |
CN113158415B (zh) | 一种基于误差分析的车辆轨迹相似度评估方法 | |
CN114154231A (zh) | 一种基于深度强化学习的车辆行驶状态估计系统及方法 | |
Lai et al. | Adaptive fuzzy neural network-aided progressive Gaussian approximate filter for GPS/INS integration navigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |