CN113879048B - 基于速度区间插值的间接式胎压监测标定方法 - Google Patents

基于速度区间插值的间接式胎压监测标定方法 Download PDF

Info

Publication number
CN113879048B
CN113879048B CN202111071594.1A CN202111071594A CN113879048B CN 113879048 B CN113879048 B CN 113879048B CN 202111071594 A CN202111071594 A CN 202111071594A CN 113879048 B CN113879048 B CN 113879048B
Authority
CN
China
Prior art keywords
speed
value
calibration
speed interval
domain analysis
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
CN202111071594.1A
Other languages
English (en)
Other versions
CN113879048A (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.)
Ruolun Automobile Technology Wuhan Co ltd
Original Assignee
Ruolun Automobile Technology Wuhan 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 Ruolun Automobile Technology Wuhan Co ltd filed Critical Ruolun Automobile Technology Wuhan Co ltd
Priority to CN202111071594.1A priority Critical patent/CN113879048B/zh
Publication of CN113879048A publication Critical patent/CN113879048A/zh
Application granted granted Critical
Publication of CN113879048B publication Critical patent/CN113879048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60C—VEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C23/00—Devices for measuring, signalling, controlling, or distributing tyre pressure or temperature, specially adapted for mounting on vehicles; Arrangement of tyre inflating devices on vehicles, e.g. of pumps or of tanks; Tyre cooling arrangements
    • B60C23/02—Signalling devices actuated by tyre pressure
    • B60C23/04—Signalling devices actuated by tyre pressure mounted on the wheel or tyre
    • B60C23/0408—Signalling devices actuated by tyre pressure mounted on the wheel or tyre transmitting the signals by non-mechanical means from the wheel or tyre to a vehicle body mounted receiver
    • B60C23/0474—Measurement control, e.g. setting measurement rate or calibrating of sensors; Further processing of measured values, e.g. filtering, compensating or slope monitoring
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明公开了一种基于速度区间插值的间接式胎压监测标定方法，将同一车辆的速度根据大小划分出不同的速度区间，并通过学习阶段的原始标定值拟合出各个速度区间之间的关系，从而根据车辆当前速度所属速度区间计算出可插入区间的标定值，并将计算的值与该速度区间的原始标定值作对比，从而判断轮胎是否出现异常。通过上述方式，本发明能够解决车辆的iTPMS(indirectTirePressureMonitoringSystem)系统在学习阶段中无法完成所有主要速度区间的学习时导致的iTPMS系统不能在未学习的速度区间提供应有的报警功能的问题，完善了当前iTPMS系统的监测报警功能，具有极大的市场前景。

Description

