CN103927895A - 一种基于车路/车车通讯的车辆弯道通行辅助系统 - Google Patents

一种基于车路/车车通讯的车辆弯道通行辅助系统 Download PDF

Info

Publication number
CN103927895A
CN103927895A CN201410169875.4A CN201410169875A CN103927895A CN 103927895 A CN103927895 A CN 103927895A CN 201410169875 A CN201410169875 A CN 201410169875A CN 103927895 A CN103927895 A CN 103927895A
Authority
CN
China
Prior art keywords
vehicle
bend
sequence
average
gamma
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
CN201410169875.4A
Other languages
English (en)
Other versions
CN103927895B (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201410169875.4A priority Critical patent/CN103927895B/zh
Publication of CN103927895A publication Critical patent/CN103927895A/zh
Application granted granted Critical
Publication of CN103927895B publication Critical patent/CN103927895B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

本发明涉及一种基于车路/车车通讯的车辆弯道通行辅助系统,它包括若干驶过弯道的车辆、信息通讯和处理基站和即将驶入弯道的车辆;每一驶过弯道的车辆上设置有车载传感系统和车载通讯模块,即将驶入弯道的车辆上也设置有车载通讯模块,信息通讯和处理基站设置在弯道一侧;车载传感系统采集驶过弯道的车辆的行驶状态信息并通过车载通讯模块发送到路侧通讯设备,路侧通讯设备将接收的驶过弯道的车辆行驶状态信息传送给到路侧存储设备,路侧信息处理单元根据不同车辆类型将数据进行分类,并对同一类型的车辆行驶状态数据进行无监督式聚类和统计学习,得到最优速度序列和轨迹序列,并将得其传送给即将驶入弯道的车辆,即将驶入弯道的车辆根据车辆类型通过最优速度序列和轨迹序列对过弯道进行车速建议或辅助控制。

Description

一种基于车路/车车通讯的车辆弯道通行辅助系统
技术领域
本发明涉及智能交通系统领域,特别是关于一种基于车路/车车通讯的车辆弯道通行辅助系统。
背景技术
随着汽车保有量增加和道路交通流量上升,道路堵塞、环境污染和交通事故等问题日益严重。在汽车行驶的路段中,弯道路段工况复杂,是行车过程中的危险路段。车辆在弯道路段行驶时,超速行驶是发生交通事故的主要原因。驾驶员难以及时准确地估计弯道的行车状况,如弯道曲率半径、道路摩擦系数、交通状况等信息,导致以过高的车速进入弯道,容易发生侧翻、驶出车道的交通事故。限制车辆在弯道路段的行驶速度,将车速控制在安全阈值,可以显著减小弯道交通事故的发生率和事故严重性。为此,交通管理部门依据不同的道路类型和道路工况制定了相应的限速标准,在道路相应路段设置限速标志牌,但多数驾驶员很少严格执行这些标准,且该限速标准也难以根据天气条件进行动态的调整。由此可看出,作为人-车-路闭环交通系统中重要的组成部分,驾驶员由于自身条件的限制,已成为这个系统中的薄弱环节。基于智能交通系统(ITS)的驾驶辅助系统(Advanced Driving Assistance Systems/ADAS)可以提升“人”对“车-路”系统的感知、决策和控制能力。驾驶辅助系统将人-车-路作为一个闭环系统,利用现代信息技术、传感技术来扩展驾驶员的环境感知能力,在路况与车况的综合信息中辨识安全隐患。因此在危险工况中,ADAS不仅为驾驶员提供危险预警信息,而且能够自动采取措施增加车辆的安全性。影响弯道安全车速的因素主要有:弯道曲率半径、道路附着系数、能见度、道路横向坡度等。早期发展的弯道安全辅助系统有车辆侧翻预警系统、车辆防侧滑预警系统、车辆电子稳定性控制系统等,这些辅助系统都依赖于对车辆和道路状态的实时感知,由于感知的准确性和实时性难以保证,始终没有特别有效的辅助系统保障车辆在弯道上行驶的安全性。
V2V(Vehicle-to-Vehicle)和V2I(Vehicle-to-Infrastructure)通讯技术的发展使得ITS技术可以有效地将人-车-路之间的信息有机地整合起来,驾驶辅助系统感知的信息内容越来越丰富、感知的范围越来越宽、感知的实时性越来越高,从而更加有效地提高车辆的行驶安全性、舒适性和经济性。因此,基于V2V和V2I通讯的协同式驾驶安全辅助系统(Cooperative Advanced Driver Assistance Systems)得到了快速发展和应用。目前,应用于弯道的驾驶安全辅助系统的基本原理是:利用路侧装置识别弯道处的路况(如交通事故)、道路条件(如路面潮湿程度)、天气条件(如雨雪)等,并将这些信息通过V2I技术发送给即将驶入弯道的车辆,使得这些车辆提前获取弯道信息,并计算出合理的通过弯道的车速,以保证行车安全。
上述现有技术存在的问题主要包括几个方面:1)获取弯道处交通状况、天气条件、道路条件等信息的手段不够成熟,获取的信息精度不高,并且实时性难以保证;2)获取弯道信息需要高精度的传感系统,成本高,难以推广应用;3)目前能够获取的信息大多是相对静态的,如弯道曲率、道路坡度等,但是对于一些动态信息如路面附着系数、弯道交通事故、天气条件等尚难以获得;4)目前弯道安全通行车速模型的建立多依据静态道路信息,难以反应动态的道路条件变化;5)目前基于车车/车路通讯的弯道通行辅助系统是从车辆动力学控制的角度建立的,没有基于历史行车数据的自学习功能,智能化程度低。
发明内容
针对上述问题,本发明的目的是提供一种基于驶过弯道的车辆的历史行车状态数据,通过无监督式学习方法聚类区分出激进、普通、保守驾驶数据集,然后采用带时间因子的加权最小二乘方法得到统计意义下最优的行车速度序列和轨迹序列,从而为不同类型的车辆通过弯道提供合理的车速建议或者辅助控制的基于车路/车车通讯的车辆弯道通行辅助系统。
为实现上述目的,本发明采取以下技术方案:一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:它包括若干驶过弯道的车辆、一信息通讯和处理基站和一即将驶入弯道的车辆;其中,每一所述驶过弯道的车辆上设置有一车载传感系统和一车载通讯模块,所述即将驶入弯道的车辆上也设置有一车载通讯模块,所述信息通讯和处理基站设置在弯道一侧,它包括一路侧通讯设备、一路侧存储设备和一路侧信息处理单元;所述车载传感系统采集所述驶过弯道的车辆的行驶状态信息并通过自车的车载通讯模块发送到所述路侧通讯设备,所述路侧通讯设备将接收的驶过弯道的车辆行驶状态信息传送给到所述路侧存储设备进行存储,所述路侧信息处理单元获取所述路侧存储设备中的数据并根据不同车辆类型将数据进行分类,并对同一类型的车辆行驶状态数据进行无监督式聚类和统计学习,得到最优速度序列和轨迹序列,所述路侧信息处理单元将得到的最优速度序列和轨迹序列发送到所述路侧存储设备进行存储,所述路侧通讯设备获取所述路侧存储设备中的最优速度序列和轨迹序列通过所述车载通讯模块传送给所述即将驶入弯道的车辆,所述即将驶入弯道的车辆根据车辆类型通过最优速度序列和轨迹序列对过弯道进行车速建议或辅助控制,使得车辆能够安全、畅通地通过弯道。
所述路侧通讯设备采用DSRC通信设备,所述DSRC通信设备采用专用短程通信技术。
所述路侧信息处理单元包括一数据预处理模块、一无监督式学习分类模块和一统计学习最优速度序列及轨迹序列模块;所述数据预处理模块用于对接收的所述驶过弯道车辆行驶状态数据进行Kalman滤波,滤除噪声;所述无监督式学习分类模块用于采用无监督式学习方法判断个体车辆中的驾驶员类型,以确定后续统计学习的权重;所述统计学习最优速度序列及轨迹序列模块用于对每一类型车辆采用统计学习的方法,学习得到统计意义上的最优速度序列和轨迹序列。
所述无监督式学习分类模块的工作过程为:1)根据滤波后的各序列数据,对于同一类型的车辆行驶状态数据序列提取其特殊数据点,并将上述特殊数据点的值组成一个向量,用于表征车辆通过弯道的行驶特征:
xT,i=[vs,i,Ls,i,Pm,i,ax,i,vI,i,vO,i,vmax,i,vmin,i,ay,imax,imax,imax,i]T    (1)
式中,vs,i为开始制动时的车速、Ls,i为开始制动时相对弯道入口的距离、Pm,i为最大制动强度、ax,i为最大制动减速度、vI,i为进入弯道的速度、vO,i为离开弯道的速度、vmax,i为通过弯道过程中的最大速度、vmin,i为通过弯道过程中的最小速度、ay,i为通过弯道过程中的最大横向加速度、ωmax,i为通过弯道过程中的最大横摆角速度、Δmax,i为过弯轨迹与实际轨迹的最大差别量、λmax,i为四个车轮出现的最大滑移率;2)对已经驶过弯道的同一车辆类型的N辆车行驶状态的特征向量进行k-means聚类,聚成三类;聚类的目标是使得式(2)的目标函数最小,式(2)中x为每个样本的特征向量,每个样本代表驶过弯道的一个车辆的状态数据,样本的特征向量x即为表征车辆通过弯道的行驶特征xT,i,Γj∈{Γ123}表示聚成的三类数据集,mj是每类数据集的均值即中心点:
J ( e ) = Σ j = 1 3 Σ x ∈ Γ j | | x - m j | | 2 - - - ( 2 )
式(2)中,J(e)表示聚成三类数据集的类内离散程度;3)统计每一类数据集中的每个样本最大纵向加速度ax,i的均值ax,averagej)、最大横向加速度ay,i的均值ay,averagej)、最大横摆角速度ωmax,i的均值以及最大过弯轨迹偏差绝对值Δmax,i的均值计算公式如下:
a x , average ( Γ j ) = 1 N j Σ i ∈ Γ j a x , i a y , average ( Γ j ) = 1 N j Σ i ∈ Γ j a y , i ω average ( Γ j ) = 1 N j Σ i ∈ Γ j ω max , i Δ average ( Γ j ) = 1 N j Σ i ∈ Γ j Δ max , i j = 1,2,3 - - - ( 3 )
式(3)中,{Γj}j=1,2,3为无监督k-means聚类的结果,Nj为子数据集Γj的样本个数,i∈Γj表示属于类别Γj的样本;4)计算聚类三个子数据集的驾驶员类型因子Index(Γj),将驾驶员类型因子{Index(Γj)}j=1,2,3由大到小排序,依次将对应的无监督的聚类结果表示为激进、普通和保守三种驾驶员类型,其中,驾驶员类型因子Index(Γj)的计算过程为:
Index ( Γ j ) = | a x , average ( Γ j ) | max { | a x , average ( Γ j ) | } + | a y , average ( Γ j ) | max { | a y , average ( Γ j ) | } + | ω average ( Γ j ) | max { | ω average ( Γ j ) | } + | Δ average ( Γ j ) | max { | Δ average ( Γ j ) | } , j = 1,2,3 .
所述统计学习最优速度序列及轨迹序列模块对于速度序列的统计学习方法如下:1)将速度时间序列和轨迹时间序列离散为弯道的空间序列,即将不同车辆的速度-时间序列转化为速度-距离序列,并将行车轨迹序列转化为轨迹-距离序列;2)对速度-距离序列采用带时间因子的加权最小二乘估计方法得到最优的速度序列;3)对轨迹-距离序列采用带时间因子的加权最小二乘估计方法得到最优的轨迹序列。
本发明由于采取以上技术方案,其具有以下优点:1、本发明通过车载传感系统采集驶过弯道的车辆的行驶状态信息并通过自车的车载通讯模块发送到路侧通讯设备,路侧通讯设备将接收的驶过弯道的车辆行驶状态信息传送给到路侧存储设备,路侧信息处理单元获取路侧存储设备的数据根据不同车辆类型将数据进行分类,并对同一类别的车辆行驶状态数据进行无监督式聚类和统计学习,得到最优速度序列和轨迹序列,并将其发送给即将驶入弯道的车辆,即将驶入弯道的车辆根据车辆类型通过最优速度序列和轨迹序列对过弯道进行车速建议或辅助控制,使得车辆能够安全、畅通地通过弯道,因此本发明首先对驶过弯道的车辆状态进行信息采集、处理、存储与共享,并对车辆通过弯道时的行驶特性、驾驶员的驾驶行为进行自动分析、智能学习,将学习结果通过车路/车车通讯技术实时发送给即将驶入弯道的车辆,在不需要采集动态交通状况信息、天气条件信息、路面条件信息等情况下,即可为不同类型的车辆安全、畅通地通过弯道提供合理的车速建议或辅助控制。2、本发明对即将进入弯道车辆的建议车速或是辅助控制由已经通过该弯道的车辆行驶状态辨识而来,具有智能学习的功能,根据通过该路段的车辆状态数据,获得合理的弯道通行车速建议或辅助控制,对即将驶入弯道的车辆的车速辅助方式一方面可以提供速度建议,使驾驶员合理地控制车速,另一方面也可以自动控制车辆,保障行车安全。3、本发明在弯道一侧设置有路侧通讯设备和路侧存储设备,路侧通讯设备可以采用DSRC通信设备,通过与车路通讯模块交互获得通过该路段的车辆状态数据,并具有数据存储、数据分析和安全车速辨识功能,随着记录数据的增加,对通过该路段车辆的辅助就更加智能化,能够为即将通过该路段的不同车型提供合理的通行速度建议,使车辆能够安全、畅通合理地通过弯道。4、本发明的路侧信息处理单元通过对驶过该路段车辆的状态数据分析,可辨识不同类型驾驶员的弯道驾驶特性,作为弯道合理通行速度的辅助依据,具有拟人化的功能;因此不需要对弯道的路况条件、天气条件进行感知,对车辆状态数据的准确性和实时性要求也没有车辆动力学控制要求的那么高,从而降低了系统成本,提高了系统的可靠性和实用性。本发明可以广泛应用于车辆弯道行驶中。
附图说明
以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
图1是本发明的车辆弯道通行辅助系统结构示意图;
图2是本发明的弯道标记示意图;
图3是本发明的学习过程流程图;
图4是本发明的k-means聚类流程示意图;
图5是本发明的时间序列转化为沿弯道中心线位移序列示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1~3所示,本发明的基于车路/车车通讯的车辆弯道通行辅助系统包括若干驶过弯道的车辆1、一信息通讯和处理基站2和一即将驶入弯道的车辆3;其中,每一驶过弯道的车辆1上设置有一车载传感系统4和一车载通讯模块5,即将驶入弯道的车辆3上也设置有一车载通讯模块6;信息通讯和处理基站2设置在弯道一侧,具体位置不作限定,只要保证信息通讯能够覆盖车辆在弯道行驶的区域即可,它包括一路侧通讯设备21、一路侧存储设备22和一路侧信息处理单元23;车载传感系统4采集驶过弯道的车辆1的行驶状态信息并通过自车的车载通讯模块5发送路侧通讯设备21,路侧通讯设备21将接收的驶过弯道的车辆1行驶状态信息传送给路侧存储设备22进行储存,路侧信息处理单元23通过获取路侧存储设备22的数据并根据驶过弯道的不同车辆类型将数据进行分类,并对同一类别的车辆行驶状态数据进行无监督式聚类和统计学习,得到最优速度序列和轨迹序列,路侧信息处理单元23将得到的最优速度序列和轨迹序列发送到路侧存储设备22,路侧通讯设备21获取路侧存储设备中的最优速度序列和轨迹序列通过车载通讯模块6传送给即将驶入弯道的车辆3,即将驶入弯道的车辆3根据所属车辆类型通过最优速度序列和轨迹序列对过弯道进行车速建议或辅助控制,使得车辆能够安全、畅通地通过弯道。
在一个优选的实施例中,路侧通讯设备21可以采用DSRC通信设备,DSRC通信设备采用专用短程通信技术(Dedicated Short Range Communications,DSRC),DSRC可以实现在小区域内与高速运动目标间的双向通信,通过车辆的车载通讯模块和DSRC通信设备能够实现对车辆的行驶状态信息的采集和发送。
在一个优选的实施例中,如图2、图5所示,在对驶过弯道的车辆的行驶状态信息采集时,需要预先对弯道进行标记,具体标记过程为:对于辅助车辆行驶的弯道路段,一般可认为由直线、缓和曲线和圆曲线组成,几何上标记如下的关键点:A为数据零点,即车辆驶过该点时,路侧通讯设备21接收车辆行驶状态数据的时间零点;B为弯道入口,沿车辆行驶的方向,道路曲率半径由无穷大变为刚小于某一设定值处定义为弯道入口;C为弯道出口,沿车辆行驶的方向,道路曲率半径由小于某一设定值变为刚大于此设定值处定义为弯道出口(弯道入口和弯道出口的定义中某一设定值可以根据具体的道路条件进行确定,本发明试验时此设定值取为2000m,但是不限于此,可以根据具体情况确定);D为数据终点,即车辆驶过该点之后,路侧通讯设备21将不接受车辆的行驶状态数据,上述标记点均在路侧通讯设备21的DSRC通信覆盖范围内。另外,从数据零点A到数据终点D路段的道路中心轨迹坐标在信息通讯和处理基站2上是已知的。
在一个优选的实施例中,驶过弯道的车辆1的车载传感系统4包括有GPS定位仪和车载自主式传感器,GPS定位仪用于对车辆的位置坐标进行定位,车载自主式传感器包括设置在车辆上的车速传感器、轮速传感器、纵/横加速度传感器、横摆角速度传感器、制动压力传感器、踏板开度传感器、方向盘转角传感器,由于上述各传感器均为现有车辆中常用传感器因此对其不作赘述;采用上述车载传感系统4采集驶过弯道的车辆1的行驶状态信息包括:车辆类型T,时间戳Time_stamp、驶过弯道时车辆的速度序列vT,i(tk)、位置序列{xT,i(tk),yT,i(tk)}、纵向加速度序列aTx,i(tk)、横向加速度序列aTy,i(tk)、横摆角速度序列ωT,i(tk)、四个车轮的轮速序列ΩT1,i(tk),ΩT2,i(tk),ΩT3,i(tk),ΩT4,i(tk)、制动压力序列PT,i(tk)、加速踏板序列αT,i(tk)、方向盘转角序列θT,i(tk),符号下标T表示车辆的类型,取值如下表所示:
表1
其中,时间戳Time_stamp表示当前通信数据的绝对时间,符号下标i表示在每一类型中通过弯道的车辆标号,该标号在路侧存储设备22中进行记载并区分;tk表示弯道行驶的时间序列,序列间隔为通信频率的倒数,时间序列以车辆经过弯道入口A时为零点,该标记点与路侧通讯设备21的距离保证在DSRC的通信范围内。速度序列vT,i(tk)、位置序列{xT,i(tk),yT,i(tk)}和纵向加速度序列aTx,i(tk)用以表征车辆在弯道的行驶状态,是用于学习的主要状态量;制动压力序列PT,i(tk)、加速踏板序列αT,i(tk)和方向盘转角序列θT,i(tk)用以表征驾驶员的驾驶行为;横向加速度序列aTy,i(tk)、横摆角速度序列ωT,i(tk)、四个车轮的轮速序列ΩT1,i(tk),ΩT2,i(tk),ΩT3,i(tk),ΩT4,i(tk)、车辆过弯的轨迹{xT,i(tk),yT,i(tk)}和弯道中心轨迹的横向偏差量ΔT,i(tk)是用以评价车辆驶过弯道的安全和舒适程度的基本参数,后两组状态用于判断驾驶员的类型:激进、普通和保守,以确定该车辆状态信息在学习过程中的权重。
在一个优选的实施例中,如图1和图3所示,路侧信息处理单元23包括一数据预处理模块、一无监督式学习分类模块和一统计学习最优速度序列及轨迹序列模块;数据预处理模块用于对接收的若干驶过弯道的车辆1的行驶状态数据进行Kalman滤波,滤除噪声,Kalman滤波为现有的滤波方法,具体滤波过程在此不再赘述;无监督式学习分类模块用于采用无监督式学习方法判断个体车辆中的驾驶员类型,以确定后续统计学习的权重;统计学习最优速度序列及轨迹序列模块用于对每一类型车辆采用统计学习的方法,学习得到统计意义下的最优速度序列和轨迹序列。
其中,无监督式学习分类模块的工作过程为:
1)根据滤波后的各序列数据,对于同一类型(I为普通轿车、II为普通客车、III为普通货车)的车辆行驶状态数据序列提取其特殊数据点,特殊数据点包括:开始制动时的车速vs,i、开始制动时相对弯道入口的距离Ls,i、最大制动强度Pm,i、最大制动减速度ax,i、进入弯道的速度vI,i、离开弯道的速度vO,i、通过弯道过程中的最大速度vmax,i、通过弯道过程中的最小速度vmin,i、通过弯道过程中的最大横向加速度ay,i、通过弯道过程中的最大横摆角速度ωmax,i、过弯轨迹与实际轨迹的最大差别量Δmax,i,四个车轮出现的最大滑移率λmax,i;在本发明中所有数据量的量纲均采用国际标准量纲,若上述特殊点并不直接对应通信时间序列tk,则采用相邻状态点进行线性插值得到;并将上述特殊数据点的值组成一个向量,用于表征车辆通过弯道的行驶特征:
xT,i=[vs,i,Ls,i,Pm,i,ax,i,vI,i,vO,i,vmax,i,vmin,i,ay,imax,imax,imax,i]T    (1)
2)对已经驶过弯道的同一车辆类型的N辆车行驶状态的特征向量进行k-means聚类,聚成三类;
聚类主要的目标是使得式(2)的目标函数最小,式(2)中x为每个样本的特征向量,本发明中每个样本代表驶过弯道的一个车辆的状态数据,样本的特征向量x即为表征车辆通过弯道的行驶特征xT,i,Γj∈{Γ123}表示聚成的三类数据集,mj是每类数据集的均值即中心点:
J ( e ) = Σ j = 1 3 Σ x ∈ Γ j | | x - m j | | 2 - - - ( 2 )
式(2)中,J(e)表示聚成三类数据集的类内离散程度。k-means聚类是要使得式(2)的类内离散程度最小,反映每一类的聚集程度最好,从而将路侧通讯设备21记录的行驶状态分成的三类聚集程度较高的子数据集。其中,如图4所示的具体的k-means算法的过程为:
step1:将总体样本随机划分三类子数据集,总体样本指路侧存储设备22记录的同一车辆类型的历史行驶状态数据集;
step2:随机选择任何一个样本的特征向量,设x∈Γj,表示x属于子数据集Γj
step3:若该子数据集Γj中样本个数Nj=1,则转step2;
step4:根据式(2)计算样本的特征向量x被调整到其他各类子数据集时的J(e);
step5:将样本的特征向量x调整到一个子数据集使J(e)最小;
step6:若step2到step5连续迭代次数Count大于等于预先设定的Num次,样本的特征向量x所处的子数据集类别不变,则迭代终止;输出聚类的结果Γ123;否则转step2。
3)统计每一类数据集中的每个样本最大纵向加速度ax,i的均值ax,averagej)、最大横向加速度ay,i的均值ay,averagej)、最大横摆角速度ωmax,i的均值ωaveragej)以及最大过弯轨迹偏差绝对值Δmax,i的均值Δaveragej),计算公式如下:
a x , average ( Γ j ) = 1 N j Σ i ∈ Γ j a x , i a y , average ( Γ j ) = 1 N j Σ i ∈ Γ j a y , i ω average ( Γ j ) = 1 N j Σ i ∈ Γ j ω max , i Δ average ( Γ j ) = 1 N j Σ i ∈ Γ j Δ max , i j = 1,2,3 - - - ( 3 )
式(3)中,{Γj}j=1,2,3为无监督k-means聚类的结果,Nj为子数据集Γj的样本个数,i∈Γj表示属于类别Γj的样本。
4)计算聚类三个子数据集的驾驶员类型因子Index(Γj),并将驾驶员类型因子{Index(Γj)}j=1,2,3由大到小排序,依次将对应的无监督的聚类结果{Γj}j=1,2,3表示为激进、普通和保守三种驾驶员类型,三种驾驶员类型在后面统计学习过程中权重分别为λ123。例如Index(Γ2)>Index(Γ1)>Index(Γ3),则数据集Γ2的驾驶员类型为激进,学习权重取λ1;数据集Γ1的驾驶员类型为普通,学习权重取λ2;数据集Γ3的驾驶员类型为保守,学习权重取λ3。其中,驾驶员类型因子的计算过程为:
Index ( Γ j ) = | a x , average ( Γ j ) | max { | a x , average ( Γ j ) | } + | a y , average ( Γ j ) | max { | a y , average ( Γ j ) | } + | ω average ( Γ j ) | max { | ω average ( Γ j ) | } + | Δ average ( Γ j ) | max { | Δ average ( Γ j ) | } , j = 1,2,3
式中,max{|ax,averagej)|}表示取集合{|ay,average1)|,|ay,average2)|,|ay,average3)|}中最大的数,max{|ay,averagej)|},max{ωaveragej)|},max{|Δaveragej)|}含义类似,在此不再赘述。
统计学习最优速度序列及轨迹序列模块对于速度序列的统计学习方法如下:
1)将速度时间序列和轨迹时间序列离散为弯道的空间序列
如图5所示,黑色的弧形点画线表示实际的弯道中心轨迹,黑色弧形实线表示某一具体的车辆驶过弯道的行驶轨迹,s表示在弯道上沿弯道中心线到弯道起始标记点A的距离,将弯道在空间上等距离的离散化,表示为在本文中符号表示“定义为”。虚线表示垂直于sj处道路中心轨迹的切线。在弯道中心线上的sj处,虚线与车辆实际行驶轨迹的交点可对应车辆的行驶状态{xT,i(tj′),yT,i(tj′),vT,i(tj′)},tj′若不能直接对应某一通信的时间序列的tj,此时则采用相近的两个状态进行线性插值得到{xT,i(tj′),yT,i(tj′),vT,i(tj′)},将此时的车辆行驶状态表示为{xT,i(sj),yT,i(sj),vT,i(sj)}。从而,将不同车辆的速度-时间序列vT,i(tk)转化为速度-距离序列vT,i(s),且将行车轨迹序列{xT,i(tk),yT,i(tk)}转化为轨迹-距离序列{xT,i(s),yT,i(s)},其中 s = ^ { s 0 , s 1 , . . . , s j , . . . , s K } .
2)对速度-距离序列{vT,i(sk)}i=1,…,N采用带时间因子的加权最小二乘估计方法得到最优的voptimal(sk),其中 s k ∈ s = ^ { s 0 , s 1 , . . . , s K } ;
对于弯道内sk点处驶过同一车辆类型的N辆车的速度序列{vT,i(sk)}i=1,…,N,采用带时间因子的加权最小二乘估计方法得到最优的voptimal(sk),即:
v optimal ( s k ) = min x ^ Σ i = 1 N t i w i [ v T . i ( s k ) - x ^ ] 2 - - - ( 4 )
式(4)中,表示待估计的量,在该式中表示弯道内sk处待估计的速度值;加权系数wi的取值是根据无监督式学习分类模块对个体车辆中的驾驶员类型的聚类结果,相应地从权重λ123中取得;ti表示时间因子系数,对于不同时间段上的数据采用不同的权重:例如可将一天分为4个时间段,加权学习时对不同时间段历史数据采用不同的权重处理,与当前驶入弯道车辆所处的时间段相同的历史数据的权重大于其他时间段的权重,比如当前驶入弯道车辆所处时间段为6:00-11:59,则令时间因子t2>ti,i=1,3,4,如下表所示:
表2
式(4)的目的是对离散化弯道的每一点sk∈{s0,s1,…,sK}处,均计算得到一个加权最小二乘意义下的最优速度量voptimal(sk)。对离散化弯道的每一点sk∈{s0,s1,…,sK}重复式(4)计算过程,得到弯道过程中该类型车辆在统计意义上最优的速度序列Voptimax(sk),k=0,1,2,…,K本发明中的统计意义指的是带时间因子的加权最小二乘方法。
3)将轨迹-距离序列{xT,i(sk),yT,i(sk)}采用带时间因子加权最小二乘估计方法得到最优的轨迹序列:
对于轨迹-距离序列{xT,i(sk),yT,i(sk)},在弯道内sk点,分别对x坐标xT,i(s)和y坐标yT,i(s)采用带时间因子加权最小二乘估计方法:
x optimal ( s k ) = min x ^ Σ i = 1 N t i w i [ x T . i ( s k ) - x ^ ] 2 y optimal ( s k ) = min x ^ Σ i = 1 N t i w i [ y T . i ( s k ) - x ^ ] 2 - - - ( 5 )
式中ti与wi的含义与式(4)中一致,从而得到最优的轨迹序列{xoptimal(sk),yoptimal(sk)},k=O,1,2,…,K。
在一个优选的实施例中,对即将驶入弯道的车辆3的辅助过程根据车辆的类型(I代表普通轿车,II代表普通客车,III代表普通货车)利用路侧存储设备22记录的与即将驶入弯道的车辆3相同车辆类型的行车状态数据,根据上述的聚类、统计学习过程得到最优的速度序列和行车轨迹序列;本发明对即将驶入弯道的车辆3进行车速建议或者辅助控制具体是指:1)若即将驶入弯道的车辆3具有速度控制器、转向控制器,能够对油门开度、制动压力、方向盘转角进行主动控制时,车辆的车载通讯模块6从路侧通讯设备21接收最优速度序列和位置序列,车辆的控制器结合当前车辆在弯道中行驶的速度和位置进行闭环反馈控制,使得车辆通过弯道的速度和位置趋向于最优速度和位置序列,从而保证车辆以安全、畅通合理的方式通过该弯道;2)若即将驶入弯道车辆3不具备速度控制器、转向控制器,不具备主动干预油门开度、制动压力、方向盘转角的功能,则根据车辆当前的行驶速度与位置,通过适当的车载预警单元对驾驶员进行车速建议,促使车辆能够安全、畅通地通过弯道。不具备车载通信模块的车辆不在本发明方案的应用范围之内;同时对于具有车载通信模块的车辆,但是既不具有速度控制器、转向控制器,也不具有预警单元车辆也不在本发明的应用范围内。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (5)

