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

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

————————————————————————————————作者: ————————————————————————————————日期:

20１１年数学建模论文

——对公交车调度问题的研究

摘要:本文根据所给的客流量及运营情况排出公交车调度时刻表，以及反映客运公司和乘客的利益有多个指标，建立了乘客的利益及公司利益两个目标函数的多目标规划数学模型。基于多目标规划分析法，进行数值计算，从而得到原问题的一个明确、完整的数学模型，并在模型扩展中运用已建的计算机模拟系统对所得的结果和我们对于调度方案的想法进行分析和评价。

首先通过数据的分析,并考虑到方案的可操作性，将一天划为；引入乘客的利益、公司利益作为两个目标函数，建立了两目标优化模型。通过运客能力与运输需求(实际客运量) 达到最优匹配、满载率高低体现乘客利益;通过总车辆数较少、发车次数最少表示公司利益建立两个目标函数。应用matｌaｂ中的ｆgoalａttain进行多目标规划求出发车数,以及时间步长法估计发车间隔和车辆数。

关键字：公交车调度；多目标规划;数据分析；数学模型；时间步长法,matlaｂ

一问题的重述：

１、路公交线路上下行方向各２４站,总共有L 辆汽车在运行，开始时段线路两端的停车场中各停放汽车m辆,每两车可乘坐S人。这些汽车将按照发车时刻表及到达次序次发车，循环往返地运行来完成运送乘客的任务。建立数学模型,根据乘客人数大小,配多少辆车、多长时间发一班车使得公交公司的盈利最高,乘客的抱怨程度最小。假设公交车在运行过程中是匀速的速度为v。

1路公交车站点客流量见下表

从新汽车站出发到市检察院

站点名称新汽车

站

汉庭花

园

天九湾电信

公司

天九湾车

场

西环小

区

步行街

上车人

数

1１３１ 1 1 2 下车人

数

1 ０0 ０ 4 等待时

间

3．8 5 2 1.5 2

站点名称实验小

学

莆一中

后门

十字街旧汽车站新街口市农行

上车人５ 1 ３ 4 8 3

下车人

数

1 0 ０１0 1 等待时

间

3．8 3 ５.３ 1 4.１ 3.８

站点名称市公交

公司

中国银

行

凤凰山

八十亩小

区

石室路

口

市公交稽

证处

上车人

数

3 1 3 ３ 2 2

下车人

数

２ 3 3 9 ２ 3

等待时

间

1.７１０．5 2．５２.2 5．5 站点名

称北磨

交通花

园

三信家园市政府

龙桥市

场

市检察院

上车人

数

２0 0 0 ０0 下车人

数

4 2 ２

5 7 １0 等待时 4

从市检察院出发到新汽车站站点名

称市检察院

龙桥市

场

市政府三信家园

交通花

园

北磨

上车人

数

1７ 3 0 １ 2 7 下车人

数

0 1 1 1 1 4 等待时

间

3．5 1.2 ２.８４.8 ２．6 ４

站点名称市公交稽

证处

石室路

口

八十亩

小区

凤凰山

中国银

行

市公交

公司

上车人

数

2 ３１５8

3 下车人

数

1 0 1

2 ２ 1 等待时

间

3.3 １．6 5 ４０9 站点名

称市农行新街口

旧汽车

站

十字街

莆一中

后门

实验小

学

上车人

数

２0 2 4 1 0 下车人

数

２7 2 5 ２２等待时

间

3.４６ 5 1

站点名

称步行街

西环小

区

天九湾

车场

天九湾电信

公司

汉庭花

园

新汽车

站

上车人

数

0 ０ 1 0 ００下车人

数

3 １ 2 ４３１

4 等待时

间

1

１已知数据及问题的提出

我们要考虑的是莆田市的一路公交线路上的车辆调度问题。现已

知该线路上行的车站总数Ｎ1 ( = ２4 )，下行的车站总数N2 ( ＝

24 ),并且给出每一个站点上下车的人数。公交线路总路程L(＝L)；

公交行驶的速度V＝2０km/ h;运营调度要求，车辆满载率不应超过r

＝１２０% ,一般也不要底于ｒ＝ 50 %。

现要我们根据以上资料和要求，为该线路设计一个公交公司发车时间的调度方案、一共需要多少辆车、公交车道路行驶过程中的速度以及公交车车型的选择的方案。并给出刻划乘客和公交公司双方利益、满意程度的指标,进行评估等。

2准备工作

我们首先来看一下上、下行线的有关客流量数据。由给出数据对数据进行处理,可简化模型求出每分钟的平均客流量。

