CN103968910B - 一种提高机动车燃油量估算精度的方法 - Google Patents
一种提高机动车燃油量估算精度的方法 Download PDFInfo
- Publication number
- CN103968910B CN103968910B CN201410196768.0A CN201410196768A CN103968910B CN 103968910 B CN103968910 B CN 103968910B CN 201410196768 A CN201410196768 A CN 201410196768A CN 103968910 B CN103968910 B CN 103968910B
- Authority
- CN
- China
- Prior art keywords
- oil mass
- data
- value
- filtering
- oil
- 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
Landscapes
- Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明公开了一种提高机动车燃油量估算精度的方法,包括步骤一:油量检测设备采集油量电压信号;步骤二:复合采样滤波;步骤三:通过GPRS向后台监控中心发送油量信息;步骤四:判断是否需要进行油量补偿,如果是,进入步骤五,否则,进入步骤六;步骤五:进行油量补偿;步骤六:通过油量检测设备上传的油量电压信号,得出此时车辆油箱中的油量值。本发明提高了油量监测精度、提高了油量信号的稳定性。
Description
技术领域
本发明涉及一种用于提高机动车燃油量估算精度的方法,特别适合对货车、客车等大容积油箱车辆的燃油量监测精度进行改善。
背景技术
油量监测系统主要由油量检测设备、GPRS传输网络和中心监控软件三部分组成。其中,油量检测设备采集车辆油量传感器的电压信号,并结合GPS信息,通过GPRS网络发送到后台监控中心软件,最终在后台软件中显示里程-油量关系图和时间-油量关系图。
基于原车油量传感器的油量监测系统,其油量监测精度只能达到70%左右,这主要是受如下两个子问题的影响。
子问题1.车辆在实际行驶过程中会受加减速、转弯、路面波动以及车身电瓶供压不稳等因素的影响,致使油量检测设备采集到的油量信息含有大量杂波,极大地影响了油量监测精度。
子问题2.现有技术在实际应用中会出现检测盲区的问题。如图1所示,当车辆油箱中的油面高于传感器的最大检测范围时,会造成油量检测失真的问题。在后台监控软件中会出现如图2所示的里程-油量图,在A到B这段里程变化中,其对应的油量值一直保持不变,此即为油量检测盲区。该问题是由原车油量传感器自身测量范围有限引起的,对油量监测精度会产生较大影响。
发明内容
本发明的目的是为了提高油量监测稳定性及精度,依次解决两个子问题。对于子问题1,本发明中的油量检测设备采用了一种复合采样滤波方法,可以有效滤除车辆在行驶过程中受到的各种油量干扰信号,最终通过GPRS网络每30秒钟向监控中心发送一条油量信号。油量检测设备所用的单片机型号为STM32F107,GPRS模块的型号为SIM900D。
在后台监控中心软件对数据处理的过程中,通过一套基于里程和百公里油耗的油量补偿方法,可以解决子问题2,进而提高油量监测的精度。里程信息通过GPS模块及单片机计算得到,所用GPS模块型号为HOLUX M-8929。
一种提高机动车燃油量估算精度的方法,包括以下几个步骤:
步骤一:油量检测设备采集油量电压信号;
油量检测设备通过油量信号线采集原车油量传感器和仪表盘之间的油量电压信号;
步骤二:复合采样滤波;
每1ms采集一次油量电压信号,然后对其进行滤波,每100ms得到一条油量电压数据;
步骤三:通过GPRS向后台监控中心发送油量信息;
对同一时刻采集到的油量电压数据、位置、海拔、速度及里程信息进行整合,通过GPRS模块每隔30秒向监控中心发送一条数据;
步骤四:判断是否需要进行油量补偿,如果是,进入步骤五,否则,进入步骤六;
监控中心对相邻两条油量电压数据中的油量值作差;设ΔYi=Yi+1-Yi,Yi表示油量值,如果从ΔYi到ΔYi+m的m+1个连续值都为0且m的值大于100,则从Yi到Yi+m的油量值保持不变,进行油量补偿,进入步骤五,否则不用补偿,进入步骤六;
步骤五:进行油量补偿;
根据该车最新的百公里油耗及待补偿油量区间段对应的里程值,进行油量补偿,得到修正油量;
步骤六:通过油量检测设备上传的油量电压信号,得出此时车辆油箱中的油量值。
本发明的优点在于:
(1)本发明提高了油量监测精度;
(2)本发明提高了油量信号的稳定性;
附图说明
图1是背景技术中油量检测盲区示意图;
图2是背景技术中里程-油量图;
图3本发明整体技术方案流程图;
图4本发明采样滤波方法流程图;
图5本发明油量补偿方法流程图;
图6本发明油量估算方法处理前的里程-油量图;
图7本发明油量估算方法处理后的里程-油量图;
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种提高机动车燃油量估算精度的方法,流程如图3所示,包括以下几个步骤:
步骤一:油量检测设备采集油量电压信号;
油量检测设备通过油量信号线采集原车油量传感器和仪表盘之间的油量电压信号。
步骤二:复合采样滤波;
每1ms采集一次油量电压信号,然后对其进行滤波,最后每100ms得到一条油量电压数据。
步骤三:通过GPRS向后台监控中心发送油量信息;
对同一时刻采集到的油量电压数据、位置、海拔、速度及里程信息进行整合,通过GPRS模块每隔30秒向监控中心发送一条数据。
步骤四:判断是否需要进行油量补偿,如果是,进入步骤五,否则,进入步骤六;
对相邻两条油量电压数据中的油量值作差。定义ΔYi=Yi+1-Yi,Yi表示油量值,如果从ΔYi到ΔYi+m的m+1个连续值都为0且m的值大于100,则认为从Yi到Yi+m这一段油量值是保持不变的,即为需要进行油量补偿的区间段,进入步骤五,反之则不用补偿,进入步骤六。
步骤五:进行油量补偿;
根据该车最新的百公里油耗及待补偿油量区间段对应的里程值,进行油量补偿,得到修正油量。
步骤六:通过油量检测设备上传的油量电压信号,得出此时车辆油箱中的油量值。
所述的步骤二复合采样滤波中,具体的复合采样滤波方法为:
针对车辆加减速、转弯、启停等因素所引起的大幅值脉冲杂波和车身供压不稳等因素所引起的高频杂波,本发明提出了一套复合采样滤波方法,可以对车辆在行驶过程中受到的各种油量信号干扰起到良好的过滤效果。
如图4所示,复合采样滤波方法包括如下步骤:
步骤1:采样
油量检测设备每1ms采集一次油量信号,然后对100ms内采集到的100条油量信号进行最大值采样。
步骤2:滤波
步骤2包括以下三个子步骤:
步骤2.1:限幅递推中位值平均滤波
通过步骤1的采样,单片机每100ms接收一条数据,把连续采集的10个数据看成一个队列,新接收的数据与上一条数据作差并取绝对值,如果结果大于等于0.2V,则认为新数据是噪声干扰,进行抛除;如果结果小于0.2V就放入队尾,同时抛除原来队首的数据。然后抛掉这10个数据中的最大值和最小值,再对剩下的八个数据求均值,最后用该平均值代替这条新接收的数据。
通过上述滤波处理,得到油量数据组A。
步骤2.2:一阶滞后滤波
经过步骤2.1的滤波处理后,再对数据组A进行一阶滞后滤波,即采用本次采样值与上次滤波输出值进行加权,得到有效滤波值,使得输出对输入有反馈作用。一阶滞后滤波表达式为:
Bn=α·An+(1-α)·Bn-1
式中,Bn表示本次滤波结果,Bn-1表示上次滤波结果,An表示对数据组A的新采样值,α为滤波系数,根据实际经验,取0.5。
最后,通过一阶滞后滤波得到油量数据组B。
步骤2.3:卡尔曼滤波
经过步骤2.2的滤波处理后,再对数组B进行卡尔曼滤波。记数据组B中的第k条油量电压值为观测值zk,由卡尔曼滤波计算得到的第k条油量电压值为xk。
根据车辆油箱内油量的液位变化特性可知,油面在短时间内是不变的,故可得状态方程为:
xk=xk-1+ωk-1
量测方程为:
zk=xk+υk
式中,xk为状态变量,zk为观测变量,ωk和υk分别为过程激励噪声和观测噪声,它们为相互独立、正态分布的白噪声,即p(ω)~N(0,Q),p(υ)~N(0,R),Q、R分别为过程激励噪声和观测噪声的方差。
定义xk -(-代表先验,^代表估计)为在已知第k步之前状态情况下第k步的先验状态估计。定义xk为已知测量变量zk时第k步的后验状态估计。由此定义先验估计误差和后验估计误差:
ek -=xk-xk -
ek=xk-xk
先验估计误差的协方差:Pk -=E[ek -ek -T]。
后验估计误差的协方差:Pk=E[ekek T]。
用K表示后验估计误差协方差的最小值,则根据卡尔曼滤波器的概率原型,先验估计xk -和加权的测量变量zk及其预测Hxk -之差的线性组合便可构成后验状态估计xk。即xk=xk -+K(zk-Hxk -)。通常我们将K称为卡尔曼增益,有
根据车辆油箱内的油位变化特性可知,油面在短时间内是不变的,故可认为油面的变化速度及变化加速度为0,进而得到时间更新方程为:
xk -=xk-1
Pk -=Pk-1+Q
测量更新方程为:
xk=xk -+Kk(zk-xk -)
Pk=(I-Kk)Pk -
根据实际经验,油箱液位的真值可近似为正态分布,其期望为油箱高度的一半,标准差为1。记油箱油量为油箱容积一半时所对应的油量电压信号值为Uh,则我们可确定滤波器的初始条件x0=Uh。由卡尔曼滤波的实际经验知,几乎任何P0≠0都会使滤波器最终收敛。这里我们令P0=1。过程激励噪声方差Q=10-5,观测噪声方差R取0.01。代入以上初始值,即可完成对油量电压数据的卡尔曼滤波处理,得到的xk即为经过卡尔曼滤波后的第k条油量电压数据。
所述的步骤五为油量补偿方法,流程如图5所示,包括如下步骤:
步骤1:计算该车的百公里油耗K。
在里程-油量图中,每个点都可以用(里程,油量)坐标进行唯一表示,用Xi表示第i条数据的里程,Yi表示第i条数据的油量,其中Xi的单位为km,Yi的单位为升。
记Xn=Xi+m+1+100,Yn为里程Xn下的平均油量,据此可求出该车的百公里油耗
步骤2:根据盲区内各点的里程及百公里油耗K进行油量修正。
在需要进行油量补偿的区间段内,第i条数据的修正油量第i+1条数据的修正油量……第i+m-1条数据的修正油量
本发明中提高油量估算精度的方法,可以有效过滤掉油量信号中掺杂的大幅值脉冲杂波,高频杂波的振动幅值平均下降了80%,同时对盲区内的油量进行了补偿,最终解决了影响油量监测精度的两个子问题,将油量监测精度从70%左右提高到95%左右。图6和图7分别为处理前后的里程-油量曲线示意图。
本发明提出了一种提高机动车燃油量估算精度的方法,包括一套可以有效过滤掉车辆在行驶过程中受到的各种信号干扰的复合采样滤波方法,以及一种针对油量检测盲区问题的油量补偿方法。通过这两个方法,对油量检测设备采集到的油量电压数据进行处理,最终提高了机动车燃油量的数据稳定性及油量估算精度。
Claims (2)
1.一种提高机动车燃油量估算精度的方法,包括以下几个步骤:
步骤一:油量检测设备采集油量电压信号;
油量检测设备通过油量信号线采集原车油量传感器和仪表盘之间的油量电压信号;
步骤二:复合采样滤波;具体包括如下步骤:
步骤1:采样
每1ms采集一次油量电压信号,然后对其进行滤波,每100ms得到一条油量电压数据;
油量检测设备每1ms采集一次油量信号,然后对100ms内采集到的100条油量信号进行最大值采样;
步骤2:滤波,包括以下三个子步骤:
步骤2.1:限幅递推中位值平均滤波;
通过步骤1的采样,每100ms接收一条数据,把连续采集的10个数据设为一个队列,新接收的数据与上一条数据作差并取绝对值,如果结果大于等于0.2V,则认为新数据是噪声干扰,进行抛除;如果结果小于0.2V就放入队尾,同时抛除原来队首的数据;然后抛掉10个数据中的最大值和最小值,对剩下的八个数据求均值,最后用该平均值代替这条新接收的数据,通过上述滤波处理,得到油量数据组A;
步骤2.2:一阶滞后滤波
经过步骤2.1的滤波处理后,再对数据组A进行一阶滞后滤波,即采用本次采样值与上次滤波输出值进行加权,得到有效滤波值;一阶滞后滤波表达式为:
Bn=αAn+(1-α)·Bn-1
式中,Bn表示本次滤波结果,Bn-1表示上次滤波结果,An表示对数据组A的新采样值,α为滤波系数,最后,通过一阶滞后滤波得到油量数据组B;
步骤2.3:卡尔曼滤波
经过步骤2.2的滤波处理后,再对数组B进行卡尔曼滤波;记数据组B中的第k条油量电压值为观测值zk,由卡尔曼滤波计算得到的第k条油量电压值为xk;
状态方程为:
xk=xk-1+ωk-1
量测方程为:
zk=xk+υk
式中,xk为状态变量,zk为观测变量,ωk和υk分别为过程激励噪声和观测噪声,p(ω)~N(0,Q),p(υ)~N(0,R),Q、R分别为过程激励噪声和观测噪声的方差;
设xk -为在已知第k步之前状态情况下第k步的先验状态估计,xk为已知测量变量zk时第k步的后验状态估计,其中-代表先验,^代表估计,则先验估计误差和后验估计误差:
ek -=xk-xk -
ek=xk-xk
先验估计误差的协方差:Pk -=E[ek -ek -T];
后验估计误差的协方差:Pk=E[ekek T];
用K表示后验估计误差协方差的最小值,则根据卡尔曼滤波器的概率原型,先验估计xk -和加权的测量变量zk及其预测Hxk -之差的线性组合便可构成后验状态估计xk;即xk=xk -+K(zk-Hxk -);K为卡尔曼增益,
设油面的变化速度及变化加速度为0,时间更新方程为:
xk -=xk-1
Pk -=Pk-1+Q
测量更新方程为:
xk=xk -+Kk(zk-xk -)
Pk=(I-Kk)Pk -
设油箱液位的真值为正态分布,其期望为油箱高度的一半,标准差为1;记油箱油量为油箱容积一半时所对应的油量电压信号值为Uh,滤波器的初始条件x0=Uh;令P0=1;过程激励噪声方差Q=10-5,观测噪声方差R取0.01;代入以上初始值,完成对油量电压数据的卡尔曼滤波处理,得到的xk即为经过卡尔曼滤波后的第k条油量电压数据;
步骤三:通过GPRS向后台监控中心发送油量信息;
对同一时刻采集到的油量电压数据、位置、海拔、速度及里程信息进行整合,通过GPRS模块每隔30秒向监控中心发送一条数据;
步骤四:判断是否需要进行油量补偿,如果是,进入步骤五,否则,进入步骤六;
监控中心对相邻两条油量电压数据中的油量值作差;设ΔYi=Yi+1-Yi,Yi表示油量值,如果从ΔYi到ΔYi+m的m+1个连续值都为0且m的值大于100,则从Yi到Yi+m的油量值保持不变,进行油量补偿,进入步骤五,否则不用补偿,进入步骤六;
步骤五:进行油量补偿;
根据该车最新的百公里油耗及待补偿油量区间段对应的里程值,进行油量补偿,得到修正油量;
步骤六:通过油量检测设备上传的油量电压信号,得出此时车辆油箱中的油量值。
2.根据权利要求1所述的一种提高机动车燃油量估算精度的方法,所述的步骤五进行油量补偿方法,包括如下步骤:
步骤1:计算该车的百公里油耗K;
在里程-油量图中,每个点用(里程,油量)坐标进行唯一表示,用Xi表示第i条数据的里程,Yi表示第i条数据的油量,其中Xi的单位为km,Yi的单位为升;
记Xn=Xi+m+1+100,Yn为里程Xn下的平均油量,则该车的百公里油耗
步骤2:根据盲区内各点的里程及百公里油耗K进行油量修正;
在需要进行油量补偿的区间段内,第i条数据的修正油量第i+1条数据的修正油量……第i+m-1条数据的修正油量
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410196768.0A CN103968910B (zh) | 2014-05-09 | 2014-05-09 | 一种提高机动车燃油量估算精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410196768.0A CN103968910B (zh) | 2014-05-09 | 2014-05-09 | 一种提高机动车燃油量估算精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103968910A CN103968910A (zh) | 2014-08-06 |
CN103968910B true CN103968910B (zh) | 2017-02-01 |
Family
ID=51238695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410196768.0A Active CN103968910B (zh) | 2014-05-09 | 2014-05-09 | 一种提高机动车燃油量估算精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103968910B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10365146B2 (en) * | 2014-12-26 | 2019-07-30 | Ford Global Technologies, Llc | Method and system for engine cooling system control |
US9719409B2 (en) * | 2014-12-26 | 2017-08-01 | Ford Global Technologies, Llc | Method and system for engine cooling system control |
CN104617914B (zh) * | 2015-02-11 | 2018-09-07 | 珠海格力电器股份有限公司 | 一种电器设备的信号滤波方法及系统 |
WO2018081945A1 (zh) | 2016-11-02 | 2018-05-11 | 深圳市汇顶科技股份有限公司 | 传感器本底噪声的检测方法及其装置 |
CN108656951A (zh) * | 2017-03-31 | 2018-10-16 | 华晨汽车集团控股有限公司 | 一种汽车燃油表的控制方法 |
CN110426101A (zh) * | 2019-05-20 | 2019-11-08 | 北京中位科技有限公司 | 一种油箱油量监测优化方法 |
CN111967950B (zh) * | 2020-10-13 | 2021-07-20 | 南京闻航汽车科技有限公司 | 与车辆交互的方法、计算设备和计算机可读存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5301113A (en) * | 1993-01-07 | 1994-04-05 | Ford Motor Company | Electronic system and method for calculating distance to empty for motorized vehicles |
CN101470017A (zh) * | 2007-12-29 | 2009-07-01 | 厦门雅迅网络股份有限公司 | 一种通过无线网络实时监测汽车油量及统计汽车油耗的方法 |
CN101865715A (zh) * | 2010-06-17 | 2010-10-20 | 天合导航通信技术有限公司 | 一种用于车载定位终端的油量信号同步补偿方法 |
CN202393445U (zh) * | 2011-12-22 | 2012-08-22 | 深圳市赛格导航科技股份有限公司 | 支持定位追踪的油量检测系统和支持定位跟踪的油量检测盒 |
CN102890858A (zh) * | 2012-09-13 | 2013-01-23 | 南京交通职业技术学院 | 客车能耗远程监测系统 |
CN103017854A (zh) * | 2011-09-28 | 2013-04-03 | 艾默生网络能源系统北美公司 | 油量计算方法、装置及监控系统 |
-
2014
- 2014-05-09 CN CN201410196768.0A patent/CN103968910B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5301113A (en) * | 1993-01-07 | 1994-04-05 | Ford Motor Company | Electronic system and method for calculating distance to empty for motorized vehicles |
CN101470017A (zh) * | 2007-12-29 | 2009-07-01 | 厦门雅迅网络股份有限公司 | 一种通过无线网络实时监测汽车油量及统计汽车油耗的方法 |
CN101865715A (zh) * | 2010-06-17 | 2010-10-20 | 天合导航通信技术有限公司 | 一种用于车载定位终端的油量信号同步补偿方法 |
CN103017854A (zh) * | 2011-09-28 | 2013-04-03 | 艾默生网络能源系统北美公司 | 油量计算方法、装置及监控系统 |
CN202393445U (zh) * | 2011-12-22 | 2012-08-22 | 深圳市赛格导航科技股份有限公司 | 支持定位追踪的油量检测系统和支持定位跟踪的油量检测盒 |
CN102890858A (zh) * | 2012-09-13 | 2013-01-23 | 南京交通职业技术学院 | 客车能耗远程监测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103968910A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103968910B (zh) | 一种提高机动车燃油量估算精度的方法 | |
CN102700551B (zh) | 一种车辆行驶过程中路面坡度的实时估计方法 | |
CN102627108B (zh) | 一种基于高频信息提取的整车质量估算方法 | |
CN104192145B (zh) | 一种车辆有限变速巡航方法 | |
CN107003430B (zh) | 路面状态预测方法和路面状态预测系统 | |
CN101539113B (zh) | 确定风力涡轮机塔架的谐振频率的装置和方法 | |
CN106926845B (zh) | 一种车辆状态参数的动态估计方法 | |
CN105223614B (zh) | 一种基于dwt_sta/lta的含噪信号p波初至峰度拾取方法 | |
RU2549598C2 (ru) | Система и способ оценки стиля вождения моторного транспортного средства | |
CN102608542A (zh) | 动力电池荷电状态估计方法 | |
CN103661395A (zh) | 动态道路坡度估计 | |
CN105651254B (zh) | 基于道路线形和频谱特征的道路坡度计算方法 | |
US11117590B2 (en) | Method and system for determining effective wind speeds for motor vehicles | |
CN104316716B (zh) | 一种用gps速度信息改进车载速度表的方法 | |
CN103661352A (zh) | 静态道路坡度估计 | |
CN102486400A (zh) | 车辆质量辨识方法和装置 | |
CN105015561B (zh) | 动车组的粘着控制系统 | |
CN113159435B (zh) | 一种新能源车辆剩余驾驶里程预测方法及系统 | |
CN103994748A (zh) | 一种采用飞行和风洞试验数据估计无人机配平迎角的方法 | |
CN103310637B (zh) | 基于车载数据的公交车速度获取装置及方法 | |
EP3155373A1 (de) | Verfahren und system zur anpassung eines navigationssystems | |
CN107064865A (zh) | 基于深度聚类的极坐标动态规划无源协同定位方法 | |
CN105627978B (zh) | 一种雨季滑坡变形监测方法及系统 | |
CN111038519B (zh) | 一种车载道路坡度实时估计方法 | |
WO2015189165A1 (de) | Verfahren und system zur anpassung eines navigationssystems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |