红绿灯亮起时间方案设计

1、论文题目：红绿灯亮起时间方案设计姓名： 专业： 姓名： 专业： 姓名： 专业： 方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。摘 要本文就典型城市道路红绿灯实时配时问题进行研究：首先，以单点红绿灯配时为设计方案，针对尽可能顺畅的目标，我们考虑用延误时间和滞留车辆数两个量作为路口顺畅程度的评价指标，由此引入延误率，滞留率及阻塞度的概念，在尽可能顺畅的要求下，极小化延误率和平均滞留车辆数，由此建立目标函数，设计相应算法。其次，将整个典型环路尽可能顺畅,定义为一个周期内从某路口流向下一路口并无阻碍通过的车辆数尽可能多。采用图解的方法，提出离开率、离开率时间积、滑动窗口的概念，建立出混合线性规划

2、模型，得出在同一时刻各十字路口交通灯的相位分布，最后作出了运作时序图。模型一中，两个目标函数都是比较复杂的非线性函数，约束条件中也因含有max等函数成为非线性，而变量约束为整数，故整个模型一是非线性规划，并且利用计算复杂性的理论，我们可以大致确定其模型属于NP难问题。对于NP难问题，现在还没有多项式算法，不能在较短的时间内得出最优解(通过将优化模型输入LINGO，LINGO不能找出全局最优解)，故我们选择采用启发式的禁忌搜索算法，禁忌搜索能够同时拥有高效性和鲁棒性，它的算法设计忌搜索是对局部领域搜索的扩展，得出了较好的结果。以所有路口的阻塞度总和为目标，得出典型环路上各十字路口的公共周期时长为

3、120s，以及在该周期时长下的红绿灯配时方案，结果见表14。由此，完成了对整个典型环路的红绿灯配时方案的优化设计。该方案给出了共用周期时长下的每个路口红绿配时方案和各个路口交通灯的相位时序图。关键词： 双目标优化规划模型 混合线性规划模型 单点红绿灯配时 滞留率阻塞度 禁忌搜索算法 LINGO软件方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。一问题提出随着城市汽车数量不断增加，每天的上下班高峰期，主要交通路口都会出现堵车现象。 为了增加道路的通过能力，减少堵车现象的发生，红绿灯能够发挥重要的作用。请建立一个数学模型，确定合理的红绿灯亮起时间。二问题分析十字路口是组成城市道路网的基本单元，

4、对其进行研究是开展城市交通系统研究的前提和基础。城市交通控制分为单交叉口控制和多交叉口协调控制(即按控制范围分为点控、线控和面控)，并非后者比前者控制效果好，它们各有自己的适用范围。单交叉口控制不但更适于那些孤立的交叉口，而且还担负着线控、面控控制方案的落实。以下引入一组概念：相位：一个交叉口中可同时获得通行权的一组互不冲突的交通车流。交通流量：单位时间内通过某个位置的车辆数。单位一般为辆/小时。周期时长：各相位车流轮流获得通行权一次所需的时间，单位为秒。相位差：相邻路口同一相位绿灯或红灯起始时间之差。绿信比：一个周期内内一个相位所获得的有效绿灯时间长度与整个周期时长的之间比值。针对典型环路设

5、计出的优化方案能提供如下信息：1) 典型环路上每个路口的红绿灯配时方案和周期时长。2) 典型环路上所有路口在同一时刻的相位分布，并按共用周期时长轮转。对某一个交通繁忙的十字路口，每一条支路按车行方向可分为三种车道：左行，右行和直行。每个车道的交通流量因时间、路口、支路、行车方向的不同而不同。每一个路口的交通灯红绿配时方案及周期时长也可能不尽相同。而且，同一时刻，各个路口的交通灯相位分布也可能不一致。这些都反映了该问题背景的复杂性。结合上述方案的设计内容，我们需要对问题做进一步的简化：第一，只考虑上下班高峰期时的交通问题，这使的交通流量相对稳定，从而进行离线方案设计。第二，结合目前国内道路交通管

6、理状况，只讨论典型四相位三车道的共用周期时长的信号联动管理，并将直行和右行视为一个车道，如图1所示：图1 路口的四个相位由前面相位的定义可知：在任一时刻，某路口的车流通行状态必处于四个相位之一。配时方案模型：模型一：点交叉口交通红绿灯配时方案方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。我们主要考虑设置红灯时长要足够长以保证行人顺利穿过人行横道线的约束，来达到每个路口交通安全的要求。通过每个路口在一个周期内的平均滞留车辆总数最少和每一辆车平均延误时间最少的双目标优化问题，达到尽可能顺畅的要求。求解时，为减少变量个数，将周期时长在一定范围内按一定步长取值，得到相应的一系列配时方案，从中选择