从新汽车站出发到市检察院

站点名称新汽车

站

汉庭花

园

天九湾电信

公司

天九湾车

场

西环小

区

步行街

上车平

均客流

量

11 0.8 0．2 0.5 ０．７１

下车平

均客流

量

０．３0 ０0 2

站点名称实验小

学

莆一中

后门

十字街旧汽车站新街口市农行

上车平

均客流

1.3 0.3 0.6 4

2.0 0.8

下车平

均客流

量

0.3 0 0 １0 ０．3

站点名称市公交

公司

中国银

行

凤凰山

八十亩小

区

石室路

口

市公交稽

证处

上车平

均客流

量

1.81６１.20．90.4

下车平

均客流

量

１．2 3 6 3.６０.9 ０.5 站点名

称北磨

交通花

园

三信家园市政府

龙桥市

场

市检察院

上车平

均客流

量

0.50 0 ０0 0

下车平

均客流

量

1 2 2 5 7 10

从市检察院出发到新汽车站站点名

称市检察院

龙桥市

场

市政府三信家园

交通花

园

北磨

上车平

均客流

量

4．9１.50 ０．2 0.8 1.8

下车平

均客流

量

０0.8 ０．4 ０.2 ０.4 １

站点名称市公交稽

证处

石室路

口

八十亩

小区

凤凰山

中国银

行

市公交

公司

上车平

均客流

量

0.6 1.90.5１．3８0．3

下车平

均客流

量

0.3 0 ０.2 ０.5 2 ０．1 站点名

称市农行新街口

旧汽车

站

十字街

莆一中

后门

实验小

学

上车平

均客流

量

0.６0 ０.3 ０．8 1 0

下车平

均客流

量

0.６7 0．3 1 2 2 站点名

称步行街

西环小

区

天九湾

车场

天九湾电信

公司

汉庭花

园

新汽车

站

上车人

数

0 0 1 0 0 0 下车人

数

３ 1 2 4 3 14 请注意这个表格是对数据的简单处理，剔除、修正一些不合理的数据,

并且以人／分钟为单位保留一位小数进行四舍五入。需要注意:下车

平均客流量是将在每个站点的下车人数进行时间划分,是在以无论等

多长时间乘客都不会离开的前提下假设的；表格中出现的0不是说此

站没有人，而是客流量太小可以忽略不计。

然后将各站的上车平均客流量和下车平均客流量之和进行比较,

大于或等于3.0人次每分钟的定为大站，认为在该站上下车需耗时2

分钟,即Δt i = 2 ;1.0人次每分钟至3.0人次定为中等站,消耗1分

钟;上下车耗时1 分钟小于50０0 人为小站,上下车耗时0.５分钟。

3 问题的初步分析及基本假设

制定公交车调度方案需要考虑的因素非常多，且很多因素都是随机的。为了抓住重点,简化模型建立及求解，必须作一定的简化假设和设定。

1) 汽车从起点站发车后,都能在额定的时间里到达终点站;

２) 汽车行驶过程都看做匀速行驶；

3）乘客在规定的时间内都可以乘车;

4) 乘客的满意程度只以他所乘的车的拥挤程度来衡量;

5) 在车站等待的人绝大多数不会离去。

6) 公交站点确定，距离的调整不考虑特殊情况

7) 根据给定数据客流量是一个平均值，则考虑调度问题只需考虑一个时间段即可。

8）车辆上行或下行到达终点时，所有的乘客必须下车；

9)在同一个时间段内,相邻两辆车发车时间间隔相等；１0）对全天而言客车公司基本把所有的顾客运完;

４模型的建立

4. 1 符号说明:

N 某时段发车次数(注:由于数据给定为平均客流量只需考虑在一个完整的周期内的车次，即从始发站到终点站的这段时间)

Ｂ某时段的平均满载率

Ｔ=L/v+ Σｔｉ一辆公交从始发站到终点站的整个时间

aｉ第ｉ站上车平均客流量

R =T*Σai

B＝R ／( c ×N) R 为某时段的总上车人数, c ＝100 人/ 车次

α供求匹配比α= ( ΣV) / ( ΣQ)

k 控制参数

Q 某时段运客能力（人×公里）

Ｑ = 某时段发车次数Ni ×每辆车标准载客量c ×单程(上行或下行) 总运行距离L 。其中,上行时， L公里; 下行时, L公里

V 某时段的需要运客量(人×公里）

V = Σj（ｘj-yj）＊Ｔ* Lj j ∈(２4 ,12 ．..,1 ,0) , 上行方向； j ∈ (0 ,2 ,3 , ．.．24），下行方向。

其中, x ｊ为某时段内A ｊ站的上车人数; ｙj为某时段内Ａｊ站的下车人数

Ｌ j 为A j 站距该单程方向上终点站的距离。

问题一:发车次数的确定

依据前面的分析,兼顾乘客与公交公司双方的利益,分别对单程的上行路线和下行路线建立如下的多目标规划模型:

目标函数：Ⅰ供求的最优匹配 min ( Q ×Ｂ- Ｖ )^ 2

Ⅱ各时段的发车车次均最小 mｉn{ N}

约束条件: ①各时段的平均满载率限制０、5 ≤Ｂ≤1、２

②供求匹配比限制α≤k

４、2 目标函数说明:

目标函数Ⅰ使某时段的运客能力Q 与运输需求(实际客运量) Ｖ达到最优匹配,β反映满载率高低的影响。

目标函数Ⅱ使所需的最大发车次,在满足约束条件下尽可能少, 以使总车辆数较少。

4、3 约束条件说明：

条件①是限制满载率满足运营调度要求,是考虑了程客的利益。

条件②是限制供求匹配比α小于常数k。

补充约束条件：为使始发站车场的每天起始时刻的车辆数保持不变,需使总发车次数与总收车次数相等,即必须使单程车次总数达到匹配( Ｎ1 ＝ N2) ，而N１不能减少(受满载率限制) ，因此我们在求解下行方向的Ｎｉ时增加约束ΣN２ｉ = Ｎ1. 在增添约束条件ΣＮ2 i = N１之后,用二次规划求得各时段发车次数N1 i 和N２i 。

问题二:发车数量及发车间隔的确定

(１）发车间隔的确定

在这部分,我们采用时间步长法，根据假设一个时段内发车间隔时间t i 相等,则ｔ i 可由N确定,从而得到发车时刻表。按此发车时刻表模拟实际运行过程, 目标是确定满足时刻表的最小车辆数n ,统计各项运营指标，搜索最优调度方案解。

(2) 模拟子程序一:确定最小车辆数目n

根据“按流发车”和“先进先出”的原则,对起点站，在发车时刻应

至少有一辆车可以发出（处于等待发车状态) 。若有多辆车,则先进站者先发车,其余车辆“排队”等候；若无车可发,则出现“间断”。完整的运营过程应保证车辆严格按时刻表发车,不发生间断。

设Ａ 23站和A 0 站分别有车场A 和B ,从车场中不断有车发出,同时接受车进场,则车场中的车的数目是随时间变化的状态量。用Na 和Nb 来描述车场A 和车场B 中要满足车流不间断所需的最小数目,分别搜索其在运行过程中的最大值,则所需最小车量数目ｎ＝ Nａ+ Nb。

(3)模拟子程序二:统计各项运营指标

确定各项运营指标，采用模拟统计的计算方法, 对不同的运营指标进行定量计算, 主要功能是通过定量分析运营指标来检验方案的可行性，以确定方案调整。

由于车次与发车时刻一一对应,而车辆的队列顺序是不发生改变,因而对所需车辆进行统一编号,则对每一车次，与其对应的车辆编号是确定的,故我们直接对第ｋ次车进行考察。

我们统计的指标及其定义如下：

平均满载率上行方向Ｂ0１＝ ( ΣkΣj1B( ｋ, j１) ／（Ｎ1 ·J1）

下行方向B０２= ( ΣkΣj2B( k , j2) / ( N2 ·J2)

满载率分布可以由B（ k , j）确定。

平均候车时间上行方向T1 = ( ΣｋΣj１T ( ｋ, ｊ1) / （ N １·J1)

下行方向T2 = ( ΣkΣj2T ( k , j２) / ( Ｎ

2 ·Ｊ2)

符号说明:

D ( k ， j) 第k 次车到第j 站时上车与下车的人数之差； (已知) Ｃ(ｋ, j) 第k次车离开第j 站时站台上的滞留人数; C( k ， j) = C( k - 1 , j) + D ( k ， j） -（120 - Ｂ（ k , j - 1)

B ( k , j) 第k 次车离开第j 站时车上的人数; B ( k , ｊ）= B ( k , ｊ - 1) + D ( ｋ， j) + C( k -１, j) - Ｃ( k , j) Ｔ( k ，ｊ）为第k 次车离开第j 站时站台上滞留者的滞留时间; Ｔ ( k , j） = Ｃ( k , j) ·ｔｉ

β( k , j) 为第ｋ次车离开第ｊ站时的满载率，β( ｋ , j) = B ( k , ｊ) / １00 ；

N1 ， N２为一天单程所发的车次总数；J1 , Ｊ２为单程站台

总数;

4．4试验过程及数据处理

程序一

ｘ＝[1１;0．5；０.２;0．5；0.7；-1;１；0.3;0.６;３;2；0.5；０.６;-2;0;-2.4;０;－0.1;－０.5;-2;－2；-５;－７；10]；

l=11．50０:-０.5:０；

T=67;

V=T＾２＊l*x

程序二

fｕｎcｔion ｏutput=mｙfun(n)

output=zｅros(2,1）;

oｕｔｐuｔ(1）=(ni＊c*ｌ*bｉ－vｉ)^２;

outpuｔ(2)=nｉ;

%其次在命令窗口中输入如下语句;

A=[0,１;０,-１]；

b=[1.2;－0.5］;

ｏｐｔions=oｐtimｓet（'MaxＦunＥvals',1０0０000,'ＭaxＩte ｒ＇,10000０0,'ToｌX',０.１);

ｎ=fgoalaｔtａin（@myｆｕn，[3,1]，[1０,１0]，［-10，－10]，Ａ,b) ;

4.5模拟结果及统计指标分析

设全程的长度为L=12kｍ,车的平均速度为v＝２0km/ｈ,假设站点间的距离为等距的。根据大小站及距离得出上行方向需要的时间为67分钟，下行方向需要的时间为6６.5分钟，上行方向乘客需求量V= 189 下行方向乘客需求量V=９０。根据程序得上行N= 1０，B＝１。下行N= ９，B= １。上行平均每5分钟发一辆车,下上行平均每６分钟发一辆车

5调度方案：

认为在该路线上运行的总车数固定不变,形成序贯流动的车流,依照“按流开车”和“先进先出”的原则,按发车时刻表发车。

6 模型的进一步讨论

1) 关于采集运营数据的讨论

由于我们假设在一个时段内乘客到站服从均匀分布, 而实际中乘客到站时间不可能都服从均匀分布。特别是数据给出的是一个公交车在每一站的经验值就会使模型结论误差较大。我们建议以下几种改进采集方式的方法：

(1）由于实际问题中统计的人数不可能服从均匀分布，而且存在一天中不同时间段分高峰期和低峰期。所以数据的统计应按照一天为单位分段统计更加合理。

(2) 采取不等的统计人数的间隔时间在高峰期的情况下,为削弱乘客到站时间不均匀分布带来的影响，可适当减小统计的间隔时间但统计时间加密应有一定限度。对客流量很小的时段，我们可适当增大统计的间隔时间。

(３）增加能反应有关滞留人数的统计数据。

（4）按相等到站人数来区分时间段的统计方法是统计达到一定到站人数时的时间点,其优点是能较为准确地反映客流量的变化情况，有利于按其分布的疏密进行车辆调度，以更好的满足乘客的需要。

(５）考虑到实际问题中,公交线路存在单车场和双车场。单车场问题相当于将重点站看作一个中间站点，即把双车场模型的方法延伸到4８个站点进行单方向运算。方法同双车场的上行运算，所以不再重复。(６)单车场调度方案与双车场调度方案的选用

由结果分析可知单车场调度方案减少了公司的前期投资成本;双车场调度方案的运营成本小,更好的兼顾到乘客与公司双方的利益。我们建议，在有双车场的条件下选取双车场调度方案更好。当需进行路线规划，需要选取单车场或双车场时, 建议根据实际所需成本来选取方

案。

7 模型的评价

本文的优点如下：

1) 模型的主体是采用时间步长法,模拟生成的发车时刻表的实际运行过程,准确性高,容量大，逻辑性严格,计算速度快,具有较强的说服力和适应能力。

2) 定义了能定量衡量我们的调度方案对乘客和公交公司双方利益满足程度的统计指标。

3）在求最少车辆数时，将两个车场看作两个发射源, 通过对两个车场的存车状态的实时模拟,形成不间断的运营过程，从而求得所需车辆数目。

本文的缺点是:

１) 对于运营数据的采集方式,只给出了一些原则和想法,没有经过

仿真验证。

２) 对于乘客到站的分布,直接将人流量假设为均匀分布，没有对其他分布的情况再作讨论。可将分布时段划分为高峰期，平常期,低峰期。

3）对于车辆行驶速度假设为匀速,在实际运营过程存在,站点的停留。未将停留时间考虑进去可能产生误差。