CN113537620A - 基于Markov模型优化和工况识别的车速预测方法 - Google Patents

基于Markov模型优化和工况识别的车速预测方法 Download PDF

Info

Publication number
CN113537620A
CN113537620A CN202110869128.1A CN202110869128A CN113537620A CN 113537620 A CN113537620 A CN 113537620A CN 202110869128 A CN202110869128 A CN 202110869128A CN 113537620 A CN113537620 A CN 113537620A
Authority
CN
China
Prior art keywords
vehicle speed
working condition
markov
acceleration
speed prediction
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
Application number
CN202110869128.1A
Other languages
English (en)
Other versions
CN113537620B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202110869128.1A priority Critical patent/CN113537620B/zh
Publication of CN113537620A publication Critical patent/CN113537620A/zh
Application granted granted Critical
Publication of CN113537620B publication Critical patent/CN113537620B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/10Machine learning using kernel methods, e.g. support vector machines [SVM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Economics (AREA)
  • Medical Informatics (AREA)
  • Strategic Management (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Human Resources & Organizations (AREA)
  • General Engineering & Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Development Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于Markov模型优化和工况识别的车速预测方法,首先建立基于标准法规工况的基础Markov车速预测模型;然后分别在四种工况下借助模拟退火算法对基础Markov车速预测模型的加速度序列和速度序列进行优化,得到四种工况对应的优化Markov车速预测模型;接着,建立基于支持向量机的工况识别模型;最后,基于工况识别模型对当前车辆的工况进行识别,并根据识别出的工况采用对应的优化Markov车速预测模型进行车速预测。本发明能够根据工况在线识别的结果调用对应的Markov车速预测模型,提高车速预测的精度。

Description

基于Markov模型优化和工况识别的车速预测方法
技术领域
本发明涉及混合动力汽车能量管理领域,尤其涉及一种基于Markov模型优化和工况识别的车速预测方法。
背景技术
能量管理策略通过协调动力需求在多种能源之间的分配方式能在不牺牲动力需求的前提下提升整车的燃油经济性和排放性能,在混合动力汽车的控制中具有举足轻重的地位。现有能量管理策略一般分为基于规则和基于优化的策略。基于规则的策略控制时效性强,但对整车参数和工况的适用性差,且规则的制定需要丰富的工程经验。基于优化的策略分为全局优化策略和瞬时优化策略。全局优化策略能提供全局最优解,但需要预知完整的工况信息且计算量大,难以实现在线应用。瞬时优化策略针对当前时刻的工况信息寻求最优的动力需求分配方式,能实现在线的应用,但无法保证全局最优的控制效果。模型预测控制策略综合了全局优化策略和瞬时优化策略的优势,通过进行局部的车速预测及优化控制,既考虑了未来的行程信息,又极大地缩短了优化的时间,有助于实现策略的在线应用。
对未来工况的预测是模型预测控制策略的一个主要组成部分,通常采用Markov模型或者神经网络模型进行车速的预测,这两种基础车速预测都已经应用成熟,研究的热点集中在如何结合工况和驾驶员等不确定因素建立更为精准的车速预测模型。
发明内容
本发明所要解决的技术问题是针对背景技术中所涉及到的缺陷,提供一种基于Markov模型优化和工况识别的车速预测方法。
本发明为解决上述技术问题采用以下技术方案:
基于Markov模型优化和工况识别的车速预测方法,包括以下步骤:
步骤1),建立基于标准法规工况的基础Markov车速预测模型;
所述标准法规工况包含高速工况、城郊工况、城市畅通工况和城市拥堵工况四种工况;所述基础Markov车速预测模型采用单阶多步车速预测模型,用于描述车辆在当前的加速度与车速状态下每个预测时刻车辆的加速度状态转移概率分布,包含加速度序列as0、速度序列Vs0、以及大小为La*La*Lv*Lp的概率矩阵,其中,as0为长度为La的amin到amax的等间距序列,Vs0为长度为Lv的0到Vmax的等间距序列,La为离散加速度状态数,Lv为离散速度状态数,Lp为预测步长,amin为预设的最小加速度阈值,amax为预设的最大加速度阈值,Vmax为预设的最大速度阈值;
步骤2),分别在四种工况下借助模拟退火算法对基础Markov车速预测模型的加速度序列和速度序列进行优化,得到四种工况对应的优化Markov车速预测模型;
步骤3),建立基于支持向量机的工况识别模型,用于识别高速工况、城郊工况、城市畅通工况和城市拥堵工况,其输入为给定时间长度tS内的平均车速、最大加速度、最小加速度和停车时间比,输出为工况类别;
步骤4),基于工况识别模型对当前车辆的工况进行识别,并根据识别出的工况采用对应的优化Markov车速预测模型进行车速预测。
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,步骤1)中基础Markov车速预测模型的训练过程如下:
步骤1.1),设置大小为La*La*Lv*Lp的频数矩阵T并将其零初始化,设置k=1;
步骤1.2),令过程预测步长q=1;
步骤1.3),获取标准法规工况的速度V(k)、加速度a(k)和a(k+q),采用最邻近法则得到V(k)在速度序列中对应的离散值Vn、a(k)在加速度序列中对应的离散值aj、a(k+q)在加速度序列中对应的离散值ai,n为离散值Vn在速度序列中的位置编号、j为离散值aj在加速度序列中的位置编号、i为离散值ai在加速度序列中的位置编号,依据位置编号n、j、i、q对频数矩阵T的相应位置累积加1,即T(i,j,n,q)累积加1,i,j∈{1,2,3,...,La},n∈{1,2,3,...,Lv};
步骤1.4),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.3);
步骤1.5),更新k=k+1,并比较k和L-Lp之间的大小,L为标准法规工况的长度,如果k≤L-Lp,跳转执行步骤1.2);
步骤1.6),依据频数矩阵完成概率矩阵,表示如下:
Figure BDA0003188435790000021
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,步骤1)中所述基础Markov车速预测模型的预测过程如下:
步骤1.a),获取当前状态的车速V(t)与加速度a(t),t是车辆的真实运行时间,采用最邻近法则得到V(t)在速度序列中对应的离散值Vn1、a(t)在加速度序列中对应的离散值aj1;n1和为离散值Vn1在速度序列中的位置编号、j1为离散值aj1在加速度序列中的位置编号,令q=1;
步骤1.b),将概率矩阵P(:,j1,n1,q)中最大值所在位置记为L1,则q步车速预测值为:
V(t+q)=V(t+q-1)+as0(L1)*3.6
步骤1.c),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.b),否则所述基础Markov车速预测模型的预测过程结束。
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,所述步骤2)包含以下具体步骤:
对于四种工况中的每一种工况:
步骤2.1),以加速度序列as和速度序列Vs中的具体数值为优化变量、以预测误差为优化目标,将优化变量记为x=[as,Vs],优化目标E表示如下:
Figure BDA0003188435790000031
其中,
Figure BDA0003188435790000032
是所有预测时刻的均方根误差的均值,
Figure BDA0003188435790000033
是所有预测时刻的平均绝对误差的均值,
Figure BDA0003188435790000034
为预测的车速,Vt为真实的车速;
步骤2.2),采用模拟退火算法对优化变量进行优化,模拟退火算法能以一定的概率接受比当前解更差的解,具有跳出局部最优陷阱的能力:
步骤2.2.1),设定初始温度TSAA,初始优化变量xold为步骤1)中的加速度与速度序列,令Q=0;
步骤2.2.2),产生新优化变量xnew,将当前优化变量xold与新优化变量xnew带入Markov车速预测模型,求解得到对应的优化目标值Eold和Enew,并计算新优化变量的可接受概率PSAA
Figure BDA0003188435790000035
步骤2.2.3),当Enew≥Eold时,产生一个随机数ε∈(0,1),当ε>PSAA时,优化目标值增加的新状态不被接受,执行步骤2.2.5),否则,接受新状态,执行步骤2.2.4);
步骤2.2.4),更新xold=xnew、TSAA=α*TSAA、Q=0,然后跳转执行步骤2.2.2),其中,α为温度参数衰减系数,α∈[0,1];
步骤2.2.5),更新Q=Q+1,判断Q是否大于预设的阈值N,如果Q小于N,跳转执行步骤2.2.2),如果Q大于等于N,优化结束,将此时的Markov车速预测模型记录为工况对应的优化Markov车速预测模型。
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,步骤3)中的工况识别模型借助MATLAB的Classification Learner工具箱进行训练,其中训练数据的生成步骤如下:
对于四种工况中的每一种工况:
步骤3.1),预先设定样本量为K,i=1,l为工况的长度;
步骤3.2),随机产生起始点n0,满足n0+ts<l;
步骤3.3),i≤K,则记录
Figure BDA0003188435790000041
为该工况的样本Si,否则,执行步骤3.5);
步骤3.4),更新n0=n0+ts,i=i+1,比较n0+ts和l之间的大小,如果n0+ts<l,跳转执行3.3),否则跳转执行步骤3.2);
步骤3.5),针对记录的样本进行特征参数提取,所述特征参数包含平均车速、最大加速度、最小加速度和停车时间比,再根据预先设定的工况类别代号对样本进行工况标记;
步骤3.6),将特征参数及工况类别组合作为该工况的训练数据。
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,所述步骤4)具体包含以下步骤:
步骤4.1),当车辆实际运行时间t<ts时,采用预设的默认工况,并调用预设默认工况对应的优化Markov车速预测模块进行车速预测;
步骤4.2),当t≥ts时,按照预设的间隔时间阈值tup对[t-ts,t-1]的工况进行特征参数提取并输入工况识别模型,并依据工况识别模型识别的工况调用对应的优化Markov车速预测模块进行车速预测。
作为本发明基于Markov模型优化和工况识别的车速预测方法进一步的优化方案,步骤4.1)中预设的默认工况为城市畅通工况。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明结合工况识别,对基础Markov车速预测模型的加速度和速度序列进行优化,能有效提高车速预测的精度,应用在基于模型预测控制的混合动力汽车能量管理策略中能够进一步提升车辆的燃油经济性。
附图说明
图1为本发明基于Markov模型优化和工况识别的车速预测方法的流程示意图;
图2为本发明中借助模拟退火算法对基础Markov车速预测模型的加速度序列和速度序列进行优化的流程示意图;
图3为本发明中建立基于支持向量机的工况识别模型过程中的模型训练数据生成流程示意图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明可以以许多不同的形式实现,而不应当认为限于这里所述的实施例。相反,提供这些实施例以便使本公开透彻且完整,并且将向本领域技术人员充分表达本发明的范围。在附图中,为了清楚起见放大了组件。
如图1所示,基于Markov模型优化和工况识别的车速预测方法,包括以下步骤:
步骤1),建立基于标准法规工况的基础Markov车速预测模型;
所述标准法规工况包含高速工况、城郊工况、城市畅通工况和城市拥堵工况四种工况;所述基础Markov车速预测模型采用单阶多步车速预测模型,用于描述车辆在当前的加速度与车速状态下每个预测时刻车辆的加速度状态转移概率分布,包含加速度序列as0、速度序列Vs0、以及大小为La*La*Lv*Lp的概率矩阵,其中,as0为长度为La的amin到amax的等间距序列,Vs0为长度为Lv的0到Vmax的等间距序列,La为离散加速度状态数,Lv为离散速度状态数,Lp为预测步长,amin为预设的最小加速度阈值,amax为预设的最大加速度阈值,Vmax为预设的最大速度阈值;
步骤2),分别在四种工况下借助模拟退火算法对基础Markov车速预测模型的加速度序列和速度序列进行优化,得到四种工况对应的优化Markov车速预测模型;
步骤3),建立基于支持向量机的工况识别模型,用于识别高速工况、城郊工况、城市畅通工况和城市拥堵工况,其输入为给定时间长度tS内的平均车速、最大加速度、最小加速度和停车时间比,输出为工况类别;
步骤4),基于工况识别模型对当前车辆的工况进行识别,并根据识别出的工况采用对应的优化Markov车速预测模型进行车速预测。
步骤1)中基础Markov车速预测模型的训练过程如下:
步骤1.1),设置大小为La*La*Lv*Lp的频数矩阵T并将其零初始化,设置k=1;
步骤1.2),令过程预测步长q=1;
步骤1.3),获取标准法规工况的速度V(k)、加速度a(k)和a(k+q),采用最邻近法则得到V(k)在速度序列中对应的离散值Vn、a(k)在加速度序列中对应的离散值aj、a(k+q)在加速度序列中对应的离散值ai,n为离散值Vn在速度序列中的位置编号、j为离散值aj在加速度序列中的位置编号、i为离散值ai在加速度序列中的位置编号,依据位置编号n、j、i、q对频数矩阵T的相应位置累积加1,即T(i,j,n,q)累积加1,i,j∈{1,2,3,...,La},n∈{1,2,3,...,Lv};
步骤1.4),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.3);
步骤1.5),更新k=k+1,并比较k和L-Lp之间的大小,L为标准法规工况的长度,如果k≤L-Lp,跳转执行步骤1.2);
步骤1.6),依据频数矩阵完成概率矩阵,表示如下:
Figure BDA0003188435790000061
步骤1)中所述基础Markov车速预测模型的预测过程如下:
步骤1.a),获取当前状态的车速V(t)与加速度a(t),t是车辆的真实运行时间,采用最邻近法则得到V(t)在速度序列中对应的离散值Vn1、a(t)在加速度序列中对应的离散值aj1;n1和为离散值Vn1在速度序列中的位置编号、j1为离散值aj1在加速度序列中的位置编号,令q=1;
步骤1.b),将概率矩阵P(:,j1,n1,q)中最大值所在位置记为L1,则q步车速预测值为:
V(t+q)=V(t+q-1)+as0(L1)*3.6
步骤1.c),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.b),否则所述基础Markov车速预测模型的预测过程结束。
如图2所示,所述步骤2)包含以下具体步骤:
对于四种工况中的每一种工况:
步骤2.1),以加速度序列as和速度序列Vs中的具体数值为优化变量、以预测误差为优化目标,将优化变量记为x=[as,Vs],优化目标E表示如下:
Figure BDA0003188435790000062
其中,
Figure BDA0003188435790000063
是所有预测时刻的均方根误差的均值,
Figure BDA0003188435790000064
是所有预测时刻的平均绝对误差的均值,
Figure BDA0003188435790000065
为预测的车速,Vt为真实的车速;
步骤2.2),采用模拟退火算法对优化变量进行优化,模拟退火算法能以一定的概率接受比当前解更差的解,具有跳出局部最优陷阱的能力:
步骤2.2.1),设定初始温度TSAA,初始优化变量xold为步骤1)中的加速度与速度序列,令Q=0;
步骤2.2.2),产生新优化变量xnew,将当前优化变量xold与新优化变量xnew带入Markov车速预测模型,求解得到对应的优化目标值Eold和Enew,并计算新优化变量的可接受概率PSAA
Figure BDA0003188435790000071
步骤2.2.3),当Enew≥Eold时,产生一个随机数ε∈(0,1),当ε>PSAA时,优化目标值增加的新状态不被接受,执行步骤2.2.5),否则,接受新状态,执行步骤2.2.4);
步骤2.2.4),更新xold=xnew、TSAA=α*TSAA、Q=o,然后跳转执行步骤2.2.2),其中,α为温度参数衰减系数,α∈[0,1];
步骤2.2.5),更新Q=Q+1,判断Q是否大于预设的阈值N,如果Q小于N,跳转执行步骤2.2.2),如果Q大于等于N,优化结束,将此时的Markov车速预测模型记录为工况对应的优化Markov车速预测模型。
步骤3)中的工况识别模型借助MATLAB的Classification Learner工具箱进行训练,其中训练数据的生成步骤如下:
如图3所示,对于四种工况中的每一种工况:
步骤3.1),预先设定样本量为K,i=1,l为工况的长度;
步骤3.2),随机产生起始点n0,满足n0+ts<l;
步骤3.3),i≤K,则记录
Figure BDA0003188435790000072
为该工况的样本Si,否则,执行步骤3.5);
步骤3.4),更新n0=n0+ts,i=i+1,比较n0+ts和l之间的大小,如果n0+ts<l,跳转执行3.3),否则跳转执行步骤3.2);
步骤3.5),针对记录的样本进行特征参数提取,所述特征参数包含平均车速、最大加速度、最小加速度和停车时间比,再根据预先设定的工况类别代号对样本进行工况标记;
步骤3.6),将特征参数及工况类别组合作为该工况的训练数据。
所述步骤4)具体包含以下步骤:
步骤4.1),当车辆实际运行时间t<ts时,采用预设的默认工况,并调用预设默认工况对应的优化Markov车速预测模块进行车速预测;
步骤4.2),当t之ts时,按照预设的间隔时间阈值tup对[t-ts,t-1]的工况进行特征参数提取并输入工况识别模型,并依据工况识别模型识别的工况调用对应的优化Markov车速预测模块进行车速预测。
步骤4.1)中预设的默认工况优先采用城市畅通工况。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.基于Markov模型优化和工况识别的车速预测方法,其特征在于,包括以下步骤:
步骤1),建立基于标准法规工况的基础Markov车速预测模型;
所述标准法规工况包含高速工况、城郊工况、城市畅通工况和城市拥堵工况四种工况;所述基础Markov车速预测模型采用单阶多步车速预测模型,用于描述车辆在当前的加速度与车速状态下每个预测时刻车辆的加速度状态转移概率分布,包含加速度序列as0、速度序列Vs0、以及大小为La*La*Lv*Lp的概率矩阵,其中,as0为长度为La的amin到amax的等间距序列,Vs0为长度为Lv的0到Vmax的等间距序列,La为离散加速度状态数,Lv为离散速度状态数,Lp为预测步长,amin为预设的最小加速度阈值,amax为预设的最大加速度阈值,Vmax为预设的最大速度阈值;
步骤2),分别在四种工况下借助模拟退火算法对基础Markov车速预测模型的加速度序列和速度序列进行优化,得到四种工况对应的优化Markov车速预测模型;
步骤3),建立基于支持向量机的工况识别模型,用于识别高速工况、城郊工况、城市畅通工况和城市拥堵工况,其输入为给定时间长度ts内的平均车速、最大加速度、最小加速度和停车时间比,输出为工况类别;
步骤4),基于工况识别模型对当前车辆的工况进行识别,并根据识别出的工况采用对应的优化Markov车速预测模型进行车速预测。
2.根据权利要求1所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,步骤1)中基础Markov车速预测模型的训练过程如下:
步骤1.1),设置大小为La*La*Lv*Lp的频数矩阵T并将其零初始化,设置k=1;
步骤1.2),令过程预测步长q=1;
步骤1.3),获取标准法规工况的速度V(k)、加速度a(k)和a(k+q),采用最邻近法则得到V(k)在速度序列中对应的离散值Vn、a(k)在加速度序列中对应的离散值aj、a(k+q)在加速度序列中对应的离散值ai,n为离散值Vn在速度序列中的位置编号、j为离散值aj在加速度序列中的位置编号、i为离散值ai在加速度序列中的位置编号,依据位置编号n、j、i、q对频数矩阵T的相应位置累积加1,即T(i,j,n,q)累积加1,i,j∈{1,2,3,...,La},n∈{1,2,3,...,Lv};
步骤1.4),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.3);
步骤1.5),更新k=k+1,并比较k和L-Lp之间的大小,L为标准法规工况的长度,如果k≤L-Lp,跳转执行步骤1.2);
步骤1.6),依据频数矩阵完成概率矩阵,表示如下:
Figure FDA0003188435780000021
3.根据权利要求2所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,步骤1)中所述基础Markov车速预测模型的预测过程如下:
步骤1.a),获取当前状态的车速V(t)与加速度a(t),t是车辆的真实运行时间,采用最邻近法则得到V(t)在速度序列中对应的离散值Vn1、a(t)在加速度序列中对应的离散值aj1;n1和为离散值Vn1在速度序列中的位置编号、j1为离散值aj1在加速度序列中的位置编号,令q=1;
步骤1.b),将概率矩阵P(:,j1,n1,q)中最大值所在位置记为L1,则q步车速预测值为:
V(t+q)=V(t+q-1)+as0(L1)*3.6
步骤1.c),更新q=q+1,比较q和Lp之间的大小,如果q≤Lp,跳转执行步骤1.b),否则所述基础Markov车速预测模型的预测过程结束。
4.根据权利要求3所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,所述步骤2)包含以下具体步骤:
对于四种工况中的每一种工况:
步骤2.1),以加速度序列as和速度序列Vs中的具体数值为优化变量、以预测误差为优化目标,将优化变量记为x=[as,Vs],优化目标E表示如下:
Figure FDA0003188435780000022
其中,
Figure FDA0003188435780000023
是所有预测时刻的均方根误差的均值,
Figure FDA0003188435780000024
是所有预测时刻的平均绝对误差的均值,
Figure FDA0003188435780000025
为预测的车速,Vt为真实的车速;
步骤2.2),采用模拟退火算法对优化变量进行优化,模拟退火算法能以一定的概率接受比当前解更差的解,具有跳出局部最优陷阱的能力:
步骤2.2.1),设定初始温度TSAA,初始优化变量xold为步骤1)中的加速度与速度序列,令Q=0;
步骤2.2.2),产生新优化变量xnew,将当前优化变量xold与新优化变量xnew带入Markov车速预测模型,求解得到对应的优化目标值Eold和Enew,并计算新优化变量的可接受概率PSAA
Figure FDA0003188435780000026
步骤2.2.3),当Enew≥Eold时,产生一个随机数ε∈(0,1),当ε>PSAA时,优化目标值增加的新状态不被接受,执行步骤2.2.5),否则,接受新状态,执行步骤2.2.4);
步骤2.2.4),更新xold=xnew、TSAA=α*TSAA、Q=0,然后跳转执行步骤2.2.2),其中,α为温度参数衰减系数,α∈[0,1];
步骤2.2.5),更新Q=Q+1,判断Q是否大于预设的阈值N,如果Q小于N,跳转执行步骤2.2.2),如果Q大于等于N,优化结束,将此时的Markov车速预测模型记录为工况对应的优化Markov车速预测模型。
5.根据权利要求4所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,步骤3)中的工况识别模型借助MATLAB的Classification Learner工具箱进行训练,其中训练数据的生成步骤如下:
对于四种工况中的每一种工况:
步骤3.1),预先设定样本量为K,i=1,l为工况的长度;
步骤3.2),随机产生起始点n0,满足n0+ts<l;
步骤3.3),i≤K,则记录
Figure FDA0003188435780000031
为该工况的样本Si,否则,执行步骤3.5);
步骤3.4),更新n0=n0+ts,i=i+1,比较n0+ts和l之间的大小,如果n0+ts<l,跳转执行3.3),否则跳转执行步骤3.2);
步骤3.5),针对记录的样本进行特征参数提取,所述特征参数包含平均车速、最大加速度、最小加速度和停车时间比,再根据预先设定的工况类别代号对样本进行工况标记;
步骤3.6),将特征参数及工况类别组合作为该工况的训练数据。
6.根据权利要求5所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,所述步骤4)具体包含以下步骤:
步骤4.1),当车辆实际运行时间t<ts时,采用预设的默认工况,并调用预设默认工况对应的优化Markov车速预测模块进行车速预测;
步骤4.2),当t≥ts时,按照预设的间隔时间阈值tup对[t-ts,t-1]的工况进行特征参数提取并输入工况识别模型,并依据工况识别模型识别的工况调用对应的优化Markov车速预测模块进行车速预测。
7.根据权利要求6所述的基于Markov模型优化和工况识别的车速预测方法,其特征在于,步骤4.1)中预设的默认工况为城市畅通工况。
CN202110869128.1A 2021-07-30 2021-07-30 基于Markov模型优化和工况识别的车速预测方法 Active CN113537620B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110869128.1A CN113537620B (zh) 2021-07-30 2021-07-30 基于Markov模型优化和工况识别的车速预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110869128.1A CN113537620B (zh) 2021-07-30 2021-07-30 基于Markov模型优化和工况识别的车速预测方法