7、最优的一组。在此基础上，以各路口总体阻塞度为目标，确定出共用周期时长和每个路口在该周期下的红绿配时方案。模型二：多交叉口交通信号灯配时方案在整个典型环路尽可能顺畅方面，我们考虑用一个周期内从某路口流向下一路口并无阻碍通过的车辆数尽可能多来量化，做出时间离开率图，并定义离开率时间积和滑动窗口等概念，结合图形，建模求解。关于方案的可行性论证，从系统模拟的角度进行。三变量说明1) ：第个路口，第相位绿灯时长。（单位：秒） 2) ：第个路口，第相位红灯时长。（单位：秒）3) ：第个路口交通灯的周期时长。（单位：秒）4) ：人穿过人行横道线的典型时间15秒。（单位：秒）5) ：第个路口在第相位时第进口道

8、的交通流量，是一个平均值。（单位：辆/）6) ：等待车队中相邻两车的车头间距。（单位：米）7) ：等待车队在绿灯放行时的离开速度。（单位：米/秒）8) ：反应时间：信号灯由红变绿时刻到等待队列中第一辆车开动时刻的时间间隔。该时间与司机的视觉灵敏度，从大脑意识到采取动作的反应速度有关，取典型值2秒。（单位：秒）四基本假设1) 模型中所涉及的原始交通数据相关的数据真实有效，不考虑特殊情况。2) 模型中研究大的车辆为标准车辆，且符合一般规律。3) 不考虑黄灯的时长，因为对司机的驾车行为而言，黄灯和红灯均表示禁止通行。4) 路口处不发生事故，所有司机遵守调度规则.5) 假设直行和右转共用一个相位，即红

9、灯时间禁止车辆右转，这种禁止措施对行人安全是有利的。6) 离线配时方案设计时，某个路口某个方向的车辆交通流量设定为常数，不同路口、不同方向的交通流量方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。不同。五模型建立5.1单点红绿灯配时的方案设计基于上述问题分析，本问题是一个双目标非线性优化问题，要求保证各路口交通安全和尽可能顺畅的前提下，确定出所有路口共用周期时长并设计出每个路口的红绿灯配时方案。我们先分别对各个路口的交通管理进行优化，确定出各个路口各自的周期时长和相应的红绿灯配时方案，然后对所有路口的交通管理进行整体优化，确定出共用周期时长和在该周期下各路口交通灯红绿配时方案。5.1.1不

10、同周期时长下的单点红绿灯配时方案设计下面省去路口下标，对某个十字路口的交通状况进行分析，如图2所示：其中分别表示四个相位,分别表示十字路口的四个进口道，表示某路口处于第相位时第进口道的交通流量（辆/小时）。从优化问题的三要素加以分析：a. 决策变量：四相位车流轮流获得通行权一次所需的时间，即周期时长该路口第个相位的绿灯时长，表示第个相位对应的车流通行时间。第个路口第相位红灯时长，表示第个相位对应的车流禁止时间。b. 约束条件：1) 结合生活中的实际情况，为保证进口道1，3上的行人能安全的穿过人行横道，一、三相位的红灯时长，必须大于行人穿过公路的典型时间，即：方案范文无法思考和涵盖全面，最好仔细

11、浏览后下载使用。2) 在任一时刻，各路口的交通流必处于四个相位之一，故每个相位的开通时间即为该相位的绿灯时长，四个相位的开启时间之和即为周期时长，因而有：3) 结合实际假定每一个相位的绿灯时长至少为秒，此处取典型值10秒，为缩小搜索范围，4) 定义1：为第相位时，第进口道的等待队列长度。其数值上应等于第相位时第进口道的交通流量与第相位的禁止时间（即红灯时长）的乘积。即：定义2：为第相位最大队列容量。表示第相位时路口所能放行的最大等待队列长度。其在数值上应等于等待队列的离开率与第相位的绿灯时间（这里假设四个进口道上等待队列的离开率都相等）。则：其中，：等待车队中相邻两车的车头间距（单位：米）。：

