2021年2月
第1期第46卷
昆明理工大学学报(自然科学版)
Joernai of Kunming University of Science and Techuoloyy(Natural Sciences)
Feb.2021
No.1Vol.46
doi：12.16112/1-chdi-53-1223/1-2221.21.20
关于远大前程的论文吴小龙4胡松6，成卫6
(1-昆明阡陌交通工程咨询有限公司，云南昆明654028；9-昆明理工大学交通工程学院，云南昆明654540)
摘要：为了使交叉口的通行效益最大化，构建以交叉口平均延误、平均停车次数最小以及最大通行能力最大化为优化目标的多目标函数模型，并提出一种基于自适应权重和levy飞行的改进鲸鱼优化算法(ALWOA)对信号控制交叉口进行配时方案优化.结合实际算例，将求解后的结果与现有方案、使用传统的Webster算法求解的方案、使用GA和标准WOA得到的方案进行对比.结果表明，使用ALWOA 得到的配时方案明显优于其他方法，从而证明了改进后的鲸鱼优化算法以及函数模型在交叉口信号配时方案优化上的有效性. 关键词：交叉口信号控制；鲸鱼优化算法；信号配时优化;levy飞行
中图分类号:U491.2文献标志码：A文章编号：O26-855X(222021-204-28
Multi-Objective Signai Timing Optimization Based on ImpreveC
Whale Optimization AlgoritUm
WU Xiaolong1,HU Song2,CHENG Wei2
(1.Kunming Qiasmo Tr—fic Engiueering Coxsul/ng Co.Ltd.,Kunming655428,China；
2.Facalty of Transporiatiox Engiueering,Kunminy University of Science and Techuolopy,Kunmine650
540,China) Abstrect:In order te yet tUo maximum benefiis of tUo comprebensive traffic intersection,wo constructeb a mulii 一okjective function mokei te achieve tUo minimum averayo delay,tUo smallest averayo nnmbyr of parking and tUo maximum capacity-Wo alse proposeb an Adaptive Wight and Lxy Flight Improveb Whaio Optimization Ai-OoritUm(ALWOA)to optimiao tUo timiny schemo io tUo sipnal control intersection-Combineb witU tUo actual x-amplos,tUo resulis and tUo schem were cempareb witU tUo schemo ixu UX Uom tUo conventional Webstor Alyo-rithm,GA and standard WOA-Tho ixu/s show that tUo ALWOA timing schemo is supxO/te otUor metUoks, which has proveb tUo ePectiveness of ALWOA and tUo function mokei on tUo intersection sipnal timing optimization-
Key wcras:zntxsxtion sipnal control-Whalo optimization alyoritUm;timing of intersections,leve Uy
o引言
我国经济的高速发展大大刺激了国民的物质需求，私家车出行因为其灵活性和方便性也因此被越来越多的人选择.由于城市化的推进和私家车保有量的快速提高，城市交通问题已成为影响城市经济发展的一个重要因素•我国的早期交通规划化与和现有的通行能力需求并不成正比，城市路网负荷超标严重，导致大量城市出现了拥堵的情况•为了在现有的交通条件下尽可能提升交通效益、减少城市的拥堵状况，对
收稿日期：2019-12-30
基金项目：国家自然科学基金项目(6136400)
作者简介：吴小龙(082-)，男，硕士，高级工程师.主要研究方向：交通系统工程.E-mail：9362648l@qq
通信作者：成卫(1277-)-男，博士，教授，博士生导师.主要研究方向：智能交通管理系统关键技术与理论方法.
E-mail：C69948233@qq
第1期吴小龙，胡松，成卫，：基于改进鲸鱼优化算法的多目标信号配时优化06
交叉口信号控制进行优化是最有效的途径之一.交叉口信号配时优化要考虑多方面的综合因素来进行交叉口信号控制方案设计以达到交叉口的交通效益最优的目的.在评价交叉口通行效益时通常用停车延误、停车率、车辆排放、通行能力等指标来进行.国内外学者在交叉口配时优化上做出了很多研究:Webster配时法是以交叉口车辆延误的估计为基础，通过对周期长度的优化计算，确定相应的一系列配时参数的经典配时算法［1］；国内外还开发了许多有关信号控制的软件系统如Syychrv等⑵；文献［3］针对交通流的时变特性，建立了以延误、排队长度和停车次数最小为目标的交叉口多目标信
号配时优化模型，并结合实例证明了模型的有效性；文献［4］使用解析法来构建函数模型，并引入了基于排队论的理论模型，结合人工操作经验来完成交叉口的信号控制，但是该方式在情况复杂和波动性较大的城市交通中适用性不高.近年来，各种启发式算法的兴起使得交叉口配时的方法有了新的发展方向，已经有大量将启发式算法应用于交叉口配时优化的研究：文献［5］构建了车辆和行人总延误最小的信号配时模型，使用遗传算法寻优获得交通实体产生的平均交通延误达到最小的方案；文献［6］改善了过饱和条件下的信号配时优化模型，以延误最小和通行量最大为目标，使用遗传算法（GA）获得更优的配时方案；文献［7］基于TSTM结合GA算法提出的GATSTM系统能够通过校准系统参数来处理和管理交通网络状况的动态变化；文献［8］通过使用改进的微粒算法对构建的函数模型求解，获得比传统Webster算法更优的配时方案；文献［9］将模拟退火算法引入自适应控制交叉口中，证明了该方法可以明显提升交叉口的通行能力-
鲸鱼优化算法（W0A）作为一种新兴的元启发式搜索算法，在2217年由澳大利亚的mirjaPli等［12］根据观察海洋中座头鲸独特的捕食方式提出，算法通过模拟鲸鱼一系列的捕食行为实现目标函数的求解寻优.W0A算法因其操作简单、参数少而且性能好的特点被大量研究人员所关注并作出了许多相关研究.当前对WOA的研究主要集中在应用和改进两方面.文献：D-17］分别通过引入自适应权重和柯西变异以及随机自适应权重和模拟退火的策略对W0A进行改进，并使用测试函数表明了改进后算法的优越性；文献［13］对WOA进行改进并将其应用到充电站选址的项目上，证明了该算法在工程应用上的有
效性.文献［14］使用WOA进行水资源的配置优化，验证WOA在该领域的适用性.文献［0］则通过对WOA进行电网无功优化调度，实例验证了WOA在解决该问题上的鲁棒性和有效性-
文章提及的遗传算法、模拟退火算法和WOA等在进行求解时都存在着传统启发式算法通有的易陷入局部最优解、收敛精度低等缺陷•本文首先通过使用自适应权重和levy飞行策略对WOA进行改进，提升了WOA的全局寻优能力、局部寻优能力以及收敛精度，并把改进后的算法应用到交叉口信号配时优化中，然后通过构造多目标寻优函数模型并使用ALWOA对其进行寻优求解，获得综合通行效益最高的配时方案-
1多目标信号配时优化的数学模型描述
进行交叉口配时优化时，优化效果的评价指标主要包括延误时间、停车次数、道路通行能力、饱和度、油耗、尾气排放等•本文选取交叉口的车辆平均延误、平均停车次数以及最大通行能力这三个参数作为优化目标，利用加权的方法将这三个优化目标联合起来构造目标函数并进行寻优，以获取交叉口的最大交通效益-
12交叉口平均延误时间
车辆延误由均匀延误和随机延误组成，本文运用韦伯斯特（Webster）提出的延误时间的计算方法来计算车辆的平均延误时间.由Webster延误计算公式可得相位i的车辆平均延误仏如下:
必=c（］-扎）0殆
2（1-打）*2伽（1-切）(0
式中：第一部分表示均匀延误，第二部分为随机延误;9表示周期9是绿信比，表示第i相位有效绿灯时间与信号周期的比值;9i为第i相位第J进口道的流量比；切为第i相位第丿进口道的饱和度为第i相位第j进口道的实际到达的当前交通量jcu/h.由式（1）可知交叉口的所有车辆的平均延误表示为：
136昆明理工大学学报（自然科学版）第46卷
d=丫也/丫％⑵
1.2交叉口平均停车次数
进入交叉口的车辆在信号控制的情况下会产生停车的总次数如式（3）所示:
排放因子h=X仏％/》弘3）
i i
式中：，表示第i相位的车辆平均停车次数.由式（3）可得一个信号周期内的交叉口的车辆平均停车次
数表示为：
h=X也/X弘（4）
i i
1.3交叉口通行能力计算
计算交叉口通行能力的方法：首先将交叉口各进口道划分为若干车道组，然后计算各车道组的通行能力，再将各相位的通行能力加起来，最后得到该交叉口一个周期的通行能力，表达式如下：
q=X a=X s入⑸
i i
式中:Q，表示第i相位的通行能力，血为车道组i的饱和流率.
13目标函数
根据实际到达交通量，通过加权的方式将信号周期内交叉口车辆的平均停车次数、平均延误以及交叉口通行能力联合转化为交叉口信号控制优化目标函数，以控制周期内的有效绿灯时间作为自变量.考虑
到路口交通量的变换，对模型在设置权重时要根据实际交通流率进行分配，权重设置为«3=2|3-y|;2=3.5|3-y=0.5 y（y为交叉口各相位关键流率比之和），目标函数如下：
r d h Q o
（6）约束条件包括以下几点:
X（（■“+4）邑*⑸
5m女=
gm m W&W&m*（8）式中：g a为相位i有效绿灯时间，厶为相位i损失时间,g lmin为相位i的最小绿灯时间，gm#为相位i的最大绿灯时间,5m为最小信号周期长度,5#为最大信号周期长度•由目标函数构成可知，要使目标函数取得最小值，就要求在约束条件下尽可能减小车辆的平均延误和停车次数而通行能力则尽量增大.
2鲸鱼优化算法
鲸鱼优化算法是一种新兴的智能优化算法，鲸鱼优化算法具有简单、调整参数少、性能高效的特点.算法的主要思路是通过模拟海洋中的座头鲸的觅食行为得到的，把食物位置作为寻优目标通过包围捕食、气泡网攻击及随机游走等方式获得最优解•
2.3包围捕食策略
鲸鱼个体能识别猎物的位置区域.由于位置的优化设计在搜索空间不是预先确定的,WOA假定当前的最佳解决方案是目标位置或接近目标最优个体位置.在定义了最佳搜索位置之后，个体开始按照一定的策略朝着当前最优位置进行游动.此行为表示如下：
D=\C•云⑸）-壬⑸）（9）
X⑸+1）=十⑸）-A•D（10）式中/表示个体位置和最优个体间的距离丿表示当前迭代次数，丈表示当前最优个体位置向量,X表示其余个体位置向量,A和c表示系数向量,A、的计算方法如下：
rcc第1期吴小龙，胡松，成卫，：基于改进鲸鱼优化算法的多目标信号配时优化07
A(11)
C=2-r(12)式中&随着迭代次数增加从2到2线性递减的向量，表示一个2到1之间的随机向量.
2.2气泡袭击阶段
描述气泡网攻击行为的数学模型用如下两种方式设计实现：
1)收缩包围机制：这种行为是通过降低式(11)中r的值来实现的，这时候A的值也会随之变化•换句话说,A是区间［-a,a］内的一个随机值，在这个区间内，呈现从2到2的线性递减迭代过程•如果A的随机值在［-1,0之间时，个体会朝着最优位置更新自身位置，更新方式通过式(0)实现.
2)螺旋更新位置：此方法首先计算鲸鱼之间的距离，个体从当前位置朝着最优个体位置进行螺旋式移动，用一个螺旋方程来模拟鲸鱼的螺旋形运动如下：
壬(+0=可•e"•cos(2n)+X*()(13)式中：o7=IX*()-x()1是鲸鱼个体与最优个体之间的距离向量,e是对数螺线形状常数S是［-1,0中的随机数•鲸鱼个体在猎物周围游动收缩包围和沿螺旋形的路径游走是同时进行的•这时为了模拟这两种同时进行的机制，使用2.0的阀值来决定是选择缩小包围机制还是螺旋模型更新个体的位置，数学模型如下:
X(t)-A•D:P<2.5
X(+0=｛一“一(14)
［D'•e•cvs(2tt/)+X*()j MO.5
心宽一寸病退一丈
2.3搜寻猎物阶段
当心生活中的核辐射
该阶段同样基于向量A的变化方法来寻猎物，个体会根据彼此的位置进行随机搜索.因此，我们使用A在大于1或小于-1的范围的随机值来偏离目标位置进行随机搜索，从而扩大搜索范围，这时就是迫使WOA执行全局搜索，从而表现出一定的全局寻优能力•数学模型如下：
D=iC・X n()-X()I(15)
X(+0=X Ud()-A•D(16)式中:X Vd是种中随机一个个体的位置.
3改进的鲸鱼优化算法
为了避免传统的WOA在求解后期容易陷入局部最优导致的算法早熟从而收敛使进度不高的问题，本文改进WOA的思路从以下两个方法入手：一是使用自适应权重方法，使得WOA的局部寻优能力得到提升；另一方法是通过引入levy飞行策略对鲸鱼位置进行更新，以提升WOA的全局寻优能力-
32自适应权重方法
由于WOA的局部搜索实现方式是以公式(17)和公式(13)进行的，当个体以公式(14)的更新方式进行局部搜索时，这种方式只能在局部最优解附近徘徊，而不能实现更好的局部寻优•惯性权重是可以用来平衡算法局部搜索能力和全局搜索能力的重要参数•本文使用了一种呈指数改变的自适应权重方法，算法前期使用较大的权重实现较强的全局搜索性能：保证搜索范围：随着迭代次数的增长，接近最优解时：
权重值呈现指数减小：使得算法的局部寻优能力大大提升.自适应权值公式如(17)所示，改进后的位置更新公式如(18所示:
5J=e-(12!/T)6(0)
138昆明理工大学学报(自然科学版)第46卷
r®•0()-A-D,p<0.5
(18)
•X*()+D'•e b l•cos(2n/),pM0.5
式中：表示当前迭代次数,T表示最大迭代次数.
汤姆叔叔视频18以上观看
3.2levy飞行策略
levy飞行这个概念的正式定义是“步长具有重尾概率分布的随机行走”.我们可以说这是一个随机游动，它的特殊性在于它表现出较大的跳跃，因为这个过程的步长来自于一个具有无限方差的分布•与任何随机过程一样,levy飞行起源于扩散过程.正因为如此，它们在随机测量和随机或伪随机自然现象的模拟中
很有用，特别是表现出一种反常的扩散，系统中存在一种“微观结构”，与混沌理论有关.levy飞行的过程是实体在进行运动的过程中进行大量的小步长移动，同时还有少量跨越式大步长移动的过程•在Matlab
图1Levy飞行二维平面示意图
Fig.2Two-dimensional diagram of Lvvy flight
研究发现许多生物的觅食行为符合levy飞行模式，目前也有大量的研究将其应用到一些仿生算法中，并取得了不错的效果•受文献［16-87］启发，本文将levy飞行应用于鲸鱼的位置更新中，在算法进行更新后再进行一次levy飞行更新个体位置，可以实现跳出局部最优解，扩大搜索能力的效果.位置更新的方式为：
X(+8)=X()+a©levy(A)(19)式中:a是步长缩放因子,levy(X)就是随机步长，㊉就是'•*'运算.2009年，Yang X Y［18_18］把levy分布函数经过简化和傅立叶变换后得到其幕次形式的概率密度函数，使用Mantegna方法生成levy分布随机步长.levy飞行概率密度函数及生成随机步长的公式如下：
levy:：二)f,lw入W3(22)
S(28)即进行更新时levy。)使用S表示；进行计算时，参数0取值为1.5，取值为1,“~N(0,/),i~N(0, 1)，取值为:
厂(1+0)•sin(n0/2)]1
pr^)•o(-1)/2}(22
3.3改进WOA的算法流程图
改进WOA通过使用levy飞行策略来跳出局部最优解，避免算法早熟.通过加上自适应权重的方式使
