公园内电动汽车配置服务中心位置规模的设计方案及其分析.docx
《公园内电动汽车配置服务中心位置规模的设计方案及其分析.docx》由会员分享，可在线阅读，更多相关《公园内电动汽车配置服务中心位置规模的设计方案及其分析.docx（26页珍藏版）》请在冰豆网上搜索。
公园内电动汽车配置服务中心位置规模的设计方案及其分析
公园内电动汽车配置、服务中心位置规模的设计方案及其分析
摘要
本文利用多动态优化的方法给出了计算最大旅客数的方法。
首先通过数据分析，考虑到了方案的可操作性。
通过建立数学模型求出提高车辆利用率的方案。
根据游玩旅客人数的分布情况，设计出每个综合服务中心的位置及规模。
针对问题1，以实际对两辆车发车顺序进行假设（即两辆同时发，先忙后快，先快后慢），再通过列出不同的发车间隔，通过excel表格进行分析分，并用matlab画出拟合图像，通过图像分析得出结论，每个服务中心辆车同时发车，并间隔13min为最佳。
则在一个服务中心的车辆配置为A车33辆，B车辆。
全程之中，A车825辆，B车925辆。
旅客总数：
。
针对问题2，提高每辆车的利用率即提高每辆车的满载率，通过动态规划再次控制发车间隔，建立下一辆车到来后关于时间的效用函数，以此求得解，
针对问题3，因为考虑到每天进入公园游览的旅客数成poisson分布，所以为了减少服务中心资源不必要的浪费，根据旅客数的分布情况来设置服务中心的位置和规模。
在人流密集区域扩大规模，增加容纳量，反之，在人流较为稀疏的区域适当缩小规模及容纳量。
关键字：
汽车配置时间间隔载客量利用率poisson分布动态规划
问题提出
大型的森林公园要对公众开放，为了保护森林生态和防止旅行者受到猛兽等的侵害，特别设置了一条半封闭的旅游线路（线路在森林公园内部穿梭而进，最后回到出发点）。
该线路全长250公里，游览公园必须乘坐带有防护装置的两种电动汽车，时速分别为5公里每小时和10公里每小时，座位数分别为20位和8位。
沿线每隔10公里设有综合服务中心，提供临时的休息和晚间的住宿（每处容纳人数的上限是1000）。
电动汽车的开行时间为早上8点到下午5点之间，旅客可中途停靠服务中心休息，但是不允许在综合服务中心之间的线路上下车。
为了使旅行者更好地享受到旅行的乐趣，应尽量减少电动汽车相遇的次数和保持一定的车距。
每个旅行者因各自身体情况和喜好的不同，每天行车的时间约为3～6个小时。
1.旅游的旺季，在整条线路上如何配置两种电动汽车，使森林公园每天可以接待的旅客数达到最大？
2.请设计方案，在保证能达到每天旅客数为最大容量的85%这个前提下，尽量提高每部车的利用率。
3.如果由你来设置综合服务中心的位置和规模（容纳人数，最大不超过1500），如何配置两种车辆，使得森林公园可以接待较多的旅客数，以及使得每部车有较高的利用率?
模型假设
1.假设每辆电车汽车上下车人数是等比例的。
2.乘客上下车时间忽略不计。
3.电动汽车速度恒定，且无特殊事件发生（如抛锚、乘客由于特殊原因中途下车）。
4.以分钟作为最小的时间单位，这对安排时刻表是合理的。
5.假设电动汽车在行驶中都是单向的。
6.速度为5km/h，,座位数为20的车设为A车，速度为10km/h，座位数8的车设为B车。
字符说明
1.
表示公园内每个10公里设置的24个服务中心
2.
3.
4.
5.
6.
7.
问题分析
公园每天对旅客的接待数即为每天滞留在公园内部的旅客总数。
需考虑24个服务中心及入口的游客数。
电动汽车在每个服务中心都有发车。
为减少相遇次数同时保证游客数达到最大，车与车之间的发车时间是有时间间隔的，我们应用excel和matlab程序，得出最优时间间隔。
且按车的种类批次发车。
为提高利用率，通过动态规划建立效用函数，用lingo求得最优解。
由于客流的分布是呈poisson分布的，根据客流的多少，建立服务中心的规模。
以尽量减少资源浪费或不足的状况发生。
模型建立
问题一：
两种电动汽车的配置
由于现在处于旅游旺季，故每个服务中心都是有足够旅客的。
我们假设每个服务中心都是相互独立的且车在行驶过程中满载。
故可先着眼于一个服务中心进行考虑。
单独针对于一个服务中心来说，假设从此站经过的车的上车人数等于下车人数，则相当于本服务中心的人都可看作是乘坐本服务中心出发的车继续旅行的，故本服务中心的一天发车所能承载的总人数应小于服务中心的人数上限（即1000人）。
考虑到每个乘客的个人喜好不同，乘坐两种车的概率是相同的，所以我们让两种车间隔出发或同时出发以满足乘客需求。
同时以分钟为单位划分时间间隔，通过excel罗列出不同时间间隔两种车的发车数量、载客量及相遇次数。
通过excel数据和matlab图像比较，找出载客量最为接近1000且满足相遇次数较少的时间间隔。
由此得知，最佳的时间间隔为13min/辆。
对于相遇次数问题，见以下四个表格和matlab图像
当两车同时发车时，可分别列出在不同时间间隔两车的发车量和相遇次数
AB同时发车
间隔时间t
A车
B车
相遇次数
10
43
49
215
15
29
33
93
20
22
25
46
30
15
17
16
40
11
13
11
50
9
10
10
60
8
9
0
表1
画出表格为
分别于总发车数和相遇次数的关系
图1
当先发B车后发A车，两车间隔发车时，可分别列出在不同时间间隔两车的发车量和相遇次数
间隔时间
A车
B车
相遇次数
A先B后
10
22
24
66
15
15
16
31
20
11
12
11
30
8
8
8
40
6
6
6
50
5
5
5
60
4
4
0
表2
画出表格为
分别于总发车数和相遇次数的关系
图2
当先发A车后发B车，两车间隔发车时，可分别列出在不同时间间隔两车的发车量和相遇次数
间隔时间
A车
B车
相遇次数
5
43
49
266
B先A后
10
21
25
63
15
14
17
28
20
11
13
11
30
7
9
6
40
5
7
4
50
4
5
3
60
4
5
0
表3
画出表格为
分别于总发车数和相遇次数的关系
图3
由以上三个图像交点位置可知，当两车同时发车时可使相遇次数尽量减少而发车总数尽量大，同时可以以根据图1中，交点以下部分斜率可知间隔时间最值在10到15分钟内，由此我们列出下表
间隔时间
A车
B车
相遇次数
10
43
49
215
11
39
45
195
12
36
41
144
13
33
37
132
14
31
35
124
15
29
33
93
表4
对于载客量的问题，由于数据庞大，详见附录表一
由表4可知，在时间间隔为13min时，载客数达到最大。
A车可载旅客660人，B车可载旅客296人。
则A车的车辆数为：
；B车的车辆数为：
。
因为每个服务中心看作是相同的，全程250公里，所以A车的总数为：
；B的总数为：
。
旅客最大数：
。
问题二：
提高电动汽车的利用率
在问题一中，我们假设每个服务中心是相互独立的，显然，这样当前一中心发出的车经过后一服务中心时，由于时间间隔问题，导致载客量小，即车辆的利用率低，为此，我们再次控制时间间隔，通过动态规划的思想建立效用函数由此达到最大利用率。
首先，我们已知乘客每天的行车时间为3到6小时，而汽车运行时间为9小时，则乘客在服务中心的滞留时间为3到6小时（不包括晚间住宿），故平均每站的滞留时间为0.5到1.5小时，计为
。
我们可以知道，在上一发车时刻已确定从而
已确定的情况下，则下一辆车到达的效用函数
和遗留人数
都是发车时间
的函数
根据问题一中数据计算得
由题意在保证能达到每天旅客数为最大容量的85%这个前提下，则要求每个服务中心的人数要保持在850至1000之间。
以此作为约束条件，效用函数作为目标函数，列出动态规划方程
目标函数
约束条件
问题三：
服务中心位置及规模的设计方案
因为考虑到每天进入公园游览的旅客数成poisson分布，所以为了减少服务中心资源不必要的浪费，根据旅客数的分布情况来设置服务中心的位置和规模。
在人流密集区域扩大规模，增加容纳量，反之，在人流较为稀疏的区域适当缩小规模及容纳量。
误差分析与稳定求解
（1）误差分析
如果不考虑模型自身的缺陷，模型的误差主要来自两个方面：
一是我们假设各个服务中心是相同的，忽略了各个服务中心之间的相互影响；二是没有考虑乘客上下车时间、上下车的人数差。
由于服务中心的人数远大于一辆车的载客量，每个服务中心之间的相互影响不大，故所得结果与真实值之间的误差保持在一个较小的范围之内。
（2）稳定性分析
我们将时间以分钟为单位进行划分，相对于全天9小时的运行时间来说较为精确和细致，车与车之间的时间间隔稳定相等，所以具有较好的稳定性；反之，若有车发生突发事件，因为此事为小概率事件且后面的车辆可通过调节来维持系统的稳定性。
综上，模型是稳定的。
参考文献
[1]甘应爱、田丰、李维铮、李梅生、陈秉正、胡运权、顾基发、郭耀煌、钱頌迪，运筹学（第三版），北京：
清华大学出版社，2005
[2]张无非、张弛、严奇琦，对于公交汽车调度问题的求解，上海交通大学，2002
附录
Matlab运行程序
Z=20*X+8*Y
t=[102030405060];
x=[4322151198];
y=[49251713109];
m=[215461611100];
Z=20*X+8*Y
Z=20.*X+8.*Y
Z=20*x+8*y
p=polyfit（x,z,m）
p=polyfit（x,Z,m）
plot3=（t,z,m）
plot3=（t,Z,m）
plot3（t,Z,m）
plot（t,x,y,m）
plot3（t,Z,m）
plot（t,m）
plot（1/t,m）
q=1/t
q=1\t
plot3（q,Z,m）
meshgrid（t,Z,m）
clc
surf（t,Z,m）
%--2012/6/919:
09--%
plot3（t,Z,m）
clc
t=[102030405060]
x=[4322151198]
y=[49251713109]
n=[215461611100]
plot3（t,x,y,n）
z=x+y
plot3（t,z,n）
surf（t,z,m）
surf（t,z,n）
surf（t,x,y,n）
z=[4322151198
49251713109]
surf（t,z,n）
%--2012/6/1011:
16--%
t=[102030405060]
y=[49251713109]
x=[4322151198]
t=[10152030405060]
x=[432922151198]
y=[4933251713109]
n=[21593461611100]
plot（t,n）
m=1.\t
m=1/t
m=1./t
plot（m,n）
holdon
z=x+y
holdon
plot（m,z）
z=20x+8y
z=20*x+8*y
z=x+y
plot（m,n）
holdon
plot（m,z）
X1=22
15
11
8
6
5
4
X1=[2215118654]
Y1=[2416128654]
N1=[6631118650]
plot（m,N1）
z1=X1+Y1
holdon
plot（m,z1）
x2=[2114117544]
y2=[2517139755]
z2=x2+y2
n2=[6328116430]
plot（m,n2）
holdon
plot（m,z2）
t3=[101112131415]
x3=[433936333129]
y3=[494541373533]
n3=[21519514413212493]
z3=x3+y3
m2=1/t3
m2=1\t3
plot（m2,z3）
holdon
plot（m2,n3）
表二
AB同时发车
A车
乘客人数
B车
乘客人数
间隔十分钟
8:
00
10:
00
20
8:
00
9:
00
8
8:
10
10:
10
40
8:
10
9:
10
16
8:
20
10:
20
60
8:
20
9:
20
24
8:
30
10:
30
80
8:
30
9:
30
32
8:
40
10:
40
100
8:
40
9:
40
40
8:
50
10:
50
120
8:
50
9:
50
48
9:
00
11:
00
140
9:
00
10:
00
56
9:
10
11:
10
160
9:
10
10:
10
64
9:
20
11:
20
180
9:
20
10:
20
72
9:
30
11:
30
200
9:
30
10:
30
80
9:
40
11:
40
220
9:
40
10:
40
88
9:
50
11:
50
240
9:
50
10:
50
96
10:
00
12:
00
260
10:
00
11:
00
104
10:
10
12:
10
280
10:
10
11:
10
112
10:
20
12:
20
300
10:
20
11:
20
120
10:
30
12:
30
320
10:
30
11:
30
128
10:
40
12:
40
340
10:
40
11:
40
136
10:
50
12:
50
360
10:
50
11:
50
144
11:
00
13:
00
380
11:
00
12:
00
152
11:
10
13:
10
400
11:
10
12:
10
160
11:
20
13:
20
420
11:
20
12:
20
168
11:
30
13:
30
440
11:
30
12:
30
176
11:
40
13:
40
460
11:
40
12:
40
184
11:
50
13:
50
480
11:
50
12:
50
192
12:
00
14:
00
500
12:
00
13:
00
200
12:
10
14:
10
520
12:
10
13:
10
208
12:
20
14:
20
540
12:
20
13:
20
216
12:
30
14:
30
560
12:
30
13:
30
224
12:
40
14:
40
580
12:
40
13:
40
232
12:
50
14:
50
600
12:
50
13:
50
240
13:
00
15:
00
620
13:
00
14:
00
248
13:
10
15:
10
640
13:
10
14:
10
256
13:
20
15:
20
660
13:
20
14:
20
264
13:
30
15:
30
680
13:
30
14:
30
272
13:
40
15:
40
700
13:
40
14:
40
280
13:
50
15:
50
720
13:
50
14:
50
288
14:
00
16:
00
740
14:
00
15:
00
296
14:
10
16:
10
760
14:
10
15:
10
304
14:
20
16:
20
780
14:
20
15:
20
312
14:
30
16:
30
800
14:
30
15:
30
320
14:
40
16:
40
820
14:
40
15:
40
328
14:
50
16:
50
840
14:
50
15:
50
336
15:
00
17:
00
860
15:
00
16:
00
344
15:
10
16:
10
352
15:
20
16:
20
360
15:
30
16:
30
368
15:
40
16:
40
376
15:
50
16:
50
384
16:
00
17:
00
392
（1）
AB同时发车
A车
B车
间隔11分钟
8:
00
10:
00
20
8:
00
9:
00
8
8:
11
10:
11
40
8:
11
9:
11
16
8:
22
10:
22
60
8:
22
9:
22
24
8:
33
10:
33
80
8:
33
9:
33
32
8:
44
10:
44
100
8:
44
9:
44
40
8:
55
10:
55
120
8:
55
9:
55
48
9:
06
11:
06
140
9:
06
10:
06
56
9:
17
11:
17
160
9:
17
10:
17
64
9:
28
11:
28
180
9:
28
10:
28
72
9:
39
11:
39
200
9:
39
10:
39
80
9:
50
11:
50
220
9:
50
10:
50
88
10:
01
12:
01
240
10:
01
11:
01
96
10:
12
12:
12
260
10:
12
11:
12
104
10:
23
12:
23
280
10:
23
11:
23
112
10:
34
12:
34
300
10:
34
11:
34
120
10:
45
12:
45
320
10:
45
11:
45
128
10:
56
12:
56
340
10:
56
11:
56
136
11:
07
13:
07
360
11:
07
12:
07
144
11:
18
13:
18
380
11:
18
12:
18
152
11:
29
13:
29
400
11:
29
12:
29
160
11:
40
13:
40
420
11:
40
12:
40
168
11:
51
13:
51
440
11:
51
12:
51
176
12:
02
14:
02
460
12:
02
13:
02
184
12:
13
14:
13
480
12:
13
13:
13
192
12:
24
14:
24
500
12:
24
13:
24
200
12:
35
14:
35
520
12:
35
13:
35
208
12:
46
14:
46
540
12:
46
13:
46
216
12:
57
14:
57
560
12:
57
13:
57
224
13:
08
15:
08
580
13:
08
14:
08
232
13:
19
15:
19
600
13:
19
14:
19
240
13:
30
15:
30
620
13:
30
14:
30
248
13:
41
15:
41
640
13:
41
14:
41
256
13:
52
15:
52
660
13:
52
14:
52
264
14:
03
16:
03
680
14:
03
15:
03
272
14:
14
16:
14
700
14:
14
15:
14
280
14:
25
16:
25
720
14:
25
15:
25
288
14:
36
16:
36
740
14:
36
15:
36
296
14:
47
16:
47
760
14:
47
15:
47
304
14:
58
16:
58
780
14:
58
15:
58
312
15:
09
17:
09
800
15:
09
16:
09
320
15:
20
16:
20
328
15:
31
16:
31
336
15:
42
16:
42
344
15:
53
16:
53
352
16:
04
17:
04
360
（2）
间隔12分钟
8:
00
10:
00
20
8:
00
9:
00
8
8:
12
10:
12
40
8:
12
9:
12
16
8:
24
10:
24
60
8:
24
9:
24
24
8:
36
10:
36
80
8:
36
9:
36
32
8:
48
10:
48
100
8:
48
9:
48
40
9:
00
11:
00
120
9:
00
10:
00
48
9:
12
11:
12
140
9:
12
10:
12
56
9:
24
11:
24
160
9:
24
10:
24
64
9:
36
11:
36
180
9:
36
10:
36
72
9:
48
11:
48
200
9:
48
10:
48
80
10:
00
12:
00
220
10:
00
11:
00
88
10:
12
12:
12
240
10:
12
11:
12
96
10:
24
12:
24
260
10:
24
11:
24
104
10:
36
12:
36
280
10:
36
11:
36
112
10:
48
12:
48
300
10:
48
11:
48
120
11:
00
13:
00
320
11:
00
12:
00
128
11:
12
13:
12
340
11:
12
12:
12
136
11:
24
13:
24
360
11:
24
12:
24
144
11:
36
13:
36
380
11:
36
12:
36
152
11:
48
13:
48
400
11:
48
12:
48
160
12:
00
14:
00
420
12:
00
13:
00
168
12:
12
14:
12
440
12:
12
13:
12
176
12:
24
14:
24
460
12:
24
13:
24
184
12:
36
14:
36
480
12:
36
13:
36
192
12:
48
14:
48
500
12:
48
13:
48
200
13:
00
15:
00
520
13:
00
14:
00
升级会员 