CN113269959A - 一种基于变尺度梯度校正的随机用户均衡交通流分配方法 - Google Patents
一种基于变尺度梯度校正的随机用户均衡交通流分配方法 Download PDFInfo
- Publication number
- CN113269959A CN113269959A CN202110442220.XA CN202110442220A CN113269959A CN 113269959 A CN113269959 A CN 113269959A CN 202110442220 A CN202110442220 A CN 202110442220A CN 113269959 A CN113269959 A CN 113269959A
- Authority
- CN
- China
- Prior art keywords
- path
- flow
- iteration
- vector
- gradient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000012937 correction Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 239000013598 vector Substances 0.000 claims description 52
- 230000006870 function Effects 0.000 claims description 28
- 230000019771 cognition Effects 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Analytical Chemistry (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于变尺度梯度校正的随机用户均衡交通流分配方法,本方法在既约变量空间中进行迭代,每步迭代通过生成一个校正矩阵来动态调整目标函数的梯度信息,该方法具有超线性的收敛速度,且无需计算任何逆矩阵,能够大大减少运算的复杂度;且相比梯度投影法,该方法能够提高迭代后期的收敛速度,并大大节省迭代后期的CPU时间,非常适合求解大规模无约束优化问题。通过运用基于变尺度梯度校正的随机用户均衡交通流分配方法,能够提高求解效率,节省运算时间。
Description
技术领域
本发明涉及交通分配技术领域,特别是一种基于变尺度梯度校正的随机用户均衡交通流分配方法。
背景技术
交通分配模型是城市交通系统分析和设计中最基础的工具之一,在城市交通规划中占有核心地位。交通分配模型用于预测在均衡状态下,交通网络中的路段流量或者路径流量。交通分配模型可分为确定性用户均衡分配模型(User Equilibrium简称UE)及随机用户均衡分配模型(Stochastic User Equilibrium简称SUE),其中,随机用户均衡分配模型考虑了现实生活中用户在路径选择中具有一定的随机性的要求,更具有现实性指导意义。随机用户均衡分配模型假设,人们对于路径出行时间的感知是有误差的,在SUE解点上,没有一个出行者可以通过单方面改变路径来降低自己的理解出行时间。随机用户均衡模型考虑了出行者感知的出行费用,能更好地描述出行者在网络中的出行选择行为,基于Logit的SUE模型由于具有结构简单,可解释性强等特点,在交通分配中的应用最为广泛。
目前,Logit型随机用户均衡路径交通流分配常用的分配方法是梯度投影法,该方法的核心思想是将当前迭代点沿目标函数的负梯度方向移动一个步长,并将结果投影到问题的可行域上,从而得到新的迭代点,但该方法只在原始变量空间中进行迭代,只具有线性的收敛速度,当迭代点接近最优解时收敛速度相对较慢。
发明内容
本发明所要解决的技术问题是克服现有技术的不足而提供一种基于变尺度梯度校正的随机用户均衡交通流分配方法,本方法在既约变量空间中进行迭代,每步迭代通过生成一个校正矩阵来动态调整目标函数的梯度信息,该方法具有超线性的收敛速度,且无需计算任何逆矩阵,能够大大减少运算的复杂度;且相比梯度投影法,该方法能够提高迭代后期的收敛速度,并大大节省迭代后期的CPU时间,非常适合求解大规模无约束优化问题。
本发明为解决上述技术问题采用以下技术方案:
根据本发明提出的一种基于变尺度梯度校正的随机用户均衡交通流分配方法,
考虑一个交通网络G(N,A),其中,N是节点集合、A是路段集合,令W表示路网中所有OD对的集合,Rw表示OD对w之间路径集合,w∈W,bw表示OD对w之间的交通需求量;设分别为OD对w之间路径r上的流量与阻抗,r∈Rw,xa、ta(xa)分别为路段a上的流量与阻抗,a∈A,ta(xa)为xa的可微单调递增函数,且其中,Ca为路段a的通行能力、为自由流行驶时间,γ和m为确定性参数;
假设为出行者的路径流量向量,为出行者的路径行驶时间向量,为出行者的路段流量向量,为出行者的路段行驶时间向量,ta为出行者的路段行驶时间;定义路段-路径关联矩阵其中为指标变量,如果路段a在连接OD对w的路径r上,否则
根据上述符号定义,随机用户均衡模型表示为如下最优化问题P1,P1包括公式(1)-公式(4):
其中,Z(f)为目标函数,x为路段流量,ta(x)为路段阻抗,θ为参数用于反映人们对网络阻抗的认知程度;式(2)表示OD需求量与路径流量之间的守恒关系,式(3)是路径流量非负约束,式(4)表示路段流量与路径流量的关系;
求解最优化问题P1的方法具体如下:
步骤1、组织交通调查,确定每个OD对之间的交通需求量{bw,w∈W}及每个OD对之间的路径集合{Rw,w∈W};
则OD对w∈W所对应的所有既约路径的下标集合Qw表示为:
其中rN为既约路径下标,则对应的既约路径流量向量fN为
对(8)式求一阶偏导,得到第k次迭代时的既约梯度向量gk;
步骤7、令第k次迭代时的搜索方向向量dk=-Hkgk,Hk是第k次迭代时的矩阵;
其中,上标T代表矩阵转置。
作为本发明所述的一种基于变尺度梯度校正的随机用户均衡交通流分配方法进一步优化方案,所述步骤2的初始路径流量向量f0由以下3个子步骤求得:
步骤2.1、令xa=0,a∈A,由ta(xa)得出各路段的自由流行驶时间;
步骤2.3、按照下式,将出行者OD需求量加载到出行路径上:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明首先运用变量消去法,将Logit型SUE模型转化为一个无约束的最优化问题;然后,运用基于变尺度梯度校正的随机用户均衡交通流分配方法求解该无约束最优化问题得到SUE最优解;最后,在SiouxFalls网络上,分别对变尺度梯度校正法与梯度投影法进行对比,以证明本发明所提方法的有效性;
(2)本发明所提出的方法只需要利用目标函数值和一阶导数构造出目标函数的曲率近似,通过前后两次迭代的梯度信息与下降方向信息进行向量点乘运算得到校正矩阵,而不需要计算任何逆矩阵,从而简化了运算的复杂度,大大降低了计算工作量;
(3)本发明提出的方法能够在既约变量空间中进行迭代,每步迭代通过生成一个校正矩阵,动态调整目标函数的梯度信息,该方法具有超线性收敛性,能够大大提高算法的收敛速度;
(4)因此,通过运用基于变尺度梯度校正的随机用户均衡交通流分配方法,能够提高求解效率,节省运算时间,是求解大规模非线性问题的一种切实可行的有效方法;本发明是对交通分配理论和实践的探索和创新,具有较强的理论研究意义和实践指导价值。
附图说明
图1是本发明流程图。
图2是本发明方法在Sioux Falls网络上迭代次数的收敛性能。
图3是本发明方法在Sioux Falls网络上CPU时间的收敛性能。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明进行详细描述。
首先,运用变量消去法,将Logit型SUE模型转化为一个无约束的最优化问题;然后,运用基于变尺度梯度校正的随机用户均衡交通流分配方法求解该无约束最优化问题得到SUE最优解;最后,在SiouxFalls网络上,分别对变尺度梯度校正法与梯度投影法进行对比,以证明本发明所提方法的有效性。
(1)模型描述
考虑一个交通网络G(N,A),其中,N是节点集合、A是路段集合,令W表示路网中所有OD对的集合,Rw表示OD对w之间路径集合,w∈W,bw表示OD对w之间的交通需求量;设分别为OD对w之间路径r上的流量与阻抗,r∈Rw,xa、ta(xa)分别为路段a上的流量与阻抗,a∈A,ta(xa)为xa的可微单调递增函数,且其中,Ca为路段a的通行能力、为自由流行驶时间,γ和m为确定性参数;
假设为出行者的路径流量向量,为出行者的路径行驶时间向量,为出行者的路段流量向量,为出行者的路段行驶时间向量,ta为出行者的路段行驶时间;定义路段-路径关联矩阵其中为指标变量,如果路段a在连接OD对w的路径r上,否则
根据上述符号定义,随机用户均衡模型表示为如下最优化问题P1,P1包括公式(1)-公式(4):
其中,Z(f)为目标函数,x为路段流量,ta(x)为路段阻抗,θ为参数用于反映人们对网络阻抗的认知程度;式(2)表示OD需求量与路径流量之间的守恒关系,式(3)是路径流量非负约束,式(4)表示路段流量与路径流量的关系。
(2)算法设计
求解出最优化问题P1,即能实现随机用户均衡交通流分配;求解最优化问题P1的方法具体如下:
步骤1、组织交通调查,确定每个OD对之间的交通需求量{bw,w∈W}及每个OD对之间的路径集合{Rw,w∈W};
则OD对w∈W所对应的所有既约路径的下标集合Qw表示为:
其中rN为既约路径下标,则对应的既约路径流量向量fN为
对(8)式求一阶偏导,得到第k次迭代时的既约梯度向量gk;
步骤7、令第k次迭代时的搜索方向向量dk=-Hkgk,Hk是第k次迭代时的矩阵;
其中,上标T代表矩阵转置。
步骤2的初始路径流量向量f0由以下3个子步骤求得:
步骤2.1、令xa=0,a∈A,由ta(xa)得出各路段的自由流行驶时间;
步骤2.3、按照下式,将出行者OD需求量加载到出行路径上:
运行上述如图1的算法,最终可以找到一个SUE最优解,值得指出的是,上述算法给出的是求解过程的总体框架,在算法实施过程中,会涉及到更多的细节问题,需要结合具体实施方式作进一步说明。
(3)算法分析
下面通过一个测试用例来对比基于变尺度梯度校正的随机用户均衡交通流分配方法与梯度投影法的性能,进而验证本发明的有效性。测试网络采用Sioux Falls网络,Sioux Falls网络是一个中型网络,它由76个路段,24个节点和528个OD对组成,该网络取自Bar-Gera。
通过算法性能测试比较变尺度梯度校正法与梯度投影法在不同迭代阶段的收敛速度。本测试用例中,使用Leurent提出的收敛性依据。
其中,Zk为第k次迭代时,最优化问题P1目标函数的值,Z*为目标函数的最优值。
本发明对比了变尺度梯度校正法与梯度投影法在不同迭代次数和CPU时间下的收敛性能。图2显示了上述方法在Sioux Falls网络上迭代次数的收敛性能。从图2中可以看出,在迭代早期,变尺度梯度校正法的斜率与梯度投影法类似,这意味着两种算法的收敛速度也类似;在迭代后期,变尺度梯度校正法的斜率远大于梯度投影法,这意味着变尺度梯度校正法的收敛速度远快于梯度投影法。这种现象验证了变尺度梯度校正法具有局部超线性的收敛速度,而梯度投影法仅具有线性的收敛速度。
图3显示了上述方法在Sioux Falls网络上CPU时间的收敛性能。从图3可以看出,在迭代早期,梯度投影法的计算效率与变尺度梯度校正法类似;在迭代后期,梯度投影法的计算效率明显低于变尺度梯度校正法。
综上所述,本发明所提出的基于变尺度梯度校正的随机用户均衡交通流分配方法比梯度投影法有更好的性能。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围内。
Claims (2)
1.一种基于变尺度梯度校正的随机用户均衡交通流分配方法,其特征在于,
考虑一个交通网络G(N,A),其中,N是节点集合、A是路段集合,令W表示路网中所有OD对的集合,Rw表示OD对w之间路径集合,w∈W,bw表示OD对w之间的交通需求量;设fr w、分别为OD对w之间路径r上的流量与阻抗,r∈Rw,xa、ta(xa)分别为路段a上的流量与阻抗,a∈A,ta(xa)为xa的可微单调递增函数,且其中,Ca为路段a的通行能力、为自由流行驶时间,γ和m为确定性参数;
假设为出行者的路径流量向量,为出行者的路径行驶时间向量,为出行者的路段流量向量,为出行者的路段行驶时间向量,ta为出行者的路段行驶时间;定义路段-路径关联矩阵其中为指标变量,如果路段a在连接OD对w的路径r上,否则
根据上述符号定义,随机用户均衡模型表示为如下最优化问题P1,P1包括公式(1)-公式(4):
其中,Z(f)为目标函数,x为路段流量,ta(x)为路段阻抗,θ为参数用于反映人们对网络阻抗的认知程度;式(2)表示OD需求量与路径流量之间的守恒关系,式(3)是路径流量非负约束,式(4)表示路段流量与路径流量的关系;
求解最优化问题P1的方法具体如下:
步骤1、组织交通调查,确定每个OD对之间的交通需求量{bw,w∈W}及每个OD对之间的路径集合{Rw,w∈W};
则OD对w∈W所对应的所有既约路径的下标集合Qw表示为:
其中rN为既约路径下标,则对应的既约路径流量向量fN为
对(8)式求一阶偏导,得到第k次迭代时的既约梯度向量gk;
步骤7、令第k次迭代时的搜索方向向量dk=-Hkgk,Hk是第k次迭代时的矩阵;
其中,上标T代表矩阵转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110442220.XA CN113269959B (zh) | 2021-04-23 | 2021-04-23 | 一种基于变尺度梯度校正的随机用户均衡交通流分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110442220.XA CN113269959B (zh) | 2021-04-23 | 2021-04-23 | 一种基于变尺度梯度校正的随机用户均衡交通流分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113269959A true CN113269959A (zh) | 2021-08-17 |
CN113269959B CN113269959B (zh) | 2022-04-15 |
Family
ID=77229248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110442220.XA Active CN113269959B (zh) | 2021-04-23 | 2021-04-23 | 一种基于变尺度梯度校正的随机用户均衡交通流分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113269959B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116311932A (zh) * | 2023-03-16 | 2023-06-23 | 东南大学 | 一种MaaS背景下考虑混合均衡的动态交通分配方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694278A (zh) * | 2018-04-27 | 2018-10-23 | 东南大学 | 一种基于道路负载均衡的城市离散交通网络设计方法 |
CN108765944A (zh) * | 2018-05-30 | 2018-11-06 | 东南大学 | 基于多重路径集的最优交通流预测方法与拥挤收费方法 |
CN108776721A (zh) * | 2018-04-27 | 2018-11-09 | 东南大学 | 一种基于目标配流的城市离散交通网络设计方法 |
CN108776722A (zh) * | 2018-04-27 | 2018-11-09 | 东南大学 | 一种基于设计速度的城市离散交通网络设计方法 |
CN108804801A (zh) * | 2018-05-25 | 2018-11-13 | 东南大学 | 基于目标配流的城市离散交通网络设计r语言实现方法 |
-
2021
- 2021-04-23 CN CN202110442220.XA patent/CN113269959B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694278A (zh) * | 2018-04-27 | 2018-10-23 | 东南大学 | 一种基于道路负载均衡的城市离散交通网络设计方法 |
CN108776721A (zh) * | 2018-04-27 | 2018-11-09 | 东南大学 | 一种基于目标配流的城市离散交通网络设计方法 |
CN108776722A (zh) * | 2018-04-27 | 2018-11-09 | 东南大学 | 一种基于设计速度的城市离散交通网络设计方法 |
CN108804801A (zh) * | 2018-05-25 | 2018-11-13 | 东南大学 | 基于目标配流的城市离散交通网络设计r语言实现方法 |
CN108765944A (zh) * | 2018-05-30 | 2018-11-06 | 东南大学 | 基于多重路径集的最优交通流预测方法与拥挤收费方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116311932A (zh) * | 2023-03-16 | 2023-06-23 | 东南大学 | 一种MaaS背景下考虑混合均衡的动态交通分配方法 |
CN116311932B (zh) * | 2023-03-16 | 2024-03-01 | 东南大学 | 一种MaaS背景下考虑混合均衡的动态交通分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113269959B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Patriksson et al. | A mathematical model and descent algorithm for bilevel traffic management | |
Chiou | Bilevel programming for the continuous transport network design problem | |
CN105678054B (zh) | 一种基于a星策略的最优多会合点路径搜索方法及装置 | |
CN109951392B (zh) | 一种基于深度学习的中大型网络智能路由选择方法 | |
Cantarella et al. | Stochastic equilibrium assignment with variable demand: theoretical and implementation issues | |
CN108647837B (zh) | 考虑出行者路径选择惰性的网络交通流预测方法 | |
CN108765944B (zh) | 基于多重路径集的最优交通流预测方法与拥挤收费方法 | |
CN113269959B (zh) | 一种基于变尺度梯度校正的随机用户均衡交通流分配方法 | |
Wu et al. | Learning to effectively estimate the travel time for fastest route recommendation | |
CN114518755A (zh) | 全局路径的规划方法、装置、设备及存储介质 | |
Zhou et al. | Research on optimal path based on Dijkstra algorithms | |
CN107300388B (zh) | 基于Q-learning算法和回声状态网络的骑行旅游路线规划方法 | |
Larsson et al. | A column generation procedure for the side constrained traffic equilibrium problem | |
CN104821086A (zh) | 一种定位大规模交通网络中低效路段组合的方法 | |
CN114969234A (zh) | 一种面向全量城市路网的向量表征方法 | |
CN116702389B (zh) | 一种针对混合交通流的嵌套流量计算方法 | |
Sonar | Classical finite volume methods | |
CN113326468B (zh) | 一种基于路径尺度修正及两阶段优化的交通流分配方法 | |
CN116341634A (zh) | 神经结构搜索模型的训练方法、装置及电子设备 | |
CN111833596B (zh) | 考虑出行者决策惯性的逐日路段流量预测方法 | |
CN105427581A (zh) | 基于浮动车数据的交通仿真方法及系统 | |
Ran et al. | Virtual SDN network embedding algorithm based on load balance | |
JP2020119108A (ja) | データ処理装置、データ処理方法、データ処理プログラム | |
Guo et al. | Safe zeroth-order optimization using linear programs | |
Aguirre et al. | A hybrid selection strategy using scalarization and adaptive ε-ranking for many-objective optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |