一种城市轨道运力优化方法、电子设备及其存储介质

摘要

一种城市轨道运力优化方法、电子设备及其存储介质，属于城市轨道交通客流管理技术领域。为提高城市轨道运力与乘客出行需求的匹配程度。本发明将车站作为节点构建城市轨道线路网络拓扑结构；统计各个时间片内的各个断面的客流数据，对各个断面乘客出行需求预测数据进行推演得到各方向起始站在各个时间片的发车需求；合并各个时间片不同断面的发车需求；对每个时间片不同断面的最大发车需求进行Fisher聚类得到动态规划转移函数，进行运营时段划分；枚举建立交路方案集合，建立最小化运能缺口模型及约束条件，建立最小化剩余运能模型及约束条件；对每个运营时段进行城市轨道运力优化得到优化后的全天的运营交路计划及各个时段的发车频率。

权利要求

1.一种城市轨道运力优化方法，其特征在于：包括如下步骤：S1、以车站作为节点构建城市轨道线路网络拓扑结构；

S2、统计各个时间片内的各个断面的客流数据，根据历史OD客流，对各个断面乘客出行需求预测数据进行推演，得到各方向起始站在各个时间片的不同断面的发车需求；

S3、合并各个时间片不同断面的发车需求；

S4、对每个时间片不同断面的最大发车需求进行Fisher聚类，得到动态规划转移函数，求解动态规划问题，最终得到城市轨道交通运营时段划分；

步骤S4的具体实现方法为：提取第d个方向、第t个时间片、第j个断面的最大发车需求Ed为d方向的断面集合，pdtj为第d个方向、第t个时间片、第j个断面的发车需求，得到乘客发车需求矩阵：

P为乘客发车需求矩阵，pt为第t个时间片的上下行方向最大发车需求，即pt＝max{p0，t，p1，t}；

对于覆盖时间片[α，β]的某时段，定义时段内发车需求差异情况指标：其中：Δ(α，β)为发车需求差异情况指标，覆盖时间片[α，β]内的最大发车需求的平均值为记b(t，k)为将t个时间片分为k个有序类别，定义损失函数L[b(t，k)]，则m个时间片分为k类的最佳分类结果 需要满足以下动态规划转移函数：其中，Δ(1，τ‑1)为时间片(1，τ‑1)的发车需求差异情况指标，Δ(τ，m)为时间片(τ，m)的发车需求差异情况指标；

使用以上动态规划递推式，能够求出使得损失函数最小的k个时间片划分，最终得到城市轨道交通的运营时段划分；

S5、城市轨道运力优化：枚举建立交路方案集合，建立最小化运能缺口模型及约束条件，建立最小化剩余运能模型及约束条件；对每个运营时段进行所述城市轨道运力优化，得到优化后的全天的运营交路计划及各个时段的发车频率。

2.根据权利要求1所述的一种城市轨道运力优化方法，其特征在于：步骤S1中城市轨道线路网络拓扑结构的方向集合为D＝{0，1}，其中0为上行，1为下行；站点集合为S＝{1，

2，…，n}，其中n为城市轨道线路站点数量，上行站台集合记为S0＝{1，2，…，n}，下行站台集合为S1＝{n+1，n+2，…，2n}；

上行断面集合为E0＝{<1，2>，<2，3>，…，}，下行断面集合为E1＝{，，…，<2n‑1，2n>}， 为d方向、第j个断面的运行时长。

3.根据权利要求2所述的一种城市轨道运力优化方法，其特征在于：步骤S2的具体实现方法包括如下步骤：S2.1、首先基于日乘客出行需求预测数据进行推演，vdtj表示第d个方向、第t个时间片、第j个断面出行需求预测数据，然后按照15分钟的时间片间隔对乘客出行需求进行统计，将一日划分为96个时间片，则时间片集合T＝[1，96]，时间片t∈T；

S2.2、设置mt为时间片t开始的时间戳，根据各个断面的运行时间，将断面的乘客出行需求推演至起点，第d个方向、第j个断面出行需求对应起点站的发车时间片t′为：为d方向、第i个断面的运行时长。

4.根据权利要求3所述的一种城市轨道运力优化方法，其特征在于：步骤S3的具体实现方法为：设置pdti为第d个方向、第t个时间片、第i个断面经过推演的发车需求，首先，初始化pdtj＝0，即 令pdtj＝0；

合并发车需求，即 更新pdt′j＝pdt′j+vdti；pdt′j为第d个方向、第t′个时间片、第j个断面经过推演的发车需求。

