CN103544850B - 车联网中基于车间距离概率分布的碰撞预测方法 - Google Patents

Abstract

本发明公开了一种基于车间距概率分布的高速路模型下车辆碰撞的预测方法，包括以下步骤：车辆周期性(10HZ)广播当前运动状态信息Beacons（速度、加速度、GPS）；动态计算周围环境的车辆密度，构建车间距分布概率模型；根据自身和前方紧相邻车辆的运动状态，动态计算两者同时紧急刹车情况下避免发生碰撞所需的最小安全距离；根据车间距概率分布，计算相邻两车发生碰撞的概率（车间距小于最小安全距离概率）；构建多车辆发生碰撞的Markov链和状态转移矩阵，评估整个路段某时刻发生车辆碰撞次数的期望。该方法创新度高，扩展性强，较好的弥补了当前基于车辆定位预测碰撞算法中GPS数据精度不足和不稳定的缺陷，在GPS卫星信号盲区本方法的作用尤为突出，应用前景十分广阔。

Description

车联网中基于车间距离概率分布的碰撞预测方法

技术领域

本发明属于智能交通系统的车联网技术领域，具体涉及一种在紧急情况下基于车间距概率分布的多车碰撞预测方法。车辆碰撞预测方法，可以用于实时计算单独车辆行驶的安全系数，也可以用于评估整条路段在遭遇突发事故时引发次生碰撞的风险。

背景技术

近年来，车联网(VehicularAd-hocNetwork,VANET)作为智能交通领域中一种新兴技术引起了众多汽车厂商和学者的广泛关注，联邦通讯委员会(FederalCommunicationCommission，FCC)专门为车联网应用分配了75MHz的通信频段。在VANET中,每个移动车辆均具备无线通信的能力，正常情况下周期性广播包含自身运动状态和GPS位置信息的Beacons消息，当车辆遭遇危险情形时会触发并广播警告消息。通过实时分享各自运动状态和及时告知危险事件达到提高行车安全的目的。

基于车联网技术，实时预测相邻车辆碰撞概率是智能交通领域中的一项强烈需求，可以有效保障和提高道路交通的安全。目前常见预测碰撞的方法大部分是基于车辆GPS定位技术，通过Beacon消息获取周围车辆的位置信息，结合自身GPS坐标计算相对车距。同时根据相邻车辆的运动状态计算是否会发生碰撞。由于当前车载GPS定位系统存在较大误差，造成此类预测车辆碰撞的方法可靠性较差。尤其在GPS信号覆盖不到或者信号较弱的路段，此类方法将会失去作用。

虽然文献[C.Garcia-Costa,etc,“Astochasticmodelforchaincollisionsofvehiclesequippedwithvehicularcommunications”,IEEETransactionsonIntelligentTransportationSystems,vol.13,no.2,Jun.2012]给出了一种预测车辆连环碰撞的随机模型，但是由于该模型中各种输入变量均假设为随机变量，无法有效应用到实际车联网中实时预测多车碰撞的问题。为了摆脱对GPS定位精度强依赖，本发明首次采用基于车间距概率分布的方法来预测车辆的碰撞事件，并充分利用车联网中的两类消息的作用，提高的预测碰撞的准确度。

发明内容

本发明目的在于提供一种在紧急情况下基于车间距概率分布的多车碰撞预测方法，有效提高车联网中行车碰撞预测的准确度，保障了道路交通安全，同时摆脱了GPS定位精度不足的制约，实用范围更广。

为了解决现有技术中的这些问题，本发明提供的技术方案是：

一种车联网中基于车间距离概率分布的碰撞预测方法，假设车联网整个路段上共有N辆编号依次为V₁,V₂,...,V_N的车相继前行，每一个移动车辆均具备无线发射和接受信息的能力，其特征在于所述方法包括以下步骤：

(1)车辆周期性的向邻居车辆广播包含自身实时运动状态的Beacons消息，其中实时运动状态信息包括速度v,加速度a,运动方向以及GPS位置向量X；

(2)车联网中的每个车辆通过解析来自邻居车辆的Beacons消息，获取周围环境中邻居车辆的实时运动状态信息，动态计算周围环境的车辆分布密度λ，构建车间距概率分布模型；