12、等待车队在绿灯放行时的离开速度（单位：米/秒），则表示等待队列的离开率（单位:辆/小时）。从避免阻塞的角度考虑，我们希望在第相位禁止时期内在该路口滞留的最大等待队列长度（单位：辆）能在第相位的开通时间被全部放行。即：由此得出约束（4）：c. 目标函数：针对尽可能顺畅的目标，我们考虑用延误时间和滞留车辆数两个量作为路口顺畅程度的评价指标，由此引入延误率，滞留率及阻塞度的概念，现说明如下：定义3：延误率（delay）：表示一个周期时长内所有经过该路口的每一辆车的平均延误时间和周期时长的百分比，无量纲。延误率的计算如下：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。图4 等待车队示意图123k

13、TH·qVc停车线l第TH·q辆a) 计算第第相位时，第进口道上第辆车的延误时间：在如图4所示的等待队列中，第辆车延误时间=红灯时长-红灯消逝时间+等待前辆车离开的时间 +反应时间。其中，红灯消逝时间=，表示相邻两车到达路口的间隔时间，为等待车队中相邻两车的车头间距。为等待车队在绿灯时的离开速度。反应时间取典型值2秒。由此可得:b) 计算第相位时，第进口道上每辆车的平均延误时间其中，分子表示第进口道上等待队列的总延迟时间，分母表示在一个周期时长内经过第进口道的总车辆数，二者之比表示在一个周期内通过第进口道的每一辆车的平均延误时间。c) 延误率的计算：我们将第相位时，第进口道

14、上每辆车的平均延误时间在一个周期各个相位及相应进口道上取平均，得到了一个周期内通过该路口的每一辆车的平均延误时间：，再除以周期时长，得到了延误周期百分比，即d) 平均滞留车辆数的计算：表示路口处于第相位时，第进口道上的交通流量，即为第相位时，第进口道红灯期间滞留的车辆数，则在一个周期时长内平均每个进口道的滞留车辆数为：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。由上面的推算，我们得到了优化目标的量化表达式，在尽可能顺畅的要求下，极小化延误率和平均滞留车辆数，即： d. 不同周期时长单点红绿灯配时优化模型：至此，我们得出了单点红绿灯配时优化模型：其中, ，考虑到延误率和平均滞留车辆数的数

15、量级和量纲均不同，因此我们将它们分别归一化为两个无量纲的量并通过与基准值相比较消除数量级的影响，即：和，其中和均为基准值(通过一个特定的解)。最后，通过加权因子对其加权并求和，引出阻塞度的概念：阻塞度（block）：。尽可能顺畅，即阻塞度最小。5.1.2共用周期下的单点红绿灯配时设计我们已经确定出使相关指标达到最优下的各个路口各自的周期时长，在此基础上，结合实际情况和便于统一管理，我们需要从总体上对各路口的交通状况进行优化，确定出一个最佳共用周期时长。因此，我们建立出以为目标的优化模型，约束条件和单点红绿灯配时优化模型相似，不再赘述。方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。5.2

16、典型环路网交通信号联动方案的设计我们已经通过极小化延误率和平均滞留车辆数的方法来保证某个十字路口尽可能顺畅，并设计出周期时长和红绿灯配时方案，在此基础上，为方便城市交通系统的管理，提出了共用周期时长的概念，从总体尽可能顺畅的角度出发，制定出共用周期时长下的各个典型环路口的红绿配时方案。现在，我们需要从宏观上对整个典型环路上的交通信号灯进行联动控制，结合前面提到的方案设计内容，以保证整个典型环路交通尽可能顺畅，避免出现阻塞为目标，设计出在初始时刻典型环路各路口的相位状态分布及各相邻路口间的相位差，进而得到典型环路各十字路口交通灯运作时序图。针对上述问题，我们考虑以一个周期内从某路口流向下一路口并

17、无阻碍通过的车辆数尽可能多为目标，确定出相邻两路口的相位差。现分析如下：图5 典型环路多交叉口为方便叙述和读者理解，我们先提出相位差和离开率的概念，相位差：是指相邻路口同一相位绿灯或红灯起始时间之差。如果知道上一路口在某一时刻所处的相位状态，则下一路口的相位状态随之而确定。离开率：是指离开该路口并驶向特定的下一路口的车辆离开速率，单位为辆/小时。离开率的计算：考虑到一个周期内某车道到达的车全部在该相位运行时离开，则第一辆车与最后一辆车的离开时间间隔为，离开率（为到达第个进口道的车辆并驶向下一路口的车流量）。并且假设队列中各辆车速度与时间间隔在行驶过程中保持不变，则这一列车的流速即为离开率。因为

18、第一、二、三相位运行时均有车辆使向下一路口，第四相位没有，则在两路口间形成连续的具有三种时间间隔的车队。 由上分析可知车队并不是泊松流，而是分段函数。根据实地考察所采集的数据和按3.2中计算出的该路口红绿配时方案，我们可以做出以时间为横轴，以离开率为纵轴的图像，如图6所示：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。图6 时间与离开率的关系示意图离开率(q)时间一相位二相位三相位四相位滑动窗口S1S2t1t2令该路口起始车即进口道1右转的第一辆车到达下一路口的时刻为0，则该路口第三相位的开始时刻为(L为两路口间距，为车辆行驶速度，取)。显然，离开率时间积在几何意义上表示该时间段上离开率

19、和时间轴所围的面积，在实际意义上表示该时间段上离开第个路口驶向第个路口的车辆数。根据我们前面提到的建模目标，要使这些车辆尽可能多的通过第个路口，第+1个路口必须处于第三相位来让该方向上的车辆无阻碍通过，且数量尽可能多。因此，第个路口的第三相位开启状态可用图7中的矩形表示：该矩形在时间轴上滑动表示第三相位在第个路口的任意一种分布状态，类似于一个长度变动(长度的具体分析见下)的滑动窗口。我们考虑的是无阻碍通过的车辆数尽可能多，而红灯时间内排队的车辆会使有效绿灯时间流失(即一个周期内保证车辆无阻碍通过的时间)，其值应是第三相位的运行时间减去排队车辆用去的时间，即，其中，2表示其前面的每一辆车离开所需

20、要的平均时间为2秒，表示在等待队列中的车辆数，也即覆盖区外的面积，如图7所示。从而该路口第三相位的开始时刻为，由此可得相位差为最优解即有效绿时所覆盖的时间区域上离开率所围面积(s)最大。下面建立出基于该方法的模型：局部变量说明：第一辆无阻碍通过下一路口的时刻：最后一辆无阻碍通过下一路口的时刻：第k个路口第个相位时间，也即绿灯时间（若是则表示第个相位时间）：第个路口的车辆平均离开率，单位辆/小时方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。结合图7，针对下一路口的第三相位时间，分类讨论如下：显然，有下面约束成立 1) 当，时，当，时， 当，时，：2) 当时：，当：，3) 当，时：综合以上可

21、得此优化模型为：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。六模型求解6.1模型一求解： 模型一中，两个目标函数都是比较复杂的非线性函数，约束条件中也因含有max等函数成为非线性，而变量约束为整数，故整个模型一是非线性规划，并且利用计算复杂性的理论，我们可以大致确定其模型属于NP难问题。对于NP难问题，现在还没有多项式算法，不能在较短的时间内得出最优解(通过将优化模型输入LINGO，LINGO不能找出全局最优解)，故我们选择采用启发式的禁忌搜索算法，得出了较好的结果6.1.1禁忌搜索算法介绍：禁忌搜索是一种全局逐步寻优的人工智能算法，它常能有效的应用于一些典型NP问题，如TSP。禁忌搜

22、索能够同时拥有高效性和鲁棒性。禁忌搜索算法设计忌搜索是对局部领域搜索的扩展。传统局部邻域搜索是基于贪婪思想在当前解的邻域中进行搜索，搜索性能完全依赖于邻域结构和初始解的选取，搜索结果容易陷入局部极小而无法保证全局最优。而禁忌搜索从一个初始可行解s开始，通过变换得解的邻域函数V(s)，按照某种选择策略从中选取一个解best，从s移动到best，把best作为一个新解，重新叠代搜索，直到满足退出机制。为避免循环和陷入局部最优，禁忌搜索使用禁忌表记录已经到达的局部最优点，也即最近进行的移动状态。在下一步的搜索中利用规定的禁忌规则，在一定搜索次数内不允许选择这些被禁的搜索点，从而可以跳出局部最优的限制

