CN113009173B - 一种高精度的车速计算方法 - Google Patents
一种高精度的车速计算方法 Download PDFInfo
- Publication number
- CN113009173B CN113009173B CN202110140353.1A CN202110140353A CN113009173B CN 113009173 B CN113009173 B CN 113009173B CN 202110140353 A CN202110140353 A CN 202110140353A CN 113009173 B CN113009173 B CN 113009173B
- Authority
- CN
- China
- Prior art keywords
- tire
- acceleration
- vehicle speed
- current moment
- value
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 37
- 230000001133 acceleration Effects 0.000 claims abstract description 162
- 238000012937 correction Methods 0.000 claims abstract description 89
- 238000005259 measurement Methods 0.000 claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 4
- 101100391172 Dictyostelium discoideum forA gene Proteins 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 206010039203 Road traffic accident Diseases 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 102100034112 Alkyldihydroxyacetonephosphate synthase, peroxisomal Human genes 0.000 description 1
- 229920002430 Fibre-reinforced plastic Polymers 0.000 description 1
- 101000799143 Homo sapiens Alkyldihydroxyacetonephosphate synthase, peroxisomal Proteins 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000000848 angular dependent Auger electron spectroscopy Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011151 fibre-reinforced plastic Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W50/0098—Details of control systems ensuring comfort, safety or stability not otherwise provided for
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
- B60W40/10—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
- B60W40/105—Speed
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/02—Knowledge representation; Symbolic representation
- G06N5/022—Knowledge engineering; Knowledge acquisition
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0019—Control system elements or transfer functions
- B60W2050/0028—Mathematical models, e.g. for simulation
- B60W2050/0031—Mathematical model of the vehicle
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0062—Adapting control system settings
- B60W2050/0075—Automatic parameter input, automatic initialising or calibrating means
- B60W2050/0083—Setting, resetting, calibration
- B60W2050/0088—Adaptive recalibration
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/26—Wheel slip
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
- G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mechanical Engineering (AREA)
- Transportation (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Human Computer Interaction (AREA)
- Measuring Fluid Pressure (AREA)
- Tires In General (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本发明涉及一种高精度的车速计算方法,包括计算当前时刻的轮胎滑移率,并判断滑移率是否小于阈值,若是,则继续对轮胎半径和加速度进行修正;否则,停止修正过程;针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式;根据关系式,建立系统预估模型;根据关系式,建立量测系统模型;将系统预估模型和量测系统模型代入卡尔曼滤波器预估模型中,以获取当前时刻的轮胎半径修正值以及加速度偏移量修正值;依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速。通过该方法获得的车速与真实车速基本匹配,误差小,精度高,对辅助驾驶和自动驾驶系统而言,具有更高参考价值,利于提高驾驶的安全性。
Description
技术领域
本发明及车辆技术领域,特别涉及一种高精度的车速计算方法。
背景技术
辅助驾驶(ADAS,Advanced Driving Assistant System)和自动驾驶系统的车辆控制演算法都需要建立在精确的车速计算与加速度两侧基础上,方能够在各种行车情景下达到精准的车辆行为控制。例如,自动跟车功能(ACC,Adaptive Cruise Control),若车速或加速度在量测或计算上存在误差,则会直接反应在本车与前车之鉴的距离控制上,由于距离控制的不精确甚至会造成碰撞,引起严重的交通事故。
目前,在量产车型中最常用于车速计算的参数为轮速感测器所量测的轮胎速度以及加速度计(Accelerometer)所量测的加速度。主流的做法有将非驱动的轮速进行优化处理或者以轮速和加速度两种参数所计算的车速进行叠加(如专利TW106125538中所描述),但无论是对轮速进行优化还是将轮速和加速度进行叠加,这两种方法都没有考虑到轮胎半径和加速度的测量值可能存在的误差。由于车辆在行驶过程中,轮胎不可避免会与地面产生摩擦,这就使得轮胎半径必然是逐渐变小的,若一致以轮胎的出厂半径计算轮速,势必就会引起一定误差。此外,加速度的测量在上坡、下坡、急加速或急减速等情况下,也会受惯性影响而产生偏差。而使用未经优化的轮胎半径以及加速度量测值来估算车速,估算结果是不准确的,参考价值较低。
为了解决现有车速估算方法存在的缺陷,本申请提供了一种高精度的车速计算方法,以对量产车型的轮胎半径以及加速度量测值进行修正,使得车速的计算精度得到提升。
发明内容
为了克服上述缺陷,本发明提供了一种高精度的车速计算方法,其特征在于,包括:
计算当前时刻的轮胎滑移率,并判断滑移率是否小于阈值,若是,则继续对轮胎半径和加速度进行修正;否则,停止修正过程;
针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式;
根据关系式,建立系统预估模型;
根据关系式,建立量测系统模型;
将系统预估模型和量测系统模型代入卡尔曼滤波器预估模型中,以获取当前时刻的轮胎半径修正值以及加速度偏移量修正值;
依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速。
进一步的,所述计算轮胎的滑移率步骤包括:
获取车辆当前时刻的轮胎线速度;
获取车辆当前时刻的中心点移动速度;
获取轮胎线速度与中心点移动速度之间的差值;
若当前时刻的车辆是处于加速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以轮胎线速度作为分母,获取两者的比值,将比值作为滑移率;
若当前时刻的车辆处于减速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以中心点移动速度作为分母,获取两者的比值,将该比值作为滑移率。
进一步的,所述获取车辆当前时刻的轮胎线速度步骤,包括:
获取当前时刻的轮胎角速度;
获取轮胎出厂半径;
以当前时刻的轮胎角速度与轮胎出厂半径的乘积作为当前时刻的轮胎线速度。
进一步的,当滑移率大于等于阈值且停止修正过程之后,还包括:
基于上一时刻的轮胎半径修正值以及加速度偏移量修正值,计算出当前时刻的车速。
进一步的,在所述针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式步骤中,所述关系式如下:
其中,所述为轮胎角加速度,所述ak(i)为加速度量测值,所述ε(i)为加速度偏移量,所述r为轮胎半径。
进一步的,在所述针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式步骤之后,在所述根据关系式,建立系统预估模型步骤之前,还包括公式改写步骤:
将建立的关系式改写成矩阵型公式,所述矩阵型公式如下:
其中,所述为轮胎角加速度,所述ak(i)为加速度量测值,所述ε(i)为加速度偏移量,所述r为轮胎半径。
进一步的,所述根据关系式,建立系统预估模型步骤,包括:
以作为系统预估模型中的预估变量X,设定/>
基于矩阵型公式,按照预估模型的标准式,建立系统预估模型,具体如下:
其中,所述为i时刻的X,所述/>则为i-1时刻的X,所述/>为预测系统杂讯。
进一步的,所述根据关系式,建立量测系统模型步骤,包括:
以作为量测系统模型中的量测值Z,设定/>
基于矩阵型公式,按照量测模型的标准式,建立量测系统模型,具体如下:
其中,所述v(i)为量测系统杂讯。
进一步的,所述预估模型的标准式如下:
X(i)=A·X(i-1)+B·U(i)+N(i),
其中,所述X(i)为i时刻X的状态,所述X(i-1)为i-1时刻的X状态,所述A和B均为系统参数,所述U(·)为i时刻输入的系统变量,所述N(·)为系统杂讯;
所述量测模型的标准式如下:
Z(i)=H·X(i)+V(·),
其中,所述Z(i)为i时刻的量测值,所述X(i)为i时刻X的状态,所述H为系统参数,所述V(i)为系统杂讯。
进一步的,所述依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速步骤,包括:
将当前时刻的加速度量测值与加速度偏移量修正值进行求和计算,并将和值作为当前时刻的加速度修正值;
将当前时刻的加速度修正值与取样周期相乘,并将乘积值与上一时刻的车速相加,所得和值即为当前时刻的车速。
与现有技术相比,本发明的有益效果是如下:
本发明提供了一种高精度的车速计算方法,在该车速计算方法中,通过对同一时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径建立关系式,并将关系式改写成矩阵型公式,然后建立系统预估模型和量测系统模型,并将系统预估模型和量测系统模型代入至卡尔曼滤波器预估模型中,得到轮胎半径修正值以及加速度偏移量修正值。得到的轮胎半径修正值以及加速度偏移量修正值与当前时刻的轮胎真实半径和真实加速度偏移量具有高度吻合性,使得基于轮胎半径修正值以及加速度偏移量修正值计算得到的车速与当前时刻的真实车速基本匹配,误差小,精度高。这对辅助驾驶和自动驾驶系统而言,具有更高的参考价值,利于提高驾驶的安全性,以降低交通事故的发生概率,具有非常高的使用价值。
附图说明
图1为实施例1中高精度的车速计算方法。
图2为实施例1中由GPS获取的加速度量测值、加速度修正值以及未修正的加速度量测值的曲线对比图。
图3为实施例1中由GPS获取的当前时刻的车速、修正后的当前车速以及未修正的当前车速之间的曲线比较图。
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际尺寸;对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的;相同或相似的标号对应相同或相似的部件;附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制。
具体实施方式
下面对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征更易被本领域技术人员理解,从而对本发明的保护范围作出更为清楚的界定。
实施例1
如图1-图3所示,本实施例提供了一种高精度的车速计算方法,该方法的实施基于安装在汽车上的轮速感测器以及设置在汽车上的加速度计。所述轮速感测器用于获取轮胎的转动角速度(由轮胎的转动角速度可以计算出轮胎线速度);而所述加速度计则是用于检测车辆的加速度值,以获取到车辆的加速度量测值。
一种高精度的车速计算方法,包括如下步骤:
101、计算当前时刻的轮胎滑移率,并判断滑移率是否小于阈值,若是,则继续对轮胎半径和加速度进行修正;否则,停止修正过程。
在具体实施过程中,需要先获取车辆当前时刻的轮胎线速度Vtire,然后再获取车辆当前时刻的中心点移动速度VCG(即车速)。随后获取轮胎线速度与中心点移动速度之间的差值(Vtire-VCG)。然后根据当前时刻车辆是处于加速状态还是减速状态,来判断具体使用哪种计算方法来获取车辆的滑移率。一般来讲,若当前时刻的车辆是处于加速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以轮胎线速度作为分母,获取两者的比值,将该比值作为滑移率。也就是说,当车辆处于加速度状态时,S(滑移率)=(Vtire-VCG)/Vtire。而如果当前时刻的车辆处于减速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以中心点移动速度作为分母,获取两者的比值,将该比值作为滑移率。也就是说,如果当前时刻的车辆处于减速状态时,则S(滑移率)=(Vtire-VCG)/VCG。
本技术方案中,为了获取车辆当前时刻的轮胎线速度,一般需要首选获取到当前时刻的轮胎角速度,轮胎角速度通过轮速感测器来获取。然后获取轮胎出厂半径,轮胎出厂半径通常由厂家直接提供。最后,再以当前时刻的轮胎角速度与轮胎出厂半径的乘积作为当前时刻的轮胎线速度。有时候,ESC(车身稳定控制系统)会直接把四个轮胎的角速度传输至CAN(控制器局域网络、Controller Area Network)上,由相关的计算模块通过四个轮的角速度以及出厂轮胎半径参数,计算出平均轮胎线速度;当然,ESC也可以直接将轮胎的角速度转换成轮胎线速度,在传输至CAN上,在此不做唯一性限制。
本技术方案主要是针对量产车的车速估计,由于量产车上通常还没有安装GPS系统,所以无法通过GPS系统直接获取到车速。此时,为了计算出车辆的滑移率,对后续修正做出准确判断,就需要我们根据轮胎线速度以及车辆加速度来推算出当前时刻的车速,这里所说的车辆加速度是由加速度车辆所得。具体的,依据轮胎线速度以及车辆加速度来估算车速的计算公式如下:
其中,VCG(i)为i时刻的车速,而VCG(i-1)为i-1时刻的车速,所述Vrr、Vlr、Vrf、Vlf分别为四个轮胎的轮胎线速度,所述Krr、Klr、Krf、Klf为各轮胎的权重值,所述a(i)为i时刻的加速度,所述T为取样周期,所述Kacc为加速度权重值,所述Ks=Krr+Klr+Krf+Klf+Kacc。基于上述估算出的车速以及推算出的的轮胎线速度,就可以完成对滑移率的计算。
作为优选的,当滑移率大于等于阈值且停止修正过程之后,还可以调用上一时刻的轮胎半径修正值以及加速度偏移量修正值,并基于上一时刻的轮胎半径修正值以及加速度偏移量修正值,计算出当前时刻的车速。本实施例中,滑移率的阈值通常设为3%-8%之间,优选5%。也就是说,当上述所计算出的滑移率S小于阈值,则继续对加速度以及轮胎半径进行修正,若滑移率S大于或等于阈值时,则会暂停对加速度的修正,并使用上一时刻的轮胎半径修正值以及加速度偏移量修正值来对车速进行估算。根据多次试验观察,一般用路环境轮胎不会处于长时间打滑,加速度还是有很多机会得到修正,而且这个算法的修正主要是为了消除坡度对加速度量测的误差,在打滑瞬间车辆所处的坡度与上一个时刻相比通常不会有太大变化;即使在停止对加速度修正的这段时间里,车辆行驶的坡度有变化,这样的计算结果依然比完全不修正来得准确。
102、针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式。
在建立关系式的过程中,需要先对轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径之间的关系做初步分析。其中,考虑到车速计算的来源主要有有两个,一个是轮胎线速度,另一个则是车辆中心的加速度。所以为了使车速估算精确,就必须分别对这两个来源参数进行算法上的优化。已知的现象是加速度的量测会受到车辆所在的坡度影响,我们要透过系统的预测模型使得加速度更符合使用需求。另一个会影响车速计算的参数是轮胎半径,通常轮胎线速度的量测是先由角速度ωtire而来,轮胎线速度Vtire=ωtire·r,由此可见r的误差也会影响轮胎线速度。其中,ε(i)=ak(i)-a(i),所述ε(i)为加速度偏移量(测量值与真实值之间存在的误差),所述ak(i)为加速度量测值,而a(i)则视为i时刻的真实加速度值,也就是我们想要获取的加速度数值。
假设车辆在坡面上行驶使,加速度获取的车辆加速度值已经收到了坡面的影响,那么此时为了修正加速度,我们就需要通过轮胎线速度的微分来获取加速度,藉此作为另一个加速度来源的参考。把由轮胎线速度微分所得到的加速度记作atire(i),由物理公式可知轮胎线性加速度atire(i)与轮胎角加速度的关系为/>然后,我们假设在车辆模型中,轮胎线加速度atire(i)等于车辆中心加速度a(i),并带入加速度偏移量ε(i),就可得如下关联公式:
进而推导出如下所建立的关系式:
其中,所述为轮胎角加速度,所述ak(i)为加速度量测值,所述ε(i)为加速度偏移量,所述r为轮胎半径。
为了便于后续建立系统预估模型和量测系统模型,以便成功代入卡尔曼滤波器预估模型中,还需要将上述建立的关系式进一步改写成矩阵型公式,所述矩阵型公式具体如下:
其中,在矩阵型公式中,所述为轮胎角加速度,所述ak(i)为加速度量测值,所述ε(i)为加速度偏移量,所述r为轮胎半径。
103、根据关系式,建立系统预估模型。
在建立系统预估模型的过程中,首先需要确立预估变量。本实施例中,以作为系统预估模型中的预估变量X,即设定/>然后,再基于矩阵型公式,并按照预估模型的标准式,建立系统预估模型,所述系统预估模型具体如下:
其中,在系统预估模型中,所述为i时刻的X,所述/>则为i-1时刻的X,所述/>为预测系统杂讯。
本技术方案中,所说的预估模型的标准式如下所示:
X(i)=A·X(i-1)+B·U(i)+N(i),
其中,在预估模型的标准式中,所述X(i)为i时刻X的状态,所述X(i-1)为i-1时刻的X状态,所述A和B均为系统参数,所述U(i)为i时刻输入的系统变量(本技术的方案中,U=0),所述N(i)为系统杂讯。
104、根据关系式,建立量测系统模型。
在建立量测系统模型的过程中,首先需要确定量测值。本实施例中,以作为量测系统模型中的量测值Z,设定/>然后,再基于矩阵型公式,按照量测模型的标准式,建立量测系统模型,具体如下:
其中,在量测系统模型中,所述v(i)为量测系统杂讯。
本技术方案中,所说的量测模型的标准式如下所示:
Z(i)=H·X(i)+V(i),
其中,所述Z(i)为i时刻的量测值,所述X(i)为i时刻X的状态,所述H为系统参数(H=[ak(i) -1]),所述V(i)为系统杂讯。
105、将系统预估模型和量测系统模型代入卡尔曼滤波器预估模型中,以获取当前时刻的轮胎半径修正值以及加速度偏移量修正值。
具体的,将步骤103和步骤104中建立的系统预估模型和量测系统模型带入至卡尔曼滤波器预估模型中。由于卡尔曼滤波器预估模型已经形成了固定的算法模式,也已被广泛使用,所以在此不做赘述。然后由系统预估模型中的预测系统杂讯以及量测系统模型中的量测系统杂讯v(i)来计算出卡尔曼滤波器预估模型中的参数Q和R,具体计算方式套用统计学中的协方差公式即可,在此不做详述。最后通过常规的演算法计算每个时刻的预估变量X。也就是说,每个时刻的轮胎半径r与加速度偏移量ε都会被估算出来,即可得到每一时刻的轮胎半径修正值和加速度偏移量修正值。
106、依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速。
当得到了当前时刻的轮胎半径修正值及加速度偏移量修正值之后,便可以计算当前时刻的车速了。具体的,将当前时刻的加速度量测值与加速度偏移量修正值进行求和(或求差)计算,即可获得该时刻较为准确的加速度修正值acorrection(i)。将得到的加速度修正值、由加速度计获取的未经修正的加速度量测值以及由GPS获取的加速度量测值(视为准确的加速度值)进行对比,如图2所示。由图2可知,由加速度计获取的未经修正的加速度量测值与由GPS获取的加速度量测值之间存在较为明显的误差,而加速度修正值与由GPS获取的加速度量测值之间则误差非常小,两条曲线具有高度吻合性,这就表明,加速度修正值与由GPS获取的加速度量测值几乎是吻合的,得到的加速度修正值是能够代表真实加速度值来估算车速的。
为了进一步获取到当前时刻的车速,还需要将当前时刻的加速度修正值acorrection(i)与取样周期T相乘,并将乘积值与上一时刻的车速相加,所得和值即为当前时刻的车速。具体如车速计算公式一所示。
车速计算公式一:VCG(i)=VCG(i-1)+acorrection(i)*T
其中,VCG(i)为i时刻的车速,而VCG(i-1)为i-1时刻的车速。
当然,对车速的计算方式并不唯一,也可以选择将轮胎半径修正值r直接带入至车速计算公式二中,以得到较为精确的车速计算值。
车速计算公式二:Vtire=ωtire·r
其中,Vtire为轮胎线速度,ωtire为轮胎角速度,而r为轮胎半径修正值。
基于轮胎半径修正值以及加速度修正值得到的车速与未修正的当前车速以及由GPS获取的当前时刻的车速进行比较,发现基于轮胎半径修正值以及加速度修正值得到的车速与GPS获取的当前时刻的车速吻合度更高,两者之间的误差更小,其参考价值更大,具体如图3所示。
本实施例提供的车速计算方法,通过对轮胎角加速度、加速度测量值、加速度偏移量以及轮胎半径之间的关系进行分析,建立了具体的关系式(在步骤102中有详细描述)。并通过对建立的关系式进行简单改写后形成矩阵型公式(在步骤102中有详细描述)。形成的矩阵型公式不仅能够改写成与预估模型的标准式相符合的格式(即系统预估模型),也能够改写成与量测模型的标准式相符的格式(即量测系统模型)。而且,最重要的是,形成的系统预估模型和量测系统模型均与卡尔曼滤波器预估模型的特性相吻合,刚好都能够把各变量放在正确的位置,并完成对各参数矩阵的合理设置。这就表明利用卡尔曼滤波器预估模型计算轮胎半径修正值以及加速度偏移量修正值是比较适合。当然,在具体实验过程中,通过采集准确的轮胎半径以及加速度偏移量,来验证卡尔曼滤波器预估模型输出的轮胎半径修正值以及加速度偏移量修正值,发现卡尔曼滤波器预估模型输出的结果的确也是可靠有效的。
值得注意的是,在本技术方案中,卡尔曼滤波器预估模型并不是计算轮胎半径修正值以及加速度偏移量修正值的唯一选择,在实际使用过程中,也可以采用最小平方公式(Least squares formula)替代卡尔曼滤波器预估模型来完成计算。由于最小平方公式(Least squares formula)也是非常成熟的计算方法,在此不再做具体赘述。
本实施例公开了一种高精度的车速计算方法,在该车速计算方法中,通过对同一时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径建立关系式,并将关系式改写成矩阵型公式,然后建立系统预估模型和量测系统模型,并将系统预估模型和量测系统模型代入至卡尔曼滤波器预估模型中,得到轮胎半径修正值以及加速度偏移量修正值。得到的轮胎半径修正值以及加速度偏移量修正值与当前时刻的轮胎真实半径和真实加速度偏移量具有高度吻合性,使得基于轮胎半径修正值以及加速度偏移量修正值计算得到的车速与当前时刻的真实车速基本匹配,误差小,精度高。这对辅助驾驶和自动驾驶系统而言,具有更高的参考价值,利于提高驾驶的安全性,以降低交通事故的发生概率,具有非常高的使用价值。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (9)
1.一种高精度的车速计算方法,其特征在于,包括:
计算当前时刻的轮胎滑移率,并判断滑移率是否小于阈值,若是,则继续对轮胎半径和加速度进行修正;否则,停止修正过程;
针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式;
将建立的关系式改写成矩阵型公式,所述矩阵型公式如下:
,
其中,所述为轮胎角加速度,所述/>为加速度量测值,所述/>为加速度偏移量,所述/>为轮胎半径;
根据关系式,建立系统预估模型;
根据关系式,建立量测系统模型;
将系统预估模型和量测系统模型代入卡尔曼滤波器预估模型中,以获取当前时刻的轮胎半径修正值以及加速度偏移量修正值;
依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速。
2.如权利要求1所述一种高精度的车速计算方法,其特征在于,所述计算轮胎的滑移率步骤包括:
获取车辆当前时刻的轮胎线速度;
获取车辆当前时刻的中心点移动速度;
获取轮胎线速度与中心点移动速度之间的差值;
若当前时刻的车辆是处于加速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以轮胎线速度作为分母,获取两者的比值,将比值作为滑移率;
若当前时刻的车辆处于减速状态,则以轮胎线速度与中心点移动速度之间的差值作为分子,以中心点移动速度作为分母,获取两者的比值,将该比值作为滑移率。
3.如权利要求2所述一种高精度的车速计算方法,其特征在于,所述获取车辆当前时刻的轮胎线速度步骤,包括:
获取当前时刻的轮胎角速度;
获取轮胎出厂半径;
以当前时刻的轮胎角速度与轮胎出厂半径的乘积作为当前时刻的轮胎线速度。
4.如权利要求1所述一种高精度的车速计算方法,其特征在于,当滑移率大于等于阈值且停止修正过程之后,还包括:
基于上一时刻的轮胎半径修正值以及加速度偏移量修正值,计算出当前时刻的车速。
5.如权利要求1所述一种高精度的车速计算方法,其特征在于,在所述针对当前时刻的轮胎角加速度、加速度量测值、加速度偏移量以及轮胎半径,建立关系式步骤中,所述关系式如下:
其中,所述为轮胎角加速度,所述/>为加速度量测值,所述/>为加速度偏移量,所述/>为轮胎半径。
6.如权利要求1所述一种高精度的车速计算方法,其特征在于,所述根据关系式,建立系统预估模型步骤,包括:
以作为系统预估模型中的预估变量X,设定X=/>;
基于矩阵型公式,按照预估模型的标准式,建立系统预估模型,具体如下:
,
其中,所述为/>时刻的X,所述/>则为/>时刻的X,所述/>为预测系统杂讯。
7.如权利要求6所述一种高精度的车速计算方法,其特征在于,所述根据关系式,建立量测系统模型步骤,包括:
以作为量测系统模型中的量测值Z,设定Z=/>;
基于矩阵型公式,按照量测模型的标准式,建立量测系统模型,具体如下:
,
其中,所述为量测系统杂讯。
8.如权利要求7所述一种高精度的车速计算方法,其特征在于,所述预估模型的标准式如下:
,
其中,所述为/>时刻X的状态,所述/>为/>时刻的X状态,所述A和B均为系统参数,所述/>为/>时刻输入的系统变量,所述/>为系统杂讯;
所述量测模型的标准式如下:
,
其中,所述为i时刻的量测值,所述/>为/>时刻X的状态,所述H为系统参数,所述为系统杂讯。
9.如权利要求1所述一种高精度的车速计算方法,其特征在于,所述依据当前时刻的轮胎半径修正值及加速度偏移量修正值,计算当前时刻的车速步骤,包括:
将当前时刻的加速度量测值与加速度偏移量修正值进行求和计算,并将和值作为当前时刻的加速度修正值;
将当前时刻的加速度修正值与取样周期相乘,并将乘积值与上一时刻的车速相加,所得和值即为当前时刻的车速。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110140353.1A CN113009173B (zh) | 2021-02-02 | 2021-02-02 | 一种高精度的车速计算方法 |
PCT/CN2021/102532 WO2022166083A1 (zh) | 2021-02-02 | 2021-06-25 | 一种高精度的车速计算方法 |
US18/270,504 US20240077510A1 (en) | 2021-02-02 | 2021-06-25 | High-precision vehicle velocity calculation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110140353.1A CN113009173B (zh) | 2021-02-02 | 2021-02-02 | 一种高精度的车速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113009173A CN113009173A (zh) | 2021-06-22 |
CN113009173B true CN113009173B (zh) | 2023-10-03 |
Family
ID=76384838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110140353.1A Active CN113009173B (zh) | 2021-02-02 | 2021-02-02 | 一种高精度的车速计算方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240077510A1 (zh) |
CN (1) | CN113009173B (zh) |
WO (1) | WO2022166083A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009173B (zh) * | 2021-02-02 | 2023-10-03 | 惠州市德赛西威汽车电子股份有限公司 | 一种高精度的车速计算方法 |
CN116039648B (zh) * | 2023-04-03 | 2023-06-27 | 成都赛力斯科技有限公司 | 一种基于权重的坡度计算方法、装置及车辆 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014167929A1 (ja) * | 2013-04-12 | 2014-10-16 | 古野電気株式会社 | 速度算出装置、及び速度算出方法 |
CN111806430A (zh) * | 2020-06-03 | 2020-10-23 | 惠州市德赛西威汽车电子股份有限公司 | 一种用于自动泊车的车速计算方法 |
CN111912426A (zh) * | 2020-08-10 | 2020-11-10 | 电子科技大学 | 一种基于mems imu的低成本里程计设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009173B (zh) * | 2021-02-02 | 2023-10-03 | 惠州市德赛西威汽车电子股份有限公司 | 一种高精度的车速计算方法 |
-
2021
- 2021-02-02 CN CN202110140353.1A patent/CN113009173B/zh active Active
- 2021-06-25 WO PCT/CN2021/102532 patent/WO2022166083A1/zh active Application Filing
- 2021-06-25 US US18/270,504 patent/US20240077510A1/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014167929A1 (ja) * | 2013-04-12 | 2014-10-16 | 古野電気株式会社 | 速度算出装置、及び速度算出方法 |
CN111806430A (zh) * | 2020-06-03 | 2020-10-23 | 惠州市德赛西威汽车电子股份有限公司 | 一种用于自动泊车的车速计算方法 |
CN111912426A (zh) * | 2020-08-10 | 2020-11-10 | 电子科技大学 | 一种基于mems imu的低成本里程计设计方法 |
Non-Patent Citations (2)
Title |
---|
丁能根 等.基于卡尔曼滤波和ABS控制输入的车速估计.北京航空航天大学学报.2011,37(01),第67-71页. * |
芦冰 等.UKF车速估计器的算法研究与仿真.吉林大学学报(信息科学版).2015,第33卷(第1期),第7-11页. * |
Also Published As
Publication number | Publication date |
---|---|
CN113009173A (zh) | 2021-06-22 |
WO2022166083A1 (zh) | 2022-08-11 |
US20240077510A1 (en) | 2024-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113009173B (zh) | 一种高精度的车速计算方法 | |
US11472236B2 (en) | Method, control device, and system for determining a profile depth of a profile of a tire | |
KR101884485B1 (ko) | 자동차의 스티어링 각도의 결정 | |
EP3115765B1 (en) | Tire sensor-based vehicle state estimation system and method | |
EP2777999B1 (en) | Dynamic tire slip angle estimation system and method | |
KR100222506B1 (ko) | 제동 차량의 속도 추정 방법 | |
US20140371990A1 (en) | Sensor system comprising a vehicle model unit | |
CN109606378A (zh) | 面向非高斯噪声环境的车辆行驶状态估计方法 | |
CN113002549B (zh) | 一种车辆状态估算方法、装置、设备及存储介质 | |
TW201940370A (zh) | 基於車輛量測資料處理之車輛操作 | |
CN111086520A (zh) | 一种适用于四轮驱动车辆多轮高滑移率时的速度估计算法 | |
US6834543B2 (en) | Underinflation detector and method of correcting output of sensor using the same | |
JP2011201404A (ja) | 車両特性情報推定装置及びそれを用いた警告装置 | |
JP3271952B2 (ja) | 車両の路面摩擦係数推定装置 | |
JP7368181B2 (ja) | 摩耗量推定システムおよび演算モデル生成システム | |
US6865456B2 (en) | Underinflation detector | |
CN116767240A (zh) | 车轮轮速修正方法、装置及电子设备 | |
Goos et al. | Model-based state estimator for an intelligent tire | |
JP5002203B2 (ja) | タイヤ動荷重半径の基準値初期化方法 | |
CN116409327A (zh) | 侧向工况下考虑轮胎瞬态特性的路面附着系数估计方法 | |
WO2022203026A1 (ja) | 運転制御装置 | |
CN114043986B (zh) | 一种考虑质量失配的轮胎路面附着系数多模型融合估计方法 | |
JP2001116569A (ja) | 距離係数算出方法と、移動距離算出装置 | |
JP2008247243A (ja) | タイヤの内圧低下判断方法における閾値設定方法 | |
CN113771865B (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 |