基于速度区间插值的间接式胎压监测标定方法
技术领域
本发明涉及轮胎胎压间接监测技术领域，特别是涉及一种基于速度区间插值的间接式胎压监测标定方法。
背景技术
随着科技水平的不断进步，各种车辆已经成为人们生活中必不可少的交通工具。近年来，车辆安全问题愈加严重，人们对于车辆安全性能的关注也越来越高。其中，以ABS/ESC为基础的主动安全概念受到越来越多汽车厂商以及消费者的关注。在汽车领域中，轮胎的安全问题一直占据及其重要的地位。统计表明，在高速行驶中，如果车速超过160公里/小时，前轮爆胎时，车内人员的生存几率几乎为零。而交通事故中爆胎原因所占的比例高达60％，这其中75％的爆胎都是由胎压不足造成，且轮胎欠压对于油耗和轮胎寿命都有影响。因此，对轮胎胎压的实时监测显得及其重要。
胎压监测技术主要包括直接式胎压监测技术dTPMS与间接式胎压监测技术iTPMS。目前，iTPMS对于胎压的分析主要有两部分，分别为时域分析TDA和频域分析FDA。其中，TDA是根据车轮之间的轮速差异对车轮的相对车轮滚动半径进行计算；FDA则是计算每个车轮的角速度频谱特征。在车辆自学习阶段完成之后，用实时监测到的TDA和FDA分析结果与学习阶段的值做比较，从而得出胎压的相对变化情况。
由于TDA和FDA的分析都与车速相关，所以学习阶段所存的值是分速度区间来存储的。在车辆自学习完成之后，实时计算的TDA和FDA分析结果也是与相同速度区间的标定值进行比较。在这种情况下iTPMS的监测报警功能只能在完成学习的速度区间中实现，这不符合实际应用时车辆对iTPMS功能的需求。在工程设计上需要根据速度区间划分标定界限，但对于iTPMS的用户来讲，这一划分规则是隐性的，客户在使用iTPMS进行自学习时，车辆极有可能在短时间内无法完成所有主要速度区间的学习，使得iTPMS不能在未学习的速度区间提供应有的监测报警功能。
因此，设计一种能够解决车辆的iTPMS系统在学习阶段无法完成所有主要速度区间的学习时导致的iTPMS系统不能在未学习的速度区间提供应有的报警功能的问题的基于速度区间插值的间接式胎压监测标定方法就很有必要。
发明内容
为了克服上述问题，本发明提供一种基于速度区间插值的间接式胎压监测标定方法，将同一车辆的速度根据大小划分出不同的速度区间，并通过学习阶段的原始标定值拟合出各个速度区间之间的关系。从而根据车辆当前速度所属速度区间计算出可插入区间的标定值，并将计算的标定值与该速度区间的原始标定值作对比，从而判断轮胎是否出现异常。这种方法能够实现实时采集的TDA及FDA跨速度区间的插值判断，解决了车辆的iTPMS系统在学习阶段无法完成所有主要速度区间的学习时导致的iTPMS系统不能在未学习的速度区间提供应有的报警功能的问题。
为实现上述的目的，本发明采用的技术方案是：
一种基于速度区间插值的间接式胎压监测标定方法，包括以下步骤：
S1、车辆进入自学习阶段，iTPMS系统采集时域分析和频域分析的分析结果以作为原始标定值；
S2、将各个原始标定值对应的速度按数值大小划分成不同速度区间，并对各个速度区间内的时域分析标定值与频域分析标定值分别进行归一化处理；
S3、将不同速度区间内的时域分析标定值与频域分析标定值进行拟合，以得出不同速度区间内的时域分析标定值及频域分析标定值之间的关系模型；
S4、判断行驶阶段车辆速度是否位于学习阶段的速度区间内，若是，将实时监测到的当前时域分析实时值与频域分析实时值直接与初始标定值比较后判断是否报警；若不是，则执行S5；
S5、判断当前时域分析实时值与频域分析实时值所属速度区间是否满足插值条件，若是，则将实时监测到的当前时域分析实时值与频域分析实时值与插值后得到的速度区间的初始标定值比较后判断是否报警。
进一步的，车辆自学习阶段的时间不小于20min，并在平整光滑的路面上进行，以保证后续实时采样数据能满足插值条件。
进一步的，在步骤S2中，将样本点划分为七个速度区间，并将各个速度区间内的时域分析标定值与频域分析标定值分别进行归一化处理：
Figure BDA0003260564600000031
式中，
Figure BDA0003260564600000032
是修正后的标定值，V是原始值，μ是所有原始值的平均值，σ是所有原始值的标准差。
进一步的，步骤S3包括以下步骤：
S31、将归一化后的时域分析标定值与频域分析标定值采用梯度下降的方法进行多项式拟合：
y＝θ0+θ1x+θ2x2+θ3x3+θ4x4+θ5x5+θ6x6+θ7x7
式中，x为每一个样本点对应的速度区间，y为每一个样本点对应的标定值，θ0至θ7为需要进行拟合的参数；
S32、采用随机数对参数θ0至θ7进行初始化，并将其整合成参数矩阵：
θ＝[θ0 θ1 θ2 θ3 θ4 θ5 θ6 θ7]
S33、将所有样本点的速度区间x整合成输入矩阵：
Figure BDA0003260564600000033
式中，m为所有样本点的个数；
S34、将所有样本点的标定值整合成输出矩阵：
Y＝[y1 … ym]
S35、将参数矩阵、输入矩阵以及输出矩阵带入拟合关系式，并采用梯度下降法求解θ：
Figure BDA0003260564600000041
将上式整合为：
Figure BDA0003260564600000042
S36、对S35中整合后的拟合公式进行误差分析，以得出最优化的拟合曲线。
进一步的，在步骤S36中，误差分析包括以下步骤：
S361、对S35中整合后的拟合公式的因变量进行误差预测，并定义损失函数：
Figure BDA0003260564600000043
式中，
Figure BDA0003260564600000044
为拟合值与实际值之间的误差；
S362、对参数矩阵θ进行梯度计算：
Figure BDA0003260564600000045
S363、对参数矩阵θ进行更新：
θ＝θ-αdθ
式中，α为学习速率，并取值0.1；
S364、重复步骤S363，直至
Figure BDA0003260564600000046
开始收敛，并将此时的参数矩阵θ记作/>
Figure BDA0003260564600000047
设此时的拟合曲线方程为：
Figure BDA0003260564600000048
S365、根据上述拟合曲线方程计算不同速度区间下标定值的差值，得出差值公式：
ΔVj＝f(j+1)-f(j)，j∈[1，6]。
进一步的，步骤S4包括以下步骤：
S41、ABS系统监测当前车辆行驶速度，并寻找其对应的速度区间；
S42、读出步骤S41的速度区间内是否存在原始标定值；
S43、若是，将实时监测到的当前时域分析实时值与频域分析实时值直接与原始标定值比较后判断是否报警；若不是，则执行S5。
进一步的，在步骤S5中，所述插值条件包括部分速度区间在学习阶段采集到满足标定个数的原始标定值、以及不是跨超过两个的速度区间进行插值。
进一步的，满足标定个数的原始标定值为单个速度区间的时域分析标定值个数大于25、以及频域分析标定值个数大于30，在此条件下可将实测值插入至其所述速度区间的相邻两个速度区间内进行比较。
进一步的，满足标定个数的原始标定值为两个速度区间的时域分析标定值个数大于10、以及频域分析标定值个数大于15，或者，两个速度区间的时域分析标定值个数之和大于30、以及频域分析标定值个数之和大于45，在此条件下可将实测值插入至其所述速度区间的相邻两个速度区间以及相隔一个速度区间的速度区间内进行比较。
进一步的，当任意一个速度区间的时域分析标定值个数大于25且频域分析标定值个数大于30时，不再使用插值后的值作为当前速度区间的标定值，而是使用当前速度区间自身样本点计算出的值作为标定值。
与现有技术相比，本发明的有益效果是：
1.本发明的基于速度区间插值的间接式胎压监测标定方法，通过将同一车辆的速度根据大小划分出不同的速度区间，并通过学习阶段的原始标定值拟合出各个速度区间之间的关系。从而根据车辆当前速度所属速度区间计算出可插入区间的标定值，并将计算的标定值与该速度区间的原始标定值作对比，从而判断轮胎是否出现异常。这种方法能够实现实时采集的TDA及FDA跨速度区间的插值判断，解决了车辆的iTPMS系统在学习阶段无法完成所有主要速度区间的学习时导致的iTPMS系统不能在未学习的速度区间提供应有的报警功能的问题，有效解决了当前iTPMS系统存在的技术短板。
2.本发明的基于速度区间插值的间接式胎压监测标定方法，通过对计算得到的拟合曲线进行误差分析，将拟合公式内的参数进行不断的更新，以得到最接近实际曲线的拟合曲线，保证插值计算公式的精准性。
附图说明
图1是本发明的基于速度区间插值的间接式胎压监测标定方法的流程框图；
图2是本发明的基于速度区间插值的间接式胎压监测标定方法的对TDA与FDA的标定值进行归一化处理前的分布示意图；
图3是本发明的基于速度区间插值的间接式胎压监测标定方法的对TDA与FDA的标定值进行归一化处理后的分布示意图；
图4是本发明的基于速度区间插值的间接式胎压监测标定方法的拟合曲线方程的曲线示意图；
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面结合附图和具体实施例对本发明进行详细描述。
在此，还需要说明的是，为了避免因不必要的细节而模糊了本发明，在附图中仅仅示出了与本发明的方案密切相关的结构和/或处理步骤，而省略了与本发明关系不大的其他细节。
另外，还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。
实施例
如图1所示，一种基于速度区间插值的间接式胎压监测标定方法100，包括以下步骤：
S1、车辆进入自学习阶段，iTPMS系统采集时域分析(Time Domain Analysis，下文简称TDA)和频域分析(Frequency Domain Analysis，下文简称FDA)的分析结果以作为原始标定值。
在本步骤中，iTPMS系统在车辆自学习阶段所采集的TDA与FDA的原始标定值是分速度区间来存储的。在车辆自学习完成之后，实时计算的TDA和FDA分析结果也是与相同速度区间的原始标定值进行比较。而本方法能够实现实时采集的TDA及FDA标定值的跨速度区间的插值判断。值得注意的是，该方法需要采集足够多的数据样本，以保证存在原始标定值的速度区间之间的间隔不会太大。基于此种情况，车辆自学习阶段是在平整光滑的路面上进行，且自学习阶段的时间不小于20min，以保证后续实时采样数据能满足插值条件。
S2、将各个原始标定值对应的速度按数值大小划分成不同速度区间，并对各个速度区间内的TDA与FDA分别进行归一化处理。
具体来讲，将样本点按对应速度大小进行归类，并以七个速度区间为例，将样本点划分成七个速度区间，并将各个速度区间内的TDA与FDA分别进行归一化处理：
Figure BDA0003260564600000071
式中，
Figure BDA0003260564600000072
是修正后的标定值，V是原始值，μ是所有原始值的平均值，σ是所有原始值的标准差。
随后，将TDA与FDA的标定值分别代入上式：
Figure BDA0003260564600000073
Figure BDA0003260564600000074
Figure BDA0003260564600000075
Figure BDA0003260564600000076
Figure BDA0003260564600000081
Figure BDA0003260564600000082
式中：
s表示速度区间，从1至7共有7个速度区间。
Figure BDA0003260564600000083
表示在速度区间s内TDA标定点的个数。
Figure BDA0003260564600000084
表示在速度区间s中的每一个TDA标定点的值，其中i是序号，取值在/>
Figure BDA0003260564600000085
内。
Figure BDA0003260564600000086
表示在速度区间s中的所有TDA标定值的平均值。
Figure BDA0003260564600000087
表示在速度区间s中的所有TDA标定值的标准差。
Figure BDA0003260564600000088
表示在速度区间s中经过归一化处理后的TDA标定值，并分别与每一个/>
Figure BDA0003260564600000089
相对应。
Figure BDA00032605646000000810
表示在速度区间s内FDA标定点的个数。
Figure BDA00032605646000000811
表示在速度区间s中的每一个FDA标定点的值，其中i是序号，取值在/>
Figure BDA00032605646000000812
内。
Figure BDA00032605646000000813
表示在速度区间s中的所有FDA标定值的平均值。
Figure BDA00032605646000000814
表示在速度区间si中的所有FDA标定值的标准差。
Figure BDA00032605646000000815
表示在速度区间si中经过归一化处理后的FDA标定值，分别与每一个/>
Figure BDA00032605646000000816
相对应。
如图2所示，图中分布点为对TDA与FDA的标定值进行归一化处理前，TDA与FDA的标定值分布示意图。从图中可以看出，此时TDA与FDA的标定值呈离散型分布，相同速度区间内的TDA与FDA的标定值之间间隔较大，整体样本点波动较大，不利于后续分析。
如图3所示，图中分布点为对TDA与FDA的标定值进行归一化处理后，TDA与FDA的标定值分布示意图。从图中可以看出，此时TDA与FDA的标定值分布较为紧凑，在相同速度区间内的TDA与FDA的标定值彼此重叠或接近，整体样本点平稳性好，后续分析结果更可靠。
S3、将不同速度区间内的TDA与FDA标定值进行拟合，以得出不同速度区间内的TDA及FDA标定值之间的关系模型。
具体来讲，步骤S3包括以下步骤：
S31、将归一化后的TDA与FDA标定值采用梯度下降的方法进行多项式拟合，将速度区间作为自变量x，将速度区间对应的标定值作为因变量y来进行特征升维，升维后的结果如下：
y＝θ0+θ1x+θ2x2+θ3x3+θ4x4+θ5x5+θ6x6+θ7x7
式中，x为每一个样本点对应的速度区间，y为每一个样本点对应的标定值，θ0至θ7为需要进行拟合的参数。该式用于反映速度区间与速度区间对应的标定值之间的关系模型，后续对该关系模型内的未知项进行分析，得出未知项后以将该关系模型补全。
S32、采用随机数对参数θ0至θ7进行初始化，得到一组数组，并将该数组整合成参数矩阵：
θ＝[θ0 θ1 θ2 θ3 θ4 θ5 θ6 θ7]
该参数矩阵可通过后续的计算进行求解，以得出参数θ的具体值。
S33、将所有样本点的速度区间x整合成输入矩阵：
Figure BDA0003260564600000091
式中，m为所有样本点的个数；
S34、将所有样本点的标定值整合成输出矩阵：
Y＝[y1 … ym]
S35、将参数矩阵、输入矩阵以及输出矩阵带入拟合关系式，并采用梯度下降法求解θ：
Figure BDA0003260564600000101
将上式整合为：
Figure BDA0003260564600000102
该式为速度区间与速度区间对应的标定值之间的拟合曲线，该拟合曲线由于采用近似算法进行拟合，其关系模型与实际模型存在一定的偏差。为了保证该拟合曲线的准确度，需要对该关系模型进行误差分析，以使拟合曲线最接近实际关系模型。
S36、对S35中整合后的拟合公式进行误差分析，以得出最优化的拟合曲线。
具体来讲，在步骤S36中，误差分析包括以下步骤：
S361、对S35中整合后的公式的因变量进行误差预测，并定义损失函数：
Figure BDA0003260564600000103
式中，
Figure BDA0003260564600000104
为拟合值与实际值之间的误差，以表示该拟合曲线所得到的速度区间的标定值与实际关系模型内的标定值的误差大小。
S362、对参数矩阵θ进行梯度计算，以θ来对
Figure BDA0003260564600000105
求偏导：
Figure BDA0003260564600000106
S363、对参数矩阵θ进行更新：
θ＝θ-αdθ
式中，α为学习速率，并取值0.1。
S364、重复步骤S363，以重复对θ进行更新，直至
Figure BDA0003260564600000107
开始收敛，并将此时的参数矩阵θ记作/>
Figure BDA0003260564600000111
由于当/>
Figure BDA0003260564600000112
开始收敛时，θ为极值点，此时求得的值最接近实际数值，拟合效果最好。设此时的拟合曲线方程为：
Figure BDA0003260564600000113
如图4所示，该图为上式关系模型得到的拟合曲线。图中，曲线平滑过渡，极点分布均匀，拟合效果好。因此，上述拟合曲线方程的拟合效果好。
S365、根据上述拟合曲线方程计算不同速度区间下标定值的差值，得出差值公式：
ΔVj＝f(j+1)-f(j)，j∈[1，6]。
根据上述公式，可以将未学习的速度区间的样本点插入至其它学习过的速度区间内，并与该速度区间内的原始标定值比较，以判断当前轮胎的胎压状态，实现了实时采集的TDA及FDA的跨速度区间的插值判断。
S4、判断行驶阶段车辆速度是否位于学习阶段的速度区间内，若是，将实时监测到的当前TDA与FDA直接与初始标定值比较后得出当前车胎的胎压状态，并判断是否报警；若不是，则执行S5。
具体来讲，步骤S4包括以下步骤：
S41、ABS系统监测当前车辆行驶速度，并寻找其对应的速度区间；
S42、读出步骤S41的速度区间内是否存在原始标定值；
S43、若是，将实时监测到的当前TDA与FDA直接与原始标定值比较后判断是否报警；若不是，则执行S5。
值得注意的是，由于该差值公式是根据原始标定值拟合出的关系曲线得到，速度区间相隔越近，其差值公式得到的关系模型越准确，速度区间相隔越远，其差值公式得到的关系模型误差越大。因而，当需要使用该公式时，需要满足相应的插值条件，才能保证该差值公式得出的结果最准确，保证胎压监测的精准度及可靠度。
S5、判断当前TDA与FDA所属速度区间是否满足插值条件，若是，则将实时监测到的当前TDA与FDA与插值后得到的速度区间的初始标定值比较后判断是否报警。
具体来讲，在步骤S5中，插值条件包括部分速度区间在学习阶段采集到满足标定个数的原始标定值、以及不是跨超过两个的速度区间进行插值。
具体分为下述三种情况：
经过试验论证得知，当满足标定个数的原始标定值为单个速度区间的TDA标定值个数大于25、以及FDA标定值个数大于30。TDA的标定值通常25秒左右会更新一个，所以在法规要求的20分钟自学习标定时间内最多会更新48个TDA的标定值。但是，考虑到车辆在标定阶段会出现起停以及刹车等无效时间，所以设定当有35个TDA标定值时自学习阶段完成。FDA的标定速度比TDA略快，在法规要求的20分钟自学习标定时间内最多会更新45个FDA标定值。当TDA标定值有25个，FDA标定值有30个时，此时虽然二者均未标定完成，但样本点数量足够多，根据统计学规律此时的标定值已经有相当高的可靠性，与标定完成时的标定值相比差别不大，此时可以使用该速度区间的标定值对相邻的速度区间进行插值。
在此条件下可将实测值插入至其速度区间的相邻两个速度区间内进行比较，此时根据差值公式得到的结果具备可靠性。
根据此种情况，设满足标定个数要求的速度区间为速度区间s，则满足该插值条件的速度区间的标定值计算如下：
当s＝1时，
Figure BDA0003260564600000121
当s＞1且s＜7时，
Figure BDA0003260564600000131
当s＝7时，
Figure BDA0003260564600000132
例如，将当前各个速度区间的标定值的个数用数组(N N Y N N N N)表示，其中Y表示该速度区间的标定值个数满足上述条件，N表示该速度区间的标定值个数不满足上述条件(下面所示的标记方法与该标记方法相同)。则能够进行采用差值公式进行插值预测的速度区间如数组(0 1 1 1 0 0 0)所示，其中1表示能采用差值公式进行插值预测，0表示不能采用差值公式进行插值预测(下面所示的标记方法与该标记方法相同)。
当满足标定个数的原始标定值为两个速度区间的TDA标定值个数大于10、以及FDA标定值个数大于15，或者，两个速度区间的TDA标定值个数之和大于30、以及FDA标定值个数之和大于45。在此条件下可将实测值插入至其速度区间的相邻速度区间以及相隔一个速度区间的速度区间内进行比较，此时根据差值公式得到的结果具备可靠性。这种情形下要求两个满足条件的速度区间至多相隔一个速度区间，即二者要么相邻，要么间隔一个速度区间。
在该情形中，以一个速度区间为相邻的速度区间的基准，这需要作为基准的速度区间的标定值足够可信。在当前情形中，是相邻两个速度区间互相做基准，所以此时单个速度区间的标定值个数可以少于情形一，但也要满足一定数量保证可靠性且两个速度区间总共的标定值个数要足够多。所以设置单个速度区间的TDA标定值个数要大于10，FDA标定值个数要大于15，两个速度区间总共的TDA标定个数要大于30，FDA标定值个数要大于45，在此条件下可以使用这两个速度区间的标定值对其他速度区间进行插值。
例如，当两个速度区间相邻，则该相邻的两个速度区间能采用差值公式进行插值预测。将当前各个速度区间的标定值的个数用数组(N N N Y Y N N)表示，则能够进行采用差值公式进行插值预测的速度区间如数组(0 0 0 1 1 0 0)所示。
根据此种情况，设满足标定个数要求的速度区间为速度区间s及速度区间s+1，则满足该插值条件的速度区间的标定值计算如下：
Figure BDA0003260564600000141
例如，当两个速度区间之间间隔一个速度区间，则该彼此之间间隔一个速度区间的两个速度区间能采用差值公式进行插值预测。将当前各个速度区间的标定值的个数用数组(N N N Y N Y N)表示，则能够进行采用差值公式进行插值预测的速度区间如数组(0 00 1 0 1 0)所示。
根据此种情况，设满足标定个数要求的速度区间为速度区间s-1及速度区间s+1，则满足该插值条件的速度区间的标定值计算如下：
Figure BDA0003260564600000151
值得注意的是，上述两种情况可以同时生效。
当任意一个速度区间的TDA标定值个数大于25且FDA标定值个数大于30时，不再使用插值后的值作为当前速度区间的标定值，而是使用当前速度区间自身样本点计算出的值作为标定值。这是由于插值算法是一种估计算法，且由于路况、车辆批次等因素会导致不同速度区间之间标定值的相关关系发生变化，这导致插值算法本身具有一定的系统误差，无法保证插值之后的分析结果一定准确，只能作为一种补充手段而无法替代真正的标定值，所以当某一速度区间的标定个数足够多时，应当首先使用其自身的标定值，而不是通过插值算法计算得到的标定值。
通过上述方法的各个步骤，本发明通过将同一车辆的速度根据大小划分出不同的速度区间，并通过学习阶段的原始标定值拟合出各个速度区间之间的关系。从而根据车辆当前速度所属速度区间计算出可插入区间的标定值，并将计算的标定值与该速度区间的原始标定值作对比，从而判断轮胎是否出现异常。这种方法能够实现实时采集的TDA及FDA跨速度区间的插值判断，解决了车辆的iTPMS系统在学习阶段无法完成所有主要速度区间的学习时导致的iTPMS系统不能在未学习的速度区间提供应有的报警功能的问题，有效解决了当前iTPMS系统存在的技术短板。
以上所述仅用以说明本发明的技术方案，而非对其进行限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (7)

