CN114677847A - 一种基于迭代优化的区域绿波协调方案求解方法 - Google Patents

Abstract

本发明公开了一种基于迭代优化的区域绿波协调方案求解方法，包括：选取第一条东西走向协调干道，计算不同公共信号周期C取值下协调干道上所有乙类交叉口不同可选相位相序组合方案对应的甲类交叉口最佳相位相序组合方案；选取第一条南北走向协调干道，计算不同公共信号周期C取值下协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案；确定当前协调干道，计算不同公共信号周期C取值下协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案，直至均为已标定协调干道；利用记录矩阵依次反推，获得最优区域绿波协调控制方案，获得理想的绿波协调控制效果。

Description

一种基于迭代优化的区域绿波协调方案求解方法

技术领域

本发明涉及交通信号控制技术领域，更具体的，涉及一种基于迭代优化的区域绿波协调方案求解方法。

背景技术

随着城市路网密度及支路通行能力的不断提升，干道协调控制已然满足不了城市交通的发展需求，故城市区域交通信号优化研究工作逐渐引起了各界的高度关注。现有的区域绿波协调控制设计通常是从模型构建与优化的角度出发，以干道绿波协调控制模型为基础，采用模型法求解最优协调控制方案。然而在模型构建过程中往往由于存在相位差闭环约束，导致区域绿波协调控制模型难以获取可行解或需要牺牲部分路段的绿波带宽，并存在模型复杂度较高、运算效率较低、以及难以确保在有限时间内获得有效设计方案等问题。数解法，如卢凯等在《非对称通行条件下的双向绿波协调控制数解算法.中国公路学报(06)，99-107.》中提出的算法，是通过解析绿波协调控制参数之间的数学关系，利用数值计算实现绿波协调控制方案的快速求解，具备可操作性强、计算量较小的优点，但目前鲜见拓展应用至包含闭环协调干道的复杂路网。

发明内容

本发明的主要目的在于克服现有区域绿波协调方案设计方法求解速度较慢、较难获取可行解或需要牺牲部分路段绿波带宽的缺点与不足，将干道数解法的求解思路拓展应用至包含闭环协调干道的复杂路网，通过提供一种基于迭代优化的区域绿波协调方案求解方法，实现对复杂路网交叉口公共信号周期、信号相位相序以及相位差的快速综合优化，使得协调区域内各协调干道均能获得较为理想的绿波协调控制效果。

为了实现本发明目的，本发明提供的一种基于迭代优化的区域绿波协调方案求解方法，协调区域内存在n条东西走向与m条南北走向构成封闭路网的相交协调干道；依次由北往南将第i条东西走向协调干道命名为A_EWi，定义由西往东为协调正向，由东往西为协调反向；由西往东将第j条南北走向的协调干道命名为A_SNj，定义由北往南为协调正向，由南往北为协调反向；若经过某交叉口的干道中仅存在一条协调干道，则定义该交叉口为非节点交叉口，若经过某交叉口的干道存在两条或两条以上协调干道，则定义该交叉口为节点交叉口；将东西走向协调干道A_EWi与南北走向协调干道A_SNj的节点交叉口记为I_(i，j)，相位相序设置方式记为P_(i，j)；将节点交叉口I_(1，1)西进口直行相位的绿灯起始时刻点作为协调区域的相位差零点，包括以下步骤：

S1、选取第一条东西走向协调干道A_EW1，计算不同公共信号周期C取值下，协调干道上所有乙类交叉口不同可选相位相序组合方案对应的甲类交叉口最佳相位相序组合方案，其中，若经过某交叉口的所有协调干道均为已标定协调干道或当前协调干道，则定义该交叉口为甲类交叉口；若经过某交叉口的所有协调干道中存在但不全为已标定协调干道或当前协调干道，则定义该交叉口为乙类交叉口；若经过某交叉口的所有协调干道中不存在已标定协调干道或当前协调干道，则定义该交叉口为丙类交叉口；

S2、选取第一条南北走向协调干道A_SN1，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案及对应的最佳相位相序组合方案计入相应的记录矩阵；

S3、选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案以及对应的最佳相位相序组合方案计入相应的记录矩阵，重复该步骤直至东西走向协调干道A_EWn与南北走向协调干道A_SNm均为已标定协调干道；

S4、利用记录矩阵依次反推获得在上一协调干道组合计算时所确定的甲类交叉口最佳相位相序组合方案，最终获得最优区域绿波协调控制方案。

与现有的技术相比，本发明的有益效果如下：

1)本发明给出了一套完整的区域绿波协调方案迭代优化求解步骤，能够实现公共信号周期、信号相位相序以及相位差的全局综合优化，可适用于不同绿波协调控制需求下的区域绿波协调设计方案求解。

2)本发明利用记录矩阵完成协调干道的逐步迭代组合计算与最优方案的逐步反推，利用步骤化的求解思路降低了求解算法的复杂度，解决了交通网络信号优化求解所面临的维数灾难问题，实现了区域绿波协调方案的快速求解。

3)本发明是干道数解法在区域绿波协调控制设计中的进一步拓展，通过对非基准节点交叉口相位差确定准则的优化设计，克服了干道数解法在拓展应用到区域闭环网络时所面临的相位差计算冲突问题，实现了东西走向与南北走向干道交叉的协同优化设计。

4)本发明给出了一种基于迭代优化的区域绿波协调方案求解方法，利用逐步迭代优化的思想，将干道数解法的求解思路拓展应用至区域协调路网，给出一种新型的区域绿波协调方案求解方法，可以提高区域绿波协调控制方案设计的求解效率，并获得更为理想的区域绿波协调控制方案。

附图说明

图1是本发明实施例提供的一种基于迭代优化的区域绿波协调方案求解方法流程图；

图2是本发明东西走向协调干道A_EW1甲类交叉口与乙类交叉口分类情况示意图；

图3是本发明节点交叉口东西走向相位差示意图；

图4是本发明南北走向协调干道A_SN1甲类交叉口与乙类交叉口分类情况示意图；

图5是本发明节点交叉口南北走向相位差示意图；

图6是本发明南北走向协调干道A_SN2为当前协调干道时甲类交叉口与乙类交叉口分类情况示意图；

图7是本发明东西走向协调干道A_EWn为终止协调干道时甲类交叉口与乙类交叉口分类情况示意图；

图8是本发明节点交叉口I_(i，j)的相位差O_(i，j)示意图；

图9是本发明实施例中的协调区域路网结构示意图；

图10是本发明实施例中东西走向协调干道A_EW1组合计算示意图；

图11是本发明实施例中南北走向协调干道A_SN1组合计算示意图；

图12是本发明实施例中南北走向协调干道A_SN2组合计算示意图；

图13是本发明实施例中东西走向协调干道A_EW2组合计算示意图；

图14是本发明实施例中南北走向协调干道A_SN3组合计算示意图；

图15是本发明实施例中东西走向协调干道A_EW3组合计算示意图；

图16是本发明实施例中南北走向协调干道A_SN4组合计算示意图；

