CN108051004B - 一种针对四轮机器人的转动瞬心估计方法 - Google Patents
一种针对四轮机器人的转动瞬心估计方法 Download PDFInfo
- Publication number
- CN108051004B CN108051004B CN201711228282.0A CN201711228282A CN108051004B CN 108051004 B CN108051004 B CN 108051004B CN 201711228282 A CN201711228282 A CN 201711228282A CN 108051004 B CN108051004 B CN 108051004B
- Authority
- CN
- China
- Prior art keywords
- terrain
- rotation
- covariance
- state estimation
- posterior
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000008569 process Effects 0.000 claims abstract description 19
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 18
- 239000013598 vector Substances 0.000 claims description 18
- 230000001133 acceleration Effects 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 10
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001052 transient 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
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- 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/005—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 with correlation of navigation data from several sources, e.g. map or contour matching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Manipulator (AREA)
Abstract
本发明公开了一种针对四轮机器人的转动瞬心估计方法,通过引入最优估计理论与地形检测算法,能够获取机器人的转动瞬心的估计值。该方法包括初始化步骤、采集传感器数据步骤、地形检测步骤、调整过程噪声方差步骤、状态预测步骤、状态更新步骤、以及输出转动瞬心的估计值的步骤。本发明优点如下:1)针对运动学方程是非线性而观测方程是线性的情况,在状态预测过程中引入无色变换,以保证在强非线性情况下的状态预测的准确性;2)引入新息协方差估计器;3)由于引入了地形检测,在地形发生明显变化的时候,本发明会调整转动瞬心的过程噪声方差,这种自适应机制能够保证转动瞬心估计的平稳性,同时缩减收敛时间,适用于地形复杂的场景。
Description
技术领域
本发明涉及机器人技术领域,特别是涉及一种针对四轮机器人的转动瞬心估计方法。
背景技术
四轮机器人是一种常见的野外机器人,它可以通过改变左右轮子的速度来控制机器人的方向,具有很好的鲁棒性和灵活性,并且能够实现零半径转向。但是,四轮机器人在移动过程中不可避免的出现打滑现象,这就导致运动学模型的不准确性。
在论文“Mandow A,Morales J,Pedraza S.Approximating Kinematics forTracked Mobile Robots[J].International Journal of Robotics Research,2005,24(10):867-878.”中,作者首次提出了转动瞬心运动学模型,做法是将左右轮子与机器人本体的转动瞬心的位置坐标引入到传统的运动学模型中,以提升运动学模型在存在打滑现象的准确度。
转动瞬心本身可以看作对四轮机器人打滑现象的一种描述方式,而且转动瞬心往往随着地形的变化而变化。因此,如何实时的获取四轮机器人的转动瞬心成为了野外机器人领域中一项重要且具有挑战性的工作。
发明内容
本发明的目的在于提供一种针对四轮机器人的转动瞬心估计方法,以实时获取四轮机器人的转动瞬心。
为此,本发明提供了一种针对四轮机器人的转动瞬心估计方法,包括以下步骤:
步骤一:对采样点序号k、后验状态估计后验误差协方差Pk、新息协方差εk、地形特征向量pk、过程噪声和观测噪声的方差Qk和Rk、采样间隔T以及车身宽度B进行初始化,其中,后验状态估计中的与分别表示东向坐标、北向坐标与航向角的后验状态估计,与为转动瞬心的3个运动学参数的后验状态估计;
步骤二:将采样点序号自增k←k+1,采集加速度计关于垂直于地面轴向的加速度数据,在一个采样周期内按照相等时间间隔采集N次,得到加速度数据集合{ak,i},i=1,…,N;利用朝向地面的摄像头拍摄地面照片,得到像素矩阵Mk;采集左右轮编码器数据,获得左右车轮的旋转速度vL,k与vR,k;采集电子罗盘数据与GPS模块,获得观测向量zk=[ze,k zn,kzθ,k]′,其中ze,k与zn,k为东向坐标与北向坐标的观测值,由GPS模块采集获得,zθ,k为航向角的观测值,由电子罗盘采集获得;
步骤三:根据步骤一中获得的地形特征向量以及步骤二中获得的加速度数据集合和地面照片像素矩阵进行地形检测,判断地形是否发生显著变化;
步骤四:若判定地形发生显著变化,在接下来五个采样点内将转动瞬心的过程噪声方差,乘以设定倍数;若地形没有发生变化,则保持原来的方差;
步骤五:根据步骤一中获得的后验状态估计、后验误差协方差、采样间隔、车轮半径以及车身宽度,步骤二中获得的左右车轮的旋转速度以及步骤四中调整后的过程噪声方差,进行状态预测,得到先验状态估计与先验误差协方差;
本发明与现有技术相比,优点在于:1)针对运动学方程是非线性而观测方程是线性的情况,在状态预测过程中引入无色变换,以保证在强非线性情况下的状态预测的准确性;2)引入新息协方差估计器,在系统噪声统计特性发生或者模型参数发生变化的时候,能够保证估计方法的稳定性;3)由于引入了地形检测,在地形发生明显变化的时候,本发明会调整转动瞬心的过程噪声方差,这种自适应机制能够保证转动瞬心估计的平稳性,同时缩减收敛时间,适用于地形复杂的场景。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明的针对四轮机器人的转动瞬心估计方法的流程图;
图2是根据本发明一实施例的地形相似度仿真结果;以及
图3是根据本发明一实施例的转动瞬心仿真结果。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
本发明提供了一种针对四轮机器人的转动瞬心估计方法,通过引入最优估计理论与地形检测算法,能够获取机器人的转动瞬心的估计值。与已有的技术相比,优点在于:1)针对运动学方程是非线性而观测方程是线性的情况,在状态预测过程中引入无色变换,以保证在强非线性情况下的状态预测的准确性;2)引入新息协方差估计器,在系统噪声统计特性发生或者模型参数发生变化的时候,能够保证估计方法的稳定性;3)由于引入了地形检测,在地形发生明显变化的时候,本发明会调整转动瞬心的过程噪声方差,这种自适应机制能够保证转动瞬心估计的平稳性,同时缩减收敛时间,适用于地形复杂的场景。
如图1所示,本发明的基于卡尔曼滤波的四轮机器人的转动瞬心估计方法的具体实施步骤如下:
S10、初始化
对采样点序号、后验状态估计、后验误差协方差、新息协方差、地形特征向量、过程噪声和观测噪声的方差、采样间隔、车轮半径以及车身宽度进行初始化。具体如下:
采样点序号k=0,后验状态估计中的六个元素需要根据实际情况确定,后验误差协方差Pk=0.1×I6×6,新息协方差εk=O3×3,地形特征向量pk=O8×1,过程噪声和观测噪声的方差Qk和Rk分别为6行6列对角阵和3行3列对角阵,需要根据实际情况确定,采样间隔T,车轮半径Φ,车身宽度B需要根据实际情况确定。其中,下标k表示采样点序号,I6×6为6行6列的单位矩阵,O8×1为8行1列的零向量,O3×3为3行3列的零矩阵,与分别表示东向坐标、北向坐标与航向角的后验状态估计,与为转动瞬心的3个运动学参数的后验状态估计。
在本发明中,单引号上标表示矩阵的转置,例如,C′为矩阵C转置。
S20、采集传感器数据
将采样点序号自增,并采集加速度计、摄像头、左右轮编码器、电子罗盘与GPS模块的数据。具体如下:
采样点序号k←k+1;采集加速度计关于垂直于地面轴向的加速度数据,在一个采样周期内按照相等时间间隔采集N次,得到加速度数据集合{ak,i},i=1,…,N;利用朝向地面的摄像头拍摄地面照片,得到像素矩阵Mk;采集左右轮编码器数据,获得左右车轮的旋转速度vL,k与vR,k;采集电子罗盘数据与GPS模块,获得观测向量zk=[ze,k zn,k zθ,k]′,其中ze,k与zn,k为东向坐标与北向坐标的观测值,由GPS模块采集获得,zθ,k为航向角的观测值,由电子罗盘采集获得。
S30、地形检测
根据步骤S10中获得的地形特征向量以及步骤S20中获得的加速度数据集合和地面照片像素矩阵进行地形检测,判断地形是否发生显著变化。具体如下:
3.4)对地形特征向量进行归一化处理。
3.5)判断地形是否变化:计算地形相似度距离D(pk,pk-1):
S40、调整转动瞬心的过程噪声方差
根据步骤S20中对地形是否发生显著变化的判断,调整转动瞬心的过程噪声方差:如果地形发生显著变化,在接下来五个采样点内将转动瞬心的过程噪声方差,即Qk主对角线第4至6个元素,乘以10倍;如果地形没有发生变化,则保持原来的方差。
S50、状态预测
根据步骤S10中获得的后验状态估计、后验误差协方差、采样间隔、车轮半径以及车身宽度,步骤S20中获得的左右车轮的旋转速度以及步骤S40中调整后的过程噪声方差,进行状态预测,得到先验状态估计与先验误差协方差,具体如下:
S60、状态更新
根据步骤S10中获得的新息协方差、观测噪声方差,步骤S20中获得的观测向量以及步骤S50中获得的先验状态估计和先验误差协方差,进行状态更新,具体如下:
6.5)计算后验误差协方差Pk如下:Pk=(I6×6-KkC)Pk,k-1,其中,矩阵I6×6为6维单位矩阵。
为了验证本发明,我们采用软件MATLAB对进行仿真实验,设置2000个采样点,采样间隔为0.4秒,轮胎半径为35厘米,车架宽度为65厘米,转动瞬心初始为21.6、-21.6、10,在第1001个采样点变为43.2、-43.2、5。同时,我们利用加速度传感器和摄像头分别采集水泥地和草地两种地形的数据各1000组,与MATLAB进行联合实验,模拟移动机器人切换地形。
地形相似度仿真结果如图2所示,可见在地形发生变化的时候,地形相似度距离会发生短暂的跳跃性变化;在地形不发生变化的时候,地形相似度距离会稳定的维持在一个相对较小的值。转动瞬心仿真结果如图3所示,初始值分别设置为65、-65、20,估计结果能够在状态突变后很快的收敛到真实值。由以上仿真可以验证本发明的有效性。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种针对四轮机器人的转动瞬心估计方法,其特征在于,包括以下步骤:
步骤一:对采样点序号k、后验状态估计后验误差协方差Pk、新息协方差εk、地形特征向量pk、过程噪声和观测噪声的方差Qk和Rk、采样间隔f以及车身宽度B进行初始化,其中,后验状态估计中的与分别表示东向坐标、北向坐标与航向角的后验状态估计,与为转动瞬心的3个运动学参数的后验状态估计;
步骤二:将采样点序号自增k←k+1,采集加速度计关于垂直于地面轴向的加速度数据,并在一个采样周期内按照相等时间间隔采集N次,得到加速度数据集合{ak,i},i=1,…,N;利用朝向地面的摄像头拍摄地面照片,得到像素矩阵Kk;采集左右轮编码器数据,获得左右车轮的旋转速度vL,k与vR,k;采集电子罗盘数据与GPS模块数据,获得观测向量zk=[ze,k zn,kzθ,k]′,其中,ze,k与zn,k为东向坐标与北向坐标的观测值,由GPS模块采集获得,zθ,k为航向角的观测值,由电子罗盘采集获得;
步骤三:根据步骤一中获得的地形特征向量以及步骤二中获得的加速度数据集合和地面照片像素矩阵进行地形检测,判断地形是否发生显著变化;
步骤四:若判定地形发生显著变化,在接下来五个采样点内将转动瞬心的过程噪声方差,乘以设定倍数;若地形没有发生变化,则保持原来的方差;
步骤五:根据步骤一中获得的后验状态估计、后验误差协方差、采样间隔以及车身宽度,步骤二中获得的左右车轮的旋转速度以及步骤四中调整后的过程噪声方差,进行状态预测,得到先验状态估计与先验误差协方差;
3.根据权利要求2所述的针对四轮机器人的转动瞬心估计方法,其特征在于,所述步骤四中的设定倍数为10倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711228282.0A CN108051004B (zh) | 2017-11-29 | 2017-11-29 | 一种针对四轮机器人的转动瞬心估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711228282.0A CN108051004B (zh) | 2017-11-29 | 2017-11-29 | 一种针对四轮机器人的转动瞬心估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108051004A CN108051004A (zh) | 2018-05-18 |
CN108051004B true CN108051004B (zh) | 2020-10-09 |
Family
ID=62121361
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711228282.0A Active CN108051004B (zh) | 2017-11-29 | 2017-11-29 | 一种针对四轮机器人的转动瞬心估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108051004B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109249429B (zh) * | 2018-09-25 | 2019-10-01 | 安徽果力智能科技有限公司 | 一种双足机器人地形分类系统 |
CN110160527B (zh) * | 2019-05-06 | 2020-08-28 | 安徽红蝠智能科技有限公司 | 一种移动机器人导航方法与装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706342A (zh) * | 2012-05-31 | 2012-10-03 | 重庆邮电大学 | 一种智能移动机器人的定位与环境建模方法 |
CN104154916A (zh) * | 2013-08-30 | 2014-11-19 | 北京航天发射技术研究所 | 一种基于激光陀螺捷联惯组的车载定位设备 |
CN105573310A (zh) * | 2014-10-11 | 2016-05-11 | 北京自动化控制设备研究所 | 一种煤矿巷道机器人定位与环境建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013034561A1 (en) * | 2011-09-06 | 2013-03-14 | Land Rover | Suspension control device |
-
2017
- 2017-11-29 CN CN201711228282.0A patent/CN108051004B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706342A (zh) * | 2012-05-31 | 2012-10-03 | 重庆邮电大学 | 一种智能移动机器人的定位与环境建模方法 |
CN104154916A (zh) * | 2013-08-30 | 2014-11-19 | 北京航天发射技术研究所 | 一种基于激光陀螺捷联惯组的车载定位设备 |
CN105573310A (zh) * | 2014-10-11 | 2016-05-11 | 北京自动化控制设备研究所 | 一种煤矿巷道机器人定位与环境建模方法 |
Non-Patent Citations (1)
Title |
---|
基于自适应卡尔曼滤波的侧滑移动机器人运动模型估计;吴耀等;《电子与信息学报》;20151231;第37卷(第12期);第3016-3204页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108051004A (zh) | 2018-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107991110B (zh) | 一种履带式机器人滑动参数检测方法 | |
CN108731670B (zh) | 基于量测模型优化的惯性/视觉里程计组合导航定位方法 | |
CN105021184B (zh) | 一种用于移动平台下视觉着舰导航的位姿估计系统及方法 | |
CN104732518B (zh) | 一种基于智能机器人地面特征的ptam改进方法 | |
CN107341814B (zh) | 基于稀疏直接法的四旋翼无人机单目视觉测程方法 | |
CN110702107A (zh) | 一种单目视觉惯性组合的定位导航方法 | |
CN105180937B (zh) | 一种mems‑imu初始对准方法 | |
CN108731700B (zh) | 一种视觉惯性里程计中的加权欧拉预积分方法 | |
CN110517324B (zh) | 基于变分贝叶斯自适应算法的双目vio实现方法 | |
CN108051004B (zh) | 一种针对四轮机器人的转动瞬心估计方法 | |
CN108036789B (zh) | 一种野外机器人航迹推算方法 | |
CN105004337B (zh) | 基于直线匹配的农用无人机自主导航方法 | |
CN112258409A (zh) | 一种用于无人驾驶的单目相机绝对尺度恢复方法及装置 | |
CN114234967A (zh) | 一种基于多传感器融合的六足机器人定位方法 | |
CN112388635B (zh) | 机器人多传感器融合感知与空间定位的方法、系统及装置 | |
CN109443353B (zh) | 基于模糊自适应ickf的视觉-惯性紧耦合组合导航方法 | |
CN108020855B (zh) | 一种滑移转向机器人的位姿和转动瞬心联合估计方法 | |
CN107123128B (zh) | 一种保证准确性的车辆运动状态估计方法 | |
CN117272525A (zh) | 一种智能电动汽车路面附着系数估计方法 | |
CN113569393B (zh) | 一种稳健的绝对定向算法 | |
CN108088498B (zh) | 一种航向角与转动阻力系统联合估计方法 | |
CN114638858B (zh) | 一种基于车载双相机系统的运动目标位置与速度估计方法 | |
CN116338719A (zh) | 基于b样条函数的激光雷达-惯性-车辆融合定位方法 | |
EP3410074B1 (en) | Method and device for improving performance of relative-position sensor, and computer storage medium | |
CN109674480A (zh) | 一种基于改进互补滤波的人体运动姿态解算方法 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220429 Address after: 241000 first floor, building 5, standardized workshop, Anhui Nanling Economic Development Zone, Nanling County, Wuhu City, Anhui Province Patentee after: Anhui Zhongfen Intelligent Technology Co.,Ltd. Address before: 230088 room C-104, National University Science Park, 602 Huangshan Road, high tech Zone, Hefei City, Anhui Province Patentee before: ANHUI 11TONG INFORMATION TECHNOLOGY Co.,Ltd. |