1.一种基于速度区间插值的间接式胎压监测标定方法，其特征在于，包括以下步骤：
S1、车辆进入自学习阶段，iTPMS系统采集时域分析和频域分析的分析结果以作为原始标定值；
S2、将各个原始标定值对应的速度按数值大小划分成七个速度区间，并对各个速度区间内的时域分析标定值与频域分析标定值分别进行归一化处理：
Figure QLYQS_1
式中，
Figure QLYQS_2
是修正后的标定值，V是原始值，μ是所有原始值的平均值，σ是所有原始值的标准差；
S3、将不同速度区间内的时域分析标定值与频域分析标定值进行拟合，得到拟合曲线方程，用于反映速度区间与速度区间对应的标定值之间的关系模型；
包括以下步骤：
S31、将归一化后的时域分析标定值与频域分析标定值采用梯度下降的方法进行多项式拟合：
y＝θ0+θ1x+θ2x2+θ3x3+θ4x4+θ5x5+θ6x6+θ7x7
式中，x为每一个样本点对应的速度区间，y为每一个样本点对应的标定值，θ0至θ7为需要进行拟合的参数；
S32、采用随机数对参数θ0至θ7进行初始化，并将其整合成参数矩阵：
θ＝[θ0 θ1 θ2 θ3 θ4 θ5 θ6 θ7]
S33、将所有样本点的速度区间x整合成输入矩阵：
Figure QLYQS_3
式中，m为所有样本点的个数；
S34、将所有样本点的标定值整合成输出矩阵：
Y＝[y1 … ym]
S35、将参数矩阵、输入矩阵以及输出矩阵带入拟合关系式，并采用梯度下降法求解θ：
Figure QLYQS_4
将上式整合为：
Figure QLYQS_5
S36、对S35中整合后的拟合公式进行误差分析，以得出最优化的拟合曲线；
S4、判断行驶阶段车辆速度是否位于学习阶段的速度区间内，若是，将实时监测到的当前时域分析实时值与频域分析实时值直接与初始标定值比较后判断是否报警；若不是，则执行S5；
S5、判断当前时域分析实时值与频域分析实时值所属速度区间是否满足插值条件，所述插值条件包括部分速度区间在学习阶段采集到满足标定个数的原始标定值、以及不是跨超过两个的速度区间进行插值；若是，则将实时监测到的当前时域分析实时值与频域分析实时值与插值后得到的速度区间的初始标定值比较后判断是否报警。
2.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，车辆自学习阶段的时间不小于20min，并在平整光滑的路面上进行，以保证后续实时采样数据能满足插值条件。
3.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，在步骤S36中，误差分析包括以下步骤：
S361、对S35中整合后的拟合公式的因变量进行误差预测，并定义损失函数：
Figure QLYQS_6
式中，L为拟合值与实际值之间的误差；
S362、对参数矩阵θ进行梯度计算：
Figure QLYQS_7
S363、对参数矩阵θ进行更新：
θ＝θ-αdθ
式中，α为学习速率，并取值0.1；
S364、重复步骤S363，直至L开始收敛，并将此时的参数矩阵θ记作
Figure QLYQS_8
设此时的拟合曲线方程为：
Figure QLYQS_9
S365、根据上述拟合曲线方程计算不同速度区间下标定值的差值，得出差值公式：
ΔVj＝f(j+1)-f(j)，j∈[1，6]；
根据上述公式，将未学习的速度区间的样本点插入至其它学习过的速度区间内，并与该速度区间内的原始标定值比较，以判断当前轮胎的胎压状态，实现实时采集的时域分析及频域分析的跨速度区间的插值判断。
4.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，步骤S4包括以下步骤：
S41、ABS系统监测当前车辆行驶速度，并寻找其对应的速度区间；
S42、读出步骤S41的速度区间内是否存在原始标定值；
S43、若是，将实时监测到的当前时域分析实时值与频域分析实时值直接与原始标定值比较后判断是否报警；若不是，则执行S5。
5.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，满足标定个数的原始标定值为单个速度区间的时域分析标定值个数大于25、以及频域分析标定值个数大于30，在此条件下可将实测值插入至其所述速度区间的相邻两个速度区间内进行比较。
6.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，满足标定个数的原始标定值为两个速度区间的时域分析标定值个数大于10、以及频域分析标定值个数大于15，或者，两个速度区间的时域分析标定值个数之和大于30、以及频域分析标定值个数之和大于45，在此条件下可将实测值插入至其所述速度区间的相邻两个速度区间以及相隔一个速度区间的速度区间内进行比较。
7.根据权利要求1所述的基于速度区间插值的间接式胎压监测标定方法，其特征在于，当任意一个速度区间的时域分析标定值个数大于25且频域分析标定值个数大于30时，不再使用插值后的值作为当前速度区间的标定值，而是使用当前速度区间自身样本点计算出的值作为标定值。
CN202111071594.1A 2021-09-14 2021-09-14 基于速度区间插值的间接式胎压监测标定方法 Active CN113879048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111071594.1A CN113879048B (zh) 2021-09-14 2021-09-14 基于速度区间插值的间接式胎压监测标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111071594.1A CN113879048B (zh) 2021-09-14 2021-09-14 基于速度区间插值的间接式胎压监测标定方法

