CN113741199B - 一种基于智能网联信息的整车经济性速度规划方法 - Google Patents

一种基于智能网联信息的整车经济性速度规划方法 Download PDF

Info

Publication number
CN113741199B
CN113741199B CN202111087524.5A CN202111087524A CN113741199B CN 113741199 B CN113741199 B CN 113741199B CN 202111087524 A CN202111087524 A CN 202111087524A CN 113741199 B CN113741199 B CN 113741199B
Authority
CN
China
Prior art keywords
vehicle
speed
terminal
time
network
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
Application number
CN202111087524.5A
Other languages
English (en)
Other versions
CN113741199A (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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN202111087524.5A priority Critical patent/CN113741199B/zh
Publication of CN113741199A publication Critical patent/CN113741199A/zh
Application granted granted Critical
Publication of CN113741199B publication Critical patent/CN113741199B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明属于智能网技术领域,具体的说是一种基于智能网联信息的整车经济性速度规划方法。包括以下步骤:步骤一、通过智能网联环境获得路段上的动态交通信息,通过传感器获得本车的状态,确定车辆的安全边界即最大的行驶速度和最大的加速度;步骤二、建立车辆运动学模型,并且分析车辆行驶过程中的能量消耗;步骤三、根据动态交通信息设定***的性能指标和约束条件;步骤四、基于滚动时域强化学习,在线求解优化速度即对速度进行规划。本发明充分利用强化学习的求解力,并且融合了预测控制解决约束及扰动的优势,在线求解非线性时变最优化问题,具有明显节能优势和工程应用能力,解决了长预测时域控制求解时间问题。

Description

一种基于智能网联信息的整车经济性速度规划方法
技术领域
本发明属于智能网技术领域,具体的说是一种基于智能网联信息的整车经济性速度规划方法。
背景技术
融合外界多源信息和汽车传统节能技术的智能网联纯电动汽车已经逐步成为新能源汽车领域的研究重点,通过优化车辆的行驶速度,减少车辆因交通环境带来的频繁的加减速,进而有效减少等待时间和燃料消耗。经济性行驶对于节能控制具有良好的应用前景,但是该技术也存在诸多挑战。因经济性行驶的速度规划是基于宏观交通信息,具有时空跨度大,预测时域长的特点,动态规划、预测控制等方法在求解过程中会遇到高维问题,无法在线求解及工程应用。
发明内容
本发明提供了一种基于智能网联信息的整车经济性速度规划方法,该规划方法充分利用强化学习的求解力,并且融合了预测控制解决约束及扰动的优势,在线求解非线性时变最优化问题,具有明显节能优势和工程应用能力,解决了长预测时域控制求解时间问题。
本发明技术方案结合附图说明如下:
一种基于智能网联信息的整车经济性速度规划方法,包括以下步骤:
步骤一、通过交通信息层的智能网联环境获得路段上的动态交通信息,并通过传感器获得本车的状态,确定车辆的安全边界即最大的行驶速度和最大的加速度;
步骤二、建立车辆运动学模型,并且分析车辆行驶过程中的能量消耗;
步骤三、根据动态交通信息设定***的性能指标和约束条件;
步骤四、在速度轨迹规划层基于滚动时域强化学习,在线求解优化速度即对速度进行规划,并通过轨迹规划层控制四轮力矩,驱动电动汽车按照需求力矩行驶。
步骤一中所述动态交通信息包括外交通灯位置和时序、限速及前车信息。
所述步骤二的具体方法如下:
在上层经济性行驶速度轨迹规划中,只考虑车辆的纵向运动学,运动学模型如下:
式中,d为行驶距离;v为车速;m为整车质量;δ为旋转质量系数;Ft为总驱动力;FR为行驶阻力,包括空气阻力,滚动阻力和坡度阻力,表示如下:
式中,fw为空气阻力;ff为滚动阻力;fi为坡度阻力;CD为阻力系数;A为迎风面积;v为车速;fr为滚动阻力系数;m为整车质量;θ为坡度;
车辆在加速行驶过程中,驱动力由轮毂电机提供,车辆在制动过程中,总驱动力Ft等于0,制动力由液压制动提供,表达式如下:
式中,Tmi是每个电机的驱动力矩i=1,2,3,4,Tm1、Tm2、Tm3、Tm4分别为左前轮、右前轮、左后轮、右后轮,rw为轮胎半径;
能源消耗总量为:
P=Ftv=(mδa+FR)v (5)
式中,Ft为总驱动力;v为车速;m为整车质量;δ为旋转质量系数;a为加速度;FR为行驶阻力。
所述步骤三的具体方法如下:
31)将车辆的行驶距离d和车速v作为状态变量,将能耗公式(1)—(2)进行离散表示如下:
v(k+1)=v(k)+a(k)Δt (7)
式中,k为当前时刻;Δt为***采样时间;v(k)为当前时刻的车速;a(k)为***的加速度,经济性行驶被描述为一个加速度的优化问题;状态变量:x=[x1,x2]=[d,v]T,控制输入:u=[a],式中,d为行驶距离;v为车速;a为加速度;
32)运动学模型表示为:
x(k+1)=x(k)+f(x(k),u(k)) (8)
式中,x(k)为***状态;u(k)为控制输入;f(x(k),u(k))为由控制输入引起的状态变化;
33)动力学***的性能指标函数为:
式中,k+N是终端约束项;L为电动汽车能量消耗,
L(x,u,k)=β1P(x,u,k)+β2||u(k)-u(k-1)||2 (10)
式中,P(x,u,k)为燃料消耗,其由公式(5)计算得到;x为***状态;u为控制输入;k表示当前时刻;||u(k)-u(k-1)||2为***控制输入的幅值约束,防止过大的加速度,导致舒适性下降,u(k)为***的控制输入;u(k-1)为上一时刻***控制输入;β1,β2为权重因子;T为设定预测时间长度;在性能指标公式(9)中,预测时间步长度N=T/Δt,且随着运行时间递减;车辆行驶的终端时间根据前方交通灯时序确定,即终端时间kf是固定值;每个时刻下的预测时间长度T=kf-k;式中,kf为终端时间;k为表示当前时刻;
φ[x(k+N)]为***的终端状态约束,定义为:
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (11)
其中,kf为终端时间;sf为终端距离,由路段上交通灯时序及位置决定;vf为终端速度;a1、a2为权重系数;v(kf)为终端时刻车速;s(kf)为终端时刻位置;
34)根据动态交通信息中的交通灯的实际情况,转化为规划问题的约束条件;当车辆接近交通灯时,具有三种情况,即加速、巡航和减速;根据道路要求1设置车辆自适应巡航初始速度为vd,而车辆在时域末端的速度vf根据是否通过红绿灯确定;具体如下:
a、当车辆在低密度路况且无前车的情况下,考虑前方红绿灯信号作用,将红绿灯信号的时间及距离转化为速度约束,当车辆保持当前车速不能通过下一个红绿灯路口时,需要选择提前加速或者减速达到节能通过红绿灯或停车的目的;通过与红绿灯距离Δd和绿灯所剩时间tgreen进行判断,当车速提高到vmax时,若仍然不能在绿灯时刻通过路口,则选择提前减速,车辆准备在下一个绿灯周期通过路口,即终端速度约束vf=0;反之,车辆定速巡航或加速通过路口,终端速度约束vf=vd;在自适应巡航控制中,设计速度浮动参数ε,即车辆的容许范围∈[1-ε,1+ε];车辆的速度上边界为:vmax=min(vd(1+ε),vlim),vlim为道路限速;
b、当与前车的距离宽裕且能满足速度及加速度需求的情况,本车仍然根据剩余绿灯时间进行跟进或者减速度到停止线停车的控制;反之,车辆的速度容许度受前车影响,车辆必须遵守跟车规则;设定前车的行驶速度为vleader,车辆速度上边界表示为:vmax=min(vleader,vd(1+ε),vlim),通过vmax*tgreen≥Δd的关系,判断是否通过路口,即而得到终端速度条件;
35)通过31)、32)、33)和34)得到***需满足的约束条件如下;
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (12a)
vmin(k)≤v(k)≤vmax(k) (12b)
amin(k)≤a(k)≤amax(k) (12c)
v0=vd,s0=0 (12d)
式中,为***的终端状态约束;a1、a2为权重系数;kf为终端时间;v(kf)为终端时刻车速;s(kf)为终端时刻位置;sf为终端距离;v(k)为终端速度约束;vmin(k)为车辆速度的下边界;vmax(k)为车辆速度的上边界;a(k)为车辆加速度;amin(k)为车辆加速度的下边界;amax(k)为车辆加速度的上边界;v0为初始速度,根据自适应控制车速设定;vd为终端速度约束;s0车辆的初始位置,通常选择为0。
所述步骤四为:
根据步骤二和步骤三获得的车辆动力学模型、性能指标和约束条件构建优化问题,基于滚动时域强化学习,得到优化后的加速度值,并选择控制序列即加速度值的前Tc(Tc≤N)个控制动作作用于***中,Tc为选择的控制序列长度;车辆根据加速度需求值控制目标车辆驱动力输出,当到达时刻k+Tc时,重复求解控制序列,实现滚动优化控制。
所述步骤四的具体方法如下:
41)根据经济性巡航任务要求,设计单步回报函数表示如下:
r(x(k),u(k))=β1P(x(k),u(k))+β2||u(k)-u(k-1)||2 (13)
式中,β1,β2为权重因子;x(k)为***状态;u(k)为控制输入;P(x(k),u(k))为车辆燃料消耗;u(k-1)为***上一时刻控制量;
在有限时域内,代价函数由终端代价和累积回报两部分组成:
式中,EU[·]为在当前策略下的期望算子;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数,根据公式(12a)确定;γ为折扣因子,且0≤γ≤1;k表示当前时刻;N为预测时域长度;x(j)为当前时刻状态;u(j)为当前时刻控制输入;
42)设计动作网络,动作网络用于逼近最优控制输入u*(k),采用模糊小脑模型神经网络作为逼近器,模糊小脑模型神经网络包括五个模块,分别为输入变量、变量模糊化、概念映射、实际映射和输出变量,具体如下:
a、输入变量及变量模糊化,将d,v的模糊集均表示为[NB,NS,ZE,PS,PB],其中NB,NS,ZE,PS,PB分别表示为负大、负小、0、正小、正大;两者的论域[-2,-1,0,1,2];隶属度函数向量λj=[λ1j,…λnj]T,其中,j=1,2,n=1,2…5;模糊化变量选择高斯隶属度函数:
式中,ci为第i个隶属度函数中心点,σi为方差;
b、模糊隶属度至输入空间U,将空间U划分为10个存储单元,每个单元对应一个向量,从空间划分的存储单独找到对应地址;
c、输入空间U到Ac的概念映射,输入空间U的10个模糊隶属度函数进行规则划分后得到52个状态,每个状态作为一个指针映射到虚拟空间Ac的C个存储单元中,并找到对应该状态的地址;
d、概念映射Ac到实际映射Ap,采用杂散编码技术中的除留余数法,将概念映射的C个单元映射到Ap的C个单元;
e、实际映射Ap到输出变量模糊小脑模型神经网络的输出表示为:
式中,am(k)为高斯隶属度函数映射的乘积;wm(k)为相应的权值;动作网络的学习目标为极小化动作网络估计值与最优值u*(k)之间的误差值,误差函数为:
权值wa(k)修正迭代公式表示为:
式中,la为动作网络学习率,它随时间的增加而逐渐减小为一个固定值,θ∈(0,1]为网络惯性系数;λi(k)为隶属度函数向量;为前一时刻权值矩阵;/>为(k-2)时刻权值矩阵;
43)设计评价网络,评价网络采用模糊小脑模型神经网络作为逼近器,评价网络的输出表示为:
式中,wc,n(k)为当前时刻权值;an(k)为状态在高斯隶属度函数映射的乘积;
根据评价网络输出的及公式(13)单步回报函数,公式(14)的代价函数表示为:
式中,x(k)为当前时刻***状态;为当前估计控制输入;/>为当前时刻的单步回报;γ为折扣因子;x(k+1)为下一刻的***状态;/>为下一时刻动作估计值;/>为下一时刻代价函数估计值;
根据Bellman最优值原理,最优值函数表示为:
式中,x(k)为***状态;u(k)为***输入;r(x(k),u(k))为当前时刻单步回报值;γ为折扣因子;J*(k+1)为下一时刻代价函数最优值;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数;
最优控制u*表示为:
式中,x(k)为***状态;;u(k)为***输入;J*(x(k),u(k))为代价函数最优值。
44)设计评价网络误差函数,评价网络误差函数用来修正评价权值,设计为累积代价误差和终端代价误差之和,根据公式(20)累积代价误差表示为:
基于模糊小脑模型神经网络的终端评价网络表示为:
式中,wc,N(k)为终端评价网络的权值矩阵;aN(k)为终端状态估计值在高斯隶属度函数映射的乘积;
终端代价表示为:
式中,为终端网络估计值;/>为终端状态估计值,根据初始状态和待优化控制序列确定;/>为终端约束;根据终端约束(12a)得到/>为终端状态估计值;QN为权值矩阵;根据公式(23)与(25),得到评价网络的误差函数:
式中,ec_N为终端误差;ec_J为累积误差;Ec为评价网络误差代价函数;
评价网络误差wc(k)更新规则如下:
式中,lc(k)为当前时刻评价网络学习率;Ec(k)为当前时刻评价网络的误差函数;wc(k)为相应的评价网络权值;λi(k)为评价网络隶属度函数向量;ec_N为终端误差;ec_J为累积误差;C为模糊小脑模型神经网络存储单元个数;为评价网络误差;
45)在每个有限时域内设置最大的取值迭代次数nmax,当迭代次数到达nmax,便停止当前预测时域内的策略优化,将所得的策略作用到时间的***当中去;在每个预测时域内的初始权值可以选择为上一预测时域已经收敛的权值。
本发明的有益效果为:
本发明采用滚动时域强化学习方法解决了长预测时域控制求解时间问题,该算法充分利用强化学习的求解力,并融合了预测控制解决约束及扰动的优势,在线求解非线性时变最优化问题,本发明在能量效率和计算效率方面具有明显优势,具有工程应用能力。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明的整体***结构示意图;
图2为本发明中步骤四的结构示意图;
图3为模糊小脑模型神经网络的结构示意图;
图4为速度在距离域的对比曲线图;
图5为多个路口下的交通灯时序及速度曲线图;
图6为电动汽车的能量消耗对比曲线图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
参阅图1-图3,一种基于智能网联信息的整车经济性速度规划方法,包括以下步骤:
步骤一、通过智能网联环境获得路段上的动态交通信息,通过传感器获得本车的状态,确定车辆的安全边界即最大的行驶速度和最大的加速度;
所述动态交通信息包括外界交通灯、限速及前车信息。
步骤二、建立车辆运动学模型,并且分析车辆行驶过程中的能量消耗,选择预测时域N内的燃料消耗作为目标函数,目标函数中同时考虑的车辆行驶过程中的动力性和安全性;具体如下:
在上层经济性行驶速度轨迹规划中,只考虑车辆的纵向运动学,运动学模型如下:
式中,d为行驶距离;v为车速;m为整车质量;δ为旋转质量系数;Ft为总驱动力;FR为行驶阻力,包括空气阻力,滚动阻力和坡度阻力,表示如下:
式中,fw为空气阻力;ff为滚动阻力;fi为坡度阻力;CD为阻力系数;A为迎风面积;v为车速;fr为滚动阻力系数;m为整车质量;θ为坡度;
车辆在加速行驶过程中,驱动力由轮毂电机提供,车辆在制动过程中,总驱动力Ft等于0,制动力由液压制动提供(暂不考虑能量回收),表达式如下:
式中,Tmi是每个电机的驱动力矩i=1,2,3,4,Tm1、Tm2、Tm3、Tm4分别为左前轮、右前轮、左后轮、右后轮,rw为轮胎半径;
能源消耗总量为:
P=Ftv=(mδa+FR)v (5)
式中,Ft为总驱动力;v为车速;m为整车质量;δ为旋转质量系数;a为加速度;FR为行驶阻力。
控制方法用于速度规划。通过寻找最优的车辆速度减少车辆因交通灯等信息引起的频繁加速和制动操作,进而达到最小化整车能耗的目的。
步骤三、根据动态交通信息设定***的性能指标和约束条件;
31)将车辆的行驶距离d和车速v作为状态变量,将能耗公式(1)—(2)进行离散表示如下:
v(k+1)=v(k)+a(k)Δt (7)
式中,k为当前时刻;Δt为***采样时间;v(k)为当前时刻的车速;a(k)为***的加速度,经济性行驶被描述为一个加速度的优化问题;状态变量:x=[x1,x2]=[d,v]T,控制输入:u=[a],式中,d为行驶距离;v为车速;a为加速度;
32)动力学***模型表示为:
x(k+1)=x(k)+f(x(k),u(k)) (8)
式中,x(k)为***状态;u(k)为控制输入;f(x(k),u(k))为由控制输入引起的状态变化;
33)动力学***的性能指标函数为:
式中,k+N是终端约束项;L为电动汽车能量消耗,
L(x,u,k)=β1P(x,u,k)+β2||u(k)-u(k-1)||2 (10)
式中,P(x,u,k)由公式(5)计算得到;x为***状态;u为控制输入;k表示当前时刻;||u(k)-u(k-1)||2为***控制输入的幅值约束,防止过大的加速度,导致舒适性下降,u(k)为***的控制输入;u(k-1)为上一时刻***控制输入;β1,β2为权重因子;T为设定预测时间长度;在性能指标公式(9)中,预测时间步长度N=T/Δt,且随着运行时间递减;车辆行驶的终端时间根据前方交通灯时序确定,即终端时间kf是固定值;每个时刻下的预测时间长度T=kf-k;式中,kf为终端时间;k表示当前时刻;
φ[x(k+N)]为***的终端状态约束,定义为:
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (11)
其中,kf为终端时间;sf为终端距离,由路段上交通灯时序及位置决定;vf为终端速度;a1、a2为权重系数;v(kf)为终端时刻车速;s(kf)为终端时刻位置;
34)根据动态交通信息中的交通灯的实际情况,转化为规划问题的约束条件;当车辆接近交通灯时,具有三种情况,即加速、巡航和减速;根据道路要求1设置车辆自适应巡航初始速度为vd,而车辆在时域末端的速度vf根据是否通过红绿灯确定;具体如下:
a、当车辆在低密度路况且无前车的情况下,考虑前方红绿灯信号作用,将红绿灯信号的时间及距离转化为速度约束,当车辆保持当前车速不能通过下一个红绿灯路口时,需要选择提前加速或者减速达到节能通过红绿灯或停车的目的;通过与红绿灯距离Δd和绿灯所剩时间tgreen进行判断,当车速提高到vmax时,若仍然不能在绿灯时刻通过路口,则选择提前减速,车辆准备在下一个绿灯周期通过路口,即终端速度约束vf=0;反之,车辆定速巡航或加速通过路口,终端速度约束vf=vd;在自适应巡航控制中,设计速度浮动参数ε,即车辆的容许范围∈[1-ε,1+ε];车辆的速度上边界为:vmax=min(vd(1+ε),vlim),vlim为道路限速;
b、当与前车的距离宽裕且能满足速度及加速度需求的情况,本车仍然根据剩余绿灯时间进行跟进或者减速度到停止线停车的控制;反之,车辆的速度容许度受前车影响,车辆必须遵守跟车规则;设定前车的行驶速度为vleader,车辆速度上边界表示为:vmax=min(vleader,vd(1+ε),vlim),通过vmax*tgreen≥Δd的关系,判断是否通过路口,即而得到终端速度条件;
35)通过31)、32)、33)和34)得到***需满足的约束条件如下;
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (12a)
vmin(k)≤v(k)≤vmax(k) (12b)
amin(k)≤a(k)≤amax(k) (12c)
v0=vd,s0=0 (12d)
式中,为***的终端状态约束;a1、a2为权重系数;kf为终端时间;v(kf)为终端时刻车速;s(kf)为终端时刻位置;sf为终端距离;v(k)为终端速度约束;vmin(k)为车辆速度的下边界;vmax(k)为车辆速度的上边界;a(k)为车辆加速度;amin(k)为车辆加速度的下边界;amax(k)为车辆加速度的上边界;v0为车辆初始速度,根据自适应控制车速设定;vd为终端速度约束;s0为车辆的初始位置,通常选择为0。
步骤四、基于滚动时域强化学习,在线求解优化速度即对速度进行规划所述步骤四为:
根据步骤一和步骤三获得的车辆动力学模型、性能指标和约束条件构建优化问题,基于滚动时域强化学习,得到优化后的加速度值,并选择控制序列即加速度值的前Tc(Tc≤N)个控制动作作用于***中,Tc为选择的控制序列长度;车辆根据加速度需求值控制目标车辆驱动力输出,当到达时刻k+Tc时,重复求解控制序列,实现滚动优化控制。具体如下:
41)根据经济性巡航任务要求,设计单步回报函数表示如下:
r(x(k),u(k))=β1P(x(k),u(k))+β2||u(k)-u(k-1)||2 (13)
式中,β1,β2为权重因子;x(k)为***状态;u(k)为控制输入;P(x(k),u(k))为车辆燃料消耗;u(k-1)为***上一时刻控制量;
在有限时域内,代价函数由终端代价和累积回报两部分组成:
式中,EU[·]为在当前策略下的期望算子;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数,根据公式(12a)确定;γ为折扣因子,且0≤γ≤1;k为当前时刻;N为预测时域长度;x(j)为当前时刻状态;u(j)为当前时刻控制输入;
42)设计动作网络,动作网络用于逼近最优控制输入u*(k),采用模糊小脑模型神经网络作为逼近器,模糊小脑模型神经网络包括五个模块,分别为输入变量、变量模糊化、概念映射、实际映射和输出变量,具体如下:
a、输入变量及变量模糊化,将d,v的模糊集均表示为[NB,NS,ZE,PS,PB],其中NB,NS,ZE,PS,PB分别表示为负大、负小、0、正小、正大;两者的论域[-2,-1,0,1,2];隶属度函数向量λj=[λ1j,…λnj]T,其中,j=1,2,n=1,2…5;模糊化变量选择高斯隶属度函数:
式中,ci为第i个隶属度函数中心点,σi为方差;
b、模糊隶属度至输入空间U,将空间U划分为10个存储单元,每个单元对应一个向量,从空间划分的存储单独找到对应地址;
c、输入空间U到Ac的概念映射,输入空间U的10个模糊隶属度函数进行规则划分后得到52个状态,每个状态作为一个指针映射到虚拟空间Ac的C个存储单元中,并找到对应该状态的地址;
d、概念映射Ac到实际映射Ap,采用杂散编码技术中的除留余数法,将概念映射的C个单元映射到Ap的C个单元;
e、实际映射Ap到输出变量模糊小脑模型神经网络的输出表示为:
式中,am(k)为高斯隶属度函数映射的乘积;wm(k)为相应的权值;动作网络的学习目标为极小化动作网络估计值与最优值u*(k)之间的误差值,误差函数为:
权值wa(k)修正迭代公式表示为:
式中,la为动作网络学习率,它随时间的增加而逐渐减小为一个固定值,θ∈(0,1]为网络惯性系数;λi(k)为隶属度函数向量;为前一时刻权值矩阵;/>为(k-2)时刻权值矩阵;
43)设计评价网络,评价网络采用模糊小脑模型神经网络作为逼近器,评价网络的输出表示为:
式中,wc,n(k)为当前时刻权值;an(k)为状态在高斯隶属度函数映射的乘积;
根据评价网络输出的及公式(13)单步回报函数,公式(14)的代价函数表示为:
式中,x(k)为当前时刻***状态;为当前估计控制输入;/>为当前时刻的单步回报;γ为折扣因子;x(k+1)为下一刻的***状态;/>为下一时刻动作估计值;/>为下一时刻代价函数估计值;
根据Bellman最优值原理,最优值函数表示为:
式中,x(k)为***状态;u(k)为***输入;r(x(k),u(k))为当前时刻单步回报值;γ为折扣因子;J*(k+1)为下一时刻代价函数最优值;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数;
最优控制u*表示为:
式中,x(k)为***状态;u(k)为***输入;J*(x(k),u(k))为代价函数最优值。
44)设计评价网络误差函数,评价网络误差函数用来修正评价权值,设计为累积代价误差和终端代价误差之和,根据公式(20)累积代价误差表示为:
基于模糊小脑模型神经网络的终端评价网络表示为:
式中,wc,N(k)为终端评价网络的权值矩阵;aN(k)为终端状态估计值在高斯隶属度函数映射的乘积;
终端代价表示为:
式中,为终端网络估计值;/>为终端状态估计值,根据初始状态和待优化控制序列确定;/>为终端约束;
根据终端约束(12a)得到为终端状态估计值;QN为权值矩阵;根据公式(23)与(25),得到评价网络的误差函数:
式中,ec_N为终端误差;ec_J为累积误差;Ec为评价网络误差代价函数;
评价网络误差wc(k)更新规则如下:
式中,lc(k)为当前时刻评价网络学习率;Ec(k)为当前时刻评价网络的误差函数;wc(k)为相应的评价网络权值;λi(k)为评价网络隶属度函数向量;ec_N为终端误差;ec_J为累积误差;C为模糊小脑模型神经网络存储单元个数;为评价网络误差;
45)在每个有限时域内设置最大的取值迭代次数nmax,当迭代次数到达nmax,便停止当前预测时域内的策略优化,将所得的策略作用到时间的***当中去;在每个预测时域内的初始权值可以选择为上一预测时域已经收敛的权值。
综上所述,本发明首先获得路段前方的交通信息,确定预测时域N及终端约束条件,根据道路要求设定车辆巡航车速vd及目标车辆速度边界和加速度边界值;在时域的起点时刻,采用强化学习的方法来解决受有约束的最优控制问题,思想是通过不断迭代,逐步降低回报函数值,当迭代次数达到设定的最大次数nmax,则停止迭代,输出当前控制量u值,通过上述方法分步求解一个长度为N的最优控制序列。
实施例
基于智能网联汽车硬件在环仿真平台用于对一种基于智能网联信息的整车经济性速度规划方法进行测试。采用某段城市道路信息及交通灯时序构建SCANeR场景,实验路段约3.6km,全程具有路口5个,分别在500m、1200m、1900m、2590m和3320m处。实验中自适应巡航***(ACC)的初始速度为60km/h,离散***的采样时间Δt=1s,滚动时域强化学习的动作网络和评价网络学习率lc(0)…lc(f)为0.3,0.25,0.2…0.05逐渐缩小序列。实验结果如图4-6所示。
速度在距离域的对比曲线如图4所示。相比于ACC控制,在600m和2590m处,MPC(预测控制)算法和RHRL(滚动时域强化学习)算法根据绿灯所剩时间,提前加速通过路口。而ACC算法则在红绿灯处制动停车。MPC算法和RHRL算法相比,RHRL算法采用根据前方信息逐渐提高车速,而MPC速度变化明显大于RHRL算法。在1200处,三种算法均不能通过路口,MPC算法和RHRL算法均采用在容许范围内怠速滑行的方式,减少能耗消耗。图5给出了多个路口下的交通灯时序及速度曲线,可以看出,与基于规则的ACC定速巡航相比,本发明提出的算法和MPC算法均能实现驾驶意图,并在交通路口规划车速,避免频繁的启停,并在节约运行时间方面有较大优势。
在经济性方面,在约3.6km且具有5个交通灯的仿真路段上,三种方法的能量消耗如图6所示。
实验结果表明,本发明可以使车辆尽可能在绿灯时间内通过路口,减少车辆由于红灯而频繁启停,同时增加电机在高效区域的概率和降低轮胎滑移耗散值,从宏观和微观层面共同降低整车的能量消耗。和ACC定速巡航控制相比能耗降低约26.49%,和MPC算法相比,本发明所提的方法在计算过程中,对目标函数进行不断学习优化,故能耗降低约2.31%。
在计算效率方面,本发明采用滚动时域的机制,在每个预测时域内均采用有限时域的强化学***均计算时间代价降低约81.7%。

