[go: up one dir, main page]

CN118114498A - Data filtering method, device, computer equipment, storage medium and program product - Google Patents

Data filtering method, device, computer equipment, storage medium and program product Download PDF

Info

Publication number
CN118114498A
CN118114498A CN202410328547.8A CN202410328547A CN118114498A CN 118114498 A CN118114498 A CN 118114498A CN 202410328547 A CN202410328547 A CN 202410328547A CN 118114498 A CN118114498 A CN 118114498A
Authority
CN
China
Prior art keywords
moment
state data
priori
data
motion state
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.)
Pending
Application number
CN202410328547.8A
Other languages
Chinese (zh)
Inventor
张乾坤
苏瑞
宋艳辉
刘勇敢
张玉龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Geely Holding Group Co Ltd
Radar New Energy Vehicle Zhejiang Co Ltd
Original Assignee
Zhejiang Geely Holding Group Co Ltd
Radar New Energy Vehicle Zhejiang Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Geely Holding Group Co Ltd, Radar New Energy Vehicle Zhejiang Co Ltd filed Critical Zhejiang Geely Holding Group Co Ltd
Priority to CN202410328547.8A priority Critical patent/CN118114498A/en
Publication of CN118114498A publication Critical patent/CN118114498A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Navigation (AREA)

Abstract

The invention relates to the technical field of permanent magnet synchronous motors, and discloses a data filtering method, a device, computer equipment, a storage medium and a program product, wherein the data filtering method is applied to a Kalman filter and comprises the following steps: acquiring motion state data of a target vehicle in running, wherein the motion state data comprise motion state data of a first moment and a second moment; determining a first prior estimated covariance of the current moment based on the motion state data; and filtering the first priori estimation data at the current moment based on the first priori estimation covariance to obtain the optimal estimation of the motion state data at the current moment. Based on the method, the process noise Q can be not predicted, but the first priori estimated covariance of the current moment is determined through the motion state data of the first moment and the second moment, so that the accuracy of the optimal estimation of the motion state data of the current moment is improved on the basis of getting rid of the dependence on the accuracy Q.

Description

数据滤波方法、装置、计算机设备、存储介质及程序产品Data filtering method, device, computer equipment, storage medium and program product

技术领域Technical Field

本发明涉及永磁同步电机技术领域,具体涉及数据滤波方法、装置、计算机设备、存储介质及程序产品。The present invention relates to the technical field of permanent magnet synchronous motors, and in particular to a data filtering method, device, computer equipment, storage medium and program product.

背景技术Background technique

在新能源汽车开发的过程中,采用的电机通常为永磁同步电机,从而在降低车辆能耗的同时,延长了电池的续航里程。然而,考虑到永磁同步电机驱动器的磁场定向矢量控制方法会在同步参考系中产生正负六次谐波电流,导致了扭矩波动和控制性能恶化。In the process of developing new energy vehicles, the motors used are usually permanent magnet synchronous motors, which reduce vehicle energy consumption and extend the battery life. However, considering that the field-oriented vector control method of the permanent magnet synchronous motor drive will generate positive and negative sixth harmonic currents in the synchronous reference system, it leads to torque fluctuations and deterioration of control performance.

基于此,在新能源汽车领域通常使用卡尔曼滤波器对谐波电流进行滤除,通常来说,卡尔曼滤波器的滤波效果依赖于针对过程噪声Q的预测,然而,在汽车行驶过程中,由于无法预测行驶路况,因此无法得到准确的过程噪声Q,从而影响了卡尔曼滤波器的预测精度。Based on this, Kalman filters are usually used in the field of new energy vehicles to filter out harmonic currents. Generally speaking, the filtering effect of the Kalman filter depends on the prediction of the process noise Q. However, during the driving of the car, the road conditions cannot be predicted, so the accurate process noise Q cannot be obtained, which affects the prediction accuracy of the Kalman filter.

发明内容Summary of the invention

有鉴于此,本发明提供了一种数据滤波方法、装置、计算机设备、存储介质及程序产品,以解决新能源汽车行驶过程中无法确定准确的过程噪声Q从而影响了卡尔曼滤波器的预测精度的问题。In view of this, the present invention provides a data filtering method, device, computer equipment, storage medium and program product to solve the problem that the accurate process noise Q cannot be determined during the driving of new energy vehicles, thereby affecting the prediction accuracy of the Kalman filter.

第一方面,本发明提供了一种数据滤波方法,应用于卡尔曼滤波器,该方法包括:In a first aspect, the present invention provides a data filtering method, applied to a Kalman filter, the method comprising:

获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据;Acquire motion state data of the target vehicle during driving, wherein the motion state data includes motion state data at a first moment and a second moment;

基于运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为所述第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差;Determine a first a priori estimation covariance at a current moment based on the motion state data, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimation covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori;

基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,第一先验估计数据用于指示对所述当前时刻运动状态数据的先验估计。The first a priori estimation data at the current moment is filtered based on the first a priori estimation covariance to obtain an optimal estimate of the motion state data at the current moment, wherein the first a priori estimation data is used to indicate a priori estimation of the motion state data at the current moment.

在一种可选的实施方式中,基于运动状态数据,确定当前时刻的第一先验估计协方差,包括:In an optional implementation, determining the first a priori estimated covariance at the current moment based on the motion state data includes:

基于运动状态数据,确定第一时刻对应的先验调整量,其中,先验调整量用于指示对第一时刻的第二先验估计协方差的误差进行修正时的调整量;Determine, based on the motion state data, a priori adjustment amount corresponding to the first moment, wherein the priori adjustment amount is used to indicate an adjustment amount when correcting an error of the second priori estimate covariance at the first moment;

根据先验调整量与所述第二先验估计协方差的和确定当前时刻的第一先验估计协方差。The first a priori estimation covariance at the current moment is determined according to the sum of the a priori adjustment amount and the second a priori estimation covariance.

在本公开实施例中,在相关的卡尔曼滤波算法中需要通过噪声Q来计算当前时刻k对应的第一先验估计协方差,而在本公开中,计算第一先验估计协方差时舍弃了噪声Q,从而提高了确定出的第一先验估计协方差的精度,进而提高了对当前时刻运动状态数据的最优估计的准确性。In an embodiment of the present disclosure, in the relevant Kalman filter algorithm, noise Q is required to calculate the first prior estimate covariance corresponding to the current moment k. However, in the present disclosure, noise Q is discarded when calculating the first prior estimate covariance, thereby improving the accuracy of the determined first prior estimate covariance, and further improving the accuracy of the optimal estimate of the motion state data at the current moment.

在一种可选的实施方式中,第一时刻的运动状态数据包括第一时刻对应的后验状态数据,其中,后验状态数据用于指示所述第一时刻对应的最优估计;In an optional implementation, the motion state data at the first moment includes a posteriori state data corresponding to the first moment, wherein the a posteriori state data is used to indicate the optimal estimate corresponding to the first moment;

基于所述运动状态数据,确定第一时刻对应的先验调整量,包括:Determining a priori adjustment amount corresponding to a first moment based on the motion state data includes:

获取预先基于第二时刻的运动状态数据确定出的第一时刻的第二先验估计协方差;Obtaining a second a priori estimated covariance at the first moment that is determined in advance based on the motion state data at the second moment;

基于后验状态数据与第二先验估计数据的差值,确定第一时刻对应的后验残差;Determine a posterior residual corresponding to the first moment based on a difference between the posterior state data and the second a priori estimated data;

基于后验残差对第二先验估计协方差进行修正,得到第一时刻对应的先验调整量。The second a priori estimated covariance is corrected based on the posterior residual to obtain a priori adjustment amount corresponding to the first moment.

在本公开实施例中,可以基于第一时刻以及第二时刻对应的运算数据来计算第一时刻对应的先验调整量,以基于该先验调整量来计算当前时刻k对应的第一先验估计协方差,从而摆脱计算过程中对噪声Q的依赖。In the embodiment of the present disclosure, the prior adjustment amount corresponding to the first moment can be calculated based on the operation data corresponding to the first moment and the second moment, so as to calculate the first prior estimated covariance corresponding to the current moment k based on the prior adjustment amount, thereby getting rid of the dependence on noise Q in the calculation process.

在一种可选的实施方式中,基于所述第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到当前时刻的最优估计,包括:In an optional implementation, filtering the first a priori estimation data at the current moment based on the first a priori estimation covariance to obtain the optimal estimate at the current moment includes:

获取目标车辆的后验信息,并根据后验信息确定所述第一时刻的运动状态数据对应的后验状态数据;Acquire a posteriori information of the target vehicle, and determine the a posteriori state data corresponding to the motion state data at the first moment according to the a posteriori information;

根据后验状态数据生成所述当前时刻的第一先验估计数据;Generate first a priori estimation data at the current moment according to the a priori state data;

根据第一先验估计协方差确定所述当前时刻对应的卡尔曼滤波矩阵,并根据所述卡尔曼滤波矩阵对所述第一先验估计数据进行滤波,得到当前时刻的最优估计。The Kalman filter matrix corresponding to the current moment is determined according to the first a priori estimate covariance, and the first a priori estimate data is filtered according to the Kalman filter matrix to obtain the optimal estimate at the current moment.

在本公开实施例中,可以根据第一时刻对应的后验状态数据生成当前时刻的第一先验估计数据,以根据当前时刻对应的卡尔曼滤波矩阵对该第一先验估计数据进行滤波,从而尽可能的滤除噪声对第一先验估计数据的影响。In the embodiment of the present disclosure, first a priori estimation data at the current moment can be generated based on the a priori state data corresponding to the first moment, and the first prior estimation data can be filtered according to the Kalman filter matrix corresponding to the current moment, so as to filter out the influence of noise on the first prior estimation data as much as possible.

在一种可选的实施方式中,根据后验状态数据生成当前时刻的第一先验估计数据,包括:In an optional implementation, generating first a priori estimation data at the current moment according to the a priori state data includes:

获取目标车辆对应的目标状态动态模型,其中,目标状态动态模型用于预测目标车辆在各个时刻的运动状态数据;Obtaining a target state dynamic model corresponding to the target vehicle, wherein the target state dynamic model is used to predict the motion state data of the target vehicle at each moment;

根据目标状态动态模型对所述第一时刻的运动状态数据进行预测,得到预测结果;Predicting the motion state data at the first moment according to the target state dynamic model to obtain a prediction result;

基于后验状态数据对预测结果进行调整,得到当前时刻的第一先验估计数据。The prediction result is adjusted based on the posterior state data to obtain the first a priori estimation data at the current moment.

在本公开实施例中,可以基于第一时刻的后验状态数据对第一时刻运动状态数据的预测结果进行调整,以得到当前时刻的第一先验估计数据,从而为后续基于该第一先验估计数据计算当前时刻的最优估计提供了技术基础。In the embodiment of the present disclosure, the prediction result of the motion state data at the first moment can be adjusted based on the posterior state data at the first moment to obtain the first a priori estimation data at the current moment, thereby providing a technical basis for subsequently calculating the optimal estimate of the current moment based on the first a priori estimation data.

在一种可选的实施方式中,根据所述卡尔曼滤波矩阵对所述第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,包括:In an optional implementation, filtering the first a priori estimation data according to the Kalman filter matrix to obtain an optimal estimation of the motion state data at the current moment includes:

获取目标车辆在当前时刻的角度信息,并计算角度信息与所述第一先验估计数据的残差;Obtaining angle information of the target vehicle at the current moment, and calculating the residual between the angle information and the first a priori estimation data;

基于卡尔曼滤波矩阵与所述残差对第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计。The first a priori estimation data is filtered based on the Kalman filter matrix and the residual to obtain an optimal estimation of the motion state data at the current moment.

在一种可选的实施方式中,最优估计包括:当前时刻对所述目标车辆的角速度信息以及角度信息的估计值。In an optional implementation, the optimal estimation includes: estimated values of angular velocity information and angle information of the target vehicle at the current moment.

在本公开实施例中,可以计算当前时刻的角度信息与第一先验估计数据的残差,以基于该残差对第一先验估计数据进行修正,以得到当前时刻对应的最优估计,从而提高了该最优估计的准确性。In the disclosed embodiment, the residual between the angle information at the current moment and the first a priori estimation data can be calculated, and the first a priori estimation data can be corrected based on the residual to obtain the optimal estimate corresponding to the current moment, thereby improving the accuracy of the optimal estimate.

第二方面,本发明提供了一种数据滤波装置,该装置包括:In a second aspect, the present invention provides a data filtering device, the device comprising:

获取模块,获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据;An acquisition module is used to acquire motion state data of the target vehicle during driving, wherein the motion state data includes motion state data at a first moment and a second moment;

确定模块,用于基于运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差;A determination module, used to determine a first a priori estimation covariance at a current moment based on the motion state data, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimation covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori;

滤波模块,用于基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。A filtering module is used to filter the first prior estimation data at the current moment based on the first prior estimation covariance to obtain the optimal estimate of the motion state data at the current moment, wherein the first prior estimation data is used to indicate the prior estimation of the motion state data at the current moment.

第三方面,本发明提供了一种计算机设备,包括:存储器和处理器,存储器和处理器之间互相通信连接,存储器中存储有计算机指令,处理器通过执行计算机指令,从而执行上述第一方面或其对应的任一实施方式的数据滤波方法。In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the data filtering method of the first aspect or any corresponding embodiment thereof by executing the computer instructions.

第四方面,本发明提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机指令,计算机指令用于使计算机执行上述第一方面或其对应的任一实施方式的数据滤波方法。In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the data filtering method of the first aspect or any corresponding embodiment thereof.

第五方面,本发明提供了一种计算机程序产品,包括计算机指令,计算机指令用于使计算机执行上述第一方面或其对应的任一实施方式的数据滤波方法。In a fifth aspect, the present invention provides a computer program product, comprising computer instructions for causing a computer to execute the data filtering method of the first aspect or any corresponding embodiment thereof.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

图1是根据本发明实施例的一种数据滤波方法的流程图;FIG1 is a flow chart of a data filtering method according to an embodiment of the present invention;

图2是根据本发明实施例的另一种数据滤波方法的流程示意图;FIG2 is a schematic flow chart of another data filtering method according to an embodiment of the present invention;

图3是根据本发明实施例的又一种数据滤波方法的流程示意图;FIG3 is a schematic flow chart of another data filtering method according to an embodiment of the present invention;

图4是根据本发明实施例的再一种数据滤波方法的卡尔曼滤波器的计算框图;FIG4 is a calculation block diagram of a Kalman filter of another data filtering method according to an embodiment of the present invention;

图5是根据本发明实施例的数据滤波装置的结构框图;FIG5 is a structural block diagram of a data filtering device according to an embodiment of the present invention;

图6是本发明实施例的计算机设备的硬件结构示意图。FIG. 6 is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention.

具体实施方式Detailed ways

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

结合数据滤波方法的执行所依赖的应用场景,在此对应用场景进行描述。In conjunction with the application scenarios on which the execution of the data filtering method depends, the application scenarios are described here.

在新能源汽车开发的过程中,采用的电机通常为永磁同步电机,从而在降低车辆能耗的同时,延长了电池的续航里程。然而,考虑到永磁同步电机驱动器的磁场定向矢量控制方法会在同步参考系中产生正负六次谐波电流,导致了扭矩波动和控制性能恶化。In the process of developing new energy vehicles, the motors used are usually permanent magnet synchronous motors, which reduce vehicle energy consumption and extend the battery life. However, considering that the field-oriented vector control method of the permanent magnet synchronous motor drive will generate positive and negative sixth harmonic currents in the synchronous reference system, it leads to torque fluctuations and deterioration of control performance.

基于此,在新能源汽车领域通常使用卡尔曼滤波器对谐波电流进行滤除,通常来说,卡尔曼滤波器的滤波效果依赖于针对过程噪声Q的预测,然而,在汽车行驶过程中,由于无法预测行驶路况,因此无法得到准确的过程噪声Q,从而影响了卡尔曼滤波器的预测精度。Based on this, Kalman filters are usually used in the field of new energy vehicles to filter out harmonic currents. Generally speaking, the filtering effect of the Kalman filter depends on the prediction of the process noise Q. However, during the driving of the car, the road conditions cannot be predicted, so the accurate process noise Q cannot be obtained, which affects the prediction accuracy of the Kalman filter.

举例来说,在相关的卡尔曼滤波器中的算法如下:For example, the algorithm in the relevant Kalman filter is as follows:

卡尔曼滤波器的时间更新方程如下:The time update equation of the Kalman filter is as follows:

卡尔曼滤波器的状态更新方程如下:The state update equation of the Kalman filter is as follows:

其中,和/>分别表示k-1时刻和k时刻的后验状态估计值,是滤波的结果之一,即更新后的结果,也叫最优估计(估计的状态,理论上不可能实时测得每时刻状态的确切结果,所以叫估计)。in, and/> They represent the posterior state estimates at time k-1 and time k respectively, which are one of the results of filtering, that is, the updated result, also called the optimal estimate (the estimated state, it is theoretically impossible to measure the exact result of the state at each moment in real time, so it is called an estimate).

表示k时刻的先验状态估计值,是滤波的中间计算结果,即根据上一时刻(k-1时刻)的最优估计预测的k时刻的结果,是预测方程的结果。 It represents the prior state estimate at time k, which is the intermediate calculation result of filtering, that is, the result of time k predicted based on the optimal estimate at the previous time (time k-1), and is the result of the prediction equation.

PK表示k时刻的后验估计协方差(即的协方差,表示状态的不确定度),是滤波的结果之一。P K represents the posterior estimated covariance at time k (i.e. The covariance of , representing the uncertainty of the state), is one of the results of filtering.

H表示状态变量到测量(观测)的转换矩阵,表示将状态和观测连接起来的关系,卡尔曼滤波里为线性关系,它负责将m维的测量值转换到n维,使之符合状态变量的数学形式,是滤波的前提条件之一。H represents the conversion matrix from state variables to measurements (observations), which represents the relationship connecting the state and observations. In Kalman filtering, it is a linear relationship. It is responsible for converting m-dimensional measurements to n-dimensional ones to conform to the mathematical form of state variables, and is one of the prerequisites for filtering.

zk表示测量值(观测值),是滤波的输入。z k represents the measurement value (observation value) and is the input of the filter.

Kk表示滤波增益矩阵,是滤波的中间计算结果,卡尔曼增益,或卡尔曼系数。K k represents the filter gain matrix, which is the intermediate calculation result of the filter, the Kalman gain, or the Kalman coefficient.

A表示状态转移矩阵,实际上是对目标状态转换的一种猜想模型。例如在机动目标跟踪中,状态转移矩阵常常用来对目标的运动建模,其模型可能为匀速直线运动或者匀加速运动。当状态转移矩阵不符合目标的状态转换模型时,滤波会很快发散。A represents the state transfer matrix, which is actually a conjecture model of the target state transition. For example, in maneuvering target tracking, the state transfer matrix is often used to model the target's motion, which may be uniform linear motion or uniformly accelerated motion. When the state transfer matrix does not conform to the target's state transition model, the filter will diverge quickly.

Q表示过程激励噪声协方差(系统过程的协方差)。该参数被用来表示状态转换矩阵与实际过程之间的误差。因为我们无法直接观测到过程信号,所以Q的取值是很难确定的。是卡尔曼滤波器用于估计离散时间过程的状态变量,也叫预测模型本身带来的噪声。Q represents the process excitation noise covariance (covariance of the system process). This parameter is used to represent the error between the state transition matrix and the actual process. Because we cannot directly observe the process signal, the value of Q is difficult to determine. It is the state variable used by the Kalman filter to estimate the discrete time process, also known as the noise brought by the prediction model itself.

R表示测量噪声协方差。滤波器实际实现时,测量噪声协方差R一般可以观测得到,是滤波器的已知条件。R represents the measurement noise covariance. When the filter is actually implemented, the measurement noise covariance R can generally be observed and is a known condition of the filter.

B表示将输入转换为状态的矩阵。B represents the matrix that transforms input into state.

表示实际观测和预测观测的残差,和卡尔曼增益一起修正先验(预测),得到后验。 It represents the residual between the actual observation and the predicted observation, and together with the Kalman gain, it corrects the prior (prediction) to obtain the posterior.

本发明实施例提供了一种数据滤波方法,首先可以获取目标车辆在行驶中的运动状态数据,其中,该运动状态数据包括第一时刻与第二时刻的运动状态数据,然后,可以基于该运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。接下来,可以基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,该第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。基于此,本公开可以无需对过程噪声Q进行预测,而是通过第一时刻与第二时刻的运动状态数据来确定当前时刻的第一先验估计协方差,从而在摆脱了对精确Q依赖的基础上提高了对当前时刻运动状态数据的最优估计的精度。The embodiment of the present invention provides a data filtering method, which can first obtain the motion state data of the target vehicle during driving, wherein the motion state data includes the motion state data at the first moment and the second moment, and then, based on the motion state data, determine the first a priori estimated covariance at the current moment, wherein the current moment is the moment after the first moment and the second moment, and the first a priori estimated covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori. Next, the first a priori estimated data at the current moment can be filtered based on the first a priori estimated covariance to obtain the optimal estimate of the motion state data at the current moment, wherein the first a priori estimated data is used to indicate the a priori estimate of the motion state data at the current moment. Based on this, the present disclosure can determine the first a priori estimated covariance of the current moment through the motion state data at the first moment and the second moment without predicting the process noise Q, thereby improving the accuracy of the optimal estimate of the motion state data at the current moment on the basis of getting rid of the dependence on accurate Q.

根据本发明实施例,提供了一种数据滤波方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。According to an embodiment of the present invention, an embodiment of a data filtering method is provided. It should be noted that the steps shown in the flowchart 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 flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

在本实施例中提供了一种数据滤波方法,可用于上述新能源汽车的驱动控制系统,具体的,可应用于上述卡尔曼滤波器,图1是根据本发明实施例的一种数据滤波方法的流程图,如图1所示,该流程包括如下步骤:In this embodiment, a data filtering method is provided, which can be used in the driving control system of the new energy vehicle. Specifically, it can be applied to the Kalman filter. FIG. 1 is a flow chart of a data filtering method according to an embodiment of the present invention. As shown in FIG. 1 , the process includes the following steps:

步骤S101,获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据。Step S101, obtaining motion state data of a target vehicle during driving, wherein the motion state data includes motion state data at a first moment and a second moment.

在本公开实施例中,上述运动状态数据可以包括目标车辆在运行过程中的角速度ωe以及角度θ,该运动状态数据x可以表示为一个状态向量其中,在确定角速度ωe时,可以获取目标车辆的行驶速度v以及转弯半径r,这里,ωe=v/r。同时,可以基于目标车辆的车轮位置安装的角度传感器获取角度θ。In the embodiment of the present disclosure, the above-mentioned motion state data may include the angular velocity ω e and the angle θ of the target vehicle during operation. The motion state data x can be expressed as a state vector When the angular velocity ω e is determined, the travel speed v and the turning radius r of the target vehicle can be obtained, where ω e =v/r. At the same time, the angle θ can be obtained based on the angle sensor installed at the wheel position of the target vehicle.

在获取上述运动状态数据时可以分别获取第一时刻与第二时刻的运动状态数据,应理解的是,本公开中的数据滤波过程是在目标车辆的行驶过程中实时进行的,相对于当前时刻k来说,第一时刻k-1是k的前一时刻,第二时刻k-2是k-1的前一时刻。具体的,时刻之间的间隔可以是秒,毫秒等,具体以实际使用需求为准,本公开对此不做限定。When acquiring the above motion state data, the motion state data of the first moment and the second moment can be acquired respectively. It should be understood that the data filtering process in the present disclosure is performed in real time during the driving process of the target vehicle. Relative to the current moment k, the first moment k-1 is the moment before k, and the second moment k-2 is the moment before k-1. Specifically, the interval between moments can be seconds, milliseconds, etc., which is subject to actual use requirements, and the present disclosure does not limit this.

步骤S102,基于上述运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为所述第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。Step S102, based on the above motion state data, determine the first a priori estimated covariance of the current moment, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimated covariance is used to indicate the covariance of the error when performing a priori estimation on the motion state data at the current moment.

在本公开实施例中,应理解的是,在卡尔曼滤波算法分为两步:状态预测和状态更新,其中,状态预测用于对当前时刻k的第一先验估计协方差Pk|k-1进行预测,状态更新用于根据目标车辆在前一时刻的运动状态来预测当前时刻的运动状态,具体包括:基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻k运动状态数据的最优估计。具体确定第一先验估计协方差的方式如下所述,此处不再赘述。In the disclosed embodiment, it should be understood that the Kalman filter algorithm is divided into two steps: state prediction and state update, wherein the state prediction is used to predict the first a priori estimated covariance P k|k-1 at the current moment k, and the state update is used to predict the motion state at the current moment based on the motion state of the target vehicle at the previous moment, specifically including: filtering the first a priori estimated data at the current moment based on the first a priori estimated covariance to obtain the optimal estimate of the motion state data at the current moment k. The specific method for determining the first a priori estimated covariance is as follows and will not be repeated here.

