560 
深圳大学学报理工版
t;  = C~ L (ωl,k~k+lll,k~k+l Vl ,k~k+l ) 
第 29 卷
(3) 
传算法优化发车间隔和松弛时间.杨兆升 [5 J 提出了
以社会效益最大化为目标的发车频率优化模型.宋
端等 [6] 运用机会约束规划方法,研究随机需求条件
下公交运营设计的优化问题.于滨等 [7] 针对整个公
交系统,以运营总成本最低为目标建立了一个双层
的公交线路发车频率优化模型.
通常一个城市公交系统中公交车辆的规模是确
定的.因此,在设计发车间隔时,若忽略公交车辆
规模的限制,则会影响研究结果在实际中应用.本
研究在车队总规模固定的前提下,以乘客和运营者
总费用最小为目标,建立了一个模型以优化公交线
路的发车间隔.当城市规模较大、公交线路较多
时,该模型较复杂,使用解析算法很难在可接受的
时间内求得结果.遗传算法 (genetic algorithm , 
GA) 是近年发展起来的一种进化算法,巳成功用
于很多领域 [8-9] 本研究采用遗传算法对该模型进
行求解,并通过双种群提高算法的求解性能 [ω
1  发车问隔优化模型
运营者和乘客是公交系统中两个相互权衡的因
素.在乘客数量不变时,运营者希望发车间隔尽可
能大,以减少其可变成本,而发车间隔变大,会增
加乘客的候车时间,无疑增加了乘客的时间费用.
只有同时考虑运营者和乘客双方利益,才能使城市
公交系统获得最大的社会效益.本研究分别建立乘
客总费用函数和运营者总费用函数,并在此基础上
构建公交发车间隔优化模型.
1.1  乘客时间费用函数
其中, tpllssen 为乘客时间费用;可和 c; 分别表示乘客
等车时间费用系数和乘车时间费用系数 Ul. k 为公
交线路 J 第 k 站等待的乘客数;可为乘客等待公交线
路 l 的期望时间,本研究设乘客到达服从均匀分布,
则可 = h/2; ψ 为惩罚系数,与乘客错过的车辆数
目成正相关,反映了由于车辆容量不足时无法上车
的乘客的额外费用,可以用大于 1 的整数或用发车
频率乘以等待的车辆数表示;岛 k 为公交线路 l 第 k
站上因无法乘车而滞留的乘客数,即 δl ,k→k+l
屿 , k
- bi.k;  hl 为公交线路 l 的发车间隔,该变量是模型
的决策变量 ;ωi.k→k+l 为公交线路 l 第 k 站到第 k + 1 
站区间上公交车内乘客的舒适系数,本研究采用拥
挤度来表示,即 ωl. k~k+l 二 Vi.k~k+l h /  (HV);  H 为优
化周期 V 为标准车辆的额定容量 Vm因为车辆的最
大载客数量 ll. k~k+l 为公交线路 1 第 k 站到第 k + 1 站
的平均运行时间 Vl.k~k+l 为公交线路 J 第 k 站到第 k
V l. k k+l  + bi.k  -
+  1 站间车内的乘客数,即 Vl.k~k+l
αi.k; αi.k 为公交线路 l 第 k 站下车的乘客数, αl.k
V l.k -l-.k  X  qu;  qi.k 为公交线路 1 第 k 站乘客下车比率,
即该站点优化周期内的下车人数占从该站点到终点
所有站点下车人数的比例 bu 为公交线路 l 第 k 站
上车的乘客数,即
b"  = [Ul ,k , 屿,k 运 IH/hl 丁 x (凡JaX
l , k 
- 11 H/hl 丁 x (Vmax  - V l. k - l
- V l .k- 1
k  + αl ,k) 
k  + αi.k) ,  otherwise 
(4) 
1. 2  运营者费用函数
运营者费用主要包括燃油消耗、车辆折旧、维
与运营调度相关的乘客时间费用( tpassen )主要
修、司机及售票员的工资等,其表达式为
由乘客等车和乘车两部分费用构成.其中,等车时
间费用 (tn 通常分为两种情况:①乘客到站后,
等待第 1 辆公交车的时间费用,通常用乘客数与平
均期望等车时间的乘积表示;②乘客到站后,未
能乘坐第 1 辆到站的公交车,必须等待下一辆车的
等待时间费用.乘车时间费用 (t; ) 主要指乘客在车
内所耗费的时间费用,包括乘客客观耗费时间和车
内拥挤程度所造成的主观时间费用.乘客时间费用
函数的具体形式为
t阿en = L (t;'  + t;) 
t;' 
c; 二 (μl,kt~ + ψδl,khl) 
(1) 
(2) 
(5) 
1阳 = L 盯
盯= C~I H/hll +  C: 三 IH/hl 丁 x tl ,k
(6) 
其中 Toper表示运营者费用;打表示公交线路 l 的运
营费用 ;c; 和 C; 分别表示公交车辆的费用系数和公
交公司的运营费用系数.
1. 3  公交发车间隔优化模型
k+l 
要使整个城市公交系统的总费用最低,必须找
到乘客总时间费用和运营者费用之间的均衡点,本
研究采用线性加权法将两种费用转化成单目标优化
问题,得到公交发车间隔优化模型为
W passen t passen  + W oper Toper 
min Ttotal 
(7) 
第 6 期
姚锦宝,等:基于双种群遗传算法的公交线路发车间隔优化
561 
s.  t.  hl  ~ 60 
L, gl 运 G
(8) 
(9) 
通过一系列 (0 , 3600J 的随机数组成.由于模型存
在车队规模的约束,因此在产生初始染色体时,如
W  passen'  W op巳r 二", 0
W passell  +  W op{~r  = 1 
(1 0) 
( 11 ) 
其中, T,o,a' 为系统总成本 Wpassen 和 Woper 分别为乘客
时间费用和运营着费用的权重值 ; gl 为线路 l 公交车
辆的规模,即 gl =  12 L, t l ,k--'k+l/hl l ,  Il 表示向下
取整,如 16. 5l = 6;  G 表示公交网络的车辆总数.
1. 4  权重确定
测试前,需要确定优化模型目标函数中乘客时
间费用和运营者费用的权重值 Wpassen 和 Woper " 本研究
以大连市公交系统所有公交线路的乘客时间费用和
运营者费用的比例作为两个权重的取值,
lz() 
果不满足车队规模约束,即 L, gl >  G , 则重新生成
染色体,直至产生一个种群.
2.2  适应值函数
标准遗传算法适宜求解目标最大的优化模型,
本研究通过引人一个常数 Q , 将目标函数 (7) 转换
为最大化函数.此外,遗传算法在进化过程中可能
会产生一些违背模型约束的染色体.通常有 2 种方
法解决这个问题,一种是通过给该染色体增加一个
惩罚因子,来降低其在后续进化中存活的几率,但
是并不强制扔掉该染色体;另一种则是修改不满足
约束的染色体的基因片段,以使其满足约束.本研
究采用第 1 种方法对目标函数进行处理,得到遗传
算法适应值函数为
三( t p…(i)  + 凡川 i))
(1 2) 
F  = Q/[T+ φ(7) ( L, g/  -的]
w叩
1
- Wpassen 
其中, tpassen (i) 和 Toper ( i) 分别表示第 i 条线路乘客
时间费用和运营者费用.假设有 2 条公交线路,它
们的乘客和运营费用分别为 (0.9 , 0.76) 和 (0.8 ,
0.78) ,则有 ω阳ssen = (0.9  + O.  8)/(0. 9  + 0.76  + 
O.  8  + O.  78)  =  O.  52 , w叩盯 =  1  - w阳osen
O.  48. 
2  双种群遗传算法
遗传算法 (GA) 是模仿自然选择、物种进化
和群体遗传学建立的一种随机搜索技术,已被成功
用于工业、经济管理、交通运输及工业设计等领
域 [11] 标准遗传算法在进化过程中,由于生成的
染色体可能重复或相似,在种群规模较小时容易发
生早熟[山3] 与传统遗传算法不同的是,双种群遗
传算法在优化前,同时生成 2 个种群, 2 个种群的
编码和规模完全相同. 2 个种群独立进行进化操作
(交叉和变异) .在进行选择操作时,双种群遗传算
法采用"混合选择"方式,在进化过程中可以保持
种群相对丰富多样,减小早熟导致收敛的可能,提
高求解质量.
2. 1  编码策略
公交发车间隔优化模型网络中每条线路的发车
间隔(秒)为决策变量,本研究采用整数编码的方
案,编码为 1 e l ,  e 2 , …,  e l , … ,  e N  1. 初始种群是
(1 3 ) 
(1 4 ) 
φ(7)  = f ‘
rß, 7β) ,g-,  >  G 
l~. 
l  0 . 
otherwise 
其中 , F 表示适应值函数 ;Q 、 βl 和 β2 是常量, βl 和
β2 用来控制对违反约束染色体的惩罚强度.
2.3  选择操作
双种群遗传算法与标准遗传算法最大区别在于
染色体的选择.双种遗传算法中染色体选择步骤
为:首先,将 2 个种群的所有染色体合并,形成一
个 2 倍于原种群规模的染色体集合;然后,分别从
这个染色体集合中选择各自的新种群.在选择过程
中,为避免新种群中包含重复的染色体,禁止同一
条染色体出现在同一种群中.本研究选用了精英策
略,将适应值最高的染色体复制到下一代,确保最
好的基因在下一代进化中不会丢失,避免退化.
2.4  交叉操作
交叉操作是通过交换两个父染色体的基因来产
生后代.具体交叉操作为
(;「 δyarenl ,
j  +  (1δ) eyurent , L l 
efIld , 2= 「 δ , leyarenl , 2  +  (1  _δ 'l)erare时 , ll , 
币", 
深圳大学学报理工版
第 29 卷
辆的额定容量、最大载客量、乘客的等车、乘车时
间费用系数、公交车辆的费用系数和公交公司的运
营费用系数等数据.根据这些实际数据,对 Wpa
和切0叩阳P严阳巳盯r 2 个目标的权重进行标定.最终模型和算法
中相关参数的设定如表 1 和表 2.
表 1 优化模型的参数取值
Table 1 
Parameters in  headway  optimization  model 
H 
3600 s 
4 
V 
80 
4 
Vmux 
120 
C~ 
2.0元/小时 8.75方车
3 元Y分
'p 
2 
Wpassen 
0.53 
CW 
P 
2. 7 元/小时
W oper 
0.47 
表 2 双种群遗传算法中的参数取值
Table 2  Parameters in DGA 
Q 
1 000 
λ 
3 
Psiz‘祖
240 
R 
p ,-
0.6 
Pm 
o.  1 
number 
β1 
β2 
Tmax 
2000 
3.1  发车间隔优化结果
为检验本研究提出的发车间隔优化模型的效
果,将优化结果与现状进行比较.优化前后整个网
络的总费用分别是 57.0 万元和 52. 2 万元,优化方
案可节省 8.4% ,优化结果分别见图 1 和图 2. 可
以看出,优化前后各条公交线路的发车间隔发生了
很大波动,同时各线路的舒适性(满载程度)也发
生了较大变化.优化前有些线路过于拥挤,有些线
路却乘坐率甚低,说明优化前的公交车辆的资源未
得到充分利用,优化后公交系统内各线路的满载率
则能够获得有效改善(现状所有线路满载率的均方
562 
因 ;e;lnhl ,
l 和 e~俨hil队l
6矶F:t 和 8矶"飞F:l 表示[阳0 , 1 ]间的随机数 ; PC 是交叉率.
2.5  变异操作
变异操作是通过改变一条父染色体中的一段基
因来产生新的染色体.它可以引人新的基因信息到
种群中.与交叉操作相似,变异操作也是通过一个
变异率 Pm 来控制变异的数量.如果父染色体的第 J
段基因被选择进行变异,则变异的结果为
'hEn 
tujq 
「
l
== 
-A1i (( ×× nn ee rr aa 
。
0
m
'
'
b
、
l
/
+
一
(( 
xx aa nn TT 
JJJ/ 
A
λ
 
「l
l
l
,, 
δ[  <  PC 
《
,
4
r
t
E
F
、
J
E
l
l
-
‘
ρ
l
v
ρ
t
v
 
已
,
ι
「
l
l
p
a
ρ
,
u
ι
P
ρ
A
U
F
ι
 
 
。
。
川m
z
e
b
、
}
/
T 
「l
i
t
otherwise 
(1 6 ) 
其中, efaM 表示父染色体的第 l 段基因 ;efzkl 表示子
染色体的第 l 段基因 ;δl 和 δl 表示 [0 , 1] 间的随机
数;凡是变异率 'Tmax 和 T 分别表示最大的进化代
数和当前的进化代数.此外,在算法最初阶段,染
色体变异的程度越大,意味着算法可以在更大范围
内搜索,能够丰富种群的多样性;而算法进化到后
期,染色体过大的变异程度则会影响收敛速度,本
研究设置了参数 λ(2 '"二 λ 运 5) 来控制进化代数对
变异程度的影响.
3 
实例研究
为验证上述模型和算法,本研究采用大连市主
城区的公交数据进行了样本研究.大连市市区面积
为 180 km2 ,目前共有的条公交线路, 3  004 个公
交站点.对 89 条线路进行随车调查,获得了早高
峰期间的客流 OD 数据、运行时间、高峰小时、车
+现状线路发车间隔
·优化线路发车间隔
500 
400 
300 
200 
Z
E
E
-恃
抖
@+ 
·+ 
@+ +· 
@+ 
@
亨
.+ • 协
+· 
l 
nu 
nu 
i
陆
• • +  ... 
+ 
.. 
+ 
Q  +. + 
+ 
+ 
...  . .. 
.. 
+. 
+ - : -
+ 
-+ 
..- +、
+ 
.. 
+ 
Q 
..  .. 
+ 
@ 
• 
• • +· 
+ 
• • +@ 
• @+ 
@ 
•  ..+  .. 
+  + 
+ 
++ 
+  .. 
+‘' 
+ 
.. .. 
+ 
@+ 
+ 
@+ 
+ 
' 
• 
.. • 
·+ 
·+ 
.. 
+@ +· 
@ +· 
+@ +· 
+@ 
+@ 
+ 
+· 
• ·+ 
-7 
·+ 
·+ 
4· 
+ • 
89 
85 
81 
+-
+曹+ + 
+ 
ø T 
~  't 
Q 
Q 
+ 
Q 
++ 
0.... 
"'Q 
Q 
-:-‘' 
8' 
-:-. 
~+ 
-'1 
.no'!' 
Q
-:- +'a 
.". 
Q 
+ 
。
5 
9 
13 
17 
21 
25 
29 
33 
37 
41  45  49 
公交线路
53  57 
61 
65 
69 
73 
77 
图 1 优化前后备线路发车间隔比较
Fig. 1  Comparison of headway  of each  route 
第 6 期
姚锦宝,等:基于双种群遗传算法的公交线路发车间隔优化
563 
2.0 
+现状线路舒适水平 ·优化线路舒适水平
1.6~  + 
+ 
+ 
+ 
\ l 2 ·+ . . . . …+ 
叶
怜明越E
+  9 .@  
0080+  ~  Q  0 
0
+  Q 
.  G+ 
4b 
0  + 
_+ 
o
品!
+ 
+ 
4· 
• 
+ 
+ 
• • +@ +· 
O 8+ .+++  + 
+ 
+ 
+ 
+ 
• 
+ 
+ 
+· +· 
• 
• • ·+ 
@ 
• 
• 
·+ ·+ 
@+ 
·+ 
• 
• 
·+ 
• 
·+ 
+· 
·+ • + 
• • 
+ 
@+ @+ 
@+ 
+ 
@+ 
@ 
+ 
+@ 
+· • 
azT • • , 
·+ 
• 
+ 
+++ 
• • • 
+· ·+ 
+ 
+ 
·+ • • +· 
+ + 
00 
+ 
• 
• 
@ 
+ 
O.4 r 
。
+ 
+ 
+  + 
+  + 
00+  
+ 
5 
9 
13 
17 
21 
25  29  33  37 
41 
45  49  53  57 
61 
65 
69  73  77 
81 
85  89 
公交线路
圄 2 优化前后各线路舒适度比较
Fig. 2  Comparison of comfortable level  of each  route 
误差为 0.35 ,优化后为 0.23) .因此如果对目前大
连市现有的公交车辆资源整合,并进行统一调度,
将会节省整个公交系统的总费用.
3.2  双种群遗传算法性能测试
为验证双种群遗传算法( dual population genetic 
algorithm ,  DGA) 的求解效率,在相同参数下,本
研究进行了 10 次试算,结果如图 3. 可以看出,前
600 代内适应度迅速增加,但从 800 代开始,算法
适应度的变化程度逐渐减小,到 1 200 代左右算法
开始趋于收敛 .10 次测算结果相差不大,最差的结
果与最好结果相差大约 2% ,说明 DGA 算法的收敛
性能良好.与此同时,为检验 DGA 与 GA 算法的性
能差异,本研究将 DGA 与 GA 在相同的情况下连续
运行 10 次,计算结果如图 4. 可以看出,双种群遗
传算法在求解公交线路发车间隔优化问题上优于标
准遗传算法,这进一步说明,双种群机制是一种保
持种群多样性的有效手段.
18 
自 15
明
12 
9 
200 
400 
800 
600 
进化代数
1000  1200  1400 
圄 3 双种群遗传算法的测试结果
Fig. 3  The result of each  calculation 
..  DGA 
'"  GA 
19.2~ 
A 
喝
' 
• 
A 
日画
~ 
也
.. 
a 
@ 
A 
• 
句画
• 
A 
生吕
因F司
18.8 
18.4 
18.0 
.. 
• 
A 
.. 
2 
3 
4 
6 
5 
测试次数
7 
8 
9 
10 
图 4 双种群与单种群遗传算法性能比较
Fig. 4  Comparison of the  perfonnances 
of single  and dual  GA 
为分析罚函数中参数 βl 、 β2 对 DGA 算法性能
的影响,本研究设计了 5 种参数组合来测试算法的
收敛速度和求解质量,具体结果如图 5 所示.组合
(β2 , β2  = 2) 的收敛速度最快,但其解的质量
不高.组合 (β1 =  0.5 , β2  =  0.5) 搜索解的质量很
高,但是其收敛的速度却非常惶组合 (β1β2
= 1) 则具有很好的搜索性能和求解质量.原因在
·适应度 ........收敛速度(进化代)
2200 
1900 
.:: 
19.5 
19.2 
刨 18.9
国
甘冒 18.6 
18.3 
T~白ld
V4P 3 
11 
CQ.. 
参数组合
圄 5 罚函数参数比较
Fig. 5  Comparison of parameters in penalty function 
564 
深圳大学学报理工版
第 29 卷
于,随着 2 个参数的增大,罚函数对违背约束的染
色体进行了更加严厉的惩罚,使其很快从种群中消
失.但是,由于搜索方案很容易产生一些不满足约
束的解方案,这些方案虽然违背约束但是其很可能
会存在部分有益的基因段,如果过快淘汰这些基因
可能会影响求解的质量.因此,对于进化过程中,
容易产生不可行方案的 DGA 应该考虑采用较大的
控制参数(乱, β2) ;反之,则应选择较小的参数.
结语
发车间隔优化对公交运营至关重要,本研究从
乘客和公交运营者双方利益出发,考虑城市公交车
辆资源的约束,开发了一个同时考虑乘客时间费用
和运营者费用最低的发车间隔优化模型.考虑到该
模型属于公交调度问题,设计了双种群遗传算法对
其进行了求解.以大连市主城区公交系统的数据对
该模型和算法进行了检验.结果显示,如果对大连
市公交车辆资源进行整合,统一调度会改善整个系
统的服务水平,同时也可以降低系统的总成本.此
外,通过与标准的遗传算法进行比较表明,双种群
遗传算法的求解效果优于标准遗传算法.
基金项目:国家自然科学基金资助项目 (51208079 , 51108053)
作者简介:姚锦宝 (1972-) ,男(汉族) .安徽省安庆市人,北京
交通大学讲师、博士. E-mail:  jbyao@bjtu.edu.cn 
引
文:姚锦宝,姚宝珍,尹智宏,等.基于双种群遗传算法的
公交线路发车间隔优化[1].深圳大学学报理工版,
2012 ,  29(6):  559δ64. 
参考文献/References:
[  1  J Scheele  S.  A  supply  model  for  public  transit  services 
[ 1].  Transportation  Research:  B ,  1980 ,  14:  133 -146. 
[  2  J Sun  C  J ,  Zhou  W ,  Wang  Y  Q.  Scheduling  combination 
and  headway  optimization  of  bus  rapid  transit  [J J.  Jour
nal  of Transportation  Systems  Engineering  and  Information 
Technology ,  2008 ,  8 (5) :  61 -67. 
[  3  ]  Tom  V  M,  Mohan  S.  Transit  route  network  design  using 
frequency  coded  genetic  algorithm  [J J.  J  Transp  Een
asce ,  2003 ,  129(2):  186-195. 
[  4  ]  Park S J.  Bus  Network  Scheduling with  Genetic  Algorithms 
and  Simulation  [D J.  Maryland  (USA):  University  of 
Marγland ,  2005. 
[  5  ]  YANG  Zhao-sheng.  Theory  and  Method  of  Urban  Intelli
gent  Public  Transportation  System  [M J.  Beiji吨 Railway
Press ,  2002:  56-76.  (  in  Chinese) 
杨兆升城市智能公共交通系统理论与方法[ MJ 
北京:中国铁道出版社, 2002:  56-76. 
[  6  J SONG  Rui ,  HE  Shi-wei ,  YANG  Hai ,  et  al.  Optimal  tran
sit  operation  plan  model  and  algorithm  with  stochastic  OD 
demands  [J J.  China  Civil  Engineering Joumal ,  2006 ,  39 
(4):  110-115.  (in Chinese) 
宋瑞,何世伟,杨海,等.基于随机需求的公交运
营设计优化模型及算法[J J. 土木工程学报, 2006 , 
39( 4):  110-115. 
[  7  ]  YU  Bin ,  Y ANG  Zhong-zhen ,  CHENG  Chun-tian ,  et  al. 
Bi-level  programming  model  for  optimizing  bus  frequencies 
and  its  algorithm  [J J.  Joumal  of Jilin  University 
Engi
neering  and  Technology  Edition ,  2006 ,  36 (5) :  664-668. 
(in Chinese) 
于 滨,杨忠振,程春田,等.公交线路发车频率优化
的双层规划模型及其解法[1] .吉林大学学报工学
版, 2006 ,  36  (5):  664-668 
[  8  J CAI  Liang-wei ,  HU  Shi-xi.  A  genetic  algorithm  based  on 
similarity  and  its  application  on  JSP  [J J. 
Joumal  of 
Shenzhen  University  Science  and  Engineering ,  2006 ,  23 
(2):  107-11 1.  (in Chinese) 
蔡良伟,胡世曦-基于相似性遗传算法及其在 JSP 中
的应用[ J J .深圳大学学报理工版, 2006 ,  23  (2): 
107 -11 1. 
[  9  J WU  Xu-yi ,  WU  Xiao-yu.  Improved  genetic  algorithm  for 
job-shop  scheduling  [J J.  Joumal  of  Shenzhen  University 
Science  and  Engineering ,  2006 ,  23  ( 3 ):  272-277.  (in 
Chinese) 
吴序一,伍晓宇.非量产模式下车间调度的改进遗传
算法[J] .深圳大学学报理工版, 2006 ,  23  ( 3 ) : 
272-277. 
[ 10 J Zhao  G  H ,  Shen  W  H ,  Wu  M L.  Ultra-wideband  antenna 
design  using  FDTD  and  double  population  genetic  algo
rithm  [J J.  Microwave  and  Optical  Technology  Letters , 
2009 ,  51 (2) :  361-364. 
[ 11  J Holland  J  H.  Adaptation  in  Natural  and  Artificial  Systems 
[ M J.  Michigan  (USA):  University  of  Michigan  Press , 
1975.234-245. 
[12J  LI  Bi ,  LIN  Tu-she吨. Modified  genetic  algorithm  based  on 
competitive  coevolution  [J J.  Joumal  of Shenzhen  Univer
sity  Science  and  Engineering ,  2009 ,  26 ( 1 ):  24-29.  (in 
Chinese) 
李 碧,林土胜.基于竞争协同进化的改进遗传算法
[J].深圳大学学报理工版, 2009 ,  26  (1):  24-29. 
[ 13 J CAI  Liang-wei.  Real-coded  adaptive  genetic  annealing  al
gorithm  based  on  distance  measurement  [J J.  Joumal  of 
Shenzhen  University  Science  and  Engineering ,  2004 ,  21 
(4):  291-294.  (in Chinese) 
蔡良伟.基于距离视~度的实数编码自适应遗传退火算
法 [JJ. 深圳大学学报理工版, 2004 , 21  (4):  291-
294. 
[中文责编:坪梓;英文责编:之幸]