图17是本发明实施例中计算所得最佳区域绿波协调方案示意图；

图18是本发明实施例中计算所得最佳区域绿波协调方案对应的空间时距图。

具体实施方式

下面结合附图和具体实施方式对本发明做详细描述。

在本发明的其中一些实施例中，如图9所示为一个三横四纵的不规则路网，把协调区域内由北往南的三条东西走向的协调干道依次命名为A_EW1、A_EW2和A_EW3；协调区域内由西往东的四条南北走向协调干道依次命名为A_SN1、A_SN2、A_SN3和A_SN4。依次将东西走向协调干道A_EWi与南北走向协调干道A_SNj的节点交叉口记为I_(i，j)，其中节点交叉口I_(1，1)、节点交叉口I_(1，4)与节点交叉口I_(2，1)均为T型交叉口。路网所有路段的绿波设计行驶速度均为12.5m/s，公共信号周期的取值范围为[100，150]s，各交叉口均处于未饱和交通状态，基于现状交通流量数据进行进口单独放行相位设计下的交叉口绿信比分配设计如表1所示。

表1各交叉口绿信比分配表

图1为本发明提供的一种基于迭代优化的区域绿波协调方案求解方法，包括步骤以下：

步骤S1、选取第一条东西走向协调干道A_EW1，计算不同公共信号周期C取值下，协调干道上所有乙类交叉口不同可选相位相序组合方案对应的甲类交叉口最佳相位相序组合方案，其中，若经过某交叉口的所有协调干道均为已标定协调干道或当前协调干道，则定义该交叉口为甲类交叉口；若经过某交叉口的所有协调干道中存在但不全为已标定协调干道或当前协调干道，则定义该交叉口为乙类交叉口；若经过某交叉口的所有协调干道中不存在已标定协调干道或当前协调干道，则定义该交叉口为丙类交叉口；

S101、选取第一条东西走向协调干道A_EW1为当前协调干道，确定当前协调干道A_EW1上甲类交叉口以及乙类交叉口分类情况；

在本发明的其中一些实施例中，当前协调干道A_EW1上没有非节点交叉口，故不存在甲类交叉口，当前协调干道A_EW1上的所有节点交叉口均为乙类交叉口，如图10所示；

S102、以节点交叉口I_(1，1)为基准交叉口，计算在不同公共信号周期C的取值下，当前协调干道A_EW1上各交叉口的相位差取值；

在本发明的其中一些实施例中，协调区域内各交叉口的公共信号周期C的优化步长取为1s，优化范围记为[C_min，C_max]，其中，C_min为允许最小公共信号周期，为100s，C_max为允许最大公共信号周期，为150s；根据节点交叉口I_(1，1)，计算当前协调干道A_EW1上的任意节点交叉口I_(1，j)(2≤j≤4)的相位差O_(1，j)，如图3所示，满足式(1)：

式中，T_GW(1，1)和T_GW(1，j)分别为节点交叉口I_(1，1)和节点交叉口I_(1，j)的西进口绿灯时长；v_(1，j)为节点交叉口I_(1，1)和节点交叉口I_(1，j)间的车辆平均行驶速度；a_(1，j)为节点交叉口I_(1，1)和节点交叉口I_(1，j)间的交叉口理想间距，满足式(2)：

式中，Δt_W→E(1，1)和Δt_W→E(1，j)分别为节点交叉口I_(1，1)和节点交叉口I_(1，j)西进口超前于东进口直行相位绿灯中心时刻点的时长；n_(1，j)为整数，其取值应使得理想间距a_(1，j)与实际间距s_(1，j)最为接近；

将东西走向协调干道A_EWi的非节点交叉口判定系数定义为0-1变量δ_EWi，若东西走向协调干道A_EWi上存在非节点交叉口，则δ_EWi等于1，东西走向协调干道A_EWi的非节点交叉口的数量记为p_i，由西往东将第k个非节点交叉口记为I_s(i，k)，相位相序设置方式记为P_s(i，k)，若东西走向协调干道A_EWi上不存在非节点交叉口，则δ_EWi等于0；

根据节点交叉口I_(1，1)，计算当前协调干道A_EW1上的任意非节点交叉口I_s(1，k)(1≤k≤p₁)的相位差O_s(1，k)，满足式(3)：

式中，T_GWs(1，k)为非节点交叉口I_s(1，k)西进口绿灯时长；v_s(1，k)为节点交叉口I_(1，1)与非节点交叉口I_s(1，k)间的车辆平均行驶速度；a_s(1，k)为节点交叉口I_(1，1)和非节点交叉口I_s(1，k)间的理想间距，满足式(4)：

式中，Δt_{W→Es(1，k)}为非节点交叉口I_s(1，k)西进口超前于东进口直行相位绿灯中心时刻点时长；n_s(1，k)为整数，其取值应使得理想间距a_s(1，k)与实际间距s_s(1，k)最为接近；

在本发明的其中一些实施例中，当前协调干道A_EW1上不存在非节点交叉口，即δ_EW1＝0，则无需计算当前协调干道A_EW1上的非节点交叉口相位差；

S103、计算并存储不同公共信号周期C取值以及所有乙类交叉口不同可选相位相序组合方案[P_(1，1)，P_(1，2)，P_(1，3)，P_(1，4)]对应的当前协调干道A_EW1甲类交叉口最佳相位相序组合方案

其中，所述最佳相位相序组合方案

