1、数学建模A题车道被占用对城市道路通行能力的影响参赛队员 (打印并签名) :1、徐胜杰2、包小红3、冯金慧指导教师或指导教师组负责人 (打印并签名): 刘志伟 摘要本文研究车道被占用对城市道路通行能力的影响的问题,根据所给的附件,运用数据统计和回归拟合方法,建立了实际道路通行能力模型,运用Excel、SPSS软件进行求解和作图,进而得出了同一横断面交通事故所占车道不同对该横断面实际通行能力的影响差异和变化过程。针对问题一,首先,使用数理统计方法,分别统计出不同车型的车辆;其次,将统计出来的车型数量换算成标准的小车数量;最后,根据换算后的小车数量来计算出道路的实际通行能力并用Excel作成折线图,
2、直观描述事故所处横断面实际通行能力的变化过程。针对问题二,采用和问题一相同的方法,用Excel绘制实际通行能力的变化图,然后与问题一的结果相比较,并绘制出折线图,通过比较得出同一横断面交通事故所占车道不同对该横断面实际通行能力影响有较大的差异。视频2(事故发生在车道一和车道二)中统计的事故所处横断面的通行能力普遍高于视频1(事故发生在车道二和车道三)中的横断面的通行能力。视频2横断面的通行能力高于视频1横断面的通行能力主要因为左车道的车流量高于右车道车流量。 针对问题三,采用线性回归与非线性回归模型相结合,通过SPSS分别判断车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间
3、的关系。因此我们得到非线性函数关系式:,即:车辆排队长度与事故横断面实际通行能力、路段上游车流量呈线性关系,与事故持续时间呈非线性关系。经过检验,该模型和计算结果均是合理的。针对问题四,首先考虑到绿灯的周期性,司机看见绿灯,行车通过十字路口到达下方车道的时间,绿灯期间,上游车流量在大于下方车道通行能力时,会出现排队现象,在一定时间内,拥堵车辆越来越多,排队长度越来越长,通过这个差值与时间、车长、车距可以建立数学模型来计算车辆排队长度。根据此数学模型计算出到达排队长度140m所需时间约为4.47min。经过检验,该模型和计算结果均是合理的。最后,我们总结了模型的优缺点,并提出了改进方法和推广。
4、关键词:通行能力 数理统计 SPSS软件 回归分析方程 数学模型 一、问题的重述1.1 问题重述车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。由于城市道路具有交通流密度大、连续性强等特点,一条车道被占用,也可能降低路段所有车道的通行能力,即使时间短,也可能引起车辆排队,出现交通阻塞。如处理不当,甚至出现区域性拥堵。车道被占用的情况种类繁多、复杂,正确估算车道被占用对城市道路通行能力的影响程度,将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。视频1(附件1)和视频2(附
5、件2)中的两个交通事故处于同一路段的同一横断面,且完全占用两条车道。请研究以下问题:1.根据视频1(附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。2.根据问题1所得结论,结合视频2(附件2),分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。3.构建数学模型,分析视频1(附件1)中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。4.假如视频1(附件1)中的交通事故所处横断面距离上游路口变为140米,路段下游方向需求不变,路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零,且
6、事故持续不撤离。请估算,从事故发生开始,经过多长时间,车辆排队长度将到达上游路口。1.2 问题分析 1、问题一的分析道路通行能力是指在一定的道路条件,交通条件和服务水平的情况下,单位时间能够通过车道上某截面处的最大交通流量。根据视频1统计出一定时间内通过车辆数,计算出道路的实际通行能力,并作图以观察道路通行能力的变化。2、问题二的分析在问题一的基础上,运用相同的方法,绘制实际通行能力的变化图。比较两幅图的区别分析对横断面实际通行能力影响的差异。3、问题三的分析通过分析可以确定的是车辆长度是因变量,而其他三个量是作为自变量的来影响因变量的变化的。采用线性回归与非线性回归模型相结合,分别判断车辆排
7、队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。4、问题四的分析首先我们考虑到红绿灯的周期性。其次,司机看见绿灯到启动的反应时间为,驾驶车辆通过十字路口到达车道所需时间,结合问题一中我们得到下方车道平均通行能力,可以由此建立数学模型。三、模型的假设与符号说明(一) 模型假设1、假设统计数据真实有效2、假设问题中车辆都各行其道,不抢占车道3、假设上游路段红路灯处来车分布均匀4、假设事故期间不重复发生事故5、假设各车辆排队时保持相等间距6、司机反应时间均为1s(二) 符号说明表1 符号说明表符号说明单位y事故影响的路段车辆排队长度 m事故横断面实际通行能力 pcu/h事故持续
8、时间h路段上游车流量pcu/h上游车流量pcu/h下方车道通行能力pcu/h车长m安全距离m司机反应时间s车辆通过十字路口到达车道所需时间sL排队长度mt到达排队长度的时间min第i条车道的修正系数四、模型的分析、建立与求解4.1 问题一4.1.1 问题一的分析问题一要求根据视频1(附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。道路通行能力是指在一定的道路条件,交通条件和服务水平的情况下,单位时间能够通过车道上某截面处的最大交通流量。首先,根据视频1从交通事故发生至撤离期间选择每个30秒左右的时间间隔来统计出通过横断面不同车型的数量;其次,将统计出来的车型数
9、量换算成标准的小车数量;最后,根据换算后的小车数量来计算出道路的实际通行能力,并通过excel软件进行作图,以便直观的反应和观察。4.1.2 模型的建立与求解首先,根据视频1从交通事故发生至撤离期间选择每个30秒左右的时间间隔来统计出通过横断面不同车型的数量;其次,将统计出来的车型数量换算成标准的小车数量;最后,根据换算后的小车数量来计算出道路的实际通行能力,并通过excel软件进行作图,以便直观的反应和观察。通行能力,指单位时间内通过断面的最大车辆数TC(traffic capacity)=n/t(n为通过车辆数,t是时间)视频1中不同时刻点各种车型的数量与横断面道路实际通行能力状况见下表:
10、表2 视频1中不同时刻点各种车型的数量与横断面道路实际通行能力时间段大车小车换算总和通行能力pcu/h42:32-43:0511011.5125443:06-43:352811126343:36-44:0601212144044:36-44:460101092444:47-45:1808892945:19-45:54099101245:55-46:341910.5105046:35-47:16099101247:17-47:500111199747:51-48:3211011.598548:33-49:0509998149:06-50:050131396250:06-50:39189.5964
11、50:40-51:1201212107651:13-51:4411112.5109851:45-52:240111199052:25-52:4606698252:47-53:2126992553:22-54:05189.587854:06-54:380121298154:39-55:10178.593755:11-55:4309998259:10-59:45271097459:46-01:2201010107401:23-01:54219222475平均:1086pcu/h根据上表绘出不同时点下事故处横断面的实际通行能力折线图如下:图1 视频1中不同时点下事故处横断面的实际通行能力由图1可知:
12、在事故刚发生时,即在16:42:32-16:43:05,事故所处横断面的实际通行能力为1254pcu/h;在16:42:32-16:44:06内,道路的实际通行能力随时间的增加呈上升趋势,在16:44:06时,通行能力达到1440pcu/h;在16:44:06-16:44:46内,道路的实际通行能力随时间的增加呈下降趋势,在16:44:46时,通行能力达到924pcu/h;在16:44:46-17:01:22内,道路的实际通行能力基本维持一个稳定的状态,在970pcu/h附近较小范围的波动;在17:01:23-17:01:54内,道路的实际通行能力呈较大幅度上升至2475pcu/h,即在事故撤
13、离时,道路的通行能力逐渐恢复。4.2 问题二4.2.1 问题二的分析分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。在问题一的基础上,运用相同的方法,绘制实际通行能力的变化图。比较两幅图的区别分析对横断面实际通行能力影响的差异。4.2.2 模型的建立与求解在问题一的基础上,运用相同的方法,绘制实际通行能力的变化图。在与问题一的结果进行比较,进而说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。视频2中不同时刻点各种车型的数量与横断面道路实际通行能力状况见附录二。根据附录二绘出不同时点下事故处横断面的实际通行能力折线图如下: 图2 视频2中不同时点下事故
14、处横断面的实际通行能力由图2可知:在事故刚发生时,即在17:35:17-17:36:47,事故所处横断面的实际通行能力由1680pcu/h下降至960pcu/h;在17:36:47-18:02:47内,事故所处横断面的实际通行能力在不断的上下波动;在18:02:47-18:03:17内,道路的实际通行能力逐步上升,即在事故即将撤离时,道路的通行能力逐渐恢复。将问题一求解结果和问题二求解结果的折线图绘制在同一坐标内,结果如下:图3 视频1、2中不同时点下事故处横断面的实际通行能力由上图可知,同一横断面交通事故所占车道不同对该横断面实际通行能力影响有较大的差异。由图3可知,视频2(事故发生在车道一
15、和车道二)中统计的事故所处横断面的通行能力普遍高于视频1(事故发生在车道二和车道三)中的横断面的通行能力。视频2横断面的通行能力高于视频1横断面的通行能力主要因为左车道的车流量高于其他两根车道。4.3 问题三4.3.1 问题三的分析分析视频1(附件1)中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。因为实际情况的复杂性以及不确定性,很难通过交通流问题来确定交通事故的车辆长度与道路通行能力,事故持续时间,上游车流量之间的关系,但是可以确定的是车辆长度是因变量,而其他三个量是作为自变量的来影响因变量的变化的。采用线性回归与非线性回归模型相结合,通过
16、spss分别判断车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。4.3.2 模型的建立与求解从问题一采集的数据中,我们选取部分数据整理如下: 表3 车辆排队长度与通行能力、持续时间、上游车流量的数据统计表路段车辆排队长度(m)通行能力(pcu/h)持续时间(min)上游车流量(pcu/h)1111332124847418002233183169232380921548424221201440525201241332625951201296726541301328827121201292926681021228102440651978112357 根据表3的数据,我们
17、采用回归分析,通过spss分别判断车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。多元线性回归的数学模型可以表示为,式中是4个待估计的参数,是相互对立且服从一正态分布的随机变量,是随机变量,是非随机变量。 设y为事故影响的路段车辆排队长度,为事故横断面实际通行能力,为事故持续时间,为路段上游车流量。以为因变量,分别为自变量,假设与都分别成线性关系,如图:图4 y与的线性拟合关系图5 y与的线性拟合关系图6 y与的线性拟合关系 图7 y与的立方关系如上图我们得到与成线性关系,与不成线性关系,所以我们假定与成三次关系,对数据进行拟合,得到D图,由此可以得到与成三次关系,
18、所以,我们设方程为,然后我们用SPSS软件确定该方程的系数图8 SPSS确定该方程的系数因此我们得到交通事故所影响车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的非线性函数关系式:所以,车辆排队长度与事故横断面实际通行能力、路段上游车流量呈线性关系,与事故持续时间呈非线性关系。这个结果与视频当中观察到的情况基本相符。4.3.3模型检验在表3车辆排队长度与通行能力、持续时间、上游车流量的数据统计表中取出两组数据,如:路段车辆排队长度(m)通行能力(pcu/h)持续时间(min)上游车流量(pcu/h)120144052520130132882712带入所得回归方程中检验,计
19、算得到车辆排队长度与实际非常接近,说明该模型合理。4.4问题四4.4.1 问题四的分析首先我们考虑到红绿灯的周期性,上方需间隔一个相位才会来车,因此,一个周期内流量为题目所给流量的二分之一,而下方车辆则一直在通行,不存在周期问题。其次,司机看见绿灯到启动的反应时间为,驾驶车辆通过十字路口到达车道所需时间,结合问题一中我们得到下方车道平均通行能力,每分钟上游来车减去下方来车为排队车辆,排队车辆乘以车辆长度与车距之和可得到一分钟排队长度,并以此建立数学模型,计算并得到结果。4.4.2 模型的建立与求解查阅资料得到城市主干道三车道修正系数,我们取车道一修正系数为1.0,车道二修正系数为0.8,车道三
20、修正系数为0.6。首先我们假设司机看见绿灯到启动的反应时间为=1s,驾驶车辆通过十字路口到达车道所需时间=5s,假设车长=5m,假设排队时两车相隔安全距离=1m,已知上游车流量=1500pcu/h,在问题一中我们得到下方车道通行能力平均为1478pcu/h,设经过时间t到达排队长度L。考虑到红绿灯的周期性,上方需间隔一个相位才会来车,因此,一个周期内流量为的二分之一,而下方车辆则一直在通行,不存在周期问题,每分钟上游来车减去下方来车为排队车辆,排队车辆乘以车辆长度与车距之和为一分钟排队长度,我们建立数学模型如下:根据此数学模型计算出到达排队长度140m所需时间约为4.47min。4.4.3 模
21、型检验在表3车辆排队长度与通行能力、持续时间、上游车流量的数据统计表中取出两组数据,如:路段车辆排队长度(m)通行能力(pcu/h)持续时间(min)上游车流量(pcu/h)8316923238092154842422带入问题四构建的数学模型中检验,计算达到一定车辆排队长度与相应持续时间,与在视频中统计的数据相差不大,说明该模型是合理的。五、模型的评价(一) 模型的特点 本模型具有较大创新性,通过图像分析数据,直观形象,简单明了,构建清晰的数学模型,方便计算,利于检验。(二) 模型的优点(1)对于问题一和问题二所构建的模型,我们分别求出事故发生前、事故发生时、事故撤离后的通行能力,且得出的结果
22、的相对误差率符合实际要求,这样更能体现建立的模型的可行性与准确性,利用图形的方式清晰的表达了道路发生故障时可流通道路的通行能力的比较,形象生动。(2)对于问题三,我们同样采用线性回归分析和非线性回归分析相结合分析,一方面对采集的数据进行拟合,使模型相对简便化,得到与实际较符合的拟合方程模型,从而确定车辆排队长度与实际通行能力、事故持续时间、路段上游车流量间的关系。推导出四者之间的具体数学关系,具备一定的严密性和准确性,两方面的相互结合达到我们需要的理想结果。(3)对于问题四首先我们根据红绿灯的周期性,得到了一个周期通过红绿灯的车流量,然后我们构建数学模型,简单明了的得出了排队长度与车流量、时间
23、等的简单关系,并得出结论。(三) 模型的缺点(1)在统计车辆数量时,人工计数会造成不同程度误差。(2)在数据的采集上比较单一,并且在视频1和视频2中,有部分画面自动跳转,导致统计有较大困难,这样会造成不可避免的误差。六、模型的改进与推广(一)改进 我们的数学模型都是建立在比较理想的条件下,而对于一些细节问题可能没法考虑,因此这与真实情况会有偏差,本模型需要改进的地方就是每次在建模前尽可能的在现实的生活中多做实验,模拟题目中的情景,统计更多的数据,这样才能与真实情况相对符合,提高模型的精度。(二)推广对于问题一建立的模型,我们可以对交通异常事件对道路通行能力进行预算,进而可以对道路上各条车道上的
24、车流量比例进行修改,从而减少事故对通行能力的影响。此种模型也适用于研究数据量较大,数据收集也比较容易的实际生活问题。对于问题三及问题四的模型,我们可以用于事故对交通流的影响的研究以及对策分析,通过对车辆排队长度与各因素之间的关系,在事故发生后,交通警察可以采取疏散上游车流量,以及缩短交通撤离时间使得事故对交通流的影响降为最低。七、参考文献1何晓群,刘文卿,应用回归分析,北京:中国人民大学出版社,2007.7。2段进宇, 杨佩昆,混合车辆的流通能力 J ,同济大学学报, 1995(1) : 37-42.3 杜强,贾丽艳,SPSS 统计分析从入门到精通M,北京:人民邮电出版社,2009,P257-
25、271,P503.4董文永,刘进,丁健利,朱福喜,最优化技术与数学建模,北京:清华大学出版社,2010.9。5姜启源,谢金星,叶俊,数学建模,北京:高等教育出版社,2011年。八、附录清单附录一:车辆折算系数附录二:视频2车辆数与横断面道路实际通行能力九、附录附录一:交通量调查车型划分及车辆折算系数车型换算系数小汽车1.0公交车2.0摩托车0.5大型货车2.0附录二: 视频2中不同时刻点各种车型的数量与横断面道路实际通行能力统计时刻实际收集数据间隔时间事故断面处大巴车数量/辆事故断面处小汽车/辆折算后的车辆数通行能力pcu/h17:35:173021114168017:35:473011213
26、.5162017:36:1730189.5114017:36:473008896017:37:173011112.5150017:37:473021013156017:38:1730189.5114017:38:473001212144017:39:17301910.5126017:39:47303610.5126017:40:173008896017:40:47302912144017:41:1730099108017:41:473001010120017:46:1730099108017:46:47302811132017:47:4730269108017:48:1730281113201
27、7:48:473001010120017:49:1730189.5114017:49:473001111132017:50:173011011.5138017:50:473001111132017:51:173001010120017:52:173001011.5138017:52:47301910.5126017:54:173001111132017:54:473011112.5150017:55:173001010120017:55:47302710120017:56:17301910.5126017:57:1730189.5114017:57:47301910.5126017:58:1730189.5114017:58:4730359.5114017:59:47302710120018:00:173001111132018:00:47301910.5126018:01:1730189.5114018:01:4730099108018:02:17302811132018:02:47302811132018:03:1730111121440平均:1273pcu/h
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1