(3)前后相邻同向行驶的车辆V_i-1和V_i根据两者的最新速度v_i-1和v_i、加速度a_i-1和a_i、最大刹车速度a_max,i-1和a_max,i，根据物理学运动模型动态计算两者同时紧急刹车情况下为避免碰撞所需保持的最小安全距离d_ms；

(4)根据车辆V_i与前方相邻车辆V_i-1的车间距概率分布模型，计算相邻两车发生碰撞的概率；所述V_i和V_i-1发生碰撞的概率即为V_i与V_i-1的车间距小于最小安全间距d_ms的概率；1<＝i<＝N；

(5)根据相邻两车发生碰撞的概率构建多车辆发生碰撞的Markov链和状态转移矩阵，评估整个路段某时刻发生车辆碰撞次数的期望。

优选的技术方案是：所述方法步骤(2)中邻居车辆的Beacons消息解析出来的GPS位置向量数据用于局部计算周围环境的实时车辆分布密度λ；通过侦听范围内最前和最后的两个邻居车辆GPS位置向量计算车辆的感知范围总长度L,各个车辆进而可以实时计算出所处局部环境的车辆分布密度λ＝N/L。

优选的技术方案是：所述方法步骤(3)中前后相继行驶车辆V_i-1和V_i的刹车延迟t_res包含接受到前方事故车辆V_j触发的告警消息传输延迟和司机采取紧急刹车的反应延迟两部分。

优选的技术方案是：相邻车辆V_i-1和V_i的最小安全行驶距离d_ms的计算方法根据刹车延迟t_res结束时刻两车的初始运动状态，分为三种情况：

i)j＝i-1,即V_i-1为车队中首先触发告警消息并向后传输的车辆；

ii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度大于V_i-1的速度；

iii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度小于等于V_i-1的速度。

优选的技术方案是：所述方法中当其中编号为V_j的车遭遇紧急事故并将此告警消息向后广播，后继的N-j辆车在收到此消息后几乎同时采取紧急刹车动作会引发连环的多车碰撞，基于后续车辆发生连环碰撞过程的所有状态，构建包含了(N-j+1)(N-j+2)/2种状态的齐次Markov链和对应的状态转移矩阵P；通过N-j次自乘运算生成新矩阵P^N-j。矩阵元素P^N-j(1,(N-j+1)(N-j+2)/2-k)表示N-j个后续车辆中最终共有k个车辆发生碰撞的概率，其中0≤k≤N-j；计算出由V_j引发的后续N-j辆车发生碰撞次数的期望。

优选的技术方案是：所述方法步骤(5)中构建Markov链的状态集合为C＝(c_0,0,c_1,0,c_0,1,...,c_N-j,0,c_N-j-1,1,c_N-j-2,2,...,c_1,N-j-1,c_0,N-j)，集合大小为(N-j+1)(N-j+2)/2。其中c_0,0表示触发告警消息的V_j，定义为初始状态，c_1,0表示从V_j往后的1辆车中有一个发生碰撞，c_N-j-2,2表示从V_j开始的后续N-j辆车中有N-j-2个发生碰撞，依次类推。

本发明提供了一种基于车间距概率分布的高速路模型下车辆碰撞的预测方法，其特征在于所述方法包括以下步骤：(1)移动车辆通过周期性的向邻居节点广播状态信息(Beacons),获取邻居车辆当前的运动状态并计算所处环境下的车辆密度λ；(2)结合自身当前的运动状态，每个车辆动态计算在紧急刹车情况下和前方相邻车辆避免发生碰撞的最小安全行车距离；(3)根据车辆分布密度，动态更新车间距的概率分布模型。结合最小安全行车距离，计算两两相邻车辆紧急情况下发生碰撞的概率；(4)针对相同车道上同向相继行驶的多个车辆，基于连环碰撞发生的所有中间状态构建Markov链和状态转移矩阵，并计算每种最终碰撞状态发生的概率。进而评估整个路段紧急情况下发生车辆连环碰撞次数的总体期望。该方法创新度高，扩展性强，较好的弥补了当前基于车辆定位预测碰撞算法中GPS数据精度不足和不稳定的缺陷，在GPS卫星信号盲区本方法的作用尤为突出，应用前景十分广阔。

