tceic.com
学霸学习网 这下你爽了
相关文章
当前位置:首页 >> 英语学习 >>

外文翻译


西南交通大学本科毕业论文(外文翻译)

利用混合整数线性规划求解固定与移动闭塞条 件下的多列车操纵优化问题
王义惠
摘要: 本文主要解决基于固定闭塞与移动闭塞条件下的多列车操纵优化问题,利用婪式 法以及同步法进行求解。对于两种求解方法,文中都将列车操纵优化问题转化为混合 整数规划问题。特别地,文中将能量消耗目标函数以及列车模型利用分段仿射模型进 行近似处理。在模型建立过程中,考虑了线路基本阻力、分段速度限制和最大牵引力 等影响因素。并且,文中将固定闭塞与移动闭塞条件下的前车运行进行离散化,利用 分段仿射处理转化为混合整数规划问题。利用同步法与婪式法进行求解,发现同步法 的仿真结果的控制更加优化,但是需要花费更多的计算时间。并且将两种方法的仿真 结果与伪谱法的仿真结果进行比较,发现伪谱法在求解优化控制中表现更佳。主要表 现为婪式混合整数规划法的求解结果比婪式伪谱法的求解结果在节能与准时方面表 现更差,但是求解时间呈现一至两个数量级的差异,婪式伪谱法需要花费更多的时间 进行求解。同样的,同步混合整数规划法与同步伪谱法存在相同的差异性。 关键词:运行曲线,列车操纵,信号系统,混合整数规划,伪谱法 1.简介 随着能源价格的增长与环境问题的严重化, 交通运输系统的节能变得越来越重要。 轨道交通因其安全,快捷,准时与舒适的服务,在交通运输系统可持续发展中扮演着 至关重要的作用。列车运行能耗占到中国铁路年运营与维护能耗的 13%到 16%,因此 将其作为主要优化目标,哪怕是小小的一个提升,也会为铁路系统节约一大笔资金。 为了达到这个目的,一些辅助驾驶系统已经被开发,例如 FreightMiser、Metromiser 以及 driving style manager。随着铁路系统的现代化,自动驾驶系统在确保准确停车, 准点,节能和驾驶舒适方面起到越来越重要的作用。铁路控制中心或者自动化控制系 统主要依据列车监控系统提供的线路基本阻力、速度限制、最大牵引力和制动力等数 据来决定列车运行曲线。 在参考文献中, 对于列车操纵优化的研究开始于 1960 年代,主要致力于解决列车 从一个站到相邻站的运行曲线。自从发现该研究对于列车运行节能有帮助,多种方法 第1页

西南交通大学本科毕业论文(外文翻译) 被应用到这个研究中,主要分成两种类别:解析法与数值优化。对于解析法,Howlett 等人将极大值原理运用到列车运行操纵优化中,并将操纵工况划分为四种情况:全力 牵引,恒速运行,惰行以及全力制动。但是在该方法中,如果将更多现场状况的约束 引入到模型中, 形成非线性约束条件, 很难得到解。 随着当今计算机计算能力的提升, 越来越多数值优化方法被应用到这个研究中,比如婪式法以及遗传算法等,但是通过 这些方法不一定能得到优化方案。另一方面,二次多项式优化被应用到求解该问题, 主要特点是通过分段仿射处理列车模型。 该理论主要由 Va?ak 提出, 在先前的文章中, 我们希望将运行曲线优化问题转化为混合整数规划问题,并通过现成的求解工具进行 全局求解。 但是先前提到的方法都没有考虑铁路信号系统的影响,例如固定闭塞信号系统与 移动闭塞信号系统。固定闭塞信号系统是通过轨旁信号机对固定闭塞区间进行保护, 列车只有在信号机显示通过状态下,才能进入对应的闭塞区间。移动闭塞信号系统需 要通过列车与区间控制者之间的沟通,得知列车位置以及速度的实时数据,对于每一 趟列车定义移动闭塞。与固定闭塞相比,移动闭塞情况下,运行列车运行间隔减小, 因此增加了区间通过能力。 Lu 和 Feng 考虑在移动闭塞区间下, 同一方向, 同一轨道, 前车对后车造成的约束下,对后车运行曲线进行优化,对前后车使用相同的遗传算法 进行运行曲线优化, 使总的能源消耗降低。 Gu,Lu 和 Tang 利用非线性优化对后车运行 曲线进行优化,在前车运行或者停车的两种情况都被考虑,并且提出对应的后车运行 策略。Ding 等人将移动闭塞区间下多车运行的约束条件考虑到列车运行节能模型中, 并且将全力牵引,惰行,全力制动三种工况结合到算法中,通过既定规则安排工况顺 序,最终得到优化的运行曲线。 在本文中,移动闭塞区间以及固定闭塞区间条件下,由于前车造成的对于后车的 约束进行建模,先将其离散化,再通过分段仿射处理转化为线性约束。因此,这些约 束可以融入到混合整数规划模型中,并且通过比其他方法更快地进行求解。文中,运 用婪式法以及同步法对问题进行求解,并且与伪谱法进行对比。在过去地十年中,伪 谱法在数值优化领域被大量使用,例如轨道转移,轨迹引导,引力控制,因此文中选 择其作为仿真比较。 接下来的论文内容如下:第二部分,列车模型以及混合整数规划在单车运行控制 的应用基于先前论文进行说明;第三部分,介绍固定闭塞区间以及移动闭塞区间的原 则;第四部分,将固定闭塞条件下,后车跟驰前车受到的约束转换到混合整数规划模 型中;第五部分,介绍了移动闭塞条件下的混合整数规划模型;第六部分,利用北京 亦庄的数据进行仿真计算;第七部分,进行总结与展望。 2.列车运行模型与混合整数规划模型 第2页