步骤S103,基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。Step S103, filtering the first a priori estimation data at the current moment based on the first a priori estimation covariance to obtain an optimal estimate of the motion state data at the current moment, wherein the first a priori estimation data is used to indicate a priori estimation of the motion state data at the current moment.

在本公开实施例中,由上可知,在卡尔曼滤波算法中包括状态预测和状态更新,其中,状态更新用于根据目标车辆在前一时刻的运动状态来预测当前时刻的运动状态。在进行状态更新时包括:基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻k运动状态数据的最优估计。In the disclosed embodiment, it can be seen from the above that the Kalman filter algorithm includes state prediction and state update, wherein the state update is used to predict the motion state of the target vehicle at the current moment according to the motion state of the target vehicle at the previous moment. When performing the state update, it includes: filtering the first a priori estimation data at the current moment based on the first a priori estimation covariance to obtain the optimal estimate of the k motion state data at the current moment.

具体的,在进行上述状态更新时,首先可以确定当前时刻k对应的第一先验估计数据应理解的是,该第一先验估计数据往往存在一定的偏差,因此,可以通过当前时刻的第一先验估计协方差Pk|k-1对该第一先验估计数据进行滤波,从而得到当前时刻运动状态数据的最优估计/>应理解的是,该/>既可以作为卡尔曼滤波器输出的当前时刻的最优估计,同时也可作为当前时刻的后验状态数据输入到卡尔曼滤波器,参与到针对下一时刻k+1的最优估计的计算过程中。Specifically, when performing the above state update, the first a priori estimated data corresponding to the current time k can be determined. It should be understood that the first a priori estimation data often has a certain deviation. Therefore, the first a priori estimation data can be filtered by the first a priori estimation covariance P k|k-1 at the current moment, so as to obtain the optimal estimation of the motion state data at the current moment. It should be understood that the /> It can be used as the optimal estimate of the current moment output by the Kalman filter, and can also be input into the Kalman filter as the posterior state data of the current moment, participating in the calculation process of the optimal estimate for the next moment k+1.

通过上述描述可知,在本公开实施例中,首先可以获取目标车辆在行驶中的运动状态数据,其中,该运动状态数据包括第一时刻与第二时刻的运动状态数据,然后,可以基于该运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。接下来,可以基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,该第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。基于此,本公开可以无需对过程噪声Q进行预测,而是通过第一时刻与第二时刻的运动状态数据来确定当前时刻的第一先验估计协方差,从而在摆脱了对精确Q依赖的基础上提高了对当前时刻运动状态数据的最优估计的准确性。It can be known from the above description that in the embodiment of the present disclosure, first, the motion state data of the target vehicle during driving can be obtained, wherein the motion state data includes the motion state data of the first moment and the second moment, and then, based on the motion state data, the first a priori estimation covariance of the current moment can be determined, wherein the current moment is the moment after the first moment and the second moment, and the first a priori estimation covariance is used to indicate the covariance of the error when the motion state data of the current moment is estimated a priori. Next, the first a priori estimation data of the current moment can be filtered based on the first a priori estimation covariance to obtain the optimal estimate of the motion state data of the current moment, wherein the first a priori estimation data is used to indicate the a priori estimate of the motion state data of the current moment. Based on this, the present disclosure can determine the first a priori estimation covariance of the current moment through the motion state data of the first moment and the second moment without predicting the process noise Q, thereby improving the accuracy of the optimal estimate of the motion state data of the current moment on the basis of getting rid of the dependence on the precise Q.

在本实施例中提供了另一种数据滤波方法,可用于上述新能源汽车的驱动控制系统,具体的,可应用于上述卡尔曼滤波器,图2是根据本发明实施例的数据滤波方法的流程图,如图2所示,该流程包括如下步骤:In this embodiment, another data filtering method is provided, which can be used for the driving control system of the new energy vehicle. Specifically, it can be applied to the Kalman filter. FIG. 2 is a flow chart of the data filtering method according to an embodiment of the present invention. As shown in FIG. 2, the process includes the following steps:

步骤S201,获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据。详细请参见图1所示实施例的步骤S101,在此不再赘述。Step S201, obtaining the motion state data of the target vehicle during driving, wherein the motion state data includes the motion state data at the first moment and the second moment. Please refer to step S101 of the embodiment shown in FIG1 for details, which will not be repeated here.

步骤S202,基于上述运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。Step S202, based on the above motion state data, determine the first a priori estimated covariance of the current moment, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimated covariance is used to indicate the covariance of the error when performing a priori estimation on the motion state data at the current moment.

具体地,上述步骤S202包括:Specifically, the above step S202 includes:

步骤S2021,基于第一时刻与第二时刻的运动状态数据,确定第一时刻对应的先验调整量,其中,先验调整量用于指示对第一时刻的第二先验估计协方差的误差进行修正时的调整量。Step S2021, based on the motion state data at the first moment and the second moment, determine the a priori adjustment amount corresponding to the first moment, wherein the a priori adjustment amount is used to indicate the adjustment amount when correcting the error of the second a priori estimate covariance at the first moment.

步骤S2022,根据先验调整量与上述第二先验估计协方差的和确定当前时刻的第一先验估计协方差。Step S2022: determining the first a priori estimated covariance at the current moment according to the sum of the a priori adjustment amount and the second a priori estimated covariance.

在本公开实施例中,上述先验调整量可以记为ΔPk-1,当前时刻的第一先验估计协方差可以记为Pk|k-1。具体的,In the embodiment of the present disclosure, the a priori adjustment amount can be recorded as ΔP k-1 , and the first a priori estimated covariance at the current moment can be recorded as P k|k-1 . Specifically,

Pk|k-1=Pk-1|k-2+ΔPk-1 P k|k-1 =P k-1|k-2 +ΔP k-1

其中,Pk-1|k-2为第一时刻k-1对应的第二先验估计协方差。Wherein, P k-1|k-2 is the second prior estimation covariance corresponding to the first moment k-1.

这里,在确定上述第二先验估计协方差Pk-1|k-2时,可以基于第二时刻对应的先验调整量来确定,该先验调整量可以基于第二时刻的运动状态数据对应的后验状态数据以及第三先验估计数据/>来确定,具体确定第二时刻的先验调整量的方式如下述确定第一时刻对应的先验调整量的实施例中所述,此处不再赘述。Here, when determining the second a priori estimated covariance P k-1|k-2 , it can be determined based on the a priori adjustment amount corresponding to the second moment, and the a priori adjustment amount can be determined based on the a posteriori state data corresponding to the motion state data at the second moment And the third prior estimation data/> The specific method of determining the prior adjustment amount at the second moment is as described in the following embodiment of determining the prior adjustment amount corresponding to the first moment, and will not be repeated here.

步骤S203,基于上述第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,第一先验估计数据用于指示对所述当前时刻运动状态数据的先验估计。详细请参见图1所示实施例的步骤S101,在此不再赘述。Step S203, based on the first a priori estimation covariance, the first a priori estimation data at the current moment is filtered to obtain the optimal estimation of the motion state data at the current moment, wherein the first a priori estimation data is used to indicate the a priori estimation of the motion state data at the current moment. For details, please refer to step S101 of the embodiment shown in FIG1 , which will not be described in detail here.

在本公开实施例中,由上可知,在相关的卡尔曼滤波算法中需要通过噪声Q来计算当前时刻k对应的第一先验估计协方差,而在本公开中,计算第一先验估计协方差时舍弃了噪声Q,从而提高了确定出的第一先验估计协方差的精度,进而提高了对当前时刻运动状态数据的最优估计的准确性。In the embodiment of the present disclosure, it can be seen from the above that in the relevant Kalman filter algorithm, noise Q is needed to calculate the first prior estimate covariance corresponding to the current moment k, while in the present disclosure, noise Q is discarded when calculating the first prior estimate covariance, thereby improving the accuracy of the determined first prior estimate covariance, and further improving the accuracy of the optimal estimate of the motion state data at the current moment.

在一些可选的实施方式中,第一时刻的运动状态数据包括第一时刻对应的后验状态数据,其中,该后验状态数据用于指示第一时刻对应的最优估计,上述步骤S2021包括:In some optional implementations, the motion state data at the first moment includes a posteriori state data corresponding to the first moment, wherein the a posteriori state data is used to indicate the optimal estimate corresponding to the first moment, and the above step S2021 includes:

步骤a1,获取预先基于第二时刻的运动状态数据确定出的第一时刻的第二先验估计协方差。Step a1, obtaining a second a priori estimated covariance at the first moment that is determined in advance based on the motion state data at the second moment.

步骤a2,基于所述后验状态数据与第二先验估计数据的差值,确定第一时刻对应的后验残差。Step a2: determining the a posteriori residual corresponding to the first moment based on the difference between the a posteriori state data and the second a priori estimated data.

步骤a3,基于所述后验残差对第二先验估计协方差进行修正,得到第一时刻对应的先验调整量。Step a3: correcting the second a priori estimated covariance based on the a priori residual to obtain a priori adjustment amount corresponding to the first moment.

在本公开实施例中,上述第一时刻的运动状态数据中的后验状态数据可以记为上述第二先验估计数据可以记为/>基于此,第一时刻对应的后验残差 In the embodiment of the present disclosure, the a posteriori state data in the motion state data at the first moment can be recorded as The second priori estimation data can be recorded as/> Based on this, the posterior residual corresponding to the first moment

接下里,可以根据上述后验参数Δxk-1对第二先验估计协方差进行修正,得到第一时刻对应的先验调整量ΔPk-1,具体的:Next, the second a priori estimated covariance can be corrected according to the above a posteriori parameter Δx k-1 to obtain the a priori adjustment amount ΔP k-1 corresponding to the first moment, specifically:

其中,为转置矩阵,Kk-1为k-1时刻对应的卡尔曼滤波矩阵,Pk-1|k-2为第一时刻k-1对应的先验估计协方差,H为状态变量到测量的转换矩阵,表示将状态和观测连接起来的关系,负责将m维的测量值转换到n维,使之符合状态变量的数学形式,为卡尔曼滤波器中已知的滤波前提条件之一。in, is the transposed matrix, K k-1 is the Kalman filter matrix corresponding to the k-1 moment, P k-1|k-2 is the prior estimated covariance corresponding to the first moment k-1, and H is the state variable to measurement conversion matrix, which represents the relationship connecting the state and observation. It is responsible for converting the m-dimensional measurement value to the n-dimensional one so that it conforms to the mathematical form of the state variable. It is one of the known filtering prerequisites in the Kalman filter.

具体的,确定上述Kk-1的方式如下述确定当前时刻k对应的卡尔曼滤波矩阵的实施例中所述,此处不再赘述。应理解的是,上述Pk-1|k-2为根据第二时刻k-2的运动状态数据对应的后验状态数据以及第三先验估计数据/>确定的。因此,在当前时刻的前序步骤中计算出该/>与/>后,可以将该/>与/>预存到缓存空间中,以便在后续运算中直接进行读取。Specifically, the method for determining the above K k-1 is as described in the embodiment of determining the Kalman filter matrix corresponding to the current time k below, and will not be repeated here. It should be understood that the above P k-1|k-2 is the posterior state data corresponding to the motion state data at the second time k-2. And the third prior estimation data/> Therefore, in the previous step at the current moment, the / > With/> After that, you can put the /> With/> Pre-stored in the cache space so that it can be directly read in subsequent operations.

在本公开实施例中,可以基于第一时刻以及第二时刻对应的运算数据来计算第一时刻对应的先验调整量,以基于该先验调整量来计算当前时刻k对应的第一先验估计协方差,从而摆脱计算过程中对噪声Q的依赖。In the embodiment of the present disclosure, the prior adjustment amount corresponding to the first moment can be calculated based on the operation data corresponding to the first moment and the second moment, so as to calculate the first prior estimated covariance corresponding to the current moment k based on the prior adjustment amount, thereby getting rid of the dependence on noise Q in the calculation process.

在本实施例中提供了又一种数据滤波方法,可用于上述新能源汽车的驱动控制系统,具体的,可应用于上述卡尔曼滤波器,图3是根据本发明实施例的数据滤波方法的流程图,如图3所示,该流程包括如下步骤:In this embodiment, another data filtering method is provided, which can be used for the driving control system of the new energy vehicle. Specifically, it can be applied to the Kalman filter. FIG3 is a flow chart of the data filtering method according to an embodiment of the present invention. As shown in FIG3, the process includes the following steps:

步骤S301,获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据。详细请参见图1所示实施例的步骤S101,在此不再赘述。Step S301, obtaining the motion state data of the target vehicle during driving, wherein the motion state data includes the motion state data at the first moment and the second moment. Please refer to step S101 of the embodiment shown in FIG1 for details, which will not be repeated here.

步骤S302,基于上述运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。详细请参见图1所示实施例的步骤S102,在此不再赘述。Step S302, based on the above motion state data, determine the first a priori estimated covariance of the current moment, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimated covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori. For details, please refer to step S102 of the embodiment shown in Figure 1, which will not be repeated here.

步骤S303,基于上述第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到当前时刻的最优估计。Step S303: filtering the first a priori estimation data at the current moment based on the first a priori estimation covariance to obtain the optimal estimation at the current moment.

具体地,上述步骤S303包括:Specifically, the above step S303 includes:

步骤S3031,获取目标车辆的后验信息,并根据该后验信息确定第一时刻的运动状态数据对应的后验状态数据。Step S3031, obtaining a posteriori information of the target vehicle, and determining the a posteriori state data corresponding to the motion state data at the first moment according to the a posteriori information.

步骤S3032,根据上述后验状态数据生成当前时刻的第一先验估计数据。Step S3032, generating first a priori estimation data at the current moment according to the a priori state data.

步骤S3033,根据上述第一先验估计协方差确定当前时刻对应的卡尔曼滤波矩阵,并根据该卡尔曼滤波矩阵对第一先验估计数据进行滤波,得到当前时刻的最优估计。Step S3033, determining the Kalman filter matrix corresponding to the current moment according to the first a priori estimate covariance, and filtering the first a priori estimate data according to the Kalman filter matrix to obtain the optimal estimate at the current moment.

在本公开实施例中,上述后验信息中可以为卡尔曼滤波器针对当前时刻之前的各个时刻进行预测时输出的最优估计,其中,每个时刻的最优估计都可以作为计算之后时刻的最优估计的后验状态数据。In the disclosed embodiment, the above-mentioned posterior information may be the optimal estimate output by the Kalman filter when predicting each moment before the current moment, wherein the optimal estimate at each moment may be used as the posterior state data for calculating the optimal estimate at subsequent moments.

在确定出上述第一先验估计数据之后,可以基于该/>确定上述当前时刻的卡尔曼滤波矩阵Kk,具体的:After determining the first a priori estimation data Afterwards, based on this Determine the Kalman filter matrix K k at the current moment, specifically:

Kk=Pk|k-1H[Pk|k-1HT+R]-1 K k = P k|k-1 H[P k|k-1 H T + R] -1

其中,H如上所述,此处不再赘述,R为测量噪声协方差,在卡尔曼滤波器实际实现时,测量噪声协方差R一般可以观测得到,是滤波器的已知条件。Wherein, H is as described above and will not be described again here. R is the measurement noise covariance. When the Kalman filter is actually implemented, the measurement noise covariance R can generally be observed and is a known condition of the filter.

在本公开实施例中,可以根据第一时刻对应的后验状态数据生成当前时刻的第一先验估计数据,以根据当前时刻对应的卡尔曼滤波矩阵对该第一先验估计数据进行滤波,从而尽可能的滤除噪声对第一先验估计数据的影响。In the embodiment of the present disclosure, first a priori estimation data at the current moment can be generated based on the a priori state data corresponding to the first moment, and the first prior estimation data can be filtered according to the Kalman filter matrix corresponding to the current moment, so as to filter out the influence of noise on the first prior estimation data as much as possible.

在一些可选的实施方式中,上述步骤S3032包括:In some optional implementations, the above step S3032 includes:

步骤b1,获取所述目标车辆对应的目标状态动态模型,其中,所述目标状态动态模型用于预测所述目标车辆在各个时刻的运动状态数据。Step b1, obtaining a target state dynamic model corresponding to the target vehicle, wherein the target state dynamic model is used to predict the motion state data of the target vehicle at each moment.

步骤b2,根据所述目标状态动态模型对所述第一时刻的运动状态数据进行预测,得到预测结果。Step b2: predicting the motion state data at the first moment according to the target state dynamic model to obtain a prediction result.

步骤b3,基于所述后验状态数据对所述预测结果进行调整,得到所述当前时刻的第一先验估计数据。Step b3, adjusting the prediction result based on the posterior state data to obtain the first a priori estimation data at the current moment.

