CN115482655A - 一种基于偏最小二乘Kriging模型的路径诱导方法 - Google Patents
一种基于偏最小二乘Kriging模型的路径诱导方法 Download PDFInfo
- Publication number
- CN115482655A CN115482655A CN202210405251.2A CN202210405251A CN115482655A CN 115482655 A CN115482655 A CN 115482655A CN 202210405251 A CN202210405251 A CN 202210405251A CN 115482655 A CN115482655 A CN 115482655A
- Authority
- CN
- China
- Prior art keywords
- road section
- road
- model
- travel time
- upstream
- 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
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
-
- 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
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/09—Arrangements for giving variable traffic instructions
- G08G1/0962—Arrangements for giving variable traffic instructions having an indicator mounted inside the vehicle, e.g. giving voice messages
- G08G1/0968—Systems involving transmission of navigation instructions to the vehicle
- G08G1/096805—Systems involving transmission of navigation instructions to the vehicle where the transmitted instructions are used to compute a route
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/09—Arrangements for giving variable traffic instructions
- G08G1/0962—Arrangements for giving variable traffic instructions having an indicator mounted inside the vehicle, e.g. giving voice messages
- G08G1/0968—Systems involving transmission of navigation instructions to the vehicle
- G08G1/096833—Systems involving transmission of navigation instructions to the vehicle where different aspects are considered when computing the route
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于偏最小二乘Kriging模型的路径诱导方法,包括:S1、确定拥堵上游路段;S2、确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的优化问题,并构建中观交通仿真模型;S3、基于拉丁超立方抽样在解空间中均匀抽取样本,并计算样本的目标函数值,建立仿真样本集合;S4、建立基于偏最小二乘Kriging模型的代理模型;S5、利用群遗传算法和改善期望函数在解空间内进行二次采样,获得最大价值采样点;S6、计算最大价值采样点的目标函数值,并更新仿真样本集合;S7、重复迭代直至达到最大迭代次数,输出当前得到的最优解。与现有技术相比,本发明具有计算复杂度低、对大规模复杂拥堵路段的改善效果显著的优点。
Description
技术领域
本发明涉及交通管控与智能交通领域,尤其是涉及一种基于偏最小二乘 Kriging模型的路径诱导方法。
背景技术
随着机动车保有量的持续增长,多数城市的交通供给能力已无法满***通需求,交通拥堵现象在城市路网中愈发常见。各类偶发性事件,如交通事故、极端天气等同样严重了影响交通***的通行效率。拥堵不仅造成了经济损失和资源浪费,还会带来严重的环境污染问题。为了缓解交通拥堵现象,先进的出行者信息***(Advanced TravelerInformation Systems,ATIS)成为提升城市路网通行效率的重要手段。
ATIS可将路网的实时交通状况和突发事件信息发送给用户(出行者),触发用户的在途路径选择行为(En-route Diversion),帮助用户及时调整出行路线,从而避开发生拥堵或交通事故的路段。ATIS的实现形式包括各类可变信息标志(Variable Message Sign,VMS):拥堵警告、速度建议、路径提示,导航***(Navigation System),以及智能网联环境下基于车载无线通信技术的路径引导(Connected Vehicle Guidance)等。大量研究表明,ATIS信息触发的在途路径选择行为不仅可以缩短个体用户的出行时间,还将显著提升整个***的运行效率。
早期关于ATIS设计的研究大都关注VMS的布设位置优化问题。后来的研究多将VMS的布设位置与VMS的种类、显示信息等变量综合考虑,并协同优化。 ATIS优化设计中需要考虑一个重要的因素:用户对实时信息的反应和遵从度 (Compliance Rate)。用户的遵从度与在途换路的用户比例(Diversion Rate)直接相关,多数VMS布设优化的研究都分析了该参数的灵敏度,结果显示在途换路的用户比例对路网通行效率有显著影响。因此若要通过ATIS发布实时的路径诱导信息以改善拥堵网络的通行效率,首先需要确定拥堵瓶颈路段上游需要重新规划路线的在途用户比例,然后发布适当的ATIS信息引导这部分用户重新选择路线。目前已有大量研究讨论了如何合理设计ATIS信息,从而影响在途用户对该信息的遵从度并达到期望的换路比例:通过动态调整ATIS信息,或引入其他交通控制手段调整受ATIS信息影响在途用户的比例。此外,确定瓶颈上游所需的最佳在途换路用户比例也是一个重要的课题,但目前对该问题的研究较少。大规模城市路网在高峰期容易多点同时发生常发性拥堵,且可能伴随着由于突发事件引起的偶发性拥堵,此时网络的通行效率将显著降低。然而该类场景的优化决策空间较大,优化难度高,目前关于该场景下ATIS优化设计的研究较少。
在方法层面,ATIS优化设计的研究中通常有两种建模方式:解析模型及仿真模型。解析模型通过排队论模型、确定性用户均衡或随机用户均衡等建立起决策变量和目标函数之间的映射关系,并基于一定的优化算法进行求解。然而解析模型很难同时准确刻画大规模城市路网中的交通动态演化过程、由于拥堵引发的流量溢出以及用户对ATIS信息的反应行为等复杂交通现象。仿真模型即利用成熟的交通仿真工具模拟实际的交通流演变过程,基于仿真结果计算目标函数值。与解析模型相比,仿真模型可以准确刻画上述的复杂交通现象,但也具有不可解析、计算成本高、仿真噪声大的难题。早期基于仿真模型的ATIS优化设计大都使用启发式算法,该方法在优化过程中需要频繁调用仿真模型计算目标函数值,在有限的计算资源内很难获得理想的优化结果。最新的研究利用基于贝叶斯优化的仿真优化方法设计 ATIS,建立高斯过程代理模型,通过代理模型评估目标函数值可以极大减少仿真模型的调用次数。然而该研究使用的代理模型计算复杂度高,在高维仿真优化问题中的计算效率较低。因此受限于仿真模型和代理模型较高的计算成本,目前仿真优化方法只能处理较低维度的ATIS设计问题。
基于以上分析,现有技术目前还存在以下问题:
1)没有考虑在大规模偶发拥堵和常发拥堵共存场景中的在途换路用户比例确定;
2)高维度ATIS仿真优化的效率低。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供了一种高效的基于偏最小二乘Kriging模型的路径诱导方法,该方法考虑了大规模偶发拥堵和常发拥堵共存场景中的在途换路用户比例确定问题,且提高了高维度ATIS仿真优化的效率。
本发明的目的可以通过以下技术方案来实现:
本发明提供了一种基于偏最小二乘Kriging模型的路径诱导方法,该方法包括以下步骤:
步骤S1、基于实际的交通供给、需求及突发交通状况信息,确定需要发送“前方拥堵警告”信息的拥堵上游路段;
步骤S2、确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的连续优化问题,并构建用于计算目标函数值的中观交通仿真模型;
步骤S3、基于拉丁超立方抽样在解空间中均匀抽取样本,并带入所述中观交通仿真模型计算样本的目标函数值,建立仿真样本集合;
步骤S4、利用仿真样本集合,建立基于偏最小二乘Kriging模型的代理模型;
步骤S5、基于代理模型,利用群遗传算法和改善期望函数在解空间内进行二次采样,获得最大价值采样点;
步骤S6、将最大价值采样点带入所述中观交通仿真模型中计算其目标函数值,并更新仿真样本集合;
步骤S7、判断是否达到最大迭代次数,若否,则返回步骤S4,否则结束输出当前得到的最优解。
优选地,所述步骤S1中的拥堵上游路段包括常发性拥堵上游路段和偶发性拥堵上游路段;
所述偶发性拥堵上游路段采用直接确定法进行确定;所述直接确定法为将偶发性拥堵上游的上下匝道设定范围内路段作为偶发性拥堵上游路段。
优选地,所述常发性拥堵上游路段采用主成分分析法进行确定;所述主成分分析法包括:路段选取、路段降维以及路段筛选。
优选地,所述路段选取具体为:
通过用于判断路段车流拥堵程度的路段速度损失率函数选择拥堵路段上游的所有路段中速度损失大于预设值的路段,得到路段选取结果Lb,1;
所述路段速度损失率函数表达式为:
优选地,所述路段降维为使用主成分分析法对路段选取结果Lb,1进行降维,具体包括以下子步骤:
1)将样本数据标准化,得到标准化阵Z,表达式为:
式中,zij为标准化阵Z的元素,sj分别为变量Yj对应样本的均值和标准差;变量Yj包括:各在途换路用户比例下***总出行时间减少量平均值Y1、***总出行时间减少量的最小值Y2、***总出行时间减少量的最大值Y3、***总出行时间减少量的上四分位数Y4、***总出行时间减少量的中位数Y5、***总出行时间减少量的下四分位数Y6、各在途换路用户比例下***总出行时间减少量的极差Y7、以及各在途换路用户比例下***总出行时间减少量总和Y8;
2)计算标准化阵Z的协方差矩阵R,表达式为
式中,rij为协方差矩阵R的元素,n为样本容量:
3)计算协方差矩阵R的特征值λ1≥λ2≥…≥λn≥0以及对应特征向量 u1,u2,…,un;
4)计算各主成分yi的贡献率bj以及y1,y2,…,yp累计贡献率αp:
选取特征值大于1且累计贡献率约达到90%的成分作为主成分;
5)获得成分得分系数矩阵,计算各成分得分,之后通过旋转后的方差百分比与各成分得分求出最后的综合评价得分,得到降维后的上游路段Lb,2。
优选地,所述路段筛选具体为:引入候选路段的最短间距限制,从降维后的上游路段Lb,2中进行筛选,将筛选后的结果作为常发性拥堵上游路段Lb。
优选地,所述步骤S2中确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的连续优化问题,并构建用于计算目标函数值的中观交通仿真模型,具体为:
1)定义决策变量为n个上游路段在途换路的用户比例,目标函数为***总出行时间的期望f(x)最小;其中,所述决策变量表达式为:
x=[x1,x2,...,xn]T
式中,1,2,...,n表示拥堵上游发布ATIS信息的路段,即步骤S1中确定的拥堵上游路段L,xi表示第i个路段上在途换路用户的比例;
2)将p次重复仿真所得***总出行时间的平均值y(x)作为***总出行时间期望f(x)的近似估计值,将原目标函数转化为近似目标函数,优化数学模型具体如下:
所述原目标函数表达式为:
式中,f(x)为***总出行时间的期望;为***总出行时间,包含一个随机变量;r,R分别为用户出行的起点节点及所有起点节点的集合;s,S分别为用户出行的终点节点及所有终点节点的集合;t,T分别为用户的出发时刻及所有出发时刻的集合;为高峰时段t时刻从起点r到终点s的总出行量;为高峰时段t 时刻从起点r到终点s的平均出行时间,与x有关;为非高峰时段t时刻从起点r 到终点s的总出行量;为非高峰时段t时刻从起点r到终点s的平均出行时间,与x无关;
所述近似目标函数表达式为:
约束条件为:
xmin≤x≤xmax
式中,xmin=[0,0,...,0]T为在途换路用户比例的下限;xmax=[1,1,...,1]T为在途换路用户比例的上限。
优选地,所述步骤S4具体为:
步骤S4.1、构建基于偏最小二乘Kriging模型的代理模型,表达式为:
fm(x)=μ+ε(x)
式中,μ是目标函数的平均值;ε(x)为估计误差,服从均值为0、方差为σ2的高斯分布,满足以下协方差关系:
Cov[ε(x(i)),ε(x(j))]=σ2ψ(x(i),x(j))
式中,x(i)为第i个路段上在途换路用户的比例;ψ为偏最小二乘Kriging基函数,表达式为:
式中,R为基于仿真样本集合D中的观测数据得到的相关矩阵,表达式为:
Ri,j=ψ(x(i),x(j));
r为协方差向量,表达式为:
r={Cov[ε(x*),ε(x(1))],...,Cov[ε(x*),ε(x(d))]}T
步骤S4.2、采用网格搜索法和交叉验证法确定代理模型的初始参数θ0的最佳取值。
优选地,所述步骤S4.2具体为:
1)将仿真样本集合D随机分为K份,得到子样本集合:D1,D2,...,DK;
2)将θ0的取值网格化,Θ={0.001,0.01,0.1,1,10};
3)将Θ中的值逐一赋给θ0,使用K折交叉验证计算代理模型的平均预测误差;
4)选择Θ中平均预测误差最小的参数作为θ0的初始参数值。
优选地,所述步骤S5包括以下子步骤:
步骤S5.1、利用代理模型辅助的多种群遗传算法进行同时寻优;
步骤S5.2、选择样本点改善函数I(x)和样本点改善期望函数E[I(x)]作为采集函数,表达式为:
I(x)=ymin-y(x)
步骤S5.3、在解空间内进行二次采样,获得最大价值采样点。
与现有技术相比,本发明具有以下优点:
1)本发明提出的基于偏最小二乘Kriging代理模型的仿真优化方法可以确定在大规模偶发拥堵和常发拥堵共存场景中最优的在途换路用户比例,跟真实场景下的交通拥堵路况,提高了模型预测的精确度;
2)本发明使用偏最小二乘Kriging模型作为中观仿真模型的代理模型,与常用的Kriging模型相比,该模型在高维空间中的计算复杂度更低,且模型的预测精度更高;
3)本发明提出的基于偏最小二乘Kriging代理模型的仿真优化算法可以高效解决较高维度的ATIS的设计问题;
4)本发明使用仿真模型作为优化方案的评价工具,能够模拟交通***的动态演变过程,并且优化方法不依赖于某种具体的仿真模型,具有较高的可移植性。
附图说明
图1为本发明的流程示意图;
图2为本发明实施例中A市徐汇区路网的示意图;
图3为本发明实施例中常发性拥堵节点或路段的上游路段示意图;
图4为本发明实施例中交通需求标定结果的示意图;
图5为本发明实施例中地面主干道交通流模型参数标定结果的示意图;
图6为本发明实施例中快速交通流模型参数标定结果的示意图;
图7为本发明实施例中优化方案与基准方案在路段平均速度上的差异图;
图8为本发明实施例中优化方案与基准方案在路段平均行程时间上的差异图;
图9为本发明实施例中优化方案与基准方案在网络排队车辆总数上的差异图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
实施例
如图1所示,本实施例提供了一种基于偏最小二乘Kriging模型的路径诱导方法,该方法包括以下步骤:
步骤S1、基于实际的交通供给、需求及突发交通状况信息,确定需要发送“前方拥堵警告”信息的拥堵上游路段,具体为:
步骤S1.1:偶发性拥堵节点或路段的拥堵上游路段直接确定,即选择拥堵节点或路段上游上下匝道附近的路段作为偶发拥堵的上游路段La。
步骤S1.2:常发性拥堵节点或路段的拥堵上游路段通过主成分分析法确定,包括三个步骤:路段选取、路段降维以及路段筛选。
路段选取具体为通过路段速度损失率函数选择拥堵节点或路段上游的所有路段中速度损失较大的路段。路段速度损失率可以判断该路段的车流拥堵程度,其形式如下,其中Li为第i条路段的速度损失率,T为总的时间片段数,为第i条路段在t时间片段的速度,vi,lim为第i条路段的限速。
选取整个网络中速度损失率较大的路段及其上游路段作为初始选取的上游路段Lb,1。使用中观仿真技术分别计算初始选取路段Lb,1上的在途换路用户比例在 10%、20%、30%、40%、50%、60%、70%、80%、90%、100%时的***总出行时间。
路段降维具体为使用主成分分析法对初始选取的路段结果Lb,1进行降维。定义 8个变量,Y1-Y8。其中,Y1为各在途换路用户比例下***总出行时间减少量平均值,Y2为***总出行时间减少量的最小值,Y3为***总出行时间减少量的最大值, Y4为***总出行时间减少量的上四分位数,Y5为***总出行时间减少量的中位数, Y6为***总出行时间减少量的下四分位数,Y7为各在途换路用户比例下***总出行时间减少量的极差,Y8为各在途换路用户比例下***总出行时间减少量总和。对该 8个变量进行主成分分析,可以将多指标转换为少数几个特征性指标,并借助综合评价得分对不同路径诱导方案的效果进行评价,步骤如下:
2)计算标准化阵Z的协方差矩阵R,其中为rij矩阵R的元素,n为样本容量:
3)计算矩阵R的特征值λ1≥λ2≥…≥λn≥0及对应特征向量u1,u2,…,un;
4)计算各主成分yi的贡献率bj以及y1,y2,…,yp累计贡献率αp:
选取特征值大于1,累计贡献率约达到90%的成分作为主成分。
5)获得成分得分系数矩阵,计算各成分得分;之后通过旋转后的方差百分比与各成分得分求出最后的综合评价得分。最终得到降维后的上游路段Lb,2。
路段筛选具体为,引入候选路段的最短间距限制,保证得到的拥堵上游路段具有一定的位置间隔。从降维后的上游路段Lb,2中进行筛选,将筛选后的结果作为常发性拥堵的上游路段Lb。
步骤S1.3、将偶发拥堵的上游路段La和常发性拥堵的上游路段Lb结合,得到整个网络的拥堵上游路段L。
步骤S2、确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的连续优化问题,并构建用于计算目标函数值的中观交通仿真模型,具体为:
决策变量为n个上游路段在途换路的用户比例。目标函数为***总出行时间的期望f(x)最小。由于仿真结果中含有较大的随机因素,且f(x)无法通过显式表达,所以本实施例将p次重复仿真所得***总出行时间的平均值y(x)作为***总出行时间期望f(x)的近似估计值,将原目标函数转化为近似目标函数。优化数学模型具体如下:
决策变量:
x=[x1,x2,...,xn]T
其中1,2,...,n表示拥堵上游发布ATIS信息的路段,即S1中确定的拥堵上游路段L,xi表示第i个路段上在途换路用户的比例。
原目标函数:
其中f(x)为***总出行时间的期望;为***总出行时间,包含随机因素,为一个随机变量;r,R分别为用户出行的起点节点及所有起点节点的集合;s,S分别为用户出行的终点节点及所有终点节点的集合;t,T分别为用户的出发时刻及所有出发时刻的集合;为高峰时段t时刻从起点r到终点s的总出行量;为高峰时段t时刻从起点r到终点s的平均出行时间,与x有关;为非高峰时段t时刻从起点r到终点s的总出行量;为非高峰时段t时刻从起点r到终点s的平均出行时间,与x无关。
近似目标函数:
约束条件:
xmin≤x≤xmax
其中xmin=[0,0,...,0]T,为在途换路用户比例的下限;xmax=[1,1,...,1]T,为在途换路用户比例的上限。
步骤S3、基于拉丁超立方抽样在解空间中均匀抽取样本,并带入所述中观交通仿真模型计算样本的目标函数值,建立仿真样本集合,具体为:
基于拉丁超立方抽样在决策变量的解空间中均匀采样,得到样本点集 {x1,x2,...,xd},利用式(3)获得每个样本点对应的目标函数值{y(x1),y(x2),...,y(xd)},由此建立初始仿真样本集合D:
D={[x1,y(x1)],[x2,y(x2)],...,[xd,y(xd)]}
步骤S4、利用仿真样本集合,建立基于偏最小二乘Kriging模型的代理模型,具体为:
步骤S4.1、利用仿真样本集合,建立仿真模型的基于偏最小二乘Kriging模型的代理模型,模型形式如下:
fm(x)=μ+ε(x)
式中,μ是目标函数的平均值;ε(x)为估计误差,是一个均值为0,方差为σ2的高斯过程,满足如下的协方差关系:
Cov[ε(x(i)),ε(x(j))]=σ2ψ(x(i),x(j))
利用偏最小二乘技术将决策空间的维度压缩到可接受的范围。整合偏最小二乘的降维结果和Kriging的相关矩阵,得到如下的偏最小二乘Kriging基函数ψ:
其中θ=[θ1,θ2,...,θh]T为缩放系数向量,可以使得n维决策变量每个维度的基函数具体不同的宽度。PLS降维后的变量个数h远小于原变量个数n,模型中需要估计的参数个数减少,仿真优化算法的效率将显著提升。
基于仿真样本集合D中的观测数据可以得到如下的相关矩阵:
Ri,j=ψ(x(i),x(j))
式中:
r={Cov[ε(x*),ε(x(1))],...,Cov[ε(x*),ε(x(d))]}T
y={y(x1),y(x2),...,y(xd)}T
步骤S4.2、代理模型的初始值设置。
建立偏最小二乘Kriging代理模型需要设置初始参数值θ0。利用网格搜索法和交叉验证法确定初始参数值的最佳取值,步骤如下:
1)将仿真样本集合D随机分为K份,得到子样本集合:D1,D2,...,DK。
2)将θ0的取值网格化,Θ={0.001,0.01,0.1,1,10}。
3)将Θ中的值逐一赋给θ0,使用K折交叉验证计算代理模型的平均预测误差。
4)选择Θ中平均预测误差最小的参数作为θ0的初始参数值。
步骤S5、基于代理模型,利用群遗传算法和改善期望函数在解空间内进行二次采样,获得最大价值采样点,具体为:
S5.1:利用代理模型辅助的进化算法进行搜索。基于仿真模型的目标函数具有非线性、计算成本高的特点,因此各种代理模型辅助的进化算法常被用来求解该类优化问题。本方法使用代理模型辅助的多种群遗传算法进行寻优,引入不同的种群同时寻优,可以提升算法在解空间中的搜索能力。
S5.2:由于代理模型是对仿真模型的有限近似,故代理模型对目标函数值的预测结果通常存在误差。采集函数可用来评估二次采样阶段得到样本点的潜在价值。选择样本点改善函数I(x)和样本点改善期望函数E[I(x)]作为采集函数,改善期望值可以描述该样本点比仿真样本集合D中最优解ymin更优的期望。
I(x)=ymin-y(x)
步骤S6、将最大价值采样点带入所述中观交通仿真模型中计算其目标函数值,并更新仿真样本集合;
步骤S7、判断是否达到最大迭代次数,若否,则返回步骤S4,否则结束输出当前得到的最优解输出当前得到的最优解,即各拥堵上游路段所需的最优在途换路用户比例。
接下来,选择A区路网作为实施例验证对象,对本发明的方法进行验证。
该路网面积约60平方千米,包含A中环高架(西南侧)、内环高架(西南侧)、沪闵高架三条快速路,及所有高架上下匝道和地面主干道,共有134个交通小区、 372个节点、962条路段(如图2)。设置如下的拥堵场景:A、B、C三处同时发生了不同严重程度的交通事故,事故所在路段的通行能力分别下降100%、50%、 25%,事故的存在时间为上午7:00~8:00;D、E、F为路网中存在的三个常发性拥堵节点,高峰时段(上午7:00~10:00)这些节点附近的通行效率将明显下降。为了引导拥堵路段上游的用户避开拥堵重新规划路线,在相关区域共设置了18个布有 ATIS的拥堵警告区(G~X),高峰时段经过该区域的部分用户将收到前方的实时路况信息。ATIS设备包括但不限于可变信息标志(VMS)、导航***、车载无线通讯***,所有可以将前方即时拥堵信息传递给上游的通讯方式均可被认为是 ATIS设备。
步骤S1中确定常发性拥堵节点或路段的拥堵上游路段过程为:
利用DynusT中观交通仿真软件进行交通仿真,可得在没有路径诱导方案的情况下各路段在仿真时间段内的车流速度。通过路段速度损失率可以判断出在仿真过程中的车流拥堵程度,选取其中损失率最大的100条路段及其上游路段(共计371 条)作为初始决策空间进一步分析,分别求出其在10%、20%、30%、40%、50%、 60%、70%、80%、90%、100%在途换路用户比例下的***总出行时间减少量。
利用KMO和巴特利特检验进行因子分析前的数据适用性验证。分析结果如表 1所示,根据分析结果可知KMO值为0.574大于0.5,巴特利特球形度检验显著性为0.000小于0.05,说明适合进行主成分分析。
表1
表2为Y1-Y8八个指标的方差贡献率可知,第一、二、三主成分方差贡献率分别为58.178%、24.967%、8.551%,累计方差贡献率达91.696%。第一主成分(Z1) 主要反映平均值、总和、均值、极差等综合性数据,第二主成分(Z2)主要反映最大值和上四分位数等数值较大数据,第三主成分(Z3)主要反映最小值和下四分位数等数值较小数据,三个主成分各自具有代表性。
表2
由主成分得分系数矩阵可得各主成分得分,同时以主成分方差贡献率为权重,构建综合评价得分函数:Z=0.534×Z1+0.233×Z2+0.233×Z3,得出各路段综合得分和排名,其中综合得分最高的前十条路段如表3所示。
表3
按照PCA处理得到的综合评价得分由高到低依次将路段加入到候选的拥堵上游路段集合中,若待加入的候选路段与已存在的候选路段间最短间距小于最小间距阈值(本实施例设为500m)时,则跳过该路段继续进行,直至候选集合大小达到目标容量。得到的常发性拥堵的上游路段Lb如图3所示。图中实线为经过路段筛选后最终确认的拥堵上游路段Lb,虚线为不满足最小间距阈值设定而被筛选掉的路段。
步骤S1中确定偶发性拥堵节点或路段的拥堵上游路段过程为:将所有偶发性拥堵节点或路段上游上下匝道附近的路段作为偶发拥堵的上游路段La。
如图2所示,将偶发拥堵的上游路段La和常发性拥堵的上游路段Lb结合,得到整个网络的拥堵上游路段L,具体共包含以下的18条路段:G、H、I、J、K、L、 M、N、O、P、Q、R、S、T、U、V、W、X。
步骤S2中定义的优化问题如下:
决策变量:
x=[x1,x2,...,xn]T
原目标函数:
近似目标函数:
约束条件:
xmin≤x≤xmax
其中n的值为18,p的值为3。
步骤S2中建立徐汇区路网仿真模型的过程如下:
利用粒子滤波算法从徐汇区卡口电警数据中重构出机动车轨迹数据,经过高架线圈数据的修正(动态OD反推),共得到约23万条机动车轨迹数据。图4为标定需求后的仿真结果与实测路段流量的比较,可以看出,需求标定结果较理想。基于上述供给和需求数据建立DynusT中观仿真模型,仿真时间段为上午6:00~11:00。 DynusT是目前应用较为广泛的基于中观交通仿真的动态交通分配软件,本实施例中仅使用其中的中观仿真功能。假设所有用户均为en-route类型,即所有车辆均支持在途选路。同时,本实施例对地面主干道和快速路的交通流模型分别进行了标定,标定结果如表4。
本实施例在DynusT中观仿真模型中选用形如下式的交通流模型:
其中vi当前速度,ki为当前密度,v0为路段的最低速度,vf为路段的自由流速度,kjam为拥挤阻塞密度,α、β分别为需要标定的模型参数。图5和图6分别为地面主干道与快速路的交通流模型的标定结果,其中散点代表观测值,虚线代表标定得到的交通流模型,可以看出,交通流模型的标定结果较理想。
表4
α | β | |
地面主干道 | 2.40 | 0.80 |
快速路 | 4.52 | 1.00 |
步骤S3中基于拉丁超立方抽样共获得了90个初始样本,将这些样本带入仿真模型进行计算其目标函数值,最终可以得到仿真样本集合。
步骤S4中基于仿真样本集合建立了偏最小二乘Kriging的代理模型。本实施例选择了以下五个指标(均方根误差(RMSE)、归一化均方根误差(NRMSE)、最大误差(MAE)、归一化最大误差(NMAE)、拟合代理模型所需的计算时间(CPU)) 来比较偏最小二乘Kriging代理模型和普通Kriging代理模型的预测精度和计算效率差异,比较结果如表5。可以看出,与普通Kriging模型相比,偏最小二乘Kriging 模型预测结果的均方根误差降低近20%,最大误差降低25%,且所消耗的计算资源仅为KRG的4%。偏最小二乘Kriging模型在预测精度和计算效率上都明显优于普通Kriging模型。
表5
代理模型 | RMSE | MAE | NRMSE | NMAE | CPU |
KRG | 40 814.05 | 68 697.42 | 1.27% | 2.58 | 63.44s |
KPLS | 32 832.63 | 42 544.07 | 1.02% | 1.93 | 2.53s |
注:CPU在i9-10900 2.81GHz的计算机上测得。
步骤S5中利用多种群遗传算法和改善期望函数在解空间内进行二次采样,获得最大价值的采样点。本实施例中遗传算法的细节如下:采用采用二进制格雷编码,决策变量的编码精度精确到小数点后2位,种群选择过程为随机选择,DNA交叉方式为两点交叉,采用二进制染色体的变异算子。
由于本发明使用了启发式算法,求解过程具有一定的随机性,因此本实施例运行了该算法三次,最终得到了三种路径诱导***的优化设计方案,如表6所示。
表6
表7整理了三种优化方案和基准方案(无路径诱导,所有的在途换路用户比例均为0)下的网络通行效率。与基准方案相比,最优方案可以使网络中***总出行时间节省约3.9×105分钟,平均每位用户在网络中的通行时间减少约1.8分钟,整个网络的通行效率提升近11%。三种方案的效果十分类似,仿真优化方法的鲁棒性较强。
表7
图7、图8、图9分别比较了在最优方案和基准方案下网络在高峰时段的宏观性能差异,其中图7为路段加权平均速度差异、图8为路段加权平均行程时间差异、图9为网络中排队车辆总数的差异。可以看出与基准方案相比,高峰时段内最优方案的平均速度更大、行程时间和网络中的排队车辆数更小。在事故发生的时段前后 (7:30~8:30),两种方案的差异最明显。可见最优方案可以显著缓解高峰期交通拥堵。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (10)
1.一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,该方法包括以下步骤:
步骤S1、基于实际的交通供给、需求及突发交通状况信息,确定需要发送“前方拥堵警告”信息的拥堵上游路段;
步骤S2、确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的连续优化问题,并构建用于计算目标函数值的中观交通仿真模型;
步骤S3、基于拉丁超立方抽样在解空间中均匀抽取样本,并带入所述中观交通仿真模型计算样本的目标函数值,建立仿真样本集合;
步骤S4、利用仿真样本集合,建立基于偏最小二乘Kriging模型的代理模型;
步骤S5、基于代理模型,利用群遗传算法和改善期望函数在解空间内进行二次采样,获得最大价值采样点;
步骤S6、将最大价值采样点带入所述中观交通仿真模型中计算其目标函数值,并更新仿真样本集合;
步骤S7、判断是否达到最大迭代次数,若否,则返回步骤S4,否则结束输出当前得到的最优解。
2.根据权利要求1所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述步骤S1中的拥堵上游路段包括常发性拥堵上游路段和偶发性拥堵上游路段;
所述偶发性拥堵上游路段采用直接确定法进行确定;所述直接确定法为将偶发性拥堵上游的上下匝道设定范围内路段作为偶发性拥堵上游路段。
3.根据权利要求2所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述常发性拥堵上游路段采用主成分分析法进行确定;所述主成分分析法包括:路段选取、路段降维以及路段筛选。
5.根据权利要求3所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述路段降维为使用主成分分析法对路段选取结果Lb,1进行降维,具体包括以下子步骤:
1)将样本数据标准化,得到标准化阵Z,表达式为:
式中,zij为标准化阵Z的元素,分别为变量Yj对应样本的均值和标准差;变量Yj包括:各在途换路用户比例下***总出行时间减少量平均值Y1、***总出行时间减少量的最小值Y2、***总出行时间减少量的最大值Y3、***总出行时间减少量的上四分位数Y4、***总出行时间减少量的中位数Y5、***总出行时间减少量的下四分位数Y6、各在途换路用户比例下***总出行时间减少量的极差Y7、以及各在途换路用户比例下***总出行时间减少量总和Y8;
2)计算标准化阵Z的协方差矩阵R,表达式为
式中,rij为协方差矩阵R的元素,n为样本容量:
3)计算协方差矩阵R的特征值λ1≥λ2≥…≥λn≥0以及对应特征向量u1,u2,…,un;
4)计算各主成分yi的贡献率bj以及y1,y2,…,yp累计贡献率αp:
选取特征值大于1且累计贡献率约达到90%的成分作为主成分;
5)获得成分得分系数矩阵,计算各成分得分,之后通过旋转后的方差百分比与各成分得分求出最后的综合评价得分,得到降维后的上游路段Lb,2。
6.根据权利要求3所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述路段筛选具体为:引入候选路段的最短间距限制,从降维后的上游路段Lb,2中进行筛选,将筛选后的结果作为常发性拥堵上游路段Lb。
7.根据权利要求1所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述步骤S2中确定以拥堵上游路段的在途换路用户比例为决策变量、以最小化***总出行时间为优化目标的连续优化问题,并构建用于计算目标函数值的中观交通仿真模型,具体为:
1)定义决策变量为n个上游路段在途换路的用户比例,目标函数为***总出行时间的期望f(x)最小;其中,所述决策变量表达式为:
x=[x1,x2,...,xn]T
式中,1,2,...,n表示拥堵上游发布ATIS信息的路段,即步骤S1中确定的拥堵上游路段L,xi表示第i个路段上在途换路用户的比例;
2)将p次重复仿真所得***总出行时间的平均值y(x)作为***总出行时间期望f(x)的近似估计值,将原目标函数转化为近似目标函数,优化数学模型具体如下:
所述原目标函数表达式为:
式中,f(x)为***总出行时间的期望;F(x,ζ)为***总出行时间,包含一个随机变量;r,R分别为用户出行的起点节点及所有起点节点的集合;s,S分别为用户出行的终点节点及所有终点节点的集合;t,T分别为用户的出发时刻及所有出发时刻的集合;为高峰时段t时刻从起点r到终点s的总出行量;为高峰时段t时刻从起点r到终点s的平均出行时间,与x有关;为非高峰时段t时刻从起点r到终点s的总出行量;为非高峰时段t时刻从起点r到终点s的平均出行时间,与x无关;
所述近似目标函数表达式为:
约束条件为:
xmin≤x≤xmax
式中,xmin=[0,0,...,0]T为在途换路用户比例的下限;xmax=[1,1,...,1]T为在途换路用户比例的上限。
8.根据权利要求7所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述步骤S4具体为:
步骤S4.1、构建基于偏最小二乘Kriging模型的代理模型,表达式为:
fm(x)=μ+ε(x)
式中,μ是目标函数的平均值;ε(x)为估计误差,服从均值为0、方差为σ2的高斯分布,满足以下协方差关系:
Cov[ε(x(i)),ε(x(j))]=σ2ψ(x(i),x(j))
式中,x(i)为第i个路段上在途换路用户的比例;ψ为偏最小二乘Kriging基函数,表达式为:
式中,R为基于仿真样本集合D中的观测数据得到的相关矩阵,表达式为:
Ri,j=ψ(x(i),x(j));
r为协方差向量,表达式为:
r={Cov[ε(x*),ε(x(1))],...,Cov[ε(x*),ε(x(d))]}T
步骤S4.2、采用网格搜索法和交叉验证法确定代理模型的初始参数θ0的最佳取值。
9.根据权利要求7所述的一种基于偏最小二乘Kriging模型的路径诱导方法,其特征在于,所述步骤S4.2具体为:
1)将仿真样本集合D随机分为K份,得到子样本集合:D1,D2,...,DK;
2)将θ0的取值网格化,Θ={0.001,0.01,0.1,1,10};
3)将Θ中的值逐一赋给θ0,使用K折交叉验证计算代理模型的平均预测误差;
4)选择Θ中平均预测误差最小的参数作为θ0的初始参数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210405251.2A CN115482655B (zh) | 2022-04-18 | 2022-04-18 | 一种基于偏最小二乘Kriging模型的路径诱导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210405251.2A CN115482655B (zh) | 2022-04-18 | 2022-04-18 | 一种基于偏最小二乘Kriging模型的路径诱导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115482655A true CN115482655A (zh) | 2022-12-16 |
CN115482655B CN115482655B (zh) | 2023-08-29 |
Family
ID=84420795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210405251.2A Active CN115482655B (zh) | 2022-04-18 | 2022-04-18 | 一种基于偏最小二乘Kriging模型的路径诱导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115482655B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502566A (zh) * | 2023-06-27 | 2023-07-28 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 基于贝叶斯优化的燃气轮机燃烧室性能多目标优化方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101633360B1 (ko) * | 2015-01-12 | 2016-06-27 | 한국과학기술원 | 구조물 상태 평가를 위한 크리깅 모델 기반 순차적 샘플링 방법 |
CN109409560A (zh) * | 2018-08-16 | 2019-03-01 | 北京交通大学 | 基于多智能体仿真的城市轨道交通客流诱导方法 |
WO2021109318A1 (zh) * | 2019-12-03 | 2021-06-10 | 东南大学 | 一种城市路网短期交通运行状态估计与预测方法 |
WO2021142916A1 (zh) * | 2020-01-15 | 2021-07-22 | 深圳大学 | 基于代理辅助进化算法的翼型优化方法及装置 |
CN114048547A (zh) * | 2021-11-18 | 2022-02-15 | 汉思科特(盐城)减震技术有限公司 | 一种基于自适应代理模型的车用空气弹簧工程优化设计方法 |
-
2022
- 2022-04-18 CN CN202210405251.2A patent/CN115482655B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101633360B1 (ko) * | 2015-01-12 | 2016-06-27 | 한국과학기술원 | 구조물 상태 평가를 위한 크리깅 모델 기반 순차적 샘플링 방법 |
CN109409560A (zh) * | 2018-08-16 | 2019-03-01 | 北京交通大学 | 基于多智能体仿真的城市轨道交通客流诱导方法 |
WO2021109318A1 (zh) * | 2019-12-03 | 2021-06-10 | 东南大学 | 一种城市路网短期交通运行状态估计与预测方法 |
WO2021142916A1 (zh) * | 2020-01-15 | 2021-07-22 | 深圳大学 | 基于代理辅助进化算法的翼型优化方法及装置 |
CN114048547A (zh) * | 2021-11-18 | 2022-02-15 | 汉思科特(盐城)减震技术有限公司 | 一种基于自适应代理模型的车用空气弹簧工程优化设计方法 |
Non-Patent Citations (2)
Title |
---|
张干锋;王德禹;: "基于分区间泛化Kriging近似模型的船舶局部结构形状优化", 海洋工程, no. 06, pages 48 - 56 * |
张敖木翰;高自友;: "突发事故下交通信息策略研究", 交通运输***工程与信息, no. 03 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502566A (zh) * | 2023-06-27 | 2023-07-28 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 基于贝叶斯优化的燃气轮机燃烧室性能多目标优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115482655B (zh) | 2023-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | A machine learning method for predicting driving range of battery electric vehicles | |
CN107610464B (zh) | 一种基于高斯混合时间序列模型的轨迹预测方法 | |
CN111301426B (zh) | 基于gru网络模型预测未来行驶过程能耗的方法 | |
CN111461410B (zh) | 一种基于迁移学习的空气质量的预测方法及装置 | |
Zheng et al. | A stochastic simulation-based optimization method for equitable and efficient network-wide signal timing under uncertainties | |
CN114781538B (zh) | 一种ga-bp神经网络耦合决策树的空气质量预测方法和*** | |
CN113112090B (zh) | 基于综合互信息度的主成分分析的空间负荷预测方法 | |
CN109376331A (zh) | 一种基于梯度提升回归树的城市公交车排放率估计方法 | |
CN114048920A (zh) | 充电设施建设的选址布局方法、装置、设备及存储介质 | |
Lee et al. | Dynamic OD estimation using three phase traffic flow theory | |
CN115482655B (zh) | 一种基于偏最小二乘Kriging模型的路径诱导方法 | |
Sun et al. | Bi-objective traffic count location model for mean and covariance of origin–destination estimation | |
CN116046008A (zh) | 基于态势感知的路线规划方法、***及效能评估装置 | |
CN113642162B (zh) | 一种城市道路交通应急预案仿真综合分析方法 | |
Fahs et al. | Traffic congestion prediction based on multivariate modelling and neural networks regressions | |
Festa et al. | Experimental analysis of different simulation models for motorway traffic flow | |
KR102601046B1 (ko) | 차량 가용도 예측 방법 및 장치 | |
Dong et al. | An identification model of urban critical links with macroscopic fundamental diagram theory | |
CN115691140A (zh) | 一种汽车充电需求时空分布的分析与预测方法 | |
CN116434566A (zh) | 交通流量的确定方法及设备 | |
Haragos et al. | Modelling road traffic using service center | |
CN117172516B (zh) | 充电桩动态调度决策方法、装置、设备及存储介质 | |
CN112215520B (zh) | 一种多目标融合通行方法、装置、计算机设备及存储介质 | |
CN117894174B (zh) | 一种道路交通安全主动管控方法、装置及*** | |
CN117723999B (zh) | 电池使用寿命预测方法、装置、设备及存储介质 |
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 |