基于差分方程和元胞自动机交通阻塞模型.docVIP

基于差分方程和元胞自动机交通阻塞模型 　　摘要：在城市化进程发展快速的当下，车道被占用会带来一系列的消极影响。通过计算修正的道路通行能力值，得到事故所处横断面实际通行能力会下降并持续一段时间。通过t检验可以得到所占车道不同对横断面实际通行能力具有显著性差异。然后建立了分段差分方程模型，得到了交通事故所影响的路段车辆排队长度与其他变量间的关系。当交通事故发生是红灯或绿灯时，估算了车辆排队到达上游路口所需的时间，最后通过元胞自动机的模拟仿真得到了很好的验证。 　　关键词：通行能力；ARMA模型；分段差分方程；元胞自动机 　　中图分类号：TB文献标识码：Adoi：10.19311/ki2018.04.088 　　1引言 　　车道被占用是指由于交通事故等各种原因，导致车道或道路横断面通行能力在单位时间内降低的现象。在城市化进程发展快速的当下，车道被占问题情况复杂，车道被占用极大程度上影响人民的正常工作生活。为此，根据2013年全国大学生数学建模竞赛A题，采用数学方法分析建立相关模型，估算了其影响程度，将为交通部门处理类似问题提供参考。 　　2模型假设 　　（1） 假设视频 1，视频2 中的车祸发生于同一横断面面，并且都完全占用两条车道。 　　（2） 假设两个小区路口出入的车辆对上游车流量的影响可忽略不计。 　　（3） 假设所有车辆严格遵守交通规则。 　　3模型的建立与求解 　　3.1数据的预处理 　　实际通行能力是指单位时间内通过的标准小汽车当量数，考虑到视频中的实际情况，在本文中只考虑四轮及四轮以上的机动车、摩托车及电瓶车。并且对公交车，摩托车等进行当量化的计算。通过公路工程技术标准，采用标准小汽车当量数计算车型折算系数如表1所示。 　　由于红绿灯以30秒为一个相位进行变化，于是本文以30秒的时间为一个单位间隔对视频1中通过车祸横截面处对车辆进行计数，由于16：42：32为车祸发生时间，所以本文从此时开始计数。观看视频1我们可以得知，部分时刻的画面存在缺失。对缺失的数据我们采用MATLAB进行拟合得到。 　　可以建立t时刻的标准小汽车当量数Ct的计算公式如下： 　　Ct=St1+0.5St2+1.5St3 　　其中Sti表示t时刻第i种车型的数量。 　　t时刻实际通行能力Yt的计算公式为： 　　Yt=120Ct 　　3.2模型的建立 　　ARMA模型（auto regressive moving average model）自回归滑动平均模型，这种方法是研究平稳随机过程有理谱的典型方法，适用于很大一类实际问题。首先对数据进行AIC检验，通过MATLAB编程得出当p=3，q=3时，AIC达到最小值即AICmin=9923776，最后进行chi2检验，h=0说明检验通过，拟合的曲线为： 　　Yt=2.48Xt-1-2.21Xt-2+0.73Xt-3+εt-2598εt-1+2.41εt-2-0.7968εt-3 　　实际通行能力是基本通行能力与各修正系数的连乘积，其计算公式如下： 　　Y=Y0*FEQ*FLT*FRT*FLA*FFR 　　其中基本通行能力Y0是理想情况下每条车道单位时间内能够通过的最大交通量，其计算公式如下： 　　Y0=3600t0=3600l0v3.6=1000vl0pcuh 　　其中l0=ld+lb+ls+lc 　　t0表示驾驶员反应时间，取0.5s， l0表示?头最小间距，ld表示驾驶员在反应时间内车辆行驶的距离，lb表示车辆的制动距离，ls表示车辆间安全距离，lc表示车辆平均长度。 　　由条件知车道宽度为3.25m，车辆平均长度为各车车长及其所占车流量比例乘积和，结合视频有v=33km/h，ld=5.56，lb=10，ls=5，lc=3.42。求得Y0=1388.66。 　　各修正系数为：主支路流量不平衡影响修正系数FEQ=1-0.32lnx，左转修正系数FLT=1-0.40x，右转修正系数FRT=1+0.1x，交叉平面口大型车汇入修正系数FLA=1+0.20x，结合实际观测出的各时段横向干扰修正系数的值FFR=0.8，综上可得到： 　　Y=Y0×[1-0.32×ln（qaq23）]×（1+0.1×q23qa）×（1+0.2×qlqa）×0.8 　　通过MATLAB编程求解得到模型值与通过视频统计的实际值做出对比如图1所示，从图1可以看出，该模型计算的实际通行能力值是真实可靠的。交通事故发生至撤离期间， 事故所处横断面实际通行能力会下降且会持续一段时间。 　　3.3差异性检验 　　为了说明Yt1，Yt2的数据之间是否存在显著性差异，我们通过SPSS软件对其进行T检验，首先假设Yt1，Yt2之间不存在显著性差异，检验结果如表2所