CN110926483B - 一种用于自动驾驶的低成本传感器组合定位***及方法 - Google Patents

一种用于自动驾驶的低成本传感器组合定位***及方法 Download PDF

Info

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
Application number
CN201911168614.XA
Other languages
English (en)
Other versions
CN110926483A (zh
Inventor
李赵
袁瑞廷
刘心刚
陈诚
张旸
刘洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AutoCore Intelligence Technology Nanjing Co Ltd
Original Assignee
AutoCore Intelligence Technology Nanjing Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AutoCore Intelligence Technology Nanjing Co Ltd filed Critical AutoCore Intelligence Technology Nanjing Co Ltd
Priority to CN201911168614.XA priority Critical patent/CN110926483B/zh
Publication of CN110926483A publication Critical patent/CN110926483A/zh
Application granted granted Critical
Publication of CN110926483B publication Critical patent/CN110926483B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/28Navigation; 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining 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)和观测值矩阵
Figure BDA0002288120700000021
其中,x表示输入卡尔曼滤波器的车辆在UTM坐标系下的横坐标,y表示输入卡尔曼滤波器的车辆在UTM坐标系下的纵坐标,v表示输入卡尔曼滤波器的车辆的速度,yaw表示输入卡尔曼滤波器的车辆的航向角,yawv表示输入卡尔曼滤波器的车辆的航向角速度,vscale表示输入卡尔曼滤波器的速度比例因子,yawbias表示输入卡尔曼滤波器的航向角偏差;zx和zy表示GPS输出的位置信息,zyaw是IMU输出的航向角,zv为ODMO输出的速度信息,zyawv是IMU输出的航向角的角速度。速度比例因子vscale解决由车辆轮胎气压、车辆载重、路面不平整等造成的测量速度和真实速度不一致的情况,引入航向角偏差yawbias解决道路磁场环境变换引起的航向角不稳定的情况,这两个状态量能够根据GPS的观测数据和卡尔曼滤波算法的预测模型动态调整。更加有效的提高对车辆定位的精度。
进一步,所述预测模型为:
Figure BDA0002288120700000031
Figure BDA0002288120700000032
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表示时间差,
Figure BDA0002288120700000033
表示t-1时刻车辆的速度比例因子,
Figure BDA0002288120700000034
表示t-1时刻车辆的航向角偏差,
Figure BDA0002288120700000035
为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:在车辆位置定位单元构建卡尔曼滤波器,车辆位置定位单元根据获得的车辆位置信息,初始化或者修改卡尔曼滤波器参数,其中卡尔曼滤波器参数包括***协方差矩阵
Figure BDA0002288120700000061
测量噪声协方差矩阵
Figure BDA0002288120700000062
***状态量矩阵X(x,y,v,yaw,yawv,vscale,yawbias)和观测值矩阵Z(zx,zy,zyaw
Figure BDA0002288120700000063
zv),其中,Px_x表示车辆在UTM坐标系下的横坐标的方差值,px_y表示车辆在UTM坐标系下的横坐标与纵坐标的协方差值,px_v表示车辆在UTM坐标系下的横坐标与车辆速度的协方差值,Px_yaw表示车辆在UTM坐标系下的横坐标与航向角的协方差值,
Figure BDA0002288120700000064
表示车辆在UTM坐标系下的横坐标与航向角速度的协方差值,
Figure BDA0002288120700000065
表示车辆在UTM坐标系下的横坐标与速度比例因子的协方差值,
Figure BDA0002288120700000066
表示车辆在UTM坐标系下的横坐标与航向角偏差的协方差值,Py_x表示车辆在UTM坐标系下的纵坐标与横坐标的协方差值,py_y表示车辆在UTM坐标系下的纵坐标的方差值,py_v表示车辆在UTM坐标系下的纵坐标与车辆速度的协方差值,py_yaw表示车辆在UTM坐标系下的纵坐标与航向角的协方差值,
Figure BDA0002288120700000067
表示车辆在UTM坐标系下的纵坐标与航向角速度的协方差值,
Figure BDA0002288120700000068
表示车辆在UTM坐标系下的纵坐标与速度比例因子的协方差值,
Figure BDA0002288120700000069
表示车辆在UTM坐标系下的纵坐标与航向角偏差的协方差值,Pv_x表示车辆的速度与在UTM坐标系下横坐标的协方差值,pv_y表示车辆的速度与在UTM坐标系下纵坐标的协方差值,pv_v表示车辆速度的方差值,pv_yaw表示车辆速度与航向角的协方差值,
Figure BDA00022881207000000610
表示车辆速度与航向角速度的协方差值,
Figure BDA00022881207000000611
表示车辆速度与速度比例因子的协方差值,
Figure BDA0002288120700000071
表示车辆速度与航向角偏差的协方差值,Pyaw_x表示车辆的航向角与车辆在UTM坐标系下的横坐标的协方差值,pyaw_y表示车辆的航向角与车辆在UTM坐标系下的纵坐标的协方差值,pyaw_v表示车辆的航向角与车辆速度的协方差值,pyaw_yaw表示车辆航向角的方差值,
Figure BDA0002288120700000072
表示车辆航向角与航向角速度的协方差值,pyaw_vsacle表示车辆航向角与速度比例因子的协方差值,
Figure BDA0002288120700000073
表示车辆航向角与航向角偏差的协方差值,
Figure BDA0002288120700000074
表示车辆航向角速度与车辆在UTM坐标系下的横坐标的协方差值,
Figure BDA0002288120700000075
表示车辆航向角速度与车辆在UTM坐标系下的纵坐标的协方差值,
Figure BDA0002288120700000076
表示车辆航向角速度与车辆速度的协方差值,
Figure BDA0002288120700000077
表示车辆航向角速度与车辆航向角的协方差值,
Figure BDA0002288120700000078
表示车辆航向角速度的方差值,
Figure BDA0002288120700000079
表示车辆航向角速度与速度比例因子的协方差值,
Figure BDA00022881207000000710
表示车辆航向角速度与航向角偏差的协方差值,
Figure BDA00022881207000000711
表示车辆的速度比例因子与车辆在UTM坐标系下的横坐标的协方差值,
Figure BDA00022881207000000712
表示车辆的速度比例因子与车辆在UTM坐标系下的纵坐标的协方差值,
Figure BDA00022881207000000713
表示车辆的速度比例因子与车辆速度的协方差值,
Figure BDA00022881207000000714
表示车辆的速度比例因子与车辆航向角的协方差值,
Figure BDA00022881207000000715
表示车辆的速度比例因子与车辆航向角速度的协方差值,
Figure BDA00022881207000000716
表示车辆的速度比例因子的方差值,
Figure BDA00022881207000000717
表示车辆的速度比例因子与车辆航向角偏差的协方差值,
Figure BDA00022881207000000718
表示车辆航向角偏差与车辆在UTM坐标系下的横坐标的协方差值,
Figure BDA00022881207000000719
表示车辆航向角偏差与车辆在UTM坐标系下的纵坐标的协方差值,
Figure BDA00022881207000000720
表示车辆航向角偏差与车辆速度的协方差值,
Figure BDA00022881207000000721
表示车辆航向角偏差与车辆航向角的协方差值,
Figure BDA00022881207000000722
表示车辆航向角偏差与车辆航向角速度的协方差值,
Figure BDA00022881207000000723
表示车辆航向角偏差与车辆的速度比例因子的协方差值,
Figure BDA00022881207000000724
表示车辆航向角偏差的方差值;rx_x表示GPS在横向位置的方差,ry_y表示GPS在纵向位置的方差,rv_v表示车辆速度的方差,ryaw_yaw表示车辆航向角的方差,
Figure BDA00022881207000000725
表示车辆航向角速度的方差,
Figure BDA00022881207000000726
表示车辆的速度比例因子的方差,
Figure BDA00022881207000000727
表示车辆航向角偏差的方差,其中,矩阵R中的元素是根据传感器的特性设定的方差,其与传感器的精度相关,
Figure BDA00022881207000000728
Figure BDA00022881207000000729
的值为0;x表示输入卡尔曼滤波器的车辆在UTM坐标系下的横坐标,y表示输入卡尔曼滤波器的车辆在UTM坐标系下的纵坐标,v表示输入卡尔曼滤波器的车辆的速度,yaw表示输入卡尔曼滤波器的车辆的航向角,yawv表示输入卡尔曼滤波器的车辆的航向角速度,vscale表示输入卡尔曼滤波器的速度比例因子,yawbias表示输入卡尔曼滤波器的航向角偏差;zx和zy分别表示GPS输出的位置信息在UTM坐标系下的横坐标和纵坐标,zyaw是IMU输出的航向角,zv为ODMO输出的速度信息,
Figure BDA0002288120700000081
是IMU输出的航向角的角速度。
步骤3:卡尔曼滤波器根据接收到的车辆的车辆的速度、车辆的航向信息和车辆的位置信息,结合其中的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出。具体包括如下步骤:
步骤301:根据卡尔曼滤波器中的CTRV预测模型分别预测下一时刻的***状态量X中每个量的值;其中,CTRV预测模型为:
Figure BDA0002288120700000082
Figure BDA0002288120700000083
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表示时间差,
Figure BDA0002288120700000084
表示t-1时刻车辆的速度比例因子,
Figure BDA0002288120700000085
表示t-1时刻车辆的航向角偏差,
Figure BDA0002288120700000086
为t-1时刻航向角的角速度。
步骤302:把车辆位置信息采集单元采集的车辆的经纬度数据转转化到UTM平面坐标,卡尔曼滤波器中的更新模型分别根据GPS、IMU、ODOM输出的位置信息、航向、速度值,更新***状态量矩阵X(x,y,v,yaw,yawv,vscale,yawbias);
其中,vscale的调节更新方法:根据预测模型和根据GPS输出的位置(zx,zy)和车辆的速度zv可以计算出vscale的值,在预测模型计算中***车辆的速度的计算为
Figure BDA0002288120700000087
其中含有上一时刻的速度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的值。在预测模型计算中***车辆的航向角的计算公式为:
Figure BDA0002288120700000091
其中含有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:
Figure BDA0002288120700000111

Claims (5)

1.一种用于自动驾驶的低成本传感器组合定位***,其特征在于:包括车辆速度信息采集单元、车辆航向信息采集单元、车辆位置信息采集单元和车辆位置定位单元,其中,车辆速度信息采集单元用于采集车辆的实时车速信息;车辆航向信息采集单元用于实时采集车辆的航向信息;车辆位置信息采集单元为GPS,用于实时采集车辆位置信息;所述车辆位置定位单元包括状态机和卡尔曼滤波器,所述状态机根据GPS反馈的参数对卡尔曼滤波器的参数进行调整,所述卡尔曼滤波器根据车辆速度信息采集单元、车辆航向信息采集单元和车辆位置信息采集单元输出的车辆的信息,结合卡尔曼滤波器的预测模型和更新模型不断的迭代得到最终的车辆定位信息并输出对车辆进行定位;所述预测模型为:
Figure FDA0002770180700000011
Figure FDA0002770180700000012
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表示时间差,
Figure FDA0002770180700000013
表示t-1时刻车辆的速度比例因子,
Figure FDA0002770180700000014
表示t-1时刻车辆的航向角偏差,
Figure FDA0002770180700000015
为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为判断车辆速度参数;
所述预测模型为:
Figure FDA0002770180700000031
Figure FDA0002770180700000032
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表示时间差,
Figure FDA0002770180700000033
表示t-1时刻车辆的速度比例因子,
Figure FDA0002770180700000034
表示t-1时刻车辆的航向角偏差,
Figure FDA0002770180700000035
为t-1时刻航向角的角速度。
5.根据权利要求4所述的定位方法,其特征在于:所述步骤2中还包括用GPS采集的航向信息修正车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息;根据公式:
yawoffset=yawgps-yawimu
yawinput=yawimu+yawoffest
对车辆航向信息采集单元输入到卡尔曼滤波器中的航向信息进行修正,式中yawgps为GPS星数大于等于12,车速大于5m/s、水平精度因子小于0.6时,RMC报文输出的航向角度,yawimu为IMU 输出的航向角,yawoffset为偏移角,yawinput为输入卡尔曼滤波器的yaw值。
CN201911168614.XA 2019-11-25 2019-11-25 一种用于自动驾驶的低成本传感器组合定位***及方法 Active CN110926483B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 众泰新能源汽车有限公司 一种车载导航车辆组合定位装置及其组合定位方法

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.