5.根据权利要求4所述的一种城市轨道运力优化方法，其特征在于：步骤S5的枚举建立交路方案集合的具体实现方法为：对于一条城市轨道线路中的所有具备折返功能的站点，任意选择两个，进行小交路方案的枚举，经过枚举之后，得到上行方向、下行方向经过各个断面的情况矩阵的总数为：其中：R为折返站数量， 为从R个折返站选出2个折返站的排列数。

6.根据权利要求5所述的一种城市轨道运力优化方法，其特征在于：步骤S5建立最小化运能缺口模型及约束条件的具体实现方法为：最小化运能缺口模型的目标函数为：

其中，min为最小化目标函数，wdtj表示第d个方向、第t个时间片、第j个断面的运能缺口变量；

最小化运能缺口模型的包括以下约束条件：发车频率范围约束、预期平均满载率约束、运能缺口约束；

所述发车频率范围约束为：

f为发车频率下限，f为发车频率， 为发车频率上限；

所述预期平均满载率约束为：

α为预期平均满载率下限， 为预期平均满载率上限，V为核载人数；

所述运能缺口约束为：

7.根据权利要求6所述的一种城市轨道运力优化方法，其特征在于：步骤S5建立最小化剩余运能模型及约束条件的具体实现方法为：最小化剩余运能模型的目标函数为：

其中：min为最小化目标函数，udtj表示第d个方向、第t个时间片、第j个断面的运能缺口变量；

最小化剩余运能模型的包括以下约束条件：预期平均满载率约束、交路方案集合约束、剩余运力约束、线路通过能力约束；

xj为小交路是否经过第j个断面，经过为1、不经过为0；

交路方案集合约束为：

其中：G为小交路集合， 为第g种小交路方案是否经过方向d的第j个断面，经过为1，不经过为0；zg为是否使用第g种小交路方案，使用为1，不使用为0；

所述剩余运力约束为：

所述线路通过能力约束为：

为方案g起点站的折返时间， 为方案g终点站的折返时间。

8.电子设备，其特征在于，包括存储器和处理器，存储器存储有计算机程序，所述的处理器执行所述计算机程序时实现权利要求1‑7任一项所述的一种城市轨道运力优化方法的步骤。

9.计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1‑7任一项所述的一种城市轨道运力优化方法。

说明书

一种城市轨道运力优化方法、电子设备及其存储介质

技术领域

[0001] 本发明属于城市轨道交通客流管理技术领域，具体涉及一种城市轨道运力优化方法、电子设备及其存储介质。

背景技术

[0002] 城市轨道交通是一种高效、绿色的客运方式，近年来在我国得到快速发展。在城市轨道交通实际运营过程中，列车运力能否满足客流需求是一个城市轨道交通系统的重要评价指标。由于乘客出行需求在时空分布上表现不均匀，不合理的列车开行方案容易造成运力不足或者运力浪费，从而导致乘客等待时间过长或者城市轨道交通客运公司盈利过少。因此，如何对运力进行优化，已经成为城市轨道交通实际运营过程中迫不及待的问题。

[0003] 黎茂盛在《城市轨道交通运营管理》中进行了运力优化的研究，此研究方案能够计算出满足最大乘客出行最大断面的需求，但是忽略了实际运营过程中的约束，如最大发车间隔约束、列车通过能力约束、最小追踪时间约束等。在高峰特大客流期间，按照这种方式计算出来的发车间隔常常是不可行的，并且该方案还忽略了乘客出行需求在空间的不均匀性，仅用最大断面的需求进行计算造成大量的城市轨道交通运能的浪费。

[0004] 公开号为CN 111882156 A、发明名称为面向随机动态客流和节能运行的列车时刻表鲁棒优化方法的研究方案为：S1、在列车时刻表鲁棒优化模型中引入有效服务决策变量、列车在首站的发车时间决策变量以及各站间运行曲线选择决策变量；S2、构建列车时刻表与列车的发车时间决策变量和各站间运行曲线选择决策变量之间的线性关联约束；S3、构建列车时刻表与有效服务决策变量之间的非线性关联约束；构建列车运行安全时间间隔约束和列车承载能力约束；S4、结合步骤S1‑S3，构建面向随机动态客流和节能运行的列车时刻表鲁棒优化模型；S5、对鲁棒优化模型进行求解，得到最优列车时刻表。本发明可以高效求解不确定客流条件下的地铁列车时刻表优化问题，提高了解的鲁棒性。但该技术方案所建立的数学模型是一个大规模的MIP问题，存在求解时间长、效率低等问题，使用CPLEX等商业求解器需要花费巨额资金进行购买；该技术方案基于随机生成的多场景进行鲁棒性优化，难以刻画乘客到达的规律，如果不能对乘客的到达特征进行刻画，得到结果不能保证乘客满意度的提升；该技术方案在模型中没有考虑小交路方案的存在，那么得到的结果将会造成运能的浪费，不利于城市轨道运营公司降低成本。

