CN110926483B - 一种用于自动驾驶的低成本传感器组合定位***及方法 - Google Patents
一种用于自动驾驶的低成本传感器组合定位***及方法 Download PDFInfo
- Publication number
- CN110926483B CN110926483B CN201911168614.XA CN201911168614A CN110926483B CN 110926483 B CN110926483 B CN 110926483B CN 201911168614 A CN201911168614 A CN 201911168614A CN 110926483 B CN110926483 B CN 110926483B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- gps
- state
- kalman filter
- entering
- 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
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/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
-
- 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
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
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)
- Navigation (AREA)
Abstract
本发明公开了一种用于自动驾驶的低成本传感器组合定位***,通过采用低成本GPS、IMU传感器,根据低成本GPS的定位特性设计一套状态机,用于动态调整卡尔曼滤波器参数,然后通过卡尔曼滤波器的预测和状态更新模块输出车辆当前的位置以及速度和航向信息。本发明还提供了一种用于自动驾驶的低成本传感器组合定位方法。本发明不仅降低定位套件的硬件成本,便于自动驾驶市场化;而且最终使整个***定位精度可达亚米级。
Description
技术领域
本发明属于自动驾驶领域,特别涉及一种用于自动驾驶的低成本传感器组合定位***及方法。
背景技术
组合导航是自动驾驶领域的的重要研究课题,目前常用价格高昂的高精度GPS和IMU。RTK(高精度GPS)套件设备价格数万元;此外其移动终端通过通信模块与基站和数据处理中心处理卫星信号,在通信环境恶劣地方,其无法有效工作。高精度IMU由磁力计、陀螺仪、加速度计组成,可获取当前车辆的姿态信息。工业级IMU的价格在数千到数万元之间。常用组合导航设备因传感器价格昂贵无法,无法实现自动驾驶技术市场化。
发明内容
发明目的:本发明针对现有技术存在的问题,提出了一种成本低,定位效果好的用于自动驾驶的低成本传感器组合定位***。
技术方案:为实现上述目的,本发明提供了一种用于自动驾驶的低成本传感器组合定位***,包括车辆速度信息采集单元、车辆航向信息采集单元、车辆位置信息采集单元和车辆位置定位单元,其中,车辆速度信息采集单元用于采集车辆的实时车速信息;车辆航向信息采集单元用于实时采集车辆的航向信息;车辆位置信息采集单元为GPS,用于实时采集车辆位置信息;所述车辆定位单元中包括状态机和卡尔曼滤波器,所述状态机根据GPS反馈的参数对卡尔曼滤波器的参数进行调整,所述卡尔曼滤波器根据车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元输入的车辆的信息对车辆进行定位。
进一步,所述状态机中包括重启状态、GPS无效状态、GPS有效状态和停车/缓慢行驶状态,其中,重启状态时,初始化卡尔曼滤波器参数中的***协方差矩阵P、测量噪声协方差矩阵R和***状态量矩阵X;GPS无效状态时,GPS采集的位置信息不输入卡尔曼滤波器;GPS有效状态时,GPS采集的位置信息输入卡尔曼滤波器;停车/缓慢行驶状态时,将GPS的位置信息每间隔1秒输入卡尔曼滤波器一次;当***上电后,进入重启状态,在重启状态时,当hdop>h_mid或者star_num<10时候,进入GPS无效状态;当hdop<=h_mid且v_car>v_mid时,进入GPS有效状态;当hdop<=h_mid且v_car<=v_mid时,进入停车/缓慢行驶状态;在GPS无效状态时,若hdop<h_good,进入GPS有效状态;若P异常,进入重启状态;在GPS有效状态时,若hdop>h_bad,进入GPS无效状态;若v_car<v_ba,则进入停车/缓慢行驶状态;若P异常,进入重启状态;在停车/缓慢行驶状态时,若v_car>v_good,进入GPS有效状态;若hdop>h_bad,进入GPS无效状态;若P异常,进入重启状态;hdop为水平精度因子、star_num为接收星数、v_car为从ODOM读取的当前车速、P为卡尔曼滤波器的***协方差矩阵,h_good、h_mid、h_bad为判断GPS定位精度参数,v_good、v_bad、v_mid为判断车辆速度参数。这样能够有效的提高定位的精度。
其中,h_good的范围为0.55-0.65,h_mid的范围为0.68-0.73,h_bad的范围为0.75-0.77。v_good的范围为3-5m/s,v_mid的范围为1-2m/s,v_bad的范围为0.3_0.5m/s。这样有效的保证了整个方法的鲁棒性,保证了状态见切换的合理性,使定位效果更好。
进一步,车辆速度信息采集单元采用设置在车辆上的ODOM,ODOM采集车辆的前进速度信息和旋转速度;车辆航向信息采集单元采用装在车辆旋转中心位置的IMU,其主要采集和提供车辆的航向信息;车辆位置信息采集单元采用设置在车顶并且位于车辆旋转中心位置的GPS,其主要采集和提供车辆的位置信息。
本发明还提供了一种基于上述用于自动驾驶的低成本传感器组合定位***的组合定位方法,其包括以下步骤:
步骤1:车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元分别将采集到车辆的速度、车辆的航向信息和车辆的位置信息发动到车辆位置定位单元;
步骤2:在车辆位置定位单元构建卡尔曼滤波器,车辆位置定位单元中的状态机根据获得的车辆位置信息,初始化或者修改卡尔曼滤波器参数;
步骤3:卡尔曼滤波器根据接收到的车辆的车辆的速度、车辆的航向信息和车辆的位置信息,结合其中的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出。
进一步,其中卡尔曼滤波器参数包括***协方差矩阵P、测量噪声协方差矩阵R、***状态量矩阵X(x,y,v,yaw,yawv,vscale,yawbias)和观测值矩阵其中,x表示输入卡尔曼滤波器的车辆在UTM坐标系下的横坐标,y表示输入卡尔曼滤波器的车辆在UTM坐标系下的纵坐标,v表示输入卡尔曼滤波器的车辆的速度,yaw表示输入卡尔曼滤波器的车辆的航向角,yawv表示输入卡尔曼滤波器的车辆的航向角速度,vscale表示输入卡尔曼滤波器的速度比例因子,yawbias表示输入卡尔曼滤波器的航向角偏差;zx和zy表示GPS输出的位置信息,zyaw是IMU输出的航向角,zv为ODMO输出的速度信息,zyawv是IMU输出的航向角的角速度。速度比例因子vscale解决由车辆轮胎气压、车辆载重、路面不平整等造成的测量速度和真实速度不一致的情况,引入航向角偏差yawbias解决道路磁场环境变换引起的航向角不稳定的情况,这两个状态量能够根据GPS的观测数据和卡尔曼滤波算法的预测模型动态调整。更加有效的提高对车辆定位的精度。
进一步,所述预测模型为:
xt=xt-1+vt*cos(yawt)*Δt;
yt=yt-1+vt*sin(yawt)*Δt;
上式中vt表示在t时刻车辆的速度,yawt表示车辆在t时刻的航向角,xt、yt分别表示在t时刻车辆在UTM坐标系下的横坐标和纵坐标,vt-1表示t-1时刻车辆的速度,yawt-1表示t-1时刻车辆的航向角,xt-1和yt-1分别代表在t-1时刻车辆在UTM坐标系下的横坐标和纵坐标,Δt表示时间差,表示t-1时刻车辆的速度比例因子,表示t-1时刻车辆的航向角偏差,为t-1时刻航向角的角速度。
进一步,所述步骤2中还包括用GPS采集的其航向信息修正车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息;根据公式:
yawoffset=yawgps-yawimu;
yawinput=yawimu+yawoffest;
对车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息进行修正,式中yawgps为GPS星数大于等于12,车速大于5m/s、水平精度因子小于0.6时,RMC报文输出的航向角度,yawimu为imu输出的航向角,yawoffset为偏移角,yawinput为输入卡尔曼滤波器的yaw值。这样能够有效的提高***的定位精度。
进一步,所述步骤2中初始化或者修改卡尔曼滤波器参数的方法为:设定重启状态、GPS无效状态、GPS有效状态和停车/缓慢行驶状态,其中,重启状态时,初始化卡尔曼滤波器参数中的***协方差矩阵P、测量噪声协方差矩阵R和***状态量矩阵X;GPS无效状态时,GPS采集的位置信息不输入卡尔曼滤波器;GPS有效状态时,GPS采集的位置信息输入卡尔曼滤波器;停车/缓慢行驶状态时,将GPS的位置信息每间隔1秒输入卡尔曼滤波器一次;当***上电后,进入重启状态,在重启状态时,当hdop>h_mid或者star_num<10时候,进入GPS无效状态;当hdop<=h_mid且v_car>v_mid时,进入GPS有效状态;当hdop<=h_mid且v_car<=v_mid时,进入停车/缓慢行驶状态;在GPS无效状态时,若hdop<h_good,进入GPS有效状态;若P异常,进入重启状态;在GPS有效状态时,若hdop>h_bad,进入GPS无效状态;若v_car<v_ba,则进入停车/缓慢行驶状态;若P异常,进入重启状态;在停车/缓慢行驶状态时,若v_car>v_good,进入GPS有效状态;若hdop>h_bad,进入GPS无效状态;若P异常,进入重启状态;hdop为水平精度因子、star_num为接收星数、v_car为从ODOM读取的当前车速、P为卡尔曼滤波器的***协方差矩阵,h_good、h_mid、h_bad为判断GPS定位精度参数,v_good、v_bad、v_mid为判断车辆速度参数。
工作原理:本发明提出的定位方案采用低成本GPS、IMU传感器,其中IMU输出的航向信息(Yaw)、GPS提供的位置信息(Position)、车辆ODOM提供的速度(Velocity)和角速度(Yaw_rate)输入扩展卡尔曼滤波器,根据低成本GPS的定位特性设计一套状态机,用于动态调整卡尔曼滤波器参数,然后通过卡尔曼滤波器的预测和状态更新模块输出车辆当前的位置以及速度和航向信息。车辆的ODOM信息可从车辆的CAN总线读取,其速度和角速度值分别由安装在左右车轮的脉冲计数器求和平均和做差求出。本发明还加入速度比例因子和航向角偏差用于实时修正车速和航向角;同时,在车辆速度较大、GPS定位稳定时,用GPS的航向信息修正IMU的航向信息。
有益效果:与现有技术相比,本发明利采用的低成本GPS和IMU,降低定位套件的硬件成本,便于自动驾驶市场化;同时,通过状态机调参机制,最终***定位精度可达亚米级。
附图说明
图1为本发明的***示意图;
图2为状态机工作流程示意图;
图3为本发明与现有技术的车辆行驶轨迹对比图。
具体实施方式
下面将结合本发明实例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种用于自动驾驶的低成本传感器组合定位***,包括车辆速度信息采集单元、车辆航向信息采集单元、车辆位置信息采集单元和车辆位置定位单元,其中车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元将采集到的信息发送到车辆位置定位单元,车辆位置定位单元采用引入了速度比例因子vscale和航向角偏差yawbias的卡尔曼滤波器对接收到的车辆的信息进行融合,从而得到车辆当前的定位信息,定位信息包括车辆当前的速度,航向角、在UTM平面坐标系下的横坐标和纵坐标。
其中,车辆速度信息采集单元采用设置在车辆上的ODOM,ODOM采集车辆的前进速度信息和旋转速度,其从车上Can总线上可读取,自动驾驶车辆都含有此信息,其通过左右轮脉冲计数器求得;车辆航向信息采集单元采用装在车辆旋转中心位置的IMU,其主要采集和提供车辆的航向信息,本实施例中采用mti-3型号的IMU;车辆位置信息采集单元采用设置在车顶并且位于车辆旋转中心位置的GPS,其主要采集和提供车辆的位置信息,即车辆所在位置的经度和纬度,本实施例中采用GN200B型号的GPS。
本实施例公开的采用上述用于自动驾驶的低成本传感器的定位***的组合定位方法,具体包括以下步骤:
步骤1:车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元分别将采集到车辆的速度、车辆的航向信息和车辆的位置信息发动到车辆位置定位单元;
步骤2:在车辆位置定位单元构建卡尔曼滤波器,车辆位置定位单元根据获得的车辆位置信息,初始化或者修改卡尔曼滤波器参数,其中卡尔曼滤波器参数包括***协方差矩阵测量噪声协方差矩阵***状态量矩阵X(x,y,v,yaw,yawv,vscale,yawbias)和观测值矩阵Z(zx,zy,zyaw,zv),其中,Px_x表示车辆在UTM坐标系下的横坐标的方差值,px_y表示车辆在UTM坐标系下的横坐标与纵坐标的协方差值,px_v表示车辆在UTM坐标系下的横坐标与车辆速度的协方差值,Px_yaw表示车辆在UTM坐标系下的横坐标与航向角的协方差值,表示车辆在UTM坐标系下的横坐标与航向角速度的协方差值,表示车辆在UTM坐标系下的横坐标与速度比例因子的协方差值,表示车辆在UTM坐标系下的横坐标与航向角偏差的协方差值,Py_x表示车辆在UTM坐标系下的纵坐标与横坐标的协方差值,py_y表示车辆在UTM坐标系下的纵坐标的方差值,py_v表示车辆在UTM坐标系下的纵坐标与车辆速度的协方差值,py_yaw表示车辆在UTM坐标系下的纵坐标与航向角的协方差值,表示车辆在UTM坐标系下的纵坐标与航向角速度的协方差值,表示车辆在UTM坐标系下的纵坐标与速度比例因子的协方差值,表示车辆在UTM坐标系下的纵坐标与航向角偏差的协方差值,Pv_x表示车辆的速度与在UTM坐标系下横坐标的协方差值,pv_y表示车辆的速度与在UTM坐标系下纵坐标的协方差值,pv_v表示车辆速度的方差值,pv_yaw表示车辆速度与航向角的协方差值,表示车辆速度与航向角速度的协方差值,表示车辆速度与速度比例因子的协方差值,表示车辆速度与航向角偏差的协方差值,Pyaw_x表示车辆的航向角与车辆在UTM坐标系下的横坐标的协方差值,pyaw_y表示车辆的航向角与车辆在UTM坐标系下的纵坐标的协方差值,pyaw_v表示车辆的航向角与车辆速度的协方差值,pyaw_yaw表示车辆航向角的方差值,表示车辆航向角与航向角速度的协方差值,pyaw_vsacle表示车辆航向角与速度比例因子的协方差值,表示车辆航向角与航向角偏差的协方差值,表示车辆航向角速度与车辆在UTM坐标系下的横坐标的协方差值,表示车辆航向角速度与车辆在UTM坐标系下的纵坐标的协方差值,表示车辆航向角速度与车辆速度的协方差值,表示车辆航向角速度与车辆航向角的协方差值,表示车辆航向角速度的方差值,表示车辆航向角速度与速度比例因子的协方差值,表示车辆航向角速度与航向角偏差的协方差值,表示车辆的速度比例因子与车辆在UTM坐标系下的横坐标的协方差值,表示车辆的速度比例因子与车辆在UTM坐标系下的纵坐标的协方差值,表示车辆的速度比例因子与车辆速度的协方差值,表示车辆的速度比例因子与车辆航向角的协方差值,表示车辆的速度比例因子与车辆航向角速度的协方差值,表示车辆的速度比例因子的方差值,表示车辆的速度比例因子与车辆航向角偏差的协方差值,表示车辆航向角偏差与车辆在UTM坐标系下的横坐标的协方差值,表示车辆航向角偏差与车辆在UTM坐标系下的纵坐标的协方差值,表示车辆航向角偏差与车辆速度的协方差值,表示车辆航向角偏差与车辆航向角的协方差值,表示车辆航向角偏差与车辆航向角速度的协方差值,表示车辆航向角偏差与车辆的速度比例因子的协方差值,表示车辆航向角偏差的方差值;rx_x表示GPS在横向位置的方差,ry_y表示GPS在纵向位置的方差,rv_v表示车辆速度的方差,ryaw_yaw表示车辆航向角的方差,表示车辆航向角速度的方差,表示车辆的速度比例因子的方差,表示车辆航向角偏差的方差,其中,矩阵R中的元素是根据传感器的特性设定的方差,其与传感器的精度相关,和的值为0;x表示输入卡尔曼滤波器的车辆在UTM坐标系下的横坐标,y表示输入卡尔曼滤波器的车辆在UTM坐标系下的纵坐标,v表示输入卡尔曼滤波器的车辆的速度,yaw表示输入卡尔曼滤波器的车辆的航向角,yawv表示输入卡尔曼滤波器的车辆的航向角速度,vscale表示输入卡尔曼滤波器的速度比例因子,yawbias表示输入卡尔曼滤波器的航向角偏差;zx和zy分别表示GPS输出的位置信息在UTM坐标系下的横坐标和纵坐标,zyaw是IMU输出的航向角,zv为ODMO输出的速度信息,是IMU输出的航向角的角速度。
步骤3:卡尔曼滤波器根据接收到的车辆的车辆的速度、车辆的航向信息和车辆的位置信息,结合其中的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出。具体包括如下步骤:
步骤301:根据卡尔曼滤波器中的CTRV预测模型分别预测下一时刻的***状态量X中每个量的值;其中,CTRV预测模型为:
xt=xt-1+vt*cos(yawt)*Δt;
yt=yt-1+vt*sin(yawt)*Δt;
上式中vt表示在t时刻车辆的速度,yawt表示车辆在t时刻的航向角,xt、yt分别表示在t时刻车辆在UTM坐标系下的横坐标和纵坐标,vt-1表示t-1时刻车辆的速度,yawt-1表示t-1时刻车辆的航向角,xt-1和yt-1分别代表在t-1时刻车辆在UTM坐标系下的横坐标和纵坐标,Δt表示时间差,表示t-1时刻车辆的速度比例因子,表示t-1时刻车辆的航向角偏差,为t-1时刻航向角的角速度。
步骤302:把车辆位置信息采集单元采集的车辆的经纬度数据转转化到UTM平面坐标,卡尔曼滤波器中的更新模型分别根据GPS、IMU、ODOM输出的位置信息、航向、速度值,更新***状态量矩阵X(x,y,v,yaw,yawv,vscale,yawbias);
其中,vscale的调节更新方法:根据预测模型和根据GPS输出的位置(zx,zy)和车辆的速度zv可以计算出vscale的值,在预测模型计算中***车辆的速度的计算为其中含有上一时刻的速度vt-1和比例因子vscale,同时将速度vt带入到公式xt=xt-1+vt*cos(yawt)*Δt和yt=yt-1+vt*sin(yawt)*Δt中;在更新模型中,传感器会输入GPS的位置信息(x,y)和速度v,此时会根据预测值和传感器的输出值和其噪声协方差矩阵R和协方差矩阵P更新***变量X中的vscale,用以修正车辆速度。
yawbias的调节更新方法:根据预测模型和根据GPS输出的位置(zx,zy)和IMU输出的航向角yaw和yawv,可以计算出yawbias的值。在预测模型计算中***车辆的航向角的计算公式为:其中含有yawbias;将航向角yawt带入到公式xt=xt-1+vt*cos(yawt)*Δt和yt=yt-1+vt*sin(yawt)*Δt中;在更新模块部分,传感器会输入GPS的位置信息(x,y)和imu输出的航向角以及角速度,此时会根据预测值和传感器的输出值和其噪声协方差矩阵R和协方差矩阵P更新***变量X中的yawbias,用以修正车辆的航向角。
步骤303:不断重复步骤301~步骤302,每0.02秒输出一次车辆的定位信息,输出的定位信息即为当前***状态量矩阵X。
其中,当GPS定位精度高、车速大的情况的下,用其航向信息修正IMU输入到卡尔曼滤波器中的航向信息。其公式如下:
yawoffset=yawgps-yawimu;
yawinput=yawimu+yawoffest;
上式中yawgps为GPS星数大于等于12,车速大于5m/s、水平精度因子小于0.6时,RMC报文输出的航向角度,yawimu为imu输出的航向角,yawoffset为偏移角,yawinput为输入卡尔曼滤波器的yaw值,即为Z_yaw的值。
步骤2中初始化或者修改卡尔曼滤波器参数的方法为:根据GPS在动态和静态情况下定位精度不同、GPS的星数star_num和水平精度因子hdop等参数以及信号特性设置一套状态机调参方法,动态设置***的观测值Z和测量噪声协方差R。主要包括重启、GPS无效、GPS有效、停车/缓慢行驶四个状态,并据此设置卡尔曼滤波中的测量噪声协方差R,动态调整算法参数。各状态分别如下:
重启状态:初始化卡尔曼滤波器参数中的***协方差矩阵P、测量噪声协方差矩阵R和***状态量矩阵X。初始化时,***状态量矩阵X中的x、y设置为GPS经纬度转换成UTM系下的横、纵坐标,其它值设为0;***协方差矩阵P对角线上元素值设为0.01,测量噪声协方差矩阵R对角线上元素值设为0.1。
GPS无效状态:当车辆长期经过桥底、隧道等信号遮挡地方时,GPS输出的水平精度因子较高、接收星数较少,此时GPS无效,GPS采集的位置信息不输入卡尔曼滤波器,用***变量X中的值和传感器输出的速度和航向值输入预测模型进行推算。
GPS有效状态:当车辆在空旷路况或城市道路正常行驶,GPS定位数据可用。但当经过高楼密集或者数目遮挡地方时,GPS的定位数据不稳定,会出现抖动情况。此时根据水平精度因子hdop动态调整GPS观测方差,hdop越小,测量噪声方差越小。即Hdop<0.55时,rx_x和ry_y设置为1.0;0.55<Hdop<0.65时,rx_x和ry_y设置为1.5;0.65<Hdop<0.68时,rx_x和ry_y设置为2.5;0.68<Hdop<0.73时,rx_x和ry_y设置为15;Hdop>0.75时,rx_x和ry_y设置为150。
停车/缓慢行驶状态:当从车辆速度较小时,GPS的定位精度较低。将GPS的位置信息每间隔1秒输入卡尔曼滤波器一次,同时设置GPS的测量噪声方差,将GPS的测量噪声协方差rx_x和ry_y均设置成2.0。其中,rx_x和ry_y分别表示GPS在横向和纵向位置的方差。
如图2所示,当***上电后,进入重启状态,在重启状态时,当hdop>h_mid或者star_num<10时候,进入GPS无效状态;当hdop<=h_mid且v_car>v_mid时,进入GPS有效状态;当hdop<=h_mid且v_car<=v_mid时,进入停车/缓慢行驶状态。
在GPS无效状态时,若hdop<h_good,说明GPS定位精度提高,此时进入GPS有效状态;若P异常,进入重启状态;
在GPS有效状态时,若hdop>h_bad,说明车辆进入隧道或者高架桥下,此时GPS信息不可以,便进入GPS无效状态;若v_car<v_ba,说明车速降低,则进入停车/缓慢行驶状态;若P异常,进入重启状态;
在停车/缓慢行驶状态时,若v_car>v_good,说明车速提高,则GPS定位精度增加,便进入GPS有效状态;若hdop>h_bad,说明GPS天线受遮挡,便进入GPS无效状态。若P异常,进入重启状态。
其中hdop为水平精度因子、star_num为接收星数、v_car为从ODOM读取的当前车速、P为卡尔曼滤波器的***协方差,h_good、h_mid、h_bad为算法中判断GPS定位精度参数,h_good的优选范围为0.55-0.65,h_mid的优选范围为0.68-0.73,h_bad的优选范围为0.75-0.77;v_good、v_bad、v_mid为算法中判断车辆速度参数,v_good的优选范围为3-5m/s,v_mid的优选范围为1-2m/s,v_bad的优选范围为0.3_0.5m/s。
其中,当***初始化或者运行过程中***的协方差P值出现非正定或其中横纵向和其它参数的协方差值存在明显非对称差异表示P异常。例如,***协方差矩阵P中的元素P_x_v和p_y_v差值大于1000则认为是横纵向和其它参数的协方差值存在明显非对称差异,属于P异常的情况。
如图3所示,在城市道路,穿越高架桥时的定位轨迹对比图,其中,黑色轨迹为本发明实现的定位结果,灰色轨迹为DJI-A3的定位结果。可明显看出黑色轨迹相对平滑,灰色轨迹在桥底出现抖动情况。
采用本发明提供的定位方法与现有技术中的DJI-A3和TrimbleDB982定位效果,以及整个***的成本等相关内容进行比较,如表一所示,可以看出本发明提供的方法在降低了整个***硬件成本的同时,保证了***的定位精度,不管在有遮挡还是没有遮挡的地方,都可以准确的定位到车辆信息。
表1:
Claims (5)
1.一种用于自动驾驶的低成本传感器组合定位***,其特征在于:包括车辆速度信息采集单元、车辆航向信息采集单元、车辆位置信息采集单元和车辆位置定位单元,其中,车辆速度信息采集单元用于采集车辆的实时车速信息;车辆航向信息采集单元用于实时采集车辆的航向信息;车辆位置信息采集单元为GPS,用于实时采集车辆位置信息;所述车辆位置定位单元包括状态机和卡尔曼滤波器,所述状态机根据GPS反馈的参数对卡尔曼滤波器的参数进行调整,所述卡尔曼滤波器根据车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元输出的车辆的信息,结合卡尔曼滤波器的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出对车辆进行定位;所述预测模型为:
xt=xt-1+vt*cos(yawt)*Δt;
yt=yt-1+vt*sin(yawt)*Δt;
上式中vt表示在t时刻车辆的速度,yawt表示车辆在t时刻的航向角,xt、yt分别表示在t时刻车辆在UTM坐标系下的横坐标和纵坐标,vt-1表示t-1时刻车辆的速度,yawt-1表示t-1时刻车辆的航向角,xt-1和yt-1分别代表在t-1时刻车辆在UTM坐标系下的横坐标和纵坐标,Δt表示时间差,表示t-1时刻车辆的速度比例因子,表示t-1时刻车辆的航向角偏差,为t-1时刻航向角的角速度;
所述状态机中包括重启状态、GPS无效状态、GPS有效状态和停车/缓慢行驶状态,其中,重启状态时,初始化卡尔曼滤波器参数中的***协方差矩阵P、测量噪声协方差矩阵R和***状态量矩阵X;GPS无效状态时,GPS采集的位置信息不输入卡尔曼滤波器;GPS有效状态时,GPS采集的位置信息输入卡尔曼滤波器;停车/缓慢行驶状态时,将GPS的位置信息每间隔1秒输入卡尔曼滤波器一次;当***上电后,进入重启状态,在重启状态时,当hdop>h_mid或者star_num<10时候,进入GPS无效状态;当hdop<=h_mid且v_car>v_mid时,进入GPS有效状态;当hdop<=h_mid且v_car<=v_mid时,进入停车/缓慢行驶状态;在GPS无效状态时,若hdop<h_good,进入GPS有效状态;若P异常,进入重启状态;在GPS有效状态时,若hdop>h_bad,进入GPS无效状态;若v_car<v_bad,则进入停车/缓慢行驶状态;若P异常,进入重启状态;在停车/缓慢行驶状态时,若v_car>v_good,进入GPS有效状态;若hdop>h_bad,进入GPS无效状态;若P异常,进入重启状态;hdop为水平精度因子、star_num为接收星数、v_car为从ODOM读取的当前车速、P为卡尔曼滤波器的***协方差矩阵,h_good、h_mid、h_bad为判断GPS定位精度参数,v_good、v_bad、v_mid为判断车辆速度参数;车辆速度信息采集单元采用设置在车辆上的ODOM,ODOM采集车辆的前进速度信息和旋转速度;车辆航向信息采集单元采用装在车辆旋转中心位置的IMU,其主要采集和提供车辆的航向信息;车辆位置信息采集单元采用设置在车顶并且位于车辆旋转中心位置的GPS,其主要采集和提供车辆的位置信息。
2.根据权利要求1所述的用于自动驾驶的低成本传感器组合定位***,其特征在于:h_good的范围为0.55-0.65,h_mid的范围为0.68-0.73,h_bad的范围为0.75-0.77。
3.根据权利要求1所述的用于自动驾驶的低成本传感器组合定位***,其特征在于:v_good的范围为3-5m/s,v_mid的范围为1-2m/s,v_bad的范围为0.3-0.5m/s。
4.一种基于权利要求1所述的用于自动驾驶的低成本传感器组合定位***的定位方法,其特征在于:包括以下步骤:
步骤1:车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元分别将采集到车辆的速度、车辆的航向信息和车辆的位置信息发送到车辆位置定位单元;
步骤2:在车辆位置定位单元构建卡尔曼滤波器,车辆位置定位单元中的状态机根据获得的车辆位置信息,初始化或者修改卡尔曼滤波器参数;
步骤3:卡尔曼滤波器根据接收到的车辆的速度、车辆的航向信息和车辆的位置信息,结合其中的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出;
所述步骤2中初始化或者修改卡尔曼滤波器参数的方法为:设定重启状态、GPS无效状态、GPS有效状态和停车/缓慢行驶状态,其中,重启状态时,初始化卡尔曼滤波器参数中的***协方差矩阵P、测量噪声协方差矩阵R和***状态量矩阵X;GPS无效状态时,GPS采集的位置信息不输入卡尔曼滤波器;GPS有效状态时,GPS采集的位置信息输入卡尔曼滤波器;停车/缓慢行驶状态时,将GPS的位置信息每间隔1秒输入卡尔曼滤波器一次;当***上电后,进入重启状态,在重启状态时,当hdop>h_mid或者star_num<10时候,进入GPS无效状态;当hdop<=h_mid且v_car>v_mid时,进入GPS 有效状态;当hdop<=h_mid且v_car<=v_mid时,进入停车/缓慢行驶状态;在GPS无效状态时,若hdop<h_good,进入GPS有效状态;若P异常,进入重启状态;在GPS有效状态时,若hdop>h_bad,进入GPS无效状态;若v_car<v_bad,则进入停车/缓慢行驶状态;若P异常,进入重启状态;在停车/缓慢行驶状态时,若v_car>v_good,进入GPS有效状态;若hdop>h_bad,进入GPS无效状态;若P异常,进入重启状态;hdop为水平精度因子、star_num为接收星数、v_car为从ODOM读取的当前车速、P为卡尔曼滤波器的***协方差矩阵,h_good、h_mid、h_bad为判断GPS定位精度参数,v_good、v_bad、v_mid为判断车辆速度参数;
所述预测模型为:
xt=xt-1+vt*cos(yawt)*Δt;
yt=yt-1+vt*sin(yawt)*Δt;
5.根据权利要求4所述的定位方法,其特征在于:所述步骤2中还包括用GPS采集的航向信息修正车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息;根据公式:
yawoffset=yawgps-yawimu;
yawinput=yawimu+yawoffest;
对车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息进行修正,式中yawgps为GPS星数大于等于12,车速大于5m/s、水平精度因子小于0.6时,RMC报文输出的航向角度,yawimu为IMU 输出的航向角,yawoffset为偏移角,yawinput为输入卡尔曼滤波器的yaw值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911168614.XA CN110926483B (zh) | 2019-11-25 | 2019-11-25 | 一种用于自动驾驶的低成本传感器组合定位***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911168614.XA CN110926483B (zh) | 2019-11-25 | 2019-11-25 | 一种用于自动驾驶的低成本传感器组合定位***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110926483A CN110926483A (zh) | 2020-03-27 |
CN110926483B true CN110926483B (zh) | 2020-12-25 |
Family
ID=69851986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911168614.XA Active CN110926483B (zh) | 2019-11-25 | 2019-11-25 | 一种用于自动驾驶的低成本传感器组合定位***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110926483B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112505718B (zh) * | 2020-11-10 | 2022-03-01 | 奥特酷智能科技(南京)有限公司 | 用于自动驾驶车辆的定位方法、***及计算机可读介质 |
CN115932924B (zh) * | 2022-07-29 | 2023-09-05 | 润芯微科技(江苏)有限公司 | 一种基于imu辅助定位的方法及*** |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2898196B1 (fr) * | 2006-03-01 | 2008-04-25 | Eurocopter France | Procede et dispositif de positionnement hybride |
US8473207B2 (en) * | 2008-10-21 | 2013-06-25 | Texas Instruments Incorporated | Tightly-coupled GNSS/IMU integration filter having calibration features |
CN102608641A (zh) * | 2012-03-30 | 2012-07-25 | 江苏物联网研究发展中心 | 单轴陀螺仪和单轴加速度计的车载组合导航***及方法 |
CN103472459A (zh) * | 2013-08-29 | 2013-12-25 | 镇江青思网络科技有限公司 | 一种基于gps伪距差分的车辆协作定位方法 |
CN105606094B (zh) * | 2016-02-19 | 2018-08-21 | 北京航天控制仪器研究所 | 一种基于mems/gps组合***的信息条件匹配滤波估计方法 |
CN107436444A (zh) * | 2017-06-23 | 2017-12-05 | 北京机械设备研究所 | 一种车载多模式组合导航***及方法 |
CN108983270A (zh) * | 2018-06-14 | 2018-12-11 | 兰州晨阳启创信息科技有限公司 | 一种基于多传感器融合的列车安全定位***及方法 |
CN108827292A (zh) * | 2018-06-27 | 2018-11-16 | 四川大学 | 基于gnss和地面基站的组合式导航精确测速定位方法及*** |
CN109343095B (zh) * | 2018-11-15 | 2020-09-01 | 众泰新能源汽车有限公司 | 一种车载导航车辆组合定位装置及其组合定位方法 |
-
2019
- 2019-11-25 CN CN201911168614.XA patent/CN110926483B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110926483A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109946732B (zh) | 一种基于多传感器数据融合的无人车定位方法 | |
CN110631593B (zh) | 一种用于自动驾驶场景的多传感器融合定位方法 | |
CN110940344B (zh) | 一种用于自动驾驶的低成本传感器组合定位方法 | |
Ahmed et al. | Accurate attitude estimation of a moving land vehicle using low-cost MEMS IMU sensors | |
CN110307836B (zh) | 一种用于无人清扫车辆贴边清扫的精确定位方法 | |
CN110160542B (zh) | 车道线的定位方法和装置、存储介质、电子装置 | |
EP3109589B1 (en) | A unit and method for improving positioning accuracy | |
Georgy et al. | Modeling the stochastic drift of a MEMS-based gyroscope in gyro/odometer/GPS integrated navigation | |
CN104061899B (zh) | 一种基于卡尔曼滤波的车辆侧倾角与俯仰角估计方法 | |
CN113945206A (zh) | 一种基于多传感器融合的定位方法及装置 | |
Fakharian et al. | Adaptive Kalman filtering based navigation: An IMU/GPS integration approach | |
CN109186597B (zh) | 一种基于双mems-imu的室内轮式机器人的定位方法 | |
CN108931244B (zh) | 基于列车运动约束的惯导误差抑制方法及*** | |
CN112505737B (zh) | 一种gnss/ins组合导航方法 | |
CN109343095B (zh) | 一种车载导航车辆组合定位装置及其组合定位方法 | |
CN107132563B (zh) | 一种里程计结合双天线差分gnss的组合导航方法 | |
CN110530361B (zh) | 基于农业机械双天线gnss自动导航***的转向角度估计器 | |
Chen et al. | A novel fusion methodology to bridge GPS outages for land vehicle positioning | |
CN110926483B (zh) | 一种用于自动驾驶的低成本传感器组合定位***及方法 | |
US20230182790A1 (en) | Method for calculating an instantaneous velocity vector of a rail vehicle and corresponding system | |
Gao et al. | An integrated land vehicle navigation system based on context awareness | |
CN112433531A (zh) | 一种自动驾驶车辆的轨迹跟踪方法、装置及计算机设备 | |
CN111678514A (zh) | 一种基于载体运动条件约束和单轴旋转调制的车载自主导航方法 | |
CN106646569B (zh) | 一种导航定位方法及设备 | |
CN115451949A (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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 210012 room 401-404, building 5, chuqiaocheng, No. 57, Andemen street, Yuhuatai District, Nanjing, Jiangsu Province Patentee after: AUTOCORE INTELLIGENT TECHNOLOGY (NANJING) Co.,Ltd. Address before: 211800 building 12-289, 29 buyue Road, Qiaolin street, Jiangbei new district, Pukou District, Nanjing City, Jiangsu Province Patentee before: AUTOCORE INTELLIGENT TECHNOLOGY (NANJING) Co.,Ltd. |