CN112362052B - 一种融合定位方法及系统 - Google Patents
一种融合定位方法及系统 Download PDFInfo
- Publication number
- CN112362052B CN112362052B CN202011163752.1A CN202011163752A CN112362052B CN 112362052 B CN112362052 B CN 112362052B CN 202011163752 A CN202011163752 A CN 202011163752A CN 112362052 B CN112362052 B CN 112362052B
- Authority
- CN
- China
- Prior art keywords
- rfid
- inertial navigation
- module
- noise
- navigation module
- 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 63
- 230000004927 fusion Effects 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 48
- 238000001914 filtration Methods 0.000 claims abstract description 40
- 238000005259 measurement Methods 0.000 claims abstract description 20
- 230000001133 acceleration Effects 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 21
- 230000004807 localization Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000007704 transition Effects 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 230000003068 static effect Effects 0.000 abstract description 8
- 238000012937 correction Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- 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
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
本发明提供一种融合定位方法及系统。本发明的方法用于对具有行进轮的目标主体进行定位,所述目标主体具有惯性导航模块和基于行进轮的速度测量模块,所述方法包括:在所述目标主体保持静止,行进轮旋转的情况下,测定惯性导航模块和速度测量模块的噪声方差;构建与所述目标主体的位置相关的卡尔曼滤波状态量和观测量;基于所述噪声方差以及所述状态量和观测量构建表征位置预测不确定性的噪声矩阵;实时测量所述观测量,并基于所测得的观测量基于卡尔曼滤波进行所述目标主体的位置预测。
Description
技术领域
本发明涉及导航定位技术领域,更具体涉及组合导航技术领域。本发明可应用于室外无人驾驶系统、室内智能小车系统等需要移动主体定位信息的系统。
背景技术
车辆导航定位是无人驾驶系统的关键技术之一,精密、安全的无人驾驶系统需在任何时间和空间获得高精度的车辆定位信息。由于基于微机电系统(Micro-Electro-Mechanical System,MEMS)技术的惯性导航模块(Inertial Navigation System,INS)的加速度计测量结果受多种噪声影响,基于惯性导航模块航位推算的车辆定位结果误差随时间迅速累积,无法长时间定位,需与其他传感器结合使用。目前,惯性导航模块与速度计(里程计)和射频识别技术(Radio Frequency Identification,RFID)融合是一种可能的定位方案。以如图1所示的定位场景为例,路面设置RFID标签,车辆上安装车载RFID阅读器,惯性导航模块和速度计安装在车辆上。RFID标签提供离散的绝对位置坐标,惯性导航模块与速度计完成相对位置推算。由于RFID标签能够对惯性导航模块与速度计进行修正,在整体上降低了惯性导航模块与速度计长距离定位难度,是一种可行的高精度定位方法。然而,针对无人驾驶系统场景,此方案相关算法缺乏这一场景下可实现的构造方法,同时还存在进一步提升定位精度的空间。
目前,现有大量方案均是基于传感器误差的卡尔曼滤波融合定位算法,将惯性导航模块、速度计和RFID传感器在各个方向的系统误差作为卡尔曼滤波的状态量,以惯性导航模块推算位置与RFID位置的差值、惯性导航模块推算速度与速度计的差值作为观测量进行卡尔曼滤波得到定位结果。比如,孙伟等人在中国专利CN107402005A中公开的“一种基于惯性/里程计/RFID的高精度组合导航方法”,以及KIM IN-JUN(KR)等在美国专利US2005143916中公开的“Positioning apparatus and method combining RFID,GPS andINS”都是基于上述方法。
但是,现有的导航定位方案主要存在三个问题:
(1)观测量使用的是惯性导航模块推算的位置和速度,即,观测量构建方式直接采用惯性导航模块单独推算的位置和速度与其他传感器对应量相减,而由于惯性导航模块误差累积速度非常快,导致观测量严重远离实际误差,难以达到高精度的定位;
(2)缺少对噪声协方差矩阵构建,而卡尔曼滤波的精度和收敛性取决于过程噪声和观测噪声协方差矩阵,因此如何设计十分重要,但是现有技术均根据经验和实际效果进行调参,未给出如何根据其物理意义构建噪声协方差矩阵,以避免盲目调参的问题;
(3)在RFID标签处,以RFID内存储的绝对位置信息直接去修正如惯性导航模块等其他传感器。而在实际中,RFID阅读器天线存在一定的覆盖范围,当车辆阅读到RFID标签时,车辆与RFID标签之间存在一定的位置偏差,直接将RFID标签存储的绝对位置信息视为车辆当前位置,误差最大可达RFID标签与阅读器的最大通信距离,大大降低位置修正精度。
因此,现有的基于传感器误差的卡尔曼滤波融合定位方法都存在上述问题中的至少一种,在定位精度上均有更进一步提高的需要。
发明内容
申请人首先从理论上简述本发明的融合定位方法的原理,卡尔曼滤波状态量和观测量的形式,并概括说明本发明的方案如何解决了背景技术中提到的三个问题。
针对上述问题(1):
现有技术中惯性导航模块可测量车辆三轴加速度,在每一个轴方向上,对加速度在时间上一次积分得到速度,对加速度在时间上二次积分得到位置,其位置推算方程为
其中Δt为惯导数据更新周期,pk为k时刻的位置,其由k-1时刻位置pk-1、速度vk-1和加速度ak-1计算得到。由于加速度计噪声的存在,由其推算出的速度和位置误差随时间累积。引入速度计,其可测量车辆线速度,结合惯性导航模块测量的偏航角分解可得各方向速度,其可对惯性导航模块一次积分结果,即速度,进行卡尔曼滤波融合;引入RFID标签,其内部存储绝对位置坐标信息,同时结合接收的信号强度指示RSSI(Received SignalStrength Indicator)解算出车辆相对于标签的位置,对惯性导航模块二次积分结果,即位置,进行修正。
本发明根据速度计和RFID标签的数据输出频率,设计了一种速度计与惯导卡尔曼滤波、离散RFID修正相结合的方法。卡尔曼滤波状态量为各轴方向的位置、速度和加速度;观测量为各轴方向离散的RFID处位置、速度计线速速度和惯导加速度;构造的状态量和观测量直接使用惯导的加速度原始值,不直接取积分结果,且结合精度较高的惯性导航模块偏航角对线速度分解,解决了问题(1)。
针对问题(2):
卡尔曼滤波中存在两个先验噪声量,过程噪声和观测噪声。过程噪声指从k-1时刻推测k时刻状态量的误差;观测噪声指传感器(惯性导航模块、速度计)观测数据的噪声。本发明方法中以协方差矩阵给出这两个噪声量。对于观测噪声,可根据实测的惯性导航模块和速度计的静态数据进行统计,得到噪声的方差情况,而不同方向、不同传感器噪声是独立的,因此可构建观测噪声协方差矩阵;对于过程噪声,除了自协方差,还存在位置、速度与加速度之间的耦合,其代表了加速度的不确定性带给速度预测和位置预测的不确定性度量,本发明中根据加速度计的噪声对其进行了推导,获得接近实际的噪声矩阵,使卡尔曼滤波具有高输出精度和快收敛性,解决了问题(2)。针对问题(3):
现有技术中,基于接收的信号强度指示的三边RFID标签定位在室内静态场景的方案有很多。但是发明人注意到在移动主体定位,比如,车辆定位中,RFID标签阅读器相对于RFID标签往往是单向运动的,基于这样的特征,本发明在基于RFID标签进行定位时,采用了一种创新的RFID标签定位方法,即提出了一种基于速度判断的双阅读器定位方法,利用双阅读器接收的信号强度指示得到两个可行解,再结合速度确定唯一解。最终在RFID标签处得到精度更高的车辆位置修正,解决了问题(3)。
针对上述问题,本申请的发明人提出了一种融合定位方法及系统,其可以实现基于惯导/RFID/速度计融合的高精度定位。
根据本发明的一个方面,提供了一种融合定位方法,所述方法用于对具有行进轮的目标主体进行定位,所述目标主体具有惯性导航模块和基于行进轮的速度测量模块,其特征在于,所述方法包括:
步骤1:利用所述惯性导航模块和速度测量模块实时测量所述目标主体的加速度和速度;
步骤2:将所测得的速度和加速度作为滤波观测量代入卡尔曼滤波方程组,进行所述目标主体的位置预测,其中,在所述卡尔曼滤波方程组中利用各个滤波状态量之间的噪声协方差矩阵对所述滤波状态量进行修正,所述噪声协方差矩阵基于所述惯性导航模块和速度测量模块的噪声方差确定。
在一种优选实现方式中,在所述目标主体以及其行进路径上设置有基于位置标签的定位模块,所述方法还包括步骤3:利用所述定位模块获取所述目标主体的位置信息,基于所述位置信息对所述位置预测的结果进行修正。
在另一种优选实现方式中,所述惯性导航模块和速度测量模块的噪声方差在所述目标主体保持静止,行进轮旋转的情况下进行测定。
在另一种优选实现方式中,所述速度测量模块为速度计,所述定位模块包括设置于所述目标主体的RFID阅读器和设置于路面内的RFID标签,所述目标主体为车辆。
在另一种优选实现方式中,噪声协方差矩阵包括:
(1)过程噪声协方差矩阵Q,
以及(2)观测噪声协方差矩阵R,R=diag(Dv,Dv,Da,Da),其中,Dv为所述速度测量模块的噪声方差。
在另一种优选实现方式中,所述卡尔曼滤波方程组为:
在另一种优选实现方式中,所述滤波状态量为X=[px,py,vx,vy,ax,ay]T,其中,px,py为所述目标主体的二维位置坐标、vx,vy为所述目标主体的二维速度,ax,ay为所述目标主体的二维加速度;所述观测量为其中,vl,k为所述目标主体的线速度,θk为所述目标主体的偏航角。
在另一种优选实现方式中,所述目标主体内沿垂直于其行进方向横向布置两个RFID阅读器
所述步骤3包括:
3.1、基于两个所述RFID阅读器获取的RFID标签的信号强度指示值,由路径损耗模型解算出所述RFID标签与所述目标主体的两个可能相对位置;
3.2、基于两个所述RFID阅读器对所述RFID标签的识别时间、车速确定所述RFID标签的相对位置;
3.3、读取所述RFID标签内存储的所述RFID标签的绝对位置信息,
计算出所述目标主体的位置信息。
根据本发明的另一方面,提供了一种融合定位系统,所述系统用于对具有行进轮的目标主体进行定位,所述融合定位系统包括惯性导航模块、基于行进轮的速度测量模块以及融合定位解算模块,其特征在于,所述融合定位解算模块采用所述的方法进行定位解算。
根据本发明的另一方面,提供了一种计算机可读存储介质,其上存储有计算机程序,其中,该程序被处理器执行时实现根据所述的方法。
技术效果
本发明通过对卡尔曼滤波参数的重新设计,对卡尔曼滤波算法中状态量、观测量、过程噪声协方差矩阵和观测噪声协方差进行创新构造,可在任意时间和空间给出定位结果,实现高精度定位输出,满足实际应用要求。
此外,本发明在RFID部分仅使用两个RFID阅读器实现了比采用三个RFID阅读器同等、甚至更加准确的位置修正。本发明引入了基于速度判断的双阅读器RFID标签定位算法进一步提升了位置修正精度。
附图说明
以下附图仅对本发明作示意性的说明和解释,并不用于限定本发明的范围,其中:
图1为现有的惯性导航模块/RFID/速度计定位方案的场景示意图;
图2为本发明实施例中的惯性导航模块/RFID/速度计融合定位的流程示意图;
图3为本发明实施例中采用RFID标签进行定位的原理示意图;
图4为采用本发明实施例中的定位方案和采用现有方案的仿真性能对比图。
具体实施方式
为了使本发明的目的、技术方案、设计方法及优点更加清楚明了,以下结合附图通过具体实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
下面首先结合附图对本发明的原理以及实施例进行详细说明。
本实施例中,融合定位系统包括惯性导航模块、RFID标签及阅读器、速度计。惯性导航模块、RFID标签及阅读器、速度计的设置与图1中所示的现有设置方式类似,若干RFID标签设置于路面内,优选沿着路面中心线布置,惯性导航模块、RFID阅读器和速度计设置于目标车辆或其他移动主体内。
图2为本实施例中的惯性导航模块/RFID/速度计融合定位方法的示意流程图。当车辆在二维平面内运动时,首先实时利用车辆上所安装的惯性导航模块进行车辆运动中的三轴加速度、三个姿态角的测量运算,并将所解算出的三轴加速度、三个姿态角的信息提供给融合定位系统(比如,其中的融合定位模块或者其他运算单元),速度计测量车辆的线速度,RFID阅读器读取路面的RFID标签内存储的位置信息,进而提供绝对位置修正。
本实施例中,利用该融合定位系统进行定位的过程包括两部分,模型预先构建以及实时数据采集和位置预测。
其中,模型预先构建部分包括:
步骤1:噪声预估计;
步骤2:卡尔曼滤波状态量和观测量构建;
步骤3:噪声矩阵构建;
实时数据采集和位置预测部分包括:
步骤4:实时数据采集和基于卡尔曼滤波的位置预测;
步骤5:基于RFID标签定位对位置预测结果进行修正。
下面分别就上述五个步骤进行展开说明。
1、噪声预估计
噪声估计前需要数据采样,数据采样时,在车辆的车身静止、但是车轮运动的状态下测量加速度计和速度计的静态数据。可以将车辆后轮悬空,设定一定的转速,测量速度计和惯性导航模块中加速度计的数据,作为静态数据,低阶均值滤波后估计出加速度计噪声方差Da和速度计噪声方差Dv。
由于惯性导航模块和速度计测量值存在惯性噪声等诸多非高斯噪声成分,一般为在局部采样点上幅值较大的高频噪声,因此,在采样频率100Hz情况下,本实施例中首先对二者提供的原始数据使用阶数为10-20阶的低阶均值滤波来滤除这些噪声的影响,具体表达为ξ(k)为k时刻采样值,Ω为阶数。
2、卡尔曼滤波状态量和观测量构建
设卡尔曼滤波状态量为X=[px,py,vx,vy,ax,ay]T,其中px,py为导航坐标系下的二维位置坐标(即平面坐标,这里不考虑z方向),vx,vy为二维速度,ax,ay为二维加速度。则设计的状态预测方程为
式中Xk为k时刻滤波状态量,A代表状态转移矩阵,表示两个相邻时刻状态推算关系,Δt为惯性导航模块更新周期,Wk-1表示由第k-1时刻预测第k时刻滤波状态量的过程噪声值,W=[n1 n2 n3 n4 n5 n6]T,服从均值向量为06×1的6维联合正态分布,即
其中Q为6*6的协方差矩阵,n1 n2 n3 n4 n5 n6为各个状态量的过程噪声。当然,本领域技术人员应该理解,当在本实施例的状态量基础上额外增加状态量时,矩阵阶数也相应的增加,本实施例中的6行、6列仅仅是个示例。
本实施例的方案中,二维的加速度和速度为可观测量,其中加速度值由惯性导航模块测得,速度由速度计测得。观测方程为
其中Zk为观测状态量,为速度计和惯性导航模块测量到的XY方向速度和加速度,H为观测矩阵,代表了观测量和实际量的关系;Vk为k时刻观测噪声,其表征惯性导航模块和速度计传感器噪声情况,即为观测值与真值之间的差异,用协方差矩阵R衡量其统计特性。
由于速度计测量到的为车辆线速度vl,k,需根据惯性导航模块输出的偏航角θk分解至XY方向,即当前惯性导航模块偏航角精度高,误差小于1°,稳定性在0.1°内,可平滑输出角度值,因此可忽略其对速度噪声的影响。
3、噪声协方差矩阵构建
对于过程噪声协方差矩阵,由于状态预测方程中速度和位置均由加速度计算得到,因此协方差矩阵Q表示了预测过程中由加速度的不确定性带给速度预测和位置预测不确定性的度量。
本实施例中,根据协方差的定义,矩阵Q中从第一行到第六行以及第一列到第六列代表的变量与滤波状态量一致,即从px到ay的顺序,每个元素为所在行和列表示的相应滤波状态量之间协方差。例如,X方向位置px过程噪声n1包含了从k-1时刻推算k时刻时XY方向加速度、XY方向速度和XY方向位置六部分的不确定性,因此Q的第一行为
[Cov(px,px) Cov(px,py) Cov(px,vx) Cov(px,vy) Cov(px,ax) Cov(px,ay)]
同理第二行为
[Cov(py,px) Cov(py,py) Cov(py,vx) Cov(py,vy) Cov(py,ax) Cov(py,ay)]
以此类推,其中由位置推算的不确定性Cov(px,px)可由式(4)计算得到
其中Da代表加速度计噪声方差,其可由预先采集的数据估计。按照上述构造方法,可得到6*6的Q矩阵中每个元素值。特别地,由于两个时刻的加速度ax,k-1、ay,k-1与ax,k、ay,k不存在推算关系,故Cov(ax,ax)=0,Cov(ay,ay)=0;由于两个正交方向XY上的各参数(加速度、速度和位置)是独立的,因此不存在推算关系。最终协方差矩阵Q为
得到协方差矩阵后,可确定过程噪声值W的概率密度函数,即公式(3),由此可生成过程噪声值。
对于预测噪声协方差矩阵,由于速度计和加速度计是独立的,并且XY方向也是独立的,因此矩阵R为对角矩阵,对角元素按照Zk所代表的参数顺序自上到下分别为X方向速度计噪声方差、Y方向速度计噪声方差、X方向加速度计噪声方差和Y方向加速度计噪声方差,即Vk协方差矩阵R为
R=diag(Dv,Dv,Da,Da) (7)
4、实时数据采集和位置预测
正式测量时,将惯性导航模块和速度计所测得的数据使用卡尔曼滤波算法融合,引入过程噪声协方差矩阵和观测噪声协方差矩阵,该过程噪声协方差矩阵和观测噪声协方差矩阵已经由上述估计的噪声方差值推导得到。对惯性导航模块和速度计所测得的数据按照公式(8)进行卡尔曼滤波后,输出以惯性导航模块更新频率为准的定位结果输出。
5、基于RFID标签定位对卡尔曼滤波结果进行修正。
在基于惯性导航模块和速度计利用卡尔曼滤波进行位置预测的同时,实时利用RFID阅读器进行RFID标签的检测。当RFID阅读器检测到RFID标签存在时,执行RFID标签定位,结合接收的信号强度指示值、车辆速度计算车辆与标签之间相对位置和RFID标签内部存储的绝对位置,修正卡尔曼滤波结果,并以此位置继续进行位置推算。RFID定位过程可以采用现有方法,也可以采用本实施例下面介绍的改进方法。
针对现有的基于RFID标签定位的缺陷,本申请的发明人进行了改进。
发明人在研发设计过程中注意到在车辆定位中,阅读器相对于标签往往是单向运动的。于是在本实施例的方案中,发明人巧妙的对该特征进行了利用,提出了创新的基于RFID标签进行定位的定位方案以弥补现有技术的缺陷。
按照现有的RFID标签使用方法,在静态的RFID标签定位场景中,使用三个阅读器可确定唯一位置解,因此必须使用三个阅读器来进行基于RFID标签的定位。
但是在本实施例的一种优选实现方案中,发明人根据车辆定位这种动态的场景,车辆运动具有单向性且速度可得这一特性,减少一个RFID阅读器,即设置两个RFID阅读器,并且也获得了唯一定位解。由于RFID标签沿车道中心线布置,其出现在车的底部区域,因此本实施例中给出了如图3所示的阅读器安装方式,将两个RFID阅读器横向设置在车体两侧,使两个阅读器间距为LR1,R2,距离地面Hc。由于在自动驾驶过程中,基本是以车道中心为基准进行车道保持,因此,RFID标签出现位置主要在左右两侧阅读器之间,如果使用两个阅读器,RFID标签坐标解有两个,分别位于阅读器连线前后。然后,获取车速信息和RFID识别时间,基于车速信息和RFID识别时间可以得到唯一解。
车辆上的RFID阅读器实时监听路面上是否存在RFID标签,设在某时刻检测到标签,并成功接收到RFID存储的绝对位置信息xt和yt以及两个阅读器相应的接收的信号强度指示测量值为RSS1和RSS2。根据图中所示几何关系,设符合自由空间路径损耗模型,标签反射功率Pt,工作频率为f,以R1为例,则在路面形成的候选位置为圆,其半径r1为
r2计算同理。由两圆相交公式,可得两个坐标解为
仿真实例:
(1)仿真参数设置
车辆从原点由静态出发,沿x正方向直线行驶约600m,随后转弯至y正方向,直线行驶约120m,转弯至x负方向直线行驶约600m,进入曲率约为1/60的弯道路段,行驶至原点处停止。整个过程涉及的参数值如表1所示。
表1仿真参数设置及说明
仿真选择了:(1)惯性导航模块单独工作、(2)惯性导航模块/速度计共同工作并利用卡尔曼滤波、(3)惯性导航模块/速度计/RFID直接修正和(4)本发明实施例中的惯性导航模块/速度计/RFID标签定位算法四种方案进行对比,以说明方案设计的有效性。
图4给出了四种方案定位误差随时间变化曲线。
(1)惯性导航模块单独工作时,由于自身测量的加速度噪声存在,定位误差随时间迅速累积,最终累积到4.5m的误差。(2)惯性导航模块/速度计卡尔曼滤波对速度量进行了优化,可以大大降低惯性导航模块的推算误差,但长距离定位时,由于惯性导航模块和速度计噪声的存在,依然会有定位误差随时间的累积,最终达到约0.5m的误差。(3)引入RFID直接修正后,以一定的间距对位置进行修正,可以消除整体上的误差累积趋势,但由于此方案直接以RFID标签存储的绝对位置信息作为当前车辆位置,因此修正时的误差是不确定的,最大可达RFID标签最大通信距离,仿真中局部定位误差可达0.3m左右。(4)本发明实施例中的方案不仅对卡尔曼滤波进行了改进,而且对RFID的定位算法进行了改进,引入了新的RFID标签定位算法,可以有效的提高RFID标签处定位精度,同时使得后续的无RFID标签位置推算结果更加精确,使得本发明的定位方法在整体和局部均可以达到误差小于0.1m高精度的定位结果。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (9)
1.一种融合定位方法,所述方法用于对具有行进轮的目标主体进行定位,所述目标主体具有惯性导航模块和基于行进轮的速度测量模块,其特征在于,所述方法包括:
步骤1:利用所述惯性导航模块和速度测量模块实时测量所述目标主体的加速度和速度;
步骤2:将所测得的速度和加速度作为滤波观测量代入卡尔曼滤波方程组,进行所述目标主体的位置预测,其中,在所述卡尔曼滤波方程组中利用各个滤波状态量之间的噪声协方差矩阵对所述滤波状态量进行修正,所述噪声协方差矩阵基于所述惯性导航模块和速度测量模块的噪声方差确定,
其中,噪声协方差矩阵包括:
(1)过程噪声协方差矩阵Q,
以及(2)观测噪声协方差矩阵R,R=diag(Dv,Dv,Da,Da),其中,Dv为所述速度测量模块的噪声方差。
2.根据权利要求1所述的融合定位方法,其特征在于,在所述目标主体以及其行进路径上设置有基于位置标签的定位模块,所述方法还包括步骤3:利用所述定位模块获取所述目标主体的位置信息,基于所述位置信息对所述位置预测的结果进行修正。
3.根据权利要求1所述的融合定位方法,其特征在于,所述惯性导航模块和速度测量模块的噪声方差在所述目标主体保持静止,行进轮旋转的情况下进行测定。
4.根据权利要求2所述的融合定位方法,其特征在于,所述速度测量模块为速度计,所述定位模块包括设置于所述目标主体的RFID阅读器和设置于路面内的RFID标签,所述目标主体为车辆。
7.根据权利要求4所述的融合定位方法,其中,所述目标主体内沿垂直于其行进方向横向布置两个RFID阅读器,其特征在于:
所述步骤3还包括:
3.1、基于两个所述RFID阅读器获取的RFID标签的信号强度指示值,由路径损耗模型解算出所述RFID标签与所述目标主体的两个可能相对位置;
3.2、基于两个所述RFID阅读器对所述RFID标签的识别时间、车速确定所述RFID标签的相对位置;
3.3、读取所述RFID标签内存储的所述RFID标签的绝对位置信息,计算出所述目标主体的位置信息,
其中,设在某时刻检测到标签,并成功接收到RFID存储的绝对位置信息xt和yt以及两个阅读器相应的接收的信号强度指示测量值为RSS1和RSS2,设符合自由空间路径损耗模型,标签反射功率Pt,工作频率为f,对于任意一个所述RFID阅读器,RFID标签在路面形成的候选位置为圆,其半径r1为
r2计算同理,由两圆相交公式,可得两个坐标解为
8.一种融合定位系统,所述系统用于对具有行进轮的目标主体进行定位,所述融合定位系统包括惯性导航模块、基于行进轮的速度测量模块以及融合定位解算模块,其特征在于,所述融合定位解算模块采用权利要求1-7中任意一项所述的方法进行定位解算。
9.一种计算机可读存储介质,其上存储有计算机程序,其中,该程序被处理器执行时实现根据权利要求1至7中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011163752.1A CN112362052B (zh) | 2020-10-27 | 2020-10-27 | 一种融合定位方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011163752.1A CN112362052B (zh) | 2020-10-27 | 2020-10-27 | 一种融合定位方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112362052A CN112362052A (zh) | 2021-02-12 |
CN112362052B true CN112362052B (zh) | 2022-09-16 |
Family
ID=74510712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011163752.1A Active CN112362052B (zh) | 2020-10-27 | 2020-10-27 | 一种融合定位方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112362052B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113257044B (zh) * | 2021-07-09 | 2022-02-11 | 中航信移动科技有限公司 | Ads-b数据的滤波方法、装置、计算机设备及存储介质 |
CN114964262A (zh) * | 2022-05-24 | 2022-08-30 | 深圳大学 | 一种基于多源融合的机器人定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007459A (zh) * | 2014-05-30 | 2014-08-27 | 北京融智利达科技有限公司 | 一种车载组合定位装置 |
CN105467382A (zh) * | 2015-12-31 | 2016-04-06 | 南京信息工程大学 | 基于svm的多传感器的目标跟踪数据融合算法及其系统 |
CN106767783A (zh) * | 2016-12-15 | 2017-05-31 | 东软集团股份有限公司 | 基于车载通信的定位修正方法和装置 |
CN107289933A (zh) * | 2017-06-28 | 2017-10-24 | 东南大学 | 基于mems传感器和vlc定位融合的双卡尔曼滤波导航装置和方法 |
CN107343052A (zh) * | 2017-07-25 | 2017-11-10 | 苏州大学 | 一种wsn与rfid结合的物流监控系统 |
CN107402005A (zh) * | 2016-05-20 | 2017-11-28 | 北京自动化控制设备研究所 | 一种基于惯性/里程计/rfid的高精度组合导航方法 |
CN109916407A (zh) * | 2019-02-03 | 2019-06-21 | 河南科技大学 | 基于自适应卡尔曼滤波器的室内移动机器人组合定位方法 |
CN110243363A (zh) * | 2019-07-03 | 2019-09-17 | 西南交通大学 | 一种基于低成本imu与rfid技术结合的agv实时定位方法 |
CN111272174A (zh) * | 2020-02-27 | 2020-06-12 | 中国科学院计算技术研究所 | 一种组合导航方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120086598A1 (en) * | 2010-10-08 | 2012-04-12 | Canadian Space Agency | Apparatus and methods for driftless attitude determination and reliable localization of vehicles |
-
2020
- 2020-10-27 CN CN202011163752.1A patent/CN112362052B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007459A (zh) * | 2014-05-30 | 2014-08-27 | 北京融智利达科技有限公司 | 一种车载组合定位装置 |
CN105467382A (zh) * | 2015-12-31 | 2016-04-06 | 南京信息工程大学 | 基于svm的多传感器的目标跟踪数据融合算法及其系统 |
CN107402005A (zh) * | 2016-05-20 | 2017-11-28 | 北京自动化控制设备研究所 | 一种基于惯性/里程计/rfid的高精度组合导航方法 |
CN106767783A (zh) * | 2016-12-15 | 2017-05-31 | 东软集团股份有限公司 | 基于车载通信的定位修正方法和装置 |
CN107289933A (zh) * | 2017-06-28 | 2017-10-24 | 东南大学 | 基于mems传感器和vlc定位融合的双卡尔曼滤波导航装置和方法 |
CN107343052A (zh) * | 2017-07-25 | 2017-11-10 | 苏州大学 | 一种wsn与rfid结合的物流监控系统 |
CN109916407A (zh) * | 2019-02-03 | 2019-06-21 | 河南科技大学 | 基于自适应卡尔曼滤波器的室内移动机器人组合定位方法 |
CN110243363A (zh) * | 2019-07-03 | 2019-09-17 | 西南交通大学 | 一种基于低成本imu与rfid技术结合的agv实时定位方法 |
CN111272174A (zh) * | 2020-02-27 | 2020-06-12 | 中国科学院计算技术研究所 | 一种组合导航方法和系统 |
Non-Patent Citations (2)
Title |
---|
An_Improved_Posteriori_Variance-Covariance_Components_Estimation_Applied_to_Unconventional_GPS_and_Multiple_Low-Cost_Imus_Integration_Strategy;MINGHONG ZHU 等;《SPECIAL SECTION ON EMERGING TRENDS, ISSUES AND CHALLENGES FOR ARRAY SIGNAL PROCEAPPLICATIONS IN SMART CITY》;20190917 * |
基于改进UKF的BDS_IMU组合列车定位方法;蔡煊 等;《西南交通大学学报》;20200430;第55卷(第2期);第393-400页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112362052A (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fan et al. | Data fusion for indoor mobile robot positioning based on tightly coupled INS/UWB | |
JP2022113746A (ja) | 判定装置 | |
CN105628026B (zh) | 一种移动物体的定位定姿方法和系统 | |
CN110702091B (zh) | 一种沿地铁轨道移动机器人的高精度定位方法 | |
US8775063B2 (en) | System and method of lane path estimation using sensor fusion | |
CN103777220B (zh) | 基于光纤陀螺、速度传感器和gps的实时精确位姿估计方法 | |
CN107132563B (zh) | 一种里程计结合双天线差分gnss的组合导航方法 | |
CN110208842A (zh) | 一种车联网环境下车辆高精度定位方法 | |
CN111426320B (zh) | 一种基于图像匹配/惯导/里程计的车辆自主导航方法 | |
CN109752725A (zh) | 一种低速商用机器人、定位导航方法及定位导航系统 | |
CN103487050A (zh) | 一种室内移动机器人定位方法 | |
CN112362052B (zh) | 一种融合定位方法及系统 | |
CN114623823B (zh) | 一种融合里程计的uwb多模态定位系统、方法及装置 | |
WO2016203744A1 (ja) | 測位装置 | |
KR20220024791A (ko) | 차량의 궤적을 결정하기 위한 방법 및 장치 | |
CN112147651B (zh) | 一种异步多车协同目标状态鲁棒估计方法 | |
CN113342059A (zh) | 基于位置和速度误差的多无人机跟踪移动辐射源方法 | |
JP2021085880A (ja) | 移動体の定位誤差の分析 | |
CN118129733A (zh) | 一种典型激光雷达退化场景下的移动机器人定位方法 | |
CN117490683A (zh) | 一种井下隧道多传感器融合算法的定位建图方法 | |
Anousaki et al. | Simultaneous localization and map building of skid-steered robots | |
CN113063441A (zh) | 里程计累计推算误差的数据源纠正方法及装置 | |
Song et al. | RFID/in-vehicle sensors-integrated vehicle positioning strategy utilising LSSVM and federated UKF in a tunnel | |
CN116047480A (zh) | 一种激光雷达到姿态传感器的外参标定方法 | |
CN114666732B (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 |