西南交通大学本科毕业论文(外文翻译) 本部分对于单车运行优化模型建立与混合整数规划模型的运用进行总结。 2.1 优化控制模型 连续质点模型经常被运用到列车运行优化控制,如下: mρ = ? ? , ,


=

1 2

,(1)

其中,m代表的列车单位质量,ρ代表动态质量参数, 代表单位质量的动能,等价于 0. 5 2 , 代表列车运行速度, 代表列车的位置, 代表控制变量即牵引力或者制动力, 上界为最大牵引力 ,下界为最大制动力 ,因此 < < , 代表 基本运行阻力,包括滚动摩擦阻力以及空气阻力, , 代表附加阻力,由坡度,曲 线,隧道等因素产生。 单位质量的动能 = 0.5 2 和时间t是相对固定的,位置参数在模型中是固定参数。 运行曲线优化问题可以建立如下: J= s. t. ≤ ≤ , 0 < ≤ , = , = , t = 0,t = ,
2 其中, 目标函数J是不考虑再生制动的列车运行能量消耗, 等价于0.5 ,

max? 0,

(2)

代表最大运行速度,主要取决于列车参数以及线路状况,是关于位置的分段 常数函数, , ,t 分别代表行程的起始位置,起始单位动能,起 始时刻, , , t 分别代表行程的结束位置, 结束单位动能, 结束时刻, 代表计划运行时间又时刻表或者重订时刻表给定。文中假设 > 0,即列车速度 大于零,在路途中不停车,这个假设在现实中可以通过定义起车与停车为极小值来实 现,而且列车运行管理系统也不希望列车在路途中途停车。 2.2 转化规则 本小节将介绍通过 Williams 提出的转化规则将非线性优化控制问题转化混合整数 规划问题。首先将非线性的列车模型转化为混合逻辑模型,然后将优化控制问题的约 束条件利用分段仿射处理转化为混合整数线性约束,最终形成混合整数规划问题。 考虑声明 ≤ 0,: → 表示仿射, ∈ ,X ∈ ,并且使 = max ∈ , = min ∈ 第3页

西南交通大学本科毕业论文(外文翻译) 如果我们引入逻辑变量δ ∈ 0,1 ,那么可以得到如下等式: 当且仅当 ≤ 1 ? δ 成立时, ≤ 0 ? = 1 成立 ≥ + ?

其中, 是一个小正整数(典型的机械精度) ,引入其将严格等式转化为非严格等式, 适合混合整数规划问题的框架。 两个逻辑变量的乘积1 2 ,可用另一个辅助逻辑变量3 代替,也就是说 3 = 1 ? 1 = 1 ∧ 2 = 1 等价于 ?1 + 3 ≤ 0, ?2 + 3 ≤ 0, 1 + 2 ? 3 ≤ 1. 乘 积 δ 可 用 辅 助 变 量 z = δ 代 替 , 并 且 满 足 = 0 ? = 0 以 及 = 1 ? = ,因此z = δ 等价于 ≤ , ≥ , ≤ ? 1 ? , ≥ ? 1 ? . 不等式集(5)-(7)均基于为仿射。 2.3 混合整数规划方法 在 Wang 的论文中将列车运行的空间连续模型离散化,将起点站与终点站之间的 位置区段 , 不包含中间站离散化成N 个间隔,并且假设线路,列车的特性 以及牵引力或者制动力为恒量在每一个区段 , +1 ,区段长度? = +1 ? , k = 1,2, … . , N。然后将空间离散模型通过 PWA 公式转化为分段仿射模型,并且利用 2.2 中描述的转化规则将分段仿射模型建立成为混合逻辑动态模型,如下: x + 1 = + + 1, + 2, + 1 +1, + 2, + 1 + ,(8) 1, + 2, + 1 + 3, + 4, + 1 ≤ 5, + 6, + 7, (9) 其中,x = ,


,δ · 和z · 包含二元变量以及辅助变量,是由于引入转化

规则。同时,等式(9)存在上下界约束变量E ,t 以及u 。为了方便起见,我 们用E 和t 分别标记 和 , 混合逻辑动态模型的系数矩阵主要由列车模型, PWA 近似值以及上下界等因素决定。 列车操作优化模型的目标为不考虑再生制动情况下的列车运行能耗,可以描述如 下:

第4页

西南交通大学本科毕业论文(外文翻译)


J=
=1

0,

?

在 Wangl 的论文中, 引入新的辅助变量w 将优化控制模型重新如下建立, 辅助变量 w 的引入主要是为了转化等式(10)中的 0, 。 min


约束于: 1 ≤ 2 1 + 3 4 = 5 1 + 6 其中, = 0 … 0 ?1 … ? , = 1 1 2 2 = , = , = ? ? + 1