发明内容

[0005] 本发明要解决的问题是提高城市轨道运力与乘客出行需求的匹配程度，发明了一种城市轨道运力优化方法、电子设备及其存储介质。

[0006] 为实现上述目的，本发明通过以下技术方案实现：

[0007] 一种城市轨道运力优化方法，包括如下步骤：

[0008] S1、以车站作为节点构建城市轨道线路网络拓扑结构；

[0009] S2、统计各个时间片内的各个断面的客流数据，根据历史OD客流，对各个断面乘客出行需求预测数据进行推演，得到各方向起始站在各个时间片的不同断面的发车需求；

[0010] S3、合并各个时间片不同断面的发车需求；

[0011] S4、对每个时间片不同断面的最大发车需求进行Fisher聚类，得到动态规划转移函数，求解动态规划问题，最终得到城市轨道交通运营时段划分；

[0012] S5、城市轨道运力优化：枚举建立交路方案集合，建立最小化运能缺口模型及约束条件，建立最小化剩余运能模型及约束条件；对每个运营时段进行所述城市轨道运力优化，得到优化后的全天的运营交路计划及各个时段的发车频率。

[0013] 进一步的，步骤S1中城市轨道线路网络拓扑结构的方向集合为 ，其中0为上行，1为下行；站点集合为 ，其中为城市轨道线路站点数量，上行站台集合记为 ，下行站台集合为 ；

[0014] 上行断面集合为 ，下行断面集合为， 为 方向、第 个断面的运
行时长。

[0015] 进一步的，步骤S2的具体实现方法包括如下步骤：

[0016] S2.1、首先基于日乘客出行需求预测数据进行推演， 表示第 个方向、第 个时间片、第个断面出行需求预测数据，然后按照15分钟的时间片间隔对乘客出行需求进行统计，将一日划分为96个时间片，则时间片集合 ，时间片 ；

[0017] S2.2、设置mt为时间片t开始的时间戳，根据各个断面的运行时间，将断面的乘客出行需求推演至起点，第 个方向、第 个断面出行需求对应起点站的发车时间片 为：

[0018]

[0019] 为 方向、第 个断面的运行时长。

[0020] 进一步的，步骤S3的具体实现方法为：设置 为第 个方向、第 个时间片、第 个断面经过推演的发车需求，首先，初始化 ，即 ， , ,令。

[0021] 合并发车需求，即 ， , ，更新 ；为第 个方向、第 个时间片、第 个断面经过推演的发车需求。

[0022] 进一步的，步骤S4的具体实现方法为：提取第 个方向、第 个时间片、第 个断面的最大发车需求 ，得到乘客发车需求矩阵：

[0023]

[0024] P为乘客发车需求矩阵，pt为第 个时间片的上下行方向最大发车需求，即；

[0025] 对于覆盖时间片 的某时段，定义时段内发车需求差异情况指标：

[0026]

[0027] 其中： 为发车需求差异情况指标，覆盖时间片 内的最大发车需求的平均值为 ；

[0028] 记 为将 个时间片分为 个有序类别，定义损失函数 ，则  个时间片分为 类的最佳分类结果 需要满足以下动态规划转移函数：

[0029]

[0030] 使用以上动态规划递推式，能够求出使得损失函数最小的 个时间片划分，最终得到城市轨道交通的运营时段划分。

[0031] 进一步的，步骤S5的枚举建立交路方案集合的具体实现方法为：对于一条城市轨道线路中的所有具备折返功能的站点，任意选择两个，进行小交路方案的枚举，经过枚举之后，得到上行方向、下行方向经过各个断面的情况矩阵的总数为：

[0032]

[0033] 其中： 为折返站数量， 为从 个折返站选出2个折返站的排列数。

[0034] 进一步的，步骤S5建立最小化运能缺口模型及约束条件的具体实现方法为：

[0035] 最小化运能缺口模型的目标函数为：