Publications (2)

Publication Number Publication Date
CN113879048A CN113879048A (zh) 2022-01-04
CN113879048B true CN113879048B (zh) 2023-06-27

Family

ID=79008889

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111071594.1A Active CN113879048B (zh) 2021-09-14 2021-09-14 基于速度区间插值的间接式胎压监测标定方法

Country Status (1)

Country Link
CN (1) CN113879048B (zh)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2125393A1 (en) * 2007-03-16 2009-12-02 Nira Dynamics AB Calibration in indirect tire pressure monitoring
JP5530395B2 (ja) * 2011-04-20 2014-06-25 富士重工業株式会社 運転評価装置
US20130151504A1 (en) * 2011-12-09 2013-06-13 Microsoft Corporation Query progress estimation
CN103832223A (zh) * 2012-11-20 2014-06-04 傅黎明 一种tpms技术缺陷改进方案
KR101373151B1 (ko) * 2012-12-18 2014-03-20 현대오트론 주식회사 차량용 타이어의 상태 탐지 장치 및 방법
US10254129B2 (en) * 2015-02-27 2019-04-09 Infineon Technologies Ag Low noise zero crossing detection for indirect tire pressure monitoring
KR101735728B1 (ko) * 2015-07-24 2017-05-15 현대오트론 주식회사 제로 크로싱을 이용한 타이어 압력 모니터링 방법 및 장치
CN109774389B (zh) * 2019-03-25 2021-08-06 宗培亮 一种基于频谱分析的间接式胎压监测方法
CN111016552B (zh) * 2019-12-26 2021-08-24 武汉理工大学 一种间接式胎压监测系统和方法
CN111224643B (zh) * 2020-01-16 2021-08-06 西安交通大学 快速收敛的样条自适应滤波方法
CN112345869A (zh) * 2020-11-25 2021-02-09 武汉光庭信息技术股份有限公司 一种汽车电子设备测试方法、系统、电子设备及存储介质

Also Published As

Publication number Publication date
CN113879048A (zh) 2022-01-04

Similar Documents

Publication Publication Date Title
CN111301426B (zh) 基于gru网络模型预测未来行驶过程能耗的方法
CN110562263B (zh) 一种基于多模型融合的轮毂电机驱动车辆车速估计方法
CN111806240B (zh) 一种基于中国工况的电动汽车续驶里程预测方法
JP5837341B2 (ja) 路面状態判定方法とその装置
CN108556682B (zh) 一种续驶里程预测方法、装置及设备
CN106956680B (zh) 一种电动汽车驾驶行为识别分析方法
CN112700156A (zh) 一种新能源汽车运行安全性能评价体系构建方法
Formentin et al. h-TPMS: a hybrid tire pressure monitoring system for road vehicles
CN113879048B (zh) 基于速度区间插值的间接式胎压监测标定方法
CN114161888B (zh) 双胎iTPMS胎压监测方法及其系统
EP3835161A1 (en) Vehicle mass estimation device
CN113505932A (zh) 一种基于大数据技术评估的动力电池容量算法
CN111483277A (zh) 基于机器学习的轮胎欠压识别方法、系统及存储介质
CN117057219A (zh) 轮胎受损预测方法、装置、设备、存储介质及车辆
Zhao et al. Indirect tire pressure monitoring method based on the fusion of time and frequency domain analysis
CN117171506A (zh) 一种基于imu的工程车辆运行状态监测方法及装置
CN115049996B (zh) 基于证据推理规则的双传感器目标检测融合方法和系统
CN116452008A (zh) 一种基于多项式建模的二手车风险预测方法及系统
CN115891520A (zh) 间接式轮胎压力监测方法、装置、车辆及存储介质
CN108407812A (zh) 轮胎更换后的轮径纠偏方法及系统
CN109345136A (zh) 基于我国机动车登记制度的车辆存活曲线模型优化方法
CN114707037A (zh) 新能源车辆节能驾驶推荐方法、装置、设备及存储介质
CN117755154B (zh) 基于新能源汽车实车工况的动力电池健康状态评估方法
CN114435378A (zh) 一种基于神经网络的纯电动汽车整车质量估算方法
CN115130060A (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