, 1 1 2 2 , = ? ? + 1 + 1

进行合适的定义矩阵和矢量1 ,2 ,3 ,4 ,5 和6 ,混合整数规划问题可以通过现 有的商业和免费软件进行求解,例如 CPLEX,Xpress-MP,GLPK。 3.区间闭塞信号系统的原则 区间闭塞信号主要是用来保持同一线路上前后相邻列车的安全距离,固定闭塞信 号系统以及移动闭塞信号系统将在接下来阐述。 3.1 固定闭塞信号系统 固定闭塞信号系统(FBS)在当今的列车控制系统中被普遍使用。在固定闭塞信 号系统中,线路被分成不同的区段,主要依据最高列车运行速度,最差列车制动率以 及信号显示的制式,信号显示的制式决定了信号可见的颜色,例如绿灯,黄灯,红灯 等。每一个区段只能被一辆列车占用,该规则通过轨道电路进行检验。此外,区段通 过轨边信号或者车载信号进行保护,轨边信号在铁路上的经典信号系统,车载信号被 使用的越来越多,尤其是在高速铁路中,主要由于在高速下,列车司机无法有效识别 轨边信号。 固定闭塞信号系统主要分成两大类型, 单区段信号系统和多区段信号系统。 在单区段信号系统中,信号的显示根据对应区段的状况,同时需要有远信号灯为接近 列车提供必要的信息。在多区段信号系统中,信号的显示需要根据对应的两个或者多 个信号决定,一个简单的例子就是双区段信号系统伴随三中显示制式的信号灯,也就 是红,黄和绿被称为三显式信号系统。我们将阐述三显式信号系统在固定闭塞信号系 统运用中的限制,但是该方法论可以被运用到其他制式的固定闭塞信号系统中。

第5页

西南交通大学本科毕业论文(外文翻译)

图一固定闭塞信号系统下的速度编码

三显示信号系统如图一所示,在装备列车自动防护系统情况下,每一个区段在轨 道电路始端携带速度编码,速度编码主要包含两部分,本区段的批准速度编码以及后 区段的指导速度编码。速度编码由轨道电路控制设备进行编译,由轨道进行传输,通 过车载天线传输到列车,如果列车想要进入零速区段或已被占用的区段,以及速度超 过批准速度,列车会进行紧急制动。图一中显示了速度编码如何工作,当某一列车进 入区段 4 并且接近信号机 3,信号机 3 将会收到 / 编码,去确定该区段的运 行 速 度 以 及 下 一 区 段 的 目 标 速 度 。 当 列 车 进 入 区 段 2 , 编 码 变 成 / ,由于下一区段被列车 1 占用,因此当列车接近区段 2 末端时,速度需 要降到 (通常情况下 = 0) 。如果列车试图进入区段 1,车载系统将会检验最 小速度编码( / ) ,并且进行紧急制动。 为了确保列车操纵不受信号系统的影响, 也就是说列车操纵不受前面列车的影响, 最小间隔被引入, 最小间隔也就是车站发车间隔。 在不被干扰的固定闭塞信号系统中, 最小间隔可以被定义为: , = 2 + +
2 2



2 + + + + 2

其中, 代表区段长度, 表示后车在反应时间 内的运行距离,反应时间主要受司 机以及列车参数影响, 表示后车的最高速度, 表示最大制动率, 表示前车的 停留时间, 表示安全间隔的长度, 表示前车的加速度。 表示前车的车体长度,

3.2 移动闭塞信号系统 随着铁路系统列车运行密度的增加,固定闭塞信号系统下的运输能力无法符合显 示需求,即使固定闭塞信号系统下的运输能力可以通过缩短区段长度来增加,但是信 号设备以及线路设备的安装使用费用并不合理。因此,移动闭塞信号系统(MBS)由 于其良好的表现被采用。 移动闭塞信号系统依赖于列车于区域控制室之间连续双向的联系,区域控制室计 算通过计算区域内列车的移动权限限制来确保列车之间保持安全距离,更加详细地阐 第6页

西南交通大学本科毕业论文(外文翻译) 述,移动权限限制包含前车的车尾以及安全距离,也就是后车允许运行的位置。在文 献(Pearson,1973)中阐述了 4 种移动闭塞信号系统:空间移动闭塞信号系统,时间 移动闭塞信号系统,纯粹的移动闭塞信号系统,相关的移动闭塞信号系统。 Takeuchi (2003)等人评估了前三种系统,并且将其与固定闭塞信号系统基于稳定性和扰动性 进行比较,纯粹的移动闭塞信号系统表现最佳。同时,其表示相关的移动闭塞信号系 统没有被铁路系统所使用,尽管其经常在公路系统中被使用。因此,本文将围绕纯粹 的移动闭塞信号系统进行展开,相关的方法论可以在其他移动闭塞信号系统中使用。 同时,纯粹的移动闭塞信号系统是现已使用的移动闭塞信号系统的基础,其原理如下 图二所示。

图二纯粹的移动闭塞信号系统的原则

在纯粹的移动闭塞信号系统中,列车最小追踪间隔是基于后车的瞬间制动距离加 上一个安全距离,利用最小追踪间隔可以避免前车突然停车时,后车对其的冲撞。在 现实中,最小追踪间隔由于司机反应时间或者列车反应时间会变得更大,并且需要考 虑列车长度。因此,实际移动闭塞信号系统下,最小追踪间隔计算公式如下: ? ≥




+ 2

2

+ +

其中, 和 分别代表前后车在时间的位置, 代表后车的速度, 代表最 大制动率, 代表安全距离, 代表在司机或列车反应时间 内的后车运行距离, 代表前车的车长,反应时间 可从历史数据中获得。最小追踪间隔可以通过相邻列车

最小追踪时间间隔近似求得:
? ≥ + + 其中, 和 分别代表前后车经过位置的时刻,后车的制动时间 以及由于 安全距离和前车车长造成的时间间隔 可以如下计算: = / = + /

其中, 代表后车在位置的速度。 第7页

西南交通大学本科毕业论文(外文翻译) 为了使列车操纵不受信号系统的影响,也就是不被前车的运行影响(原则上缓行 或者停车) ,最小追踪时间间隔这个概念被引入。车站发车间隔如下定义:
, = + ? = + + , + 同时, ? = + , + ,其中 , 表示后车从最高速度制动到停止所需 的时间, 也就是说 表示前车完全驶出安全部分 (特殊的部分, , = / , 用来保护前车) ,其中包含安全距离, = 在计算式被视为常量。 2 + + / ,前车的加速度

4.FBS 的约束条件以及转化为 MILP 问题 在固定闭塞信号系统中,由于前车造成的约束条件首先被离散化到每一个网格点 , k = 1,2, … , N + 1, 然后这些逻辑约束被转化线性约束, 这些线性约束可以被 MILP 方法利用两种求解方法被利用到固定闭塞信号系统下多车运行操纵优化问题的求解。 4.1 离散化 FBS 约束条件 我们假设对于每个m ∈ 1,2, … , , 存在一个指数 ∈ 1,2, … , + 1 , 下式成立: , = 并且定义一个分段常数等式,如下: ? = ,对于 < ≤ +1 ,m ∈ 1,2, … , 其中, 位于固定区间 ? , ? +1 ,在三显示固定闭塞信号系统中,如图一所示, 由于前车运行造成的限制条件可以转化到每个网格点 ,如下所诉: 如果前车与后车位于同一个闭塞区间,也就是说 ∈ ? , ? = ②如果前车在后车前一个闭塞区间,也就是说 ∈ ? 是说: ≤ ? = ? +1 在间隔 ? , ? +1 内,假设列车减速度为常量,根据位置,速度以及加速度之间的 关系,我们可以得到: 2? 2? ? +1 ? ? ? ?
2 2 = ? +1 +1

,这种

情况在现实中是不允许发生的,那么后车的速度必须降低到最低速度,也就是说: , ? ,那么

+2

后车在位置? 的速度需小于或等于 ,在位置? +1 的速度需等于 ,也就

2 2 = , ?

第8页

西南交通大学本科毕业论文(外文翻译) 其中,? 表示减加速度, , 表示在位置 的最高速度,通过消除式子(23) 与式子(24)中的? ,得到: , =
2 2 2 + ( ? )

? ? ? +1 ? ?

因此, , 为常量,等式右边均为常量,可得到如下约束: ≤ , 如果前车在后车前两个闭塞区间,也就是说 ∈ ? 也就是说: ≤ ? ≤ ? +1 与求解 , 相似,我们可以得到: , =
2 2 2 + ( ? ) +2

, ?

+3

,那么后

车在位置? 的速度需小于或等于 ,在位置 ? +1 的速度需小于或等于 ,

? ? ? +1 ? ?

其中, , 表示列车在位置 的最高速度, , 同样也是一个常量,因此,我们可 以得到如下约束: ≤ , 4.2 将 FBS 中的约束条件引入 MILP 求法 为了将上诉的逻辑约束转化为线性约束,引入以下二进制变量: () ≤ (? () ≤ (?
+1 ) +2 ) +3 )

? 1 () = 1 , ? 2 () = 1 , ? 3 () = 1

() ≤ (? 以及必要的额外约束:

> 上式意味着后车经过位置 晚于前车经过该位置,以避免列车冲突。同时,基于变量 1 (),2 (),3 ()的定义,可以得到下面的逻辑式: 1 = 1 ? 2 = 3 = 1 , 2 = 1 ? 3 = 1 三显示固定闭塞信号系统环境下,由前车造成的约束可以如下重建: δ1 () () ≤ 1 ? δ1 ()2 () () ≤ , 1 ? δ2 ()3 () () ≤ , 式 (33) 中小于等于情况, 当δ1 = 0, δ1 () ()等 第9页

西南交通大学本科毕业论文(外文翻译) 于 0,并不等于 。 通过定义
= ? (? + )

≥ max? ( ? (?
+ )) ≥ ? (?

+ )), + ),

= min? ( ? (? ? ? ? ?

通过转化特性式子(5)逻辑约束式(30)等价于下式:
+

≤ (1 ? δ ()), ≥ + ( ? )δ ()

+

对于 = 1,2,3,其中, 表示前车到达终点站的时刻, 表示后侧滑离开始发站

的时刻 , 表示精 度。 同时, 我们 定义 二进 制变量 δ4 = δ1 δ2 和 δ5 = δ2 δ3 去处理非线性条件δ1 δ2 和δ2 δ3 。 根据第 2.2 节的转化特性式子, δ4 和δ5 的定义与式子(6)的等价。同时,辅助变量 ()被引入去解决非线性 因素 () (),被定义如下: = 化为混合整数规划的条件。 4.3 基于 FBS 系统下的多车优化控制问题 现在有两种发放可以用于求解基于 FBS 系统下的多车优化控制问题,为了方便起 见,我们将问题考虑为两车的运行曲线优化控制问题,但是方法论可以被应用到多车 问题中。其中一种为婪式算法(时序法) ,前车的运行速度曲线首先被确定,然后后 车的速度曲线优化基于前车的结果。另外一种方法是同时计算前后车的运行速度曲线。 使用婪式算法,前车的运行速度曲线转化为 MILP 问题进行求解,然后后车的控 制优化求解与 Wang(2011)中的求解相同,只是加入 FBS 系统下前车产生的约束条 件。式(8)和式( 9 )中的协同矩阵由前车确定,前车的运行速度曲线确定之后, ?
+1

= 1 ? 5

这个定义与式(7)等价,因此固定闭塞信号系统下由于前车造成的约束条件可以转

变为已知值,因此 和 变为常量,那么式(36)变为线性约束。

当前车与后车的优化同时进行的时候,两车的模型可由形式(8)和(9)确定, 优化控制问题可以被写成形式(11)和(12) ,但是其中包含由于列车以及 FBS 系统 产生的约束条件。但是当前车的运行速度曲线被优化后, (? 约束条件(36)可以被写成:
? ? ? ? + + ≤ ? + , + +1 )也变为已知值,

? ≤ ? ? ?

为了解决非线性问题(38) ,我们定义: = ?
+



= 1 ? 3与 相似,辅助变量 等同于线性约束,依

据第 2.2 节阐述的转化规则。 第 10 页

西南交通大学本科毕业论文(外文翻译)

5.MBS 的约束条件以及转化为 MILP 问题 5.1 离散化 MBS 约束条件 回忆第 2.3 节中,将列车操纵优化的混合逻辑动态模型离散化成 N 个空间,伴随 着网格点 ,k = 1, … , N + 1。这里我们将 MBS 系统的约束在网格点 离散化如下:
≥ + + + ≥ + + , + +

= 1,2, … , = + 1

同时,一些调解约束被引入,为了确保网格点减的位置能够满足由于 MBS 系统造成 的约束。根据式子(15) ,得到对于s ∈ , +1 满足以下约束:
? ? ? () ≥ ()

如果我们假设式子(42)左边为仿射方程,在区间 , +1 内,我们可以得到: 1?α
? ? ?

+ α 1 +

? ? + 1 ? + 1 ≥ ( + ? )

其中,α在有限子集 ∈ 0, 1 ,举例来说 = 0.1,0.2, … ,0.9 ,如果前车是固定的, 那么 ( + ? )为已知量,对于α = 0和α = 1是检索的(除了如果k = N ? 1) 。但是 如果前车的运行速度曲线是未知的,那么将同时优化前后车的运行速度曲线,在这种 情况下,条件 ( + ? )是未知的。如果我们假设式子(42)右边为仿射方程,也 就是: + ? = 1 ? + + 1 , 那么需要检验式子 (42) 在点和 + 1 (也就是 = 0和 = 1) 是否成立, 如果成立, 那么线性式子(42)将在任何中间点成立。 5.2 将 MBS 中的约束条件引入 MILP 求法 声明约束(40) , (41)和(43)在 (), ( + 1)和 ()是线性的,但是因为 安全时间间隔(17)对于列车速度 ()是非线性,所有约束在 ()和 ( + 1)是非 线性的。更多的,单位质量的动能是 ()而非 = 0.5( ())2 (具体参考第 2
章) ,因此对于 (),制动时间 ()和安全时间间隔 ()是非线性方程,如下: =

1 2 1 2

= ( + )

非线性方程f1 ? : → 2 和f2 ? : → 1/ 2 可以通过 PWA 方程进行近似化处理, 存在很多种近似化处理方法。在本文中,我们首先选择采用基于复杂性和精确性的多 层次 PWA 方法,然后用最小二乘法对f1 ? 对应的仿射方程的间隔长度和特征量进行 第 11 页

西南交通大学本科毕业论文(外文翻译) 优化,然后对f2 ? 进行相同的处理。如果我们考虑二段式,那么方程 f1 ? 和f2 ? 可以 被 PWA 近似处理如下: f1, f2, = = 1 + 1 ≤ ≤ 1 2 + 2 1 ≤ ≤ 1 + 1 ≤ ≤ 1 2 + 2 1 ≤ ≤

其中,1 ,2 ,1 ,2 ,1 ,2 ,1 ,2 和1 优化参数。现在约束(40)可以被近似 化如下线性方程:
≥ + +

1 1 + 1 + + 1 + 1 ≤ ≤ 1 1 2 + 2 + + 2 + 2

≥ + +

1 ≤ ≤ 同样的,约束条件(41)和(42)可以同样被转化为线性方程,这 些由于 MBS 系统造成的近似约束可以转入到 MILP 问题。 5.3 基于 MBS 系统下的多车优化控制问题 第 4.3 节中提到的婪式算法和同步算法可以被应用到 MBS 系统中。 运用婪式算法, 当前车的运行速度曲线已知时,变量 ()和 + ? 变成已知量对于后车,后车 的运行速度曲线的求解和 Wang(2011)的文章中阐述的相同,但是需要考虑第 5.2 节中提到的关于 MBS 系统形成的约束条件。当同时优化前后车的运行速度曲线时, 采用形式(8)和形式(9) ,然后写成 MILP 问题的格式(11)和(12) 。与单车问题 相比较,多车问题的参数,二进制变量,辅助变量以及约束条件呈现线性增加,因此 问题的规模变得更大,需要更多的求解时间。因为同步算法为全局优化算法,其结果 将比婪式算法更加优秀。 6 案例仿真 为了体现婪式算法和同步算法在解决基于固定闭塞信号系统和移动闭塞信号系统 的多车运行速度曲线优化的表现,以北京亦庄地铁线路为基础进行仿真。 6.1 建立 我们将 MILP 方法与伪谱法进行比较,同时考虑婪式算法和同步算法,为了简单 起见,我们仅仅考虑北京亦庄铁路线路的两个站:宋家庄站和肖村站,两个站之间的 线路长度伪 1332m,速度限制和线路坡度如下图三所示:

第 12 页

西南交通大学本科毕业论文(外文翻译)

图三宋家庄站与肖村站间速度限制和线路坡度

列车与线路参数如下表一所示:
表一列车与线路参数

在文献(Hansen,2008)中,旋转质量参数取 1.06,本文中也取相同的值。根据第 2.1 节中的假设单位质量动能稍大于正临界值 。在本次仿真中,我们取最小单位质量 的动能为0.1J ,对于列车的速度,最大牵引力为非线性方程,方程的最大值我们取 315kN。列车的最大加速度和最大减加速度,取 0.8m/ 2 ,为了确保低驾驶水平的实 际可以按照速度曲线进行操作。当取列车质量为2.78 × 105 时,最大牵引力和最大 制动力分别为222.4kN和?222.4kN。本文中考虑的目标函数为,不存在再生制动情况 下的列车能量消耗。 在本次仿真中,两趟列车从宋家庄站向肖村站发车,发车间隔为75s。我们考虑两 种情况:FBS 系统和 MBS 系统,假设前车在整个仿真过程中有故障,其最高速度为 40km/h,也就是11.1m/s,并且假设前后车在肖村站进入不同的站台股道,后车将在 该站超车。计算最小追踪间隔的参数如表二所示:
表二计算最小追踪间隔的参数

列车的长度取 90m,司机的反应时间取 1s。对于 FBS 系统,我们假设两站之间包含 第 13 页

西南交通大学本科毕业论文(外文翻译) 四个等长的闭塞区间, 也就是长度为 333m, 根据式子 (13) 计算最小追踪间隔为 98.4s。 对于 MBS 系统,根据表二中的数据,求得 ? = 44.6,最小追踪间隔为 69.6s。 指出初始的最小追踪间隔为 75s,大于 MBS 系统的最小追踪间隔,小于 FBS 系统的 最小追踪间隔。 在 MILP 问题中, 区间 , +1 对应的间隔长度? 的取决于限制速度, 线路坡度, 固定闭塞区间长度等因素, 如果分割的空间数量 N 越大, 那么 MILP 算法的精度越高, 但是需要更长的计算时间。对于这种情况,我们取 N 为 20,40 和 60,那么间隔长度 对应分别为 66.6m, 33.3m 和 22.2m。 运行速度曲线优化问题中的非线性条件如列车模 型中的微分方程,通过 PWA 方程进行近似处理,如果处理 PWA 方程时,取更多的子 功能, 那么近似处理的精度将更高, 同时对两段化和三段化处理进行比。 我们在 Matlab 环境下,用 CPLEX 对问题进行求解。 伪谱法是解决优化控制问题表现最好的算法,利用基于正交配置点,如高斯点, 正交多项式对状态方程和控制方程进行近似处理,微分方程可以用代数方程组进行近 似处理。伪谱法近似处理的误差可以用过增加正交配置点的数量来降低,我们分别取 高斯点数量为 20,40,80 和 120。现存在很多工具箱可以求解伪谱法,其中一个是 PROPT,支持代数微分方程,具有多个求解工具,如 MINOS 和 SNOPT,来求解非线 性问题。本次仿真中,我们在 Matlab 和 SNOPT 环境下,通过 Tomlab 进行求解。 6.2 基于 FBS 系统的仿真结果 表三展示了基于 FBS 系统的 MILP 和伪谱法用婪式算法和同步算法的结果表现对比, 结果表现中包含的例如能量消耗和时间差异,均通过非线性方程进行求解。总能量消 耗包含了前后车的能量消耗,时间差异表示前后车实际运行时间与计划运行时间的差 值。每趟车的能量消耗与实际运行时间和计划运行时间相关,如果实际运行时间比计 划运行时间长,那么列车的平均运行速度比较低,那么总能量消耗就比较少。CPU 计 算时间表示问题求解的总用时,对于婪式算法,是前后车分别计算时间,对于同步算 法,是前后车同步计算完成所需的时间。

第 14 页

西南交通大学本科毕业论文(外文翻译)
表三基于 FBS 系统的 MILP 和伪谱法用婪式算法和同步算法的结果表现对比

更多的,图 4,5 展示了,基于 FBS 系统,两种不同方法求得的优化控制输入量 和速度位置曲线,其中,MILP 方对应的空间数量 N 和伪谱法对应的正交配置点都为 40,并且 MILP 方法对应的 PWA 法采用三段式进行近似处理。速度位置曲线由优化 控制输入量求得,优化控制输入量通过求解非线性列车模型得到。可以从图中看出, 在 FBS 系统中,后车的行驶受到前车的影响,如在 666m 位置处信号机对后车现实黄 灯,主要是由于求解中采用的最小追踪间隔伪 75s,小于基于 FBS 系统计算所得的最 小追踪间隔 98.4s,因此后车必须降速,来满足黄灯对应的速度限制,40km/h,也就 是 11.1m/s。

图四 基于 FBS 系统的 75s 最小追踪间隔下采用婪式算法求解 MILP 和伪谱法所得的运行速度曲线

第 15 页

西南交通大学本科毕业论文(外文翻译)

图五 基于 FBS 系统的 75s 最小追踪间隔下采用同步算法求解 MILP 和伪谱法所得的运行速度曲线

在图四中,利用婪式法求解 MILP 问题对应的能量消耗,时间误差,计算时间分 别为 112.59MJ,2.21s,96.56s,利用婪式法求解伪谱法问题对应的能量消耗,时间误 差,计算时间分别为 110.04MJ,1.43s,1381.71s,显而易见建立 MILP 模型求解所得 能量消耗和时间误差比建立伪谱法模型求解所得的值要大,但是计算时间短的多。 在图五中,利用同步法求解 MILP 问题对应的能量消耗,时间误差,计算时间分 别为 106.58MJ,1.73s,184.76s,利用同步法求解伪谱法问题对应的能量消耗,时间 误差,计算时间分别为 106.44MJ,1.68s,1521.30s,当与婪式法求解 MILP 问题对应 的能量消耗,时间误差相比,变得更小,主要是由于前后车同时进行优化,并且前车 的优化过程中考虑后车的运行,但是计算时间变得更长,是因为优化问题的规模变成 增长了一倍。当与同步法求解伪谱法问题作比较,呈现相同的差异。 可以从表三看出,在问题离散化程度相同的情况下,利用婪式法求解 MILP 问题 对应的能量消耗,时间误差比利用婪式法求解伪谱法问题对应的要大,但是伪谱法所 需的求解时间是 MILP 法求解时间的一到两个数量级。当增加 MILP 方对应的空间数 量 N 和伪谱法对应的正交配置点时,求得的能量消耗,时间误差更小,但是两种方法 之间仍然呈现相同的规律。还可以看出利用同步法求解 MILP 问题对应的能量消耗比 用婪式法求解 MILP 问题对应的能量消耗要小,但是求解时间更长,对于伪谱法的求 解结果呈现相同的规律。 6.3 基于 MBS 系统的仿真结果 对于 MBS 系统,婪式算法和同步算法的计算结果如表四所示,图六和图七中展 示了速度位置曲线和优化控制输入量, 基于与图四和图五中相同的参数设置, 同样的, 速度位置曲线以及结束时间误差等参数是基于优化控制输入量求得的,从图中可以看 出后车并没有受到前车的干扰, 主要是由于设置的最小追踪间隔为 75s 大于基于 MBS 求得的最小追踪间隔。

第 16 页

西南交通大学本科毕业论文(外文翻译)
表四基于 MBS 系统的 MILP 和伪谱法用婪式算法和同步算法的结果表现对比

图六基于 MBS 系统 75s 最小追踪间隔下采用婪式算法求解 MILP 和伪谱法所得的运行速度曲线

图七基于 MBS 系统 75s 最小追踪间隔下采用同步算法求解 MILP 和伪谱法所得的运行速度曲线

第 17 页

西南交通大学本科毕业论文(外文翻译) 在图六中,利用婪式法求解 MILP 问题对应的能量消耗,时间误差,计算时间分 别为 67.89MJ,1.68s,75.98s,利用婪式法求解伪谱法问题对应的能量消耗,时间误 差,计算时间分别为 67.42MJ,1.66s,729.24s,利用婪式法求解 MILP 问题对应的能 量消耗,时间误差比利用婪式法求解伪谱法问题对应的稍大,但是利用婪式法求解伪 谱法问题对应的计算时间是利用婪式法求解伪谱法问题对应的值大一个数量级。同样 的,利用同步算法进行求解,两种方法呈现相同的规律。当增加 MILP 方对应的空间 数量 N 和伪谱法对应的正交配置点时,求得的能量消耗,时间误差更小,但是计算时 间迅速增加。对于婪式法或者同步算法求解 MILP 问题,当 PWA 近似处理式,使用 多端式,结束时间误差也会变小。伪谱法相较于 MILP 方法,在优化控制方面表现的 更加优秀,但是由于优化问题规模的增加,其计算负担更大。 6.4 小结 上述仿真结果显示,与婪式伪谱法相比,婪式 MILP 法得到的能量消耗和时间误 差更大,但是计算时间却要小一到两个数量级,同样地,同步法下两种方法呈现相同 的差异性。更多的,婪式 MILP 法求得的能量消耗要比同步 MILP 法求得的要大,但 是计算时间更短,对于婪式伪谱法和同步伪谱法呈现相同的差异性。 7 结论和展望 在本文中,我们提出了婪式法和同步法对多车运行优化控制进行求解。对于婪式 算法,前车的运行速度曲线先被求得,基于前车的数据,对后车的运行速度曲线进行 求解。对于同步算法,对前后车的运行速度曲线同时进行求解。在固定闭塞信号系统 和移动闭塞信号系统中,由于前车造成的约束条件被引入到运行速度曲线优化问题中, 列车模型和约束条件中的非线性项目通过分段仿射公式进行近似处理为线性项目,然 后可以将问题重新建立为混合整数规划模型。还对 MILP 法与伪谱法进行两种不同求 解实际问题结果表现的比较,结果显示 MILP 法与伪谱法有差不多的表现,并且需要 的计算时间更短。 将来需要做的工作是 MILP 法与伪谱法进行更加详细及广泛的比较(利用其他非 线性求解方法,如 MINOS 和 KNITRO) ,以及动态算法,其他文献中提出的方法和框 架(例如 AMPL,APMonitor 和 ASCEBND) ,并且进行其他实际实例和情况的仿真以 及对比。

第 18 页


推荐相关:

外文翻译技巧.doc

外文翻译技巧 - 五分钟搞定 5000 字-外文文献翻译【你想要的工具都在这里】

本科毕业设计外文翻译需要注意的问题.doc

本科毕业设计外文翻译需要注意的问题 - 外文翻译需要注意的问题, 1 外文文献的

外文翻译规范要求及模版格式.doc

外文翻译规范要求及模版格式 - (二)外文翻译 每位学生必须阅读 2 篇以上(10000 字符左右)的外文材料, 应完成 5000 汉字以上 的英译汉翻译。加“外文翻译”封面...

本科毕业设计(论文)外文翻译基本规范.doc

二、毕业论文(设计)外文翻译 (一)毕业论文(设计)外文翻译的内容要求 外文翻译内容必须与所选课题相关,外文原文不少于 6000 个印刷符号。译文 末尾要用外文注明外文...

外文翻译意见表.doc

外文翻译意见表 - 指导教师意见: 该同学翻译的外文资料与毕业设计题目有一定的联

单片机系统-毕业论文外文翻译.doc

单片机系统-毕业论文外文翻译 - 西安文理学院本科毕业设计(论文) 附录 5 外文翻译 Microcomputer System Electronic systems are used fo...

毕业论文绩效管理外文翻译.doc

标签: 管理学| 绩效管理| 外文翻译|毕业论文绩效管理外文翻译_管理学_高等

如何做毕业设计的外文翻译.ppt

如何做毕业设计的外文翻译 - 如何做毕业设计的外文翻译 y2t0d 长春翻译公司

毕设外文翻译.doc

毕设外文翻译 - 对php和mysql的外文翻译... 标签: 外文翻译| 毕设外文翻译_计算机软件及应用_IT/计算机_专业资料。对php和mysql的外文翻译 西安邮电大学 毕业设计(...

外文翻译(英文).doc

外文翻译(英文) - 土木工程道路桥梁方向毕业设计外文翻译(英文部分),是个人一

外文翻译_图文.doc

外文翻译 - 西南交通大学本科毕业论文(外文翻译) 利用混合整数线性规划求解固定

毕业论文外文翻译范例_图文.pdf

毕业论文外文翻译范例 - 外文原文(一) Savigny and his Ang

外文翻译.doc

本科毕业论文外文翻译外文译文题目 (中文) : 营运能力比率分析 学专学 院:

软件开发外文翻译.doc

软件开发外文翻译 - Requirements Phase The chance

外文翻译模板.doc

外文翻译模板_IT/计算机_专业资料。杭州电子科技大学 毕业论文外文文献翻译要求

外文翻译.doc

外文翻译 - 本科毕业设计(论文) ( 2014 届) (外文翻译) 题学专班姓

外文翻译及外文原文(参考格式).doc

外文翻译及外文原文(参考格式) - 外文翻译要求: 1、外文资料与毕业设计(论文

毕业论文和外文翻译要求.doc

毕业论文和外文翻译要求 - 沈阳农业大学本科生毕业论文(设计)撰写要求与格式规范

外文翻译_图文.doc

外文翻译 - 毕业设计(论文) 外文翻译 题学专学学 目院业生号 S 形无碳小车

毕业论文外文翻译--MySQL数据库管理中心.doc

毕业论文外文翻译--MySQL数据库管理中心 - 外文原文 Management

网站首页 | 网站地图
All rights reserved Powered by 学霸学习网 www.tceic.com
copyright ©right 2010-2021。
文档资料库内容来自网络,如有侵犯请联系客服。zhit325@126.com