本发明利用车联网中周围车辆的实时运动状态信息(Beacons),动态维护相邻车间距的概率分布模型。基于相邻车辆的运动状态计算紧急状况下保证两车避免碰撞所需的最小安全距离，车间距小于最小安全距离的概率即为碰撞概率，并基于Markov链评估整个路段上多车碰撞的危险系数，提高了行车的安全度，该方法具体可以按照如下步骤进行：

步骤一：移动车辆周期性的向周围的车辆广播Beacons消息，同一车辆广播Beacon消息的间隔周期为0.1s。Beacon消息包含了当前车辆最新的运动状态信息，具体包括速度v,加速度a,运动方向以及GPS位置向量X。通过侦听来自周围邻居车辆的Beacons消息，可以获取周围车辆的上述信息。

步骤二：计算各自感知范围总长度L＝|X_head-X_last|和周围环境中车辆分布密度λ＝N/L。其中X_head和X_last表示车辆侦听范围内的最前和最后两个邻居车辆的GPS位置向量，均从邻居车辆的Beacon消息中解析出来，N表示侦听范围内邻居总数；

步骤三：假设前后相邻同向行驶的车辆V_i-1和V_i同时收到来前方危险车辆V_j,j≤i-1触发的警告消息均立即采取紧急刹车，V_i根据自身和其前驱车辆V_i-1的最新速度v_i和v_i-1、加速度a_i-1和a_i、最大刹车速度a_max,i和a_max,i-1，按照物理学运动模型(如图3所示)，计算两者同时紧急刹车情况下为避免碰撞所需保持的最小安全间距d_ms。为了便于分析和计算，本方法中所有车辆的最大刹车加速度均为a_max。根据刹车延迟t_res结束时刻前后两车的初始运动状态，d_ms的计算方法具体分为三种情况：

i)j＝i-1,即V_i-1为车队中首先触发告警消息并向后传输的车辆V_j。此种情况下认为V_j(即V_i-1)立即停止，其后继车辆V_i是否与其相撞完全取决于V_i的刹车距离。此时的最小安全行车距离为：

$\begin{array}{c}{d}_{ms,i}=signal\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)\&CenterDot;\left(v_{0,i}\&CenterDot;t_{res}+\frac{1}{2}{a}_i\&CenterDot;t_{res}^2+\frac{{\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)}^2}{2a_{\max }}\right)\\ +signal\left(-\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)\right)\&CenterDot;\frac{v_{0,i}^2}{2a_i}+s_0,i=j+1,\end{array}$