[0036]

[0037] 其中，min 为最小化目标函数， 表示第 个方向、第 个时间片、第 个断面的运能缺口变量；

[0038] 最小化运能缺口模型的包括以下约束条件：发车频率范围约束、预期平均满载率约束、运能缺口约束；

[0039] 所述发车频率范围约束为：

[0040]

[0041] 为发车频率下限，为发车频率， 为发车频率上限；

[0042] 所述预期平均满载率约束为：

[0043]

[0044] 为预期平均满载率下限， 为预期平均满载率上限， 为核载人数；

[0045] 所述运能缺口约束为：

[0046]

[0047]

[0048] 进一步的，步骤S5建立最小化剩余运能模型及约束条件的具体实现方法为：

[0049] 最小化剩余运能模型的目标函数为：

[0050]

[0051] 其中：min 为最小化目标函数， 表示第 个方向、第 个时间片、第 个断面的运能缺口变量；

[0052] 最小化剩余运能模型的包括以下约束条件：预期平均满载率约束、交路方案集合约束、剩余运力约束、线路通过能力约束；

[0053]

[0054] 为小交路是否经过第 个断面，经过为1、不经过为0；

[0055] 交路方案集合约束为：

[0056]

[0057]

[0058] 其中： 为小交路集合， 为第 种小交路方案是否经过方向 的第 个断面，经过为1，不经过为0； 为是否使用第 种小交路方案，使用为1，不使用为0；

[0059] 所述剩余运力约束为：

[0060]

[0061]

[0062] 所述线路通过能力约束为：

[0063]

[0064]

[0065] 为方案 起点站的折返时间， 为方案 终点站的折返时间。

[0066] 电子设备，包括存储器和处理器，存储器存储有计算机程序，所述的处理器执行所述计算机程序时实现所述的一种城市轨道运力优化方法的步骤。

[0067] 进一步的，计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现所述的一种城市轨道运力优化方法。

[0068] 本发明的有益效果：

[0069] 本发明所述的一种城市轨道运力优化方法，根据城市轨道历史特征日客流预测数据，使用Fisher聚类算法进行时段划分；根据断面出行需求特征，给出合理的城市轨道运营时段划分，为城市轨道的运力配置奠定基础。

[0070] 本发明所述的一种城市轨道运力优化方法，根据历史客流数据得到预测的出行需求，本专利能够得到与之匹配的城市轨道运力配置，包括交路方案、发车频率等信息。

[0071] 本发明所述的一种城市轨道运力优化方法，针对划分好的时段，根据分时间片的出行OD需求，分阶段对城市轨道运营的发车频率、交路方案分阶段进行优化，从而提升城市轨道运力与乘客出行需求的匹配程度。

[0072] 本发明所述的一种城市轨道运力优化方法，针对数学优化问题，提出一种满足计算效率要求的程序架构，能有效的在多种平台上实施。

附图说明

[0073] 图1为本发明所述的一种城市轨道运力优化方法的流程图；

[0074] 图2为本发明所述的一种城市轨道运力优化方法的城市轨道线路网络拓扑；

[0075] 图3为本发明所述的一种城市轨道运力优化方法的运营时段划分柱状图；

[0076] 图4为本发明所述的一种城市轨道运力优化方法的枚举交路方案示意图；

[0077] 图5为本发明所述的一种城市轨道运力优化方法的运能优化前后对比图。

具体实施方式

[0078] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施方式，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅用以解释本发明，并不用于限定本发明，即所描述的具体实施方式仅仅是本发明一部分实施方式，而不是全部的具体实施方式。通常在此处附图中描述和展示的本发明具体实施方式的组件可以以各种不同的配置来布置和设计，本发明还可以具有其他实施方式。

[0079] 因此，以下对在附图中提供的本发明的具体实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定具体实施方式。基于本发明的具体实施方式，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他具体实施方式，都属于本发明保护的范围。

[0080] 为能进一步了解本发明的发明内容、特点及功效，兹例举以下具体实施方式，并配合附图1‑5详细说明如下：

[0081] 具体实施方式一：

[0082] 一种城市轨道运力优化方法，包括如下步骤：

[0083] S1、以车站作为节点构建城市轨道线路网络拓扑结构；

[0084] 进一步的，步骤S1中城市轨道线路网络拓扑结构的方向集合为  ，其中0为上行，1为下行；站点集合为 ，其中 为城市轨道线路站点数量，上行站台集合记为 ，下行站台集合为 ；

