公交车调度问题数学建模论文

2011年数学建模论文——对公交车调度问题的研究摘要：本文根据所给的客流量及运营情况排出公交车调度时刻表，以及反映客运公司和乘客的利益有多个指标，建立了乘客的利益及公司利益两个目标函数的多目标规划数学模型。基于多目标规划分析法，进行数值计算，从而得到原问题的一个明确、完整的数学模型，并在模型扩展中运用已建的计算机模拟系统对所得的结果和我们对于调度方案的想法进行分析和评价。首先通过数据的分析，并考虑到方案的可操作性，将一天划为；引入乘客的利益、公司利益作为两个目标函数，建立了两目标优化模型。通过运客能力与运输需求(实际客运量)达到最优匹配、满载率高低体现乘客利益；通过总车辆数较少、发车次数最少表示公司利益建立两个目标函数。应用matlab中的fgoalattain进行多目标规划求出发车数，以及时间步长法估计发车间隔和车辆数。关键字：公交车调度；多目标规划；数据分析；数学模型；时间步长法，matlab一问题的重述：1、路公交线路上下行方向各24站，总共有L辆汽车在运行，开始时段线路两端的停车场中各停放汽车m辆，每两车可乘坐S人。这些汽车将按照发车时刻表及到达次序次发车，循环往返地运行来完成运送乘客的任务。建立数学模型，根据乘客人数大小，配多少辆车、多长时间发一班车使得公交公司的盈利最高，乘客的抱怨程度最小。假设公交车在运行过程中是匀速的速度为v。1路公交车站点客流量见下表从新汽车站出发到市检察院站点名称新汽车站汉庭花园天九湾电信公司天九湾车场西环小区步行街上车人数1131112下车人数10004等待时间3.8521.52站点名称实验小学莆一中后门十字街旧汽车站新街口市农行上车人数513483下车人数100101等待时间3.835.314.13.8站点名称市公交公司中国银行凤凰山八十亩小区石室路口市公交稽证处上车人数313322下车人数233923等待时间1.710.52.52.25.5站点名称北磨交通花园三信家园市政府龙桥市场市检察院上车人数200000下车人数4225710等待时间4从市检察院出发到新汽车站站点名称市检察院龙桥市场市政府三信家园交通花园北磨上车人数1730127下车人数011114等待时间3.51.22.84.82.64站点名称市公交稽证处石室路口八十亩小区凤凰山中国银行市公交公司上车人数231583下车人数101221等待时间3.31.65409站点名称市农行新街口旧汽车站十字街莆一中后门实验小学上车人202410数下车人数272522等待时间3.4651站点名称步行街西环小区天九湾车场天九湾电信公司汉庭花园新汽车站上车人数001000下车人数3124314等待时间11已知数据及问题的提出我们要考虑的是莆田市的一路公交线路上的车辆调度问题。现已知该线路上行的车站总数N1(=24)，下行的车站总数N2(=24)，并且给出每一个站点上下车的人数。公交线路总路程L（=L）；公交行驶的速度V=20km/h；运营调度要求,车辆满载率不应超过r=120%,一般也不要底于r=50%。现要我们根据以上资料和要求,为该线路设计一个公交公司发车时间的调度方案、一共需要多少辆车、公交车道路行驶过程中的速度以及公交车车型的选择的方案。并给出刻划乘客和公交公司双方利益、满意程度的指标,进行评估等。2准备工作我们首先来看一下上、下行线的有关客流量数据。由给出数据对数据进行处理，可简化模型求出每分钟的平均客流量。从新汽车站出发到市检察院站点名称新汽车站汉庭花园天九湾电信公司天九湾车场西环小区步行街上车平均客流量110.80.20.50.71下车平均客流量0.30002站点名称实验小学莆一中后门十字街旧汽车站新街口市农行上车平均客流量1.30.30.642.00.8下车平均客流量0.300100.3站点名称市公交公司中国银行凤凰山八十亩小区石室路口市公交稽证处上车平均客流量1.8161.20.90.4下车平均客流量1.2363.60.90.5站点名称北磨交通花园三信家园市政府龙桥市场市检察院上车平均客流量0.500000下车平均客流量1225710从市检察院出发到新汽车站站点名称市检察院龙桥市场市政府三信家园交通花园北磨上车平均客流量4.91.500.20.81.8下车平均客流量00.80.40.20.41站点名称市公交稽证处石室路口八十亩小区凤凰山中国银行市公交公司上车平均客流量0.61.90.51.380.3下车平均客流量0.300.20.520.1站点名称市农行新街口旧汽车站十字街莆一中后门实验小学上车平0.600.30.810均客流量下车平均客流量0.670.3122站点名称步行街西环小区天九湾车场天九湾电信公司汉庭花园新汽车站上车人数001000下车人数3124314请注意这个表格是对数据的简单处理，剔除、修正一些不合理的数据,并且以人/分钟为单位保留一位小数进行四舍五入。需要注意：下车平均客流量是将在每个站点的下车人数进行时间划分，是在以无论等多长时间乘客都不会离开的前提下假设的；表格中出现的0不是说此站没有人，而是客流量太小可以忽略不计。然后将各站的上车平均客流量和下车平均客流量之和进行比较,大于或等于3.0人次每分钟的定为大站,认为在该站上下车需耗时2分钟,即Δti=2;1.0人次每分钟至3.0人次定为中等站,消耗1分钟；上下车耗时1分钟小于5000人为小站,上下车耗时0.5分钟。3问题的初步分析及基本假设制定公交车调度方案需要考虑的因素非常多,且很多因素都是随机的。为了抓住重点,简化模型建立及求解,必须作一定的简化假设和设定。1)汽车从起点站发车后,都能在额定的时间里到达终点站;2)汽车行驶过程都看做匀速行驶；3)乘客在规定的时间内都可以乘车;4)乘客的满意程度只以他所乘的车的拥挤程度来衡量;5)在车站等待的人绝大多数不会离去。6)公交站点确定，距离的调整不考虑特殊情况7)根据给定数据客流量是一个平均值，则考虑调度问题只需考虑一个时间段即可。8）车辆上行或下行到达终点时，所有的乘客必须下车；9)在同一个时间段内，相邻两辆车发车时间间隔相等；10）对全天而言客车公司基本把所有的顾客运完；4模型的建立4.1符号说明:N某时段发车次数（注：由于数据给定为平均客流量只需考虑在一个完整的周期内的车次，即从始发站到终点站的这段时间）B某时段的平均满载率T=L/v+Σti一辆公交从始发站到终点站的整个时间ai第i站上车平均客流量R=T*ΣaiB=R/(c×N)R为某时段的总上车人数,c=100人/车次α供求匹配比α=(ΣV)/(ΣQ)k控制参数Q某时段运客能力(人×公里)Q=某时段发车次数Ni×每辆车标准载客量c×单程(上行或下行)总运行距离L。其中,上行时,L公里;下行时,L公里V某时段的需要运客量(人×公里)V=Σj(xj-yj)*T*Ljj∈(24,12...,1,0),上行方向;j∈(0,2,3,...24),下行方向。其中,xj为某时段内Aj站的上车人数;yj为某时段内Aj站的下车人数Lj为Aj站距该单程方向上终点站的距离。问题一:发车次数的确定依据前面的分析,兼顾乘客与公交公司双方的利益,分别对单程的上行路线和下行路线建立如下的多目标规划模型:目标函数:Ⅰ供求的最优匹配min(Q×B-V)^2Ⅱ各时段的发车车次均最小min{N}约束条件:①各时段的平均满载率限制0、5≤B≤1、2②供求匹配比限制α≤k4、2目标函数说明:目标函数Ⅰ使某时段的运客能力Q与运输需求(实际客运量)V达到最优匹配,β反映满载率高低的影响。目标函数Ⅱ使所需的最大发车次,在满足约束条件下尽可能少,以使总车辆数较少。4、3约束条件说明:条件①是限制满载率满足运营调度要求,是考虑了程客的利益。条件②是限制供求匹配比α小于常数k。补充约束条件:为使始发站车场的每天起始时刻的车辆数保持不变,需使总发车次数与总收车次数相等,即必须使单程车次总数达到匹配(N1=N2),而N1不能减少(受满载率限制),因此我们在求解下行方向的Ni时增加约束ΣN2i=N1.在增添约束条件ΣN2i=N1之后,用二次规划求得各时段发车次数N1i和N2i。问题二：发车数量及发车间隔的确定(1)发车间隔的确定在这部分,我们采用时间步长法,根据假设一个时段内发车间隔时间ti相等,则ti可由N确定,从而得到发车时刻表。按此发车时刻表模拟实际运行过程,目标是确定满足时刻表的最小车辆数n,统计各项运营指标,搜索最优调度方案解。(2)模拟子程序一:确定最小车辆数目n根据“按流发车”和“先进先出”的原则,对起点站,在发车时刻应至少有一辆车可以发出(处于等待发车状态)。若有多辆车,则先进站者先发车,其余车辆“排队”等候;若无车可发,则出现“间断”。完整的运营过程应保证车辆严格按时刻表发车,不发生间断。设A23站和A0站分别有车场A和B,从车场中不断有车发出,同时接受车进场,则车场中的车的数目是随时间变化的状态量。用Na和Nb来描述车场A和车场B中要满足车流不间断所需的最小数目,分别搜索其在运行过程中的最大值,则所需最小车量数目n=Na+Nb。(3)模拟子程序二:统计各项运营指标确定各项运营指标,采用模拟统计的计算方法,对不同的运营指标进行定量计算,主要功能是通过定量分析运营指标来检验方案的可行性,以确定方案调整。由于车次与发车时刻一一对应,而车辆的队列顺序是不发生改变,因而对所需车辆进行统一编号,则对每一车次,与其对应的车辆编号是确定的,故我们直接对第k次车进行考察。我们统计的指标及其定义如下:平均满载率上行方向B01=(ΣkΣj1B(k,j1)/(N1·J1)下行方向B02=(ΣkΣj2B(k,j2)/(N2·J2)满载率分布可以由B(k,j)确定。平均候车时间上行方向T1=(ΣkΣj1T(k,j1)/(N1·J1)下行方向T2=(ΣkΣj2T(k,j2)/(N2·J2)符号说明:D(k,j)第k次车到第j站时上车与下车的人数之差;(已知)C(k,j)第k次车离开第j站时站台上的滞留人数;C(k,j)=C(k-1,j)+D(k,j)-(120-B(k,j-1)B(k,j)第k次车离开第j站时车上的人数;B(k,j)=B(k,j-1)+D(k,j)+C(k-1,j)-C(k,j)T(k,j)为第k次车离开第j站时站台上滞留者的滞留时间;T(k,j)=C(k,j)·tiβ(k,j)为第k次车离开第j站时的满载率,β(k,j)=B(k,j)/100;N1,N2为一天单程所发的车次总数;J1,J2为单程站台总数;4.4试验过程及数据处理程序一x=[11;0.5;0.2;0.5;0.7;-1;1;0.3;0.6;3;2;0.5;0.6;-2;0;-2.4;0;-0.1;-0.5;-2;-2;-5;-7;10];l=11.500:-0.5:0;T=67;V=T^2*l*x程序二functionoutput=myfun(n)output=zeros(2,1);output(1)=(ni*c*l*bi-vi)^2;output(2)=ni;%其次在命令窗口中输入如下语句;A=[0,1;0,-1];b=[1.2;-0.5];options=optimset('MaxFunEvals',1000000,'MaxIter',1000000,'TolX',0.1);n=fgoalattain(@myfun,[3,1],[10,10],[-10,-10],A,b);4.5模拟结果及统计指标分析设全程的长度为L=12km,车的平均速度为v=20km/h，假设站点间的距离为等距的。根据大小站及距离得出上行方向需要的时间为67分钟，下行方向需要的时间为66.5分钟，上行方向乘客需求量V=189下行方向乘客需求量V=90。根据程序得上行N=10,B=1。下行N=9,B=1。上行平均每5分钟发一辆车，下上行平均每6分钟发一辆车5调度方案:认为在该路线上运行的总车数固定不变,形成序贯流动的车流,依照“按流开车”和“先进先出”的原则,按发车时刻表发车。6模型的进一步讨论1)关于采集运营数据的讨论由于我们假设在一个时段内乘客到站服