Claims (4)

1.一种基于智能网联信息的整车经济性速度规划方法,其特征在于,包括以下步骤:
步骤一、通过交通信息层的智能网联环境获得路段上的动态交通信息,并通过传感器获得本车的状态,确定车辆的安全边界即最大的行驶速度和最大的加速度;
步骤二、建立车辆运动学模型,并且分析车辆行驶过程中的能量消耗;
步骤三、根据动态交通信息设定***的性能指标和约束条件;
步骤四、在速度轨迹规划层基于滚动时域强化学习,在线求解优化速度即对速度进行规划,并通过轨迹规划层控制四轮力矩,驱动电动汽车按照需求力矩行驶;
所述步骤二的具体方法如下:
在上层经济性行驶速度轨迹规划中,只考虑车辆的纵向运动学,运动学模型如下:
式中,d为行驶距离;v为车速;m为整车质量;δ为旋转质量系数;Ft为总驱动力;FR为行驶阻力,包括空气阻力,滚动阻力和坡度阻力,表示如下:
式中,fw为空气阻力;ff为滚动阻力;fi为坡度阻力;CD为阻力系数;A为迎风面积;v为车速;fr为滚动阻力系数;m为整车质量;θ为坡度;
车辆在加速行驶过程中,驱动力由轮毂电机提供,车辆在制动过程中,总驱动力Ft等于0,制动力由液压制动提供,表达式如下:
式中,Tmi是每个电机的驱动力矩i=1,2,3,4,Tm1、Tm2、Tm3、Tm4分别为左前轮、右前轮、左后轮、右后轮,rw为轮胎半径;
能源消耗总量为:
P=Ftv=(mδa+FR)v (5)
式中,Ft为总驱动力;v为车速;m为整车质量;δ为旋转质量系数;a为加速度;FR为行驶阻力;
所述步骤三的具体方法如下:
31)将车辆的行驶距离d和车速v作为状态变量,将能耗公式(1)—(2)进行离散表示如下:
v(k+1)=v(k)+a(k)Δt (7)
式中,k为当前时刻;Δt为***采样时间;v(k)为当前时刻的车速;a(k)为***的加速度,经济性行驶被描述为一个加速度的优化问题;状态变量:x=[x1,x2]=[d,v]T,控制输入:u=[a],式中,d为行驶距离;v为车速;a为加速度;
32)运动学模型表示为:
x(k+1)=x(k)+f(x(k),u(k)) (8)
式中,x(k)为***状态;u(k)为控制输入;f(x(k),u(k))为由控制输入引起的状态变化;
33)动力学***的性能指标函数为:
式中,k+N是终端约束项;L为电动汽车能量消耗,
L(x,u,k)=β1P(x,u,k)+β2||u(k)-u(k-1)||2 (10)
式中,P(x,u,k)为车辆燃料消耗,其由公式(5)计算得到;x为***状态;u为控制输入;k表示当前时刻;||u(k)-u(k-1)||2为***控制输入的幅值约束,防止过大的加速度,导致舒适性下降,u(k)为***的控制输入;u(k-1)为上一时刻***控制输入;β1,β2为权重因子;T为设定预测时间长度;在性能指标公式(9)中,预测时间步长度N=T/Δt,且随着运行时间递减;车辆行驶的终端时间根据前方交通灯时序确定,即终端时间kf是固定值;每个时刻下的预测时间长度T=kf-k;式中,kf为终端时间;k表示当前时刻;
φ[x(k+N)]为***的终端状态约束,定义为:
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (11)
其中,kf为终端时间;sf为终端距离,由路段上交通灯时序及位置决定;vf为终端速度;a1、a2为权重系数;v(kf)为终端时刻车速;s(kf)为终端时刻位置;
34)根据动态交通信息中的交通灯的实际情况,转化为规划问题的约束条件;当车辆接近交通灯时,具有三种情况,即加速、巡航和减速;根据道路要求1设置车辆自适应巡航初始速度为vd,而车辆在时域末端的速度vf根据是否通过红绿灯确定;具体如下:
a、当车辆在低密度路况且无前车的情况下,考虑前方红绿灯信号作用,将红绿灯信号的时间及距离转化为速度约束,当车辆保持当前车速不能通过下一个红绿灯路口时,需要选择提前加速或者减速达到节能通过红绿灯或停车的目的;通过与红绿灯距离Δd和绿灯所剩时间tgreen进行判断,当车速提高到vmax时,若仍然不能在绿灯时刻通过路口,则选择提前减速,车辆准备在下一个绿灯周期通过路口,即终端速度约束vf=0;反之,车辆定速巡航或加速通过路口,终端速度约束vf=vd;在自适应巡航控制中,设计速度浮动参数ε,即车辆的容许范围∈[1-ε,1+ε];车辆的速度上边界为:vmax=min(vd(1+ε),vlim),vlim为道路限速;
b、当与前车的距离宽裕且能满足速度及加速度需求的情况,本车仍然根据剩余绿灯时间进行跟进或者减速度到停止线停车的控制;反之,车辆的速度容许度受前车影响,车辆必须遵守跟车规则;设定前车的行驶速度为vleader,车辆速度上边界表示为:vmax=min(vleader,vd(1+ε),vlim),通过vmax*tgreen≥Δd的关系,判断是否通过路口,即而得到终端速度条件;
35)通过31)、32)、33)和34)得到***需满足的约束条件如下;
φ[x(k+N)]=a1(v(kf)-vf)2+a2(s(kf)-sf)2 (12a)
vmin(k)≤v(k)≤vmax(k) (12b)
amin(k)≤a(k)≤amax(k) (12c)
v0=vd,s0=0 (12d)
式中,为***的终端状态约束;a1、a2为权重系数;kf为终端时间;v(kf)为终端时刻车速;s(kf)为终端时刻位置;sf为终端距离;v(k)为终端速度约束;vmin(k)为车辆速度的下边界;vmax(k)为车辆速度的上边界;a(k)为车辆加速度;amin(k)为车辆加速度的下边界;amax(k)为车辆加速度的上边界;v0为车辆初始速度,根据自适应控制车速设定;vd为终端速度约束;s0为车辆的初始位置,通常选择为0。
2.根据权利要求1所述的一种基于智能网联信息的整车经济性速度规划方法,其特征在于,步骤一中所述动态交通信息包括外界交通灯位置和时序、限速及前车信息,所述动态交通信息通过车载传感器获得。
3.根据权利要求1所述的一种基于智能网联信息的整车经济性速度规划方法,其特征在于,所述步骤四为:
根据步骤一和步骤三获得的车辆运动学模型、性能指标和约束条件构建优化问题,基于滚动时域强化学习,得到优化后的加速度值,并选择控制序列即加速度值的前Tc(Tc≤N)个控制动作作用于***中,Tc为选择的控制序列长度;车辆根据加速度需求值控制目标车辆驱动力输出,当到达时刻k+Tc时,重复求解控制序列,实现滚动优化控制。
4.根据权利要求3所述的一种基于智能网联信息的整车经济性速度规划方法,其特征在于,所述步骤四的具体方法如下:
41)根据经济性巡航任务要求,设计单步回报函数表示如下:
r(x(k),u(k))=β1P(x(k),u(k))+β2||u(k)-u(k-1)||2 (13)
式中,β1,β2为权重因子;x(k)为***状态;u(k)为控制输入;P(x(k),u(k))为车辆燃料消耗;u(k-1)为***上一时刻控制量;
在有限时域内,代价函数由终端代价和累积回报两部分组成:
式中,EU[·]为在当前策略下的期望算子;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数,根据公式(12a)确定;γ为折扣因子,且0≤γ≤1;k表示当前时刻;N为预测时域长度;x(j)为当前时刻状态;u(j)为当前时刻控制输入;
42)设计动作网络,动作网络用于逼近最优控制输入u*(k),采用模糊小脑模型神经网络作为逼近器,模糊小脑模型神经网络包括五个模块,分别为输入变量、变量模糊化、概念映射、实际映射和输出变量,具体如下:
a、输入变量及变量模糊化,将d,v的模糊集均表示为[NB,NS,ZE,PS,PB],其中NB,NS,ZE,PS,PB分别表示为负大、负小、0、正小、正大;两者的论域[-2,-1,0,1,2];隶属度函数向量λj=[λ1j,…λnj]T,其中,j=1,2,n=1,2…5;模糊化变量选择高斯隶属度函数:
式中,ci为第i个隶属度函数中心点,σi为方差;
b、模糊隶属度至输入空间U,将空间U划分为10个存储单元,每个单元对应一个向量,从空间划分的存储单独找到对应地址;
c、输入空间U到Ac的概念映射,输入空间U的10个模糊隶属度函数进行规则划分后得到52个状态,每个状态作为一个指针映射到虚拟空间Ac的C个存储单元中,并找到对应该状态的地址;
d、概念映射Ac到实际映射Ap,采用杂散编码技术中的除留余数法,将概念映射的C个单元映射到Ap的C个单元;
e、实际映射Ap到输出变量模糊小脑模型神经网络的输出表示为:
式中,am(k)为高斯隶属度函数映射的乘积;wm(k)为相应的权值;动作网络的学习目标为极小化动作网络估计值与最优值u*(k)之间的误差值,误差函数为:
权值wa(k)修正迭代公式表示为:
式中,la为动作网络学习率,它随时间的增加而逐渐减小为一个固定值,θ∈(0,1]为网络惯性系数;λi(k)为隶属度函数向量;为前一时刻权值矩阵;/>为(k-2)时刻权值矩阵;
43)设计评价网络,评价网络采用模糊小脑模型神经网络作为逼近器,评价网络的输出表示为:
式中,wc,n(k)为当前时刻权值;an(k)为状态在高斯隶属度函数映射的乘积;
根据评价网络输出的及公式(13)单步回报函数,公式(14)的代价函数表示为:
式中,x(k)为当前时刻***状态;为当前估计控制输入;/>为当前时刻的单步回报;γ为折扣因子;x(k+1)为下一刻的***状态;/>为下一时刻动作估计值;为下一时刻代价函数估计值;
根据Bellman最优值原理,最优值函数表示为:
式中,x(k)为***状态;u(k)为***输入;r(x(k),u(k))为当前时刻单步回报值;γ为折扣因子;J*(k+1)为下一时刻代价函数最优值;x(k+N)为终端状态;φ[x(k+N)]为终端代价函数;
最优控制u*表示为:
式中,x(k)为***状态;u(k)为***输入;J*(x(k),u(k))为代价函数最优值;
44)设计评价网络误差函数,评价网络误差函数用来修正评价权值,设计为累积代价误差和终端代价误差之和,根据公式(20)累积代价误差表示为:
基于模糊小脑模型神经网络的终端评价网络表示为:
式中,wc,N(k)为终端评价网络的权值矩阵;aN(k)为终端状态估计值在高斯隶属度函数映射的乘积;
终端代价表示为:
式中,为终端网络估计值;/>为终端状态估计值,根据初始状态和待优化控制序列确定;/>为终端约束;根据终端约束(12a)得到/> 为终端状态估计值;QN为权值矩阵;根据公式(23)与(25),得到评价网络的误差函数:
式中,ec_N为终端误差;ec_J为累积误差;Ec为评价网络误差代价函数;
评价网络误差wc(k)更新规则如下:
式中,lc(k)为当前时刻评价网络学习率;Ec(k)为当前时刻评价网络的误差函数;wc(k)为相应的评价网络权值;λi(k)为评价网络隶属度函数向量;ec_N为终端误差;ec_J为累积误差;C为模糊小脑模型神经网络存储单元个数;为评价网络误差;
45)在每个有限时域内设置最大的取值迭代次数nmax,当迭代次数到达nmax,便停止当前预测时域内的策略优化,将所得的策略作用到时间的***当中去;在每个预测时域内的初始权值选择为上一预测时域已经收敛的权值。
CN202111087524.5A 2021-09-16 2021-09-16 一种基于智能网联信息的整车经济性速度规划方法 Active CN113741199B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111087524.5A CN113741199B (zh) 2021-09-16 2021-09-16 一种基于智能网联信息的整车经济性速度规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111087524.5A CN113741199B (zh) 2021-09-16 2021-09-16 一种基于智能网联信息的整车经济性速度规划方法

Publications (2)

Publication Number Publication Date
CN113741199A CN113741199A (zh) 2021-12-03
CN113741199B true CN113741199B (zh) 2024-04-26

Family

ID=78739357

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111087524.5A Active CN113741199B (zh) 2021-09-16 2021-09-16 一种基于智能网联信息的整车经济性速度规划方法

Country Status (1)

Country Link
CN (1) CN113741199B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114684165B (zh) * 2022-03-04 2023-03-24 阿波罗智能技术(北京)有限公司 车辆控制方法、装置、设备、存储介质和自动驾驶车辆
CN114859893A (zh) * 2022-04-06 2022-08-05 吉林大学 智能网联电动汽车在不同交通状态下的节能学习控制方法
CN114475244B (zh) * 2022-04-14 2022-09-23 禾多科技(北京)有限公司 定速巡航控制方法、整车控制器及定速巡航控制***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109591659A (zh) * 2019-01-14 2019-04-09 吉林大学 一种智能学习的纯电动汽车能量管理控制方法
CN110264757A (zh) * 2019-05-21 2019-09-20 吉林大学 基于连续信号灯信息的智能网联汽车分层速度规划方法
CN111791887A (zh) * 2020-07-03 2020-10-20 北京理工大学 一种基于分层式车速规划的车辆节能驾驶方法
CN112896161A (zh) * 2021-02-08 2021-06-04 杭州电子科技大学 一种基于强化学习的电动汽车生态自适应巡航控制***
CN113192347A (zh) * 2021-03-16 2021-07-30 东南大学 动态交通场景下智能网联车辆经济性驾驶控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110936824B (zh) * 2019-12-09 2021-06-04 江西理工大学 一种基于自适应动态规划的电动汽车双电机控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109591659A (zh) * 2019-01-14 2019-04-09 吉林大学 一种智能学习的纯电动汽车能量管理控制方法
CN110264757A (zh) * 2019-05-21 2019-09-20 吉林大学 基于连续信号灯信息的智能网联汽车分层速度规划方法
CN111791887A (zh) * 2020-07-03 2020-10-20 北京理工大学 一种基于分层式车速规划的车辆节能驾驶方法
CN112896161A (zh) * 2021-02-08 2021-06-04 杭州电子科技大学 一种基于强化学习的电动汽车生态自适应巡航控制***
CN113192347A (zh) * 2021-03-16 2021-07-30 东南大学 动态交通场景下智能网联车辆经济性驾驶控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
未知环境下智能汽车轨迹规划滚动窗口优化;张琳;章新杰;郭孔辉;王超;刘洋;刘涛;;吉林大学学报(工学版);20171116(第03期);全文 *
生态驾驶在智能网联汽车中的应用综述;吴官朴;杨昱;赵阳;马芳武;;汽车文摘;20200630(第07期);全文 *