在本公开实施例中,可以将上述目标状态动态模型记为f(x),该目标状态动态模型用于对第一时刻k-1的运动状态数据进行预测,得到预测结果然后,可以基于第一时刻的后验状态数据Δxk-1对/>进行调整,以得到第一先验估计数据/>具体的:In the embodiment of the present disclosure, the target state dynamic model can be recorded as f(x), and the target state dynamic model is used to predict the motion state data at the first time k-1 to obtain the prediction result Then, based on the a posteriori state data Δx k-1 at the first moment, Adjust to obtain the first prior estimation data/> specific:

其中,T为转置矩阵,B用于将输入uk-1转换为状态的矩阵。Among them, T is the transposed matrix, and B is used to convert the input u k-1 into the state matrix.

在本公开实施例中,可以基于第一时刻的后验状态数据对第一时刻运动状态数据的预测结果进行调整,以得到当前时刻的第一先验估计数据,从而为后续基于该第一先验估计数据计算当前时刻的最优估计提供了技术基础。In the embodiment of the present disclosure, the prediction result of the motion state data at the first moment can be adjusted based on the posterior state data at the first moment to obtain the first a priori estimation data at the current moment, thereby providing a technical basis for subsequently calculating the optimal estimate of the current moment based on the first a priori estimation data.

在一些可选的实施方式中,上述步骤S3033包括:In some optional implementations, the above step S3033 includes:

步骤c1,获取目标车辆在当前时刻的角度信息,并计算角度信息与第一先验估计数据的残差。Step c1, obtaining the angle information of the target vehicle at the current moment, and calculating the residual between the angle information and the first prior estimation data.

步骤c2,基于卡尔曼滤波矩阵与所述残差对第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计。Step c2, filtering the first priori estimation data based on the Kalman filter matrix and the residual to obtain an optimal estimation of the motion state data at the current moment.

在本公开实施例中,角度信息y可以用于指示目标车辆在运行过程中的角度,其中,y=θ,由上可知,角度信息θ可以基于目标车辆的车轮位置安装的角度传感器获取。具体计算角度信息与第一先验估计数据的残差的过程如下所述:In the embodiment of the present disclosure, the angle information y can be used to indicate the angle of the target vehicle during operation, where y=θ. As can be seen from the above, the angle information θ can be obtained based on the angle sensor installed at the wheel position of the target vehicle. The specific process of calculating the residual of the angle information and the first prior estimation data is as follows:

其中,该残差可以用于表示实际数据观测与预测数据观测的残差。The residual can be used to represent the residual between the actual data observation and the predicted data observation.

接下来,可以基于该参数与卡尔曼滤波矩阵一起修正第一先验估计数据以得到当前时刻k对应的最优估计/>具体的修正过程如下:Next, the first prior estimation data can be corrected based on this parameter together with the Kalman filter matrix To obtain the optimal estimate corresponding to the current time k/> The specific correction process is as follows:

在一些可选的实施方式中,上述最优估计包括:当前时刻对目标车辆的角速度信息以及角度信息的估计值,该估计值可以表示为一个向量/> In some optional implementations, the above optimal estimate Includes: the estimated value of the angular velocity information and angle information of the target vehicle at the current moment, which can be expressed as a vector/>

在本公开实施例中,可以计算当前时刻的角度信息与第一先验估计数据的残差,以基于该残差对第一先验估计数据进行修正,以得到当前时刻对应的最优估计,从而提高了该最优估计的准确性。In the disclosed embodiment, the residual between the angle information at the current moment and the first a priori estimation data can be calculated, and the first a priori estimation data can be corrected based on the residual to obtain the optimal estimate corresponding to the current moment, thereby improving the accuracy of the optimal estimate.

在本实施例中提供了再一种数据滤波方法,可用于上述新能源汽车的驱动控制系统,具体的,可应用于上述卡尔曼滤波器,如下所述为该卡尔曼滤波器中的算法:In this embodiment, another data filtering method is provided, which can be used in the driving control system of the new energy vehicle. Specifically, it can be applied to the Kalman filter. The algorithm in the Kalman filter is as follows:

Kk=Pk|k-1H[Pk|k-1HT+R]-1 (16)K k =P k|k-1 H[P k|k-1 H T +R] -1 (16)

具体的,上述算法对应的再一种数据滤波方法的卡尔曼滤波器的计算框图如图4所示,其中,该卡尔曼滤波器的输入数据为上述第一时刻k-1的后验状态数据目标车辆的角度信息y以及第一时刻对应的第二先验估计协方差Pk-1|k-2,输出数据为当前时刻k的最优估计/> Specifically, the calculation block diagram of the Kalman filter of another data filtering method corresponding to the above algorithm is shown in FIG4 , wherein the input data of the Kalman filter is the a posteriori state data at the first moment k-1 above. The target vehicle's angle information y and the second a priori estimated covariance P k-1|k-2 corresponding to the first moment, the output data is the optimal estimate of the current moment k/>

综上,在本公开实施例中,首先可以获取目标车辆在行驶中的运动状态数据,其中,该运动状态数据包括第一时刻与第二时刻的运动状态数据,然后,可以基于该运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差。接下来,可以基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,该第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。基于此,本公开可以无需对过程噪声Q进行预测,而是通过第一时刻与第二时刻的运动状态数据来确定当前时刻的第一先验估计协方差,从而在摆脱了对精确Q依赖的基础上提高了对当前时刻运动状态数据的最优估计的准确性。In summary, in the embodiment of the present disclosure, the motion state data of the target vehicle during driving can be first obtained, wherein the motion state data includes the motion state data at the first moment and the second moment, and then, based on the motion state data, the first a priori estimated covariance at the current moment can be determined, wherein the current moment is the moment after the first moment and the second moment, and the first a priori estimated covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori. Next, the first a priori estimated data at the current moment can be filtered based on the first a priori estimated covariance to obtain the optimal estimate of the motion state data at the current moment, wherein the first a priori estimated data is used to indicate the a priori estimate of the motion state data at the current moment. Based on this, the present disclosure does not need to predict the process noise Q, but determines the first a priori estimated covariance at the current moment through the motion state data at the first moment and the second moment, thereby improving the accuracy of the optimal estimate of the motion state data at the current moment on the basis of getting rid of the dependence on the precise Q.

在本实施例中还提供了一种数据滤波装置,该装置用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。如以下所使用的,术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。In this embodiment, a data filtering device is also provided, which is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and/or hardware of a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

本实施例提供一种数据滤波装置,如图5所示,包括:This embodiment provides a data filtering device, as shown in FIG5 , including:

获取模块501,用于获取目标车辆在行驶中的运动状态数据,其中,运动状态数据包括第一时刻与第二时刻的运动状态数据;An acquisition module 501 is used to acquire motion state data of a target vehicle during driving, wherein the motion state data includes motion state data at a first moment and a second moment;

确定模块502,用于基于运动状态数据,确定当前时刻的第一先验估计协方差,其中,当前时刻为第一时刻与第二时刻之后的时刻,第一先验估计协方差用于指示对当前时刻的运动状态数据进行先验估计时误差的协方差;A determination module 502 is used to determine a first a priori estimation covariance at a current moment based on the motion state data, wherein the current moment is a moment after the first moment and the second moment, and the first a priori estimation covariance is used to indicate the covariance of the error when the motion state data at the current moment is estimated a priori;

滤波模块503,用于基于第一先验估计协方差对当前时刻的第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计,其中,第一先验估计数据用于指示对当前时刻运动状态数据的先验估计。The filtering module 503 is used to filter the first prior estimation data at the current moment based on the first prior estimation covariance to obtain the optimal estimation of the motion state data at the current moment, wherein the first prior estimation data is used to indicate the prior estimation of the motion state data at the current moment.

在一些可选的实施方式中,确定模块502包括:In some optional implementations, the determining module 502 includes:

第一确定单元,用于基于所运动状态数据,确定第一时刻对应的先验调整量,其中,先验调整量用于指示对第一时刻的第二先验估计协方差的误差进行修正时的调整量;A first determining unit, configured to determine a priori adjustment amount corresponding to a first moment based on the motion state data, wherein the priori adjustment amount is used to indicate an adjustment amount when correcting an error of a second priori estimate covariance at the first moment;

第二确定单元,用于根据先验调整量与第二先验估计协方差的和确定当前时刻的第一先验估计协方差。The second determining unit is used to determine the first a priori estimated covariance at a current moment according to the sum of the a priori adjustment amount and the second a priori estimated covariance.

