数学建模深圳杯分析
发布时间:2018-10-13 14:04:16
发布时间:2018-10-13 14:04:16
“禁摩限电”政策效果综合分析
摘要
首先,对于交通资源总量来说,即道路通行能力,是指在特定的交通条件、道路条件及人为度量标准下单位时间通过的最大交通量。根据道路方式不同,分为机动车道、非机动车道、人行道的通行能力,而机动车道又根据道路等级划分为高速公路、快速路、主干道、支路等等。我们在下面的文章中主要讨论深圳机动车道的通行能力,大致可分为无平面交叉口的路段通行能力以及交叉口的通行能力。通过对不同路况的交通通行能力的研究,可以得出部分仿真条件。
对于交通需求结构,我们可以根据虚拟网络结构,建立超量需求模型,得出部分仿真数据。结合对交通资源总量的研究,我们可以利用MATLAB软件对某一路段进行交通状况的模拟仿真,从而得出不同情况下道路通行状况的变化情况。结合道路交通总量和交通需求结构这两个条件,我们得出结论:随着车辆驶入概率的不断增加,车辆自由通行的概率逐渐下降,车辆拥堵的概率明显上升。这就说明随着现实中车辆总数的日益上升,很可能导致城市道路无法承受现有的交通总量,出现普遍的交通拥堵状况。所以,禁摩限电政策有助减少交通总量,从而在不改变总体道路承载能力的情况下缓解交通拥堵问题。
对于交通工具效率,考虑到效率受交通工具速度、总里程数、交通工具年龄、天气、 路况等许多因素影响,所以可以采取多个指标衡量深圳的各种交通工具效率,由此我们可以建立基于多指标的综合评价模型。利用MATLAB软件进行层次分析后可得各个交通工具的效率影响。
对于交通安全性,与对于交通工具效率的建模类似,我们需要对各种交通工具的安全性进行评价,考虑到安全性受总里程数、交通工具年龄、乘客配合程度等许多因素影响,故可以采取多指标衡量深圳的各种交通工具安全性,由此我们可以建立基于多指标的综合评价模型。同样利用MATLAB软件进行层次分析后可得各个交通工具的效率影响。
对于环境影响因素,可采用碳排放值来衡量各种交通工具对安全和环境的影响。利用MATLAB软件对获得的数据进行分析后,得到公共交通在综合评分上更占 优势,从而我们得出结论:禁摩限电政策有助于减少非公共交通类的交通工具,从而促进公共交通的发展,从而保证城市各方面的发展。这也说明了禁摩限电政策的正确性。
关键字:道路交通 Matlab 层次分析 元胞机仿真
一、问题的重述
城市交通系统是城市赖以生存和发展的保证,交通的顺畅程度直接影响着城市的发展。近年来,随着城镇化进程的不断加快和汽车工业的快速发展,近年来我国城市机动车拥有量得以大幅度增加。所以城市道路交通资源是有限与交通需求的高速增长产生了越来越大的矛盾。在短时间内无法产生交通革命的情况下,限制部分交通工具成为缓解城市出行压力的主要措施。不少城市采取的限牌、限号措施取得了一定成效,然而也有不少人无法理解。现在,我们针对“禁摩限电”政策,对其进行科学的、不带意识形态的论证,并对其评价。
二、问题分析
为了使“禁摩限电”这一政策得到大多数人的支持,我们可以从深圳的交通资源总量(即道路通行能力)、交通需求结构、各种交通工具的效率及对安全和环境的影响等因素和指标出发,建立数学模型对其进行定量分析。
1.针对交通资源总量,道路通行能力是指在特定的交通条件、道路条件及人为度量标准下单位时间通过的最大交通量。根据道路方式不同,分为机动车道、非机动车道、人行道的通行能力,而机动车道又根据道路等级划分为高速公路、快速路、主干道、支路等等。下文主要讨论深圳机动车道的通行能力,大致可分为无平面交叉口的路段通行能力以及交叉口的通行能力。通过对不同路况的交通通行能力的研究,可以得出部分仿真条件。
2.针对交通需求结构,我们可以根据虚拟网络结构,建立超量需求模型,得出部分仿真数据。结合对交通资源总量的研究,我们可以利用matlab软件对某一路段进行交通状况的模拟仿真,从而得出不同情况下道路通行状况的变化情况。
3.针对交通工具效率,考虑到效率受交通工具速度、总里程数、交通工具年龄、天气、路况等许多因素影响,故可以采取多指标衡量深圳的各种交通工具效率,由此我们建立基于多指标的综合评价模型。利用matlab软件进行层次分析后可得各个交通工具的效率影响。
4.针对交通安全性,与对于交通工具效率的建模类似,我们需要对各种交通工具的安全性进行评价,考虑到安全性受总里程数、交通工具年龄、乘客配合程度等许多因素影响,故可以采取多指标衡量深圳的各种交通工具安全性,由此我们建立基于多指标的综合评价模型。同样利用matlab软件进行层次分析后可得各个交通工具的效率影响。
5.针对环境影响因素,可采用碳排放值来衡量各种交通工具对安全和环境的影响。
3、模型假设
1.假设查询所获得的数据真实可靠。
2.不考虑matlab软件内部的计算误差。
3.不考略matlab软件计算量大小对结果的影响。
4.假设题目各种类车辆各方面数据基本相同,没有明显的相关性。
5.不考虑不同价位车辆间存在的差异,只考虑车辆类型上的不同。
6.忽略模型计算过程中存在的微小误差。
7.假设不考虑其他无关因素对模型及结果的影响。
四、符号说明
五、模型的建立与求解
为了使“禁摩限电”这一政策得到大多数人的支持,我们从深圳的交通资源总量(即道路通行能力)、交通需求结构、各种交通工具的效率及对安全和环境的影响等因素和指标出发,建立数学模型对其进行定量分析。
1.模型一
5.1.1交通资源总量
道路通行能力是指在特定的交通条件、道路条件及人为度量标准下单位时间通过的最大交通量。根据道路方式不同,分为机动车道、非机动车道、人行道的通行能力,而机动车道又根据道路等级划分为高速公路、快速路、主干道、支路等等。下文主要讨论深圳机动车道的通行能力,大致可分为无平面交叉口的路段通行能力以及交叉口的通行能力。
无平面交叉口的路段通行能力
基本通行能力
基本通行能力是指道路与交通处于理想情况下,每一条车道(或每一条道路)在单位时间内能够通过的最大交通量。
基本通行能力:
word/media/image6_1.png
其中,word/media/image7_1.png为行车速度(word/media/image8_1.png),word/media/image9_1.png为车头最小时距(word/media/image10_1.png),word/media/image11_1.png为车头最小间隔(word/media/image12_1.png)。
可能通行能力
可能进行能力是指在通常的道路交通条件下,单位时间内通过道路一条车道或某断面的最大可能车辆数。
可能通行能力:
word/media/image13_1.png
其中,word/media/image14_1.png为连续小客车车流平均车头时距(word/media/image15_1.png)。
设计通行能力
设计通行能力是指道路交通的运行状态保持在某一设计的服务水平时,道路上某一路段的通行能力。
单条车道设计通行能力:
其中,word/media/image4_1.png表示道路分类系数,不同道路的分类系数见下表1:
表1:不同道路的分类系数
多车道的通行能力还要考虑变换车道的影响,因此公式为:
其中word/media/image5_1.png为机动车道单向通行能力折减系数。不同车道数目的折减系数见下表:
表2:不同车道数目的折减系数
有平面交叉口的路段通行能力
此处通行能力是指设计通行能力。
十字交叉路口通行能力计算方法如下。当出现可插间隙时间word/media/image18_1.png时,次要方向的车流可以相继通过的随车时距为word/media/image19_1.png,推导出计算公式:
。
其中,word/media/image21_1.png为非优先通行次干道上可以通过的交通量(word/media/image22_1.png),word/media/image23_1.png为主干道优先通行的双向交通量(word/media/image22_1.png),word/media/image24_1.png,word/media/image25_1.png为临界间隙时间word/media/image26_1.png(6~8s或7~8s),word/media/image27_1.png为次干道上车辆间的最小车头距(3s或5s)。
环形交叉口的通行能力计算方法如下。在常规环交的通行能力计算中较著名和使用较广泛的公式为沃尔卓普公式:
其中,word/media/image29_1.png为交织段上最大通行能力(word/media/image22_1.png),word/media/image30_1.png为交织段长度(word/media/image31_1.png),word/media/image32_1.png为交织端宽度(word/media/image31_1.png),word/media/image33_1.png为环形道路口引道的平均宽度(word/media/image31_1.png),word/media/image34_1.png为交织段内交织车辆与全部车辆之比(%)。
5.1.2交通需求结构
虚拟网络的构建
图1 a)中的网络代表某一城市道路网系统,其着重表示了某一word/media/image35_1.png对word/media/image36_1.png。定义一条虚拟的附加路段,该路段的位置为链接每一word/media/image35_1.png对word/media/image36_1.png的起终点,附加路段的流量定义为word/media/image37_1.png。
word/media/image38_1.png
a)原路网 b)增设附加路段
图1:附加虚拟网络图
假设word/media/image35_1.png对word/media/image36_1.png的交通需求函数为word/media/image39_1.png,需求函数是有上限的,假设起点word/media/image40_1.png到终点word/media/image41_1.png的上限是word/media/image42_1.png,则word/media/image43_1.png,它表示起点word/media/image44_1.png到终点word/media/image45_1.png间不能容纳的出行量,即超量需求量。令表示逆需求的超量变量函数,即word/media/image47_1.png,word/media/image48_1.png,则word/media/image49_1.png具有路段特性函数的所有性质,即为该虚拟附加路段的阻抗函数。
超量需求均衡模型
基于Wardrop均衡原理以及构建的附加虚拟网络,建立超量需求均衡模型的数学形式描述如下:
其中,。模型中所用变量和参数的含义:word/media/image3_1.png为路段word/media/image52_1.png上的交通流量,它们组成的向量为word/media/image53_1.png;word/media/image54_1.png为路段word/media/image55_1.png的交通阻抗;word/media/image56_1.png为路段word/media/image55_1.png的以流量为自变量的阻抗函数;为交通需求函数的反函数;word/media/image58_1.png为起点word/media/image44_1.png到终点word/media/image45_1.png总出行量的上限;word/media/image59_1.png为附加路段的交通流量;word/media/image60_1.png为点对word/media/image61_1.png间的第word/media/image62_1.png条路径的交通流量,其向量word/media/image63_1.png;word/media/image64_1.png为关系变量,如果路段word/media/image65_1.png在word/media/image61_1.png间的第word/media/image62_1.png条路径上,,否则,。
利用matlab软件模拟一段路的通行情况后,绘制图像如下:
word/media/image68_1.png
随着车辆驶入该路的增加,车辆自由通行概率显著下降。
word/media/image69_1.png
随着车辆驶入该路的增加,车辆拥堵概率显著上升。
通过分析Matlab绘制的图像,可知,在总体交通承受能力不变的情况下,随着车辆驶入概率的不断增加,车辆自由通行概率迅速下降,车辆拥堵概率明显上升。这说明随着现实中车辆总数的日益上升,很可能导致城市道路无法承受现有交通总量,出现普遍的交通拥堵状况。“禁摩限电”政策有助减少交通总量,从而在不改变总体道路承载能力的情况下缓解交通拥堵问题。
5.2模型二的建立与求解
5.2.1各种交通工具的效率
我们需要对各种交通工具的效率进行评价,考虑到效率受交通工具速度、总里程数、交通工具年龄、天气、路况等许多因素影响,故可以采取多指标衡量深圳的各种交通工具效率,由此我们建立基于多指标的综合评价模型。
首先进行数据的无量纲化处理,使各项指标具有可比性,这是一种通过数学手段来消除原始变量量纲影响的方法,本文采用极差化法对指标值进行无量纲化处理。
假设一个综合评价问题中共有word/media/image70_1.png个被评价对象,分别记为:
word/media/image71_1.png
系统有word/media/image72_1.png个评价指标,分别记为:
word/media/image73_1.png
word/media/image74_1.png表示第word/media/image75_1.png个评价对象中第word/media/image76_1.png个指标对应的指标值,令:
word/media/image77_1.png
其中:
word/media/image78_1.png
则word/media/image79_1.png是无量纲的标准观测值。
对于各个交通工具效率的定量综合评价,可以通过层次分析法来建立数学模型,即给出各个交通工具(地铁、公共汽车、小轿车、摩托车、电动三轮车、电动自行车)效率的排名,这是层次分析法模型的目标层;准则层为影响交通工具效率的因素,即速度、总里程数、交通工具年龄、天气、路况;该模型方案层为选择的行驶路线。
建立层次分析模型后,我们就可以在各层元素中两两比较,构造出比较判断矩阵。首先,考虑准则层的各因素对目标层的影响情形,用word/media/image80_1.png分别表示下面五个因素:速度、总里程数、交通工具年龄、天气、路况;用word/media/image81_1.png分别表示地铁、公共汽车、小轿车、摩托车、电动三轮车、电动自行车这六个交通工具,这里我们假设word/media/image80_1.png对各交通工具效率的影响是相同的,则设word/media/image82_1.png与word/media/image83_1.png对目标层的影响之比为word/media/image84_1.png(word/media/image85_1.png),则准则层对目标层的判断矩阵为:
word/media/image86_1.png
由判断矩阵word/media/image87_1.png,可以求出特征值和特征向量,通过计算可得word/media/image87_1.png的权向量为:
下面考虑方案层对准则层每个因素word/media/image89_1.png的影响权重,前面已经将数据进行了标准化处理,因此我们在构造第3层对第2层每一个因素word/media/image90_1.png的判断矩阵时,只要利用城市word/media/image91_1.png的因素word/media/image92_1.png的指标word/media/image93_1.png与城市word/media/image94_1.png的因素word/media/image92_1.png的指标word/media/image95_1.png的比值即可,可以得到判断矩阵:
word/media/image96_1.png
5.2.2各种交通工具的安全性
与对于交通工具效率的建模类似,我们需要对各种交通工具的安全性进行评价,考虑到安全性受总里程数、交通工具年龄、乘客配合程度等许多因素影响,故可以采取多指标衡量深圳的各种交通工具安全性,由此我们建立基于多指标的综合评价模型。
首先进行数据的无量纲化处理,使各项指标具有可比性,这是一种通过数学手段来消除原始变量量纲影响的方法,本文采用极差化法对指标值进行无量纲化处理。
假设一个综合评价问题中共有word/media/image70_1.png个被评价对象,分别记为:
word/media/image71_1.png
系统有word/media/image72_1.png个评价指标,分别记为:
word/media/image73_1.png
word/media/image74_1.png表示第word/media/image75_1.png个评价对象中第word/media/image76_1.png个指标对应的指标值,令:
word/media/image77_1.png
其中:
则是无量纲的标准观测值。
对于各个交通工具安全性的定量综合评价,可以通过层次分析法来建立数学模型,即给出各个交通工具(地铁、公共汽车、小轿车、摩托车、电动三轮车、电动自行车)安全性的排名,这是层次分析法模型的目标层;准则层为影响交通工具安全性的因素,即:行驶总里程数、交通工具年龄、乘客配合程度、驾驶员能力、旅客客运量;该模型的方案层为各个交通工具选择的行驶路线。
建立层次分析模型后,我们就可以在各层元素中两两比较,构造出比较判断矩阵。首先,考虑准则层的各因素对目标层的影响情形,用word/media/image80_1.png分别表示下面五个因素:行驶总里程数、交通工具年龄、乘客配合程度、驾驶员能力、旅客客运量;用word/media/image81_1.png分别表示地铁、公共汽车、小轿车、摩托车、电动三轮车、电动自行车这六个交通工具,这里我们假设word/media/image80_1.png对各交通工具安全性的影响是相同的,则设word/media/image82_1.png与word/media/image83_1.png对目标层的影响之比为word/media/image84_1.png(word/media/image85_1.png),则准则层对目标层的判断矩阵为:
word/media/image86_1.png
由判断矩阵word/media/image87_1.png,可以求出特征值和特征向量,通过计算可得word/media/image87_1.png的权向量为:
下面考虑方案层对准则层每个因素word/media/image89_1.png的影响权重,前面已经将数据进行了标准化处理,因此我们在构造第3层对第2层每一个因素word/media/image90_1.png的判断矩阵时,只要利用城市word/media/image91_1.png的因素的指标与城市的因素的指标的比值即可,可以得到判断矩阵:
5.2.3各种交通工具对环境的影响
可采用碳排放值来衡量各种交通工具对安全和环境的影响。目前,计算碳排放值一般采用能耗折算的方法,即不同的能源使用具有不同的碳排放系数word/media/image106_1.png,用公式表述为:
交通排放量=
其中,word/media/image109_1.png为能源使用量值。根据国际标准,常见能源使用量与word/media/image110_1.png排放量之间折算系数见下表。
常见能源使用量与word/media/image110_1.png排放量之间折算系数
我们平时计算时,通常会从家庭的角度来考虑,一般采用“碳足迹”的计算公式。“碳足迹”来源于一个英语单词“Carbon Footprint”,是指一个人或团体的“碳耗用量”。“碳”就是煤炭、石油、木材等由碳元素构成的自然资源,碳消耗的最多,制造二氧化碳就越多,碳足迹就越大。因此,所谓的碳足迹就是指一个人或组织团体的能源消耗行为对自然产生的影响。“碳足迹”的计算如下表:
家庭word/media/image110_1.png排放量计算表
利用mtalab软件进行数据分析后可以得出,公共交通在交通效率方面略低于其他交通工具,但是公共交通在对环境的污染方面明显低于其他交通工具,且公共交通的安全性也明显优于其他佳通工具。经过交通效率,安全性和环境污染三方面的综合分析,得出公共交通的推行对各方面的发展更为有利。
禁摩限电政策有助于减少非公共交通类的交通工具,从而促进公共交通的发展,从而保证城市各方面的发展。这也说明了禁摩限电政策的正确性。
六、参考文献
[1]米粮川; 杨洪澜; 王世刚,Matlab基础的教学思想,高师理科学刊。
[2]邵伟,蒙特卡洛方法及在一些统计模型中的应用。
[3]陆东鑫,计算机工程与应用,浙江大学,2011。
七、模型评价
7.1模型的优点
此模型对问题进行系统性分析,层次分明,定性定量综合分析了各因素“禁摩限电”的影响,并且所需要的定量数据信息较少,所得结果简单明确,容易为决策者了解和掌握。
7.2 模型的缺点
指标过多时数据统计量大,且权重难以确定,特征值和特征向量的精确求法比较复杂,不容易发现指标的相对重要性的取值里到底是哪个有问题,哪个没问题。
八、附录
function [ v d p ] = multi_driveway( nl,nc,fp,dt,nt )
% 在某一特定车流密度下的(车流密度由fp决定)单、双车道仿真模型
% nc:车道数目(1或2),nl:车道长度——输入参数
% v:平均速度,d:换道次数(1000次)p:车流密度——输出参数
% dt:仿真步长时间,nt:仿真步长数目——输入参数
% fp:车道入口处新进入车辆的概率——输入参数
% test:
% nl = 400;fp = 0.5;
% nc = 2;dt=0.01;nt=500;
%构造元胞矩阵
B=ones(2*nc+1,nl+2);
%奇数行为不可行车道
B(1:2:(2*nc+1),:)=1.2;
%初始化仿真元胞状态(1为无车,0为有车)
bb=B(2:2:2*nc,:);bb(bb~=0)=1;B(2:2:2*nc,:)=bb;B(2:2:2*nc,end)=0;
%显示初始交通流图
figure(1);
H=imshow(B,[]);
set(gcf,'position',[241 132 560 420]) ;%241 132 560 420
set(gcf,'doublebuffer','on'); %241
title('cellular-automation to traffic modeling','color','b');
%初始化化存储元胞上车辆状态的矩阵
S(1:nc,nl) = 0;
Q(1:nc,1:2) = 0;
Acc(1:nc,1:(nl+2))=0;
%初始化换道频率、平均速度、车流密度相关变量
ad = 0;
av(1:nt) = 0;
ap(1:nt) = 0;
c = 1;
for n = 1:nt
A=B(2:2:2*nc,:);
%确定前n-2个车辆的状态
S(:,:) = 0;
S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==1)=2;%加速的车
S(A(:,1:end-2)==0&A(:,2:end-1)==0)=3;%停车的车
S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==0)=1;%减速行驶的车
%确定最后2两个元胞的状态
Q(:,:) = 0;
Q(A(:,end-1)==0&A(:,end)==0) = 1;
Q(A(:,end-1)==0&A(:,end)==1) = 2;
Q(A(:,end-1)==1&A(:,end)==0) = 2;
Q(:,end) = 1;
%获得所有元胞上车辆的状态
Acc = [ S Q ];
%换路规则
if(nc>1&&n>nl/2)
%遍历每一个元胞
for g = 1:length(Acc(1,:))
%停车状态车辆如另一条路有2空位则换路
if( Acc(1,g)==3&&Acc(2,g)==0&&Acc(2,g+1)==0)
A(1,g)=1;
A(2,g)=0;
ad=ad+1;
elseif( Acc(2,g)==3&&Acc(1,g)==0&&Acc(1,g+1)==0 )
A(1,g)=0;
A(2,g)=1;
ad=ad+1;
%均速行驶车辆如另一条路有3空位则换路
elseif( Acc(1,g)==1&&Acc(2,g)==0&&Acc(2,g+1)==0&&Acc(2,g+1)==0 )
A(1,g)=1;
A(2,g)=0;
ad =ad+1;
elseif( Acc(2,g)==1&&Acc(1,g)==0&&Acc(1,g+1)==0&&Acc(1,g+1)==0 )
A(1,g)=0;
A(2,g)=1;
ad=ad+1;
end
end
%换路后重新设置元胞上的车辆状态
S(:,1:end) = 0;
S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==1)=2;%寻找加速的车
S(A(:,1:end-2)==0&A(:,2:end-1)==0)=3;%寻找停车的车
S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==0)=1;%寻找减速行驶的车
%确定最后2两个元胞的状态
Q(:,1:end) = 0;
Q(A(:,end-1)==0&A(:,end)==0) = 1;%
Q(A(:,end-1)==0&A(:,end)==1) = 2;
Q(A(:,end-1)==1&A(:,end)==0) = 2;
Q(:,end) = 1;
%获得所有元胞状态
Acc = [ S Q ];
end
%根据当前状态改变元胞位置
%匀速运行车辆向前走1格
A( Acc(:,1:end)==1 ) = 1;
A( [ zeros(nc,1) Acc(:,1:end-1)]==1 ) = 0;
%高速运行车辆向前走2格
A( Acc(:,1:end)==2) = 1;
A( [ zeros(nc,2) Acc(:,1:end-2)]==2) = 0;
%计算平均速度、换道频率、车流密度等参数
%获得运行中的车辆数目N
matN = A<1;
N = sum(sum(matN));
%获得运行中的车辆速度之和V
E = S((S==1)|(S==2));
V = sum(E);
%计算此时刻的车流密度并保存
ap(n) = N/( nc*(nl+2) );
%计算此时刻的平均速率并保存
if(N~=0&&n>nl/2)
av(c) = V/N;
c = c+1;
end
%在车道入口处随机引入新的车辆
A = [ round(fp*rand(nc,1))&A(1:nc,1) A(:,2:end)];
A(A~=0)=1;
%将新的车辆加入元胞矩阵中
B(2:2:2*nc,:)=A;
%显示交通流图
set(H,'CData',B);
%仿真步长
pause(dt);
end
%仿真结束,计算结果
d = ad;
p = mean(ap);
v = sum(av)/c;
end
%初始化仿真参量
v = 0;p=0;d=0;
nl = 100;
nc = 2;dt=0.01;nt=1000;
n=1;
fp = 1.2;
[ v d p] = multi_driveway( nl,nc,fp,dt,nt );
% 车流密度不变下的多车道仿真(包括单车道)
% nc:车道数目(1或2),nl:车道长度
% v:平均速度,d:换道次数(1000次)p:车流密度
% dt:仿真步长时间,nt:仿真步长数目
% fp:车道入口处新进入车辆的概率
v = 0;p=0;d=0;
nl = 100;nc = 1;
dt=0.01;nt=1000;
fp = 0.5;
[ v d p ] = multi_driveway( nl,nc,fp,dt,nt );
% 车流密度变化时的单车道仿真程序
% nc:车道数目(1或2),nl:车道长度
% v:平均速度,d:换道次数(1000次)p:车流密度
% dt:仿真步长时间,nt:仿真步长数目
% fp:车道入口处新进入车辆的概率
v = 0;p=0;d=0;
nl = 100;
nc = 1;dt=0.01;nt=1000;
n=1;
for fp = 2.5:-0.25:0.5
[ v d p ] = multi_driveway( nl,nc,fp,dt,nt );
va(n) = v;
pa(n) = p;
da(n) = d;
n=n+1;
%绘制平均速率-车流密度(v-p)曲线
figure(2)
plot(pa,va);
%绘制车流密度的变化曲线
figure(3)
plot(pa);
end