CN108985500A - 基于改进模拟退火算法的城际列车开行方案优化方法 - Google Patents
基于改进模拟退火算法的城际列车开行方案优化方法 Download PDFInfo
- Publication number
- CN108985500A CN108985500A CN201810687399.3A CN201810687399A CN108985500A CN 108985500 A CN108985500 A CN 108985500A CN 201810687399 A CN201810687399 A CN 201810687399A CN 108985500 A CN108985500 A CN 108985500A
- Authority
- CN
- China
- Prior art keywords
- train
- inter
- solution
- city
- passenger
- 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.)
- Pending
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000002922 simulated annealing Methods 0.000 claims abstract description 16
- 230000008901 benefit Effects 0.000 claims abstract description 11
- 238000013459 approach Methods 0.000 claims abstract description 7
- 238000012790 confirmation Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 56
- 238000012546 transfer Methods 0.000 claims description 17
- 230000003044 adaptive effect Effects 0.000 claims description 14
- 230000000977 initiatory effect Effects 0.000 claims description 9
- 230000006386 memory function Effects 0.000 claims description 8
- 238000001816 cooling Methods 0.000 claims description 7
- 238000011160 research Methods 0.000 claims description 7
- 238000003672 processing method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 7
- 238000000137 annealing Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 2
- 238000005496 tempering Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000155 melt Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
本发明公开了基于改进模拟退火算法的城际列车开行方案优化方法,主要包括:步骤S1:构建多目标函数,并建立开行方案优化模型;步骤S2:采用线性加权法将多目标函数转化为单目标函数;步骤S3:确认城际列车开行的相关影响因子;步骤S4:利用改进的模拟退火算法进行列车开行方案的优化计算;步骤S5:根据计算结果,提取最优目标函数值。本发明的基于改进模拟退火算法的城际列车开行方案优化方法,以实现以期为城际列车开行方案的优化提供一种新的方法,同时改善传统模拟退火算法求解时时间过长、初始温度较难确定以及搜索过程中因执行概率接受准则而导致遗失当前最优解等问题,可以更容易地找出全局最优解的优点。
Description
技术领域
本发明涉及城际列车开行方案优化方法领域,具体地,涉及基于改进模拟退火算法的城际列车开行方案优化方法。
背景技术
随着近年来我国经济的快速发展,城镇化进程的日益加快,城市已经不再是一个独立的个体,各城镇间物资、人才交流日益密切,城市的功能定位越来越明确,城市的发展也越来越依赖于城市群的推动。城市群的发展所带来的交通需求特征变化主要体现在:城市对外交通组织具有明显的区域化特征;城镇间的交通需求具有明显趋向于城市交通的特征;城市群旅客出行空间的分布体现了区域智能化的趋势。
城际快速轨道交通系统是我国铁路客运专线规划和建设的一个重要内容,许多城市群的城际铁路已经开通运营。城际快速轨道交通是建设在人口稠密、经济发达地区的公交化开行的,专门服务于相邻城市或者城市群的旅客运输专线铁路,它是区域性的重要交通运输基础设施。城际轨道交通所包涵的不只是城铁,还有地铁、轻轨、单轨等。多种出行方式相互补充,共同提高了城市群综合交通系统的服务质量和运行效率。
列车开行方案是客运组织的核心,是运行图编制的基础。为了更好的服务旅客和分配运力资源,制定科学高效的列车开行方案就显得尤为重要。开行方案旨在确定列车的开行数量、运行径路以及停站序列等要素。合理制定列车开行方案能最大化利用铁路既有设备,提高运输效益和旅客服务水平。
从既有的相关研究来看,现阶段关于城际铁路列车开行方案和地铁列车开行方案的研究已较为成熟。但是,目前既有研究在制定列车开行方案时均只考虑了铁路系统内部因素,未考虑旅客到站后与城市交通的匹配问题。实际上,对于时间价值较重的城际客流,单独考虑城际铁路系统内部制定开行方案可能会导致旅客到站后无法及时疏散,从而影响整体的服务质量和铁路运输效益。
模拟退火算法是80年代初期发展起来的一种求解大规模组合优化问题的随机性方法,能解决传统的优化方法难于解决的某些问题,在列车开行方案的优化问题中同样适用。模拟退火算法以优化问题的求解与物理系统退火过程的相似性为基础,利用Metropolis算法并适当地控制温度的下降过程实现模拟退火,从而达到求解全局优化问题的目的。
本发明即是在对既有研究进行分析的基础上,发现当前研究存在的问题,提出了一种基于改进模拟退火算法的城际列车开行方案优化方法。
发明内容
本发明的目的在于,针对上述问题,提出基于改进模拟退火算法的城际列车开行方案优化方法,以实现以期为城际列车开行方案的优化提供一种新的方法,同时改善传统模拟退火算法求解时时间过长、初始温度较难确定以及搜索过程中因执行概率接受准则而导致遗失当前最优解等问题,可以更容易地找出全局最优解的优点。
为实现上述目的,本发明采用的技术方案是:基于改进模拟退火算法的城际列车开行方案优化方法,主要包括:
步骤S1:构建多目标函数,并建立开行方案优化模型;
步骤S2:采用线性加权法将多目标函数转化为单目标函数;
步骤S3:确认城际列车开行的相关影响因子;
步骤S4:利用改进的模拟退火算法进行列车开行方案的优化计算;
步骤S5:根据计算结果,提取最优目标函数值。
进一步地,步骤S1中所述构建多目标函数,具体包括:
以f(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的客流量,c(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均广义出行费用;第一个优化目标总广义出行费用最小可以表示为:
票价费用和旅行时间费用共同构成旅客平均广义出行费用;OD(OriginalDestination)为起始站至终到站;
以c1表示旅客出行成本中票价的平均权重因子;c2表示旅行时间转化为费用的平均权重因子;p(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均票价;dsisj表示OD对(Si,Sj)间的路段距离;以vTu来表示考虑列车停站和起停附加时分在内的列车Tu的平均旅行速度;为0-1变量,列车Tu分担OD对(Si,Sj)间的客流量记作否则平均广义出行费用c(Si,Sj,Tu)表示为:
铁路运输部门的总收入主要来源于客票收入;运营成本由固定成本和可变成本两部分构成,将固定成本视作开行一列车必须产生的费用,可变成本由列车公里费用、车辆公里费用和停站费用三部分组成;
以为表示列车停站的0-1变量,若列车Tu在车站Sm停车记作否则cfTu表示开行列车Tu的固定成本;cpTu表示列车Tu的停站费用;cuTu为列车Tu的车公里费用;cvTu为车辆公里费用;mTu为列车Tu的编组数量;dSoSd为列车自始发站So运行到终到站Sd走行的总公里数;f T为T类列车的开行数量;
运营总收入可表示为:
运营总成本可表示为:
第二个优化目标铁路运输部门的效益最大可以表示为:
maxz2=CI-Cs。
进一步地,步骤S1中所述建立模型,具体包括:
以列车载客量、客流需求的满足、通过能力及运能匹配度作为约束条件建立模型;
约束条件的限制主要包括:
1)列车Tu输送的总客流量不大于此列车的载客量的限制:
式中,Au表示列车Tu的定员;Ψu表示列车Tu的满载率;
2)服务于OD对(Si,Sj)的所有列车的载客能力要满足此OD对间的客流需求:
式中,f(Si,Sj)表示OD对(Si,Sj)间的客流需求;
3)研究时段内所有列车开行频率之和不超过区段通过能力Nei和车站通过能力Nsi能力的限制:
4)所有列车在区段两端站,即列车的始发站和终到站必须停车:
5)若城际铁路与城市轨道交通的运能可以有效衔接,城市轨道交通可以提供的运能应满足及时疏散换乘客流的需求;
假设在客流换乘的某时段h内,T类城际列车到达终到站的列车数为[fT·h/H]([]表示取整),H表示城际列车总运营时段长度;城际列车终到站可换乘城市轨道交通的所有线路共有nr个方向;城市轨道交通列车在研究时段h内的发车间隔为tr;列车编组为mr;车厢最大载客量为cr;到达终到站后的换乘客流所占比例为α;时段h内城市轨道交通列车可供换乘旅客乘坐的容量比例为β;
经上述分析,在时段h内的换乘客流量可表示为:
若要保证城市轨道交通列车的运能能够及时疏散换乘旅客,必须满足:
进一步地,步骤S2中所述采用线性加权法将多目标函数转化为单目标函数,具体包括:
为了便于后续求解,以ε表示目标函数1的权重因子,1-ε表示目标函数2经最小化处理后的权重因子,即可将多目标函数转化为单目标函数。
进一步地,步骤S3中所述相关影响因子,包括:
城际线路的区间里程、列车开行区段、通过能力、票价率以及与之衔接的城市轨道交通车站和列车的相关信息。
进一步地,步骤S4中所述利用改进的模拟退火算法进行列车开行方案的优化计算,具体包括:
步骤S41:设置初始化参数:
步骤S42:随机产生扰动,得到新解X2并计算新解的目标函数值f(X2);计算目标函数值的增量Δf=f(X2)-f(X1);
步骤S43:判断新解X2是否满足约束,若满足则转步骤S44;否则转步骤S46;
步骤S44:根据Metropolis准则判断是否接受新解:
步骤S45:执行记忆功能:
步骤S46:令k=k+1,若k≤L,则转D2;否则转D7;
步骤S47:采用多普勒型降温函数:
T=T0αk(cos(π/(2(1-k/K))))+cos(π/(2T0(1-k/K))),
将温度进行衰减;
步骤S48:算法终止检验。
进一步地,步骤S41中所述设置初始化参数,具体包括:
设置初始温度T0,终止温度Te,每个温度T时的迭代次数L,记忆矩阵M,令当前温度T=T0,温度下降比例θ,令k=1,p=0;
初始温度T0的设定,提出一种自适应的处理方法:随机产生100个满足约束的解,并判断有无可行解;如果有可行解,则设置初始温度为5000,从可行解中选出函数值最小的作为模拟退火算法的初始解X0;否则,初始温度设置为50000,从不可行解中选出罚函数值最小的作为初始解X0;令当前解X1=X0,计算当前解的目标函数值f(X1);
引入一种自适应的函数变换方法:在初始化过程中,记录并计算函数的平均值和最小值,确定一个函数的自适应因子K:有可行解时,
否则,
这样就可以通过自适应因子K来调节解的接受准则,从而增强了算法的健壮性。
进一步地,步骤S44中所述根据Metropolis准则判断是否接受新解,具体包括:
若Δf<0或Δf≥0且满足exp(-Δf/T)>ξ,ξ为区间(0,1)上产生的均匀分布的随机数,则接受X2作为新的当前解,即X1=X2;否则保留当前解X1。
进一步地,步骤S45中所述执行记忆功能,具体包括:
若迭代次数为1,则将当前产生的最好解记录在记忆矩阵M之中;否则将当前产生的最好解与记忆矩阵M中的解进行比较,如果当前解优于记忆矩阵中的解,则将当前解写入记忆矩阵M中,否则保留记忆矩阵中的解。
进一步地,步骤S48中所述算法终止检验,具体包括:
若T<Te,则输出当前解X1,结束程序;否则令k=1,转D2。
本发明的有益技术效果:
1)能够兼顾旅客和铁路运输企业双方的利益,取得更好的社会经济效益;
2)考虑了城际列车到站后的客流与城市轨道交通运能的匹配度问题,就城市群区域范围内而言使列车开行方案更加的合理;
3)对于初始温度的确定,提出一种自适应的处理方法,解决了初始温度较难确定的问题,使算法具有较强的适应性和通用性;
4)在算法搜索过程中记住中间最优解并及时更新,增加记忆功能,可以避免算法在降温过程中遗失掉最优解;
5)降温函数采用多普勒型温度递减函数,可以较好的结合指数降温方式和快速降温方式的优点并消除各自的缺陷,使趋于低温的速度保持适中;同时不断进行的回火升温过程,可以使算法在优化的过程中多次跳出局部最优解,有效避免了传统模拟退火算法极易陷入局部极小值的缺陷。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明所述基于改进模拟退火算法的城际列车开行方案优化方法的流程图;
图2为本发明实施例中改进的模拟退火算法流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
如图1所示,基于改进模拟退火算法的城际列车开行方案优化方法,主要包括:
步骤S1:构建多目标函数,并建立开行方案优化模型;
步骤S2:采用线性加权法将多目标函数转化为单目标函数;
步骤S3:确认城际列车开行的相关影响因子;
步骤S4:利用改进的模拟退火算法进行列车开行方案的优化计算;
步骤S5:根据计算结果,提取最优目标函数值。
进一步地,步骤S1中所述构建多目标函数,具体包括:
以f(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的客流量,c(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均广义出行费用;第一个优化目标总广义出行费用最小可以表示为:
票价费用和旅行时间费用共同构成旅客平均广义出行费用;OD(OriginalDestination)为起始站至终到站;
以c1表示旅客出行成本中票价的平均权重因子;c2表示旅行时间转化为费用的平均权重因子;p(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均票价;dsisj表示OD对(Si,Sj)间的路段距离;以vTu来表示考虑列车停站和起停附加时分在内的列车Tu的平均旅行速度;为0-1变量,列车Tu分担OD对(Si,Sj)间的客流量记作否则平均广义出行费用c(Si,Sj,Tu)表示为:
铁路运输部门的总收入主要来源于客票收入;运营成本由固定成本和可变成本两部分构成,将固定成本视作开行一列车必须产生的费用,可变成本由列车公里费用、车辆公里费用和停站费用三部分组成;
以为表示列车停站的0-1变量,若列车Tu在车站Sm停车记作否则cfTu表示开行列车Tu的固定成本;cpTu表示列车Tu的停站费用;cuTu为列车Tu的车公里费用;cvTu为车辆公里费用;mTu为列车Tu的编组数量;dSoSd为列车自始发站So运行到终到站Sd走行的总公里数;f T为T类列车的开行数量;
运营总收入可表示为:
运营总成本可表示为:
第二个优化目标铁路运输部门的效益最大可以表示为:
maxz2=CI-Cs。
进一步地,步骤S1中所述建立模型,具体包括:
以列车载客量、客流需求的满足、通过能力及运能匹配度作为约束条件建立模型;
约束条件的限制主要包括:
1)列车Tu输送的总客流量不大于此列车的载客量的限制:
式中,Au表示列车Tu的定员;Ψu表示列车Tu的满载率;
2)服务于OD对(Si,Sj)的所有列车的载客能力要满足此OD对间的客流需求:
式中,f(Si,Sj)表示OD对(Si,Sj)间的客流需求;
3)研究时段内所有列车开行频率之和不超过区段通过能力Nei和车站通过能力Nsi能力的限制:
4)所有列车在区段两端站,即列车的始发站和终到站必须停车:
5)若城际铁路与城市轨道交通的运能可以有效衔接,城市轨道交通可以提供的运能应满足及时疏散换乘客流的需求;
假设在客流换乘的某时段h内,T类城际列车到达终到站的列车数为[fT·h/H]([]表示取整),H表示城际列车总运营时段长度;城际列车终到站可换乘城市轨道交通的所有线路共有nr个方向;城市轨道交通列车在研究时段h内的发车间隔为tr;列车编组为mr;车厢最大载客量为cr;到达终到站后的换乘客流所占比例为α;时段h内城市轨道交通列车可供换乘旅客乘坐的容量比例为β;
经上述分析,在时段h内的换乘客流量可表示为:
若要保证城市轨道交通列车的运能能够及时疏散换乘旅客,必须满足:
进一步地,步骤S2中所述采用线性加权法将多目标函数转化为单目标函数,具体包括:
为了便于后续求解,以ε表示目标函数1的权重因子,1-ε表示目标函数2经最小化处理后的权重因子,即可将多目标函数转化为单目标函数。
进一步地,步骤S3中所述相关影响因子,包括:
城际线路的区间里程、列车开行区段、通过能力、票价率以及与之衔接的城市轨道交通车站和列车的相关信息。
如图2所示,步骤S4中所述利用改进的模拟退火算法进行列车开行方案的优化计算,具体包括:
D1步骤S41:设置初始化参数:
D2步骤S42:随机产生扰动,得到新解X2并计算新解的目标函数值f(X2);计算目标函数值的增量Δf=f(X2)-f(X1);
D3步骤S43:判断新解X2是否满足约束,若满足则转步骤S44;否则转步骤S46;
D4步骤S44:根据Metropolis准则判断是否接受新解:
D5步骤S45:执行记忆功能:
D6步骤S46:令k=k+1,若k≤L,则转D2;否则转D7;
D7步骤S47:采用多普勒型降温函数:
T=T0αk(cos(π/(2(1-k/K))))+cos(π/(2T0(1-k/K))),
将温度进行衰减;
D8步骤S48:算法终止检验。
进一步地,步骤S41中所述设置初始化参数,具体包括:
设置初始温度T0,终止温度Te,每个温度T时的迭代次数L,记忆矩阵M,令当前温度T=T0,温度下降比例θ,令k=1,p=0;
初始温度T0的设定,提出一种自适应的处理方法:随机产生100个满足约束的解,并判断有无可行解;如果有可行解,则设置初始温度为5000,从可行解中选出函数值最小的作为模拟退火算法的初始解X0;否则,初始温度设置为50000,从不可行解中选出罚函数值最小的作为初始解X0;令当前解X1=X0,计算当前解的目标函数值f(X1);
引入一种自适应的函数变换方法:在初始化过程中,记录并计算函数的平均值和最小值,确定一个函数的自适应因子K:有可行解时,
否则,
这样就可以通过自适应因子K来调节解的接受准则,从而增强了算法的健壮性。
进一步地,步骤S44中所述根据Metropolis准则判断是否接受新解,具体包括:
若Δf<0或Δf≥0且满足exp(-Δf/T)>ξ,ξ为区间(0,1)上产生的均匀分布的随机数,则接受X2作为新的当前解,即X1=X2;否则保留当前解X1。
进一步地,步骤S45中所述执行记忆功能,具体包括:
若迭代次数为1,则将当前产生的最好解记录在记忆矩阵M之中;否则将当前产生的最好解与记忆矩阵M中的解进行比较,如果当前解优于记忆矩阵中的解,则将当前解写入记忆矩阵M中,否则保留记忆矩阵中的解。
进一步地,步骤S48中所述算法终止检验,具体包括:
若T<Te,则输出当前解X1,结束程序;否则令k=1,转D2。
至少可以达到以下有益效果:
1)能够兼顾旅客和铁路运输企业双方的利益,取得更好的社会经济效益;
2)考虑了城际列车到站后的客流与城市轨道交通运能的匹配度问题,就城市群区域范围内而言使列车开行方案更加的合理;
3)对于初始温度的确定,提出一种自适应的处理方法,解决了初始温度较难确定的问题,使算法具有较强的适应性和通用性;
4)在算法搜索过程中记住中间最优解并及时更新,增加记忆功能,可以避免算法在降温过程中遗失掉最优解;
5)降温函数采用多普勒型温度递减函数,可以较好的结合指数降温方式和快速降温方式的优点并消除各自的缺陷,使趋于低温的速度保持适中;同时不断进行的回火升温过程,可以使算法在优化的过程中多次跳出局部最优解,有效避免了传统模拟退火算法极易陷入局部极小值的缺陷。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,主要包括:
步骤S1:构建多目标函数,并建立开行方案优化模型;
步骤S2:采用线性加权法将多目标函数转化为单目标函数;
步骤S3:确认城际列车开行的相关影响因子;
步骤S4:利用改进的模拟退火算法进行列车开行方案的优化计算;
步骤S5:根据计算结果,提取最优目标函数值。
2.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S1中所述构建多目标函数,具体包括:
以f(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的客流量,c(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均广义出行费用;第一个优化目标总广义出行费用最小可以表示为:
票价费用和旅行时间费用共同构成旅客平均广义出行费用;OD(OriginalDestination)为起始站至终到站;
以c1表示旅客出行成本中票价的平均权重因子;c2表示旅行时间转化为费用的平均权重因子;p(Si,Sj,Tu)表示OD对(Si,Sj)间乘坐列车Tu的平均票价;dsisj表示OD对(Si,Sj)间的路段距离;以vTu来表示考虑列车停站和起停附加时分在内的列车Tu的平均旅行速度;为0-1变量,列车Tu分担OD对(Si,Sj)间的客流量记作否则平均广义出行费用c(Si,Sj,Tu)表示为:
铁路运输部门的总收入主要来源于客票收入;运营成本由固定成本和可变成本两部分构成,将固定成本视作开行一列车必须产生的费用,可变成本由列车公里费用、车辆公里费用和停站费用三部分组成;
以为表示列车停站的0-1变量,若列车Tu在车站Sm停车记作否则cfTu表示开行列车Tu的固定成本;cpTu表示列车Tu的停站费用;cuTu为列车Tu的车公里费用;cvTu为车辆公里费用;mTu为列车Tu的编组数量;dSoSd为列车自始发站So运行到终到站Sd走行的总公里数;fT为T类列车的开行数量;
运营总收入可表示为:
运营总成本可表示为:
第二个优化目标铁路运输部门的效益最大可以表示为:
maxz2=CI-Cs。
3.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S1中所述建立模型,具体包括:
以列车载客量、客流需求的满足、通过能力及运能匹配度作为约束条件建立模型;约束条件的限制主要包括:
1)列车Tu输送的总客流量不大于此列车的载客量的限制:
式中,Au表示列车Tu的定员;Ψu表示列车Tu的满载率;
2)服务于OD对(Si,Sj)的所有列车的载客能力要满足此OD对间的客流需求:
式中,f(Si,Sj)表示OD对(Si,Sj)间的客流需求;
3)研究时段内所有列车开行频率之和不超过区段通过能力Nei和车站通过能力Nsi能力的限制:
4)所有列车在区段两端站,即列车的始发站和终到站必须停车:
5)若城际铁路与城市轨道交通的运能可以有效衔接,城市轨道交通可以提供的运能应满足及时疏散换乘客流的需求;
假设在客流换乘的某时段h内,T类城际列车到达终到站的列车数为[fT·h/H]([]表示取整),H表示城际列车总运营时段长度;城际列车终到站可换乘城市轨道交通的所有线路共有nr个方向;城市轨道交通列车在研究时段h内的发车间隔为tr;列车编组为mr;车厢最大载客量为cr;到达终到站后的换乘客流所占比例为α;时段h内城市轨道交通列车可供换乘旅客乘坐的容量比例为β;
经上述分析,在时段h内的换乘客流量可表示为:
若要保证城市轨道交通列车的运能能够及时疏散换乘旅客,必须满足:
4.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S2中所述采用线性加权法将多目标函数转化为单目标函数,具体包括:
为了便于后续求解,以ε表示目标函数1的权重因子,1-ε表示目标函数2经最小化处理后的权重因子,即可将多目标函数转化为单目标函数。
5.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S3中所述相关影响因子,包括:
城际线路的区间里程、列车开行区段、通过能力、票价率以及与之衔接的城市轨道交通车站和列车的相关信息。
6.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S4中所述利用改进的模拟退火算法进行列车开行方案的优化计算,具体包括:
步骤S41:设置初始化参数:
步骤S42:随机产生扰动,得到新解X2并计算新解的目标函数值f(X2);计算目标函数值的增量Δf=f(X2)-f(X1);
步骤S43:判断新解X2是否满足约束,若满足则转步骤S44;否则转步骤S46;
步骤S44:根据Metropolis准则判断是否接受新解:
步骤S45:执行记忆功能:
步骤S46:令k=k+1,若k≤L,则转D2;否则转D7;
步骤S47:采用多普勒型降温函数:
T=T0αk(cos(π/(2(1-k/K))))+cos(π/(2T0(1-k/K))),
将温度进行衰减;
步骤S48:算法终止检验。
7.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S41中所述设置初始化参数,具体包括:
设置初始温度T0,终止温度Te,每个温度T时的迭代次数L,记忆矩阵M,令当前温度T=T0,温度下降比例θ,令k=1,p=0;
初始温度T0的设定,提出一种自适应的处理方法:随机产生100个满足约束的解,并判断有无可行解;如果有可行解,则设置初始温度为5000,从可行解中选出函数值最小的作为模拟退火算法的初始解X0;否则,初始温度设置为50000,从不可行解中选出罚函数值最小的作为初始解X0;令当前解X1=X0,计算当前解的目标函数值f(X1);
引入一种自适应的函数变换方法:在初始化过程中,记录并计算函数的平均值和最小值,确定一个函数的自适应因子K:有可行解时,否则,
这样就可以通过自适应因子K来调节解的接受准则,从而增强了算法的健壮性。
8.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S44中所述根据Metropolis准则判断是否接受新解,具体包括:
若Δf<0或Δf≥0且满足exp(-Δf/T)>ξ,ξ为区间(0,1)上产生的均匀分布的随机数,则接受X2作为新的当前解,即X1=X2;否则保留当前解X1。
9.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S45中所述执行记忆功能,具体包括:
若迭代次数为1,则将当前产生的最好解记录在记忆矩阵M之中;否则将当前产生的最好解与记忆矩阵M中的解进行比较,如果当前解优于记忆矩阵中的解,则将当前解写入记忆矩阵M中,否则保留记忆矩阵中的解。
10.根据权利要求1所述的基于改进模拟退火算法的城际列车开行方案优化方法,其特征在于,步骤S48中所述算法终止检验,具体包括:
若T<Te,则输出当前解X1,结束程序;否则令k=1,转D2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810687399.3A CN108985500A (zh) | 2018-06-28 | 2018-06-28 | 基于改进模拟退火算法的城际列车开行方案优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810687399.3A CN108985500A (zh) | 2018-06-28 | 2018-06-28 | 基于改进模拟退火算法的城际列车开行方案优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108985500A true CN108985500A (zh) | 2018-12-11 |
Family
ID=64539291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810687399.3A Pending CN108985500A (zh) | 2018-06-28 | 2018-06-28 | 基于改进模拟退火算法的城际列车开行方案优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108985500A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112308259A (zh) * | 2020-10-29 | 2021-02-02 | 合肥工业大学 | 一种基于俄罗斯方块下落的列车分段式座位售票方法 |
CN112749776A (zh) * | 2021-01-12 | 2021-05-04 | 南京信息工程大学 | 一种基于改进混合遗传算法的作业车间调度方法 |
CN112819316A (zh) * | 2021-01-29 | 2021-05-18 | 西南交通大学 | 一种综合客运枢纽轨道交通系统的枢纽运能识别方法 |
CN113788045A (zh) * | 2021-11-16 | 2021-12-14 | 中国铁道科学研究院集团有限公司通信信号研究所 | 基于动态多目标优化控制的有轨电车信号控制系统与方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678425A (zh) * | 2016-01-29 | 2016-06-15 | 中南大学 | 基于多节拍组合的城际铁路列车开行的优化方法 |
CN105857350A (zh) * | 2016-03-17 | 2016-08-17 | 中南大学 | 一种基于区间断面客流的高铁列车开行方法 |
CN107516147A (zh) * | 2017-08-10 | 2017-12-26 | 中南大学 | 一种高速铁路本线列车开行方案优化方法及其系统 |
-
2018
- 2018-06-28 CN CN201810687399.3A patent/CN108985500A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678425A (zh) * | 2016-01-29 | 2016-06-15 | 中南大学 | 基于多节拍组合的城际铁路列车开行的优化方法 |
CN105857350A (zh) * | 2016-03-17 | 2016-08-17 | 中南大学 | 一种基于区间断面客流的高铁列车开行方法 |
CN107516147A (zh) * | 2017-08-10 | 2017-12-26 | 中南大学 | 一种高速铁路本线列车开行方案优化方法及其系统 |
Non-Patent Citations (4)
Title |
---|
宁德圣等: "基于模拟退火算法的改进型退火策略研究", 《东华理工大学学报(自然科学版)》 * |
蒲松等: "基于改进退火算法的高速列车开行方案研究", 《计算机仿真》 * |
陈鹏: "城际轨道交通与城市交通换乘衔接研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
顾元宪等: "桁架结构截面优化设计的改进模拟退火算法", 《计算力学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112308259A (zh) * | 2020-10-29 | 2021-02-02 | 合肥工业大学 | 一种基于俄罗斯方块下落的列车分段式座位售票方法 |
CN112308259B (zh) * | 2020-10-29 | 2022-09-13 | 合肥工业大学 | 一种基于俄罗斯方块下落的列车分段式座位售票方法 |
CN112749776A (zh) * | 2021-01-12 | 2021-05-04 | 南京信息工程大学 | 一种基于改进混合遗传算法的作业车间调度方法 |
CN112749776B (zh) * | 2021-01-12 | 2023-08-15 | 南京信息工程大学 | 一种基于改进混合遗传算法的作业车间调度方法 |
CN112819316A (zh) * | 2021-01-29 | 2021-05-18 | 西南交通大学 | 一种综合客运枢纽轨道交通系统的枢纽运能识别方法 |
CN112819316B (zh) * | 2021-01-29 | 2022-06-10 | 西南交通大学 | 一种综合客运枢纽轨道交通系统的枢纽运能识别方法 |
CN113788045A (zh) * | 2021-11-16 | 2021-12-14 | 中国铁道科学研究院集团有限公司通信信号研究所 | 基于动态多目标优化控制的有轨电车信号控制系统与方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108985500A (zh) | 基于改进模拟退火算法的城际列车开行方案优化方法 | |
CN105205557B (zh) | 一种城市常规公交线网设计方法 | |
CN110070218B (zh) | 面向多制式区域轨道交通的动态服务网络优化设计方法 | |
CN110580404B (zh) | 一种基于城市多模式交通网络的网络运能确定方法 | |
CN107392360A (zh) | 一种电动公交车充电站的规划方法 | |
CN109934382B (zh) | 一种基于轴幅式布局的多级地铁货运网络系统及建构方法 | |
CN109774750A (zh) | 一种基于虚拟耦合模式的动态调度时空决策方法 | |
Zhan et al. | Energy-efficient high-speed train rescheduling during a major disruption | |
CN113988371B (zh) | 基于客流直达的城市轨道交通跨站停开行方案优化方法 | |
CN109109913B (zh) | 一种用于轨道交通系统能效行车组织的信息处理方法 | |
CN110019569A (zh) | 一种获取城市轨道交通运营状态信息的方法 | |
CN110598971A (zh) | 一种基于蚁群算法的响应式公交服务规划方法 | |
CN112434969B (zh) | 一种区域多式轨道交通运力资源调配方法 | |
CN109522584A (zh) | 一种有轨电车发车时刻表设计方法 | |
Rui et al. | Simulated annealing algorithm for solving a bi-level optimization model on high-speed railway station location | |
Lai et al. | Collaborative optimization model for the design and operation of feeder bus routes based on urban metro | |
Fayman et al. | The regional impact of the Channel Tunnel: qualitative and quantitative analysis | |
Ye et al. | A mixed decision strategy for freight and passenger transportation in metro systems | |
Jiang et al. | Overview of intelligent railway transportation systems in China | |
Wang et al. | Energy-efficient operation of multi-trains for metro systems with considering the regenerative energy | |
Qin et al. | Express/local train plan optimization for urban rail transit in condition of full-length and short-turn modes | |
Hongzhou et al. | Research on customized passenger route optimization based on improved genetic algorithms | |
Zhu et al. | A Bilevel Programming Model for Designing a Collaborative Network for Regional Railway Transportation and Logistics: The Case of the Beijing‐Tianjin‐Hebei Region in China | |
Burckhart et al. | Urban intermodality: potentials for connecting the cities’ public transport | |
Yang et al. | A Coordination Optimization for Train Operation and Energy Infrastructure Control in a Metro System |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181211 |