在一些可选的实施方式中,第一时刻的运动状态数据包括第一时刻对应的后验状态数据,其中,后验状态数据用于指示所述第一时刻对应的最优估计,第一确定单元包括:In some optional implementations, the motion state data at the first moment includes a posteriori state data corresponding to the first moment, wherein the a posteriori state data is used to indicate the optimal estimate corresponding to the first moment, and the first determining unit includes:

获取子单元,用于获取预先基于第二时刻的运动状态数据确定出的第一时刻的第二先验估计协方差;An acquisition subunit, used to acquire a second a priori estimated covariance at the first moment determined in advance based on the motion state data at the second moment;

第一确定子单元,用于基于后验状态数据与第二先验估计数据的差值,确定第一时刻对应的后验残差;A first determination subunit, configured to determine a posterior residual corresponding to a first moment based on a difference between the posterior state data and the second a priori estimation data;

修正子单元,用于基于后验残差对第二先验估计协方差进行修正,得到第一时刻对应的先验调整量。The correction subunit is used to correct the second a priori estimated covariance based on the posterior residual to obtain the a priori adjustment amount corresponding to the first moment.

在一些可选的实施方式中,滤波模块503包括:In some optional implementations, the filtering module 503 includes:

获取单元,用于获取目标车辆的后验信息,并根据后验信息确定第一时刻的运动状态数据对应的后验状态数据;An acquisition unit, used to acquire a posteriori information of the target vehicle, and determine the a posteriori state data corresponding to the motion state data at the first moment according to the a posteriori information;

生成单元,用于根据所述后验状态数据生成当前时刻的第一先验估计数据;A generating unit, configured to generate first a priori estimation data at a current moment according to the a priori state data;

滤波单元,用于根据第一先验估计协方差确定当前时刻对应的卡尔曼滤波矩阵,并根据卡尔曼滤波矩阵对第一先验估计数据进行滤波,得到当前时刻的最优估计。The filtering unit is used to determine the Kalman filter matrix corresponding to the current moment according to the first prior estimation covariance, and filter the first prior estimation data according to the Kalman filter matrix to obtain the optimal estimate at the current moment.

在一些可选的实施方式中,生成单元包括:In some optional implementations, the generating unit includes:

生成子单元,用于获取所述目标车辆对应的目标状态动态模型,其中,所述目标状态动态模型用于预测所述目标车辆在各个时刻的运动状态数据;A generating subunit, used for acquiring a target state dynamic model corresponding to the target vehicle, wherein the target state dynamic model is used for predicting the motion state data of the target vehicle at each moment;

预测子单元,用于根据所述目标状态动态模型对所述第一时刻的运动状态数据进行预测,得到预测结果;A prediction subunit, configured to predict the motion state data at the first moment according to the target state dynamic model to obtain a prediction result;

调整子单元,用于基于所述后验状态数据对所述预测结果进行调整,得到所述当前时刻的第一先验估计数据。The adjustment subunit is used to adjust the prediction result based on the posterior state data to obtain the first a priori estimation data at the current moment.

在一些可选的实施方式中,滤波单元包括:In some optional implementations, the filtering unit includes:

计算子单元,用于获取所述目标车辆在当前时刻的角度信息,并计算所述角度信息与所述第一先验估计数据的残差;A calculation subunit, used for obtaining the angle information of the target vehicle at the current moment, and calculating the residual between the angle information and the first a priori estimation data;

滤波子单元,用于基于所述卡尔曼滤波矩阵与所述残差对所述第一先验估计数据进行滤波,得到对当前时刻运动状态数据的最优估计。A filtering subunit is used to filter the first prior estimation data based on the Kalman filter matrix and the residual to obtain an optimal estimate of the motion state data at a current moment.

在一些可选的实施方式中,最优估计包括:当前时刻对目标车辆的角速度信息以及角度信息的估计值。In some optional implementations, the optimal estimate includes: estimated values of angular velocity information and angle information of the target vehicle at the current moment.

上述各个模块和单元的更进一步的功能描述与上述对应实施例相同,在此不再赘述。The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

本实施例中的数据滤波装置是以功能单元的形式来呈现,这里的单元是指ASIC(Application Specific Integrated Circuit,专用集成电路)电路,执行一个或多个软件或固定程序的处理器和存储器,和/或其他可以提供上述功能的器件。The data filtering device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and/or other devices that can provide the above functions.

本发明实施例还提供一种计算机设备,具有上述图5所示的数据滤波装置。An embodiment of the present invention further provides a computer device having the data filtering device shown in FIG. 5 .

请参阅图6,图6是本发明可选实施例提供的一种计算机设备的结构示意图,如图6所示,该计算机设备包括:一个或多个处理器10、存储器20,以及用于连接各部件的接口,包括高速接口和低速接口。各个部件利用不同的总线互相通信连接,并且可以被安装在公共主板上或者根据需要以其它方式安装。处理器可以对在计算机设备内执行的指令进行处理,包括存储在存储器中或者存储器上以在外部输入/输出装置(诸如,耦合至接口的显示设备)上显示GUI的图形信息的指令。在一些可选的实施方式中,若需要,可以将多个处理器和/或多条总线与多个存储器和多个存储器一起使用。同样,可以连接多个计算机设备,各个设备提供部分必要的操作(例如,作为服务器阵列、一组刀片式服务器、或者多处理器系统)。图6中以一个处理器10为例。Please refer to Figure 6, which is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention. As shown in Figure 6, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. The various components are connected to each other using different buses for communication, and can be installed on a common motherboard or installed in other ways as needed. The processor can process instructions executed in the computer device, including instructions stored in or on the memory to display graphical information of the GUI on an external input/output device (such as a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and/or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). In Figure 6, a processor 10 is taken as an example.

处理器10可以是中央处理器,网络处理器或其组合。其中,处理器10还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路,可编程逻辑器件或其组合。上述可编程逻辑器件可以是复杂可编程逻辑器件,现场可编程逻辑门阵列,通用阵列逻辑或其任意组合。The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

其中,所述存储器20存储有可由至少一个处理器10执行的指令,以使所述至少一个处理器10执行实现上述实施例示出的方法。The memory 20 stores instructions executable by at least one processor 10, so that the at least one processor 10 executes the method shown in the above embodiment.

存储器20可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储根据计算机设备的使用所创建的数据等。此外,存储器20可以包括高速随机存取存储器,还可以包括非瞬时存储器,例如至少一个磁盘存储器件、闪存器件、或其他非瞬时固态存储器件。在一些可选的实施方式中,存储器20可选包括相对于处理器10远程设置的存储器,这些远程存储器可以通过网络连接至该计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

存储器20可以包括易失性存储器,例如,随机存取存储器;存储器也可以包括非易失性存储器,例如,快闪存储器,硬盘或固态硬盘;存储器20还可以包括上述种类的存储器的组合。The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive; the memory 20 may also include a combination of the above types of memory.

该计算机设备还包括输入装置30和输出装置40。处理器10、存储器20、输入装置30和输出装置40可以通过总线或者其他方式连接,图6中以通过总线连接为例。The computer device further includes an input device 30 and an output device 40. The processor 10, the memory 20, the input device 30 and the output device 40 may be connected via a bus or other means, and FIG6 takes the connection via a bus as an example.

输入装置30可接收输入的数字或字符信息,以及产生与该计算机设备的用户设置以及功能控制有关的键信号输入,例如触摸屏、小键盘、鼠标、轨迹板、触摸板、指示杆、一个或者多个鼠标按钮、轨迹球、操纵杆等。输出装置40可以包括显示设备、辅助照明装置(例如,LED)和触觉反馈装置(例如,振动电机)等。上述显示设备包括但不限于液晶显示器,发光二极管,显示器和等离子体显示器。在一些可选的实施方式中,显示设备可以是触摸屏。The input device 30 can receive input digital or character information, and generate key signal input related to the user settings and function control of the computer device, such as a touch screen, a keypad, a mouse, a track pad, a touch pad, an indicator bar, one or more mouse buttons, a trackball, a joystick, etc. The output device 40 may include a display device, an auxiliary lighting device (e.g., an LED) and a tactile feedback device (e.g., a vibration motor), etc. The above-mentioned display device includes but is not limited to a liquid crystal display, a light emitting diode, a display and a plasma display. In some optional embodiments, the display device can be a touch screen.