其中v_0,i表示前方事故车辆V_j触发警告消息那一刻车辆V_i的初速度。signal是信号量函数： $signal\left(x\right)=\left\{\begin{array}{ccc}0,& if& x<0,\\ \frac{1}{2},& if& x=0,\\ 1,& if& x>0.\end{array}\right.,$ s₀表示两辆车发生碰撞的临界距离。

ii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度大于V_i-1的速度，此时车辆V_i的最小安全行车距离d_ms,i为：

$d_{ms,i}=\frac{{(v_{0,i}+a_i\&CenterDot;t_{res})}^2-{(v_{0,i-1}+a_{i-1}\&CenterDot;t_{res})}^2}{2a_{max}}+s_0,i=j+2,...,N,$

iii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度小于等于V_i-1的速度,此时车辆V_i的最小安全行车距离只要满足大于等于碰撞发生的临界距离s₀即可。

步骤四：一维高速路模型下车辆间距严格符合参数为λ的指数分布。按照步骤二中λ的计算方法，辆车V_i与前驱车辆V_i-1的车间距d_i的概率模型表示为 $F(d_i;\mathrm{\&lambda;})=1-e^{-{\mathrm{\&lambda;d}}_i},d_i\&GreaterEqual;0.$

步骤五：车辆V_i和其前驱车辆V_i-1发生碰撞的概率转化为V_i和V_i-1的间距小于最小安全距离d_ms,i的概率。因此车辆V_i与V_i-1发生碰撞的概率p_i可以表示为：

$\begin{array}{c}{p}_i=P\left(d_i<d_{ms,i}\right)\\ =1-e^{-{\mathrm{\&lambda;d}}_{ms,i}},i=j+1,\mathrm{...},N\end{array}$

步骤六：针对整个路段上事故车辆V_j后续的N-j辆车的碰撞状态构建Markov链和对应的状态转移矩阵P。Markov链对应的状态集合C是C＝(c_0,0,c_1,0,c_0,1,...,c_N-j,0,c_N-j-1,1,c_N-j-2,2,...,c_1,N-j-1,c_0,N-j),连续的两个状态转移概率为步骤(3)中的两两相邻车辆发生碰撞的概率p_i和非碰撞概率1-p_i，其中状态集合C大小为(N-j+1)(N-j+2)/2,c_0,0表示触发告警消息的V_j，定义为初始状态，c_1,0表示从V_j往后的1辆车中有一个发生碰撞，c_N-j-2,2表示从V_j开始的后续N-j辆车中有N-j-2个发生碰撞，依次类推,具体构建过程如图4所示。

步骤七：对步骤六中构建好的状态转移矩阵P做N-j次幂运算，即P^N-j。矩阵元素P^N-j(1,(N-j+1)(N-j+2)/2-k)表示N-j个后续车辆中最终共有k个车辆发生碰撞的概率，也就是状态c_k,N-j-k出现的概率，其中0≤k≤N-j。

步骤八：由车辆V_j引发的后续车辆碰撞的次数N_colli最终为：

$N_{colli}=\underset{k=0}{\overset{N-j}{\&Sigma;}}k\&CenterDot;P^{N-j}(1,(N-j+1\left)\right(N-j+2)/2-k),j=1,2...,N$

相对于现有技术中的方案，本发明的优点是：

本发明技术方案预测车辆碰撞的准确率高：由于当前GPS数据定位精度不够，基于GPS定位计算车间距的预测碰撞的方法准确率较低，本发明首次采用动态维护相邻车距的实时概率模型的方法预测车辆碰撞，取得了较高的准确率，尤其是在没有GPS信号覆盖或者信号不稳定的路段。

本发明技术方案可扩展性强：针对一维高速路交通模型，车间距较好地服从参数为λ指数分布。如果随着车辆移动，交通道路模型发生改变时，可以根据需要将车间距分布的概率模型替换成合适的类型(泊松分布等)，其他步骤不需要做任何改动，具有较好的可扩展性和灵活性。

本发明技术方案对行车安全具有现实意义：本方法除了可以实时预测单个行车紧急情况下与前车碰撞的概率，还可以针对整个路段上的多个相继行驶车辆评估紧急情况下多车碰撞的风险。

附图说明

下面结合附图及实施例对本发明作进一步描述：

图1为本发明说明框图；

图2为本发明中一维高速路碰撞模型图；

图3为本发明中计算最小安全距离的示意图；

图4为本发明建立多车碰撞状态的马尔可夫链示意图；

图5为本发明中3辆车碰撞状态的马尔科夫链和状态转移矩阵图；

图6为本发明实施例仿真实验的预测碰撞准确性率图。

具体实施方式

以下结合具体实施例对上述方案做进一步说明。应理解，这些实施例是用于说明本发明而不限于限制本发明的范围。实施例中采用的实施条件可以根据具体厂家的条件做进一步调整，未注明的实施条件通常为常规实验中的条件。

实施例

如图1所示，本发明提供的车联网中一种基于车间距概率分布的碰撞预测方法，所述车间距概率分布模型随行车周围邻居车辆的运动状态实时动态维护，所述方法包括以下步骤：

(1)假设整个路段上共有N辆编号依次为V₁,V₂,...,V_N的车相继前行，每一个移动车辆均周期性的向邻居车辆广播包含自身实时运动状态Beacons消息。通过解析来自邻居车辆的Beacons消息，获取周围环境中邻居节点的实时运动状态和车辆分布密度λ。其中运动状态信息具体包括速度v,加速度a,运动方向以及GPS位置向量X；

(2)前后相邻同向行驶的车辆V_i-1和V_i根据两者的最新速度v_i-1和v_i、加速度a_i-1和a_i、最大刹车速度a_max,i-1和a_max,i，根据物理学运动模型，计算两者同时紧急刹车情况下为避免碰撞所需保持的最小安全间距d_ms；

(3)车辆V_i动态更新与前方相邻车辆V_i-1的车间距分布概率模型，此处在高速路交通模型下，车间距概率上较好地符合指数分布。根据指数分布模型，V_i计算与V_i-1的车间距小于最小安全间距d_ms的概率，即为V_i和V_i-1发生碰撞的概率；

(4)如果编号为V_j的车遭遇紧急事故将触发告警消息向后广播(如图2所示)，后继的N-j辆车在收到此消息采取紧急刹车。根据此N-j辆车中的每一个是否与其前方相邻车辆发生碰撞，后续车辆中发生碰撞的最终数目对应了N-j+1种情况。根据计算出的两两相邻车辆发生碰撞的概率，基于后续车辆发生连环碰撞过程的所有状态，构建包含了(N-j+1)(N-j+2)/2种状态的齐次Markov链和对应的状态转移矩阵P。通过N-j次自乘运算生成新矩阵P^N-j。矩阵元素P^N-j(1,(N-j+1)(N-j+2)/2-k)表示N-j个后续车辆中最终共有k个车辆发生碰撞的概率，其中0≤k≤N-j。进而计算出由V_j引发的后续N-j辆车发生碰撞次数的期望。

所述方法步骤(1)中邻居车辆的Beacons消息解析出来的GPS位置向量数据仅用于局部计算周围环境的实时车辆密度λ。每辆车的感知范围总长度L通过L＝|X_head-X_last|计算，周围环境中车辆分布密度λ通过λ＝N/L计算。其中X_head和X_last表示侦听范围内最前和最后的两个邻居车辆的GPS位置向量，N表示侦听范围内邻居总数。

所述方法步骤(2)中前后相继行驶车辆V_i-1和V_i的刹车延迟t_res包含接受到前方事故车辆V_j触发的告警消息传输延迟和司机采取紧急刹车的反应延迟两部分。其中告警消息的传输延迟通常采用0.1s,人的反应时间采用0.9s,即总的刹车延迟t_res＝1s。

所述方法步骤(3)中相邻车辆V_i-1和V_i的最小安全行驶距离d_ms的计算方法根据刹车延迟t_res结束时刻两车的初始运动状态，具体分为三种情况：

i)j＝i-1,即V_i-1为车队中首先触发告警消息并向后传输的车辆；