23、。6.1.2禁忌搜索算法设计：初始解：初始解的选取直接影响最终结果的好坏，初始解选取得好则有利于进行全局搜索，避免陷入局部最优。通过对模型进行分析，我们选取周期时长的下限同变量约束进行规划，得出在此条件下，以此作为以后T变动时的初始解，因为直观上，对某一个T值，它是最优的。邻域的选择：领域的选取是通过对解进行一定的变化调整进行的。一般而言，变化形式分为解的直接变化、向量分量的变化和目标值的变化。鉴于解的数据结构，采用解的直接变化形式，即选取一个步长h，对当前解做进行调整。这样领域中满足约束的解最多共有24个。对h，为了防止陷入局部最优，我们采取双步长的方式，即。目标函数：在禁忌搜索中，目标函数

24、一般作为选取领域中的解的规则。对于模型一，我们选取通过与基准值比较并已加权求和的阻塞度作为目标值，即：，和为加权系数，在求解时考虑到前者的在主观上是更重要的因素，我们去，。候选集：候选集用于储存当前解的领域，由于最多有24个，故候选集设置为矩阵，其中前4列储存解的4个分量，最后一列储存对应的。禁忌表的确定：，禁忌对象的选择通常也有三种形式：解的直接变化、分量对换的变化和目标值的变化。由于分量对换在当前的数据结构和领域选择上难以实现，而目标值变化的对象过多，难以得到全局最优，故我们选择解的直接变化，但只取其中已经在领域选取的迭代中出现过的解。方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。算

25、法终止规则：在算法迭代中，当前的领域最优解与已经得到的全局最优解之间的迭代次数之差不能超过一定值，此值我们取500。6.1.3禁忌搜索算法的伪代码：初始化：; ; %流量矩阵，初始解; %设置禁忌表，实为栈 %从领域中选出的最优解 %当前解、最终解 %当前解迭代步数、最优解迭代步数; %候选集，前4列存储解，最后一列储存目标值开始： 当 %当目标值没变化的迭代步数不太多 %更新迭代步数生成的候选集;在中选择使目标函数最小的;将后4行依次上移一行将加入最后一行，第一行溢出 %更新禁忌表(栈) 若 %从领域中选出的最优解比当前解更优 %更新全局变量=; 继续6.1.4候选集生成函数的伪代码：初始化

26、：； %步长 ; %输入当前解 %选择矩阵124循环： ; %行下标 %生成第一个步长下的候选集方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。; %生成第二个步长下的候选集判断中前4列的解是否满足约束并且是非禁忌的;若不是，则将此行设为0；若是，则对此行：; 继续6.1.5模型求解结果： 表1 第一交叉口与第二交叉口红绿灯方案及其效果表：相位编号1234相位时间(秒)32.721.832.721.8最佳周期(秒)109总延误百分比0.4442平均排队长(辆)4.618表2 第一交叉口与第三交叉口红绿灯方案及其效果表：相位编号1234相位时间(秒)34263931最佳周期(秒)130总延误

27、百分比0.4927平均排队长(辆)4.214表3 第一交叉口与第四交叉口红绿灯方案及其效果表：相位编号1234相位时间(秒)22.425.618.421.6最佳周期(秒)88总延误百分比0.4156平均排队长(辆)3.108对不同的交叉口通过以上的计算对于最佳周期、总延误百分比、平均排队长三个参数的变化情况如下图所示：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。图7 最佳周期、总延误百分比、平均排队长三个参数的变化情况表4 各路口公共周期及其红绿灯相位表： 相位支路1234总延误百分比平均等待队长第二交叉口35.123.435.128.40.42534.96第三交叉口30.123.43

28、5.128.40.49833.8第四交叉口35.128.430.123.40.39524.17公共周期1206.2模型二求解：考虑到与分别是、的分段函数：鉴于分段函数属于非线性，难以处理，我们考虑再引入变量将其化为线性模型。记轴上的分点为：。当属于第1个小分区时，记。因为在上是线性的，所以方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。同样，当属于第2个小分区时，。最后当属于第3个小分区时，。为了表示在哪个小区间，引入0-1变量，当在第个小区间时，否则。这样，应满足：此时和可以统一地表示为：，同理，分析与，可表示为：，其中，这样满足：综合上述，通过引入新的变量，模型二可化为以下混合线性规划

29、模型：方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。将以上模型输入到LINGO中求得：，利用相位差计算公式的s。带入各项数据，得到一下各个节点的相位运行时序图，如图9。图8 各节点的相位分布时序图11234234123341七结果分析和检验通过该模型可求解出典型环路各路口的红绿灯配置方案。仅考虑第一交叉口与第三交叉口时，其红绿灯方案及其效果如下表所示：相位编号1234相位时间(秒)34263931最佳周期(秒)130总延误百分比0.4927平均排队长(辆)4.214从上面的结果可以看出平均每个车道在一周期内的排队长仅为4辆，比实际在该路口观测到的11辆小很多，总延时百分比也相对较小，从而

30、所明该红绿灯配置具有一定的可行性。将典型环路上所有节点统一考虑，得出各路口红绿灯方案。对于单个路口而言，阻塞度有略微上升，周期时长为120s。对于第一交叉口与第三交叉口，各相位时间分别为30.1，23.4，35.1，28.4，相比于上表中数据，各相位时间均有所减少，这主要是由邻近路口的流量造成的。总延误百分比有所下降，但排队长有所增加，总的阻塞度有略微上升。但还是在单点最优的附近，表明这样得出的结果较好。方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。通过模型二的求解，得出第一交叉口与第二交叉口与第一交叉口与第三交叉口的第三相位差为64.5s。且通过系统模拟表明这一相位差具有优越性(具体模

31、拟结果见可行性报告)。八模型评价与改进本文在成都市典型环路交通为背景，以四相位，二车道的十字路口为研究对象，建立了单点红绿配时优化模型和典型环路网交通信号联动管理模型。优点：1) 适用范围：单点红绿配时优化模型中所采用阻塞度概念和相应的评价指标及约束条件同样可以应用于任意的i相位j车道十字路口，从而得到一般情形下的单点红绿配时方案，因此该模型具有一定的普适性和通用性。而典型环路网交通信号联动管理模型中所采用的图解法，分别将几何，代数和实际意义上的面积，离开率时间积，离开车辆数结合起来，建立出信号联动管理模型，虽然思想简单却较好的解决了问题，这正是数学建模所提倡的，这在一定程度上也反映了模型的普

32、适性。2) 可行性：通过系统模拟的方法，得出的平均每辆车的延误时间，一个周期内的最大排队队长两个指标均与实际相符，说明了单点红绿配时模型的可行性（详见可行性论证报告部分）。缺点：在建模的过程中，为简化问题，我们没有考虑黄灯的时长，这对设计出的优化方案造成了一定的误差。改进：在各个路口采用各自的最优周期时长，以相应的红绿灯优化配时方案运作。十参考文献1 李尚志.数学建模竞赛教程.江苏教育出版社.19962 刘建军：交通工程学，人民交通出版社2001年版；3 吴兵，李晖等：交通管理与控制，人民交通出版社2000年版；4 谢金星，优化建模与LINDO/LINGO软件，北京：清华大学出版社，20055

33、 寿纪麟.数学建模方法与案例.西安交通大学出版社.19936 刘承平.数学建模方法.高等教育出版社.2002第一版7 李建斌, 高成修, 城市道路网络多交叉路口交通信号实时优化控制模型与算法, 系统工程, 2004 年10 期8 熊烈强，王富，李杰，彭少民，混合交通流参数之间的关系，华中科技大学学报(自然科学版)，2005年7月，第33卷第7期方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。9 王学堂城市交通信号整体优化理论的概念与方法 J 长安大学学报： 自然科学版，2 0 0 2,2 2 ( 5 ) , 9 2 9 4 10 张伯生， 任 翔 城市单交叉路 口多相位智能控制方法与计算机

34、仿真 J 系统工程, 2 0 0 1,( 6 ) ,5 2 5 6 11 王 炜，过秀成 交通工程学 M .南京:东南大学出版社， 十一附录附录1：原始流量数据：第一交叉口与第二交叉口车流表路口j相位i12341160-77-2380-30-3-330-4384-175-55第一交叉口与第三交叉口车流表路口j相位i1234160-30-240-300-3-317-3234-20-180第一交叉口与第四交叉口车流表路口j相位i12341196-176-2300-80-3-89-1244-200-200方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。附录2：计算单个节点的最佳周期和配时方案的禁

35、忌搜索算法的MATLAB实现%以第一交叉口与第三交叉口车流表为例：q=60 0 30 0;40 0 300 0;0 317 0 323;0 20 0 180/3600;global D0 L0 %全局变量 作为基准D0 L0=ami(60 90 60 90,q,100);%计算基准目标函数的值best_x=zeros(1,4); %记录最终解best_T=0; %记录最终周期时长best_fx=10;%记录最佳的目标函数评价值for T=50:5:150 %周期时长 从50开始循环 best=opt(q,T); if fx(best,q,T)<best_fx %在此T值下的最优目标评价值

