CN112541632A - 一种基于多属性决策的驾驶行为安全性评估方法 - Google Patents

一种基于多属性决策的驾驶行为安全性评估方法 Download PDF

Info

Publication number
CN112541632A
CN112541632A CN202011472843.3A CN202011472843A CN112541632A CN 112541632 A CN112541632 A CN 112541632A CN 202011472843 A CN202011472843 A CN 202011472843A CN 112541632 A CN112541632 A CN 112541632A
Authority
CN
China
Prior art keywords
index
driver
evaluation
lane
matrix
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
CN202011472843.3A
Other languages
English (en)
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.)
Jiangsu University
Original Assignee
Jiangsu 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 Jiangsu University filed Critical Jiangsu University
Priority to CN202011472843.3A priority Critical patent/CN112541632A/zh
Publication of CN112541632A publication Critical patent/CN112541632A/zh
Pending legal-status Critical Current

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
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • Development Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Educational Administration (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于多属性决策的驾驶行为安全性评估方法,通过实车数据采集平台采用不同驾驶人的驾驶数据,基于驾驶人的驾驶数据,初步选取并计算多个评价指标;对评价指标进行Pearson相关性分析,保留最能体现评估结果的关键指标;用层次分析法确定各个评估指标的主观权重pj;用熵值法确定各个评估指标的客观权重qj;基于得到的主观权重pj和客观权重qj通过乘法集成法得到指标综合权重wj,利用改进型TOPSIS方法进行多属性决策,计算驾驶人与理想解的相对贴近度Ci,根据相对贴近度Ci评估对象的驾驶行为安全等级。

Description

一种基于多属性决策的驾驶行为安全性评估方法
技术领域
本发明属于车辆主动安全领域,尤其是一种基于多属性决策的驾驶行为安全性评估方法。
背景技术
世卫组织在《道路安全全球现状报告2018》指出,世界范围内每年超过120万人丧生于道路交通事故,超过2000万人因交通事故致残致伤。仅2018年,我国交通事故发生总数24.49万起,造成25.85万人受伤,其中机动车交通事故数21.62万起,造成5.81万人死亡,22.74万人受伤,机动车事故受伤人数占比达到88%,所造成的直接经济损失达13.8456亿元。
驾驶人作为道路交通系统的主体,是道路交通安全的核心问题,对于驾驶人行为的研究已成为各国研究人员的研究热点。面对相同的道路交通环境,不同的驾驶人踩制动踏板和油门踏板,转方向盘的方式和跟车距离也各不相同,所以不同驾驶行为导致的事故风险性也各不相同。
研究驾驶行为与事故倾向之间的关系,可以量化评估个人驾驶行为的安全性并将其作为事故风险的评价指标,有效预测驾驶风险,有助于区分安全的驾驶人和不安全的驾驶人。由于驾驶行为实时数据采集方式的多样性和当前道路交通智能化的发展,基于大数据和便携设备的实时驾驶行为安全性评价研究将会是一个趋势,有利于提高车辆主动安全性。
发明内容
为了解决现有技术中的不足,本发明提出了一种基于多属性决策的驾驶行为安全性评估方法,本发明的评估方法能够对不同驾驶行为的安全等级进行较为客观的评估,有效提高该评估方法的可行性和准确性。
本发明所采用的技术方案如下:
S1、通过实车数据采集平台采集不同驾驶人的驾驶数据,基于驾驶人的驾驶数据,初步选取并计算多个评价指标;
S2、对评价指标进行Pearson相关性分析,保留最能体现评估结果的关键指标;
S3、对各个关键指标进行赋权,具体为:
S3.1、用层次分析法确定各个评估指标的主观权重pj
S3.2、用熵值法确定各个评估指标的客观权重qj
S3.3、基于得到的主观权重pj和客观权重qj通过乘法集成法得到指标综合权重wj
Figure BDA0002836446160000021
由综合权重wj构成综合权重矩阵w,表示为:
Figure BDA0002836446160000022
S4、利用改进型TOPSIS方法进行多属性决策,计算驾驶人与理想解的相对贴近度Ci,根据相对贴近度Ci评估对象的驾驶行为安全等级。
进一步,S4中利用改进型TOPSIS方法进行多属性决策的方法为:
S4.1、根据原始数据建立决策矩阵:
Figure BDA0002836446160000023
利用极值法对指标原始数据进行规范化处理,得到规范化矩阵σ:
Figure BDA0002836446160000024
对各指标赋权,将规范化指标数据矩阵乘以综合权重矩阵就得到评估矩阵V。
Figure BDA0002836446160000025
确定正理想解
Figure BDA0002836446160000026
和负理想解
Figure BDA0002836446160000027
为第n个正理想解,
Figure BDA0002836446160000028
为第n个负理想解;若评估指标Cj为效益型指标,则
Figure BDA0002836446160000029
若评估指标Cj为成本型指标,则
Figure BDA00028364461600000210
Vij为评估矩阵V中第i行第j列的参数;
S4.2、设第i个驾驶人的所有指标向量为:Xi=(Xi1,Xi2,…,Xin),则第i个驾驶人到正负理想解的马氏距离分别为:
Figure BDA0002836446160000031
Figure BDA0002836446160000032
其中,Ai为第i个评估对象即第i个驾驶人,∑表示相关系数矩阵;
S4.3、驾驶人Ai到正理想解和负理想解的联系向量距离分别表示如下:
Figure BDA0002836446160000033
Figure BDA0002836446160000034
其中,
Figure BDA0002836446160000035
为正同一度系数;
Figure BDA0002836446160000036
为正差异度系数;
Figure BDA0002836446160000037
为正对立度系数;
S4.4、通过马氏距离求得的d(Ai,S+)、d(Ai,S-)和联系向量距离求得的
Figure BDA0002836446160000038
进行无量纲化处理后,定义各待评估对象与正负理想解的合成贴近度为:
Figure BDA0002836446160000039
Figure BDA00028364461600000310
式中,系数α、β满足α+β=1;则各驾驶人与理想解的相对贴近度为:
Figure BDA00028364461600000311
据相对贴近度Ci的大小对驾驶人进行驾驶行为的安全等级排序,Ci值越趋近于1,说明该评估对象的驾驶行为安全等级越高。
进一步,S3.1中确定各个评估指标的主观权重pj的过程为:
S3.1.1、基于关键指标建立递阶层次;
S3.1.2、构造判断矩阵,若上一层指标支配的下一层指标有n个,n个指标就构成一个两两比较的判断矩阵A=(aij)n×n,aij表示指标i与指标j的重要性之比;
S3.1.3、若向量ω=[ω12,…,ωn]T满足:Aω=λmaxω,λmax为判断矩阵A的最大特征根,ω为对应的特征向量,归一化的特征向量可以作为主观权重pj
S3.1.4、在计算得到主观权重pj之后,进行一致性检验,一致性指标:
Figure BDA00028364461600000312
其中,n为矩阵阶数;随机一致性比例:
Figure BDA0002836446160000041
RI为平均随机一致性指标,当随机一致性比例CR<0.1时,认为判断矩阵具有一致性,若一致性检验通过,说明该主观权重pj符合要求,否则需要调整判断矩阵中的aij重新计算主观权重pj
进一步,S3.2、采用熵值法确定各个评估指标的客观权重qj的过程为:
S3.2.1、计算第j项指标下,第i个驾驶人的特征比重:
Figure BDA0002836446160000042
其中,i=1,2,…,m,m为驾驶人总数,j=1,2,…,n,n为评估指标总数;
S3.2.2、计算第j项指标的熵值:
Figure BDA0002836446160000043
Figure BDA0002836446160000044
k为系数,m为驾驶人总数;
S3.2.3、计算第j项指标的差异性系数:gj=1-ej
S3.2.4、确定第j项指标的客观权重:
Figure BDA0002836446160000045
进一步,S2中相关性分析的过程如下:根据相关系数公式计算不同指标之间的相关性:
Figure BDA0002836446160000046
Sxy为协方差,表示为
Figure BDA0002836446160000047
Sx为某一种指标的标准差,表示为
Figure BDA0002836446160000048
Sy为另一种指标的标准差,表示为
Figure BDA0002836446160000049
Xi表示某一个指标的实际值,
Figure BDA00028364461600000410
表示该指标的平均值,Yi表示另一个指标的实际值,
Figure BDA00028364461600000411
表示另一个指标的平均值,n为采样的数据样本数;当相关系数rxy为1时,称为完全正相关;当相关系数为-1时,称为完全负相关;相关系数的绝对值越大,相关性越强;相关系数越接近于0,相关性越弱。
进一步,用于安全性评估的关键指标分别是:变换车道视觉感知行为安全性指标LCTR、前方障碍物视觉感知安全性指标OSIA、制动幅值变化率BACV、方向盘转角熵值SE、纵向速度标准差Vstd、纵向加速度变化率LACR、方向盘转角标准差θwstd、横向加速度标准差Lastd和横向车道保持安全性指标LKSR。
进一步,获得变换车道视觉感知行为安全性指标LCTR的方法为:将驾驶人换道过程的注视区域划分为5块:车道正前方、车道左前方、车道右前方、左后视镜和右后视镜;通过眼动仪采集的不同驾驶人注视点数据,且利用K-means方法对通过眼动仪采集的不同驾驶人注视点数据进行聚类;最终选择左侧换道驾驶人注视左后视镜和车道正前方的平均次数,右侧换道注视右后视镜和车道正前方的平均次数作为分析驾驶人变换车道视觉感知行为安全性的指标;每一位驾驶人换道过程观测不同位置的次数与样本中规范性要求下变换车道注视行为水平均值进行对比,得到驾驶人变换车道视觉感知行为安全性指标
Figure BDA0002836446160000051
F为各驾驶人换道过程观测不同区域的总次数;A为规范性要求下变换车道注视行为水平均值。
本发明的有益效果:
1、本发明基于实车试验方式采集驾驶人视觉行为参数、驾驶人操纵行为数据和车辆运行状态数据,完善了驾驶行为安全性评估指标体系。其次,通过Pearson相关性分析确定关键指标,避免了因为选取的评价指标之间关联性而造成指标的冗余以及增大计算量等问题,可以提高本评估方法的效率、降低计算复杂度。
2、本发明对关键指标进行赋权时,然后通过主观的层次分析法和客观的熵值法结合确定各个评估指标的综合权重矩阵;可以使得到的指标权重更合理。
3、由于驾驶行为安全性表征指标之间存在一定的相关性,本发明通过一种改进的TOPSIS多属性决策方法对不同驾驶人行为的安全等级进行评估。本发明通过马氏距离和联系向量距离结合代替传统计算方法中的欧氏距离,可以考虑评估指标存在相关性时带来的影响,提高评估结果的准确性,同时对评估对象的相对优劣进行客观合理的评估。
附图说明
图1为搭建的实车数据采集平台;
图2为驾驶人换道过程注视区域划分示意图;
图3为世界坐标系下驾驶人前视预瞄距离示意图;
图4为层次分析法确定指标主观权重流程图;
图5为本发明的流程框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用于解释本发明,并不用于限定本发明。
一种基于多属性决策的驾驶行为安全性评估方法,包括如下步骤:
S1、初步选取多个评价指标;
在道路交通条件相同的情况下,由于不同驾驶人在跟车距离和注视区域等方面都有很大差异,最终影响车辆行驶稳定性和安全性。本发明利用K-means方法对通过眼动仪采集的不同驾驶人注视点数据进行聚类,如图2所示,将驾驶人换道过程的注视区域划分为5块:车道正前方,车道左前方,车道右前方,左后视镜和右后视镜。
通过对驾驶人注视点位置的聚类分析,左侧换道过程与右侧换道过程的视点分布均呈现出一定的规律。因此,本发明对驾驶人在要求规范性变换车道情况下,注视不同位置的平均次数进行统计分析,发现存在明显的差异性。本发明选择左侧换道驾驶人注视左后视镜和车道正前方的平均次数,右侧换道注视右后视镜和车道正前方的平均次数作为分析驾驶人变换车道视觉感知行为安全性的指标。
经过分析,本发明提出利用每一位驾驶人换道过程观测不同位置的次数与样本中规范性要求下变换车道注视行为水平均值进行对比,得到驾驶人变换车道视觉感知行为安全性指标(LCTR)。
Figure BDA0002836446160000061
式中,F为各驾驶人换道过程观测不同区域的总次数;A为规范性要求下变换车道注视行为水平均值。若LCTR=1,说明该驾驶人换道视觉行为特性与规范性要求一致,而LCTR>1或LCTR<1都说明该驾驶人换道视觉感知行为安全性降低。
驾驶人前视预瞄行为可以指导汽车未来一段时间的运动轨迹,也就是说预瞄点可以落在汽车未来行驶轨迹上。如图3所示,可利用驾驶人视线点在世界坐标系下的投影,通过三角函数关系计算得到驾驶人前视预瞄距离:
Figure BDA0002836446160000062
式中,GazeOrigin.y为视点位置的y方向坐标;h为地面与世界坐标系水平面的垂直距离;GazeDirection.y为世界坐标系下视线在y方向的单位分量;GazeDirection.z为世界坐标系下视线在z方向的单位分量。
在车辆正常行驶过程中,如果前方有障碍物出现,而道路条件又无法进行转向操作时,紧急制动就十分必要,能否在一定距离内停车直接关系到驾驶人的生命安全。本发明采用的安全制动距离表示为:
Figure BDA0002836446160000071
式中,S为安全制动距离;v为自车车速;vrel为相对速度;t1为驾驶人的反应时间(取0.1s);t2为制动器延迟时间(取0.6s);a1为自车最大减速度(取6m/s2);a2为前方障碍物最大减速度(取8m/s2);d0为最小停车距离(取3m)。
通过安全制动距离与驾驶人前视预瞄距离的比值(OSIA)来表征驾驶人直行过程对前方障碍物视觉感知安全指标:
Figure BDA0002836446160000072
若OSIA<1,说明前视预瞄距离大于安全制动距离,进行紧急制动所需距离足够,安全性较高;若OSIA>1,则说明前视预瞄距离小于安全制动距离,进行紧急制动已没有足够距离,安全性较低。
纵向操纵行为安全性指标主要表现为制动幅值变化率(BACV),在采样周期T内的制动幅值变化率可以反映出一段时间周期内制动行为的急缓程度。计算式如下:
Figure BDA0002836446160000073
式中,N为采样周期内的采样次数,b(i)为第i次制动幅值采样样本。
横向操纵行为安全性指标主要表现为方向盘转角熵值(SE),方向盘转角熵值越大,说明驾驶人的横向操纵行为越不稳定,驾驶人的心理负荷增大,进一步影响行车安全性。SE是根据方向盘转角的预测偏差出现的概率来计算的,θ(n)是t(n)时刻方向盘转角的实际值,θp(n)是根据泰勒公式展开计算得到的方向盘转角预测值。
具体计算式如下:
Figure BDA0002836446160000074
方向盘转角的预测偏差为:e(n)=θ(n)p(n) (7)
为方便进行计算和分析,将方向盘转角预测误差值划分为9个区域(P1到P9)。其中有90%的预测误差值落在[-α,+α]范围内,α确定为1°,根据误差值落在各个区间的频数确定各区间的概率分布值Pi,方向盘转角熵值计算式如下:
Figure BDA0002836446160000081
SE反映出驾驶人在任意时刻转动方向盘角度预测误差值分布情况的不确定性,熵值越小说明驾驶人横向操纵行为安全性越高,熵值越大则说明安全性越低。
纵向速度标准差体现出车速在所采集时间段内的变化情况,计算式如下:
Figure BDA0002836446160000082
式中,N为采样周期内的采样次数,Vstd为纵向速度标准差;Vi为纵向速度实际值;
Figure BDA0002836446160000083
为纵向速度平均值。
纵向加速度标准差反映出加速度在所采集时间段内的变化情况,也是一个与车辆纵向运行状态安全性有关的重要指标,计算式如下:
Figure BDA0002836446160000084
式中,N为采样周期内的采样次数,Vastd为纵向加速度标准差;αi为纵向加速度实际值;
Figure BDA0002836446160000085
为纵向加速度平均值。
急加速行为也是一种影响驾驶安全的因素,过多猛踩加速踏板行为更容易造成追尾,而采样时间段内纵向加速度的变化率可以反映出一段时间内加速的急缓程度,纵向加速度变化率的值越大,说明纵向驾驶安全性越低。计算式如下:
Figure BDA0002836446160000086
式中,N为采样周期内的采样次数,a(i)为第i次纵向加速度采样样本;采样时间段长度为T。
车辆的横向运动是由于驾驶人对方向盘的操作产生的,因此要表征车辆横向运动状态,方向盘转角是一个重要指标。
方向盘转角标准差计算式如下:
Figure BDA0002836446160000091
式中,N为采样周期内的采样次数,θwstd为方向盘转角标准差;dgi为方向盘转角实际值;
Figure BDA0002836446160000092
为方向盘转角在所测时间周期内的平均值。
横向加速度标准差计算式如下:
Figure BDA0002836446160000093
式中,N为采样周期内的采样次数,Lastd为横向加速度标准差;βi为横向加速度实际值;
Figure BDA0002836446160000094
为所测时间段内横向加速度平均值。
在驾驶过程中,如果车辆的横向偏移量过大,会造成压线行驶,形成安全隐患,严重影响行车安全。而车辆的横向偏移距离均值与实际车道宽度的百分比反映出车辆横向偏移相对于所在车道宽度的偏差程度。本发明通过横向偏移距离均值占所在车道宽度百分比指标(LKSR)来表征驾驶人横向车道保持安全性。
Figure BDA0002836446160000095
设dl(i)≤0,dr(i)≥0
式中,N为采样周期内的采样次数,dl(i)为车辆中心轴到左侧车道线的距离;dr(i)为车辆中心轴到右侧车道线的距离。
为了得到上述所选取的评价指标,本发明通过搭建的实车数据采集平台采集不同驾驶人的各类型数据;具体地,如图1所示,试验车辆集成了CAN总线车辆运行状态数据采集和驾驶人眼动信息采集系统,对相同道路交通环境下,不同驾驶人行为的安全性进行分析。所选择的试验路段为某城市快速公路,试验车辆装备有毫米波雷达、摄像头、行车记录仪、Smart-eyes眼动仪、方向盘转角传感器、轮速传感器和GPS等。且通过工控机记录上述试验车辆的各类数据,并将采集的数据集中转换,实现数据的同步输出。
S2、并通过Pearson相关性分析,保留最能体现评估结果的关键指标。
在一个评价集中,评价指标之间会存在某种联系,利用Pearson相关性分析可得到指标之间是否存在过强的关联性,如果指标间的关联性过强,会造成指标的冗余,计算量增加。因此本发明通过Pearson相关性分析,保留最能体现评估结果的关键指标,来减少计算的复杂度。
因此,对S1中初步选择的评价指标进行相关性分析,相关性分析的过程如下:
根据相关系数公式计算不同指标之间的相关性:
Figure BDA0002836446160000101
其中,Sxy为协方差,表示为
Figure BDA0002836446160000102
Sx为某一种指标的标准差,表示为
Figure BDA0002836446160000103
Sy为另一种指标的标准差,表示为
Figure BDA0002836446160000104
Xi表示某一个指标的实际值,
Figure BDA0002836446160000105
表示该指标的平均值,Yi表示另一个指标的实际值,
Figure BDA0002836446160000106
表示另一个指标的平均值,n为采样的数据样本数。当相关系数rxy为1时,称为完全正相关;当相关系数为-1时,称为完全负相关;相关系数的绝对值越大,相关性越强;相关系数越接近于0,相关性越弱。
通过相关性分析可以发现,纵向速度标准差和纵向加速度标准差的相关性较高,又属于同一类指标,因此保留纵向速度标准差指标,去除纵向加速度标准差指标,最后得到9个用于安全性评估的关键指标如下:变换车道视觉感知行为安全性指标(LCTR)、前方障碍物视觉感知安全性指标(OSIA)、制动幅值变化率(BACV)、方向盘转角熵值(SE)、纵向速度标准差Vstd、纵向加速度变化率(LACR)、方向盘转角标准差θwstd、横向加速度标准差Lastd和横向车道保持安全性指标(LKSR)。
S3、得到9个关键指标之后,需要对各个关键指标进行赋权;用层次分析法确定各个评估指标的主观权重pj;采用熵值法确定各个评估指标的客观权重qj;基于得到的主观权重pj和客观权重qj通过乘法集成法得到指标综合权重。
S3.1、如图4用层次分析法确定各个评估指标的主观权重pj的过程为:
S3.1.1、建立递阶层次;一般采用9级标度来表示不同指标之间的相对重要程度:1表示两个指标之间重要程度相同,9表示前者比后者重要很多。如下表所示:
Figure BDA0002836446160000107
Figure BDA0002836446160000111
S3.1.2、构造判断矩阵,若上一层指标支配的下一层指标有n个,那这n个指标就构成一个两两比较的判断矩阵A=(aij)n×n,其中,aij表示指标i与指标j的重要性之比。
S3.1.3、若向量ω=[ω12,…,ωn]T满足
Aω=λmaxω (16)
式中,λmax为判断矩阵A的最大特征根,ω为对应的特征向量,归一化的特征向量可以作为主观权重pj
S3.1.4、在计算得到主观权重pj之后,必须进行一致性检验。一致性指标计算式如下:
Figure BDA0002836446160000112
其中,n为矩阵阶数;
随机一致性比例计算式如下:
Figure BDA0002836446160000113
其中,RI为平均随机一致性指标,对于1~9阶判断矩阵具体取值如下表所示:
阶数 1 2 3 4 5 6 7 8 9
RI 0 0 0.58 0.9 1.12 1.24 1.32 1.41 1.45
当随机一致性比例CR<0.1时,认为判断矩阵具有一致性。若一致性检验通过,说明该主观权重pj符合要求,否则需要调整判断矩阵中的aij重新计算主观权重pj
S3.2、熵值法是一种客观赋权法,根据各项指标原始数据所提供的信息量的大小来确定指标权重。设Xij为第i个驾驶人的第j项指标数据,对于给定的指标,Xij的差异性越大,说明该指标对系统的比较作用就越大,所包含的信息就越多。本发明采用熵值法确定各个评估指标的客观权重qj的过程为:
S3.2.1、计算第j项指标下,第i个驾驶人的特征比重:
Figure BDA0002836446160000121
其中,i=1,2,…,m,m为驾驶人总数,j=1,2,…,n,n为评估指标总数;
S3.2.2、计算第j项指标的熵值:
Figure BDA0002836446160000122
其中,k为系数,m为驾驶人总数;
S3.2.3、计算第j项指标的差异性系数:
gj=1-ej (21)
S3.2.4、确定第j项指标的客观权重:
Figure BDA0002836446160000123
其中,m为驾驶人总数;
S3.3、基于主观权重pj和客观权重qj计算综合权重wj的过程为:
Figure BDA0002836446160000124
式中,wj为综合权重,
Figure BDA0002836446160000125
由综合权重构成综合权重矩阵w表示为:
Figure BDA0002836446160000126
pj和qj分别为通过主观赋权方法和客观赋权方法得到的指标权重。
S4、TOPSIS的原理是通过计算各评估对象与理想化目标的相对贴近度来进行相对优劣排序,其中理想化目标包含了正理想解S+和负理想解S-,正理想解是待评估方案集中的最优解,它表示各个单一指标值都达到了方案集中的最优,而负理想解则是所有待评估方案中的最劣解。传统的TOPSIS多属性决策方法所使用的距离计算公式是欧氏距离,不能考虑评估指标存在相关性时带来的影响,而驾驶行为安全性表征指标之间存在一定的相关性,因此会造成评估结果不准确。同时,还会出现某评估对象与正负理想解的距离都近的问题,所以也不能对评估对象的相对优劣进行客观合理的评估。因此本发明利用改进型TOPSIS多属性决策方法进行安全性评估,计算驾驶人与理想解的相对贴近度Ci,根据相对贴近度Ci的大小对不同驾驶人行为的安全等级进行排序。具体过程为:
S4.1、设有m个驾驶人分别是A1,A2,…,Am,n个评估指标分别是C1,C2,…,Cn,Xij是第i个驾驶人的第j项指标值。
(1)根据原始数据建立决策矩阵:
Figure BDA0002836446160000131
(2)由于不同指标数据的量纲存在差异,本发明先利用极值法对指标原始数据进行规范化处理,得到规范化矩阵σ:
Figure BDA0002836446160000132
(3)对各指标赋权,将规范化指标数据矩阵乘以综合权重矩阵就得到评估矩阵V。
Figure BDA0002836446160000133
(4)确定正理想解S+和负理想解S-
Figure BDA0002836446160000134
若指标Cj为效益型指标
Figure BDA0002836446160000135
若指标Cj为成本型指标
Figure BDA0002836446160000136
S4.2、分别计算各驾驶人到正负理想解的距离;
设第i个驾驶人的所有指标向量为:Xi=(Xi1,Xi2,…,Xin) (30)
Figure BDA0002836446160000141
分别为正负理想解所对应的空间坐标,
Figure BDA0002836446160000142
为第n个正理想解,
Figure BDA0002836446160000143
为第n个负理想解;则第i个驾驶人到正负理想解的马氏距离分别为:
Figure BDA0002836446160000144
Figure BDA0002836446160000145
其中,Ai为第i个评估对象即第i个驾驶人,∑表示相关系数矩阵,S+为正理想解和S-为负理想解。
S4.3、定义1.设集合A0和B1的联系度为(A0~B1)=a+bi+cj,其中,a为同一度,b为差异度,c为对立度,则向量(a,b,c)为集合A0和B1的联系向量。向量记为μ=(a,b,c),i、j为空间单位向量;
定义2.设集合A0和Bs的联系向量为μs=(as,bs,cs),集合A0和Bk的联系向量为μk=(ak,bk,ck),则集合Bs和Bk的联系向量距离为:
Figure BDA0002836446160000146
根据集对分析思想,把正负理想解看成是确定不确定系统中相互对立的关系。所以需要计算各个驾驶人与正理想解的联系度_k+、负理想解的联系度_k-
若评估指标为效益型指标,由评估对象Ai和正理想解S+组成集对H+=(Ai,S+),则联系度计算方式如下:
Figure BDA0002836446160000147
其中,
Figure BDA0002836446160000148
为正同一度系数;
Figure BDA0002836446160000149
为正差异度系数;
Figure BDA00028364461600001410
为正对立度系数;
Figure BDA00028364461600001411
为第n个驾驶人到正理想解的联系度;
Figure BDA00028364461600001412
为第t个驾驶人到正理想解的联系度,表示为
Figure BDA00028364461600001413
为第t个驾驶人的正同一度系数,
Figure BDA00028364461600001414
为第t个驾驶人的正差异度系数,
Figure BDA00028364461600001415
为第t个驾驶人的正对立度系数。
Figure BDA00028364461600001416
时,
Figure BDA00028364461600001417
Figure BDA00028364461600001418
时,
Figure BDA00028364461600001419
同理可计算评估对象Ai与负理想解S-的联系度,由评估对象Ai与负理想解组成集对H-=(Ai,S-),则有:
Figure BDA0002836446160000151
其中,
Figure BDA0002836446160000152
为负同一度系数;
Figure BDA0002836446160000153
为负差异度系数;
Figure BDA0002836446160000154
为负对立度系数;
Figure BDA0002836446160000155
为第n个驾驶人到负理想解的联系度;
Figure BDA0002836446160000156
为第t个驾驶人到负理想解的联系度,表示为
Figure BDA0002836446160000157
为第t个驾驶人的负同一度系数,
Figure BDA0002836446160000158
为第t个驾驶人的负差异度系数,
Figure BDA0002836446160000159
为第t个驾驶人的负对立度系数。
Figure BDA00028364461600001510
时,
Figure BDA00028364461600001511
Figure BDA00028364461600001512
时,
Figure BDA00028364461600001513
若评估指标为成本型指标,则联系度计算方式与效益型指标恰好相反。
根据前面的定义,驾驶人Ai到正负理想解的联系向量距离计算式如下:
Figure BDA00028364461600001514
Figure BDA00028364461600001515
S4.4、通过以上对马氏距离和联系向量距离改进原理的介绍,对通过马氏距离求得的d(Ai,S+)、d(Ai,S-)和联系向量距离求得的
Figure BDA00028364461600001516
进行无量纲化处理后,定义各待评估对象与正负理想解的合成贴近度为:
Figure BDA00028364461600001517
Figure BDA00028364461600001518
式中,系数α+β=1,具体数值由决策者决定,本发明取系数α=β=0.5
则各驾驶人与理想解的相对贴近度为:
Figure BDA00028364461600001519
这样计算得到的新相对贴近度克服了传统TOPSIS多属性决策方法在距离计算方式上的两个缺陷,使得评估结果更加科学合理。
最后,根据相对贴近度Ci的大小对驾驶人进行驾驶行为的安全等级排序,Ci值越趋近于1,说明该评估对象的驾驶行为安全等级越高。
以上实施例仅用于说明本发明的设计思想和特点,其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施,本发明的保护范围不限于上述实施例。所以,凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰,均在本发明的保护范围之内。

Claims (7)

1.一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,包括如下步骤:
S1、通过实车数据采集平台采集不同驾驶人的驾驶数据,基于驾驶人的驾驶数据,初步选取并计算多个评价指标;
S2、对评价指标进行Pearson相关性分析,保留最能体现评估结果的关键指标;
S3、对各个关键指标进行赋权,具体为:
S3.1、用层次分析法确定各个评估指标的主观权重pj
S3.2、用熵值法确定各个评估指标的客观权重qj
S3.3、基于得到的主观权重pj和客观权重qj通过乘法集成法得到指标综合权重wj
Figure FDA0002836446150000011
由综合权重wj构成综合权重矩阵w,表示为:
Figure FDA0002836446150000012
S4、利用改进型TOPSIS方法进行多属性决策,计算驾驶人与理想解的相对贴近度Ci,根据相对贴近度Ci评估不同驾驶人的驾驶行为安全等级。
2.根据权利要求1所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,S4中利用改进型TOPSIS方法进行多属性决策的方法为:
S4.1、根据原始数据建立决策矩阵:
Figure FDA0002836446150000013
利用极值法对指标原始数据进行规范化处理,得到规范化矩阵σ:
Figure FDA0002836446150000014
对各指标赋权,将规范化指标数据矩阵乘以综合权重矩阵就得到评估矩阵V。
Figure FDA0002836446150000021
确定正理想解
Figure FDA0002836446150000022
和负理想解
Figure FDA0002836446150000023
Figure FDA0002836446150000024
为第n个正理想解,
Figure FDA0002836446150000025
为第n个负理想解;若评估指标Cj为效益型指标,则
Figure FDA0002836446150000026
若评估指标Cj为成本型指标,则
Figure FDA0002836446150000027
Vij为评估矩阵V中第i行第j列的参数;
S4.2、设第i个驾驶人的所有指标向量为:Xi=(Xi1,Xi2,…,Xin),则第i个驾驶人到正负理想解的马氏距离分别为:
Figure FDA0002836446150000028
Figure FDA0002836446150000029
其中,Ai为第i个评估对象即第i个驾驶人,∑表示相关系数矩阵;
S4.3、驾驶人Ai到正理想解和负理想解的联系向量距离分别表示如下:
Figure FDA00028364461500000210
Figure FDA00028364461500000211
其中,
Figure FDA00028364461500000212
为正同一度系数;
Figure FDA00028364461500000213
为正差异度系数;
Figure FDA00028364461500000214
为正对立度系数;
S4.4、通过马氏距离求得的d(Ai,S+)、d(Ai,S-)和联系向量距离求得的
Figure FDA00028364461500000215
进行无量纲化处理后,定义各待评估对象与正负理想解的合成贴近度为:
Figure FDA00028364461500000216
Figure FDA00028364461500000217
式中,系数α、β满足α+β=1;则各驾驶人与理想解的相对贴近度为:
Figure FDA00028364461500000218
据相对贴近度Ci的大小对驾驶人进行驾驶行为的安全等级排序,Ci值越趋近于1,说明该评估对象的驾驶行为安全等级越高。
3.根据权利要求2所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,S3.1中确定各个评估指标的主观权重pj的过程为:
S3.1.1、基于关键指标建立递阶层次;
S3.1.2、构造判断矩阵,若上一层指标支配的下一层指标有n个,n个指标就构成一个两两比较的判断矩阵A=(aij)n×n,aij表示指标i与指标j的重要性之比;
S3.1.3、若向量ω=[ω12,…,ωn]T满足:Aω=λmaxω,λmax为判断矩阵A的最大特征根,ω为对应的特征向量,归一化的特征向量可以作为主观权重pj
S3.1.4、在计算得到主观权重pj之后,进行一致性检验,一致性指标:
Figure FDA0002836446150000031
其中,n为矩阵阶数;随机一致性比例:
Figure FDA0002836446150000032
RI为平均随机一致性指标,当随机一致性比例CR<0.1时,认为判断矩阵具有一致性,若一致性检验通过,说明该主观权重pj符合要求,否则需要调整判断矩阵中的aij重新计算主观权重pj
4.根据权利要求2所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,S3.2、采用熵值法确定各个评估指标的客观权重qj的过程为:
S3.2.1、计算第j项指标下,第i个驾驶人的特征比重:
Figure FDA0002836446150000033
其中,i=1,2,…,m,m为驾驶人总数,j=1,2,…,n,n为评估指标总数;
S3.2.2、计算第j项指标的熵值:
Figure FDA0002836446150000034
Figure FDA0002836446150000035
k为系数,m为驾驶人总数;
S3.2.3、计算第j项指标的差异性系数:gj=1-ej
S3.2.4、确定第j项指标的客观权重:
Figure FDA0002836446150000036
5.根据权利要求1、2、3或4所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,S2中相关性分析的过程如下:根据相关系数公式计算不同指标之间的相关性:
Figure FDA0002836446150000037
Sxy为协方差,表示为
Figure FDA0002836446150000038
Sx为某一种指标的标准差,表示为
Figure FDA0002836446150000041
Sy为另一种指标的标准差,表示为
Figure FDA0002836446150000042
Xi表示某一个指标的实际值,
Figure FDA0002836446150000043
表示该指标的平均值,Yi表示另一个指标的实际值,
Figure FDA0002836446150000044
表示另一个指标的平均值,n为采样的数据样本数;当相关系数的绝对值越大,相关性越强;相关系数越接近于0,相关性越弱。
6.根据权利要求5所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,用于安全性评估的关键指标分别是:变换车道视觉感知行为安全性指标LCTR、前方障碍物视觉感知安全性指标OSIA、制动幅值变化率BACV、方向盘转角熵值SE、纵向速度标准差Vstd、纵向加速度变化率LACR、方向盘转角标准差θwstd、横向加速度标准差Lastd和横向车道保持安全性指标LKSR。
7.根据权利要求6所述的一种基于多属性决策的驾驶行为安全性评估方法,其特征在于,获得变换车道视觉感知行为安全性指标LCTR的方法为:将驾驶人换道过程的注视区域划分为5块:车道正前方、车道左前方、车道右前方、左后视镜和右后视镜;通过眼动仪采集的不同驾驶人注视点数据,且利用K-means方法对通过眼动仪采集的不同驾驶人注视点数据进行聚类;最终选择左侧换道驾驶人注视左后视镜和车道正前方的平均次数,右侧换道注视右后视镜和车道正前方的平均次数作为分析驾驶人变换车道视觉感知行为安全性的指标;每一位驾驶人换道过程观测不同位置的次数与样本中规范性要求下变换车道注视行为水平均值进行对比,得到驾驶人变换车道视觉感知行为安全性指标
Figure FDA0002836446150000045
F为各驾驶人换道过程观测不同区域的总次数;A为规范性要求下变换车道注视行为水平均值。
CN202011472843.3A 2020-12-15 2020-12-15 一种基于多属性决策的驾驶行为安全性评估方法 Pending CN112541632A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011472843.3A CN112541632A (zh) 2020-12-15 2020-12-15 一种基于多属性决策的驾驶行为安全性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011472843.3A CN112541632A (zh) 2020-12-15 2020-12-15 一种基于多属性决策的驾驶行为安全性评估方法

Publications (1)

Publication Number Publication Date
CN112541632A true CN112541632A (zh) 2021-03-23

Family

ID=75018730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011472843.3A Pending CN112541632A (zh) 2020-12-15 2020-12-15 一种基于多属性决策的驾驶行为安全性评估方法

Country Status (1)

Country Link
CN (1) CN112541632A (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113379333A (zh) * 2021-07-22 2021-09-10 交通运输部公路科学研究所 一种辅助驾驶营运车辆智能化等级评价方法
CN113487197A (zh) * 2021-07-12 2021-10-08 吉林大学 一种基于topsis评价模型的车联网信息筛选方法
CN113722566A (zh) * 2021-08-18 2021-11-30 南斗六星系统集成有限公司 一种自动驾驶车辆功能稳定性评定的方法
CN113942520A (zh) * 2021-10-27 2022-01-18 昆明理工大学 一种驾驶人可靠度计算方法
CN114648241A (zh) * 2022-03-31 2022-06-21 西安热工研究院有限公司 一种基于topsis算法的矿用卡车行驶安全性评价方法
CN114662897A (zh) * 2022-03-17 2022-06-24 Oppo广东移动通信有限公司 数据处理方法、装置、电子设备及存储介质
CN114688981A (zh) * 2022-03-15 2022-07-01 中南大学 一种复合钢筋覆层厚度识别及厚度均匀性评价方法、设备及介质
CN115240800A (zh) * 2022-09-26 2022-10-25 北京泽桥医疗科技股份有限公司 一种基于大数据平台的医疗数据智能分析执行方法
CN115481906A (zh) * 2022-09-22 2022-12-16 中航机载系统共性技术有限公司 一种生产制造产线评估方法、装置、设备及介质
CN117152411A (zh) * 2023-11-01 2023-12-01 安徽蔚来智驾科技有限公司 视线校准方法、控制装置及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109711691A (zh) * 2018-12-17 2019-05-03 长安大学 一种基于熵权模糊综合评价模型的驾驶风格评价方法
CN110826848A (zh) * 2019-09-19 2020-02-21 安徽百诚慧通科技有限公司 基于层次分析法的驾驶人风险评估方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109711691A (zh) * 2018-12-17 2019-05-03 长安大学 一种基于熵权模糊综合评价模型的驾驶风格评价方法
CN110826848A (zh) * 2019-09-19 2020-02-21 安徽百诚慧通科技有限公司 基于层次分析法的驾驶人风险评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李扬: "驾驶行为安全性多属性评价方法及应用研究", 中国博士学位论文全文数据库 工程科技Ⅱ辑, no. 2017, pages 16 - 106 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113487197A (zh) * 2021-07-12 2021-10-08 吉林大学 一种基于topsis评价模型的车联网信息筛选方法
CN113379333A (zh) * 2021-07-22 2021-09-10 交通运输部公路科学研究所 一种辅助驾驶营运车辆智能化等级评价方法
CN113722566A (zh) * 2021-08-18 2021-11-30 南斗六星系统集成有限公司 一种自动驾驶车辆功能稳定性评定的方法
CN113942520A (zh) * 2021-10-27 2022-01-18 昆明理工大学 一种驾驶人可靠度计算方法
CN114688981A (zh) * 2022-03-15 2022-07-01 中南大学 一种复合钢筋覆层厚度识别及厚度均匀性评价方法、设备及介质
CN114688981B (zh) * 2022-03-15 2023-01-31 中南大学 一种复合钢筋覆层厚度识别及厚度均匀性评价方法、设备及介质
CN114662897A (zh) * 2022-03-17 2022-06-24 Oppo广东移动通信有限公司 数据处理方法、装置、电子设备及存储介质
CN114648241A (zh) * 2022-03-31 2022-06-21 西安热工研究院有限公司 一种基于topsis算法的矿用卡车行驶安全性评价方法
CN115481906A (zh) * 2022-09-22 2022-12-16 中航机载系统共性技术有限公司 一种生产制造产线评估方法、装置、设备及介质
CN115481906B (zh) * 2022-09-22 2024-04-16 中航机载系统共性技术有限公司 一种生产制造产线评估方法、装置、设备及介质
CN115240800A (zh) * 2022-09-26 2022-10-25 北京泽桥医疗科技股份有限公司 一种基于大数据平台的医疗数据智能分析执行方法
CN117152411A (zh) * 2023-11-01 2023-12-01 安徽蔚来智驾科技有限公司 视线校准方法、控制装置及存储介质

Similar Documents

Publication Publication Date Title
CN112541632A (zh) 一种基于多属性决策的驾驶行为安全性评估方法
CN113401143B (zh) 一种基于驾驶风格和意图的个性化自适应轨迹预测方法
Lyu et al. Using naturalistic driving data to identify driving style based on longitudinal driving operation conditions
Guo et al. Improved car-following strategy based on merging behavior prediction of adjacent vehicle from naturalistic driving data
CN112258830B (zh) 车辆编队行驶的可靠性评估方法及其应用
CN112991685A (zh) 考虑驾驶人疲劳状态影响的交通系统风险评估及预警方法
Wang et al. ARIMA Model and Few‐Shot Learning for Vehicle Speed Time Series Analysis and Prediction
CN113761715A (zh) 一种基于高斯混合和隐马尔可夫的个性化车辆跟驰模型的建立方法
CN112990563A (zh) 一种高速公路追尾事故风险实时预测方法
CN116383678B (zh) 一种营运客车异常变速行为常发性路段识别方法
CN116859352A (zh) 一种路侧毫米波雷达设备运行状态评估方法及系统
CN111222740A (zh) 基于三角模糊分析法和熵值法的制动踏板感觉评价方法
CN116975671A (zh) L3级自动驾驶车辆驾驶员信任度评估方法及系统
Anil et al. Driver behavior analysis using K-means algorithm
Zhong et al. Evaluating the Efficacy of Real-Time Connected Vehicle Basic Safety Messages in Mitigating Aberrant Driving Behaviour and Risk of Vehicle Crashes: Preliminary Insights from Highway Scenarios
Song et al. Vehicle lane change trajectory learning and prediction model considering vehicle interactions and driving styles in highway scene
CN110610611B (zh) 面向混行交通流中智能网联车辆的行车安全性评价方法
CN113252057A (zh) 一种基于高德导航数据的驾驶倾向性辨识方法和系统
Wang et al. An Improved Grey Relational Theory Evaluation Method: Considering the Comprehensive Performance of Autonomous Vehicles in Virtual Test
Wen et al. Analysis of vehicle driving styles at freeway merging areas using trajectory data
CN118197095B (zh) 一种交通事故的安全预警方法
CN118035691B (zh) 智能驾驶感知质量的量化评价方法、设备和存储介质
CN109118034A (zh) 驾驶行为三维定量综合评价方法
Han et al. The Critical Scenario Extraction and Identification Method for ICV Testing
Ruiwen et al. Research on a Comprehensive Evaluation System for Driving Behavior by Introducing Fuzzy Evaluation Based on Natural Driving Data

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