1、1关于公交车调度问题摘要随着国民生活水平的提高,公共交通问题也日益重要起来,而公交车调度是制约公共交通的重要因素。根据题中所给的数据,建立数学模型对公交车调度问题进行分析。对于问题一:首先,根据城市中某条公交线路各个时段的客流信息,得出了公交车公司的最大客容量,发车车次,发车时间间隔。运用MATLAB 编程,计算出各个时段的最大客容量,在满足公交满载率的情况下得出日最少发车车次为460次,其中上行线230车次,下行线230车次,用LINGO计算出发车时间间隔,并给出公交车发车时刻调整表。基于公交车从起始站运行到终点站的用时为44分钟,且时间间隔应为整分间隔,可算出早高峰所需最少车辆为58辆。其
2、次,一个合理的公交车调度方案应该考虑公交公司的最大利益和乘客的满意度两个方面。故建立了满意度分析模型,在此模型中,运用了层次分析法。对满意度进行了分析计算。结合整数规划模型中的结果可求得满意的分析模型中公交公司与乘客双方之间满意度,并且使二者和达到最大,同时双方满意度之差最小,得到上下行的最优满意度(0.8688,0.8688) 。最后,综合了公交车公司的最大客容量、发车车次、公交公司满意度等方面因素,且以公交公司所发的车次最小为目标,乘客的等待时间和公交载客率为约束条件提出了整数规划模型。此模型是把公交车调度问题抽象成数学模型来表达,从考虑发车车次最小出发,满足各项约束条件,寻求最优解。运用
3、LINGO编程,可计算出公交公司日发车车次最小值为461次。因此该解法是在满足乘客的情况下求的最优解。乘客的等待时间的满意度为100%,但是从舒适度考虑,上行和下行分别有11和9人不满意。这个结果为满意度模型和整数规划模型的中间情况,故此模型的建立是合理的。关键词:整数规划 满意度 MATLAB LINGO2一 问题的重述公共交通是城市交通的重要组成部分,作好公交车的调度对于完善城市交通环境、改进市民出行状况、提高公交公司的经济和社会效益,都具有重要意义。下面考虑一条公交线路上公交车的调度问题,其数据来自我国一座特大城市某条公交路线情况,一个工作日两个方向上下车的乘客数量统计表如表 1、表 2
4、 所示。已知调度要求如下:该条公交线路上行方向共 14 站,下行方向共 13 站,公交公司配给该线路同一型号的大客车,每辆标准载客 100 人,据统计客车在该线路上运行的平均速度为 20 公里/小时。运营调度要求,乘客候车时间一般不要超过 10 分钟,早高峰时一般不要超过 5分钟,车辆满载率不应超过 120%,一般也不要低于 50%。需要解决的问题:(1)为该线路设计一个便于操作全天(工作日)的公交车调度方案,包括两个起点的发车时刻表;一共需要多少辆车;这个方案以怎样的程度照度到了乘客和公交公司双方的利益:等等。(2)如何将这个调度问题抽象成一个明确、完整的数学模型,指出求解模型的方法。二 问
5、题的分析本问题需要从顾客与公交公司两方面利益下手,从而制定使双方都满意的公交调度方案。从顾客角度考虑,顾客做希望的是公交等待时间较少、公交上超载程度越小越好;从公交公司的角度考虑,公交公司希望在能满足交通需求的条件下,公交公司所安排的公交数量越少越好,从而满足公交公司利润最大。仅考虑提高公交公司的利益,只要提高公交车的满载率即可,运用数据分析法很容易得到其分配方案;仅考虑方便顾客出行,增加公交车数量即可,运用统计法,我们可以很容易的得到其调度方案。在对这两方面进行分析以后,我们考虑公交公司和顾客的满意度,在保障运行通畅的3情况下,我们选用最合理的调度方案。显然这两方案是对立,于是本题将分成两部
6、分进行分析:(1)公司经济利益用公司满意度表示,记为: ;mg(2)顾客等待时间及乘车舒适度用顾客满意度表示,记为: 。c三 符号说明符号说明如表 1:表 1 符号说明表i上行运动(此时 )14,2k2下行运动(此时 )318,j公交车从 5:00 到 23:00 运行的各个时间段ijl上行或下行第 时段最大载客量jijz上行或下行第 时段平均载客量ijc上行或下行第 时段的整车次j日所需总发车次D高峰期所需要的最大车辆数ijs上行或下行第 时段平均发车时差jijF上行或下行第 时段发车时差为小数时,向下取整数ijsC上行或下行第 时段发车时差为小数时,向上取整数jimc上行或下行乘客的日平均
7、满意度ij上行或下行第 时段乘客满意度jijt上行或下行第 时段乘客等车时间tmc乘客对等车时间的满意度w乘客对乘车舒适度的满意度4img上行或下行公交公司日平均满意度ij上行或下行第 时段公交公司的满意度j四 模型的假设(1) 交通情况、路面状况良好,不出现意外交通事故,公交车之间无超车现象;(2) 公交车车速以理想车速运行即:20 公里/小时;(3) 发车时间间隔取整数分钟数,公交车之间发车时间间隔不超过 20 分钟;(4) 乘客按顺序依次上车,不允许插队。五 模型的建立与求解5.1 模型一5.1.1 计算最大客容量(1)本文已经把数据分成上行方向和下行方向 18 个时段进行了处理,考虑到
8、每个时段乘客量不同,建立以下模型:(51)jmibalmkijkijkijkijij 18,2,18,2,)(x),1(1 将式子(1)用 MATLAB 编程得(程序如附录一)上下行各个时间段的最大客容量如表 2表 2 上行、下行每个时段的最大客容量表上行 701 2943 5018 2705 1528 1193 1355 1200 1040下行 27 1039 2752 3223 1822 1093 986 830 891上行 881 871 2133 2772 897 464 410 275 19下行 1017 1302 2196 3612 2417 1091 781 714 3375为了
9、更直观的观察最大客容量的数据以及早高峰的时间段,由此绘制各个时间段最大客容量的直方图(程序如附录二)如图 1: 图 1 上下行各时间段内最大客容量直方图由图 1 可看出:(1)在上行路线,上午的乘客要比下午的多,而下行路线,上午与下午的乘客差不多;(2)结合上下行的图可得出,早高峰在 7:008:00 这段时间,且客容量为 5000多;5.1.2 计算各个时间段最少发车次数由于公交车标准载客为 100 人,车辆满载率在 50%120%之间,这里求的是最小发车次数,所以取车辆满载率为 120%,即 人,由模型:120ijz(其中 是正整数) (52)Z120l,lllcijij ijijij12
10、0, Z(53)18ijicC6可以计算出各时间段的发车次数 ,还要注意一点,公交车之间发车时间间隔不ijc超过 20 分钟一趟,于是得到发车的车次,如表 3(程序如附录三):表 3 上、下行各时段的发车次数表时间段 上行路 下行路 时间段 上行路 下行路5:006:00 6 3 14:0015:00 8 96:007:00 25 9 15:0016:00 8 117:008:00 42 23 16:0017:00 18 198:009:00 23 27 17:0018:00 24 319:0010:00 13 16 18:0019:00 8 2110:0011:00 10 10 19:002
11、0:00 4 1011:0012:00 12 9 20:0021:00 4 712:0013:00 10 7 21:0022:00 3 713:0014:00 9 8 22:0023:00 3 3结合上下行发车车次,以及式子(3)得到总的最少发车次数为: 。46023182jiicC5.1.3 安排发车时间间隔取每个时段 60 除以车次数,得到各个时段的平均发车时间间隔:(54)csijij60以公式(4)求得上下行各个时段的平均时间间隔依次如表 4 所示:表 4 各时刻平均发车时间间隔表上行 10 2.4 1.4 2.6 4.6 6 5 6 6.7下行 20 6.7 2.6 2.2 3.8
12、6 6.7 8.6 7.5上行 7.5 7.5 3.3 2.5 7.5 15 15 20 20下行 6.7 5.5 3.2 1.9 2.9 6 8.6 8.6 20观察表中数据发现有小数存在,与模型假设不相符,应该调整为整分间隔。当7取整数时,安排的车次为 ;而当 不为整数时,不妨设以 , 为时间间隔ijsijcijs ijsFijC的发车次数分别为 ,建立模型:ijnm,(55)18,2;,160 jicsCsFijijij ij由公式(5)得到 之后,为了使得安排在同时段线路的车辆不宜过多,我们ijn,对调整的整分发车间隔对应发车量的先后顺序作调整,将相邻时间段内发车间隔相等的班次尽量安排
13、在一起,对数据进行处理。可以得到上下行各时段发车时间间隔调整表(如附录表 5)5.1.4 日需车辆数由汽车平均速度 20 千米/小时和 A0A13 的距离 14.61 公里、A13 A0 的距离14.58 公里,可求得车辆从起点站运行到终点站平均用时为 44 分钟;又由假设可知车辆到达终点后立即掉头返回。由于早高峰乘客数最多,故此时车辆实际占用数应是当日的上限,若公交公司日派车最少时能达到这个用车上限,则能满足日需车辆数。考虑到最少车辆数应满足上下行的公交车发车要求,上行方向比下行方向车辆要多发车,我们根据各时段的发车车次 cij,公交车单程运行时间 44 分钟,以及高峰期的时间段 7:008
14、:00,可以得出高峰期的占用车数为:=44+14=582314sD得到高峰期实际占用车数为 58 辆,其中 A13 站需要 44 辆,A0 站需要 14 辆,也就说明公交公司日需车辆数最少为 58 辆。5.2 模型二5.2.1 满意度的分析在问题分析中,已提到满意度的问题,下面本论文将对 , 进行分析。mgc(1)公交公司满意度对于公交公司来说,车辆的利用率越高,期望利用率越好,公交公司发车车次就少,对公交公司利益就越大。在乘客源一定的情况下,影响 的主要因素是车上的ij8乘客数即载客量 ,一般情况下 。各个时间段的公司满意度表示为ijz1205ijz从而建立模型如下:,ijmg(56)120
15、ijijzmg(57))1,2.8,;j(icg18jiijjii (2)乘客满意度对于乘客来说,影响 的主要因素是乘客的等车时间 与车上的平均载客量 ijmc ijt ijz。设 , 分别是各时段乘客因 与 的影响而产生的满意度,则 即可以表itjmciwj ijtijz ijmc示为:(58)Acciwjitjij,其中,A 是关于因素 与 的权重集。ijtijz对于乘客, 的影响不是相等的,上下车的乘客都在动态地变化,但ijiwjitjmcc对,对于车辆而言,车辆的满载率达 120%时,最大超载的 20%由于缺少座位,而注重舒适度的影响,无暇过分顾及等待时间的影响;100%的乘客因为有座
16、,而无需过分考虑舒适,更多的是考虑等车时间的影响。又设 A= ,其中, 分别是因素 的重要程度,用层次分析中的成对比witawita, ijzt,较法,可知:5201wita同时,A 应满足归一性和非负性,即0,1witwiti a可以解得 (59)6,5witia9因此 (510)iwjitjwitijitjij mcamcc615),( 我们把 满意度函数看着是常见的降半梯形分布ijitj,(511)10510ttcitj(512)22wmcjw由每时段的乘客满意度 ,每时段的乘客最大客容量 ,一天最大客容量人数为ij ijl,可以算出乘客平均日满意度为各时段的满意度的加权平均值:18ji
17、l(513))18,2;,1(18 jilmccjiijjii(3)合理调度情况分析对于公交公司,当满载 120 人时公交公司最满意,人数越少,满意度越来越低。对于乘客,可知当等车时间不超过 5 分钟,车辆满载率不超过 100%时,乘客满意度为1,随着等待时间增加和车载率的上升,乘客满意度会逐渐下降。取当公交车平均载客人数分别为 120 人,100 人,50 人时作分析。5.3 模型三在问题分析过程中,题中给了如下约束条件:(1) 乘客的等待时间一般不超过 10 分钟;(2) 早高峰时间段内乘客的等待时间不超过 5 分钟;(3) 各个时间段内的最大乘车率不超过 120%;(4) 各个时间段内的
18、最小乘车率不低于 50%。可以以公交公司所发的车次为最小,列出下面的目标函数:(514)Zczjjj117)(min10s.t. (515)%cll)j,(cjjjjj 12010541062511这个模型是整数规划模型,在满足各种约束条件的情形下,寻求全天发车车次的最小值,将式子(13) 、 (14)用 LINGO 编程求解(程序如附录五) ,算法流程图如下:开始输入数据找出可能答案目标函数是否为最小值输出数据结束是否满足各约束条件六、模型检验模型是把公交车调度问题抽象成数学模型来表达,从考虑发车车次最小出发,满足各项约束条件,寻求最优解,于是可以利用这个模型来分析此问题,对条件分析可11知
19、,约束条件满足两方面,一方面要满足乘客的等车时间早高峰不超过 5 分钟,其余时段不超过 10 分钟。对于公交公司方面,也要满足客车的载客率在 50%120%之。对于题中的客流量,我们筛选出不合要求的时段,如:第 18 时段、下行第 1 时段。于是我们利用 LINGO 编程(程序如附录三) 。得到的发车车次情况如表 6:表 6 上下行各时间段发车车次情况时间段 上行路 下行路 时间段 上行路 下行路5:006:00 6 3 14:0015:00 8 96:007:00 25 9 15:0016:00 8 117:008:00 42 23 16:0017:00 18 198:009:00 23 2
20、7 17:0018:00 24 319:0010:00 13 16 18:0019:00 8 2110:0011:00 10 10 19:0020:00 4 1011:0012:00 12 9 20:0021:00 4 712:0013:00 10 8 21:0022:00 3 713:0014:00 9 8 22:0023:00 3 3一天总发车车次为 461 次,因此该解法是在满足乘客的情况下求的最优解。乘客的等待时间的满意度为 100%,但是从舒适度考虑,上行和下行分别有 11 和 9 人不满意。此模型的结果为模型和的中间情况,故此模型的建立是合理的。七 模型的评价与推广7.1 评价(1
21、)该模型有效地求解出最优解;(2)计算早高峰的最少占用车数,很有代表性;(3)本模型运用了加权平均,对满意度进行分析,得出结果比较合理;(4)本模型中没讨论交通堵塞,红路灯,公交故障等现实状况的发生,有点脱离实际情况。7.2 推广根据前面的模型所建立的运输系统,可以很好地解决公交线上公交车的调度问题,12然而,在建模过程中,简化了许多因素,因而与实际问题有偏差。因此,要想建立更好的调度方案,可以对一条实际运营的公共汽车的运行过程进行计算机模拟,将调查得到的实际数据输入计算机程序,便可以得出更优的调度方案。八 参考文献吴建国等,公交车调度方案的优化模型,建模案例精编,中国水利水电出版社,2005 年第一版。