[0085] 上行断面集合为 ，下行断面集合为， 为 方向、第 个断面的运
行时长；

[0086] S2、统计各个时间片内的各个断面的客流数据，根据历史OD客流，对各个断面乘客出行需求预测数据进行推演，得到各方向起始站在各个时间片的不同断面的发车需求；

[0087] 步骤S2的具体实现方法包括如下步骤：

[0088] S2.1、首先基于日乘客出行需求预测数据进行推演， 表示第 个方向、第 个时间片、第 个断面出行需求预测数据，然后按照15分钟的时间片间隔对乘客出行需求进行统计，将一日划分为96个时间片，则时间片集合 ，时间片 ；

[0089] S2.2、设置mt为时间片t开始的时间戳，根据各个断面的运行时间，将断面的乘客出行需求推演至起点，第 个方向、第 个断面出行需求对应起点站的发车时间片为：

[0090]

[0091] 为方向、第 个断面的运行时长；

[0092] 进一步的，记 为时间片 开始的时间戳，根据各个断面的运行时间，可以将断面的乘客出行需求推演至起点，则发车的时间戳应为 对应的时间片 发车，才能满足第 个方向、第 个时间片、第 个断面出行需求。

[0093] S3、合并各个时间片不同断面的发车需求；

[0094] 进一步的，步骤S3的具体实现方法为：设置 为第 个方向、第 个时间片、第个断面经过推演的发车需求，首先，初始化 ，即 ， , ,令。

[0095] 合并发车需求，即 ， , ，更新 ；为第 个方向、第 个时间片、第 个断面经过推演的发车需求；

[0096] S4、对每个时间片不同断面的最大发车需求进行Fisher聚类，得到动态规划转移函数，求解动态规划问题，最终得到城市轨道交通运营时段划分；划分结果如图3所示；

[0097] 进一步的，步骤S4的具体实现方法为：提取第 个方向、第 个时间片、第 个断面的最大发车需求 ，得到乘客发车需求矩阵：

[0098]

[0099] P为乘客发车需求矩阵，pt为第 个时间片的上下行方向最大发车需求，即；

[0100] 对于覆盖时间片 的某时段，定义时段内发车需求差异情况指标：

[0101]

[0102] 其中： 为发车需求差异情况指标，覆盖时间片 内的最大发车需求的平均值为 ；

[0103] 记 为将 个时间片分为 个有序类别，定义损失函数 ，则  个时间片分为 类的最佳分类结果 需要满足以下动态规划转移函数：

[0104]

[0105] 使用以上动态规划递推式，能够求出使得损失函数最小的 个时间段的开始时间片、结束时间片，从而最终得到城市轨道交通的 个运营时段划分结果。

[0106] S5、城市轨道运力优化：枚举建立交路方案集合，建立最小化运能缺口模型及约束条件，建立最小化剩余运能模型及约束条件；对每个运营时段进行所述城市轨道运力优化，得到优化后的全天的运营交路计划及各个时段的发车频率；

[0107] 进一步的，步骤S5的枚举建立交路方案集合的具体实现方法为：对于一条城市轨道线路中的所有具备折返功能的站点，任意选择两个，进行小交路方案的枚举，经过枚举之后，得到上行方向、下行方向经过各个断面的情况矩阵的总数为：

[0108]

[0109] 其中：为折返站数量， 为从 个折返站选出2个折返站的排列数；

[0110] 枚举建立交路方案集合如附图4所示，实心圆表示可折返站点，5个站点中，有4个折返站，即 ，可知小交路方案有5种，即<1,2>、<1,4>、<2,4>、<2,5>、<4,5>，根据 定义可以得到上行方向矩阵 、下行方向矩阵 如下：

[0111]

[0112] 。

[0113] 进一步的，步骤S5建立最小化运能缺口模型及约束条件的具体实现方法为：

[0114] 最小化运能缺口模型的目标函数为：

[0115]

[0116] 其中，min 为最小化目标函数， 表示第 个方向、第 个时间片、第 个断面的运能缺口变量；

[0117] 最小化运能缺口模型的包括以下约束条件：发车频率范围约束、预期平均满载率约束、运能缺口约束；

[0118] 所述发车频率范围约束为：

[0119]

[0120] 为发车频率下限，为发车频率， 为发车频率上限；

[0121] 所述预期平均满载率约束为：

[0122]

[0123] 为预期平均满载率下限， 为预期平均满载率上限， 为核载人数；

