CN106874559B - 一种用于车轮力传感器的交互式多模型滤波方法 - Google Patents
一种用于车轮力传感器的交互式多模型滤波方法 Download PDFInfo
- Publication number
- CN106874559B CN106874559B CN201710030224.0A CN201710030224A CN106874559B CN 106874559 B CN106874559 B CN 106874559B CN 201710030224 A CN201710030224 A CN 201710030224A CN 106874559 B CN106874559 B CN 106874559B
- Authority
- CN
- China
- Prior art keywords
- model
- vehicle wheel
- dynamic model
- wheel forces
- forces
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Toys (AREA)
Abstract
本发明公开一种用于车轮力传感器的交互式多模型滤波方法,包括如下步骤:1、根据车轮力信号随机性强的特点和车轮力传感器直接输出与真实车轮力之间的数学关系,建立基于Singer模型的车轮力通用动态模型;2、根据上述建立的车轮力通用动态模型,通过选取不同的模型参数,分别建立车轮力高动态模型和车轮力低动态模型;3、根据上述建立的车轮力高动态模型和车轮力低动态模型,采用交互式多模型方法,实现对于全动态范围车轮力信号的建模,并采用卡尔曼滤波器对车轮力信号进行实时滤波。本发明能够涵盖各种动态范围的车轮力,解决随机车轮力信号难以建模的问题,从而实现基于卡尔曼滤波器的车轮力信号实时滤波。
Description
技术领域
本发明属于汽车测控技术领域,特别涉及一种用于车轮力传感器的交互式多模型滤波方法。
背景技术
汽车工业在发达国家的经济发展中,起着重要支柱的作用,从某种程度上讲,汽车行业发展水平能够体现出一个国家的整体制造业水平。在我国,近年来汽车的保有量持续增长势头,良好的市场情况刺激了汽车工业的发展,同时也暴露出我国车辆自主研发和测试方面的不足。车辆运动是由地面对车轮作用力、车辆对车轮作用力以及空气作用力等共同作用的结果,其中地面轮胎附着作用是主要因素,并受路面因素、轮胎因素、车辆因素和车辆行驶工况因素等影响,并最终在车轮受力上得到体现。因此,实时检测车辆在各种状态下的车轮力是进行车辆研发和测试的基础。由于车轮力具有很强的随机性,因此对于其进行处理的难度较大。我国在这项领域比较薄弱,尤其是对车轮力传感器输出数据的滤波方面的研究基本处于起步和探索阶段。
本发明专利所应用的对象为东南大学仪器科学与工程学院自主研发的多维轮力传感器(一种测量车轮六维力的传感器ZL201210071761.7),该多维轮力传感器能够感知车辆行驶过程中地面对车轮的作用力,包括牵引力和正压力。为了得到高精度的车轮力信号,必须对车轮力传感器的输出信号进行滤波降噪。传统的低通滤波器在去除噪声的同时也会损失车轮力信号中的高频分量,而根据不同频率分量阈值进行噪声剔除的小波滤波技术,虽然能够保留车轮力中的高频信号,但小波滤波是一种离线数据处理方法,不能满足高精度车轮力实时输出的要求。
为了解决车轮力信号实时滤波的问题,将卡尔曼滤波器引入到车轮力传感器的数据处理当中。在车轮力处理中使用卡尔曼滤波器需要对车轮力信号进行建模,然而由于车轮力的随机性和不确定性,很难对其精确建模。想要设计出适用于车轮力传感器的实时滤波器,就必须攻克这一技术难关。
发明内容
发明目的:针对现有技术中存在的问题,本发明提供一种用于车轮力传感器的交互式多模型滤波方法,能够涵盖各种动态范围的车轮力,解决随机车轮力信号难以建模的问题,从而实现基于卡尔曼滤波器的车轮力信号实时滤波。
技术方案:为解决上述技术问题,本发明提供一种用于车轮力传感器的交互式多模型滤波方法,其特征在于,包括如下步骤:
步骤一:根据车轮力信号随机性强的特点和车轮力传感器直接输出与真实车轮力之间的数学关系,建立基于Singer模型的车轮力通用动态模型;
步骤二:根据步骤一中所建立的车轮力通用动态模型,通过选取不同的模型参数,分别建立车轮力高动态模型和车轮力低动态模型;
步骤三:根据步骤二中建立的车轮力高动态模型和车轮力低动态模型,采用交互式多模型方法,实现对于全动态范围车轮力信号的建模,并采用卡尔曼滤波器对车轮力信号进行实时滤波。
进一步的,所述步骤一中建立车轮力通用动态模型的具体步骤如下:
步骤1.1:定义k时刻的模型状态量为:
其中和分别为k时刻车轮纵向力、纵向力微分以及二次微分;和分别为k时刻车轮正压力、正压力微分以及二次微分;θk、和分别为k时刻车轮旋转角度、旋转角度微分以及二次微分;
步骤1.2:根据步骤1.1中定义的参数构建车轮力通用动态模型为:
XS,k=fS(T,α)XS,k-1+wS,k-1
其中fS(T,α)=diag[fS(T,αFxw)fS(T,αFzw)fS(T,αθ)],T为采样时间,Fx w为牵引力和Fz w为正压力,αFxw、αFzw和αθ分别为Fx w、Fz w和θ的变化率,wS,k-1为系统噪声;其中fS(T,αFxw)、fS(T,αFzw)和fS(T,αθ)的具体表达式为:
进一步的,所述步骤二中建立车轮力高动态模型和车轮力地动态模型的具体步骤如下:选择αFxw=1/10,αFzw=1/10,αθ=1/10,建立车轮力高动态模型;选择αFxw=1/40,αFzw=1/40,αθ=1/30,建立车轮力低动态模型。
进一步的,所述步骤三中建立全动态范围车轮力信号模型,并进行实时滤波的具体步骤如下:根据步骤二中所建立的车轮力高动态模型和车轮力低动态模型,分别对牵引力和正压力进行卡尔曼滤波,采用交互式多模型方法,实时计算上述两个模型权重,利用经过归一化权重求取两个模型滤波结果的加权和,并将该值作为轮力滤波的最终输出,这种将两个模型滤波数据相融合的方法,解决了单一模型建模精度较低的问题,改善了车轮力信号的滤波效果。
与现有技术相比,本发明的优点在于:
通过交互式多模型方法,构建能够涵盖所有车轮力动态范围的车轮力模型,并利用卡尔曼滤波器实时对车轮力传感器输出的车轮力数据进行滤波,有效地解决了随机车轮力信号难以建模的问题,实现了车轮力数据的高精度实时滤波。该方法能够应用于东南大学仪器科学与工程学院研发的车轮力传感器的数据滤波,具有很强的实用性。
附图说明
图1为本发明的总体流程图;
图2为图1中步骤三的方法流程图;
图3为实施例中真实车轮力与车轮力传感器直接输出关系示意图。
具体实施方式
下面结合附图和具体实施方式,进一步阐明本发明。
如图1所示:一种适用于车轮力传感器的交互式多模型滤波方法,包括以下步骤:
步骤一:根据车轮力信号随机性强的特点和车轮力传感器直接输出与真实车轮力之间的数学关系,建立基于Singer模型的车轮力通用动态模型;
为了清楚阐述真实车轮力和车轮力传感器输出之间的关系,分别建立真实车轮力坐标系和车轮力传感器坐标系,如图3所示。真实车轮力坐标系定义如下:选取车轮中心为坐标原点ow,坐标轴owzw竖直向上,坐标轴owxw指向车轮前进方向;车轮力传感器坐标系定义如下:选取车轮中心为坐标原点oe,坐标轴oexe与传感器牵引力敏感梁平行,坐标轴oeze与传感器正压力敏感梁平行;θ为真实车轮力坐标系与车轮力传感器坐标系之间的偏转角度,当θ=0时,真实车轮力坐标系与车轮力传感器坐标系重合。
真实车轮力(包括牵引力和正压力)和真实车轮力坐标系与车轮力传感器坐标系之间偏转角度θ的变化都可以被表述为随机过程,由于它们都是真实存在的物理量,并且是连续变化的,且前后两个数据点之间存在一个有限的增量,可以被视为Markov过程,因此采用基于Markov过程的Singer模型来构建真实车轮力模型是合理的。Singer模型离散状态空间形式如下:
式(1)中:T为采样周期;α=1/τ表示变量的变化率,它是机动时间τ的倒数,大小取决于机动的持续时间;wk-1的方差为α的函数。
将Singer模型引入到车轮力信号处理中时,综合考虑牵引力正压力和偏转角度θ,因此选取模型的状态量为:
在如式2所示的状态的基础上,基于Singer模型构建车轮力通用动态模型为:
XS,k=fS(T,α)XS,k-1+wS,k-1 (3)
其中fS(T,α)=diag[fS(T,αFxw)fS(T,αFzw)fS(T,αθ)],T为采样时间,αFxw、αFzw和αθ分别为和θ的变化率。fS(T,αFxw)、fS(T,αDzw)和fS(T,αθ)的具体表达式为:
步骤二:根据步骤一所建立的车轮力通用动态模型,通过选取不同的模型参数,分别建立车轮力高动态模型和车轮力低动态模型;
选择αFxw、αFzw和αθ,使得车轮力通用动态模型分别适应高动态车轮力和低动态车轮力。
选择αFxw=1/10,αFzw=1/10,αθ=1/10,建立车轮力高动态模型。
选择αFxw=1/40,αFzw=1/40,αθ=1/30,建立车轮力低动态模型。
步骤三:根据步骤二所建立的车轮力高动态模型和车轮力低动态模型,采用交互式多模型方法,采用交互式多模型算法,将车轮力高动态模型和车轮力低动态模型的滤波输出相互融合,其中车轮力高动态模型和车轮力低动态模型的观测方程均为:
从而实现对于全动态范围车轮力信号的建模,并采用卡尔曼滤波器对车轮力信号进行实时滤波。
通过建立的车轮力高动态模型和车轮力低动态模型,设计基于交互式多模型方法的滤波器,分别对牵引力和正压力进行卡尔曼滤波。具体做法如图2所示:首先利用k-1时刻的滤波结果和经过归一化的模型权重,分别为轮力高动态模型和车轮力低动态模型分配模型输入,接着在k时刻,轮力高动态模型和车轮力低动态模型根据各自的模型输入与k时刻的观测量分别更新各自的模型输出与模型权重,最后利用经过归一化权重求取两个模型输出的加权和并将该值作为轮力滤波的最终输出,这种将两个模型滤波数据相融合的方法,解决了单一模型建模精度较低的问题,改善了车轮力信号的滤波效果。
本发明揭示了一种适用于车轮力传感器的交互式多模型滤波方法,通过交互式多模型方法,构建能够涵盖所有车轮力动态范围的车轮力模型,并利用卡尔曼滤波器实时对车轮力传感器输出的车轮力数据进行滤波,有效地解决了随机车轮力信号难以建模的问题,实现了车轮力数据的高精度实时滤波。该方法能够应用于东南大学仪器科学与工程学院研发的车轮力传感器的数据滤波,具有很强的实用性。
以上所述仅为本发明的实施例子而已,并不用于限制本发明。凡在本发明的原则之内,所作的等同替换,均应包含在本发明的保护范围之内。本发明未作详细阐述的内容属于本专业领域技术人员公知的已有技术。
Claims (1)
1.一种用于车轮力传感器的交互式多模型滤波方法,其特征在于,包括如下步骤:
步骤一:根据车轮力信号随机性强的特点和车轮力传感器直接输出与真实车轮力之间的数学关系,建立基于Singer模型的车轮力通用动态模型;
步骤二:根据步骤一中所建立的车轮力通用动态模型,通过选取不同的模型参数,分别建立车轮力高动态模型和车轮力低动态模型;
步骤三:根据步骤二中建立的车轮力高动态模型和车轮力低动态模型,采用交互式多模型方法,实现对于全动态范围车轮力信号的建模,并采用卡尔曼滤波器对车轮力信号进行实时滤波;
所述步骤一中建立车轮力通用动态模型的具体步骤如下:
步骤1.1:定义k时刻的模型状态量为:
其中和分别为k时刻车轮纵向力、纵向力微分以及二次微分;和分别为k时刻车轮正压力、正压力微分以及二次微分;θk、和分别为k时刻车轮旋转角度、旋转角度微分以及二次微分;
步骤1.2:根据步骤1.1中定义的参数构建车轮力通用动态模型为:
XS,k=fS(T,α)XS,k-1+wS,k-1
其中fS(T,α)=diag[fS(T,αFxw)fS(T,αFzw)fS(T,αθ)],T为采样时间,Fx w为牵引力和Fz w为正压力,αFxw、αFzw和αθ分别为Fx w、Fz w和θ的变化率,wS,k-1为系统噪声;其中fS(T,αFxw)、fS(T,αFzw)和fS(T,αθ)的具体表达式为:
其中α=1/τ表示变量的变化率,它是机动时间τ的倒数,θ为偏转角度;
所述步骤二中建立车轮力高动态模型和车轮力地动态模型的具体步骤如下:选择αFxw=1/10,αFzw=1/10,αθ=1/10,建立车轮力高动态模型;选择αFxw=1/40,αFzw=1/40,αθ=1/30,建立车轮力低动态模型;
所述步骤三中建立全动态范围车轮力信号模型,并进行实时滤波的具体步骤如下:根据步骤二中所建立的车轮力高动态模型和车轮力低动态模型,分别对牵引力Fx w和正压力Fz w进行卡尔曼滤波,采用交互式多模型方法,实时计算上述两个模型权重,利用经过归一化权重求取两个模型滤波结果的加权和,并将该值作为轮力滤波的最终输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710030224.0A CN106874559B (zh) | 2017-01-17 | 2017-01-17 | 一种用于车轮力传感器的交互式多模型滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710030224.0A CN106874559B (zh) | 2017-01-17 | 2017-01-17 | 一种用于车轮力传感器的交互式多模型滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106874559A CN106874559A (zh) | 2017-06-20 |
CN106874559B true CN106874559B (zh) | 2019-05-31 |
Family
ID=59157686
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710030224.0A Active CN106874559B (zh) | 2017-01-17 | 2017-01-17 | 一种用于车轮力传感器的交互式多模型滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106874559B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875245B (zh) * | 2018-06-29 | 2021-05-04 | 电子科技大学 | 基于glmb滤波器的多传感器自适应角度控制方法 |
CN109606378B (zh) * | 2018-11-19 | 2020-06-09 | 江苏大学 | 面向非高斯噪声环境的车辆行驶状态估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721397A (zh) * | 2012-06-07 | 2012-10-10 | 江苏科技大学 | 基于垂直动载荷的现代时间序列提取路面特征参数的方法 |
CN102834700A (zh) * | 2010-04-01 | 2012-12-19 | 基斯特勒控股公司 | 校准wim传感器的方法 |
CN104442830A (zh) * | 2014-11-05 | 2015-03-25 | 山东省计算中心(国家超级计算济南中心) | 四轮转向/驱动车辆的速度估计和摩擦辨识方法 |
CN105653728A (zh) * | 2016-01-26 | 2016-06-08 | 大连理工大学 | 基于双层模型体系的无线传感网的数据采集方法 |
CN105675312A (zh) * | 2016-01-25 | 2016-06-15 | 中国汽车技术研究中心 | 一种模拟整车状态下车轮力传递函数测试方法及装置 |
-
2017
- 2017-01-17 CN CN201710030224.0A patent/CN106874559B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102834700A (zh) * | 2010-04-01 | 2012-12-19 | 基斯特勒控股公司 | 校准wim传感器的方法 |
CN102721397A (zh) * | 2012-06-07 | 2012-10-10 | 江苏科技大学 | 基于垂直动载荷的现代时间序列提取路面特征参数的方法 |
CN104442830A (zh) * | 2014-11-05 | 2015-03-25 | 山东省计算中心(国家超级计算济南中心) | 四轮转向/驱动车辆的速度估计和摩擦辨识方法 |
CN105675312A (zh) * | 2016-01-25 | 2016-06-15 | 中国汽车技术研究中心 | 一种模拟整车状态下车轮力传递函数测试方法及装置 |
CN105653728A (zh) * | 2016-01-26 | 2016-06-08 | 大连理工大学 | 基于双层模型体系的无线传感网的数据采集方法 |
Non-Patent Citations (2)
Title |
---|
"Real-time Signal Denoising Algorithm in Wheel Force Transducer";Dong Wang等;《Applied Mechanics and Materials》;20140811;第244-247页 |
"静态小波变换在汽车道路模拟试验中的应用";李忠国等;《数据采集与处理》;20070630;第22卷(第2期);第238-241页 |
Also Published As
Publication number | Publication date |
---|---|
CN106874559A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110320504B (zh) | 一种基于激光雷达点云统计几何模型的非结构化道路检测方法 | |
CN103318180B (zh) | 一种车辆路面不平度自动识别系统及方法 | |
CN102009654B (zh) | 一种全轮电驱动车辆的纵向车速估计方法 | |
CN108715166B (zh) | 基于深度学习的车辆稳定性指标估计方法 | |
CN106874559B (zh) | 一种用于车轮力传感器的交互式多模型滤波方法 | |
CN104050681A (zh) | 一种基于视频图像的道路消失点检测方法 | |
CN108007417A (zh) | 一种农机自动驾驶控制系统角度传感器自动标定方法 | |
CN102009653A (zh) | 融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法 | |
CN105716604A (zh) | 基于地磁序列的移动机器人室内定位方法及系统 | |
Reinstein et al. | Terrain adaptive odometry for mobile skid-steer robots | |
CN109466558A (zh) | 一种基于ekf和bp神经网络的路面附着系数估计方法 | |
CN105716617B (zh) | 基于车辆数据绘制驾驶轨迹的系统及方法 | |
DE102012200445A1 (de) | Autonomes Arbeitsgerät | |
CN104809292B (zh) | 一种高速列车非线性动力学模型参数的在线辨识方法 | |
CN107238500B (zh) | 汽车操纵稳定性试验快速评价系统建立方法 | |
CN102700551A (zh) | 一种车辆行驶过程中路面坡度的实时估计方法 | |
CN104062642B (zh) | 一种对激光雷达波形数据进行高斯回波分解的方法 | |
CN104021310A (zh) | 基于UKF与修正Dugoff轮胎模型路面峰值附着系数估算方法 | |
CN106127682A (zh) | 一种顾及海岸线地理特征约束的可控分形插值方法及系统 | |
CN105509771A (zh) | 一种发动机滑油金属颗粒在线监测的信号降噪方法 | |
CN109685935A (zh) | 一种整车系统状态量的估计方法和装置 | |
CN117633696B (zh) | 一种基于bim模型的水利工程施工监测方法及系统 | |
CN103792846A (zh) | 基于Skinner操作条件反射原理的机器人避障导航方法 | |
CN203283230U (zh) | 一种车辆路面不平度自动识别系统 | |
CN108120452A (zh) | Mems陀螺仪动态数据的滤波方法 |
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 |