可以使得当前协调干道A_EW1干道偏移绿信比Δλ_EW1的绝对值|Δλ_EW1|取得最小值，将所述可选相位相序组合方案[P_(1，1)，P_(1，2)，P_(1，2，P_(1，4)]和所述最佳相位相序组合方案

计入记录矩阵M_(C，1)，

此时东西走向协调干道A_EW1转为已标定协调干道；

东西走向协调干道A_EW1的干道偏移绿信比Δλ_EW1满足式(5)：

Δλ_EW1＝|Δλ_W→E1|+|Δλ_E→W1| (5)

其中，Δλ_W→E1为东西走向协调干道A_EW1的协调正向偏移绿信比，Δλ_E→W1为协调反向偏移绿信比，分别满足式(6)和式(7)：

其中，任意节点交叉口I_(1，j)东西向的正向偏移绿信比Δλ_W→E(1，j)和反向偏移绿信比Δλ_E→W(1，j)分别满足式(8)与式(9)：

式中，ΔT_(1，j)为节点交叉口I_(1，1)和节点交叉口I_(1，j)理想间距与实际间距的差值D_(1，j)与车辆平均行驶速度v_(1，j)的比值；

任意非节点交叉口I_s(1，k)东西向的正向偏移绿信比Δλ_{sW→E(1，k)}和反向偏移绿信比Δλ_{sE→W(1，k)}分别满足式(10)与式(11)：

式中，ΔT_s(1，k)为节点交叉口I_(1，1)与非节点交叉口I_s(1，k)理想间距与实际间距的差值D_s(1，k)与车辆平均行驶速度v_s(1，k)的比值；

由于在本实施例中当前协调干道A_EW1上不存在甲类交叉口，如图10所示，则记录矩阵M_(C，1)不存在。

步骤S2、选取第一条南北走向协调干道A_SN1，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案及对应的最佳相位相序组合方案计入相应的记录矩阵；

S201、选取第一条南北走向协调干道A_SN1为当前协调干道，确定当前协调干道A_SN1上甲类交叉口以及乙类交叉口分类情况；

在本发明的其中一些实施例中，如图11所示，当前协调干道A_SN1与已标定协调干道A_EW1的节点交叉口I_(1，1)以及当前协调干道A_SN1上的所有非节点交叉口均为甲类交叉口，当前协调干道A_SN1上的其余节点交叉口均为乙类交叉口；

S202、以节点交叉口I_(1，1)为基准交叉口，计算在不同公共信号周期C的取值下，当前协调干道A_SN1其余交叉口的相位差取值；

根据节点交叉口I_(1，1)，计算当前协调干道A_SN1上的任意节点交叉口I_(i，1)(2≤i≤3)的相位差O_(i，1)，如图5所示，满足式(12)：

式中，T_GW(i，1)为节点交叉口I_(i，1)西进口绿灯时长；v_(i，1)为节点交叉口I_(1，1)和节点交叉口I_(i，1)间的车辆平均行驶速度；Δt_W→N(1，1)与Δt_W→N(i，1)分别为节点交叉口I_(1，1)和节点交叉口I_(i，1)西进口超前于北进口直行相位绿灯中心时刻点时长；a_(i，1)为节点交叉口I_(1，1)和节点交叉口I_(i，1)间的理想间距，满足式(13)：

式中，Δt_N→S(1，1)和Δt_N→S(i，1)分别为节点交叉口I_(1，1)和节点交叉口I_(i，1)北进口超前于南进口直行相位绿灯中心时刻点时长；n_(i，1)为整数，其取值应使得理想间距a_(i，1)与实际间距s_(i，1)最为接近；

将南北走向协调干道A_SNj的非节点交叉口判定系数定义为0-1变量δ_SNj，若南北走向协调干道A_SNj上存在非节点交叉口，则δ_SNj等于1，南北走向协调干道A_SNj的非节点交叉口的数量记为q_j，由北往南将第l个非节点交叉口记为I_a(j，l)，相位相序设置方式记为P_a(j，l)，若南北走向协调干道A_SNj上不存在非节点交叉口，则δ_SNj等于0；

根据节点交叉口I_(1，1)，计算当前协调干道A_SN1上的任意非节点交叉口I_a(1，l)(1≤l≤q₁)的相位差O_a(1，l)，满足式(14)：

式中，T_GWa(1，l)为非节点交叉口I_a(1，l)西进口绿灯时长；v_a(1，l)为节点交叉口I_(1，j)与非节点交叉口I_a(1，l)间的车辆平均行驶速度；Δt_{W→Na(1，l)}为非节点交叉口I_a(1，l)西进口超前于北进口直行相位绿灯中心时刻点时长；a_a(1，l)为节点交叉口I_(1，1)与非节点交叉口I_a(1，l)间的理想间距，满足式(15)：

式中，Δt_{N→Sa(1，l)}为非节点交叉口I_a(1，l)北进口超前于南进口直行相位绿灯中心时刻点时长；n_a(1，l)为整数，其取值应使得理想间距a_a(1，l)与实际间距s_a(1，l)最为接近；

在本发明的其中一些实施例中，由于在本实施例中当前协调干道A_SN1上不存在非节点交叉口，δ_SN1＝0，则无需计算当前协调干道A_SN1上的非节点交叉口的相位差；

S203、计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1，2)，P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)]对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案[P_(1，1)]，其中，所述最佳相位相序组合方案[P_(1，1)]可以使得当前协调干道A_SN1干道偏移绿信比Δλ_SN1与已标定协调干道A_EW1干道偏移绿信比Δλ_EW1的绝对值之和|Δλ_EW1|+|Δλ_SN1|取得最小值，将所述可选相位相序组合方案[P_(1，2)，P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)]和所述最佳相位相序组合方案[P_(1，1)]计入记录矩阵M_(C，2)，M_(C，2)＝|[P_(1，2)，P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)，P_(1，1)]，此时南北走向协调干道A_SN1转为已标定协调干道；

南北走向协调干道A_SN1的干道偏移绿信比Δλ_SN1满足式(16)：

Δλ_SN1＝|Δλ_N→S1|+|Δλ_S→N1| (16)

其中，南北走向协调干道A_SN1的协调正向偏移绿信比Δλ_N→S1与协调反向偏移绿信比Δλ_S→N1分别满足式(17)和式(18)：

其中，任意节点交叉口I_(i，1)南北向的正向偏移绿信比Δλ_N→S(i，1)和反向偏移绿信比Δλ_S→N(i，1)分别满足式(19)与式(20)：

式中，ΔT_(i，1)为节点交叉口I_(1，1)和节点交叉口I_(i，1)理想间距与实际间距的差值D_(i，1)与车辆平均行驶速度v_(i，1)的比值；

任意非节点交叉口I_a(1，l)南北向的正向偏移绿信比Δλ_{aN→S(1，l)}和反向偏移绿信比Δλ_{aS→N(1，l)}分别满足式(21)和式(22)：

式中，ΔT_a(1，l)为节点交叉口I_(1，1)与非节点交叉口I_a(1，l)理想间距与实际间距的差值D_a(j，l)与车辆平均行驶速度v_a(1，l)的比值。

步骤S3、选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案以及对应的最佳相位相序组合方案计入相应的记录矩阵，重复该步骤直至东西走向协调干道A_EWn与南北走向协调干道A_SNm均为已标定协调干道；

S301、根据已标定协调干道情况，选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道；

其中，确认东西走向协调干道A_EWi和南北走向协调干道A_SNj均为已标定协调干道，东西走向协调干道A_EWi+1和南北走向协调干道A_SNj+1均为未标定协调干道，则选取东西走向协调干道A_EWi+1作为当前协调干道将有j个节点交叉口转化为新增甲类交叉口，新增乙类交叉口m-j个；而选取南北走向协调干道A_SNj+1作为当前协调干道将有i个节点交叉口转化为新增甲类交叉口，新增乙类交叉口n-i个，则若满足式(23)，当前协调干道为东西走向协调干道A_EWi+1，否则为南北走向协调干道A_SNj+1：

n-i≥m-j (23)

在本发明其中一个实施例中，此时东西走向协调干道A_EW1和南北走向协调干道A_SN1均为已标定协调干道，东西走向协调干道A_EW2和南北走向协调干道A_SN2均为未标定协调干道，则选取东西走向协调干道A_EW2作为当前协调干道将新增3个乙类交叉口；而选取南北走向协调干道A_SN2作为当前协调干道将新增2个乙类交叉口；因此，当前协调干道为南北走向协调干道A_SN2；