[0124] 所述运能缺口约束为：

[0125]

[0126] ；

[0127] 最小化运能缺口模型是线性规划模型，可以使用SCIP、CPLEX、GUROBI进行求解，得到使各个断面运能缺口总和最小的发车频率。

[0128] 进一步的，步骤S5建立最小化剩余运能模型及约束条件的具体实现方法为：

[0129] 最小化剩余运能模型的目标函数为：

[0130]

[0131] 其中：min 为最小化目标函数， 表示第 个方向、第 个时间片、第 个断面的运能缺口变量；

[0132] 最小化剩余运能模型的包括以下约束条件：预期平均满载率约束、交路方案集合约束、剩余运力约束、线路通过能力约束；

[0133]

[0134] 为小交路是否经过第 个断面，经过为1、不经过为0；

[0135] 交路方案集合约束为：

[0136]

[0137]

[0138] 其中： 为小交路集合， 为第 种小交路方案是否经过方向 的第 个断面，经过为1，不经过为0； 为是否使用第 种小交路方案，使用为1，不使用为0；

[0139] 所述剩余运力约束为：

[0140]

[0141]

[0142] 所述线路通过能力约束为：

[0143]

[0144]

[0145] 为方案 起点站的折返时间， 为方案 终点站的折返时间；

[0146] 最小化剩余运能模型是线性规划模型，模型通过数学规划求解器，如SCIP、CPLEX、GUROBI等求解出使目标函数最小的交路方案、发车频率。

[0147] 对全天划分后的每一个运营时段通过进行运能优化，得到全天的运营交路计划及各个时段的发车频率。

[0148] 图5为本发明所述的一种城市轨道运力优化方法的优化后的柱状图，从图5可知：（1）优化前，上行断面从塘坑至草埔区间十个断面出现运能不足的情况，而优化后， 运能不足的断面仅剩两个，表明本实施方式能够有效地改善运能不足的情况。（2）原来的方案中没有交路，经过优化之后，本实施方式产生了一段小交路交路（塘坑‑益田）。表明本实施方式的方法能够从小交路方案集合中选择一种最优的小交路作为城市轨道交通交路方案，从而减少运能的浪费。

[0149] 具体实施方式二：

[0150] 电子设备，包括存储器和处理器，存储器存储有计算机程序，所述的处理器执行所述计算机程序时实现具体实施方式一所述的一种城市轨道运力优化方法的步骤。

[0151] 本发明的计算机装置可以是包括有处理器以及存储器等装置，例如包含中央处理器的单片机等。并且，处理器用于执行存储器中存储的计算机程序时实现上述的基于CREO软件的可修改由关系驱动的推荐数据的推荐方法的步骤。

[0152] 所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器  (Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列 (Field‑Programmable Gate Array，FPGA) 或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

[0153] 所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序（比如声音播放功能、图像播放功能等）等；存储数据区可存储根据手机的使用所创建的数据（比如音频数据、电话本等）等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡（Smart Media Card, SMC），安全数字（Secure Digital, SD）卡，闪存卡（Flash Card）、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。

[0154] 具体实施方式三：

[0155] 计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现具体实施方式一所述的一种城市轨道运力优化方法。

[0156] 本发明的计算机可读存储介质可以是被计算机装置的处理器所读取的任何形式的存储介质，包括但不限于非易失性存储器、易失性存储器、铁电存储器等，计算机可读存储介质上存储有计算机程序，当计算机装置的处理器读取并执行存储器中所存储的计算机程序时，可以实现上述的基于CREO软件的可修改由关系驱动的建模数据的建模方法的步骤。 所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器（ROM，Read‑Only Memory）、随机存取存储器（RAM，Random Access Memory）、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。

[0157] 本发明的技术关键点和欲保护点为：

[0158] (1)提出的城市轨道运力优化技术路线；

[0159] (2)所提出的数学优化模型；

[0160] (3)提出的数学模型可由其他精确数值解法或启发式求解，无论求解算法如何，应该属于专利保护范围。

[0161] 需要说明的是，术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

[0162] 虽然在上文中已经参考具体实施方式对本申请进行了描述，然而在不脱离本申请的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，本申请所披露的具体实施方式中的各项特征均可通过任意方式相互结合 起来使用，在本说明书中未对这些组合的情况进行穷举性的描述仅仅是出于省略篇幅和节 约资源的考虑。因此，本申请并不局限于文中公开的特定具体实施方式，而是包括落入权利要求的范围内的所有技术方案。