ii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度大于V_i-1的速度；

iii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度小于等于V_i-1的速度；

所述方法步骤(4)中构建Markov链的状态集合为C＝(c_0,0,c_1,0,c_0,1,...,c_N-j,0,c_N-j-1,1,c_N-j-2,2,...,c_1,N-j-1,c_0,N-j),根据步骤(3)中计算出来的两两相邻车辆发生碰撞的概率p和非碰撞概率1-p，建立Markov链对应的状态转移矩阵P。其中状态集合C大小为(N-j+1)(N-j+2)/2,c_0,0表示触发告警消息的V_j，定义为初始状态，c_1,0表示从V_j往后的1辆车中有一个发生碰撞，c_N-j-2,2表示从V_j开始的后续N-j辆车中有N-j-2个发生碰撞，依次类推。

本实施例的高速路车辆分布碰撞模型如图2所示。编号为V₁,V₂,...,V_N的N辆车相继朝同一方向行驶。每辆车周期性的向周围邻居节点广播Beacons消息，其中主要包括各自最新的速度v,加速度a,运动方向以及GPS位置向量X等信息。当其中一车(假设编号为V_j)遭遇突发情况时立即停止，并紧急触发警告消息并将其向后继车辆广播。后继的N-j辆车接受到此警告消息后立即采取紧急刹车。本实施例给出预测相邻车辆发生碰撞的方法以及计算后续N-j辆车中发生碰撞总数N_colli的过程。

具体的预测方法过程如下：

步骤一：移动车辆广播Beacons消息周期设置为0.1s。Beacon消息包含了当前车辆最新的速度v,加速度a,运动方向以及GPS位置向量X。通过侦听来自周围邻居车辆的Beacons消息，可以获取周围车辆的上述信息。

步骤二：每个移动车辆动态计算各自感知范围的总长度L＝|X_head-X_last|和周围的车辆密度λ＝N/L。其中X_head和X_last表示车辆各自侦听范围内最前和最后两个邻居车辆的GPS位置向量，N表示侦听范围内邻居总数；

步骤三：V_i-1和V_i表示位于V_j后面两相邻车辆。V_i根据自身和其前驱车辆V_i-1的最新速度v_i和v_i-1、加速度a_i和a_i-1、最大刹车速度a_max,i和a_max,i-1，按照物理学运动模型，计算两者同时紧急刹车情况下为避免碰撞所需保持的最小安全间距d_ms。为了便于分析和计算，本实施例中所有车辆具有相同的最大刹车加速度a_max。根据刹车延迟t_res结束时刻V_i和V_i-1的运动状态，d_ms,i的计算方法具体分为三种情况：

i)j＝i-1,V_i紧随事故车辆V_j。此种情况下认为V_j(即V_i-1)立即停止，此时V_i的最小安全行车距离d_ms,i为：

$\begin{array}{c}{d}_{ms,i}=signal\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)\&CenterDot;\left(v_{0,i}\&CenterDot;t_{res}+\frac{1}{2}{a}_i\&CenterDot;t_{res}^2+\frac{{\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)}^2}{2a_{\max }}\right)\\ +signal\left(-\left(v_{0,i}+a_i\&CenterDot;t_{res}\right)\right)\&CenterDot;\frac{v_{0,i}^2}{2a_i}+s_0,i=j+1,\end{array}$

其中刹车延迟t_res＝1s,相邻车辆的发生碰撞的临界车间距离s₀＝01.m，加速度a∈[4,8]m/s²,最大刹车加速度为a_max＝8m/s²，初速度v₀∈[15,32]m/s；

ii)j＜i-1,且刹车延迟t_res结束初速度满足v_0,i＞v_0i,-1，此时车辆V_i的最小安全行车距离d_ms,i为：

$d_{ms,i}=\frac{{(v_{0,i}+a_i\&CenterDot;t_{res})}^2-{(v_{0,i-1}+a_{i-1}\&CenterDot;t_{res})}^2}{2a_{max}}+s_0,$

其中的变量取值与i)中相同；

iii)j＜i-1,且v_0,i≤v_0i,-1,此时车辆V_i的最小安全行车距离d_ms,i＝s₀＝0.1m。

步骤四：根据步骤二中计算出来的λ值，更新辆车V_i与前驱车辆V_i-1的车间距d_i的概率模型 $F(d_i;\mathrm{\&lambda;})=1-e^{-{\mathrm{\&lambda;d}}_i},d_i\&GreaterEqual;0.$

步骤五：计算车辆V_i和其前驱车辆V_i-1发生碰撞的概率p_i为：

$\begin{array}{c}{p}_i=P\left(d_i<d_{ms,i}\right)\\ =1-e^{-{\mathrm{\&lambda;d}}_{ms,i}},i=j+1,\mathrm{...},N\end{array}$

步骤六：为车辆V_j后续N-j辆车的碰撞状态构建Markov链和对应的状态转移矩阵P。本实施例以V_j＝V_N-2为例，即触发事故的车辆V_j后面仅有2辆车V_N-1和V_N。此时的状态集合C＝(c_0,0,c_1,0,c_0,1,c_2,0,c_1,1,c_0,2),基于上述状态集合构建Markov链和状态转移矩阵P如图5所示。其中P为：

$P=\left(\begin{array}{cccccc}0& p_{N-1}& 1-p_{N-1}& 0& 0& 0\\ 0& 0& 0& p_N& 1-p_N& 0\\ 0& 0& 0& 0& p_N& 1-p_N\\ 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0\end{array}\right)$

步骤七：对步骤六中构建好的状态转移矩阵P做2次幂运算，即P²。矩阵元素P²(1,2-k)表示2个后续车辆中最终共有k个车辆发生碰撞的概率，也就是状态c_k,2-k出现的概率，其中0≤k≤2。

步骤八：由车辆V_j引发的后续车辆碰撞的次数N_colli最终为：

$N_{colli}=\underset{k=0}{\overset{N-j}{\&Sigma;}}k\&CenterDot;P^{N-j}(1,(N-j+1\left)\right(N-j+2)/2-k),j=1,2...,N,$

以本实施例中的V_j＝V_N-2以为例，则

为验证本发明预测车辆碰撞的准确性，本发明使用VanetMobiSim工具仿真模拟高速路模型下的车辆碰撞实验。首先使用VanetMobiSim生成记录长达一小时车辆运动轨迹的trace文件，然后在trace文件中设置100次紧急消息触发事件，按照本发明中碰撞预测方法计算每次理论碰撞期望值N_colli，并统计实际trace文件中车辆碰撞数目。并通过多次修改生成trace文件时使用的车道长度模拟多种不同的交通负载情景。仿真实验中用到的参数如表1所示：