本发明实施例还提供了一种计算机可读存储介质,上述根据本发明实施例的方法可在硬件、固件中实现,或者被实现为可记录在存储介质,或者被实现通过网络下载的原始存储在远程存储介质或非暂时机器可读存储介质中并将被存储在本地存储介质中的计算机代码,从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件的存储介质上的这样的软件处理。其中,存储介质可为磁碟、光盘、只读存储记忆体、随机存储记忆体、快闪存储器、硬盘或固态硬盘等;进一步地,存储介质还可以包括上述种类的存储器的组合。可以理解,计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件,当软件或计算机代码被计算机、处理器或硬件访问且执行时,实现上述实施例示出的方法。The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

本发明的一部分可被应用为计算机程序产品,例如计算机程序指令,当其被计算机执行时,通过该计算机的操作,可以调用或提供根据本发明的方法和/或技术方案。本领域技术人员应能理解,计算机程序指令在计算机可读介质中的存在形式包括但不限于源文件、可执行文件、安装包文件等,相应地,计算机程序指令被计算机执行的方式包括但不限于:该计算机直接执行该指令,或者该计算机编译该指令后再执行对应的编译后程序,或者该计算机读取并执行该指令,或者该计算机读取并安装该指令后再执行对应的安装后程序。在此,计算机可读介质可以是可供计算机访问的任意可用的计算机可读存储介质或通信介质。A part of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and/or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the existence of computer program instructions in computer-readable media includes, but is not limited to, source files, executable files, installation package files, etc., and accordingly, the way in which computer program instructions are executed by a computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium accessible to the computer.

虽然结合附图描述了本发明的实施例,但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型,这样的修改和变型均落入由所附权利要求所限定的范围之内。Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims (11)

1. A method of filtering data applied to a kalman filter, the method comprising:
acquiring motion state data of a target vehicle in running, wherein the motion state data comprise motion state data of a first moment and a second moment;
Determining a first priori estimated covariance of a current moment based on the motion state data, wherein the current moment is a moment after the first moment and the second moment, and the first priori estimated covariance is used for indicating covariance of errors when the motion state data of the current moment is estimated a priori;
And filtering the first priori estimation data of the current moment based on the first priori estimation covariance to obtain optimal estimation of the motion state data of the current moment, wherein the first priori estimation data is used for indicating the prior estimation of the motion state data of the current moment.
2. The method of claim 1, wherein the determining a first a priori estimated covariance of a current time instant based on the motion state data comprises:
determining a priori adjustment amount corresponding to a first moment based on the motion state data, wherein the priori adjustment amount is used for indicating an adjustment amount when correcting an error of a second priori estimated covariance of the first moment;
and determining a first prior estimation covariance of the current moment according to the sum of the prior adjustment quantity and the second prior estimation covariance.
3. The method of claim 2, wherein the motion state data at the first time comprises posterior state data corresponding to the first time, wherein the posterior state data is used to indicate an optimal estimate corresponding to the first time;
The determining, based on the motion state data, a priori adjustment corresponding to a first time includes:
acquiring a second prior estimated covariance of the first moment, which is determined in advance based on the motion state data of the second moment;
Determining a posterior residual corresponding to the first moment based on the difference value between the posterior state data and the second prior estimation data;
and correcting the second priori estimated covariance based on the posterior residual error to obtain the priori adjustment corresponding to the first moment.
4. A method according to any one of claims 1 to 3, wherein filtering the first prior estimate data at the current time based on the first prior estimate covariance to obtain an optimal estimate at the current time comprises:
acquiring posterior information of a target vehicle, and determining posterior state data corresponding to the motion state data at the first moment according to the posterior information;
Generating first priori estimated data of the current moment according to the posterior state data;
and determining a Kalman filter matrix corresponding to the current moment according to the first priori estimation covariance, and filtering the first priori estimation data according to the Kalman filter matrix to obtain the optimal estimation of the current moment.
5. The method of claim 4, wherein the generating the first prior estimate data for the current time from the posterior state data comprises:
Acquiring a target state dynamic model corresponding to the target vehicle, wherein the target state dynamic model is used for predicting motion state data of the target vehicle at each moment;
predicting the motion state data at the first moment according to the target state dynamic model to obtain a prediction result;
And adjusting the prediction result based on the posterior state data to obtain first priori estimated data of the current moment.
6. The method of claim 4, wherein filtering the first a priori estimate data based on the kalman filter matrix to obtain an optimal estimate of motion state data at a current time comprises:
Acquiring angle information of the target vehicle at the current moment, and calculating residual errors of the angle information and the first priori estimated data;
and filtering the first priori estimation data based on the Kalman filtering matrix and the residual error to obtain the optimal estimation of the motion state data at the current moment.
7. The method of claim 1, wherein the optimal estimation comprises: and estimating the angular velocity information and the angle information of the target vehicle at the current moment.
8. A data filtering apparatus, the apparatus comprising:
The system comprises an acquisition module, a control module and a control module, wherein the acquisition module acquires motion state data of a target vehicle in running, and the motion state data comprises motion state data of a first moment and a second moment;
The determining module is used for determining a first priori estimated covariance of the current moment based on the motion state data, wherein the current moment is a moment after the first moment and the second moment, and the first priori estimated covariance is used for indicating covariance of errors when the motion state data of the current moment is estimated a priori;
The filtering module is used for filtering the first priori estimation data of the current moment based on the first priori estimation covariance to obtain the optimal estimation of the motion state data of the current moment, wherein the first priori estimation data is used for indicating the prior estimation of the motion state data of the current moment.
9. A computer device, comprising:
a memory and a processor in communication with each other, the memory having stored therein computer instructions, the processor executing the computer instructions to perform the data filtering method of any of claims 1 to 7.
10. A computer-readable storage medium having stored thereon computer instructions for causing a computer to perform the data filtering method of any of claims 1 to 7.
11. A computer program product comprising computer instructions for causing a computer to perform the data filtering method of any one of claims 1 to 7.
CN202410328547.8A 2024-03-21 2024-03-21 Data filtering method, device, computer equipment, storage medium and program product Pending CN118114498A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410328547.8A CN118114498A (en) 2024-03-21 2024-03-21 Data filtering method, device, computer equipment, storage medium and program product

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410328547.8A CN118114498A (en) 2024-03-21 2024-03-21 Data filtering method, device, computer equipment, storage medium and program product

Publications (1)

Publication Number Publication Date
CN118114498A true CN118114498A (en) 2024-05-31

Family

ID=91208494

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410328547.8A Pending CN118114498A (en) 2024-03-21 2024-03-21 Data filtering method, device, computer equipment, storage medium and program product

Country Status (1)

Country Link
CN (1) CN118114498A (en)

Similar Documents

Publication Publication Date Title
CN108068882A (en) The method for controlling steering
CN115062571B (en) Method, system, device and computer-readable storage medium for dynamic selection of data area applied to integrated circuit devices
CN118114498A (en) Data filtering method, device, computer equipment, storage medium and program product
CN109631892B (en) Integrated navigation data processing method for IMU data interruption
CN114987468A (en) Method and device for controlling vehicle target speed, electronic equipment and storage medium
US20130080492A1 (en) Processing kalman filter
CN112801193B (en) Positioning data processing method and device, electronic equipment and medium
CN114090054A (en) Smart device upgrade method, device, storage medium and electronic device
KR20230041183A (en) Method for generating road topology information and system thereof
CN110320800B (en) Compensation method and device of control system, medium and intelligent equipment
CN118913273A (en) Underwater submarine navigation device positioning system, method and device, computer equipment and medium
CN111984047A (en) Multi-module system current sharing control method and device and terminal equipment
CN118244874A (en) Method, system, equipment and computer storage medium for managing and controlling power consumption of server
US8793526B2 (en) Firmware management in a computing system
CN114954021B (en) Method and device for acquiring driving range, electronic equipment and vehicle
CN110968925B (en) Robot, time alignment method and device, electronic equipment and storage medium
CN112388623B (en) Steering engine position control method and device, terminal equipment and medium
CN114665757A (en) Method and apparatus for calculating absolute position of infinite rotation of motor, and storage medium
CN112328963A (en) Method and device for calculating effective value of signal
CN119245619B (en) Kalman filter parameter setting method, device, electronic equipment and storage medium
CN118746913A (en) Visual laser time synchronization method, device and storage medium based on cross-correlation
JP7112885B2 (en) Information processing equipment
CN110737278A (en) Robot stabilization control method and device, storage medium and inspection robot
CN118611494A (en) Compressor control method, device and computer readable storage medium
CN118701059A (en) Friction torque compensation method and device, 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