1.一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:它包括若干驶过弯道的车辆、一信息通讯和处理基站和一即将驶入弯道的车辆;其中,每一所述驶过弯道的车辆上设置有一车载传感系统和一车载通讯模块,所述即将驶入弯道的车辆上也设置有一车载通讯模块,所述信息通讯和处理基站设置在弯道一侧,它包括一路侧通讯设备、一路侧存储设备和一路侧信息处理单元;所述车载传感系统采集所述驶过弯道的车辆的行驶状态信息并通过自车的车载通讯模块发送到所述路侧通讯设备,所述路侧通讯设备将接收的驶过弯道的车辆行驶状态信息传送给到所述路侧存储设备进行存储,所述路侧信息处理单元获取所述路侧存储设备中的数据并根据不同车辆类型将数据进行分类,并对同一类型的车辆行驶状态数据进行无监督式聚类和统计学习,得到最优速度序列和轨迹序列,所述路侧信息处理单元将得到的最优速度序列和轨迹序列发送到所述路侧存储设备进行存储,所述路侧通讯设备获取所述路侧存储设备中的最优速度序列和轨迹序列通过所述车载通讯模块传送给所述即将驶入弯道的车辆,所述即将驶入弯道的车辆根据车辆类型通过最优速度序列和轨迹序列对过弯道进行车速建议或辅助控制,使得车辆能够安全、畅通地通过弯道。
2.如权利要求1所述的一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:所述路侧通讯设备采用DSRC通信设备,所述DSRC通信设备采用专用短程通信技术。
3.如权利要求1或2所述的一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:所述路侧信息处理单元包括一数据预处理模块、一无监督式学习分类模块和一统计学习最优速度序列及轨迹序列模块;所述数据预处理模块用于对接收的所述驶过弯道车辆行驶状态数据进行Kalman滤波,滤除噪声;所述无监督式学习分类模块用于采用无监督式学习方法判断个体车辆中的驾驶员类型,以确定后续统计学习的权重;所述统计学习最优速度序列及轨迹序列模块用于对每一类型车辆采用统计学习的方法,学习得到统计意义上的最优速度序列和轨迹序列。
4.如权利要求3所述的一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:所述无监督式学习分类模块的工作过程为:
1)根据滤波后的各序列数据,对于同一类型的车辆行驶状态数据序列提取其特殊数据点,并将上述特殊数据点的值组成一个向量,用于表征车辆通过弯道的行驶特征:
xT,i=[vs,i,Ls,i,Pm,i,ax,i,vI,i,vO,i,vmax,i,vmin,i,ay,i,ωmax,imax,imax,i]T   (1)
式中,vs,i为开始制动时的车速、Ls,i为开始制动时相对弯道入口的距离、Pm,i为最大制动强度、ax,i为最大制动减速度、vI,i为进入弯道的速度、vO,i为离开弯道的速度、vmax,i为通过弯道过程中的最大速度、vmin,i为通过弯道过程中的最小速度、ay,i为通过弯道过程中的最大横向加速度、ωmax,i为通过弯道过程中的最大横摆角速度、Δmax,i为过弯轨迹与实际轨迹的最大差别量、λmax,i为四个车轮出现的最大滑移率;
2)对已经驶过弯道的同一车辆类型的N辆车行驶状态的特征向量进行k-means聚类,聚成三类;
聚类的目标是使得式(2)的目标函数最小,式(2)中x为每个样本的特征向量,每个样本代表驶过弯道的一个车辆的状态数据,样本的特征向量x即为表征车辆通过弯道的行驶特征xT,i,Γj∈{Γ123}表示聚成的三类数据集,mj是每类数据集的均值即中心点:
J ( e ) = Σ j = 1 3 Σ x ∈ Γ j | | x - m j | | 2 - - - ( 2 )
式(2)中,J(e)表示聚成三类数据集的类内离散程度;
3)统计每一类数据集中的每个样本最大纵向加速度ax,i的均值ax,averagej)、最大横向加速度ay,i的均值ay,averagej)、最大横摆角速度ωmax,i的均值ωaveragej)以及最大过弯轨迹偏差绝对值Δmax,i的均值Δaveragej),计算公式如下:
a x , average ( Γ j ) = 1 N j Σ i ∈ Γ j a x , i a y , average ( Γ j ) = 1 N j Σ i ∈ Γ j a y , i ω average ( Γ j ) = 1 N j Σ i ∈ Γ j ω max , i Δ average ( Γ j ) = 1 N j Σ i ∈ Γ j Δ max , i j = 1,2,3 - - - ( 3 )
式(3)中,{Γj}j=1,2,3为无监督k-means聚类的结果,Nj为子数据集Γj的样本个数,i∈Γj表示属于类别Γj的样本;
4)计算聚类三个子数据集的驾驶员类型因子Index(Γj),将驾驶员类型因子{Index(Γj)}j=1,2,3由大到小排序,依次将对应的无监督的聚类结果表示为激进、普通和保守三种驾驶员类型,其中,驾驶员类型因子Index(Γj)的计算过程为:
Index ( Γ j ) = | a x , average ( Γ j ) | max { | a x , average ( Γ j ) | } + | a y , average ( Γ j ) | max { | a y , average ( Γ j ) | } + | ω average ( Γ j ) | max { | ω average ( Γ j ) | } + | Δ average ( Γ j ) | max { | Δ average ( Γ j ) | } , j = 1,2,3 .
5.如权利要求3所述的一种基于车路/车车通讯的车辆弯道通行辅助系统,其特征在于:所述统计学习最优速度序列及轨迹序列模块对于速度序列的统计学习方法如下:
1)将速度时间序列和轨迹时间序列离散为弯道的空间序列,即将不同车辆的速度-时间序列转化为速度-距离序列,并将行车轨迹序列转化为轨迹-距离序列;
2)对速度-距离序列采用带时间因子的加权最小二乘估计方法得到最优的速度序列;
3)对轨迹-距离序列采用带时间因子的加权最小二乘估计方法得到最优的轨迹序列。
CN201410169875.4A 2014-04-25 2014-04-25 一种基于车路/车车通讯的车辆弯道通行辅助系统 Active CN103927895B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410169875.4A CN103927895B (zh) 2014-04-25 2014-04-25 一种基于车路/车车通讯的车辆弯道通行辅助系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410169875.4A CN103927895B (zh) 2014-04-25 2014-04-25 一种基于车路/车车通讯的车辆弯道通行辅助系统

Publications (2)

Publication Number Publication Date
CN103927895A true CN103927895A (zh) 2014-07-16
CN103927895B CN103927895B (zh) 2015-11-11

Family

ID=51146104

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410169875.4A Active CN103927895B (zh) 2014-04-25 2014-04-25 一种基于车路/车车通讯的车辆弯道通行辅助系统

Country Status (1)

Country Link
CN (1) CN103927895B (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104916152A (zh) * 2015-05-19 2015-09-16 苏州大学 基于车路协同的交叉口车辆右转引导系统及其引导方法
CN105844936A (zh) * 2016-03-25 2016-08-10 乐视控股(北京)有限公司 一种车辆转弯提示方法、设备及系统
CN106203703A (zh) * 2016-07-12 2016-12-07 武汉理工大学 一种用于汽车安全转向的动态数据的制作方法
CN107564280A (zh) * 2017-08-22 2018-01-09 王浩宇 基于环境感知的驾驶行为数据采集分析系统和方法
CN107830865A (zh) * 2017-10-16 2018-03-23 东软集团股份有限公司 一种车辆目标分类方法、装置、系统及计算机程序产品
CN108153854A (zh) * 2017-12-22 2018-06-12 东软集团股份有限公司 弯道分类方法、路侧单元、车载终端及电子设备
CN108280991A (zh) * 2017-01-05 2018-07-13 大唐高鸿信息通信研究院(义乌)有限公司 车载短距离通信网的车辆交通事故预测方法
CN108340923A (zh) * 2017-01-24 2018-07-31 曼卡车和巴士股份公司 用于使车辆、尤其商用车运行的方法
CN108648447A (zh) * 2018-05-08 2018-10-12 清华大学 基于行车安全场的车辆智能安全决策方法
CN110047287A (zh) * 2019-04-23 2019-07-23 临沂大学 一种道路状况提示方法及其装置
CN110136254A (zh) * 2019-06-13 2019-08-16 吉林大学 基于动态概率驾驶地图的驾驶辅助信息显示方法
CN110660232A (zh) * 2019-09-27 2020-01-07 深圳市金溢科技股份有限公司 可读存储介质、交通调度的实现方法、系统及其设备
CN111402630A (zh) * 2020-03-11 2020-07-10 浙江吉利汽车研究院有限公司 一种道路预警方法、装置及存储介质
CN112133120A (zh) * 2020-08-07 2020-12-25 北京万集科技股份有限公司 一种车辆提示处理方法及装置
CN112232254A (zh) * 2020-10-26 2021-01-15 清华大学 一种考虑行人激进度的行人风险评估方法
CN112562315A (zh) * 2020-11-02 2021-03-26 鹏城实验室 一种获取车流信息的方法、终端及存储介质
WO2021098222A1 (zh) * 2019-11-20 2021-05-27 华为技术有限公司 一种车辆行驶控制方法及相关设备
CN115713866A (zh) * 2022-10-11 2023-02-24 悉地(苏州)勘察设计顾问有限公司 一种基于车辆运行特征的道路静态信息主动服务方法
CN116013101A (zh) * 2023-03-27 2023-04-25 之江实验室 一种基于网联环境的无信号路口速度建议系统和方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002109684A (ja) * 2000-09-16 2002-04-12 Korea Telecommun 短距離高速通信装置を利用した車両航法サービス装置及び車両航法サービス方法
JP2002269678A (ja) * 2001-03-06 2002-09-20 Natl Inst For Land & Infrastructure Management Mlit 路側処理装置における検出データ入力チェック方式
CN201707794U (zh) * 2010-04-22 2011-01-12 深圳市卓尔芯电子技术有限公司 交通信息显示装置
CN103247185A (zh) * 2012-02-14 2013-08-14 厦门金龙联合汽车工业有限公司 一种车辆入弯防侧翻提醒系统和方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002109684A (ja) * 2000-09-16 2002-04-12 Korea Telecommun 短距離高速通信装置を利用した車両航法サービス装置及び車両航法サービス方法
JP2002269678A (ja) * 2001-03-06 2002-09-20 Natl Inst For Land & Infrastructure Management Mlit 路側処理装置における検出データ入力チェック方式
CN201707794U (zh) * 2010-04-22 2011-01-12 深圳市卓尔芯电子技术有限公司 交通信息显示装置
CN103247185A (zh) * 2012-02-14 2013-08-14 厦门金龙联合汽车工业有限公司 一种车辆入弯防侧翻提醒系统和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张德兆: "基于风险状态预估的弯道防侧滑超速预警系统", 《公路交通科技》, vol. 26, no. 1, 31 December 2009 (2009-12-31) *

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104916152B (zh) * 2015-05-19 2017-09-19 苏州大学 基于车路协同的交叉口车辆右转引导系统及其引导方法
CN104916152A (zh) * 2015-05-19 2015-09-16 苏州大学 基于车路协同的交叉口车辆右转引导系统及其引导方法
CN105844936A (zh) * 2016-03-25 2016-08-10 乐视控股(北京)有限公司 一种车辆转弯提示方法、设备及系统
CN106203703A (zh) * 2016-07-12 2016-12-07 武汉理工大学 一种用于汽车安全转向的动态数据的制作方法
CN108280991B (zh) * 2017-01-05 2020-12-01 大唐高鸿信息通信(义乌)有限公司 车载短距离通信网的车辆交通事故预测方法
CN108280991A (zh) * 2017-01-05 2018-07-13 大唐高鸿信息通信研究院(义乌)有限公司 车载短距离通信网的车辆交通事故预测方法
CN108340923A (zh) * 2017-01-24 2018-07-31 曼卡车和巴士股份公司 用于使车辆、尤其商用车运行的方法
CN108340923B (zh) * 2017-01-24 2023-09-05 曼卡车和巴士股份公司 用于使车辆、尤其商用车运行的方法
CN107564280A (zh) * 2017-08-22 2018-01-09 王浩宇 基于环境感知的驾驶行为数据采集分析系统和方法
CN107830865A (zh) * 2017-10-16 2018-03-23 东软集团股份有限公司 一种车辆目标分类方法、装置、系统及计算机程序产品
CN108153854A (zh) * 2017-12-22 2018-06-12 东软集团股份有限公司 弯道分类方法、路侧单元、车载终端及电子设备
CN108153854B (zh) * 2017-12-22 2021-08-17 东软集团股份有限公司 弯道分类方法、路侧单元、车载终端及电子设备
CN108648447A (zh) * 2018-05-08 2018-10-12 清华大学 基于行车安全场的车辆智能安全决策方法
CN108648447B (zh) * 2018-05-08 2019-04-05 清华大学 基于行车安全场的车辆智能安全决策方法
CN110047287A (zh) * 2019-04-23 2019-07-23 临沂大学 一种道路状况提示方法及其装置
CN110047287B (zh) * 2019-04-23 2021-04-02 临沂大学 一种道路状况提示方法及其装置
CN110136254A (zh) * 2019-06-13 2019-08-16 吉林大学 基于动态概率驾驶地图的驾驶辅助信息显示方法
CN110660232A (zh) * 2019-09-27 2020-01-07 深圳市金溢科技股份有限公司 可读存储介质、交通调度的实现方法、系统及其设备
WO2021098222A1 (zh) * 2019-11-20 2021-05-27 华为技术有限公司 一种车辆行驶控制方法及相关设备
CN111402630A (zh) * 2020-03-11 2020-07-10 浙江吉利汽车研究院有限公司 一种道路预警方法、装置及存储介质
CN112133120A (zh) * 2020-08-07 2020-12-25 北京万集科技股份有限公司 一种车辆提示处理方法及装置
CN112232254B (zh) * 2020-10-26 2021-04-30 清华大学 一种考虑行人激进度的行人风险评估方法
CN112232254A (zh) * 2020-10-26 2021-01-15 清华大学 一种考虑行人激进度的行人风险评估方法
CN112562315A (zh) * 2020-11-02 2021-03-26 鹏城实验室 一种获取车流信息的方法、终端及存储介质
CN115713866A (zh) * 2022-10-11 2023-02-24 悉地(苏州)勘察设计顾问有限公司 一种基于车辆运行特征的道路静态信息主动服务方法
CN115713866B (zh) * 2022-10-11 2023-08-22 悉地(苏州)勘察设计顾问有限公司 一种基于车辆运行特征的道路静态信息主动服务方法
CN116013101A (zh) * 2023-03-27 2023-04-25 之江实验室 一种基于网联环境的无信号路口速度建议系统和方法
CN116013101B (zh) * 2023-03-27 2023-06-20 之江实验室 一种基于网联环境的无信号路口速度建议系统和方法

Also Published As

Publication number Publication date
CN103927895B (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
CN103927895B (zh) 一种基于车路/车车通讯的车辆弯道通行辅助系统
US11648946B2 (en) Intelligent vehicle platoon lane change performance evaluation method
CN107993453B (zh) 一种基于车路协同的弯道安全车速计算方法
CN102295004B (zh) 一种车道偏离预警方法
CN107945574B (zh) 一种车辆碰撞预警方法、装置及设备
CN103921719B (zh) 驾驶员交互式商用车辆侧翻预警方法和系统
CN107161143A (zh) 一种采用人工势场法的车辆主动避撞方法
CN106056972A (zh) 一种基于车辆行驶速度和位置信息融合的安全防撞预警方法
CN102795225B (zh) 利用驾驶员侧纵向控制模型检测驾驶员干扰状态的方法
KR101575298B1 (ko) 차량 제어 시스템 및 방법
CN105151043A (zh) 一种无人驾驶汽车紧急避让的系统和方法
CN109544992A (zh) 一种基于GPS和ZigBee的车辆碰撞预警系统
CN105667509A (zh) 一种用于汽车自适应巡航控制系统的弯道控制系统及方法
CN107633705B (zh) 一种基于路面反馈的车道偏离智能预警系统及方法
CN104210489A (zh) 车路协同环境下车辆与行人碰撞规避方法与系统
CN103395419A (zh) 基于安全间距策略的车辆队列行驶控制系统及其控制方法
CN110544377A (zh) 基于车路协同的交叉口行人避撞方法
CN102837703A (zh) 估算车辆速度曲线的方法和系统及包括该系统的机动车辆
CN106205148A (zh) 一种危险货物罐车弯道安全车速获取方法及超速警示系统
CN104192144A (zh) 一种汽车主动防撞弯道虚警消除方法
CN105551282A (zh) 一种超车提示方法和装置
CN105857294A (zh) 一种汽车换道避撞控制方法
CN103164962A (zh) 一种山区公路急弯路段实时车速预警方法
CN110085057A (zh) 一种基于车路协同的车辆弯道防侧滑侧翻安全控制系统及控制方法
CN113192331A (zh) 一种网联环境下面向骑行安全的智能预警系统及预警方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant