CN107117178A - 考虑换挡和道路坡度因素的车辆质量估计方法 - Google Patents
考虑换挡和道路坡度因素的车辆质量估计方法 Download PDFInfo
- Publication number
- CN107117178A CN107117178A CN201710369680.8A CN201710369680A CN107117178A CN 107117178 A CN107117178 A CN 107117178A CN 201710369680 A CN201710369680 A CN 201710369680A CN 107117178 A CN107117178 A CN 107117178A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- mtd
- 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.)
- Granted
Links
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
- 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/12—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 parameters of the vehicle itself, e.g. tyre models
- B60W40/13—Load or weight
-
- 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/0043—Signal treatments, identification of variables or parameters, parameter estimation or state estimation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Physics (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Control Of Transmission Device (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种考虑换挡和道路坡度因素的车辆质量估计方法,包括如下步骤:步骤1:采集数据:步骤11:利用数据采集装置采集车辆行驶状态数据;步骤12:结合车辆固有参数与车辆行驶状态数据,计算得到模型所需参数;步骤2:基于已获得的车辆行驶状态数据,建立车辆纵向动力学模型;步骤3:建立考虑换挡因素的车辆质量估计模型;步骤4:建立传动系机械效率估计模型,并与步骤3提出的考虑换挡因素的车辆质量估计模型共同构成同时考虑换挡和道路坡度因素的车辆质量估计模型;步骤5:确定车辆质量实时估计***的使用条件;步骤6:将采集到的车辆行驶状态数据和相关模型参数输入车辆质量实时估计***,估计实时车辆质量。
Description
技术领域
本发明涉及车辆质量估计方法,具体的为一种考虑换挡和道路坡度因素的车辆质量估计方法。
背景技术
车辆动力学模型参数的实时估计是车辆控制的基础,而车辆质量参数是车辆动力学模型中的重要参数,准确实时的车辆质量估计参数能够有效提高车辆主动安全***的性能,如电子稳定控制***(ESP)依赖车辆质量参数实时计算质心侧偏角和横摆角,从而实现对车辆横向滑移情况的控制。同时,有研究表明,根据质量参数的实时变化情况调整换挡规律能够为汽车节油约2.62%-4.86%。因此,准确估计车辆质量参数对于提高车辆安全性和燃油经济性均具有重要意义。
现有文献中,有关车辆质量参数估计的方法可以分为两类,一类是基于传感器的质量估计方法,这类方法需要在车内加装相应的传感器,占用车内空间的同时还提高了车辆的生产成本,难以满足实际应用需求;另一类是基于车辆CAN总线数据的车辆质量估计方法,该方法基于车辆纵向动力学模型,无需附加传感器,具有较好的应用便利性。然而,车辆实际行驶的道路不可避免存在一定的坡度,在车辆纵向动力学模型中,坡度和车辆质量存在一定的耦合关系,如果在质量估计过程中不考虑道路坡度的影响,必然会降低质量估计的精度;同时,又由于车辆在正常行驶过程中需要经常换挡,而换挡过程车辆的动力学状态无法用车辆纵向动力学描述,且换挡前后车辆纵向动力学模型的旋转质量换算系数会发生改变,从而导致换挡行为发生后车辆质量估计结果存在较大的偏差,严重影响车辆质量估计结果,进而降低车辆质量估计的准确性。
发明内容
有鉴于此,为了解决现有技术中存在的不足,本发明的目的在于提供一种考虑换挡和道路坡度因素的车辆质量估计方法。
为达到上述目的,本发明提供如下技术方案:
一种考虑换挡和道路坡度因素的车辆质量估计方法,包括如下步骤:
步骤1:采集数据
步骤11:利用数据采集装置采集车辆行驶状态数据;
步骤12:结合车辆固有参数与车辆行驶状态数据,计算得到模型所需参数;
步骤2:基于已获得的车辆行驶状态数据,建立车辆纵向动力学模型,获取不同档位下转矩与车速关系;
步骤3:基于车辆纵向动力学模型,建立考虑换挡因素的车辆质量估计模型;
步骤4:基于车辆纵向动力学模型,建立传动系机械效率估计模型,并与步骤3提出的考虑换挡因素的车辆质量估计模型共同构成同时考虑换挡和道路坡度因素的车辆质量估计模型;
步骤5:确定车辆质量实时估计***的使用条件;
步骤6:将采集到的车辆行驶状态数据和相关模型参数输入车辆质量实时估计***,估计实时车辆质量。
进一步,所述步骤11中,车辆行驶状态数据包括发动机转矩Tq、车速v、发动机转速n;
所述步骤12中,车辆固有参数包括轮胎滚动半径r、主减速器传动比i0、道路滚动阻力系数f、传动系机械效率η、车辆空气阻力系数Cd、车辆正向迎风面积A;
模型所需参数包括车辆加速度a和变速器传动比ig;
其中,加速度a可由速度v对时间差分获得,第k时刻的加速度表示为:
T为数据采集装置的采样周期;
变速器传动比ig的计算方法为:
进一步,车辆纵向动力学模型为:
Ft=Ff+Fi+Faero+Fj
其中,Ft为车辆驱动力,且
Ff为滚动阻力,且Ff=mgf;
Fi为坡度阻力,且Fi=mgi;
Faero为空气阻力,且
Fj为加速阻力,且Fj=ma;
其中,Tq表示发动机转矩;ig表示变速器传动比;i0表示主减速器传动比;η表示传动系效率;r表示车轮半径;m表示车辆质量;f表示滚动阻力系数;Cd表示车辆的空气阻力系数;A表示迎风面积;v表示车辆速度;i为道路坡度;
代入后,得到不同档位下转矩与车速关系:
进一步,所述步骤3包括如下步骤:
步骤31:推导车辆纵向动力学模型的最小二乘形式;
步骤32:建立带多遗忘因子的加权最小二乘递推估计模型;
步骤33:建立带多遗忘因子的加权最小二乘递推质量估计模型。
进一步,所述步骤31中,将车辆纵向动力学模型转化成如下最小二乘形式:
其中,a表示加速度;δ表示旋转质量换算系数。
进一步,所述步骤32中,设***的输入输出关系可以描述成如下的最小二乘形式:
z(k)=hT(k)θ+n(k)
其中,z(k)是***的输出,h(k)是可观测数据向量,n(k)为白噪声,θ为待估计参数;
当模型中存在两个待估计参数时,定义准则函数为:
其中Λ(i)为加权函数,λ1和λ2分别为两个模型待估参数θ1和θ2对应的遗忘因子;
利用序列{z(k)}和{h(k)},极小化准则函数,即对θ求导,可求得参数θ的最小二乘估计值第k时刻的参数估计值可以表示为:
将上述估计结果转化为递推形式,得到带多遗忘因子的加权最小二乘递推估计模型如下:
其中:
进一步,所述步骤33中,将步骤31中的车辆纵向动力学模型最小二乘形式应用于步骤32中的最小二乘递推估计模型,则有:
将上式带入最小二乘可得带多遗忘因子的加权最小二乘质量估计模型的递推形式为:
其中,λ1和λ2分别为两个待估计参数m和δm对应的遗忘因子,取值范围为[0,1),且λ1>λ2;
Λ(k)为数据的加权因子,为保证非换挡条件下的车辆行驶状态数据能够占据更高的数据权重,设当前时刻为t(k),tgear为最近一次换挡时间,换挡前的加权因子恒为1,换挡后的加权因子可表示为:
进一步,所述步骤4包括如下步骤:
步骤41:建立传动系机械效率估计模型;
将基于车辆纵向动力学模型推导出相应的最小二乘形式应用于带多遗忘因子的加权最小二乘估计模型中,其中:
基于车辆纵向动力学模型推导出相应的最小二乘形式为:
设:
则可得到传动系机械效率的估计模型如下:
其中:
步骤42:基于传动系机械效率估计模型和考虑换挡因素的车辆质量估计模型,共同构建同时考虑换挡和道路坡度因素的车辆质量估计模型。
进一步,确定车辆质量实时估计***的使用条件的方法为
当同时满足以下条件时,质量估计***启动:
1)在车辆启动并行驶距离超过300m;
2)车辆速度大于20km/h;
3)车辆加速度为正值;
4)油门开度大于0;
5)车辆未处于制动状态;
6)车辆未处于换挡状态;
7)方向盘转角不大于30°;
8)当前行驶路段坡度小于2°;
9)车辆档位不为空挡;
当满足以下任意条件时,质量估计***暂时关闭:
1)车辆处于制动状态;
2)车辆方向盘转角小于90°;
3)车辆档位为空挡;
4)车辆处于停车状态。
进一步,所述步骤6中,估计实时车辆质量主要包括以下三个阶段:
1)启动阶段:
将质量估计模型待估计参数的估计初值均设置为0,并设置带多遗忘因子的加权最小二乘递推质量估计模型中信息增益矩阵的初值:
P(0)=a2I
其中a为充分大的正数;
随后基于当前时刻的车辆行驶状态数据,开始进行第一次质量估计,获得m和δm的第一次估计值;
2)运行阶段
将当前时刻更新的车辆行驶状态数据和上一时刻的m和δm的估计值带入加权最小二乘递推质量估计模型中,进行新一轮质量估计,并重复此步骤,获取实时的质量估计值;
3)停止阶段
当车辆质量的估计值收敛在合理范围之内,停止进行递推估计,但依旧接收车辆行驶状态数据,以便意外情况出现时,可以再次进行质量估计。
本发明的有益效果在于:
本发明考虑换挡和道路坡度因素的车辆质量估计方法,使车辆在平直和带有坡度的路面上,均能够获得较为理想的质量估计结果;同时考虑换挡因素带来的传动比突变和换挡前后车辆旋转质量换算系数的变化,采用带有多遗忘因子的加权最小二乘递推估计方法,实现对车辆质量的实时估计;即在正常行驶和车辆换挡时均可实时、准确地估计车辆质量,在降低换挡因素对于质量估计造成的影响的同时也提高了车辆质量估计的准确性和鲁棒性。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明考虑换挡和道路坡度因素的车辆质量估计方法的流程图;
图2为考虑换挡和道路坡度因素的车辆质量估计模型;
图3为车辆纵向动力学模型受力示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
本实施例考虑换挡和道路坡度因素的车辆质量估计方法,包括如下步骤:
步骤1:采集数据
步骤11:利用OpenXC数据采集装置采集车辆行驶状态数据
具体的,本实施例利用OpenXC数据采集装置采集车辆行驶状态数据,且车辆行驶状态数据包括发动机转矩Tq、车速v、发动机转速n、方向盘转角φw、油门开度、刹车信号、档位信息等,并基于数据采集装置内置传感器测算得出的实时道路坡度i;通过蓝牙接收来自车辆CAN总线上的实时车辆行驶状态数据并存储在移动数据处理终端上。
步骤12:结合车辆固有参数与车辆行驶状态数据,计算得到模型所需参数
具体的,本实施例的车辆固有参数包括轮胎滚动半径r、主减速器传动比i0、道路滚动阻力系数f、传动系机械效率η、车辆空气阻力系数Cd、车辆正向迎风面积A等;模型所需参数包括车辆加速度a和变速器传动比ig。
其中,加速度a可由速度v对时间差分获得,第k时刻的加速度表示为:
T为OpenXC数据采集装置的采样周期;
变速器传动比ig的计算方法为:
步骤2:基于已获得的车辆行驶状态数据,建立车辆纵向动力学模型,获取不同档位下转矩与车速关系。具体的,车辆纵向动力学模型为:
Ft=Ff+Fi+Faero+Fj
其中,Ft为车辆驱动力,且
Ff为滚动阻力,且Ff=mgf;
Fi为坡度阻力,且Fi=mgi;
Faero为空气阻力,且
Fj为加速阻力,且Fj=ma;
其中,Tq表示发动机转矩;ig表示变速器传动比;i0表示主减速器传动比;η表示传动系效率;r表示车轮半径;m表示车辆质量;f表示滚动阻力系数;Cd表示车辆的空气阻力系数;A表示迎风面积;v表示车辆速度;i为道路坡度;
代入后,得到不同档位下转矩与车速关系:
步骤3:基于车辆纵向动力学模型,建立考虑换挡因素的车辆质量估计模型
步骤31:推导车辆纵向动力学模型的最小二乘形式
本实施例的目标是在考虑换挡因素的条件实时估计车辆质量,而在车辆换挡前后,车辆的旋转质量换算系δ数会发生较为明显的变化,影响车辆质量估计。仅基于车辆CAN总线的数据无法实时计算车辆的旋转质量换算系数,为保证换挡条件下仍能够获得较为准确的车辆质量估计值,需要将车辆的旋转质量换算系数δ和车辆质量m一同作为车辆质量估计模型的待估计参数。因此,将车辆纵向动力学模型转化成如下最小二乘形式:
其中,a表示加速度;δ表示旋转质量换算系数。
步骤32:建立带多遗忘因子的加权最小二乘递推估计模型
在车辆启动后,车辆质量m几乎不变,是一个慢变量;而旋转质量换算系数δ会随着换挡等因素发生快速变化,是一个快变量。为保证这种快变量和慢变量共同存在的估计模型中慢变量的估计值能够快速且准确的收敛,同时快变量的估计结果能够具备较好的跟踪性能,需要为不同的待估计参数设置不同的遗忘因子,以提升模型估计结果的准确性。
又由于车辆在换挡过程中的运动状态无法用车辆纵向动力学描述,换挡过程中产生的车辆行驶状态数据会增大车辆质量估计结果的误差,为了抑制这一过程中产生的误差,需采用加权方法,减弱换挡过程中车辆行驶状态数据的权重,使质量估计方法在换挡过程中依然能够输出准确的质量估计结果。
设***的输入输出关系可以描述成如下的最小二乘形式:
z(k)=hT(k)θ+n(k)
其中,z(k)是***的输出,h(k)是可观测数据向量,n(k)为白噪声,θ为待估计参数;
当模型中存在两个待估计参数时,定义准则函数为:
其中Λ(i)为加权函数,λ1和λ2分别为两个模型待估参数θ1和θ2对应的遗忘因子;
利用序列{z(k)}和{h(k)},极小化准则函数,即对θ求导,可求得参数θ的最小二乘估计值第k时刻的参数估计值可以表示为:
为了保证车辆质量估计***能够实时更新估计结果,需将上述估计结果转化为递推形式,得到带多遗忘因子的加权最小二乘递推估计模型如下:
其中:
步骤33:建立带多遗忘因子的加权最小二乘递推质量估计模型
将步骤31中的车辆纵向动力学模型最小二乘形式应用于步骤32中的最小二乘递推估计模型,则有:
将上式带入最小二乘可得带多遗忘因子的加权最小二乘质量估计模型的递推形式为:
其中,λ1和λ2分别为两个待估计参数m和δm对应的遗忘因子,取值范围为[0,1),考虑到δm为快变量,为保证δm的估计结果具有较好的跟踪性能,应有λ1>λ2,默认情况下λ1取为0.95,λ2取为0.5。
Λ(k)为数据的加权因子,为保证非换挡条件下的车辆行驶状态数据能够占据更高的数据权重,设当前时刻为t(k),tgear为最近一次换挡时间,换挡前的加权因子恒为1,换挡后的加权因子可表示为:
步骤4:基于车辆纵向动力学模型,建立传动系机械效率估计模型,并与步骤3提出的考虑换挡因素的车辆质量估计模型共同构成同时考虑换挡和道路坡度因素的车辆质量估计模型;
步骤41:建立传动系机械效率估计模型
在带有坡度的路面上进行车辆质量估计时,仅仅依靠为质量估计模型提供准确的道路坡度信息,仍然无法准确进行质量估计。这主要是由于随着道路坡度的变化,车辆传动系效率η也在不断变化,而该参数无法通过车辆CAN总线或外部传感器获得,因此必须建立能够实时获取传动系机械效率的估计模型,并将其应用于带多遗忘因子的加权最小二乘递推质量估计模型中,从而确保车辆在存在一定坡度的路面上估计质量时的准确度。
首先基于车辆纵向动力学模型推导出相应的最小二乘形式
通过上式可以看出,最小二乘形式将传动系机械效率η和车辆质量m视为待估计参数,其中车辆质量在行驶过程中几乎不变,属于慢变量;传动系机械效率随着车辆传动系的工作状态的变化而快速变化,属于快变量。因此,需要为车辆质量赋予一个较大的遗忘因子,以保证估计的连续和稳定性;为传动系机械效率赋予一个较小的遗忘因子,以保证估计结果具有良好的跟踪性能。由于在估计过程中有可能出现换挡,为了抑制换挡时车辆行驶状态数据对估计结果造成的不良影响,也需要加入和考虑换挡因素质量估计模型一样的加权因子Λ(k),从而保证估计结果的稳定性和连续性。
将车辆纵向动力学模型转换得到的最小二乘形式应用于带多遗忘因子的加权最小二乘估计模型中,设:
则可得到传动系机械效率的估计模型如下:
其中:
步骤42:基于传动系机械效率估计模型和考虑换挡因素的车辆质量估计模型,共同构建同时考虑换挡和道路坡度因素的车辆质量估计模型。
将步骤3得到的带多遗忘因子的加权最小二乘递推质量估计模型与车辆传动系效率模型结合在一起,便构成了考虑换挡和道路坡度因素的车辆质量估计模型,如图2所示。
在车辆质量估计开始之前,根据车辆的实际情况,预先设定传动系效率初始值和旋转质量换算系数的初始值。质量估计开始时,OpenXC数据采集装置提供的车辆行驶状态数据同时输入带多遗忘因子的加权最小二乘递推质量估计模型和车辆传动系效率模型之中。车辆行驶状态数据与传动系效率初始值作为模型输入共同完成带多遗忘因子的加权最小二乘递推质量估计模型的第一次估计,同时车辆传动系效率估计模型也将基于车辆行驶状态数据和旋转质量换算系数的初值,完成第一次估计。随后将第一次获得的旋转质量换算系数和传动系机械效率的估计值分别替代其预设的初值,作为模型输入,分别进行两个模型的第二次估计。然后继续将上一时刻模型输出的旋转质量换算系数和传动系机械效率估计值替代原先的模型输入,与最新的车辆行驶状态数据一同进行下一轮递推估计,重复这一过程,直至获得质量估计曲线收敛。
考虑换挡和道路坡度因素的车辆质量估计模型克服了换挡过程带来的传动比突变问题,并通过实时估计车辆传动系效率的办法保证了车辆在坡路上估计车辆质量的精度,在驾驶员存在换挡行为以及道路存在坡度的条件下依然能够获得比较准确的质量估计结果。
步骤5:确定车辆质量实时估计***的使用条件
当同时满足以下条件时,质量估计***启动:
1)在车辆启动并行驶距离超过300m;
2)车辆速度大于20km/h;
3)车辆加速度为正值;
4)油门开度大于0;
5)车辆未处于制动状态;
6)车辆未处于换挡状态;
7)方向盘转角不大于30°;
8)当前行驶路段坡度小于2°;
9)车辆档位不为空挡;
当满足以下任意条件时,质量估计***暂时关闭:
1)车辆处于制动状态;
2)车辆方向盘转角小于90°;
3)车辆档位为空挡;
4)车辆处于停车状态。
本实施例的质量估计方法一般用于估计当前行程下的车辆质量,当车辆熄火后,车辆质量估计***会将历史数据清除,并于下次行程重新进行质量估计。
步骤6:将采集到的车辆行驶状态数据和相关模型参数输入车辆质量实时估计***,估计实时车辆质量。本实施例估计实时车辆质量主要包括以下三个阶段:
1)启动阶段:
将质量估计模型待估计参数的估计初值均设置为0,并设置带多遗忘因子的加权最小二乘递推质量估计模型中信息增益矩阵的初值:
P(0)=a2I
其中a为充分大的正数,默认值取为1000000;
随后基于当前时刻的车辆行驶状态数据,开始进行第一次质量估计,获得m和δm的第一次估计值;
2)运行阶段
将当前时刻更新的车辆行驶状态数据和上一时刻的m和δm的估计值带入加权最小二乘递推质量估计模型中,进行新一轮质量估计,并重复此步骤,获取实时的质量估计值;
3)停止阶段
当车辆质量的估计值收敛在合理范围之内,停止进行递推估计,但依旧接收车辆行驶状态数据,以便意外情况出现时,可以再次进行质量估计。
当车辆熄火后,停止进行递推估计,质量估计***中的历史数据也将清零,质量估计***停止运行。
本实施例考虑换挡和道路坡度因素的车辆质量估计方法,使车辆在平直和带有坡度的路面上,均能够获得较为理想的质量估计结果;同时考虑换挡因素带来的传动比突变和换挡前后车辆旋转质量换算系数的变化,采用带有多遗忘因子的加权最小二乘递推估计方法,实现对车辆质量的实时估计;即在正常行驶和车辆换挡时均可实时、准确地估计车辆质量,在降低换挡因素对于质量估计造成的影响的同时也提高了车辆质量估计的准确性和鲁棒性。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (10)
1.一种考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:包括如下步骤:
步骤1:采集数据
步骤11:利用数据采集装置采集车辆行驶状态数据;
步骤12:结合车辆固有参数与车辆行驶状态数据,计算得到模型所需参数;
步骤2:基于已获得的车辆行驶状态数据,建立车辆纵向动力学模型,获取不同档位下转矩与车速关系;
步骤3:基于车辆纵向动力学模型,建立考虑换挡因素的车辆质量估计模型;
步骤4:基于车辆纵向动力学模型,建立传动系机械效率估计模型,并与步骤3提出的考虑换挡因素的车辆质量估计模型共同构成同时考虑换挡和道路坡度因素的车辆质量估计模型;
步骤5:确定车辆质量实时估计***的使用条件;
步骤6:将采集到的车辆行驶状态数据和相关模型参数输入车辆质量实时估计***,估计实时车辆质量。
2.根据权利要求1所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:
所述步骤11中,车辆行驶状态数据包括发动机转矩Tq、车速v、发动机转速n;
所述步骤12中,车辆固有参数包括轮胎滚动半径r、主减速器传动比i0、道路滚动阻力系数f、传动系机械效率η、车辆空气阻力系数Cd、车辆正向迎风面积A;
模型所需参数包括车辆加速度a和变速器传动比ig;
其中,加速度a可由速度v对时间差分获得,第k时刻的加速度表示为:
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mrow>
T为数据采集装置的采样周期;
变速器传动比ig的计算方法为:
<mrow>
<msub>
<mi>i</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.377</mn>
<mi>r</mi>
<mi>n</mi>
</mrow>
<mrow>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mi>v</mi>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:车辆纵向动力学模型为:
Ft=Ff+Fi+Faero+Fj
其中,Ft为车辆驱动力,且
Ff为滚动阻力,且Ff=mgf;
Fi为坡度阻力,且Fi=mgi;
Faero为空气阻力,且
Fj为加速阻力,且Fj=ma;
其中,Tq表示发动机转矩;ig表示变速器传动比;i0表示主减速器传动比;η表示传动系效率;r表示车轮半径;m表示车辆质量;f表示滚动阻力系数;Cd表示车辆的空气阻力系数;A表示迎风面积;v表示车辆速度;i为道路坡度;
代入后,得到不同档位下转矩与车速关系:
<mrow>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<msub>
<mi>i</mi>
<mi>g</mi>
</msub>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mo>=</mo>
<mi>m</mi>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>m</mi>
<mi>g</mi>
<mi>i</mi>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<msup>
<mi>Av</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
<mo>+</mo>
<mi>m</mi>
<mi>a</mi>
<mo>.</mo>
</mrow>
4.根据权利要求3所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤3包括如下步骤:
步骤31:推导车辆纵向动力学模型的最小二乘形式;
步骤32:建立带多遗忘因子的加权最小二乘递推估计模型;
步骤33:建立带多遗忘因子的加权最小二乘递推质量估计模型。
5.根据权利要求4所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤31中,将车辆纵向动力学模型转化成如下最小二乘形式:
<mrow>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<msub>
<mi>i</mi>
<mi>g</mi>
</msub>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<msup>
<mi>Av</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
</mrow>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>m</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>m</mi>
<mi>&delta;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,a表示加速度;δ表示旋转质量换算系数。
6.根据权利要求5所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤32中,设***的输入输出关系可以描述成如下的最小二乘形式:
z(k)=hT(k)θ+n(k)
其中,z(k)是***的输出,h(k)是可观测数据向量,n(k)为白噪声,θ为待估计参数;
当模型中存在两个待估计参数时,定义准则函数为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&lambda;</mi>
<mn>1</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&lambda;</mi>
<mn>2</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中Λ(i)为加权函数,λ1和λ2分别为两个模型待估参数θ1和θ2对应的遗忘因子;
利用序列{z(k)}和{h(k)},极小化准则函数,即对θ求导,可求得参数θ的最小二乘估计值第k时刻的参数估计值可以表示为:
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msubsup>
<mi>&lambda;</mi>
<mn>1</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&lambda;</mi>
<mn>1</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<msubsup>
<mi>&lambda;</mi>
<mn>2</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>&Lambda;</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&lambda;</mi>
<mn>2</mn>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
将上述估计结果转化为递推形式,得到带多遗忘因子的加权最小二乘递推估计模型如下:
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<msup>
<mi>h</mi>
<mi>T</mi>
</msup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
2
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<msubsup>
<mi>h</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<msubsup>
<mi>h</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>h</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>h</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>.</mo>
</mrow>
7.根据权利要求6所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤33中,将步骤31中的车辆纵向动力学模型最小二乘形式应用于步骤32中的最小二乘递推估计模型,则有:
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<msub>
<mi>i</mi>
<mi>g</mi>
</msub>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mi>A</mi>
<mi>v</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
将上式带入最小二乘可得带多遗忘因子的加权最小二乘质量估计模型的递推形式为:
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<msub>
<mi>i</mi>
<mi>g</mi>
</msub>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mi>A</mi>
<mi>v</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>(</mo>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
其中,λ1和λ2分别为两个待估计参数m和δm对应的遗忘因子,取值范围为[0,1),且λ1>λ2;
Λ(k)为数据的加权因子,为保证非换挡条件下的车辆行驶状态数据能够占据更高的数据权重,设当前时刻为t(k),tgear为最近一次换挡时间,换挡前的加权因子恒为1,换挡后的加权因子可表示为:
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mn>0.95</mn>
<mfrac>
<mrow>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mi>T</mi>
</mfrac>
</msup>
<mo>.</mo>
</mrow>
8.根据权利要求7所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤4包括如下步骤:
步骤41:建立传动系机械效率估计模型;
将基于车辆纵向动力学模型推导出相应的最小二乘形式应用于带多遗忘因子的加权最小二乘估计模型中,其中:
基于车辆纵向动力学模型推导出相应的最小二乘形式为:
<mrow>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<msup>
<mi>Av</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mi>v</mi>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>+</mo>
<mi>&delta;</mi>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&eta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>m</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
设:
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mi>A</mi>
<mi>v</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
则可得到传动系机械效率的估计模型如下:
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mi>A</mi>
<mi>v</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mn>21.15</mn>
</mfrac>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>0.377</mn>
<msup>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>(</mo>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>0.377</mn>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<mi>&Lambda;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>g</mi>
<mi>f</mi>
<mo>+</mo>
<mi>g</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>;</mo>
</mrow>
步骤42:基于传动系机械效率估计模型和考虑换挡因素的车辆质量估计模型,共同构建同时考虑换挡和道路坡度因素的车辆质量估计模型。
9.根据权利要求1所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:确定车辆质量实时估计***的使用条件的方法为
当同时满足以下条件时,质量估计***启动:
1)在车辆启动并行驶距离超过300m;
2)车辆速度大于20km/h;
3)车辆加速度为正值;
4)油门开度大于0;
5)车辆未处于制动状态;
6)车辆未处于换挡状态;
7)方向盘转角不大于30°;
8)当前行驶路段坡度小于2°;
9)车辆档位不为空挡;
当满足以下任意条件时,质量估计***暂时关闭:
1)车辆处于制动状态;
2)车辆方向盘转角小于90°;
3)车辆档位为空挡;
4)车辆处于停车状态。
10.根据权利要求1所述的考虑换挡和道路坡度因素的车辆质量估计方法,其特征在于:所述步骤6中,估计实时车辆质量主要包括以下三个阶段:
1)启动阶段:
将质量估计模型待估计参数的估计初值均设置为0,并设置带多遗忘因子的加权最小二乘递推质量估计模型中信息增益矩阵的初值:
P(0)=a2I
其中a为充分大的正数;
随后基于当前时刻的车辆行驶状态数据,开始进行第一次质量估计,获得m和δm的第一次估计值;
2)运行阶段
将当前时刻更新的车辆行驶状态数据和上一时刻的m和δm的估计值带入加权最小二乘递推质量估计模型中,进行新一轮质量估计,并重复此步骤,获取实时的质量估计值;
3)停止阶段
当车辆质量的估计值收敛在合理范围之内,停止进行递推估计,但依旧接收车辆行驶状态数据,以便意外情况出现时,可以再次进行质量估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710369680.8A CN107117178B (zh) | 2017-05-23 | 2017-05-23 | 考虑换挡和道路坡度因素的车辆质量估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710369680.8A CN107117178B (zh) | 2017-05-23 | 2017-05-23 | 考虑换挡和道路坡度因素的车辆质量估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107117178A true CN107117178A (zh) | 2017-09-01 |
CN107117178B CN107117178B (zh) | 2019-06-25 |
Family
ID=59729928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710369680.8A Active CN107117178B (zh) | 2017-05-23 | 2017-05-23 | 考虑换挡和道路坡度因素的车辆质量估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107117178B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107826124A (zh) * | 2017-11-02 | 2018-03-23 | 潍柴动力股份有限公司 | 一种基于发动机制动的整车下坡提示方法及*** |
CN108189843A (zh) * | 2017-12-28 | 2018-06-22 | 天津清智科技有限公司 | 一种汽车重量的估算方法 |
CN110727994A (zh) * | 2019-10-28 | 2020-01-24 | 吉林大学 | 参数解耦的电动汽车质量与坡度估计方法 |
CN111311782A (zh) * | 2020-02-12 | 2020-06-19 | 北京经纬恒润科技有限公司 | 载重估计方法及装置 |
CN111497859A (zh) * | 2020-06-30 | 2020-08-07 | 北京主线科技有限公司 | 一种结合重量参数辨识的车辆纵向控制方法 |
CN111579037A (zh) * | 2020-04-29 | 2020-08-25 | 北理新源(佛山)信息科技有限公司 | 一种检测车辆超载的方法和*** |
CN112379348A (zh) * | 2020-11-17 | 2021-02-19 | 安徽四创电子股份有限公司 | 基于卡尔曼滤波的多模型估计目标机动识别方法 |
CN112590803A (zh) * | 2020-12-16 | 2021-04-02 | 北理慧动(常熟)车辆科技有限公司 | 一种单轴并联混合动力商用车整车质量在线估计方法 |
CN112926140A (zh) * | 2021-03-25 | 2021-06-08 | 交通运输部公路科学研究所 | 一种基于车路协同与tbox的货运车辆质量估计方法 |
CN113392518A (zh) * | 2021-06-08 | 2021-09-14 | 阿波罗智联(北京)科技有限公司 | 估计车辆重量的方法和装置 |
WO2021223334A1 (zh) * | 2020-05-06 | 2021-11-11 | 北京理工大学 | 基于mmrls和sh-stf的车辆质量与道路坡度迭代型联合估计方法 |
CN114564791A (zh) * | 2022-02-24 | 2022-05-31 | 广东工业大学 | 一种基于车辆运行数据的公交车总重量测量方法 |
CN114954494A (zh) * | 2022-06-14 | 2022-08-30 | 广西玉柴机器股份有限公司 | 一种重型商用车辆载重快速估计方法 |
CN114987510A (zh) * | 2022-06-17 | 2022-09-02 | 东风悦享科技有限公司 | 一种自动驾驶车辆质量参数在线估计方法及装置 |
WO2024138355A1 (zh) * | 2022-12-27 | 2024-07-04 | 采埃孚商用车***(青岛)有限公司 | 一种用于估算车辆质量的方法、装置和电子设备以及计算机可读介质和车辆 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102486400A (zh) * | 2010-12-06 | 2012-06-06 | 罗伯特·博世有限公司 | 车辆质量辨识方法和装置 |
CN102951158A (zh) * | 2012-11-02 | 2013-03-06 | 浙江吉利汽车研究院有限公司杭州分公司 | 车辆质量估算方法 |
CN103264669A (zh) * | 2013-05-31 | 2013-08-28 | 吉林大学 | 一种基于can信息和功能原理的重型车质量实时辨识方法 |
EP2862772A2 (de) * | 2013-10-18 | 2015-04-22 | MAN Truck & Bus AG | Verfahren zur Überprüfung eines Beladungszustandes eines Aufliegers oder Anhängers eines Nutzfahrzeuges |
-
2017
- 2017-05-23 CN CN201710369680.8A patent/CN107117178B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102486400A (zh) * | 2010-12-06 | 2012-06-06 | 罗伯特·博世有限公司 | 车辆质量辨识方法和装置 |
CN102951158A (zh) * | 2012-11-02 | 2013-03-06 | 浙江吉利汽车研究院有限公司杭州分公司 | 车辆质量估算方法 |
CN103264669A (zh) * | 2013-05-31 | 2013-08-28 | 吉林大学 | 一种基于can信息和功能原理的重型车质量实时辨识方法 |
EP2862772A2 (de) * | 2013-10-18 | 2015-04-22 | MAN Truck & Bus AG | Verfahren zur Überprüfung eines Beladungszustandes eines Aufliegers oder Anhängers eines Nutzfahrzeuges |
Non-Patent Citations (2)
Title |
---|
冯源等: "基于分段递推最小二乘估计的汽车质量辨识试验", 《同济大学学报(自然科学版)》 * |
许嘉等: "参数估计算法在汽车动态称重中的应用", 《自动化技术与应用》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107826124A (zh) * | 2017-11-02 | 2018-03-23 | 潍柴动力股份有限公司 | 一种基于发动机制动的整车下坡提示方法及*** |
CN108189843A (zh) * | 2017-12-28 | 2018-06-22 | 天津清智科技有限公司 | 一种汽车重量的估算方法 |
CN110727994A (zh) * | 2019-10-28 | 2020-01-24 | 吉林大学 | 参数解耦的电动汽车质量与坡度估计方法 |
CN111311782A (zh) * | 2020-02-12 | 2020-06-19 | 北京经纬恒润科技有限公司 | 载重估计方法及装置 |
CN111579037B (zh) * | 2020-04-29 | 2021-06-04 | 北理新源(佛山)信息科技有限公司 | 一种检测车辆超载的方法和*** |
CN111579037A (zh) * | 2020-04-29 | 2020-08-25 | 北理新源(佛山)信息科技有限公司 | 一种检测车辆超载的方法和*** |
WO2021223334A1 (zh) * | 2020-05-06 | 2021-11-11 | 北京理工大学 | 基于mmrls和sh-stf的车辆质量与道路坡度迭代型联合估计方法 |
CN111497859A (zh) * | 2020-06-30 | 2020-08-07 | 北京主线科技有限公司 | 一种结合重量参数辨识的车辆纵向控制方法 |
CN112379348A (zh) * | 2020-11-17 | 2021-02-19 | 安徽四创电子股份有限公司 | 基于卡尔曼滤波的多模型估计目标机动识别方法 |
CN112379348B (zh) * | 2020-11-17 | 2024-04-09 | 安徽四创电子股份有限公司 | 基于卡尔曼滤波的多模型估计目标机动识别方法 |
CN112590803A (zh) * | 2020-12-16 | 2021-04-02 | 北理慧动(常熟)车辆科技有限公司 | 一种单轴并联混合动力商用车整车质量在线估计方法 |
CN112590803B (zh) * | 2020-12-16 | 2022-02-22 | 北理慧动(常熟)车辆科技有限公司 | 一种单轴并联混合动力商用车整车质量在线估计方法 |
CN112926140A (zh) * | 2021-03-25 | 2021-06-08 | 交通运输部公路科学研究所 | 一种基于车路协同与tbox的货运车辆质量估计方法 |
CN112926140B (zh) * | 2021-03-25 | 2023-07-28 | 交通运输部公路科学研究所 | 一种基于车路协同与tbox的货运车辆质量估计方法 |
CN113392518A (zh) * | 2021-06-08 | 2021-09-14 | 阿波罗智联(北京)科技有限公司 | 估计车辆重量的方法和装置 |
WO2022257310A1 (zh) * | 2021-06-08 | 2022-12-15 | 阿波罗智联(北京)科技有限公司 | 估计车辆重量的方法和装置 |
CN113392518B (zh) * | 2021-06-08 | 2024-05-17 | 阿波罗智联(北京)科技有限公司 | 估计车辆重量的方法和装置 |
CN114564791B (zh) * | 2022-02-24 | 2022-10-28 | 广东工业大学 | 一种基于车辆运行数据的公交车总重量测量方法 |
CN114564791A (zh) * | 2022-02-24 | 2022-05-31 | 广东工业大学 | 一种基于车辆运行数据的公交车总重量测量方法 |
CN114954494A (zh) * | 2022-06-14 | 2022-08-30 | 广西玉柴机器股份有限公司 | 一种重型商用车辆载重快速估计方法 |
CN114954494B (zh) * | 2022-06-14 | 2024-03-26 | 广西玉柴机器股份有限公司 | 一种重型商用车辆载重快速估计方法 |
CN114987510A (zh) * | 2022-06-17 | 2022-09-02 | 东风悦享科技有限公司 | 一种自动驾驶车辆质量参数在线估计方法及装置 |
WO2024138355A1 (zh) * | 2022-12-27 | 2024-07-04 | 采埃孚商用车***(青岛)有限公司 | 一种用于估算车辆质量的方法、装置和电子设备以及计算机可读介质和车辆 |
Also Published As
Publication number | Publication date |
---|---|
CN107117178B (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107117178A (zh) | 考虑换挡和道路坡度因素的车辆质量估计方法 | |
CN106740870B (zh) | 一种考虑换挡因素的车辆质量估计方法 | |
CN107247824A (zh) | 考虑刹车和转弯影响的汽车质量‑道路坡度联合估计方法 | |
CN103264669B (zh) | 一种基于can信息和功能原理的重型车质量实时辨识方法 | |
CN105644565B (zh) | 一种混合动力汽车负载的测量方法 | |
US8768536B2 (en) | Method for determining the driving resistance of a vehicle | |
CN106840097B (zh) | 一种基于自适应扩展卡尔曼滤波的道路坡度估计方法 | |
CN102398599B (zh) | 车辆控制***和车辆控制方法 | |
CN108437998B (zh) | 基于纵向动力学的纯电动汽车坡度识别方法 | |
US11987151B2 (en) | Control system and method for controlling an electric motor | |
CN111605559B (zh) | 一种整车质量估算方法、扭矩控制方法及装置 | |
CN103209858B (zh) | 电驱动车辆 | |
WO2013123832A1 (zh) | 一种基于车辆运行工况的发动机自适应***及节油方法 | |
CN102649433A (zh) | 用于提高燃料经济性指标计算的道路坡度估计的方法 | |
CN108068802A (zh) | 一种车辆蠕行控制方法及利用其的自动泊车方法 | |
CN102358288A (zh) | 一种车辆acc驱动工况路面峰值附着系数识别方法 | |
CN109760682B (zh) | 一种纯电动车辆爬坡扭矩估算方法和控制方法及其*** | |
CN105539449A (zh) | 一种制动工况下的路面附着系数实时估算方法 | |
CN104976337A (zh) | 基于离合器传递力矩估计的车辆起步过程优化控制方法 | |
CN111252072A (zh) | 针对纯电动汽车实现坡道检测处理的方法 | |
US20100204895A1 (en) | Method and Device for Controlling the Traction Slip of the Driven Wheels of a Vehicle With the Engine Torque as a Set Parameter | |
JP6217723B2 (ja) | 駆動力制御装置 | |
CN115140070A (zh) | 一种基于amt换挡控制的整车质量估算方法 | |
JP2012192834A (ja) | 車両制御装置 | |
CN114987509A (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230301 Address after: 213399 room 5025, building B, 218 Hongkou Road, Kunlun Street, Liyang City, Changzhou City, Jiangsu Province Patentee after: Liyang Smart City Research Institute of Chongqing University Address before: 400044 No. 174, positive street, Shapingba District, Chongqing Patentee before: Chongqing University |
|
TR01 | Transfer of patent right |