43 公交车调度方案的优化设计 彭广政Word文档下载推荐.docx
- 文档编号:17148828
- 上传时间:2022-11-28
- 格式:DOCX
- 页数:21
- 大小:113.27KB
43 公交车调度方案的优化设计 彭广政Word文档下载推荐.docx
《43 公交车调度方案的优化设计 彭广政Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《43 公交车调度方案的优化设计 彭广政Word文档下载推荐.docx(21页珍藏版)》请在冰豆网上搜索。
ρ:
乘客的抱怨度;
N:
一共需要的车辆数;
S:
一天总的发车次数;
η:
平均每车次的载客率;
在第i时间段内上车的总人数;
一天的乘车总人数
四、模型分析与求解
4.1模型一
本题要求设计全天(工作日)的公交车调度方案,这里需要考虑乘客和公交公司两方面的利益,是一个多目标的优化问题。
其中可以供选择的目标函数主要有:
1.乘客候车时间要尽量短;
2.候车时间超过5分钟乘客数要尽量少;
3.公交公司所需的总车辆数尽量少;
4.全天范围内,发车的总次数尽量少;
5.平均每车次的载客率尽量高等等。
以上的目标可以用乘客利益和公司利益分为两类,这两类目标是相互冲突的,不可能同时达到最大。
工作日的早高峰正是多数乘客上班的时间,也是一天中乘坐公交车人数的高峰期,所以这段时间里所需的车辆数也是最多的。
从乘客的方面考虑,早上上班迟到对他的利益的损失相当大,因此乘客希望候车时间一般不要超过5分钟。
这时应以乘客的抱怨程度尽量小为主要目标,求得公交公司在早高峰期间的所需的最少车辆数。
在其余时间段里,乘客候车时间一般不要超过10分钟,这时考虑到公交公司的利益使其在这段时间内所发的总共发车的车次总数最少,以及提高每车次的载客率为主要目标。
因此我们首先确定出早高峰期,针对早高峰期的数据,在一定的乘客抱怨水平下,求出共需多少辆车,然后再根据全天其它时段的数据,并综合其它指标求出两个起点站的发车时刻表。
由于题目中所给出的仅是各站一个小时上下车人数的数据,对于我们的计算而言太过粗糙。
首先想到的是运用题中的数据对每一车站在各时段上车和下车的人数进行分布拟合,但这样做也有很大的缺点,因为各时段每个站点上下车人数受上下班时间以及道路沿线工厂等因素影响很大,从而导致各时段前后相关性很小。
而对各时段上车和下车的人数进行分布拟合就人为的增加了各时间段的上下班人数的相关性,与实际情况不符。
实际中如果把统计做的更细致或者知道那些影响上下车人数分布的因素,就可以较好的求出这些分布;
由于缺乏我们对这些情况得了解,所以我们假设各站的上下车人数在各个时间段(一小时)内分布是均匀的,即
E[numik(t,Δti)]=λ(Δti)=numik/Δti
其中numik(t,Δti)表示在[t,Δti]内上车的人数。
4.2模型二
为了更好的建立模型,首先要明确下面几个问题:
1.时间段的划分:
假设在题目中给出的各时段(一小时)内,各车站上下车的乘客人数分布均匀,这样就可将全天分为18个时间段,分别对每一个时间段进行考虑,并认为每个时间段内的发车间隔时间Δti上和Δti下分别为常数,但两者不一定相等。
2.对下行方向的处理:
从题中数据可以看出上行方向比下行方向多一个车站A1,我们对此的处理是在下行方向同样也补上一个车站A1,并且令这个车站在任何时段上车和下车的人数均为0。
3.对乘客平均抱怨度的定义:
考虑到一个人的抱怨程度是一个模糊的表述,它与候车时间的长短有关,候车时间越长,抱怨程度越大,但候车时间足够短时又不会抱怨。
经过分析可以定义第i个时间段上行的(或下行的)第j个乘客的抱怨度为:
或
式中
,
表示当此乘客不同等待时间wij对应不同的抱怨度。
可以看出抱怨度不仅与等待时间的长短有关,而且还会与所在的时间段i有关。
很明显,早高峰期间和平时时段里等待同样长的时间,前者给乘客造成的损失可能更大些,因此抱怨度也会相应大一些。
有了每个乘客抱怨度的定义,第i时间段的平均抱怨度为:
一天内的平均抱怨度为:
式中ωi表示第i段时间内区间的平均抱怨度对总平均抱怨度的权重,可取
。
其中wij是由Δti上或Δti下及numik(t,Δti)决定的,其中numik(t,Δti)是一个随机量,故wij也是一个随机量,从而一天内的平均抱怨度ρ也是一个随机量,可表示为f(Δti上,Δti下)。
4.总车辆数的确定:
一天所需的总车辆数N等于各时段所需的总车辆数Ni中的最大值,即N=g(Δti上,Δti下)=max{N1,N2,…,N18},而每一时段所需的总车辆数由上行车辆数、下行车辆数,加开车辆数三部分组成,有Ni=Ni上+Ni下+Ni加
其中
这里[·
]表示对括号内的数取整。
5.对平均每车次的载客率的定义:
考虑到每车次的运营成本基本不变,这样
平均每车次的载客数目的多少就能反映公司的利益。
于是我们定义平均每车次载客率定义为:
η=num总/S即:
代表第i时段在第j车站到上车人数(包括上行和下行);
代表第时段的时间间隔;
第i时段上行或下行的发车时间间隔。
平均每车次的载客率的高低直接反映了一个调度方案对于公交公司的收益率。
一般地乘坐公交车是按次计费的,所以总上车人数即反映了公交公司一天的收入,而总发车次数则反映了公交公司一天的支出。
6.据以上分析,我们建立如下模型:
目标:
minE[ρ]=E[f(Δti上,Δti下)]
minN=g(Δti上,Δti下)
min
这里[•]表示对括号内的数取整
maxη=num总/S
调度要求:
1.每辆车上承载的人数不超过120人;
2.在给定时间段Ti=60(分)内Δti上,Δti下为定值。
解法分析
在我们建立的模型中的多个目标中,总共需要的车辆数N涉及到公司建立一条公交线路的初始投资,每辆车所需的资金巨大,应被首先考虑。
而要确定总共需要的车辆数,只需求出早高峰期(我们根据题中给出的数据,假设早高峰期为7:
00—8:
00和8:
00—9:
00两个时段)内所需的车辆数即可。
考虑到实际求解过程中,对于前面模型中所定义的抱怨度在各个时间段内对于不同等待时间长度取值问题,可以通过实际的调查数据得到;
简化地想,如果对应所有的区间,顾客等待时间长度大于5分钟时都取1,而小于5时都取0,那么这是所定义的抱怨度直观意义就是指所有时间内等待时间超过5分钟的人数占总人数的比值,但显然着这种定义太粗糙;
由于缺乏实际的调查,我们在以下求解过程中对抱怨度在各个时间段内对于不同等待时间长度取值作以下假定:
早高峰期间γi1--γi4分别取[00.311.52.4];
而其它时间段内γi1--γi4分别取[00.150.50.751.2]。
解法一:
时间步长法
总体思路:
在给定的假设原则下,通过逐步改变发车时间间隔用计算机模拟各个时间段期间的系统运行状态,确定最优的发车时间间隔。
为简化计算,可以设定每个车站单位时间内上下车人数分别正比于该车站在这个时间段上下车的总人数,即numik(t,Δti)=λ(Δti)=numik/Δti。
运行步骤:
1.初始发车时间间隔
,i=1,2,…,18;
2.设置初始状态。
模拟时钟
、终点时间
3.设置每个车站等车人数以及待发车辆的最初状态,将每个车站的等车人按其等待时间的长短分为以下几类:
、
、…;
4.判断是否需要发车以及是否有车到站,并更改一次各车站等车人数及运行车辆的状态。
其中运行车辆的状态包括实际承载人数和空余座位数。
判断是否到了终点时间?
——若是,转第5步;
否则
转第2步;
5.统计整个过程中各区间等车时间超过5分钟以上的人数,对这些人加权求和,然后除以该过程总的上车人数得到平均抱怨程度。
判断是否大于给定的目标抱怨度?
——若大于给定值则转第1步,并改变发车时间间隔Δti;
否则,给出结果并转第6步;
6.结束。
解法二:
等效法
由于不同上车规则所对应的所有乘车人员的等待时间之和T总是相等的,可以把先到先乘车这个规则所对应的T总等效成后到先乘车规则的情况求出来。
由T总除以总乘车人数可以得到乘客的平均等待时间T平,对乘客而言,T平当然越小越好,即可以把T平作为一个目标;
当然也可以这样解释:
通过对乘客等待时间的统计,可得到等待时间位于各个时间区间内的乘客人数。
通过拟合有如下图所示的分布,图中虚线代表T平,那么图中阴影所示的面积就反映了等车时间超过5分钟的比重。
由图可看出,T平的值越小,抱怨度也会越小。
图一
等效法求解T平的计算原理及方法:
利用等效的概念求解是基于这样一种结论:
只要每个乘客到达的车站的时间和发车的间隔确定,那么先到先上车的规则和后到先上车的规则两种情况下所有人员的等车时间总和总是相等的。
例如:
甲在一个车站等待,过了一个周期来了一辆车和乙,但车只有一个座位;
又过一个周期,又来了一辆车和丙,也是只能一个人上。
那么,这段时间里如果按照先来先上规则,甲和乙的等车时间都是一个周期;
如果按照后来先上规则,则甲的等车时间是两个周期,而乙可以不等待,
但是两人的等待时间之和是一样的。
如果沒有更多的空位如果没有更多的空位,甲将被“滞留”在那里。
这种现象可解释为:
如果出现等待,先到先上规则是每个乘客都得等一段时间,而后到先上规则却是先到的人员一直在等,后来得人反倒可以即来即上。
相当于后来人员的等待时间被折合到先到人员身上的缘故。
先到的乘客担当了全体人员的“替罪羊”,从而形成了在整个时间隔内永远上不了车的滞留情况。
这样就可以只通过计算这些滞留人员的滞留时间之和,得到全体的平均等待时间。
定义第i个时段内,第j个车站的净上车人数aij为该车站在这一时段内的上车人数减去下去车人数。
用以上定义对数据进行预处理,即得到每个车站在各时间段内的的净上车人数aij。
考虑在Δti内,一辆本来已经满载的车经过第j个车站的情况:
如果这个车站的净上车人数大于0,则这个站滞留的乘客人数为Δti内的净上车人数aijΔti/Ti;
若净上车人数小于等于0,则这个车站可以为后面的车站提供空位置,从而使得下面站点的滞留人数减少,即可等效为它能提供的滞留人数为一负值,用来对消后面车站的滞留人数;
但是如果其后面所有车站的净上车人数之和为负的话,它所提供的空位永远也不会有人来坐,负的滞留人数就没有实际意义。
这时就可令它和它之后的所有车站所能提供的滞留人数都为0。
按照以上方法可求得每个时段上只有前面若干个车站会出现滞留人数,即可简化计算;
同时对第i时段的始发站而言,它所提供的滞留人数为:
图二给出的是7:
00段各站上行方向上的净上车人数)。
图二
按上述等效原理,算出各时间段内各站的净上车人数,即滞留人数,而滞留的人从到站起就一直等待到时间段结束,其等待时间成一个等差数列。
由此可得总的等待时间的计算公式为:
而平均等待时间Ti平=Ti总/
下面给出上行方向各段平均等待时间Ti平的目标,求得结果如下:
表一:
上行方向各段平均等待时间Ti平的目标
Ti平
7
3
2
4
Δti
9.7
2.5
1.46
2.58
4.5
5.6
4.8
5.5
6.1
6.9
3.2
2.65
11.3
12.3
15.3
37.3
误差讨论:
该方案计算时,当一个人的等待时间
时,他的等待时间就被忽略了,导致结果偏小,这里应对平均等待时间T平加上一个Δti/2的修正;
解法三:
等效时间步长法
时间步长法虽然将全天的数据作为一个整体来处理,充分考虑到各时间段的数据对相邻时间段的影响,可以较好的模拟出全天公交车的运营情况,并且给出对乘客抱怨度较为精确的描述;
但Δti是一个18维的向量,我们对于初值的确定缺乏依据,而导致大量盲目的搜索。
另一方面,等效法可方便快捷的给出在平均等待时间T平约束下各段的发车时间间隔;
但其根本缺陷在于只对每个时间段内的数据进行处理,而没有考虑到上一时间段遗留下来的人对本时间段的影响及本时间段遗留的人对下一时间段的影响。
因此算出的结果对于全天来说是就不见得特别好了。
结合前面两种算法,我们想到可以首先用等效法算出几组Δti的初值,然后将
这些初值带入时间步长法中进行计算,得出平均抱怨度最小的一组Δti作为我们的结果。
模型的结果:
按解法三求得结果如下:
1.两个起点站的发车时刻表如下:
表二:
公交车全天调度方案
发车时间间隔
上行
下行
5:
00-6:
00
6:
00-7:
7:
00-8:
9′30″
30′00″
2′30″
6′50″
1′30″
2′40″
8:
00-9:
9:
00-10:
2′20″
4′30″
4′20″
10:
00-11:
11:
00-12:
5′40″
4′50″
7′40″
8′50″
12:
00-13:
13:
00-14:
5′30″
6′10″
8′20″
14:
00-15:
7′20″
15:
00-16:
7′10″
5′50″
3′30″
16:
00-17:
3′10″
17:
00-18:
2′10″
18:
00-19:
7′00″
3′20″
9′20″
19:
00-20:
11′20″
20:
00-21:
12′20″
21:
00-22:
15′00″
22:
00-23:
20′00″
2.总共需要N=49辆车;
3.一天内总的发车次数为S=440次,其中下行方向为202次,上行方向为238次;
4.平均每次车载客率P=246人/车次
5.抱怨度ρ=0.2831
表三:
各等车时间区间等车人数
等待时间(分)
4~5
5~6
6~7
7~8
8~9
9~10
>
10
等待人数
8112
7894
5739
1557
1431
1789
3374
方案模拟:
不同工作日同一时段乘客到达总数基本相同,但由于种种随机因素的影响,它总是有一定的上下波动。
于是在单位时间内到达乘客数均值m上引入随机量
,它服从均值为0正态分布,其均方差与m成正比关系,即
~N(0,(
)2),于是实际到达人数
,其中可以通过调整
来控制波动的程度。
在
不同水平下的对上行方向调度方案仿真,同一水平下作100次仿真然后取均值得到一天不同候车时间w占全天乘车人数的比值及最后总的抱怨度ρ。
表四:
不同随机水平
上行方向调度方案仿真结果
5
6
8
9
ρ
=0
0.0688
0.05
0.0137
0.0127
0.0158
0.0297
0.2826
=0.01
0.0128
0.0156
0.2827
=0.02
0.0501
0.014
0.0129
0.2830
=0.05
0.0687
0.0139
0.0151
0.0302
0.2835
=0.1
0.0508
0.0147
0.0295
0.2841
可以看出,调度方案对
即对于数据的波动不敏感。
最大随机水平
=0.1时,抱怨度相对变化只有0.53%,此方案有较大的适用范围。
五、模型评价与改进
5.1模型的优点
1.针对题目中给出的数据,充分考虑到各方面的利益建立了多个目标下的优化模型;
2.对模型对求解结合了时间步长法和等效法的优点,求解简单,而得到的方案较为理想;
3.对人到达各车站的时间的随机变化进行了模拟,检验了方案对扰动的敏感性;
4.易操作性,一方面公交公司的时刻表比较合理可行,另一方面驾驶员能容易记住自己的上班时间,以避免时间表混乱而引起误车现象。
5.2模型的缺点
1.对于两个起点站发车次数的均衡性,对调度方案的影响未做很好的分析,只进行了定性的分析;
2.用光滑曲线拟合的方法无法模拟真实的客流量曲线。
5.3模型的改进
通过对上面所得到的结果进行分析可以看出高峰期及高峰期前各时间段等车时间较长的人数较多,而高峰期后各时间段等车时间较长的人几乎没有。
其原因是我们的模型是在第i段的起点时刻才开始启用第i段的时间间隔
,而在第i段发出第一辆车的时刻,各站的上下车人数均已按第i段的数据开始出现了,但已在路上运行的车还是按第i-1段的时间间隔
发出的,这样车就明显的与乘客的需求不相符合。
当第i段的上下车人数大于第i-1段的上下车人数时,乘客的等待时间显然会变长,增加乘客的抱怨;
我们的结果中高峰期及高峰期前各时间段等车时间较长的人数较多就是这种情况;
反之,乘客的等待时间显然会变短,同时也使车辆的满载率过低,不符合公司的利益。
对此我们可对我们的模型作如下调整:
在第i段的起始时刻之前Tr时刻开始的Ti=60(分)时段内按第i段的时间间隔Δti发车,针对不同的Tr值进行搜索,得出使得平均抱怨度ρ最小的Tr。
参考文献:
[1]吴建国.数学建模案例精编[M].北京:
中国水利水电出版社,2005
[2]杨启帆,何勇,谈之奕.数学建模竞赛——浙江大学学生获奖论文点评[M].杭州:
浙江大学出版社,2006
[3]朱道元.数学建模案例精选[M].北京:
科学出版社,2003
[4]283840255.公交车调度方案的优化设计.XX文库(
附录(可以含程序、图、表格、证明过程)
某路公交汽车各时组每站上下车人数统计表上行方向:
A13开往A0
站名
A13
A12
A11
A10
A9
A8
A7
A6
A5
A4
A3
A2
A1
A0
站间距(公里)
1.6
0.5
1
0.73
2.04
1.26
2.29
1.2
0.4
1.03
0.53
5:
00-6:
上
371
60
52
43
76
90
48
83
85
26
45
11
下
13
20
81
32
18
24
25
57
6:
00-7:
1990
376
333
256
589
594
315
622
510
176
308
307
68
99
105
164
239
588
542
800
407
208
300
288
921
615
7:
00-8:
3626
634
528
447
948
868
523
958
904
259
465
454
205
227
272
461
1058
1097
1793
801
469
560
636
1871
1459
8:
00-9:
2064
322
305
235
477
549
271
486
439
157
275
234
106
123
169
621
971
440
245
339
408
1132
759
9:
00-10
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 43 公交车调度方案的优化设计 彭广政 公交车 调度 方案 优化 设计