Also Published As

Publication number Publication date
CN113741199A (zh) 2021-12-03

Similar Documents

Publication Publication Date Title
CN113741199B (zh) 一种基于智能网联信息的整车经济性速度规划方法
CN106740846B (zh) 一种双模式切换的电动汽车自适应巡航控制方法
He et al. An improved MPC-based energy management strategy for hybrid vehicles using V2V and V2I communications
CN111845701B (zh) 一种跟车环境下基于深度强化学习的hev能量管理方法
CN107688343B (zh) 一种混合动力车辆的能量控制方法
CN110936824B (zh) 一种基于自适应动态规划的电动汽车双电机控制方法
CN112668799A (zh) 基于行驶大数据的phev的智能能量管理方法和存储介质
CN111267831A (zh) 一种混合动力车辆智能变时域模型预测能量管理方法
Lü et al. Hybrid electric vehicles: A review of energy management strategies based on model predictive control
CN113788021B (zh) 一种结合前车速度预测的自适应跟车巡航控制方法
CN112389436B (zh) 基于改进lstm神经网络的安全性自动驾驶换道轨迹规划方法
CN111923897B (zh) 一种插电式混合动力汽车能量智能管理方法
CN113635879B (zh) 一种车辆制动力分配方法
CN115257724A (zh) 一种插电式混动汽车安全节能决策控制方法和***
CN108313057B (zh) 基于mpc和凸优化算法的纯电动汽车自适应巡航控制方法
CN112896161A (zh) 一种基于强化学习的电动汽车生态自适应巡航控制***
CN113276829B (zh) 一种基于工况预测的车辆行驶节能优化变权重方法
CN108749816B (zh) 运用能量耗散理论进行智能车辆速度调控的方法
Yang et al. Variable optimization domain-based cooperative energy management strategy for connected plug-in hybrid electric vehicles
Zhang et al. An optimal vehicle speed planning algorithm for regenerative braking at traffic lights intersections based on reinforcement learning
CN113232651B (zh) 一种基于adhdp的燃油汽车生态自适应巡航控制方法
CN113561976B (zh) 基于反馈优化的车辆节能预测自适应巡航控制方法和装置
Zhou et al. Energy optimization for intelligent hybrid electric vehicles based on hybrid system approach in a car‐following process
CN114523967B (zh) 一种基于神经网络的预见巡航控制方法
Pan et al. Energy-optimized adaptive cruise control strategy design at intersection for electric vehicles based on speed planning

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