[发明专利]一种基于多传感器信息融合的智能车定位方法有效

专利信息
申请号: 201710131878.2 申请日: 2017-03-07
公开(公告)号: CN106840179B 公开(公告)日: 2019-12-10
发明(设计)人: 祝辉;邓炯;余彪;梁华为 申请(专利权)人: 中国科学院合肥物质科学研究院
主分类号: G01C21/32 分类号: G01C21/32;G01C21/34;G01S19/45
代理公司: 34126 合肥市科融知识产权代理事务所(普通合伙) 代理人: 陈思聪
地址: 230031 *** 国省代码: 安徽;34
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 基于 传感器 信息 融合 智能 定位 方法
【权利要求书】:

1.一种基于多传感器信息融合的智能车定位方法,其特征在于:包括以下步骤:

(1)、基于智能车的光电编码器和激光雷达采集的数据,并利用同步定位与地图构建方法SLAM得到智能车在SLAM方法构建的局部地图中的位置信息;

步骤(1)的具体过程如下:

(1.1)、以智能车启动时车体中心为原点Oo,沿车前进方向为Xo轴方向,垂直Xo轴指向车体左侧为Yo轴方向,Zo轴垂直向上建立里程计坐标系;由于只考虑智能车在水平面内运动且没有俯仰和翻滚动作,因此Z轴坐标始终为0且只有偏航角的变化,因此可以得到智能车的位姿模型(x,y,θ);通过智能车后轮左、右轮的光电编码器,推算智能车里程计信息,具体过程如下:

(1.1.1)、0时刻,智能车在里程计坐标系中的位姿(x0,y00)为(0,0,0);

(1.1.2)、单位采样时间内,在里程计坐标系下,车轮位置的增量ΔS通过以下公式得到:

其中ΔQ为光电码盘输出的脉冲增量,D为车轮直径,Ts为采样时间,N为光电码盘上的光栅总数,K为光电编码器减速率;

(1.1.3)、左右轮之间的距离为ω,单位采样时间内左、右轮位置的增量分别为ΔSL和ΔSR,则汽车从k–1时刻Sk-1=(xk-1,yk-1k-1)到k时刻Sk=(xk,ykk)的位姿变化为:

xk=xk-1+ΔDk cos(θk-1+Δθk),

yk=yk-1+ΔDksin(θk-1+Δθk),

θk=θk-1+Δθk

(1.1.4)、使用上述方法累积得到时刻k智能车的里程计信息;

(1.2)、通过车载激光雷达获得并记录时刻k智能车周边环境的观测信息;

(1.3)、建立与里程计坐标系同原点同轴的局部地图坐标系,并在局部地图坐标系下通过SLAM算法来进行构图和定位;通过步骤(1.1)得到的里程计信息作为SLAM算法的控制输入信息,通过步骤(1.2)得到的雷达扫描信息作为SLAM算法观测信息,使用一种基于粒子滤波器的SLAM算法得到智能车在地图坐标系下的优化位姿,同时根据优化后的位姿将k时刻智能车周边的障碍物信息加入到局部地图中;

使用粒子滤波器进行优化以得到智能车的更优位姿与地图时,智能车当前时刻的位姿根据前一时刻的位姿来推算得到的;知道k-1时刻智能车的位姿xk-1,k时刻的控制信息即里程计信息uk,k时刻的观测信息zk,估计智能车的当前位姿的xk的后验概率;

在粒子滤波器中采用序贯重要性采样算法和采样重要性重采样算法来获得后验概率分布,通过对参考分布的采样粒子进行加权来近似得到后验概率分布,

具体包括以下步骤:

(1.3.1)、初始化系统,设置N个粒子,每个粒子代表智能车当前的位姿,采样计算权重

(1.3.2)、重要性采样,在时刻k,根据k时刻智能车的控制信息uk推算k时刻智能车位姿的分布,对每个粒子进行采样

(1.3.3)、根据k时刻的观测信息zk和粒子前一时刻的权值,计算k时刻粒子重要性权值并进行归一化权值;

(1.3.4)、重采样,增加重要性权值高的粒子,删除重要性权值低的粒子,计算后验概率;

(2)、通过智能车的GPS定位获得智能车在WGS-84坐标系中的位置信息;

(3)、将步骤(1)和步骤(2)得到的智能车在局部地图坐标系中、WGS-84坐标系中的位置信息通过坐标转换统一在大地坐标系中;

(4)、基于步骤(3)的坐标转换过程,采用卡尔曼滤波技术将步骤(1)得到的智能车的位置信息、步骤(2)得到的智能车的位置信息进行融合,得到智能车精准的位置估计;

卡尔曼滤波模型中,步骤(1)通过SLAM算法得到的智能车在大地坐标系下的位姿作为预测值;步骤(2)通过GPS得到的智能车在大地坐标系下的位姿作为观测值;

定义状态空间转移方程和状态空间的观测方程如下:

X(k)=AX(k-1)+W(k),

Z(k)=HX(k)+V(k),

其中X(k)为k时刻的位姿状态向量,A为过程矩阵,H为测量矩阵,W(k)为过程噪声矩阵,其协方差为Q,V(k)为测量噪声,其协方差为R;卡尔曼滤波的过程如下:

预测方程组:

X(k|k-1)=AX(k-1|k-1),

P(k|k-1)=AP(k-1|k-1)A′+Q,

更新方程组:

X(k|k)=X(k|k-1)+Kg(k)(Z(k)-HX(k|k-1)),

P(k|k)=(1-Kg(k)H)P(k|k-1),

其中:

X(k|k-1)是k时刻的预估位姿状态向量;

X(k-1|k-1)是k-1时刻的最优位姿状态;

Kg(k)是k时刻的卡尔曼增益;

X(k|k)是k时刻的最优位姿状态估计;

Q是系统过程噪声W(k)的协方差;

R是系统测量噪声V(k)的协方差。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院合肥物质科学研究院,未经中国科学院合肥物质科学研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201710131878.2/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top