表一.仿真实验参数列表

实验结果如图6所示，实验结果表明，本发明基于车间距概率分布的预测碰撞准确率接近95％，对于评估行车安全性具有重大现实意义。

上述实例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人是能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明精神实质所做的等效变换或修饰，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种车联网中基于车间距离概率分布的碰撞预测方法，假设车联网整个路段上共有N辆编号依次为V₁,V₂,...,V_N的车相继前行，每一个移动车辆均具备无线发射和接受信息的能力，其特征在于所述方法包括以下步骤：

(1)车辆周期性的向邻居车辆广播包含自身实时运动状态的Beacons消息，其中实时运动状态信息包括速度v,加速度a,运动方向以及GPS位置向量X；

(2)车联网中的每个车辆通过解析来自邻居车辆的Beacons消息，获取周围环境中邻居车辆的实时运动状态信息，动态计算周围环境的车辆分布密度λ，构建车间距概率分布模型；

(3)前后相邻同向行驶的车辆V_i-1和V_i根据两者的最新速度v_i-1和v_i、加速度a_i-1和a_i、最大刹车速度a_max,i-1和a_max,i，根据物理学运动模型动态计算两者同时紧急刹车情况下为避免碰撞所需保持的最小安全距离d_ms；

(4)根据车辆V_i与前方相邻车辆V_i-1的车间距概率分布模型，计算相邻两车发生碰撞的概率；所述V_i和V_i-1发生碰撞的概率即为V_i与V_i-1的车间距小于最小安全间距d_ms的概率；1≤i≤N；

(5)根据相邻两车发生碰撞的概率构建多车辆发生碰撞的Markov链和状态转移矩阵，评估整个路段某时刻发生车辆碰撞次数的期望。

2.根据权利要求1所述的方法，其特征在于所述方法步骤(2)中邻居车辆的Beacons消息解析出来的GPS位置向量数据用于局部计算周围环境的实时车辆分布密度λ；通过侦听范围内最前和最后的两个邻居车辆GPS位置向量计算车辆的感知范围总长度L,各个车辆进而可以实时计算出所处局部环境的车辆分布密度λ＝N/L。

3.根据权利要求1所述的方法，其特征在于所述方法步骤(3)中前后相继行驶车辆V_i-1和V_i的刹车延迟t_res包含接受到前方事故车辆V_j触发的告警消息传输延迟和司机采取紧急刹车的反应延迟两部分。

4.根据权利要求3所述的方法，其特征在于相邻车辆V_i-1和V_i的最小安全行驶距离d_ms的计算方法根据刹车延迟t_res结束时刻两车的初始运动状态，分为三种情况：

i)j＝i-1,即V_i-1为车队中首先触发告警消息并向后传输的车辆；

ii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度大于V_i-1的速度；

iii)j＜i-1,且刹车延迟t_res结束时刻V_i的速度小于等于V_i-1的速度。

5.根据权利要求1所述的方法，其特征在于所述方法中当其中编号为V_j的车遭遇紧急事故并将此告警消息向后广播，后继的N-j辆车在收到此消息后几乎同时采取紧急刹车动作会引发连环的多车碰撞，基于后续车辆发生连环碰撞过程的所有状态，构建包含了(N-j+1)(N-j+2)/2种状态的齐次Markov链和对应的状态转移矩阵P；通过N-j次自乘运算生成新矩阵P^N-j，矩阵元素P^N-j(1,(N-j+1)(N-j+2)/2-k)表示N-j个后续车辆中最终共有k个车辆发生碰撞的概率，其中0≤k≤N-j；计算出由V_j引发的后续N-j辆车发生碰撞次数的期望。

6.根据权利要求1所述的方法，其特征在于所述方法步骤(5)中构建Markov链的状态集合为C＝(c_0,0,c_1,0,c_0,1,...,c_N-j,0,c_N-j-1,1,c_N-j-2,2,...,c_1,N-j-1,c_0,N-j)，集合大小为(N-j+1)(N-j+2)/2；其中c_0,0表示触发告警消息的V_j，定义为初始状态，c_1,0表示从V_j往后的1辆车中有一个发生碰撞，c_N-j-2,2表示从V_j开始的后续N-j辆车中有N-j-2个发生碰撞，依次类推。