36、 best_fx=fx(best,q,T); best_T=T; best_x=best; %更新全局变量 endenddisp('最佳周期时长为:')best_Tdisp('此时各个相位的绿灯时间为:')TL=best_T-best_xdisp('此时的总平均延误时间占周期的百分比和排队队长为:')D L=ami(best_x,q,best_T)各子程序如下：1.1 function best=opt(q,T)%某个节点在特定周期T和流量q下的最优方案%best为各个相位的红灯时间x=7*T/10 4*T/5 7*T/10 4*T/5; %初始

37、解taboo=zeros(5,4); %设置禁忌表 数据结构为栈s_best=x; s_now=x; best=x; %从候选集中选出的最优解、当前解、最终解k=0; best_k=0; %当前的迭代步数、最优解所在迭代步数while k-best_k<200 %当迭代步数超过最大允许迭代步数时停止运 k=k+1; %迭代步数增加 V=proV(s_now,q,T,taboo); %生成当前解的候选集 if V=0 s_best=choose(V); %从生成的候选集中选出目标函数评价值最小的解 end taboo=renew_taboo(s_now,taboo); %更新禁忌表 if

38、fx(s_best,q,T)<fx(best,q,T) %当从候选集中选出的最优解比当前最优解更优方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。 best=s_best; best_k=k; end s_now=s_best;end1.2function D L=ami(x,q,T)%计算目标函数的表达式的值%x是解向量 表示各相位的红灯时间%q是此节点的流量的矩阵%T是周期时长%D是对每辆车而言的总平均延误时间%L是单位时间内总等候车队的平均车长D=0;L=0;for i=1:4 for j=1:4 if q(i,j)>0 D=0.5*x(i)2+0.5*x(i)/q(i,

39、j)+x(i)2*q(i,j)+D; else continue end endendD=D/T/T/8;for i=1:4 for j=1:4 L=x(i)*q(i,j)+L; endendL=L/8;1.3function F=fx(x,q,T)%计算归一化后的目标函数评价值D L=ami(x,q,T);global D0 L0 %D0 L0是全局常量f1=0.85;f2=0.15; %目标函数的加权系数F=f1*(D-D0)/D0+f2*(L-L0)/L0;1.4function bool=iscon(x,q,T)%验证输入的解向量x是否满足约束条件%q是此节点的流量的矩阵方案范文无法思

40、考和涵盖全面，最好仔细浏览后下载使用。%T是周期时长bool=1;if x(1)>T-15|x(3)>T-15|x(1)<35|x(3)<35 bool=0;endif x(2)>T-10|x(4)>T-10|x(2)<40|x(4)<40 bool=0;endfor i=1:4 if max(q(i,:)*x(i)>(T-x(i)/2 bool=0; endend1.5function b=istaboo(x,taboo)%验证输入的解向量x是否被禁%q是此节点的流量的矩阵%T是周期时长b=1;for i=1:5 if taboo(i,:

41、)=x b=0; endend1.6function V=proV(s,q,T,taboo)%生成当前解的候选集V=zeros(24,5); %设置候选集 前四列储存解 后一列储存对应的目标值h=1,fix(T.5/2); %步长A=1 -1 0 0;1 0 -1 0;1 0 0 -1;0 1 -1 0;0 1 0 -1;0 0 1 -1;-1 1 0 0;-1 0 1 0;-1 0 0 1;0 -1 1 0;0 -1 0 1;0 0 -1 1; %选择矩阵for i=1:12 V(i,1:4)=s+A(i,:)*h(1); V(i+12,1:4)=s+A(i,:)*h(2);endfor i

42、=1:24 temp1=iscon(V(i,1:4),q,T); temp2=istaboo(V(i,1:4),taboo); if temp1=1&&temp2=1 %当此行的解满足约束且没有被禁忌 V(i,5)=fx(V(i,1:4),q,T); %算出此解的目标函数评价值 else V(i,:)=0; %若不满足约束或者被禁忌 则将此行全部设为0方案范文无法思考和涵盖全面，最好仔细浏览后下载使用。 endend1.7function tab=renew_taboo(s,taboo)%更新禁忌表%将taboo后四行依次上移一行 将s加入taboo最后一行 taboo第一行溢出if taboo(5,:)=s