CN115447587A - Data processing method, device, non-volatile storage medium and computer equipment - Google Patents
Data processing method, device, non-volatile storage medium and computer equipment Download PDFInfo
- Publication number
- CN115447587A CN115447587A CN202211254981.3A CN202211254981A CN115447587A CN 115447587 A CN115447587 A CN 115447587A CN 202211254981 A CN202211254981 A CN 202211254981A CN 115447587 A CN115447587 A CN 115447587A
- Authority
- CN
- China
- Prior art keywords
- state information
- vehicle
- sampling points
- moment
- calculating
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 76
- 230000001133 acceleration Effects 0.000 claims abstract description 57
- 238000005070 sampling Methods 0.000 claims description 164
- 230000015654 memory Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 20
- 230000009466 transformation Effects 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims 3
- 239000013598 vector Substances 0.000 description 37
- 238000004364 calculation method Methods 0.000 description 22
- 230000008569 process Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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/02—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 ambient conditions
- B60W40/06—Road conditions
- B60W40/076—Slope angle of the road
-
- 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
-
- 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
- B60W2050/0052—Filtering, filters
-
- 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
- B60W2420/00—Indexing codes relating to the type of sensors based on the principle of their operation
- B60W2420/40—Photo, light or radio wave sensitive means, e.g. infrared sensors
- B60W2420/408—Radar; Laser, e.g. lidar
-
- 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
-
- 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/10—Longitudinal speed
-
- 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/10—Longitudinal speed
- B60W2520/105—Longitudinal acceleration
-
- 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
- B60W2552/00—Input parameters relating to infrastructure
- B60W2552/15—Road slope, i.e. the inclination of a road segment in the longitudinal direction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Physics (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种数据处理方法、装置、非易失性存储介质和计算机设备。其中,该方法包括:获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。本发明解决了现有技术中坡度估算方法精度低的技术问题。
The invention discloses a data processing method, device, nonvolatile storage medium and computer equipment. Wherein, the method includes: acquiring state information of the vehicle, wherein the state information includes first state information of the vehicle at a first moment and second state information of the vehicle at a second moment, and any one of the state information includes: Velocity, longitudinal acceleration, slope, driving force and wheel angular acceleration, the direction of longitudinal acceleration is from the rear of the vehicle to the front; according to the first state information and the second state information, the unscented Kalman filter method is used to calculate the vehicle at the The target slope at the second moment, where the unscented Kalman filter method uses the kinematics formula and dynamics formula of the vehicle to estimate the target slope. The invention solves the technical problem of low precision of the slope estimation method in the prior art.
Description
技术领域technical field
本发明涉及车辆辅助驾驶领域,具体而言,涉及一种数据处理方法、装置、非易失性存储介质和计算机设备。The present invention relates to the field of vehicle assisted driving, in particular to a data processing method, device, non-volatile storage medium and computer equipment.
背景技术Background technique
坡度信号是汽车功能软件开发的重要信号之一,无论是传统燃油车还是新能源车辆,节油、节能等软件算法的开发离不开该信号;对于越野车或者SUV等车辆,坡度信号影响全地形软件或者驾驶模式软件的状态切换;对于自动驾驶车辆,坡度信号影响车辆行驶安全算法模块,控制算法模块等。简言之,坡度信号是非常重要的基础信号之一,该信号的精确与否直接影响上层功能软件的结果,最终影响整车的表现。现有技术中对坡度的估算精度低,估算范围小。The slope signal is one of the important signals for the development of automobile function software. Whether it is a traditional fuel vehicle or a new energy vehicle, the development of software algorithms such as fuel saving and energy saving is inseparable from this signal; for vehicles such as off-road vehicles or SUVs, the slope signal affects the entire State switching of terrain software or driving mode software; for self-driving vehicles, slope signals affect vehicle driving safety algorithm modules, control algorithm modules, etc. In short, the slope signal is one of the very important basic signals. The accuracy of the signal directly affects the results of the upper-level functional software, and ultimately affects the performance of the vehicle. The estimation accuracy of the slope in the prior art is low, and the estimation range is small.
针对上述的问题,目前尚未提出有效的解决方案。For the above problems, no effective solution has been proposed yet.
发明内容Contents of the invention
本发明实施例提供了一种数据处理方法、装置、非易失性存储介质和计算机设备,以至少解决现有技术中坡度估算方法精度低的技术问题。Embodiments of the present invention provide a data processing method, device, non-volatile storage medium and computer equipment, so as to at least solve the technical problem of low accuracy of the slope estimation method in the prior art.
根据本发明实施例的一个方面,提供了一种数据处理方法,包括:获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。According to an aspect of an embodiment of the present invention, there is provided a data processing method, including: acquiring state information of the vehicle, wherein the state information includes the first state information of the vehicle at the first moment and the first state information of the vehicle at the second moment Two state information, any one of the state information includes: speed, longitudinal acceleration, slope, driving force and wheel angular acceleration, the direction of longitudinal acceleration is from the rear to the front of the vehicle; according to the first state information and the second state information , using the unscented Kalman filter method to calculate the target slope of the vehicle at the second moment, wherein the unscented Kalman filter method uses the vehicle's kinematics formula and dynamics formula to estimate the target slope.
可选地,根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,包括:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,其中,状态信息包括第三状态信息;根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,其中,状态信息包括第四状态信息;根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益;根据第二状态信息、第三状态信息、第四状态信息和卡尔曼增益,计算得到目标坡度。Optionally, according to the first state information and the second state information, the unscented Kalman filter method is used to calculate the target slope of the vehicle at the second moment, including: according to the first state information and the kinematic formula, calculating the slope of the vehicle at the second moment The third state information at the second moment, wherein the state information includes the third state information; according to the third state information and the dynamic formula, calculate the fourth state information of the vehicle at the second moment, wherein the state information includes the fourth state information; calculate the Kalman gain at the second moment according to the third state information and the fourth state information; calculate the target slope according to the second state information, the third state information, the fourth state information and the Kalman gain.
可选地,根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,包括:对第一状态信息进行无迹变换,生成多个第一采样点及多个第一采样点的加权值;根据运动学公式、多个第一采样点和多个第一采样点的加权值进行计算,得到车辆在第二时刻下的多个第二采样点和多个第二采样点的加权值;根据多个第二采样点和多个第二采样点的加权值,计算得到第三状态信息。Optionally, according to the first state information and the kinematic formula, calculating the third state information of the vehicle at the second moment includes: performing unscented transformation on the first state information to generate multiple first sampling points and multiple first sampling points A weighted value of a sampling point; calculated according to the kinematic formula, multiple first sampling points and the weighted values of multiple first sampling points, to obtain multiple second sampling points and multiple second sampling points of the vehicle at the second moment The weighted value of the sampling point; the third state information is obtained by calculating according to the multiple second sampling points and the weighted values of the multiple second sampling points.
可选地,根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,包括:对第三状态信息进行无迹变换,生成多个第三采样点及多个第三采样点的加权值;根据动力学公式、多个第三采样点及多个第三采样点的加权值,计算多个第四采样点及多个第四采样点的加权值;根据多个第四采样点及多个第四采样点的加权值,计算得到第四状态信息。Optionally, according to the third state information and the dynamic formula, calculating the fourth state information of the vehicle at the second moment includes: performing unscented transformation on the third state information to generate multiple third sampling points and multiple first The weighted value of three sampling points; according to the dynamic formula, the weighted value of multiple third sampling points and multiple third sampling points, calculate the weighted value of multiple fourth sampling points and multiple fourth sampling points; according to multiple The fourth sampling point and weighted values of the plurality of fourth sampling points are calculated to obtain fourth state information.
可选地,根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益,包括:根据多个第三采样点和多个第四采样点,计算多个采样点和多个第四采样点之间的交叉协方差矩阵;根据交叉协方差矩阵,计算得到卡尔曼增益。Optionally, calculating the Kalman gain at the second moment according to the third state information and the fourth state information includes: calculating a plurality of sampling points and a plurality of fourth sampling points according to a plurality of third sampling points and a plurality of fourth sampling points The cross-covariance matrix between the fourth sampling points; according to the cross-covariance matrix, the Kalman gain is calculated.
可选地,获取车辆的第二状态信息,包括:在第二时刻下,接收传感器测得的物理量,其中,传感器测得的物理量包括车辆的纵向加速度、驱动力、速度和车轮角加速度;根据传感器测得的物理量和动力学公式,计算得到第二状态信息包括的第二坡度。Optionally, acquiring the second state information of the vehicle includes: at the second moment, receiving the physical quantity measured by the sensor, wherein the physical quantity measured by the sensor includes the longitudinal acceleration, driving force, speed and wheel angular acceleration of the vehicle; according to The physical quantity measured by the sensor and the dynamic formula are calculated to obtain the second slope included in the second state information.
可选地,接收传感器测得的车轮角加速度,包括:接收传感器测得的车辆的车轮角速度和传感器测量车轮角速度的采样周期;根据车轮角速度和采样周期,计算角加速度。Optionally, receiving the wheel angular acceleration measured by the sensor includes: receiving the wheel angular velocity measured by the sensor and the sampling period of the wheel angular velocity measured by the sensor; and calculating the angular acceleration according to the wheel angular velocity and the sampling period.
根据本发明实施例的另一方面,还提供了一种数据处理装置,包括:获取模块,用于获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;计算模块,用于根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。According to another aspect of the embodiments of the present invention, there is also provided a data processing device, including: an acquisition module, configured to acquire state information of the vehicle, wherein the state information includes the first state information of the vehicle at the first moment and the vehicle In the second state information at the second moment, any one of the state information includes: speed, longitudinal acceleration, slope, driving force and wheel angular acceleration, the direction of longitudinal acceleration is from the rear to the front of the vehicle; the calculation module, It is used to calculate the target slope of the vehicle at the second moment by using the unscented Kalman filter method according to the first state information and the second state information, wherein the unscented Kalman filter method uses the kinematics formula and the dynamics formula of the vehicle Estimate the target slope.
根据本发明实施例的又一方面,还提供了一种非易失性存储介质,非易失性存储介质包括存储的程序,其中,在程序运行时控制非易失性存储介质所在设备执行上述中任意一项数据处理方法。According to yet another aspect of the embodiments of the present invention, a non-volatile storage medium is also provided, and the non-volatile storage medium includes a stored program, wherein, when the program is running, the device where the non-volatile storage medium is located is controlled to execute the above-mentioned Any one of the data processing methods.
根据本发明实施例的再一方面,还提供了一种计算机设备,计算机设备包括处理器,处理器用于运行程序,其中,程序运行时执行上述中任意一项数据处理方法。According to still another aspect of the embodiments of the present invention, there is also provided a computer device, the computer device includes a processor, and the processor is used to run a program, wherein, when the program is running, any one of the above data processing methods is executed.
在本发明实施例中,通过获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计,达到了将车辆状态信息数据进行融合计算坡度的目的,从而实现了扩大坡度计算范围且提高坡度计算精度的技术效果,进而解决了现有技术中坡度估算方法精度低的技术问题。In the embodiment of the present invention, by acquiring the state information of the vehicle, the state information includes the first state information of the vehicle at the first moment and the second state information of the vehicle at the second moment, any one of the state information Including: speed, longitudinal acceleration, slope, driving force and wheel angular acceleration. The direction of longitudinal acceleration is from the rear of the vehicle to the front of the vehicle; according to the first state information and the second state information, the unscented Kalman filter method is used to calculate the vehicle The target slope at the second moment, where the unscented Kalman filter method uses the kinematics formula and dynamics formula of the vehicle to estimate the target slope, and achieves the purpose of merging vehicle state information data to calculate the slope, thereby realizing The technical effect of expanding the slope calculation range and improving the slope calculation accuracy further solves the technical problem of low precision of the slope estimation method in the prior art.
附图说明Description of drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The accompanying drawings described here are used to provide a further understanding of the present invention and constitute a part of the application. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations to the present invention. In the attached picture:
图1示出了一种用于实现数据处理方法的计算机终端的硬件结构框图;Fig. 1 shows a kind of hardware structure block diagram for realizing the computer terminal of data processing method;
图2是根据本发明实施例提供的数据处理方法的流程示意图;Fig. 2 is a schematic flowchart of a data processing method provided according to an embodiment of the present invention;
图3是根据本发明实施例提供的数据处理装置的结构框图。Fig. 3 is a structural block diagram of a data processing device provided according to an embodiment of the present invention.
具体实施方式detailed description
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to enable those skilled in the art to better understand the solutions of the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only It is an embodiment of a part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts shall fall within the protection scope of the present invention.
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first" and "second" in the description and claims of the present invention and the above drawings are used to distinguish similar objects, but not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used are interchangeable under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having", as well as any variations thereof, are intended to cover a non-exclusive inclusion, for example, a process, method, system, product or device comprising a sequence of steps or elements is not necessarily limited to the expressly listed instead, may include other steps or elements not explicitly listed or inherent to the process, method, product or apparatus.
根据本发明实施例,提供了一种数据处理的方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。According to an embodiment of the present invention, an embodiment of a data processing method is provided. It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and, although A logical order is shown in the flowcharts, but in some cases the steps shown or described may be performed in an order different from that shown or described herein.
本申请实施例一所提供的方法实施例可以在移动终端、计算机终端或者类似的运算装置中执行。图1示出了一种用于实现数据处理方法的计算机终端的硬件结构框图。如图1所示,计算机终端10可以包括一个或多个(图中采用102a、102b,……,102n来示出)处理器(处理器可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)、用于存储数据的存储器104。除此以外,还可以包括:显示器、输入/输出接口(I/O接口)、通用串行总线(USB)端口(可以作为BUS总线的端口中的一个端口被包括)、网络接口、电源和/或相机。本领域普通技术人员可以理解,图1所示的结构仅为示意,其并不对上述电子装置的结构造成限定。例如,计算机终端10还可包括比图1中所示更多或者更少的组件,或者具有与图1所示不同的配置。The method embodiment provided in Embodiment 1 of the present application may be executed in a mobile terminal, a computer terminal, or a similar computing device. Fig. 1 shows a block diagram of the hardware structure of a computer terminal for realizing the data processing method. As shown in Figure 1, the computer terminal 10 may include one or more (shown by 102a, 102b, ..., 102n in the figure) processors (processors may include but not limited to microprocessors MCU or programmable logic devices Processing device such as FPGA),
应当注意到的是上述一个或多个处理器和/或其他数据处理电路在本文中通常可以被称为“数据处理电路”。该数据处理电路可以全部或部分的体现为软件、硬件、固件或其他任意组合。此外,数据处理电路可为单个独立的处理模块,或全部或部分的结合到计算机终端10中的其他元件中的任意一个内。如本申请实施例中所涉及到的,该数据处理电路作为一种处理器控制(例如与接口连接的可变电阻终端路径的选择)。It should be noted that the one or more processors and/or other data processing circuits described above may generally be referred to herein as "data processing circuits". The data processing circuit may be implemented in whole or in part as software, hardware, firmware or other arbitrary combinations. In addition, the data processing circuit can be a single independent processing module, or be fully or partially integrated into any of the other components in the computer terminal 10 . As mentioned in the embodiment of the present application, the data processing circuit is used as a processor control (for example, the selection of the terminal path of the variable resistor connected to the interface).
存储器104可用于存储应用软件的软件程序以及模块,如本发明实施例中的数据处理方法对应的程序指令/数据存储装置,处理器通过运行存储在存储器104内的软件程序以及模块,从而执行各种功能应用以及数据处理,即实现上述的应用程序的数据处理方法。存储器104可包括高速随机存储器,还可包括非易失性存储器,如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中,存储器104可进一步包括相对于处理器远程设置的存储器,这些远程存储器可以通过网络连接至计算机终端10。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。The
显示器可以例如触摸屏式的液晶显示器(LCD),该液晶显示器可使得用户能够与计算机终端10的用户界面进行交互。The display can be, for example, a touch-screen liquid crystal display (LCD), which enables a user to interact with the user interface of the computer terminal 10 .
图2是根据本发明实施例提供的数据处理方法的流程示意图,如图2所示,该方法包括如下步骤:Fig. 2 is a schematic flow chart of a data processing method provided according to an embodiment of the present invention. As shown in Fig. 2, the method includes the following steps:
步骤S202,获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头。Step S202, acquire the state information of the vehicle, wherein the state information includes the first state information of the vehicle at the first moment and the second state information of the vehicle at the second moment, any one of the state information includes: speed, longitudinal Acceleration, slope, driving force and wheel angular acceleration, the direction of longitudinal acceleration is from the rear of the vehicle to the front of the vehicle.
在车辆行驶过程中,设置在车辆上的传感器可以测量车辆行驶时的状态参数,并将测量得到的状态参数进行数据处理后得到车辆行驶时的状态信息,其中,状态信息可以包括速度、纵向加速度、车辆所处的坡度、车辆驱动力和车轮角加速度。需要说明的是,车辆所处的坡度不是可以直接测量得到的状态参数,需要由其他状态参数通过推理计算进行估算。状态信息相较于状态参数来说,更方便参与无迹卡尔曼滤波计算,在计算时状态信息可以使用向量进行表示。During the running of the vehicle, the sensors installed on the vehicle can measure the state parameters when the vehicle is running, and process the measured state parameters to obtain the state information when the vehicle is running. The state information can include speed, longitudinal acceleration , the slope of the vehicle, the driving force of the vehicle and the angular acceleration of the wheels. It should be noted that the slope of the vehicle is not a state parameter that can be directly measured, and needs to be estimated by other state parameters through inferential calculation. Compared with state parameters, state information is more convenient to participate in unscented Kalman filter calculation, and state information can be represented by vector during calculation.
估算车辆当前所处位置的坡度值可以是一个动态循环的过程,对于两次相邻的估算过程来说,上一次的状态信息估算结果在下一次估算过程中将起到很大作用。因此,在本次估算过程中,可以获取两个时刻下的车辆状态信息,分别为在第一时刻下车辆的第一状态信息和在第二时刻下车辆的第二状态信息,第一时刻是前一个估算过程开始的时刻,车辆的第一状态信息是前一个估算过程得到的估算结果,第二时刻是本次估算过程开始的时刻,车辆的第二状态信息是在第二时刻时测量得到的状态参数经过数据处理后得到的状态信息。Estimating the slope value of the vehicle's current location can be a dynamic cycle process. For two adjacent estimation processes, the state information estimation result of the last time will play a significant role in the next estimation process. Therefore, in this estimation process, the vehicle state information at two moments can be obtained, which are the first state information of the vehicle at the first moment and the second state information of the vehicle at the second moment. The first moment is At the moment when the previous estimation process starts, the first state information of the vehicle is the estimation result obtained in the previous estimation process; the second moment is the moment when this estimation process starts, and the second state information of the vehicle is obtained from the measurement at the second moment The status information of the status parameters obtained after data processing.
步骤S204,根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。Step S204, according to the first state information and the second state information, the unscented Kalman filter method is used to calculate the target slope of the vehicle at the second moment, wherein the unscented Kalman filter method uses the kinematics formula and dynamics of the vehicle The formula estimates the target slope.
本步骤中,在获取到第一状态信息和第二状态信息之后,可以采用无迹卡尔曼滤波方法,根据合适的运动学公式和动力学公式对第一状态信息和第二状态信息进行处理,估算得出车辆在第二时刻下的目标坡度。需要说明的是,由于在使用无迹卡尔曼滤波方法对第一状态信息和第二状态信息进行数据处理的过程中,状态信息的表示方法可以为使用向量进行表示,而且状态信息中包括车辆所处的坡度,所以实际计算时,可以通过向量之间的运算得到一个向量,得到的向量中存在表示目标坡度的数字,即第一状态信息的向量表示和第二状态信息的向量表示得到一个包含目标坡度的状态信息的向量表示,这个得到的状态信息中表示所处坡度的位置处的数字即为目标坡度的值。In this step, after the first state information and the second state information are acquired, the unscented Kalman filter method may be used to process the first state information and the second state information according to appropriate kinematics formulas and dynamics formulas, Estimate the target slope of the vehicle at the second moment. It should be noted that during the data processing process of the first state information and the second state information using the unscented Kalman filter method, the representation method of the state information can be expressed by using a vector, and the state information includes Therefore, in the actual calculation, a vector can be obtained through the operation between vectors, and there is a number representing the target slope in the obtained vector, that is, the vector representation of the first state information and the vector representation of the second state information obtain a vector representation containing The vector representation of the status information of the target slope, the number at the position indicating the slope in the obtained status information is the value of the target slope.
通过上述步骤,可以达到了将车辆在不同时刻下的状态信息数据进行融合计算坡度的目的,从而实现了扩大坡度计算范围且提高坡度计算精度的技术效果,进而解决了现有技术中坡度估算方法精度低的技术问题。Through the above steps, the purpose of merging the state information data of the vehicle at different times to calculate the slope can be achieved, thereby achieving the technical effect of expanding the calculation range of the slope and improving the calculation accuracy of the slope, and then solving the problem of the slope estimation method in the prior art. Technical problems with low accuracy.
作为一种可选的实施例,根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,可以通过以下步骤实现:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,其中,状态信息包括第三状态信息;根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,其中,状态信息包括第四状态信息;根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益;根据第二状态信息、第三状态信息、第四状态信息和卡尔曼增益,计算得到目标坡度。As an optional embodiment, according to the first state information and the second state information, the target gradient of the vehicle at the second moment is calculated by using the unscented Kalman filter method, which can be realized by the following steps: according to the first state information and the kinematics formula to calculate the third state information of the vehicle at the second moment, wherein the state information includes the third state information; according to the third state information and the dynamics formula, calculate the fourth state information of the vehicle at the second moment , wherein the state information includes the fourth state information; according to the third state information and the fourth state information, calculate the Kalman gain at the second moment; according to the second state information, the third state information, the fourth state information and the Kalman gain Gain, calculated to get the target slope.
无迹卡尔曼滤波方法是一种根据上一时刻的状态向量估计这一时刻下状态向量的估算方法,并且根据观测值来更新或校正这一时刻的状态向量,得到最终的估算值。也就是说,无迹卡尔曼滤波方法在估算这一时刻的状态向量时,需要一个初始的状态向量,还需要一个观测值。The unscented Kalman filter method is an estimation method that estimates the state vector at this moment based on the state vector at the previous moment, and updates or corrects the state vector at this moment according to the observed value to obtain the final estimated value. That is to say, the unscented Kalman filter method needs an initial state vector and an observation value when estimating the state vector at this moment.
需要说明的是,运动学公式为基于车辆行驶过程中的运动学规律而生成的公式,动力学公式为基于车辆行驶中的动力学规律而生成的公式。状态信息中的坡度由于无法直接采用传感器测量得到,因此可以由传感器能够测量得到的其他物理参数根据动力学公式或运动学公式进行推算,得到车辆在某一个时刻下所处的坡度,然而,这样直接的估算并不准确,因此可以采用卡尔曼滤波的方式对车辆在当前时刻下所处的坡度进行更精确的估计。It should be noted that the kinematics formula is a formula generated based on the kinematics law during the running of the vehicle, and the dynamics formula is a formula generated based on the dynamics law during the running of the vehicle. Since the slope in the state information cannot be directly measured by sensors, other physical parameters that can be measured by sensors can be calculated according to dynamic formulas or kinematic formulas to obtain the slope of the vehicle at a certain moment. However, this The direct estimation is not accurate, so the Kalman filter can be used to estimate the slope of the vehicle at the current moment more accurately.
可选地,采用无迹卡尔曼滤波方法对第一状态信息和第二状态信息进行估算得到目标坡度时,第一状态信息和第二状态信息的表示方法可以均为向量表示,并且所有状态信息的向量均包括五个参量,五个参量依次表示车轮角加速度、纵向加速度、车辆所处的坡度、车辆驱动力和车辆速度,在多个状态信息之间进行向量运算时,不改变状态信息的向量格式。第一状态信息是无迹卡尔曼滤波方法中需要的上一时刻的状态向量,通过第一状态信息和车辆在行驶过程中遵循的运动学公式,可以估算车辆在第二时刻下的第三状态信息,再根据第三状态信息和车辆行驶过程中遵循的动力学公式,可以得到第二时刻下的第四状态信息和卡尔曼增益,卡尔曼增益即卡尔曼增益矩阵;根据卡尔曼增益矩阵和作为卡尔曼滤波方法中的观测值的第二状态信息,可以对第二时刻下的第三状态信息进行校正,计算得到包括了目标坡度的最终的状态信息,即车辆在第二时刻下经过卡尔曼滤波方法进行修正后的状态信息。Optionally, when the unscented Kalman filter method is used to estimate the first state information and the second state information to obtain the target slope, the representation methods of the first state information and the second state information can be vector representations, and all state information The vectors of each include five parameters, and the five parameters represent the angular acceleration of the wheel, the longitudinal acceleration, the slope of the vehicle, the driving force of the vehicle, and the vehicle speed in turn. When vector operations are performed between multiple state information, the state information does not change. vector format. The first state information is the state vector at the previous moment required in the unscented Kalman filter method. Through the first state information and the kinematics formula followed by the vehicle during driving, the third state of the vehicle at the second moment can be estimated information, and according to the third state information and the dynamic formula followed by the vehicle during driving, the fourth state information and Kalman gain at the second moment can be obtained, and the Kalman gain is the Kalman gain matrix; according to the Kalman gain matrix and As the second state information of the observed value in the Kalman filter method, the third state information at the second moment can be corrected, and the final state information including the target slope can be calculated, that is, the vehicle passes through Karl at the second moment Mann filter method to modify the state information.
作为一种可选的实施例,根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,可以通过以下步骤实现:对第一状态信息进行无迹变换,生成多个第一采样点及多个第一采样点的加权值;根据运动学公式、多个第一采样点和多个第一采样点的加权值进行计算,得到车辆在第二时刻下的多个第二采样点和多个第二采样点的加权值;根据多个第二采样点和多个第二采样点的加权值,计算得到第三状态信息。As an optional embodiment, the calculation of the third state information of the vehicle at the second moment according to the first state information and the kinematic formula can be realized by the following steps: performing unscented transformation on the first state information to generate multiple The weighted value of a first sampling point and a plurality of first sampling points; calculate according to kinematics formula, a plurality of first sampling points and the weighted value of a plurality of first sampling points, obtain a plurality of vehicle under the second moment The second sampling point and the weighted values of the multiple second sampling points; the third state information is obtained by calculating according to the multiple second sampling points and the weighted values of the multiple second sampling points.
可选地,对第一状态信息进行无迹变换是无迹卡尔曼滤波的一种算法,第一状态信息可以为一个采样点,根据第一状态信息中物理量的误差,可以根据无迹变换的方法得到以第一状态信息为中心的多个第一采样点,并得到多个第一采样点一一对应的多个加权值。将多个第一采样点分别带入运动学公式进行计算,得到多个第二采样点,多个第二采样点对应的加权值与多个第一采样点对应的加权值相同,对应关系也不改变,此时多个第二采样点表示第二时刻下的多个可能的取值。根据多个第二采样点和多个第二采样点一一对应的加权值,可以计算得到多个第二采样点的均值,即第三状态信息。Optionally, performing unscented transformation on the first state information is an algorithm of unscented Kalman filtering, the first state information can be a sampling point, and according to the error of the physical quantity in the first state information, it can be calculated according to the unscented transformation The method obtains a plurality of first sampling points centered on the first state information, and obtains a plurality of weighted values corresponding to each of the plurality of first sampling points. Bring multiple first sampling points into the kinematics formula for calculation to obtain multiple second sampling points. The weighted values corresponding to the multiple second sampling points are the same as the weighted values corresponding to the multiple first sampling points, and the corresponding relationship is also does not change, the multiple second sampling points represent multiple possible values at the second moment. According to the one-to-one correspondence between the plurality of second sampling points and the weighted values of the plurality of second sampling points, an average value of the plurality of second sampling points, that is, the third state information may be obtained.
可选地,本实施例中选取第一采样点的个数为11个,使用的运动学公式如下:Optionally, in this embodiment, the number of the first sampling points is selected as 11, and the kinematic formula used is as follows:
式中:aG,x为车辆质心处的纵向加速度,其中,纵向为车辆坐标系下x轴指向方向,当车辆在水平地面上且处于静止状态下,车辆坐标系的坐标原点与车辆质心重合,x轴平行于地面指向前方,z轴通过质心指向上方,y轴指向驾驶员的左侧;为车辆行驶的纵向位移的二阶导数;g为重力加速度,取值为9.81m/s2。其中,对于位移的二阶导数可以根据如下公式,根据从动轮的车轮角速度进行计算:In the formula: a G, x is the longitudinal acceleration at the center of mass of the vehicle, where the longitudinal direction is the direction of the x-axis in the vehicle coordinate system. When the vehicle is on a level ground and in a static state, the coordinate origin of the vehicle coordinate system coincides with the center of mass of the vehicle , the x-axis is parallel to the ground and points forward, the z-axis points upward through the center of mass, and the y-axis points to the left of the driver; is the second derivative of the longitudinal displacement of the vehicle; g is the acceleration due to gravity, which is 9.81m/s 2 . where, for the second derivative of the displacement It can be calculated according to the wheel angular velocity of the driven wheel according to the following formula:
式中:rstat为车辆在静止状态下滚动半径,单位为m;ω为车轮滚动角速度,单位为rad/s。In the formula: r stat is the rolling radius of the vehicle at rest, in m; ω is the rolling angular velocity of the wheel, in rad/s.
作为一种可选的实施例,接收传感器测得的车轮角加速度,可以通过以下步骤实现:接收传感器测得的车辆的车轮角速度和传感器测量车轮角速度的采样周期;根据车轮角速度和采样周期,计算角加速度。As an optional embodiment, receiving the wheel angular acceleration measured by the sensor can be achieved through the following steps: receiving the wheel angular velocity of the vehicle measured by the sensor and the sampling period of the sensor measuring the wheel angular velocity; according to the wheel angular velocity and the sampling period, calculate angular acceleration.
可选地,动力学公式中,对于离散状态下车轮的角速度求导计算,可以采用如下公式:Optionally, in the dynamic formula, the following formula can be used for the derivative calculation of the angular velocity of the wheel in a discrete state:
式中:ωk为k时刻下车轮的角速度,ωk-1为k-1时刻下车轮的速度;Δt为车轮信号的采样周期,单位为s。In the formula: ω k is the angular velocity of the wheel at time k, ω k-1 is the speed of the wheel at time k-1; Δt is the sampling period of the wheel signal, and the unit is s.
在计算车辆所处坡度时,相关技术中采用小角度估计方法,即在坡度角χroad小于10°时,采用如下公式近似:When calculating the slope of the vehicle, a small angle estimation method is used in the related art, that is, when the slope angle χ road is less than 10°, the following formula is used for approximation:
χroad=sin(χroad)=tan(χroad);χ road = sin(χ road ) = tan(χ road );
cos(χroad)=1。cos(χ road )=1.
本实施例中不采取上述近似方案,直接使用弦函数进行坡度估计,扩大了坡度信号的非线性范围。In this embodiment, the above approximation scheme is not adopted, and the gradient estimation is performed directly using the chord function, which expands the nonlinear range of the gradient signal.
作为一种可选的实施例,根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,可以通过以下步骤实现:对第三状态信息进行无迹变换,生成多个第三采样点及多个第三采样点的加权值;根据动力学公式、多个第三采样点及多个第三采样点的加权值,计算多个第四采样点及多个第四采样点的加权值;根据多个第四采样点及多个第四采样点的加权值,计算得到第四状态信息。As an optional embodiment, according to the third state information and the dynamic formula, the fourth state information of the vehicle at the second moment can be calculated through the following steps: perform unscented transformation on the third state information to generate multiple The weighted value of a third sampling point and a plurality of third sampling points; according to the weighted value of a kinetic formula, a plurality of third sampling points and a plurality of third sampling points, calculate a plurality of fourth sampling points and a plurality of fourth sampling points The weighted value of the sampling point; the fourth state information is obtained by calculating according to the multiple fourth sampling points and the weighted values of the multiple fourth sampling points.
可选地,对第三状态信息也进行无迹变换,以第三状态信息作为中心采样点得到多个第三采样点,以及多个第三采样点一一对应的加权值。将多个第三采样点分别带入动力学公式进行计算,得到多个第四采样点,多个第四采样点对应的加权值与多个第三采样点对应的加权值相同,对应关系也不改变。根据多个第四采样点和多个第四采样点一一对应的加权值,可以计算得到多个第四采样点的均值,即第四状态信息。优选的,本发明选取第一采样点的个数为11个。Optionally, unscented transformation is also performed on the third state information, using the third state information as a central sampling point to obtain a plurality of third sampling points and weighted values corresponding to each of the plurality of third sampling points. Bring multiple third sampling points into the dynamic formula for calculation, and obtain multiple fourth sampling points. The weighted values corresponding to multiple fourth sampling points are the same as the weighted values corresponding to multiple third sampling points, and the corresponding relationship is also do not change. According to the one-to-one correspondence between the multiple fourth sampling points and the weighted values of the multiple fourth sampling points, the average value of the multiple fourth sampling points, that is, the fourth state information, can be calculated. Preferably, the present invention selects 11 first sampling points.
作为一种可选的实施例,获取车辆的第二状态信息,可以通过以下步骤实现:在第二时刻下,接收传感器测得的物理量,其中,传感器测得的物理量包括车辆的纵向加速度、驱动力、速度和车轮角加速度;根据传感器测得的物理量和动力学公式,计算得到第二状态信息包括的第二坡度。As an optional embodiment, obtaining the second state information of the vehicle may be achieved through the following steps: at the second moment, receiving the physical quantity measured by the sensor, wherein the physical quantity measured by the sensor includes the longitudinal acceleration of the vehicle, the driving force, speed and wheel angular acceleration; according to the physical quantity measured by the sensor and the dynamic formula, the second slope included in the second state information is calculated.
可选地,获取车辆的第二状态信息可以根据第二时刻下车辆传感器测得的多个物理量,直接通过如下动力学公式计算得到。Optionally, the acquisition of the second state information of the vehicle may be directly calculated through the following dynamic formula according to multiple physical quantities measured by the vehicle sensors at the second moment.
对于行驶在有坡度的路面,建立汽车的受力平衡方程,即动力学公式如下:For driving on a sloped road, the force balance equation of the car is established, that is, the dynamic formula is as follows:
式中:m为车辆质量,单位为kg;Fx为地面作用于轮胎驱动力,单位为N;cair为空气阻力系数,通过风阻试验可以获得;Al为迎风面积,是车辆固有参数,单位为m2;ρa为空气密度,单位为kg/m3;vG,x为车辆质心速度,单位为m/s。其中,对于驱动力,在车辆处于滑移率较低的情况下时,可以通过车辆自身通过发动机提供的旋转扭矩转换根据如下公式获得,在电动车上,也可以通过电机提供的旋转扭矩根据如下公式获得:In the formula: m is the mass of the vehicle, the unit is kg; F x is the driving force acting on the tires on the ground, the unit is N; c air is the air resistance coefficient, which can be obtained through the wind resistance test; A l is the windward area, which is the inherent parameter of the vehicle, The unit is m 2 ; ρ a is the air density, the unit is kg/m 3 ; v G, x is the velocity of the center of mass of the vehicle, the unit is m/s. Among them, for the driving force, when the vehicle is in the case of low slip ratio, it can be obtained by converting the rotational torque provided by the vehicle itself through the engine according to the following formula. On an electric vehicle, the rotational torque provided by the motor can also be obtained according to the following formula The formula is obtained:
式中:Ttq为发动机转矩,单位N*m;ig为变速器传动比;i0为主减速器传动比;ηT为传动系机械效率;对于驱动力,也可以通过地面提供给车轮的驱动力等于垂向力乘以地面附着系数获得,但考虑到地面附着系数在不同工况下系数大小不等,故不建议使用该方式获得驱动力。In the formula: T tq is the engine torque, unit N*m; i g is the gear ratio of the transmission; i 0 is the gear ratio of the main reducer; η T is the mechanical efficiency of the drive train; the driving force can also be provided to the wheels through the ground The driving force is equal to the vertical force multiplied by the ground adhesion coefficient, but considering that the ground adhesion coefficient varies in different working conditions, it is not recommended to use this method to obtain the driving force.
作为一种可选的实施例,根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益,可以通过以下步骤实现:根据多个第三采样点和多个第四采样点,计算多个第三采样点和多个第四采样点之间的交叉协方差矩阵;根据交叉协方差矩阵,计算得到卡尔曼增益。As an optional embodiment, according to the third state information and the fourth state information, calculating the Kalman gain at the second moment can be realized through the following steps: according to multiple third sampling points and multiple fourth sampling points , to calculate the cross-covariance matrix between the multiple third sampling points and the multiple fourth sampling points; according to the cross-covariance matrix, calculate the Kalman gain.
可选地,根据多个第三采样点和多个第四采样点,可以根据无迹卡尔曼滤波方法,因为多个第四采样点分别由多个第三采样点得到,所以多个第四采样点和多个第三采样点之间存在一一对应关系,根据多个第四采样点和多个第三采样点之间存在一一对应关系,可以计算得到多个第三采样点和多个第四采样点之间的交叉协方差矩阵,采用交叉协方差矩阵乘交叉协方差矩阵的逆矩阵,得到卡尔曼增益矩阵。Optionally, according to the plurality of third sampling points and the plurality of fourth sampling points, the unscented Kalman filtering method can be used, because the plurality of fourth sampling points are respectively obtained from the plurality of third sampling points, so the plurality of fourth sampling points There is a one-to-one correspondence between the sampling points and the multiple third sampling points. According to the one-to-one correspondence between the multiple fourth sampling points and the multiple third sampling points, the multiple third sampling points and the multiple The cross-covariance matrix between the fourth sampling points is obtained by multiplying the cross-covariance matrix by the inverse matrix of the cross-covariance matrix to obtain the Kalman gain matrix.
作为一个具体的实施例,aG,x,Fx三个量是已知量,均可以通过车辆传感器采集或车辆系统中存储的相关数据带入上述给出的公式得到。As a specific example, The three quantities of a G, x , and F x are known quantities, which can be obtained by bringing the relevant data collected by the vehicle sensor or stored in the vehicle system into the formula given above.
对于上述运动学公式和动力学公式,设置状态向量X如下:For the above kinematics formula and dynamics formula, set the state vector X as follows:
对于三个已知量,假设其误差服从正态分布,误差公式如下:For the three known quantities, assuming that their errors obey the normal distribution, the error formula is as follows:
vω,k~N(0,Qk);v ω, k ~ N(0, Q k );
Va,k~N(0,Rk);V a, k ~ N (0, R k );
vF,k~N(0,Tk);v F, k ~ N (0, T k );
上述误差公式中,vω,k为aG,x,Fx的误差,va,k为aG,x的误差,vF,k为Fx的误差,Qk,Rk,Tk分别为vω,k,va,k,vF,k服从的正态分布的方差。In the above error formula, v ω, k is a G, x , the error of F x , v a, k is the error of a G, x , v F, k is the error of F x , Q k , R k , T k are respectively v ω, k , v a, k , v F, the variance of the normal distribution that k obeys.
离散系统中,k+1时刻的状态向量和k时刻的状态向量的关系如下:In a discrete system, the relationship between the state vector at time k+1 and the state vector at time k is as follows:
初始状态下,Qk,Rk,Tk取值为0.01。In the initial state, the values of Q k , R k , and T k are 0.01.
无迹卡尔曼滤波计算过程基本计算流程是:The basic calculation process of the unscented Kalman filter calculation process is:
首先,获取第一状态信息的向量表示,并确定采样点的数量和权值,本可选实施例中,采样点个数是11个;First, obtain the vector representation of the first state information, and determine the number and weight of sampling points. In this optional embodiment, the number of sampling points is 11;
其次,根据第一状态信息中各个物理量的误差,计算第一状态信息的协方差矩阵的科列斯基分解;Secondly, according to the error of each physical quantity in the first state information, calculate the Kolesky decomposition of the covariance matrix of the first state information;
其次,根据上一步中计算得到的协方差矩阵的科列斯基分解、第一状态信息和第一步中确定的采样点的数量和权值,计算多个第一采样点;Secondly, according to the Kolesky decomposition of the covariance matrix calculated in the previous step, the first state information and the number and weight of the sampling points determined in the first step, calculate a plurality of first sampling points;
其次,根据k+1时刻的状态向量和k时刻的状态向量的关系,和第一状态信息的向量表示,计算多个第一采样点在第二时刻下的数值,并分别将多个第一采样点在第二时刻下的数值带入运动学公式中计算,更新多个第一采样点中的aG,x,χroad的值,得到多个第二采样点;Secondly, according to the relationship between the state vector at time k+1 and the state vector at time k, and the vector representation of the first state information, calculate the values of the multiple first sampling points at the second time, and divide the multiple first sampling points into The value of the sampling point at the second moment is brought into the kinematics formula for calculation, and the values in multiple first sampling points are updated a G, x , the value of χ road , obtain a plurality of second sampling points;
其次,根据多个第二采样点和最初确定的权值,求解运动学公式算出的多个第二采样点的均值和协方差矩阵,其中,多个第二采样点的均值即为第三状态向量;Secondly, according to the plurality of second sampling points and the initially determined weights, solve the mean value and covariance matrix of the plurality of second sampling points calculated by the kinematic formula, wherein the mean value of the plurality of second sampling points is the third state vector;
其次,根据上一步中得到的协方差矩阵进行科列斯基分解;Secondly, perform Koleski decomposition according to the covariance matrix obtained in the previous step;
其次,根据第一步确定的采样个数和上一步中计算得到的协方差矩阵的科列斯基分解,对第三状态向量进行采样,得到多个第三采样点;Secondly, according to the number of samples determined in the first step and the Kolesky decomposition of the covariance matrix calculated in the previous step, the third state vector is sampled to obtain a plurality of third sampling points;
其次,将多个第三采样点分别代入动力学公式计算,更新多个第三采样点中的aG,x,χroad,Fx,vG,x的值,得到多个第四采样点,并计算多个第四采样点的均值和协方差矩阵,其中,多个第四采样点的均值即为第四状态向量;Secondly, a plurality of third sampling points are respectively substituted into the kinetic formula calculation, and the values of a G, x , χ road , F x , v G, x in the plurality of third sampling points are updated to obtain a plurality of fourth sampling points , and calculate the mean value and covariance matrix of a plurality of fourth sampling points, wherein, the mean value of a plurality of fourth sampling points is the fourth state vector;
其次,根据第三状态向量和第四状态向量,计算运动学公式和动力学公式分别得到的结果的交叉协方差矩阵和卡尔曼增益;Secondly, according to the third state vector and the fourth state vector, calculate the cross-covariance matrix and Kalman gain of the results obtained by the kinematic formula and the dynamic formula respectively;
其次,获取第二时刻下车辆传感器的测量值,直接根据动力学公式计算得到第二状态向量作为观测值;Secondly, the measured value of the vehicle sensor at the second moment is obtained, and the second state vector is directly calculated according to the dynamic formula as the observed value;
最后,第二状态向量减去第四状态向量的值,乘以卡尔曼增益,再加第三状态向量,即可得到最后的更新的状态向量,最后的更新的状态向量中的坡度即为目标坡度。Finally, subtract the value of the fourth state vector from the second state vector, multiply by the Kalman gain, and add the third state vector to obtain the final updated state vector, and the slope in the final updated state vector is the target slope.
需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明并不受所描述的动作顺序的限制,因为依据本发明,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本发明所必须的。It should be noted that for the foregoing method embodiments, for the sake of simple description, they are expressed as a series of action combinations, but those skilled in the art should know that the present invention is not limited by the described action sequence. Because of the present invention, certain steps may be performed in other orders or simultaneously. Secondly, those skilled in the art should also know that the embodiments described in the specification belong to preferred embodiments, and the actions and modules involved are not necessarily required by the present invention.
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到根据上述实施例的数据处理方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。Through the description of the above embodiments, those skilled in the art can clearly understand that the data processing method according to the above embodiments can be implemented by means of software plus a necessary general-purpose hardware platform, and of course also by hardware, but in many cases The former is a better implementation. Based on such an understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products are stored in a storage medium (such as ROM/RAM, disk, CD) contains several instructions to enable a terminal device (which may be a mobile phone, a computer, a server, or a network device, etc.) to execute the methods described in various embodiments of the present invention.
根据本发明实施例,还提供了一种用于实施上述数据处理方法的装置,图3是根据本发明实施例提供的数据处理装置的结构框图,如图3所示,该数据处理装置包括:获取模块32和计算模块34,下面对该数据处理装置进行说明。According to an embodiment of the present invention, a device for implementing the above data processing method is also provided. FIG. 3 is a structural block diagram of a data processing device provided according to an embodiment of the present invention. As shown in FIG. 3 , the data processing device includes: The acquisition module 32 and the calculation module 34, the data processing device will be described below.
此处需要说明的是,上述获取模块32和计算模块34对应于实施例中的步骤S202至步骤S204,二个模块与对应的步骤所实现的实例和应用场景相同,但不限于上述实施例所公开的内容。需要说明的是,上述模块作为装置的一部分可以运行在实施例提供的计算机终端10中。It should be noted here that the acquisition module 32 and the calculation module 34 correspond to steps S202 to S204 in the embodiment, and the examples and application scenarios realized by the two modules are the same as those of the corresponding steps, but they are not limited to those in the above embodiment. public content. It should be noted that, as a part of the device, the above modules can run in the computer terminal 10 provided in the embodiment.
本发明的实施例可以提供一种计算机设备,可选地,在本实施例中,上述计算机设备可以位于计算机网络的多个网络设备中的至少一个网络设备。该计算机设备包括存储器和处理器。An embodiment of the present invention may provide a computer device. Optionally, in this embodiment, the above computer device may be located in at least one network device among multiple network devices in a computer network. The computer device includes memory and a processor.
其中,存储器可用于存储软件程序以及模块,如本发明实施例中的数据处理方法和装置对应的程序指令/模块,处理器通过运行存储在存储器内的软件程序以及模块,从而执行各种功能应用以及数据处理,即实现上述的数据处理方法。存储器可包括高速随机存储器,还可以包括非易失性存储器,如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中,存储器可进一步包括相对于处理器远程设置的存储器,这些远程存储器可以通过网络连接至计算机终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。Among them, the memory can be used to store software programs and modules, such as the program instructions/modules corresponding to the data processing method and device in the embodiment of the present invention, and the processor executes various functional applications by running the software programs and modules stored in the memory. And data processing, that is, realizing the above-mentioned data processing method. The memory may include high-speed random access memory, and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some instances, the memory may further include a memory located remotely from the processor, and these remote memories may be connected to the computer terminal through a network. Examples of the aforementioned networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.
处理器可以通过传输装置调用存储器存储的信息及应用程序,以执行下述步骤:获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。The processor can call the information and application programs stored in the memory through the transmission device to perform the following steps: acquire the state information of the vehicle, wherein the state information includes the first state information of the vehicle at the first moment and the state information of the vehicle at the second moment. The second state information of the state information, any one of the state information includes: speed, longitudinal acceleration, slope, driving force and wheel angular acceleration, the direction of the longitudinal acceleration is from the rear of the vehicle to the front of the vehicle; according to the first state information and the second For state information, an unscented Kalman filter method is used to calculate the target slope of the vehicle at the second moment, wherein the unscented Kalman filter method uses the vehicle's kinematics formula and dynamics formula to estimate the target slope.
可选的,上述处理器还可以执行如下步骤的程序代码:根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,包括:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,其中,状态信息包括第三状态信息;根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,其中,状态信息包括第四状态信息;根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益;根据第二状态信息、第三状态信息、第四状态信息和卡尔曼增益,计算得到目标坡度。Optionally, the above-mentioned processor can also execute the program code of the following steps: according to the first state information and the second state information, calculate the target slope of the vehicle at the second moment by using the unscented Kalman filter method, including: according to the first A state information and kinematics formula, calculating the third state information of the vehicle at the second moment, wherein the state information includes the third state information; according to the third state information and the dynamics formula, calculating the third state information of the vehicle at the second moment Four state information, wherein the state information includes the fourth state information; according to the third state information and the fourth state information, calculate the Kalman gain at the second moment; according to the second state information, the third state information, the fourth state information and Kalman gain to calculate the target slope.
可选的,上述处理器还可以执行如下步骤的程序代码:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,包括:对第一状态信息进行无迹变换,生成多个第一采样点及多个第一采样点的加权值;根据运动学公式、多个第一采样点和多个第一采样点的加权值进行计算,得到车辆在第二时刻下的多个第二采样点和多个第二采样点的加权值;根据多个第二采样点和多个第二采样点的加权值,计算得到第三状态信息。Optionally, the above-mentioned processor may also execute the program code of the following steps: calculating the third state information of the vehicle at the second moment according to the first state information and the kinematic formula, including: performing unscented transformation on the first state information , generate multiple first sampling points and weighted values of multiple first sampling points; calculate according to the kinematics formula, multiple first sampling points and multiple first sampling points’ weighted values, and obtain the vehicle at the second moment The plurality of second sampling points and the weighted values of the plurality of second sampling points; according to the plurality of second sampling points and the weighted values of the plurality of second sampling points, the third state information is obtained through calculation.
可选的,上述处理器还可以执行如下步骤的程序代码:根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,包括:对第三状态信息进行无迹变换,生成多个第三采样点及多个第三采样点的加权值;根据动力学公式、多个第三采样点及多个第三采样点的加权值,计算多个第四采样点及多个第四采样点的加权值;根据多个第四采样点及多个第四采样点的加权值,计算得到第四状态信息。Optionally, the above-mentioned processor can also execute the program code of the following steps: calculating the fourth state information of the vehicle at the second moment according to the third state information and the dynamic formula, including: performing unscented transformation on the third state information , generate multiple third sampling points and weighted values of multiple third sampling points; calculate multiple fourth sampling points and multiple weighted values according to the dynamic formula, multiple third sampling points and multiple third sampling points A weighted value of a fourth sampling point; the fourth state information is obtained by calculating according to the multiple fourth sampling points and the weighted values of the multiple fourth sampling points.
可选的,上述处理器还可以执行如下步骤的程序代码:根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益,包括:根据多个第三采样点和多个第四采样点,计算多个采样点和多个第四采样点之间的交叉协方差矩阵;根据交叉协方差矩阵,计算得到卡尔曼增益。Optionally, the above-mentioned processor may also execute the program code of the following steps: calculating the Kalman gain at the second moment according to the third state information and the fourth state information, including: according to multiple third sampling points and multiple first sampling points Four sampling points, calculating the cross-covariance matrix between multiple sampling points and multiple fourth sampling points; according to the cross-covariance matrix, calculate the Kalman gain.
可选的,上述处理器还可以执行如下步骤的程序代码:获取车辆的第二状态信息,包括:在第二时刻下,接收传感器测得的物理量,其中,传感器测得的物理量包括车辆的纵向加速度、驱动力、速度和车轮角加速度;根据传感器测得的物理量和动力学公式,计算得到第二状态信息包括的第二坡度。Optionally, the above-mentioned processor may also execute the program code of the following steps: obtaining the second state information of the vehicle, including: receiving the physical quantity measured by the sensor at the second moment, wherein the physical quantity measured by the sensor includes the longitudinal direction of the vehicle Acceleration, driving force, speed and wheel angular acceleration; according to the physical quantity measured by the sensor and the dynamics formula, the second slope included in the second state information is calculated.
可选的,上述处理器还可以执行如下步骤的程序代码:接收传感器测得的车轮角加速度,包括:接收传感器测得的车辆的车轮角速度和传感器测量车轮角速度的采样周期;根据车轮角速度和采样周期,计算角加速度。Optionally, the above-mentioned processor can also execute the program code of the following steps: receiving the wheel angular acceleration measured by the sensor, including: receiving the wheel angular velocity of the vehicle measured by the sensor and the sampling period of the sensor measuring the wheel angular velocity; period, to calculate the angular acceleration.
采用本发明实施例,提供了一种数据处理的方案。通过获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计,达到了将车辆状态信息数据进行融合计算坡度的目的,从而实现了扩大坡度计算范围且提高坡度计算精度的技术效果,进而解决了现有技术中坡度估算方法精度低的技术问题。By adopting the embodiment of the present invention, a data processing solution is provided. By acquiring the state information of the vehicle, wherein the state information includes the first state information of the vehicle at the first moment and the second state information of the vehicle at the second moment, any one of the state information includes: speed, longitudinal acceleration, The direction of the slope, driving force, wheel angular acceleration, and longitudinal acceleration is from the rear to the front; according to the first state information and the second state information, the target of the vehicle at the second moment is calculated by using the unscented Kalman filter method Slope, among them, the unscented Kalman filter method uses the kinematics formula and dynamics formula of the vehicle to estimate the target slope, and achieves the purpose of merging the vehicle state information data to calculate the slope, thereby realizing the expansion of the slope calculation range and improving the slope. The technical effect of the calculation accuracy further solves the technical problem of low accuracy of the slope estimation method in the prior art.
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令终端设备相关的硬件来完成,该程序可以存储于一非易失性存储介质中,存储介质可以包括:闪存盘、只读存储器(Read-Only Memory,ROM)、随机存取器(RandomAccess Memory,RAM)、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing the hardware related to the terminal device through a program, and the program can be stored in a non-volatile storage medium, the storage medium It may include: a flash disk, a read-only memory (Read-Only Memory, ROM), a random access device (Random Access Memory, RAM), a magnetic disk or an optical disk, and the like.
本发明的实施例还提供了一种非易失性存储介质。可选地,在本实施例中,上述非易失性存储介质可以用于保存上述实施例所提供的数据处理方法所执行的程序代码。The embodiment of the present invention also provides a non-volatile storage medium. Optionally, in this embodiment, the above-mentioned non-volatile storage medium may be used to store the program code executed by the data processing method provided by the above-mentioned embodiment.
可选地,在本实施例中,上述非易失性存储介质可以位于计算机网络中计算机终端群中的任意一个计算机终端中,或者位于移动终端群中的任意一个移动终端中。Optionally, in this embodiment, the above-mentioned non-volatile storage medium may be located in any computer terminal in the group of computer terminals in the computer network, or in any mobile terminal in the group of mobile terminals.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:获取车辆的状态信息,其中,状态信息包括车辆在第一时刻下的第一状态信息和车辆在第二时刻下的第二状态信息,状态信息中的任意之一包括:速度、纵向加速度、所处坡度、驱动力和车轮角加速度,纵向加速度的方向为由车尾指向车头;根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,其中,无迹卡尔曼滤波方法采用车辆的运动学公式和动力学公式对目标坡度进行估计。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: acquiring state information of the vehicle, wherein the state information includes the first state of the vehicle at the first moment Information and the second state information of the vehicle at the second moment, any one of the state information includes: speed, longitudinal acceleration, slope, driving force and wheel angular acceleration, and the direction of the longitudinal acceleration is from the rear of the vehicle to the front of the vehicle; According to the first state information and the second state information, the unscented Kalman filter method is used to calculate the target slope of the vehicle at the second moment. The slope is estimated.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:根据第一状态信息和第二状态信息,采用无迹卡尔曼滤波方法计算得到车辆在第二时刻下的目标坡度,包括:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,其中,状态信息包括第三状态信息;根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,其中,状态信息包括第四状态信息;根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益;根据第二状态信息、第三状态信息、第四状态信息和卡尔曼增益,计算得到目标坡度。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: According to the first state information and the second state information, the vehicle is calculated by using the unscented Kalman filter method The target slope at the second moment includes: calculating the third state information of the vehicle at the second moment according to the first state information and the kinematic formula, wherein the state information includes the third state information; according to the third state information and The dynamics formula calculates the fourth state information of the vehicle at the second moment, wherein the state information includes the fourth state information; calculates the Kalman gain at the second moment according to the third state information and the fourth state information; The second state information, the third state information, the fourth state information and the Kalman gain are used to calculate the target slope.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:根据第一状态信息和运动学公式,计算车辆在第二时刻下的第三状态信息,包括:对第一状态信息进行无迹变换,生成多个第一采样点及多个第一采样点的加权值;根据运动学公式、多个第一采样点和多个第一采样点的加权值进行计算,得到车辆在第二时刻下的多个第二采样点和多个第二采样点的加权值;根据多个第二采样点和多个第二采样点的加权值,计算得到第三状态信息。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: calculate the third state of the vehicle at the second moment according to the first state information and the kinematics formula Information, including: performing unscented transformation on the first state information to generate a plurality of first sampling points and weighted values of the plurality of first sampling points; Calculate the weighted value of the vehicle to obtain the weighted values of multiple second sampling points and multiple second sampling points of the vehicle at the second moment; according to the weighted values of multiple second sampling points and multiple second sampling points, calculate Get the third state information.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:根据第三状态信息和动力学公式,计算车辆在第二时刻下的第四状态信息,包括:对第三状态信息进行无迹变换,生成多个第三采样点及多个第三采样点的加权值;根据动力学公式、多个第三采样点及多个第三采样点的加权值,计算多个第四采样点及多个第四采样点的加权值;根据多个第四采样点及多个第四采样点的加权值,计算得到第四状态信息。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: calculate the fourth state of the vehicle at the second moment according to the third state information and the dynamics formula information, including: performing unscented transformation on the third state information to generate a plurality of third sampling points and weighted values of the plurality of third sampling points; according to a dynamic formula, a plurality of third sampling points and a plurality of third sampling points Calculate a plurality of fourth sampling points and weighted values of the plurality of fourth sampling points; calculate and obtain fourth state information according to the plurality of fourth sampling points and the weighted values of the plurality of fourth sampling points.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:根据第三状态信息和第四状态信息,计算第二时刻下的卡尔曼增益,包括:根据多个第三采样点和多个第四采样点,计算多个采样点和多个第四采样点之间的交叉协方差矩阵;根据交叉协方差矩阵,计算得到卡尔曼增益。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: calculating the Kalman gain at the second moment according to the third state information and the fourth state information, The method includes: calculating a cross-covariance matrix between the multiple sampling points and the multiple fourth sampling points according to the multiple third sampling points and the multiple fourth sampling points; and calculating a Kalman gain according to the cross-covariance matrix.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:获取车辆的第二状态信息,包括:在第二时刻下,接收传感器测得的物理量,其中,传感器测得的物理量包括车辆的纵向加速度、驱动力、速度和车轮角加速度;根据传感器测得的物理量和动力学公式,计算得到第二状态信息包括的第二坡度。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: obtaining the second state information of the vehicle, including: at the second moment, receiving the The physical quantity, wherein the physical quantity measured by the sensor includes the vehicle's longitudinal acceleration, driving force, speed and wheel angular acceleration; according to the physical quantity measured by the sensor and the dynamics formula, the second slope included in the second state information is calculated.
可选地,在本实施例中,非易失性存储介质被设置为存储用于执行以下步骤的程序代码:接收传感器测得的车轮角加速度,包括:接收传感器测得的车辆的车轮角速度和传感器测量车轮角速度的采样周期;根据车轮角速度和采样周期,计算角加速度。Optionally, in this embodiment, the non-volatile storage medium is configured to store program codes for performing the following steps: receiving the wheel angular acceleration measured by the sensor, including: receiving the wheel angular velocity and The sensor measures the sampling period of the wheel angular velocity; according to the wheel angular velocity and the sampling period, the angular acceleration is calculated.
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。The serial numbers of the above embodiments of the present invention are for description only, and do not represent the advantages and disadvantages of the embodiments.
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above-mentioned embodiments of the present invention, the descriptions of each embodiment have their own emphases, and for parts not described in detail in a certain embodiment, reference may be made to relevant descriptions of other embodiments.
在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed technical content can be realized in other ways. Wherein, the device embodiments described above are only illustrative. For example, the division of the units may be a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or may be Integrate into another system, or some features may be ignored, or not implemented. In another point, the mutual coupling or direct coupling or communication connection shown or discussed may be through some interfaces, and the indirect coupling or communication connection of units or modules may be in electrical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple units. Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个非易失性取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。If the integrated unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a non-volatile storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium , including several instructions to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage media include: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes. .
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only a preferred embodiment of the present invention, it should be pointed out that, for those of ordinary skill in the art, without departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications can also be made. It should be regarded as the protection scope of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211254981.3A CN115447587B (en) | 2022-10-13 | 2022-10-13 | Data processing method, device, non-volatile storage medium and computer equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211254981.3A CN115447587B (en) | 2022-10-13 | 2022-10-13 | Data processing method, device, non-volatile storage medium and computer equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115447587A true CN115447587A (en) | 2022-12-09 |
CN115447587B CN115447587B (en) | 2025-04-01 |
Family
ID=84308581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211254981.3A Active CN115447587B (en) | 2022-10-13 | 2022-10-13 | Data processing method, device, non-volatile storage medium and computer equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115447587B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116985825A (en) * | 2023-09-26 | 2023-11-03 | 安徽蔚来智驾科技有限公司 | Method, system and storage medium for determining vehicle state |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140067155A1 (en) * | 2012-08-31 | 2014-03-06 | Ford Global Technologies, Llc | Dynamic road gradient estimation |
US20160332633A1 (en) * | 2013-07-11 | 2016-11-17 | C.R.F. Societa' Consortile Per Azioni | Automotive control unit programmed to estimate road slope and vehicle mass, vehicle with such a control unit and corresponding program product |
CN109684677A (en) * | 2018-12-04 | 2019-04-26 | 西安法士特汽车传动有限公司 | A kind of gradient evaluation method based on Kalman filtering algorithm |
CN111114551A (en) * | 2018-10-31 | 2020-05-08 | 广州汽车集团股份有限公司 | Vehicle ramp gradient identification method and device |
CN111391857A (en) * | 2019-12-24 | 2020-07-10 | 北京理工大学 | A method and system for slope estimation of distributed drive electric vehicle |
CN112429010A (en) * | 2020-12-02 | 2021-03-02 | 东风商用车有限公司 | Method for estimating vehicle mass and road gradient |
CN113978473A (en) * | 2021-11-25 | 2022-01-28 | 中山大学 | A Method for Estimating Vehicle Mass and Road Slope |
CN114750769A (en) * | 2021-12-14 | 2022-07-15 | 北京理工大学 | Joint estimation method and system for vehicle mass and road gradient |
-
2022
- 2022-10-13 CN CN202211254981.3A patent/CN115447587B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140067155A1 (en) * | 2012-08-31 | 2014-03-06 | Ford Global Technologies, Llc | Dynamic road gradient estimation |
US20160332633A1 (en) * | 2013-07-11 | 2016-11-17 | C.R.F. Societa' Consortile Per Azioni | Automotive control unit programmed to estimate road slope and vehicle mass, vehicle with such a control unit and corresponding program product |
CN111114551A (en) * | 2018-10-31 | 2020-05-08 | 广州汽车集团股份有限公司 | Vehicle ramp gradient identification method and device |
CN109684677A (en) * | 2018-12-04 | 2019-04-26 | 西安法士特汽车传动有限公司 | A kind of gradient evaluation method based on Kalman filtering algorithm |
CN111391857A (en) * | 2019-12-24 | 2020-07-10 | 北京理工大学 | A method and system for slope estimation of distributed drive electric vehicle |
CN112429010A (en) * | 2020-12-02 | 2021-03-02 | 东风商用车有限公司 | Method for estimating vehicle mass and road gradient |
CN113978473A (en) * | 2021-11-25 | 2022-01-28 | 中山大学 | A Method for Estimating Vehicle Mass and Road Slope |
CN114750769A (en) * | 2021-12-14 | 2022-07-15 | 北京理工大学 | Joint estimation method and system for vehicle mass and road gradient |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116985825A (en) * | 2023-09-26 | 2023-11-03 | 安徽蔚来智驾科技有限公司 | Method, system and storage medium for determining vehicle state |
CN116985825B (en) * | 2023-09-26 | 2024-03-19 | 安徽蔚来智驾科技有限公司 | Method, system and storage medium for determining vehicle state |
Also Published As
Publication number | Publication date |
---|---|
CN115447587B (en) | 2025-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109532555B (en) | Method and device for calculating endurance mileage | |
US20090157365A1 (en) | Simulation device and simulation method | |
CN110596606B (en) | A method, system and device for estimating remaining power of lithium battery | |
CN104890526A (en) | Vehicle system and method | |
CN106515740A (en) | Distributed electrically driven automobile travelling status parameter estimation algorithm based on ICDKF | |
CN113459896B (en) | Method and device for determining endurance mileage of electric vehicle | |
KR20150020270A (en) | Estimating the state of charge of a battery | |
CN110884499B (en) | A method and system for determining vehicle center of mass sideslip angle | |
CN114572224B (en) | Estimation method and terminal for maximum adhesion coefficient of road surface | |
CN116552548B (en) | A state estimation method for four-wheel distributed electric drive vehicles | |
WO2022024847A1 (en) | Computing system, battery deterioration predicting method, and battery deterioration predicting program | |
CN116702096B (en) | Method and device for measuring and calculating road sliding resistance of vehicle plateau environment | |
CN114206664B (en) | Vehicle-mounted notification device, recording medium, and computing device | |
CN110288154A (en) | Speed predicting method, device, equipment and medium | |
CN115447587A (en) | Data processing method, device, non-volatile storage medium and computer equipment | |
CN110361081A (en) | A kind of vehicular gross combined weight dynamic measurement method based on recurrent least square method | |
CN114778141A (en) | Method, device and equipment for formulating durability test of automobile electric drive assembly rack | |
CN111198032A (en) | Real-time estimation method for automobile mass | |
CN108287934B (en) | Vehicle mass center slip angle robust estimation method based on longitudinal force observer | |
CN111942399A (en) | Vehicle speed estimation method and system based on unscented Kalman filtering | |
CN114701508A (en) | Vehicle mass estimation method, apparatus, device, storage medium, and program product | |
CN102706364A (en) | Online calibration method of scaling factors of micromachining gyroscope for vehicle | |
WO2023035234A1 (en) | Vehicle state parameter estimation method and apparatus | |
CN112977082B (en) | Torque distribution method and device, modeling method and automobile | |
CN116161046B (en) | Method and device for determining slope estimation value, storage medium and electronic device |
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 |