S302、确定当前协调干道上甲类交叉口以及乙类交叉口分类情况；

当前协调干道与所有已标定协调干道的节点交叉口以及当前协调干道上的所有非节点交叉口均为甲类交叉口，当前协调干道上的其余节点交叉口均为乙类交叉口，如图12所示；

S303、计算在不同公共信号周期C的取值下，当前协调干道上新增乙类交叉口与非节点交叉口的相位差取值；

计算当前协调干道上的任意新增乙类交叉口I_(i，j)(2≤i≤n，2≤j≤m)的相位差O_(i，j)，将节点交叉口I_(i，1)与新增乙类交叉口I_(i，j)间的理想间距与实际间距的差值D_ew(i，j)与车辆平均行驶速度v_ew(i，j)的比值记为ΔT_ew(i，j)；将节点交叉口I_(1，j)与新增乙类交叉口I_(i，j)间的理想间距与实际间距的差值D_sn(i，j)与车辆平均行驶速度v_sn(i，j)的比值记为ΔT_sn(i，j)；为了同时保证东西方向与南北方向的绿波带宽，新增乙类交叉口I_(i，j)的相位差O_(i，j)应该在东西方向理想相位差O_ew(i，j)与南北方向理想相位差O_sn(i，j)之间取值，相位差O_(i，j)的取值范围为[min(O_ew(i，j)，O_sn(i，j))，max(O_ew(i，j)，O_sn(i，j))]，将东西方向理想相位差O_ew(i，j)与南北方向理想相位差O_sn(i，j)之差的绝对值记为ΔO_(i，j)＝abs(O_ew(i，j)-O_sn(i，j))，如图8所示，其中，a_ew(i，j)与s_ew(i，j)分别为节点交叉口I_(i，1)与新增乙类交叉口I_(i，j)间的理想间距与实际间距；a_sn(i，j)与s_sn(i，j)分别为节点交叉口I_(1，j)与新增乙类交叉口I_(i，j)间的理想间距与实际间距；ΔO_ew(i，j)为新增乙类交叉口I_(i，j)的相位差O_(i，j)与东西方向理想相位差O_ew(i，j)的差值；ΔO_sn(i，j)为新增乙类交叉口I_(i，j)的相位差O_(i，j)与南北方向理想相位差O_sn(i，j)的差值；

为确保相位差O_(i，j)取值尽量不给东西方向与南北方向上下行绿波总带宽带来影响，同时保证东西方向与南北方向绿波带宽均衡，确定相位差O_(i，j)的取值满足式(24)：

若当前协调干道为东西走向协调干道，则计算当前协调干道A_EWi上的任意非节点交叉口I_s(i，k)(2≤i≤n，1≤k≤p_i)的相位差O_s(i，k)，满足式(25)：

式中，T_GWs(i，k)为非节点交叉口I_s(i，k)西进口绿灯时长；v_s(i，k)为节点交叉口I_(i，1)与非节点交叉口I_s(i，k)间的车辆平均行驶速度；a_s(i，k)为节点交叉口I_(i，1)与非节点交叉口I_s(i，k)间的理想间距，满足式(26)：

式中，Δt_{W→Es(i，k)}为非节点交叉口I_s(i，k)西进口超前于东进口直行相位绿灯中心时刻点时长；n_s(i，k)为整数，其取值应使得理想间距a_s(i，k)与实际间距s_s(i，k)最为接近；

若当前协调干道为南北走向协调干道，则计算当前协调干道A_SNj上的任意非节点交叉口I_a(j，l)(2≤j≤m，1≤l≤q_j)的相位差O_a(j，l)，满足式(27)：

式中，T_GWa(j，l)为非节点交叉口I_a(j，l)西进口绿灯时长；v_a(j，l)为节点交叉口I_(1，j)与非节点交叉口I_a(j，l)间的车辆平均行驶速度；Δt_{W→Na(j，l)}为非节点交叉口I_a(j，l)西进口超前于北进口直行相位绿灯中心时刻点时长；a_a(j，l)为节点交叉口I_(1，j)与非节点交叉口I_a(j，l)间的理想间距，满足式(28)：

式中，Δt_{N→Sa(j，l)}为非节点交叉口I_a(j，l)北进口超前于南进口直行相位绿灯中心时刻点时长；n_a(j，l)为整数，其取值应使得理想间距a_a(j，l)与实际间距s_a(j，l)最为接近；

由于在本实施例中所有协调干道上均不存在非节点交叉口，则无需计算当前协调干道上的非节点交叉口相位差；

S304、计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口的最佳相位相序组合方案，其中，所述最佳相位相序组合方案可以使得当前协调干道的干道偏移绿信比与所有已标定协调干道的干道偏移绿信比的绝对值之和取得最小值，将所述可选相位相序组合方案和所述最佳相位相序组合方案计入相应的记录矩阵，此时当前协调干道转为已标定协调干道；

当前协调干道为南北走向协调干道A_SN2，则应计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)，P_(2，2)，P_(3，2)]对应的当前协调干道A_SN2甲类交叉口的最佳相位相序组合方案[P_(1，2)]，其中，所述最佳相位相序组合方案[P_(1，2)]可以使得当前协调干道A_SN2的干道偏移绿信比Δλ_SN2与所有已标定协调干道的干道偏移绿信比的绝对值之和|Δλ_EW1|+|Δλ_SN1|+|Δλ_SN2|取得最小值，将所述可选相位相序组合方案[P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)，P_(2，2)，P_(3，2)]和所述最佳相位相序组合方案[P_(1，2)]计入相应的记录矩阵M_(C，3)，M_(C，3)＝[P_(1，3)，P_(1，4)，P_(2，1)，P_(3，1)，P_(2，2)，P_(3，2)，P_(1，2)]，此时当前协调干道A_SN2转为已标定协调干道；

对于任意东西走向协调干道A_EWi，其干道偏移绿信比Δλ_EWi满足式(29)：

Δλ_EWi＝|Δλ_W→Ei|+|Δλ_E→Wi| (29)

其中，东西走向协调干道A_EWi的协调正向偏移绿信比Δλ_W→Ei与协调反向偏移绿信比Δλ_E→Wi分别满足式(30)和式(31)：

其中，任意节点交叉口I_(i，j)东西向的正向偏移绿信比Δλ_W→E(i，j)和反向偏移绿信比Δλ_E→W(i，j)分别满足式(32)和式(33)：

其中，ΔO_ew(i，j)为节点交叉口I_(i，j)的东西方向相位差偏移量，是节点交叉口I_(i，j)的实际相位差O_(i，j)与东西方向理想相位差O_ew(i，j)的差值，ΔO_ew(i，j)＝O_(i，j)-O_ew(i，j)；

对于任意南北走向协调干道A_SNj，其干道偏移绿信比Δλ_SNj满足式(34)：

Δλ_SNj＝|Δλ_N→Sj|+|Δλ_S→Nj| (34)

其中，南北走向协调干道A_SNj的协调正向偏移绿信比Δλ_N→Sj与协调反向偏移绿信比Δλ_S→Nj分别满足式(35)和式(36)：

其中，任意节点交叉口I_(i，j)南北向的正向偏移绿信比Δλ_N→S(i，j)和反向偏移绿信比Δλ_S→N(i，j)分别满足式(37)和式(38)：

其中，ΔO_sn(i，j)为节点交叉口I_(i，j)的南北方向相位差偏移量，是节点交叉口I_(i，j)的实际相位差O_(i，j)与南北方向理想相位差O_sn(i，j)的差值，ΔO_sn(i，j)＝O_(i，j)-O_sn(i，j)；

S305、重复步骤S301至S304直至所有协调干道均为已标定协调干道；对于最后一条协调干道，即终止协调干道而言，由于在终止协调干道上不再存在乙类交叉口，终止协调干道上的所有交叉口将全部转化为甲类交叉口，如图7所示，因此不同公共信号周期C取值下对应的记录矩阵M_(C，7)仅存储一组使得所有协调干道偏移绿信比绝对值之和最小的终止协调干道上所有交叉口的最佳相位相序组合方案；

南北走向协调干道A_SN2为已标定协调干道，则当前协调干道为东西走向协调干道A_EW2，如图13所示；计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1，3)，P_(1，4)，P_(3，1)，P_(3，2)，P_(2，3)，P_(2，4)]对应的当前协调干道A_EW2甲类交叉口的最佳相位相序组合方案[P_(2，1)，P_(2，2)]，其中，所述最佳相位相序组合方案[P_(2，1)，P_(2，2)]可以使得当前协调干道A_EW2的干道偏移绿信比Δλ_EW2与所有已标定协调干道的干道偏移绿信比的绝对值之和|Δλ_EW1|+|Δλ_SN1|+|Δλ_SN2|+|Δλ_EW2|取得最小值，将所述可选相位相序组合方案[P_(1，3)，P_(1，4)，P_(3，1)，P_(3，2)，P_(2，3)，P_(2，4)]和所述最佳相位相序组合方案[P_(2，1)，P_(2，2)]计入相应的记录矩阵M_(C，4)，M_(C，4)＝[P_(1，3)，P_(1，4)，P_(3，1)，P_(3，2)，P_(2，3)，P_(2，4)，P_(2，1)，P_(2，2)]，此时当前协调干道A_EW2转为已标定协调干道；

东西走向协调干道A_EW2为已标定协调干道，则当前协调干道为南北走向协调干道A_SN3，如图14所示；计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1，4)，P_(3，1)，P_(3，2)，P_(2，4)，P_(3，3)]对应的当前协调干道A_SN3甲类交叉口的最佳相位相序组合方案[P_(1，3)，P_(2，3)]，其中，所述最佳相位相序组合方案[P_(1，3)，P_(2，3)]可以使得当前协调干道A_SN3的干道偏移绿信比Δλ_SN3与所有已标定协调干道的干道偏移绿信比的绝对值之和|Δλ_EW1|+|Δλ_SN1|+|Δλ_SN2|+|Δλ_EW2|+|Δλ_SN3|取得最小值，将所述可选相位相序组合方案[P_(1，4)，P_(3，1)，P_(3，2)，P_(2，4)，P_(3，3)]和所述最佳相位相序组合方案[P_(1，3)，P_(2，3)]计入相应的记录矩阵M_(C，5)，M_(C，5)＝[P_(1，4)，P_(3，1)，P_(3，2)，P_(2，4)，P_(3，3)，P_(1，3)，P_(2，3)]，此时当前协调干道A_SN3转为已标定协调干道；

南北走向协调干道A_SN3为已标定协调干道，则当前协调干道为东西走向协调干道A_EW3，如图15所示；计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1，4)，P_(2，4)，P_(3，4)]对应的当前协调干道A_EW3甲类交叉口的最佳相位相序组合方案[P_(3，1)，P_(3，2)，P_(3，3)]，其中，所述最佳相位相序组合方案[P_(3，1)，P_(3，2)，P_(3，3)]可以使得当前协调干道A_EW3的干道偏移绿信比Δλ_EW3与所有已标定协调干道的干道偏移绿信比的绝对值之和|Δλ_EW1|+|Δλ_SN1|+|Δλ_SN2|+|Δλ_FW2|+|Δλ_SN3|+|Δλ_EW3|取得最小值，将所述可选相位相序组合方案[P_(1，4)，P_(2，4)，P_(3，4)]和所述最佳相位相序组合方案[P_(3，1)，P_(3，2)，P_(3，3)]计入相应的记录矩阵M_(C，6)，M_(C，6)＝[P_(1，4)，P_(2，4)，P_(3，4)，P_(3，1)，P_(3，2)，P_(3，3)]，此时当前协调干道A_EW3转为已标定协调干道；

东西走向协调干道A_EW3为已标定协调干道，则当前协调干道为南北走向协调干道A_SN4，如图16所示；此时协调路网中仅剩下甲类交叉口，则计算并存储不同公共信号周期C取值对应的当前协调干道A_SN4所有交叉口的最佳相位相序组合方案[P_(1，4)，P_(2，4)，P_(3，4)]，其中，所述最佳相位相序组合方案[P_(1，4)，P_(2，4)，P_(3，4)]可以使得所有协调干道偏移绿信比的绝对值之和取得最小值，将所述最佳相位相序组合方案[P_(1，4)，P_(2，4)，P_(3，4)]计入相应的记录矩阵M_(C，7)，M_(C，7)＝[P_(1，4)，P_(2，4)，P_(3，4)]。

步骤S4、利用记录矩阵依次反推获得在上一协调干道组合计算时所确定的甲类交叉口最佳相位相序组合方案，最终获得最优区域绿波协调控制方案；

S401、公共信号周期C的优化步长取为1s，以区域干道总偏移绿信比值最小作为优选原则，从允许最小公共信号周期C_min开始遍历至允许最大公共信号周期C_max，比选出终止协调干道组合计算时可以使得所有协调干道偏移绿信比绝对值之和最小的最佳公共信号周期C_B为104s；

S402、根据最佳公共信号周期C_B，利用记录矩阵M_(104，7)，确定使得所有协调干道偏移绿信比绝对值之和最小的终止协调干道上所有交叉口的最佳相位相序组合方案，再依次利用记录矩阵M_(104，6)，…，M_(104，1)，依次反推获得在上一协调干道组合计算时确定的甲类交叉口最佳相位相序组合方案，直至完成起始协调干道A_EW1的反推，最终获得最优区域绿波协调控制方案。

通过计算246s后求解得到最优区域绿波协调方案，如图17所示，绘制区域绿波空间时距图如图18所示，可见计算所得方案绿波协调优化效果较优。定义绿灯利用率为绿波带宽度与理想绿波带宽度的比值，则65％左右路段的绿灯利用率高于90％，74％以上路段的绿灯利用率大于80％，最大路段绿灯利用率达100％。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (10)

1.一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：协调区域内存在n条东西走向与m条南北走向构成封闭路网的相交协调干道；依次由北往南将第i条东西走向协调干道命名为A_EWi，定义由西往东为协调正向，由东往西为协调反向；由西往东将第j条南北走向的协调干道命名为A_SNj，定义由北往南为协调正向，由南往北为协调反向；若经过某交叉口的干道中仅存在一条协调干道，则定义该交叉口为非节点交叉口，若经过某交叉口的干道存在两条或两条以上协调干道，则定义该交叉口为节点交叉口；将东西走向协调干道A_EWi与南北走向协调干道A_SNj的节点交叉口记为I_(i,j)，相位相序设置方式记为P_(i,j)；将节点交叉口I_(1,1)西进口直行相位的绿灯起始时刻点作为协调区域的相位差零点，所述方法包括以下步骤：

S1、选取第一条东西走向协调干道A_EW1，计算不同公共信号周期C取值下，协调干道上所有乙类交叉口不同可选相位相序组合方案对应的甲类交叉口最佳相位相序组合方案，其中，若经过某交叉口的所有协调干道均为已标定协调干道或当前协调干道，则定义该交叉口为甲类交叉口；若经过某交叉口的所有协调干道中存在但不全为已标定协调干道或当前协调干道，则定义该交叉口为乙类交叉口；若经过某交叉口的所有协调干道中不存在已标定协调干道或当前协调干道，则定义该交叉口为丙类交叉口；

S2、选取第一条南北走向协调干道A_SN1，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案及对应的最佳相位相序组合方案计入相应的记录矩阵；

S3、选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案，并将所述可选相位相序组合方案以及对应的最佳相位相序组合方案计入相应的记录矩阵，重复该步骤直至东西走向协调干道A_EWn与南北走向协调干道A_SNm均为已标定协调干道；

S4、利用记录矩阵依次反推获得在上一协调干道组合计算时所确定的甲类交叉口最佳相位相序组合方案，最终获得最优区域绿波协调控制方案。

2.根据权利要求1所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S1中，所述选取第一条东西走向协调干道A_EW1，计算不同公共信号周期C取值下，协调干道上所有乙类交叉口不同可选相位相序组合方案对应的甲类交叉口最佳相位相序组合方案，具体为：

S101、选取第一条东西走向协调干道A_EW1为当前协调干道，确定当前协调干道A_EW1上甲类交叉口以及乙类交叉口分类情况；

其中，当前协调干道A_EW1上的所有非节点交叉口均为甲类交叉口，当前协调干道A_EW1上的所有节点交叉口均为乙类交叉口；

S102、以节点交叉口I_(1,1)为基准交叉口，计算在不同公共信号周期C的取值下，当前协调干道A_EW1上各交叉口的相位差取值；

S103、计算并存储不同公共信号周期C取值以及所有乙类交叉口不同可选相位相序组合方案[P_(1，1)，P_(1，2)，...，P_(1，m)]对应的当前协调干道A_EW1甲类交叉口最佳相位相序组合方案

其中，所述最佳相位相序组合方案

可以使得当前协调干道A_EW1干道偏移绿信比Δλ_EW1的绝对值

取得最小值，将所述可选相位相序组合方案[P_(1,1),P_(1,2),...,P_(1,m)]和所述最佳相位相序组合方案

计入记录矩阵M_(C,1)，

此时东西走向协调干道A_EW1转为已标定协调干道。

3.根据权利要求2所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S102中，所述不同公共信号周期C的取值下，当前协调干道A_EW1上各交叉口的相位差的确定方式为：

协调区域内各交叉口的公共信号周期C的优化范围记为[C_min,C_max]，其中，C_min为允许最小公共信号周期，C_max为允许最大公共信号周期；根据节点交叉口I_(1,1)，计算当前协调干道A_EW1上的任意节点交叉口I_(1,j)(2≤j≤m)的相位差O_(1,j)，满足式(1)：

式中，T_GW(1,1)和T_GW(1,j)分别为节点交叉口I_(1,1)和节点交叉口I_(1,j)的西进口绿灯时长；v_(1,j)为节点交叉口I_(1,1)和节点交叉口I_(1,j)间的车辆平均行驶速度；a_(1,j)为节点交叉口I_(1,1)和节点交叉口I_(1,j)间的交叉口理想间距，满足式(2)：

式中，Δt_W→E(1,1)和Δt_W→E(1,j)分别为节点交叉口I_(1,1)和节点交叉口I_(1,j)西进口超前于东进口直行相位绿灯中心时刻点的时长；n_(1,j)为整数；

将东西走向协调干道A_EWi的非节点交叉口判定系数定义为0-1变量δ_EWi，若东西走向协调干道A_EWi上存在非节点交叉口，则δ_EWi等于1，东西走向协调干道A_EWi的非节点交叉口的数量记为p_i，由西往东将第k个非节点交叉口记为I_s(i,k)，相位相序设置方式记为P_s(i,k)，若东西走向协调干道A_EWi上不存在非节点交叉口，则δ_EWi等于0；

根据节点交叉口I_(1,1)，计算当前协调干道A_EW1上的任意非节点交叉口I_s(1,k)(1≤k≤p₁)的相位差O_s(1,k)，满足式(3)：

式中，T_GWs(1,k)为非节点交叉口I_s(1,k)西进口绿灯时长；v_s(1,k)为节点交叉口I_(1,1)与非节点交叉口I_s(1,k)间的车辆平均行驶速度；a_s(1,k)为节点交叉口I_(1,1)和非节点交叉口I_s(1,k)间的理想间距，满足式(4)：

式中，Δt_W→Es(1,k)为非节点交叉口I_s(1,k)西进口超前于东进口直行相位绿灯中心时刻点时长；n_s(1,k)为整数。

4.根据权利要求2所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S103中，所述当前协调干道A_EW1干道偏移绿信比Δλ_EW1的确定方式为：

东西走向协调干道A_EW1的干道偏移绿信比Δλ_EW1满足式(5)：

Δλ_EW1＝|Δλ_W→E1|+|Δλ_E→W1| (5)

其中，Δλ_W→E1为东西走向协调干道A_EW1的协调正向偏移绿信比，Δλ_E→W1为协调反向偏移绿信比，分别满足式(6)和式(7)：

其中，任意节点交叉口I_(1,j)东西向的正向偏移绿信比Δλ_W→E(1,j)和反向偏移绿信比Δλ_E→W(1,j)分别满足式(8)与式(9)：

式中，ΔT_(1,j)为节点交叉口I_(1,1)和节点交叉口I_(1,j)理想间距与实际间距的差值D_(1,j)与车辆平均行驶速度v_(1,j)的比值；

任意非节点交叉口I_s(1,k)东西向的正向偏移绿信比Δλ_sW→E(1,k)和反向偏移绿信比Δλ_sE→W(1,k)分别满足式(10)与式(11)：

式中，ΔT_s(1,k)为节点交叉口I_(1,1)与非节点交叉口I_s(1,k)理想间距与实际间距的差值D_s(1,k)与车辆平均行驶速度v_s(1,k)的比值。

5.根据权利要求1所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S2中，所述选取第一条南北走向协调干道A_SN1，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案，具体为：

S201、选取第一条南北走向协调干道A_SN1为当前协调干道，确定当前协调干道A_SN1上甲类交叉口以及乙类交叉口分类情况；

当前协调干道A_SN1与已标定协调干道A_EW1的节点交叉口I_(1,1)以及当前协调干道A_SN1上的所有非节点交叉口均为甲类交叉口，当前协调干道A_SN1上的其余节点交叉口均为乙类交叉口；

S202、以节点交叉口I_(1,1)为基准交叉口，计算在不同公共信号周期C的取值下，当前协调干道A_SN1其余交叉口的相位差取值；

根据节点交叉口I_(1,1)，计算当前协调干道A_SN1上的任意节点交叉口I_(i,1)(2≤i≤n)的相位差O_(i,1)，满足式(12)：

式中，T_GW(i,1)为节点交叉口I_(i,1)西进口绿灯时长；v_(i,1)为节点交叉口I_(1,1)和节点交叉口I_(i,1)间的车辆平均行驶速度；Δt_W→N(1,1)与Δt_W→N(i,1)分别为节点交叉口I_(1,1)和节点交叉口I_(i,1)西进口超前于北进口直行相位绿灯中心时刻点时长；a_(i,1)为节点交叉口I_(1,1)和节点交叉口I_(i,1)间的理想间距，满足式(13)：

式中，Δt_N→S(1,1)和Δt_N→S(i,1)分别为节点交叉口I_(1,1)和节点交叉口I_(i,1)北进口超前于南进口直行相位绿灯中心时刻点时长；n_(i,1)为整数；

将南北走向协调干道A_SNj的非节点交叉口判定系数定义为0-1变量δ_SNj，若南北走向协调干道A_SNj上存在非节点交叉口，则δ_SNj等于1，南北走向协调干道A_SNj的非节点交叉口的数量记为q_j，由北往南将第l个非节点交叉口记为I_a(j,l)，相位相序设置方式记为P_a(j,l)，若南北走向协调干道A_SNj上不存在非节点交叉口，则δ_SNj等于0；

根据节点交叉口I_(1,1)，计算当前协调干道A_SN1上的任意非节点交叉口I_a(1,l)(1≤l≤q₁)的相位差O_a(1,l)，满足式(14)：

式中，T_GWa(1,l)为非节点交叉口I_a(1,l)西进口绿灯时长；v_a(1,l)为节点交叉口I_(1,j)与非节点交叉口I_a(1,l)间的车辆平均行驶速度；Δt_W→Na(1,l)为非节点交叉口I_a(1,l)西进口超前于北进口直行相位绿灯中心时刻点时长；a_a(1,l)为节点交叉口I_(1,1)与非节点交叉口I_a(1,l)间的理想间距，满足式(15)：

式中，Δt_N→Sa(1,l)为非节点交叉口I_a(1,l)北进口超前于南进口直行相位绿灯中心时刻点时长；n_a(1,l)为整数；

S203、计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案[P_(1,2),...,P_(1,m),P_(2,1),...,P_(n,1)]对应的当前协调干道A_SN1甲类交叉口最佳相位相序组合方案

其中，所述最佳相位相序组合方案

可以使得当前协调干道A_SN1干道偏移绿信比Δλ_SN1与已标定协调干道A_EW1干道偏移绿信比Δλ_EW1的绝对值之和|Δλ_EW1|+|Δλ_SN1|取得最小值，将所述可选相位相序组合方案[P_(1,2),...,P_(1,m),P_(2,1),...,P_(n,1)]和所述最佳相位相序组合方案

计入记录矩阵M_(C,2)，

此时南北走向协调干道A_SN1转为已标定协调干道。

6.根据权利要求5所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S203中，所述当前协调干道A_SN1干道偏移绿信比Δλ_SN1的确定方式为：

南北走向协调干道A_SN1的干道偏移绿信比Δλ_SN1满足式(16)：

Δλ_SN1＝|Δλ_N→S1|+|Δλ_S→N1| (16)

其中，南北走向协调干道A_SN1的协调正向偏移绿信比Δλ_N→S1与协调反向偏移绿信比Δλ_S→N1分别满足式(17)和式(18)：

其中，任意节点交叉口I_(i,1)南北向的正向偏移绿信比Δλ_N→S(i,1)和反向偏移绿信比Δλ_S→N(i,1)分别满足式(19)与式(20)：

式中，ΔT_(i,1)为节点交叉口I_(1,1)和节点交叉口I_(i,1)理想间距与实际间距的差值D_(i,1)与车辆平均行驶速度v_(i,1)的比值；

任意非节点交叉口I_a(1,l)南北向的正向偏移绿信比Δλ_aN→S(1,l)和反向偏移绿信比Δλ_aS→N(1,l)分别满足式(21)和式(22)：

式中，ΔT_a(1,l)为节点交叉口I_(1,1)与非节点交叉口I_a(1,l)理想间距与实际间距的差值D_a(j,l)与车辆平均行驶速度v_a(1,l)的比值。

7.根据权利要求1-6任一所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S3中，所述选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道，计算不同公共信号周期C取值下，协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口最佳相位相序组合方案，重复该步骤直至东西走向协调干道A_EWn与南北走向协调干道A_SNm均为已标定协调干道，具体为：

S301、根据已标定协调干道情况，选取使得新增乙类交叉口数量尽可能少的协调干道作为当前协调干道；

其中，确认东西走向协调干道A_EWi和南北走向协调干道A_SNj均为已标定协调干道，东西走向协调干道A_EWi+1和南北走向协调干道A_SNj+1均为未标定协调干道，则选取东西走向协调干道A_EWi+1作为当前协调干道将有j个节点交叉口转化为新增甲类交叉口，新增乙类交叉口m-j个；而选取南北走向协调干道A_SNj+1作为当前协调干道将有i个节点交叉口转化为新增甲类交叉口，新增乙类交叉口n-i个，则若满足式(23)，当前协调干道为东西走向协调干道A_EWi+1，否则为南北走向协调干道A_SNj+1：

n-i≥m-j (23)

S302、确定当前协调干道上甲类交叉口以及乙类交叉口分类情况；

当前协调干道与所有已标定协调干道的节点交叉口以及当前协调干道上的所有非节点交叉口均为甲类交叉口，当前协调干道上的其余节点交叉口均为乙类交叉口；

S303、计算在不同公共信号周期C的取值下，当前协调干道上新增乙类交叉口与非节点交叉口的相位差取值；

S304、计算并存储不同公共信号周期C取值以及协调路网中所有乙类交叉口不同可选相位相序组合方案对应的当前协调干道甲类交叉口的最佳相位相序组合方案，其中，所述最佳相位相序组合方案可以使得当前协调干道的干道偏移绿信比与所有已标定协调干道的干道偏移绿信比的绝对值之和取得最小值，将所述可选相位相序组合方案和所述最佳相位相序组合方案计入相应的记录矩阵，此时当前协调干道转为已标定协调干道；

S305、重复步骤S301至S304直至所有协调干道均为已标定协调干道；对于最后一条协调干道，即终止协调干道而言，由于在终止协调干道上不再存在乙类交叉口，终止协调干道上的所有交叉口将全部转化为甲类交叉口，因此不同公共信号周期C取值下对应的记录矩阵M_(C,n+m)仅存储一组使得所有协调干道偏移绿信比绝对值之和最小的终止协调干道上所有交叉口的最佳相位相序组合方案。

8.根据权利要求7所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S303中，所述不同公共信号周期C的取值下，当前协调干道上新增乙类交叉口与非节点交叉口的相位差的确定方式为：

计算当前协调干道上的任意新增乙类交叉口I_(i,j)(2≤i≤n,2≤j≤m)的相位差O_(i,j)，将节点交叉口I_(i,1)与新增乙类交叉口I_(i,j)间的理想间距与实际间距的差值D_ew(i,j)与车辆平均行驶速度v_ew(i,j)的比值记为ΔT_ew(i,j)；将节点交叉口I_(1,j)与新增乙类交叉口I_(i,j)间的理想间距与实际间距的差值D_sn(i,j)与车辆平均行驶速度v_sn(i,j)的比值记为ΔT_sn(i,j)；为了同时保证东西方向与南北方向的绿波带宽，新增乙类交叉口I_(i,j)的相位差O_(i,j)应该在东西方向理想相位差O_ew(i,j)与南北方向理想相位差O_sn(i,j)之间取值，相位差O_(i,j)的取值范围为[min(O_ew(i,j),O_sn(i,j)),max(O_ew(i,j),O_sn(i,j))]，将东西方向理想相位差O_ew(i,j)与南北方向理想相位差O_sn(i,j)之差的绝对值记为ΔO_(i,j)＝abs(O_ew(i,j)-O_sn(i,j))；

为确保相位差O_(i,j)取值尽量不给东西方向与南北方向上下行绿波总带宽带来影响，同时保证东西方向与南北方向绿波带宽均衡，确定相位差O_(i,j)的取值满足式(24)：

若当前协调干道为东西走向协调干道，则计算当前协调干道A_EWi上的任意非节点交叉口I_s(i,k)(2≤i≤n,1≤k≤p_i)的相位差O_s(i,k)，满足式(25)：

式中，T_GWs(i,k)为非节点交叉口I_s(i,k)西进口绿灯时长；v_s(i,k)为节点交叉口I_(i,1)与非节点交叉口I_s(i,k)间的车辆平均行驶速度；a_s(i,k)为节点交叉口I_(i,1)与非节点交叉口I_s(i,k)间的理想间距，满足式(26)：

式中，Δt_W→Es(i,k)为非节点交叉口I_s(i,k)西进口超前于东进口直行相位绿灯中心时刻点时长；n_s(i,k)为整数；

若当前协调干道为南北走向协调干道，则计算当前协调干道A_SNj上的任意非节点交叉口I_a(j,l)(2≤j≤m,1≤l≤q_j)的相位差O_a(j,l)，满足式(27)：

式中，T_GWa(j,l)为非节点交叉口I_a(j,l)西进口绿灯时长；v_a(j,l)为节点交叉口I_(1,j)与非节点交叉口I_a(j,l)间的车辆平均行驶速度；Δt_W→Na(j,l)为非节点交叉口I_a(j,l)西进口超前于北进口直行相位绿灯中心时刻点时长；a_a(j,l)为节点交叉口I_(1,j)与非节点交叉口I_a(j,l)间的理想间距，满足式(28)：

式中，Δt_N→Sa(j,l)为非节点交叉口I_a(j,l)北进口超前于南进口直行相位绿灯中心时刻点时长；n_a(j,l)为整数。

9.根据权利要求7所述的一种基于迭代优化的区域绿波协调方案求解方法，其特征在于：步骤S304中，所述当前协调干道的干道偏移绿信比的确定方式为：

对于任意东西走向协调干道A_EWi，其干道偏移绿信比Δλ_EWi满足式(29)：

Δλ_EWi＝|Δλ_W→Ei|+|Δλ_E→Wi| (29)

其中，东西走向协调干道A_EWi的协调正向偏移绿信比Δλ_W→Ei与协调反向偏移绿信比Δλ_E→Wi分别满足式(30)和式(31)：

其中，任意节点交叉口I_(i,j)东西向的正向偏移绿信比Δλ_W→E(i,j)和反向偏移绿信比Δλ_E→W(i,j)分别满足式(32)和式(33)：

其中，ΔO_ew(i,j)为节点交叉口I_(i,j)的东西方向相位差偏移量，是节点交叉口I_(i,j)的实际相位差O_(i,j)与东西方向理想相位差O_ew(i,j)的差值，ΔO_ew(i,j)＝O_(i,j)-O_ew(i,j)；

对于任意南北走向协调干道A_SNj，其干道偏移绿信比Δλ_SNj满足式(34)：

Δλ_SNj＝|Δλ_N→Sj|+|Δλ_S→Nj| (34)

其中，南北走向协调干道A_SNj的协调正向偏移绿信比Δλ_N→Sj与协调反向偏移绿信比Δλ_S→Nj分别满足式(35)和式(36)：

其中，任意节点交叉口I_(i,j)南北向的正向偏移绿信比Δλ_N→S(i,j)和反向偏移绿信比Δλ_S→N(i,j)分别满足式(37)和式(38)：

其中，ΔO_sn(i,j)为节点交叉口I_(i,j)的南北方向相位差偏移量，是节点交叉口I_(i,j)的实际相位差O_(i,j)与南北方向理想相位差O_sn(i,j)的差值，ΔO_sn(i,j)＝O_(i,j)-O_sn(i,j)。

10.根据权利要求7所述的一种基于迭代优化的区域绿波协调方案优化方法，其特征在于，步骤S4中，所述利用记录矩阵依次反推获得在上一协调干道组合计算时所确定的甲类交叉口最佳相位相序组合方案，最终获得最优区域绿波协调控制方案，具体为：

S401、设定公共信号周期C的优化步长，以区域干道总偏移绿信比值最小作为优选原则，从允许最小公共信号周期C_min开始遍历至允许最大公共信号周期C_max，比选出终止协调干道组合计算时可以使得所有协调干道偏移绿信比绝对值之和最小的最佳公共信号周期C_B；

S402、根据最佳公共信号周期C_B，利用记录矩阵

确定使得所有协调干道偏移绿信比绝对值之和最小的终止协调干道上所有交叉口的最佳相位相序组合方案，再依次利用记录矩阵

依次反推获得在上一协调干道组合计算时确定的甲类交叉口最佳相位相序组合方案，直至完成起始协调干道A_EW1的反推，最终获得最优区域绿波协调控制方案。

Images