Publications (2)

Publication Number Publication Date
CN113537620A true CN113537620A (zh) 2021-10-22
CN113537620B CN113537620B (zh) 2024-04-16

Family

ID=78089848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110869128.1A Active CN113537620B (zh) 2021-07-30 2021-07-30 基于Markov模型优化和工况识别的车速预测方法

Country Status (1)

Country Link
CN (1) CN113537620B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116756590A (zh) * 2023-08-17 2023-09-15 天津德科智控股份有限公司 Eps系统车速信号干扰识别及处理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107909179A (zh) * 2017-09-29 2018-04-13 北京理工大学 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法
CN109927709A (zh) * 2017-12-15 2019-06-25 郑州宇通客车股份有限公司 一种车辆行驶路线工况确定方法、能量管理方法及系统
US20190202435A1 (en) * 2017-12-29 2019-07-04 Hyundai Motor Company Method for identifying driving pattern of hybrid vehicle
CN111459995A (zh) * 2020-03-11 2020-07-28 南京航空航天大学 一种基于驾驶数据的多模态车速预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107909179A (zh) * 2017-09-29 2018-04-13 北京理工大学 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法
CN109927709A (zh) * 2017-12-15 2019-06-25 郑州宇通客车股份有限公司 一种车辆行驶路线工况确定方法、能量管理方法及系统
US20190202435A1 (en) * 2017-12-29 2019-07-04 Hyundai Motor Company Method for identifying driving pattern of hybrid vehicle
CN111459995A (zh) * 2020-03-11 2020-07-28 南京航空航天大学 一种基于驾驶数据的多模态车速预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
丁峰;王伟达;项昌乐;何;齐蕴龙;: "基于行驶工况分类的混合动力车辆速度预测方法与能量管理策略", 汽车工程, no. 11, 25 November 2017 (2017-11-25), pages 1123 - 1131 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116756590A (zh) * 2023-08-17 2023-09-15 天津德科智控股份有限公司 Eps系统车速信号干扰识别及处理方法
CN116756590B (zh) * 2023-08-17 2023-11-14 天津德科智控股份有限公司 Eps系统车速信号干扰识别及处理方法

Also Published As

Publication number Publication date
CN113537620B (zh) 2024-04-16

Similar Documents

Publication Publication Date Title
WO2021103625A1 (zh) 一种基于前车与自车互动的短期车速工况实时预测方法
CN109017809B (zh) 一种基于越野工况预测的能量分配方法
Lin et al. An ensemble learning velocity prediction-based energy management strategy for a plug-in hybrid electric vehicle considering driving pattern adaptive reference SOC
CN108520155B (zh) 基于神经网络的车辆行为模拟方法
CN107909179B (zh) 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法
CN111047085B (zh) 一种基于元学习的混合动力车辆工况预测方法
CN111284489B (zh) 智能网联汽车随机预测巡航控制系统
CN107506865B (zh) 一种基于lssvm优化的负荷预测方法及系统
CN112101684B (zh) 一种插电式混合动力车辆实时能量管理方法和系统
CN110837915B (zh) 一种基于混合集成深度学习的电力系统低压负荷点预测及概率预测方法
CN110674965A (zh) 基于动态特征选取的多时间步长风功率预测方法
CN112242059A (zh) 基于动机与风险评估的无人驾驶车辆智能决策方法
CN113479187B (zh) 一种插电式混合动力汽车分层异步长能量管理方法
CN111680413B (zh) 基于双层算法的有轨电车定时节能运行优化方法及其系统
Montazeri-Gh et al. Driving condition recognition for genetic-fuzzy HEV control
CN113135113B (zh) 一种全局soc规划方法及装置
CN115107733A (zh) 一种混合动力汽车的能量管理方法及系统
CN113537620A (zh) 基于Markov模型优化和工况识别的车速预测方法
CN116244647A (zh) 一种无人机集群的运行状态估计方法
Billert et al. A method of developing quantile convolutional neural networks for electric vehicle battery temperature prediction trained on cross-domain data
Gan et al. Intelligent learning algorithm and intelligent transportation-based energy management strategies for hybrid electric vehicles: a review
CN111516702B (zh) 一种混合动力车辆在线实时分层能量管理方法和系统
CN117465301A (zh) 基于数据驱动的燃料电池汽车实时能量管理方法
CN116968721A (zh) 一种混合动力汽车预测式能量管理方法、系统和存储介质
CN116901981A (zh) 一种在线自学习型的马尔可夫车辆速度预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant