CN105452085B - 运转曲线制作装置、运转曲线制作方法以及运转曲线控制程序 - Google Patents
运转曲线制作装置、运转曲线制作方法以及运转曲线控制程序 Download PDFInfo
- Publication number
- CN105452085B CN105452085B CN201380078342.9A CN201380078342A CN105452085B CN 105452085 B CN105452085 B CN 105452085B CN 201380078342 A CN201380078342 A CN 201380078342A CN 105452085 B CN105452085 B CN 105452085B
- Authority
- CN
- China
- Prior art keywords
- curve
- run curve
- train
- time
- energy
- 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.)
- Active
Links
- 238000002360 preparation method Methods 0.000 title claims description 12
- 230000001133 acceleration Effects 0.000 claims description 23
- 230000037396 body weight Effects 0.000 claims description 6
- 238000004134 energy conservation Methods 0.000 claims description 6
- 238000013139 quantization Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims 2
- 230000007717 exclusion Effects 0.000 claims 1
- 230000011218 segmentation Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 45
- 238000012545 processing Methods 0.000 description 46
- 238000003860 storage Methods 0.000 description 35
- 230000004069 differentiation Effects 0.000 description 15
- 230000008569 process Effects 0.000 description 14
- 230000033001 locomotion Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 9
- 230000015654 memory Effects 0.000 description 9
- 238000013138 pruning Methods 0.000 description 7
- 238000013500 data storage Methods 0.000 description 6
- 238000005096 rolling process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 230000036651 mood Effects 0.000 description 4
- 238000011002 quantification Methods 0.000 description 4
- 230000005945 translocation Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000010422 painting Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000979 retarding effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 208000000058 Anaplasia Diseases 0.000 description 1
- UNPLRYRWJLTVAE-UHFFFAOYSA-N Cloperastine hydrochloride Chemical compound Cl.C1=CC(Cl)=CC=C1C(C=1C=CC=CC=1)OCCN1CCCCC1 UNPLRYRWJLTVAE-UHFFFAOYSA-N 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000008140 language development Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L15/00—Indicators provided on the vehicle or train for signalling purposes
- B61L15/0058—On-board optimisation of vehicle or vehicle train operation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/10—Operations, e.g. scheduling or time tables
- B61L27/16—Trackside optimisation of vehicle or train operation
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/0011—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots associated with a remote control arrangement
- G05D1/0016—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots associated with a remote control arrangement characterised by the operator's input device
-
- 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/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Train Traffic Observation, Control, And Security (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Feedback Control In General (AREA)
Abstract
【课题】高速地制作更接近最优解的运转曲线。【解决方法】实施方式的运转曲线制作装置对具有规定的车辆特性的列车制作能够在规定的出发车站与规定的到达车站之间以规定的行驶时间使累计消耗能量更少地行驶的运转曲线。最短时间运转曲线制作部根据预先设定的所述列车的车辆信息以及地面信息来制作列车在所述车站间以最短时间行驶的最短时间运转曲线,节能运转曲线制作部从最短时间运转曲线、与从出发车站出发起针对每个规定的经过时间的列车的状态相当的多个候补解中选择在规定的行驶时间累计消耗能量比较少的候补解,根据所选择的候补解来制作节能运转曲线。
Description
相关申请的引用
本申请以2013年07年19日申请的在先的日本国专利申请第2013-151104号的优先权的权益为基础并要求该权益,其内容整体通过引用被包含在本文中。
技术领域
本实施方式涉及运转曲线制作装置、运转曲线制作方法以及运转曲线控制程序。
背景技术
由于近年来的电力不足,对能量效率高的铁路也要求进一步的节约能源(以下简称为节能)。即,铁路中的节能是将空调、照明等辅助设备变更为能量效率高的设备、削减列车行驶时的能量等。特别是,从过去就研究了列车行驶时的节能,以下举出的参考文献等中展开了各种讨论。
专利文献1:日本特开平8-156794号公报
专利文献2:日本特开平7-165080号公报
非专利文献1:“市川:列車運転の最適化,日本機械学会論文集(第1部)34卷258号(1968)”(“市川:列车运转的最优化,日本机械学会论文集(第1部)34卷258号(1968)”)
非专利文献2:“山崎,国方:Dynamic Programmingによる最適なランカーブ計算,サイバネティツクス(1968)”(“山崎,国方:基于动态编程的最优的运转曲线计算,控制论(1968)”)
非专利文献3:“川島,古村,松本,高岡:最適な列車運転方式の作成および修正方法の提案,電気学会(1990)”(“川岛,古村,松本,高冈:最优的列车运转方式的制作以及修正方法的提案,电气学会(1990)”)
非专利文献4:“高木,曽根:列車群制御のための省エネルギー運転パタ一ンに関する研究,電気学会(1992)”(“高木,曾根:用于列车群控制的节能运转模式有关的研究,电气学会(1992)”)
非专利文献5:“古関,姚:省エネルギ一運転曲線最適化ヘの動的計画法の応用と数值計算上の問題点,電気学会(2002)”(“古关,姚:针对节能运转曲线最优化的动态规划法的应用和数值计算上的问题,电气学会(2002)”)
非专利文献6:“宫武,高:動的計画法を用ぃた列車運転曲線最適化問题の求解法,電気学会(2004)”(“宫武,高:使用了动态规划法的列车运转曲线最优化问题的求解方法,电气学会(2004)”)
一般与列车的运转有关的节能大多采取以下的3种方针。
(1)在从出发车站出发时以最大加速度进行加速。
(2)在限制速度内尽可能增加惰性运行(惯性行驶)。
(3)在到达车站时以最大减速度进行减速。
针对这些基本方针,目前为止提出了各种节能运转曲线制作法。为了根据这些基本方针制作运转曲线,提出了很多通过启发式(heuristic)来制作运转曲线的方法。但是,在车站间最高速度只有一个的速度限制简单的情况下、行驶期间大于最短期间的情况下,即使通过启发式也能够获得接近最优的运转曲线的运转曲线。但是,当在车站间的速度限制复杂的情况下以启发式来制作运转曲线时,有时没有解或者只得到与最优解相比消耗能量相当大的解。因此,希望有如下方法:即使被给予的速度限制、坡度复杂,也能够获得接近最优解的解。
附图说明
图1是运转曲线制作装置的概要结构框图。
图2是由个人计算机构成运转曲线制作装置的情况下的概要结构框图。
图3是地面信息数据的格式以及数据例的说明图。
图4是车辆信息数据的格式以及数据例的说明图。
图5是表示组合最优化问题中的生成操作和剪枝操作的图。
图6是实施方式的处理流程图。
图7是最短时间运转曲线的制作处理流程图。
图8是最大加速度以及最大减速度的可利用条件的说明图。
图9是加速极限曲线的制作处理的处理流程图。
图10是减速极限曲线的制作处理的处理流程图。
图11是最短时间运转曲线的制作处理的处理流程图。
图12是节能运转曲线的制作处理的处理流程图。
图13是横轴取时刻、纵轴取列车的位置而成的(x-t)图形。
图14是速度条件的说明图。
图15是量化处理的说明图。
图16是运转曲线制作结果存储部的存储数据的格式以及数据例的说明图。
图17是将基于运转曲线制作结果存储部的存储数据的运转曲线进行视觉化而成的图。
图18是表示显示在显示器的运转曲线的例子的图。
(符号说明)
10:运转曲线制作装置;11:地面信息存储部;12:车辆信息存储部;13:运转曲线输入处理部;14:最短时间运转曲线制作部;15:节能运转曲线制作部;16:运转曲线制作结果存储部;17:运转曲线制作结果显示部;21:MPU;22:ROM;23:RAM;24:外部存储装置;25:IC卡;26:IC卡R/W;27:键盘;28:显示器;30:地面信息数据;31:地面ID数据;32:车站关联数据;33:速度限制数据;34:坡度数据;35:车站间信息数据;36:车站间距离数据;37:起点数据;38:终点数据;39:速度数据;40:起点数据;41:终点数据;42:坡度值数据;50:车辆信息数据;51:车辆ID数据;52:种类数据;53:车体重量数据;54:乘车率数据;55:编制车辆数数据;56:列车长度数据;57:档位特性ID数据;T:行驶时间;Tmin:最短行驶时间。
具体实施方式
一个实施方式的运转曲线制作装置对具有规定的车辆特性的列车制作能够在规定的出发车站和规定的到达车站之间以规定的行驶时间进一步减少累计消耗能量地行驶的运转曲线。
另外,最短时间运转曲线制作部根据预先设定的所述列车的车辆信息以及地面信息来制作所述列车以最短时间(最短行驶期间)在所述车站间行驶的最短时间运转曲线。
由此,节能运转曲线制作部从最短时间运转曲线以及与在从出发车站出发起针对每个规定的经过时间的列车状态相当的多个候补解中选择在规定的行驶时间累计消耗能量比较少的候补解,根据所选择的候补解来制作节能运转曲线。
一个实施方式的运转曲线制作装置能够提供一种能够高速地制作更接近最优解的运转曲线的运转曲线制作装置。
接着,进一步参照附图详细地说明详细的实施方式。
图1是运转曲线制作装置的概要结构框图。
运转曲线制作装置10具备有:地面信息存储部11,预先存储运转曲线制作所需的后述的地面信息数据;车辆信息存储部12,预先存储运转曲线制作所需的后述的车辆信息数据;运转曲线输入处理部13,输入车辆信息、车站间信息、行驶时间信息等各种数据,并且根据所输入的数据从地面信息存储部11获取运转曲线制作所需的地面信息数据,并从车辆信息存储部12获取车辆信息数据;最短时间运转曲线制作部14,制作用于通过最短时间在以根据从运转曲线输入处理部13输入的各种数据所指定的车辆条件而指定的车站间行驶的最短时间运转曲线;节能运转曲线制作部15,根据最短时间运转曲线来制作在所指定的行驶时间在所指定的车站间行驶并且能够实现节能的节能运转曲线;运转曲线制作结果存储部16,存储最短时间运转曲线制作部14制作出的最短时间运转曲线以及节能运转曲线制作部15制作出的运转曲线;以及运转曲线制作结果显示部17,显示节能运转曲线制作部15制作出的运转曲线制作结果。
在本实施方式中,假定了用通用的个人计算机等来进行运转曲线制作。
图2是由个人计算机构成运转曲线制作装置的情况下的概要结构框图。
如图2所示,运转曲线制作装置10具备有:MPU 21,构成为控制运转曲线制作装置整体的微计算机,作为运转曲线输入处理部13、最短时间运转曲线制作部14以及节能运转曲线制作部15发挥功能;ROM 22,非易失性地存储包含控制程序的各种数据;RAM 23,用作MPU 21的工作区,暂时地存储各种数据;外部存储装置24,构成为硬盘驱动装置(HDD)等,作为地面信息存储部11、车辆信息存储部12以及运转曲线制作结果存储部16发挥功能;IC卡读取器/写入器(R/W)26,能够对作为半导体存储装置的IC存储器进行各种数据的读出以及写入,能够作为运转曲线输入处理部13发挥功能;键盘27,设为操作者能够输入各种信息,能够作为运转曲线输入处理部13发挥功能;显示器28,构成为液晶显示器、有机EL显示器等,能够作为运转曲线制作结果显示部17发挥功能;以及通信接口(I/F)29,进行经由未图示的因特网、LAN等通信网络的通信的接口动作,能够作为运转曲线输入处理部13发挥功能。
这里,说明地面信息存储部11以及车辆信息存储部12中的存储数据的格式以及数据例。
首先说明存储在地面信息存储部11中的地面信息数据的格式以及数据例。
图3是地面信息数据的格式以及数据例的说明图。
存储在地面信息存储部11中的地面信息数据30大致区分为用于确定地面信息数据的地面ID数据31、包含与该地面信息数据30相对应的车站关联的信息的车站关联数据32、作为与车站关联数据32相对应的与车站间的速度限制有关的信息的速度限制数据33、以及作为与车站关联数据相对应的与车站间的坡度有关的信息的坡度数据34。
车站关联数据32包含有与地面ID数据31相对应的车站间信息数据35、以及与地面ID数据31相对应的车站间距离数据36。
更具体地说,在图3的例子的情况下,车站间信息数据35表示从A车站向B车站,车站间距离数据36表示从A车站向B车站时的车站间距离为1500m。
速度限制数据33与由车站关联数据32所确定的车站间的速度的限制有关,包含有表示施加了速度限制的速度限制区间的起点的起点数据37、表示速度限制区间的终点的终点数据38以及表示该速度限制区间中的上限速度的速度数据39。
起点数据37表示在将出发车站(在上述的例子的情况下为A车站)设为0m的情况下从出发车站至速度限制区间的起点位置为止的距离。
另外终点数据38表示在将出发车站设为0m的情况下从出发车站至速度限制区间的终点位置为止的距离。
另外,速度数据39保存了相对应的速度限制区间中的上限速度值。
具体地说,在图3的例子的情况下,速度限制区间在从作为出发车站的A车站至作为到达车站的B车站之间具有4个区间,第1区间为从作为出发车站的A车站的位置(距离0m)至500m的位置为止的区间,上限速度值为25m/sec。
同样地,第2区间为从离出发车站500m的位置至750m的位置为止的区间,上限速度值为23m/sec。第3区间为从离出发车站1000m的位置至1300m的位置为止的区间,上限速度值为20m/sec。第4区间为从离出发车站1300m的位置至1500m的位置(作为终点站的B车站的位置)为止的区间,上限速度值为20m/sec。
坡度数据34包含从出发车站至到达车站为止的全部的坡度信息,具备有表示坡度连续地固定的区间的起点的起点数据40、表示该坡度连续地固定的区间的终点的终点数据41以及表示该坡度固定区间中的坡度值的坡度值数据42。
具体地说,在图3的例子的情况下,坡度固定区间在从A车站至B车站之间具有3个区间,第1区间为从作为出发车站的A车站的位置(距离0m)至1000m的位置为止的区间,坡度值为0‰。
同样地,第2区间为从离出发车站1000m的位置至1300m的位置为止的区间,坡度值为2‰。第3区间为从离出发车站1300m的位置至1500m的位置(作为到达车站的B车站的位置)为止的区间,坡度值为0‰。
接着说明存储在车辆信息存储部12中的车辆信息数据的格式以及数据例。
图4是车辆信息数据的格式以及数据例的说明图。
存储在车辆信息存储部12中的车辆信息数据50包含有用于确定铁路车辆的车辆ID数据51、表示铁路车辆的驱动种类的种类数据52、表示车体重量的车体重量数据53、表示在运转曲线制作时假定的乘车率的乘车率数据54、表示铁路车辆的编制车辆数的编制车辆数数据55、表示列车长度的列车长度数据56、以及用于确定搭载于该铁路车辆的主控制器的档位(notch)特性的档位特性ID数据57。
具体地说,在图4的例子的情况下,车辆ID数据51=5,种类数据52表示铁路车辆的驱动种类为VVVF方式,车体重量数据53表示车体重量为10t,乘车率数据54表示在运转曲线制作时假定的乘车率为乘车率50%,编制车辆数数据55表示编制车辆数为10辆,列车长度数据56表示列车长度为200m,档位特性ID数据57表示搭载于该铁路车辆的主控制器的档位特性ID=5。
接着在说明实施方式的动作之前,说明实施方式的原理性的思路。
列车的节能运转曲线具有作为制约包含行驶阻力以及坡度阻力等的运动方程式、和由从出发起的经过时间t、速度v及离出发车站的距离x表示的将出发车站处的状态(初始状态)设为(t,v,x)=(0,0,0)、将到达车站处的状态(终端状态)设为(t,v,x)=(T,0,X)的2点边界值。
这里,T是从出发车站至到达车站所需的经过时间,X是出发车站与到达车站之间的距离。
另一方面,评价函数通过求解设为在出发车站与到达车站之间消耗的能量总和E(累计消耗能量)的最优控制问题来求得。此时,指定列车的运动的控制变量为指定加速度、减速度的档位。当前,广泛普及的列车通过具有多个的档位来实现列车的加减速,运转曲线能够视为在各时刻中选择哪个档位来运转的多级决策问题。
另外,一般求解多级决策问题时成为课题的是候补解的组合爆炸。
以下,简单说明组合爆炸。
例如,在档位的选项有10个、行驶时间为100秒的运转曲线以最小控制时间单位1秒来制作的情况下,不得不从10100个候补解中找出满足制约的解。从中一个不漏地找出满足制约的解在现实的时间的范围内是不可能的。
因此,作为解决该组合爆炸的课题的方法,提出了使用了动态规划法的最优解的搜索方法。
使用了动态规划法的运转曲线的制作是将列车的(时间、速度、位置)的状态定义在格子点上,并在该格子点上搜索累计消耗能量成为最小的路径。此时,在有到达某状态(t,v,x)的多个路径的情况下,只存储累计的消耗能量最小的路径,并继续搜索。在使用了动态规划法的情况下,预先将状态空间分割为格子状,但是一般不太可能从格子点向格子点移动。因此,在不到格子点的情况下施以某种校正来移动到格子点。在该方法中,存在当向格子点的校正变多时,与最优解相比误差变大的问题。
因此,在本实施方式中提出了如下方法:不对多级决策问题的候补解进行全枚举,而是在各级(时刻)中对是否能够成为允许解的候补进行检查(剪枝),同时每前进一级则生成候补解。
这里,剪枝是用于求解运筹学的领域中公知的组合最优化问题的通用方法。该方法是如下方法:在生成操作中生成候补解,在剪枝操作中对不能成为最优解的候补解进行剪枝,并从候补解的列表删除。重复生成操作、剪枝操作,同时搜索最优解。是接近全枚举的思路,但是不能成为候补解的部分不进行搜索这点比全枚举法效率高。
图5是表示组合最优化问题中的生成操作和剪枝操作的图。
具体地说,在经过时间t=0时状态为速度V0、位置X0的情况下,在经过最小控制时间单位后的经过时间t=1时所认为的状态是第1状态为速度V1、位置X1、第2状态为速度V2、位置X2、第3状态为速度V3、位置X3,但是在进行了该第3状态或者从该第3状态迁移了的经过时间t=2时的三个状态(速度V10、位置X10)、(速度V11、位置X11)、(速度V12、位置X12)的评价的情况下,当得不到值得进行迁移的评价值(不会有候补解)时,进行剪枝而不进行以后的搜索。但是,剪枝操作依赖于求解的问题而变化,必须根据问题来适当地设定。
以下,说明包含生成操作和剪枝操作的具体动作。
图6是实施方式的处理流程图。
首先,从运转曲线输入处理部13获取地面信息、车辆信息以及行驶时间T(步骤S101)。
即,在制作运转曲线时,需要输入在路线的哪个车站间(需要区别从A车站向B车站前进的情况和从B车站向A车站前进的情况)用哪个车辆以多少行驶时间T行驶。
这里,关于车辆信息、车站间信息、行驶时间信息等各种数据的输入,考虑通过IC卡R/W26来读入存储在IC卡25中的设定文件,或者通过经由通信I/F 29接收到的设定文件来输入、或经由键盘27等直接输入等。
当车辆信息、车站间信息、行驶时间信息等各种数据被输入到运转曲线输入处理部13中时,运转曲线输入处理部13参照地面信息存储部11以及车辆信息存储部12,根据所输入的车辆信息以及车站间信息,获取运转曲线制作所需的数据。
接着,作为最短时间运转曲线制作部14发挥功能的MPU 21从运转曲线输入处理部13接受运转曲线制作所需的数据,制作所输入的车站间的最短时间的运转曲线(步骤S102)。
图7是最短时间运转曲线的制作处理流程图。
最短时间运转极性制作处理大致区分为具备加速极限曲线制作处理(步骤S201)、减速极限曲线制作处理(步骤S202)以及基于剪枝算法的最短时间运转曲线制作处理(步骤S203)而构成。
以下,具体地进行说明。
为了制作最短时间运转曲线,将求解具有初始条件和终端条件的最短时间控制的解。一般,最短时间控制通过被称为bang-bang(乓乓)控制的最大加速度或者最大减速度来控制。然而,在列车的运转时被施加主要基于位置的速度限制,所以不一定只能够通过最大加速度和最大减速度进行控制。因此求出能够以最大加速度和最大减速度使列车移动的范围。
图8是最大加速度以及最大减速度的可利用条件的说明图。
在列车的运转时,能够使用最大加速度和最大减速度的情况是如图8所示那样情形case1~case4(情形1~情形4)这4种情况。
Case1:从出发车站出发时使用最大加速度
Case2:到达到达车站时使用最大减速度
Case3:当限制速度在行进方向增加时使用最大加速度
Case4:当限制速度在行进方向减少时使用最大减速度
在图8中,向右上升的线是以最大加速度使列车移动的情况,向左下降的线是以最大减速度使列车移动的情况。
在一个车站间,Case1的状态和Case2的状态必然只各存在一个。
与此相对,Case3的状态和Case4的状态有可能在一个车站间存在多个,因此需要对它们的所有计算加速极限曲线和减速极限曲线。
首先,使用图8说明图9和图10所示的变量ST。首先,以如下的结构体数组来给出速度限制。
struct velocity_limitation{
double position;//速度限制的开始位置
double velocity;//上述开始位置以后的限制速度
};
在图8中成为
struct velocity_limitation ST[3]
={{0,vmax},{x1,v2},{x2,vmax}};
这里,如ST[0].position<ST[1].position<ST[2].position<ST[3].position那样,设为结构体的指针(index)越小则开始位置越小。
接着检测成为加速极限曲线和减速极限曲线的目标的位置。默认设为
加速极限曲线的目标位置=0;//Case1
加速极限曲线的目标速度=0;//Case1
减速极限曲线的目标位置=X;//Case2(从出发车站至到达车站为止的距离)
减速极限曲线的目标速度=0;//Case2
除此之外将与以下的条件一致的位置设定为加速极限曲线和减速极限曲线的目标位置。
在图9中,针对加速极限曲线制作所需的目标位置和目标速度,将从默认的车站出发的时间点的Case1和有多个速度限制的情况下满足上述条件的所有目标位置、目标速度设定在数组ST。
同样地在图10中,针对减速极限曲线制作所需的目标位置和目标速度,将到达默认车站的时间点的Case:2和有多个速度限制的情况下满足上述条件的所有目标位置、目标速度设定在数组ST。
具体地说,以图8为例,
在加速极限曲线的情况下,成为
struct velocity_limitation ST[2]
={{0,0},{x2,v2}};
在减速极限曲线的情况下,成为
struct velocity_limitation ST[2]
={{x1,v2},{{X,0}};
在加速极限曲线制作处理(步骤S201)中,首先计算出最初使用该最大加速度和最大减速度时的(速度、位置)。
将图8所示的情形case1~case4中的各曲线分别称为加速极限曲线、减速极限曲线。
另外,加速极限曲线能够通过提供使用最大加速的速度和位置的初始条件来容易地求出。
图9是加速极限曲线的制作处理的处理流程图。
首先,初始化变量s,设为初始值0(步骤S401)。
接着设定加速极限曲线的目标位置、即加速经过时刻i=0、加速极限曲线中的初始速度v(0)=ST[s].velocity以及初始位置x(0)=ST[s].position(步骤S402)。
接着对时刻i累加1(步骤S403)。
接着用最大加速档位(与能够获得最大的加速度的档位位置相当)求解运动方程式来计算时刻i时的状态(v(i),x(i))(步骤S404)。
接着判别时刻i时的速度v(i)是否超过限制速度vmax、即
v(i)>vmax
或者判别时刻i时的位置x(i)是否超过目标位置ST[s].position(步骤S405)。
在步骤S405的判别中时刻i时的速度v(i)不超过限制速度vmax且时刻i时的位置x(i)不超过目标位置ST[s].position的情况下(步骤S405;否),使处理再次跳转到步骤S403,以下重复相同的处理。
在步骤S405的判别中时刻i时的速度v(i)超过限制速度vmax、或者时刻i时的位置x(i)超过目标位置ST[s].position的情况下(步骤S405;是),判别是否到达加速极限曲线的目标位置、即判别是否成为
s=sizeof(ST[s])/sizeof(ST[0])(步骤S406)。
在步骤S406的判别中未到达加速极限曲线的目标位置的情况下(步骤S406;否),对变量s累加1(步骤S407),使处理再次跳转到步骤S402,以下重复相同的处理。
在步骤S406的判别中到达了加速极限曲线的目标位置的情况下(步骤S406;是),结束加速极限曲线的制作处理。
接着与加速极限曲线的计算同样地计算减速极限曲线(步骤S202)。
在本实施方式中,在减速极限曲线的情况下用迭代方法来搜索满足终端条件的初始条件。一般,与终端条件完全一致是几乎不可能的,因此对终端条件给予允许范围(ε以及δ),求出满足该范围的初始条件。
图10是减速极限曲线的制作处理的处理流程图。
首先,将变量s初始化为初始值0、将变量a初始化为初始值0、设为变量b=ST[s].position、变量Distance=ST[s].position(步骤S501)。
接着设定减速极限曲线的目标位置、即减速经过时刻i=0、减速极限曲线中的初始速度v(0)=vmax以及初始位置x(0)=(a+b)/2.0(步骤S502)。
接着对时刻i累加1(步骤S503)。
然后,用最大减速档位求解运动方程式来计算时刻i中的列车的状态(步骤S504)。
接着判别时刻i时的速度v(i)是否满足与允许值ε相对应的下式(步骤S505)。
0≤v(i)<ε
在步骤S505的判别中,在
ε≤v(i)
的情况下(步骤S505;否),再次使处理跳转到步骤S503,以下进行相同的处理。
在步骤S505的判别中,在
0≤v(i)<ε
的情况下(步骤S505;是),判别时刻i时的位置x(i)是否满足与允许值δ相对应的下式(步骤S506)。
|x(i)-Distance|<6
在步骤S506的判别中,在
|x(i)-Distance|≥δ
的情况下(步骤S506;否),判别是否为
x(i)<Distance
(步骤S507)。
在步骤S507的判别中,在
x(i)<Distance
的情况下(步骤S507;是),设为
a=(a+b)/2.0
使处理再次跳转到步骤S502,以下进行相同的处理(步骤S508)。
在步骤S507的判别中,在
x(i)≥Distance
的情况下(步骤S507;否),设为
b=(a+b)/2.0
使处理再次跳转到步骤S502,以下进行相同的处理(步骤S509)。
另一方面,在步骤S506的判别中,在
|x(i)-Distance|<δ
的情况下(步骤S506;是),判别是否到达减速极限曲线的目标位置、即判别是否成为
s=sizeof(ST[s])/sizeof(ST[0])
(步骤S510)。
在步骤S510的判别中未到达减速极限曲线的目标位置的情况下(步骤S510;否),对变量s累加1(步骤S511),使处理再次跳转到步骤S502,以下重复相同的处理。
在步骤S510的判别中到达了减速极限曲线的目标位置的情况下(步骤S510;是),结束减速极限曲线的制作处理。
另外,为了计算减速极限曲线,有时使用将时间反向地推进来求解运动方程式的反向索引的方法,但是在存在坡度的情况下有时在反向索引的方法中与对时间正向地求解的一般的微分方程式的数值解法相比解有偏差。
与此相对,根据本实施方式的方法,使用用于更正确地求解微分方程式的方法,与反向索引的方法相比能够获得更高精度的减速极限曲线。
接着根据剪枝算法来制作最短时间运转曲线(步骤S203)
图11是最短时间运转曲线的制作处理的处理流程图。
这里,说明图11所示的变量previous_state、变量current_state以及变量middle_state。
在本实施方式中,作为实现运转曲线制作处理的控制程序,假定通过C++、JAVA(注册商标)等高级语言开发的控制程序。
特别是在C++、JAVA中,能够使用根据要素数来改变数组的大小的可变长数组。因此,在本实施方式中,设为使用可变长数组来表现各时刻时的候补解。
具体地说,在C++等中,以
struct phase{
double position;//位置
double velocity;//速度
double energy;//瞬时的能耗量
double total_energy;//累计的能耗量
int notch[500];//各时刻时的档位编号
int change_count;//档位切换次数
int notch_continuous_time;//当前的档位的持续时间};
那样的结构体来表现候补解。
另外,只要设为
vector<struct phase>previous_state;
vector<struct phase>current_state;
vector<struct phase>middle_state;
来表现候补解即可。
首先,作为出发车站的初始条件将经过时间t=0、速度v=0、位置x=0设定于变量previous_state(步骤S601)。
接着针对经过时间t+1的各个候补解,求解所有档位中的运动方程式(微分方程式),制作(经过时间t的候补解个数)×(档位选择数)个经过时间t的候补解,将该候补解保存为变量current_state,将变量previous_state从存储器中释放(清空)(步骤S602)。
这里,具体地说明经过时间t+1时的候补解的制作方法。
首先,设为在经过时间t时有N个候补解。
如果将N个中的第n个候补解的状态记述为(t,vn(t),xn(t)),则将第n个的状态(t,vn(t),xn(t))设为数值计算上的初始值,使用第k个档位NOTCH(k)所给予的加速度或者减速度,通过4次的龙格-库塔法(Runge-Kutta method)等数值计算等来计算经过时间t+1的状态。这里将4次的龙格-库塔法举为数值计算例,但是运动方程式的数值解法不限于此。
由此,经过时间t+1的候补解生成NK个。
接着,针对变量current_state的各个候补解的速度,检查是否小于加速极限曲线和减速极限曲线以及基于位置的速度限制的值,只将小的值保存到变量middle_state,变量current_state从存储器中释放(清空)(步骤S603)。
对变量middle_state的位置以降序(从大到小的顺序)排序,将高位K个保存到变量previous_state,将变量middle_state从存储器中释放(步骤S604)。
这里,剩余高位K个是为了防止变得没有候补解。K的值越大,候补解变得没有的可能性越低。然而K的值过大时计算时间变大,因此实用上认为可以是K=20左右。
因此,在以下的说明中考虑为K=20。
进而与减速极限曲线同样地,直到终端条件具有允许幅度的条件成立为止(步骤S605;是)进行重复,针对出现满足该条件的候补解的经过时间t,将t+1设为最短行驶期间(步骤S607)。
在步骤S605的判别中不满足条件的情况下(步骤S605;否),对经过时间t累加(增加)1(步骤S606),使处理再次跳转到步骤S602,以下重复相同的处理。
接着说明节能运转曲线的制作处理。
图12是节能运转曲线的制作处理的处理流程图。
首先,为了制作由运转曲线输入处理部13输入的行驶时间T的节能运转曲线,获取最短时间运转曲线制作部14的结果,将最短时间运转曲线的结果设置为上限值、将最短时间运转曲线的结果在时间方向上移动了行驶时间T的结果设置为下限值(步骤S700)。
首先,作为出发车站的初始条件将经过时间t=0、速度v=0、位置x=0设定于变量previous_state(步骤S701)。
接着针对经过时间t+1的各个候补解,求解所有档位中的运动方程式(微分方程式),制作(经过时间t的候补解个数)×(档位选择数)个经过时间t的候补解,将该候补解保存为变量current_state,将变量previous_state从存储器中释放(清空)(步骤S702)。
接着,将在变量current_state中满足解的条件(位置、速度、档位切换次数、档位持续时间)的候补解保存到变量middle_state,变量current_state从存储器中释放(清空)(步骤S703)。
将经过时间t+1时的变量middle_state以状态(量化的位置、量化的速度、档位切换次数)进行排序。然后将在视为相同状态的状态下累计能耗量最小的值选择为候补解之一,并保存到变量previous_state,将变量middle_state从存储器中释放(步骤S704)。
接着判别是否为经过时间t=行驶时间T(步骤S705)。
在步骤S705的判别中还是经过时间t<行驶时间T的情况下(步骤S705;否),对经过时间t累加(增加)1,使处理再次跳转到步骤S702,以下重复相同的处理。
在步骤S705的判别中经过时间t=行驶时间T的情况下(步骤S705;是),将保存在变量previous_state中的候补解中累计能耗量最小的候补解选择为最优解(步骤S707)。
接着详细地说明检查候补解的条件。
为了简化,这里将计算在某车站间以行驶时间Tmin+10行驶的情况下的运转曲线举为例子进行说明。这里最短行驶时间Tmin是在该车站间最短地行驶的情况下的行驶时间(最短行驶时间)。
考虑以相对于最短行驶时间Tmin长出10秒的行驶时间Tmin+10在车站间行驶。
图13是横轴取时刻、纵轴取列车的位置而成的(x-t)图形。
上面的图形表示以最短行驶时间Tmin行驶的情况下的位置x和经过时间t的关系,下面的图形表示从出发车站等待10秒后以最短行驶时间Tmin行驶的情况下的位置x和经过时间t的关系。
这里,最短行驶时间Tmin是将车辆的性能最大限度地运用而行驶的时间,如果是相同的车辆特性,则表示行驶不能比它更快。
即,如果车辆特性被决定,则最短行驶时间Tmin被唯一地求出。
因而,在某个时刻,相比于下面的图形,位置x更接近出发车站的情况下,表示在行驶时间Tmin+10不能到达到达车站(下一站)。
因此,可以说如下内容。
当在某个车站间以行驶时间Tmin+10行驶时,将在该车站间以最短行驶时间Tmin行驶的情况下的经过时间t时的位置设为Upper[t]、将最短行驶时间下的(t,Upper[t])在时间轴方向平行移动了10秒的图形设为(t,Lower[t])。其中,在t<10时设为Lower[t]=0。
如果将行驶时间Tmin+10下的能耗量变得最小的解的轨道设为(t,v(t),x(t)),则下式成立
Lower[t]≤x(t)≤Upper[t]
这里,t:0,1,2,...,Tmin+10。
图14是速度条件的说明图。
一般在铁路业界中,将该图形称为运转曲线。这里,在(v-x)图形中,将与行驶时间T相对应的位置x处的速度V记为VT[x]。由此,从先行研究的结果可知具有如下特性。
0≤VT[x]≤VTmin[x]
这里,0≤x≤X。
即,在以(v-x)图形来看的情况下,行驶时间T的运转曲线存在于比最短时间运转曲线更靠内侧。
但是,最短时间运转曲线的结果也只知道离散经过时间t时的列车的速度、位置,难以与VTmin[x]同样地计算。
因此,代替上述的VT[x]和VTmin[x]的关系式而使用下式。首先,针对行驶时间Tmin的最优运转曲线,将经过时间t时的速度以及位置分别设为VTmin(t)、XTmin(t)。针对经过时间t时的候补解(t,v(t),x(t)),找出满足XTmin(α-1)≤x(t)<XTmin(α)的α(自然数)。对找出的α判定以下式子是否成立。
[式1]
这是使用给予基于位置的速度上限的行驶时间Tmin的结果来对VTmin(x)进行线性插值而近似得到的。
另外,在单纯地用上述的方法求出运转曲线的情况下,生成累计能耗量少的运转曲线,因此不一定生成对于驾驶员而言操作性良好的运转曲线。
例如,为了沿着所生成的运转曲线进行车辆的运转,有可能必须在各时刻细微地操作的档位,有可能不实用。
因此,在本实施方式中,作为候补解的选择条件,对至该时刻为止的档位切换次数进行计数,如果达到固定值以上则从候补解中排除。例如,如果将针对某个车站间的行驶到达下一站为止档位的切换次数设为不高于10次,则只将满足存储档位切换次数的变量current_state[i].change_count≤10的条件的变量current_state[i]保存到变量middle_state中。
另外,与档位切换次数同样地,档位持续条件是用于在短的时间间隔不切换档位的条件。如果将档位持续条件设为5秒以上,则只将满足保存档位持续时间的变量current_state[i].notch_continuous_time≤5的条件的变量current_state[i]保存到变量middle_state中。
因此,在上述的步骤S703中,对所有的候补解使用位置条件、速度条件、档位切换次数条件、档位持续时间条件来进行判定,只将满足这些条件的变量保存到变量middle_state中。
进而在步骤S704中,将在步骤S703中保存了满足各个条件的候补解的变量middle_state中经过时间t+1中的状态(t,v(t),x(t))与速度和位置相关地进行量化,将被量化的速度和位置、档位切换次数、档位持续时间相同的候补解进行分组,只将进行了分组的值中累计能耗量最小的值留作候补解。
另外,在本实施方式中,也能够容易地进行关于瞬时的能耗量不超出上限值地选择档位这样的扩展。
这里,更详细地说明量化处理。
图15是量化处理的说明图。
位置和速度是连续值,因此为了进行量化,有如下方法:例如关于位置通过1[m]、而对于速度通过0.1[m/sec]这样的固定幅度来对值进行量化的方法,或关于位置,时刻t时的上限值Upper[t]和下限值Lower[t]确定,因此通过对Upper[t]和Lower[t]进行N分割来进行量化的方法。这里以按固定幅度量化位置和速度为前提,进行以后的说明。
在本实施方式中,使用位置条件、速度条件、档位切换次数条件、档位持续时间条件这4个,因此以4个轴来定义了候补解的状态,但是难以图示4维以上的状态,因此这里说明位置条件、速度条件、档位切换次数条件这3个。
对某个经过时间t时的候补解进行了剪枝操作的结果设为如图15(a)所示的状态。
量化处理在(v(t),x(t),n(t))空间上如图15(a)所示分割为格子状(矩形状),相同区域内视为相同状态(条件)。这里n(t)设为经过时间t时的档位切换次数。
另外,如果在相同区域内存在多个候补解,则在同一区域内只将累计能耗量最小的1个留作候补解。其结果,能够对更多的进行剪枝,能够进一步缩短处理时间。
在图15(b)中,黑点表示在各区域内留作候补解的值。
在步骤S705中,判定经过时间t是否与行驶时间T一致,如果不一致(步骤S705;否),则进入到步骤S706,使经过时间t增加,再次使处理跳转到步骤S702。
在步骤S705中,如果经过时间t和行驶时间T一致(步骤S705;是),则从候补解中将累计能耗量最小的候补解选择为最优解(步骤S707)。
然后,在步骤S707的处理中,为了最终选择累计能耗量小的候补解,只要按照升序(从小到大的顺序)排列累计能耗量并从中选择即可。
另外,在结束时间点有多个候补解的情况下,将累计能耗量+penalty(v,x)最小的解设为最优解。
这里penalty(v(T),x(T))表示经过时间t时的状态(T,v(T),x(T))和终端状态(T,0,X)的误差的罚(penalty)。一般,用运动方程式计算速度和位置的情况下,使用浮点型的变量,因此几乎没有速度和位置与终端状态完全一致的情况。因此,准备函数penalty(v(T),x(T)),包含误差在内地评价累计能耗量。
penalty(v(T),x(T))的决定方法与非专利文献6同样地采用下式。
[式2]
penalty(v(T),x(T))={v(T)}2+{x(T)-X}2
接着说明运转曲线制作结果存储部16的存储数据的数据格式。
图16是运转曲线制作结果存储部的存储数据的格式以及数据例的说明图。
在运转曲线制作结果存储部16中,针对确定车辆(车辆信息)的车辆ID、确定地面信息的地面ID和行驶时间T的每个组合,保存如图16所示的数据。
更详细地说,除了车辆信息、地面信息以及行驶期间的信息以外,记载有该车辆ID和地面ID中的最短行驶期间、相对于所设定的行驶期间中的时刻的速度、位置、瞬时电力、累计电力以及档位号。假如是输入到运转曲线输入处理部13的车辆ID、地面ID以及行驶时间T的数据已经保存在运转曲线制作结果存储部16中的状态,则即使不制作运转曲线也可以向运转曲线制作结果显示部17转交该数据,在运转曲线制作结果显示部17中显示该结果。
图17是对基于运转曲线制作结果存储部的存储数据的运转曲线进行了视觉化的图。
另外,作为运转曲线制作结果显示部17发挥功能的显示器28显示在节能运转曲线制作部15中制作出的行驶时间T的运转曲线的结果。
图18表示显示于显示器的运转曲线的例子。
铁路业界中的运转曲线是横轴取位置、纵轴描绘了速度和位置的曲线。
在本实施方式中,也如图18所示地显示运转曲线。
另外,在作为运转曲线制作结果显示部17发挥功能的显示器28中,还与其它行驶时刻下的累计能耗量的结果相比较,还能够比较累计能耗量的削减率。在本实施方式中,为了简化,只关注某路线的1个车站间而显示了运转曲线的结果,但是能够容易地从各车站间的结果制作路线整体的运转曲线。
此外,在所输入的行驶时间T小于最短行驶时间Tmin的情况下,在物理上行驶是不可能的,因此只要将在该输入的行驶时刻下不可能行驶的情况显示于运转曲线制作结果显示部17即可。
另外,在所输入的行驶时间T长于最短行驶时间Tmin的情况下,当在运转曲线制作结果存储部16中没有运转曲线的数据时,也可以将运转曲线制作结果存储部16没有数据的情况通知给运转曲线输入处理部13,从而运转曲线输入处理部13开始运转曲线制作处理。
如以上所说明,根据本实施方式,能够高速地制作更接近最优解的运转曲线。
进而能够自动地制作驾驶员的操作性更高的运转曲线。
在以上的说明中,将运转曲线制作装置10实现为个人计算机,但是还能够构成为实现为服务器、或实现为网络上的所谓的云。
在上述结构中,运转曲线制作结果显示部不需要一体地设置,还能够构成为配置在用户侧。
由本实施方式的运转曲线制作装置所执行的控制程序以可安装的形式或者可执行的形式的文件记录在CD-ROM、软盘(FD)、CD-R、DVD(Digital Versatile Disk:数字多功能盘)等计算机可读取的记录介质中来提供。
另外,也可以构成为通过将由本实施方式的运转曲线制作装置执行的控制程序保存在连接到因特网等网络的计算机上并经由网络下载来提供。另外,也可以构成为经由因特网等网络来提供或者发布由本实施方式的运转曲线制作装置执行的控制程序。
另外,也可以构成为将本实施方式的运转曲线制作装置的控制程序预先载入在ROM等中来提供。
由本实施方式的运转曲线制作装置执行的控制程序成为包含上述的各部(运转曲线输入处理部、最短时间运转曲线制作部、节能运转曲线制作部、运转曲线制作结果显示部)的模块结构,作为实际的硬件,CPU(处理器)从上述存储介质读出控制程序来执行,由此上述各部被加载到主存储装置上,在主存储装置上生成运转曲线输入处理部、最短时间运转曲线制作部、节能运转极性制作部、运转曲线制作结果显示部。
说明了本发明的几个实施方式,但是这些实施方式是作为例子提示的,并不意图限定发明的范围。这些新的实施方式能够以其它的各种方式来实施,能够在不脱离发明的精神的范围内进行各种省略、替换、变更。这些实施方式、其变形包含在发明的范围、精神中,并且包含在权利要求书所记载的发明及其等同的范围内。
Claims (14)
1.一种运转曲线制作装置,对具有规定的车辆特性的列车制作能够在规定的出发车站与规定的到达车站之间以规定的行驶时间使累计消耗能量减少地行驶的运转曲线,该运转曲线制作装置具备:
最短时间运转曲线制作部,根据预先设定的所述列车的车辆信息以及地面信息来制作所述列车在所述车站间以最短时间行驶的最短时间运转曲线;以及
节能运转曲线制作部,从所述最短时间运转曲线以及与从所述出发车站出发起针对每个规定的经过时间的所述列车的状态相当的用于定义列车的状态的多个候补解中选择在所述规定的行驶时间所述累计消耗能量少的候补解,根据所选择的候补解来制作节能运转曲线,
所述节能运转曲线制作部从由所述最短时间运转曲线制作部制作出的最短时间运转曲线的结果求出经过时间和位置的关系的上限值,从将所述最短时间运转曲线的结果在时间方向移动了规定行驶时间量的结果求出为经过时间和位置的关系的下限值,从该上限值和下限值中选择所述累计消耗能量少的候补解。
2.根据权利要求1所述的运转曲线制作装置,其中,
所述列车的状态是针对每个所述规定的经过时间的离所述出发车站的距离以及所述列车的速度,
所述节能运转曲线制作部对于所述候补解,针对从所述出发车站起的规定的距离范围以及规定的速度范围的每个组合,将属于各组合的任一个的状态作为代表该组合的所述候补解,来制作所述节能运转曲线。
3.根据权利要求1或者2所述的运转曲线制作装置,其中,
通过档位切换来进行所述列车的加速以及减速,
所述节能运转曲线制作部将超过规定的档位切换数的所述列车的状态从所述候补解中排除。
4.根据权利要求1或者2所述的运转曲线制作装置,其中,
通过档位切换来进行所述列车的加速以及减速,
所述节能运转曲线制作部将与同一档位的持续时间小于规定的持续时间相当的所述列车的状态从所述候补解中排除。
5.根据权利要求1或者2所述的运转曲线制作装置,其中,
所述列车的车辆信息包含车辆的种类、车体重量、乘车率、编制车辆数、列车长度。
6.根据权利要求1或者2所述的运转曲线制作装置,其中,
所述地面信息包含:车站间距离信息、包含区间位置及与该区间位置相对应的上限速度的速度限制信息以及相对于行进方向的坡度信息。
7.一种运转曲线制作方法,对具有规定的车辆特性的列车制作能够在出发车站与到达车站之间以规定的行驶时间使累计消耗能量减少地行驶的运转曲线,该运转曲线制作方法具备:
最短时间运转曲线制作步骤,根据所述列车的车辆信息以及地面信息来制作所述列车在所述出发车站与所述到达车站之间以最短时间行驶的最短时间运转曲线;以及
节能运转曲线制作步骤,从所述最短时间运转曲线以及与从所述出发车站出发起针对每个规定的经过时间的所述列车的状态相当的用于定义列车的状态的多个候补解中选择在所述规定的行驶时间所述累计消耗能量少的候补解,根据所选择的候补解来制作节能运转曲线,
在所述节能运转曲线制作步骤中,从在所述最短时间运转曲线制作步骤中制作出的最短时间运转曲线的结果求出经过时间和位置的关系的上限值,从将所述最短时间运转曲线的结果在时间方向移动了规定行驶时间量的结果求出为经过时间和位置的关系的下限值,从该上限值和下限值中选择所述累计消耗能量少的候补解。
8.根据权利要求7所述的运转曲线制作方法,其中,
还具备获取步骤,在该获取步骤中获取所述地面信息、所述车辆信息以及规定的行驶时间,
所述最短时间运转曲线制作步骤具有:加速极限曲线制作步骤,根据包含使用最大加速的速度和位置的初始条件的信息制作表示使用了最大加速度时的速度和位置的加速度极限曲线;以及减速极限曲线制作步骤,搜索满足终端条件的条件,制作表示使用了最大减速度时的速度和位置的减速度极限曲线。
9.根据权利要求7所述的运转曲线制作方法,其中,
在所述节能运转曲线制作步骤中,
在所述出发车站与所述到达车站之间以规定行驶时间行驶时,将在该车站间以最短行驶时间Tmin行驶的情况下的经过时间t时的位置设为Upper[t]、将最短行驶时间下的(t,Upper[t])在时间轴方向平行移动了规定行驶时间与最短行驶时间Tmin的差分而得到的关系设为(t,Lower[t])的情况下,
如果将规定行驶时间中的能耗量成为最小的解的轨道设为(t,v(t),x(t)),则满足如下式子:
Lower[t]≤x(t)≤Upper[t]。
10.根据权利要求7所述的运转曲线制作方法,其中,
所述节能运转曲线制作步骤还具有量化步骤,在该量化步骤中通过将位置和速度以固定幅度进行分割来量化,将在该分割的区域中累计能耗量最小的值选择为候补解之一。
11.根据权利要求7所述的运转曲线制作方法,其中,
在解的条件中除了位置、速度之外,还包含档位切换次数、档位持续时间。
12.根据权利要求7所述的运转曲线制作方法,其中,
通过档位切换来进行所述列车的加速以及减速,
在所述节能运转曲线制作步骤中将超过规定的档位切换数的所述列车的状态从所述候补解中排除。
13.根据权利要求7所述的运转曲线制作方法,其中,
通过档位切换来进行所述列车的加速以及减速,
在所述节能运转曲线制作步骤中将与同一档位的持续时间小于规定的持续时间相当的所述列车的状态从所述候补解中排除。
14.一种计算机可读取的记录介质,记录有控制程序,所述控制程序用于通过计算机来控制运转曲线制作装置,该运转曲线制作装置对具有规定的车辆特性的列车制作能够在规定的出发车站与规定的到达车站之间以规定的行驶时间使累计消耗能量减少地行驶的运转曲线,所述控制程序使所述计算机作为如下单元发挥功能:
最短时间运转曲线制作单元,根据预先设定的所述列车的车辆信息以及地面信息来制作所述列车在所述车站间以最短时间行驶的最短时间运转曲线;以及
节能运转曲线制作单元,从所述最短时间运转曲线以及与从所述出发车站出发起针对每个规定的经过时间的所述列车的状态相当的用于定义列车的状态的多个候补解中选择在所述规定的行驶时间所述累计消耗能量少的候补解,根据所选择的候补解来制作节能运转曲线,
所述节能运转曲线制作单元从由所述最短时间运转曲线制作单元制作出的最短时间运转曲线的结果求出经过时间和位置的关系的上限值,从将所述最短时间运转曲线的结果在时间方向移动了规定行驶时间量的结果求出为经过时间和位置的关系的下限值,从该上限值和下限值中选择所述累计消耗能量少的候补解。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013-151104 | 2013-07-19 | ||
JP2013151104A JP6296716B2 (ja) | 2013-07-19 | 2013-07-19 | 運転曲線作成装置、運転曲線作成装置の制御方法及び制御プログラム |
PCT/JP2013/005289 WO2015008318A1 (ja) | 2013-07-19 | 2013-09-06 | 運転曲線作成装置、運転曲線作成方法及び運転曲線制御プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105452085A CN105452085A (zh) | 2016-03-30 |
CN105452085B true CN105452085B (zh) | 2018-01-02 |
Family
ID=52345814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380078342.9A Active CN105452085B (zh) | 2013-07-19 | 2013-09-06 | 运转曲线制作装置、运转曲线制作方法以及运转曲线控制程序 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10137911B2 (zh) |
EP (1) | EP3023314A4 (zh) |
JP (1) | JP6296716B2 (zh) |
CN (1) | CN105452085B (zh) |
WO (1) | WO2015008318A1 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8522134B1 (en) * | 2012-08-13 | 2013-08-27 | Volusion, Inc. | Methods and apparatus for in-line editing of web page content with reduced disruption of logical and presentational structure of content |
JP6324157B2 (ja) * | 2014-03-27 | 2018-05-16 | インターナショナル・ビジネス・マシーンズ・コーポレーションInternational Business Machines Corporation | 情報処理装置、情報処理方法、およびプログラム |
JP6404177B2 (ja) * | 2015-05-20 | 2018-10-10 | 公益財団法人鉄道総合技術研究所 | プログラム及び運転曲線作成装置 |
US10279823B2 (en) * | 2016-08-08 | 2019-05-07 | General Electric Company | System for controlling or monitoring a vehicle system along a route |
JP6707041B2 (ja) * | 2017-02-23 | 2020-06-10 | 株式会社日立製作所 | 車両制御システム、方法及び車両 |
CN111319657B (zh) * | 2018-12-13 | 2022-08-12 | 中车株洲电力机车研究所有限公司 | 城市轨道交通仿真系统 |
CN110941277B (zh) * | 2019-12-13 | 2023-02-17 | 华南智能机器人创新研究院 | 一种小车路线规划方法及系统 |
CN113968263B (zh) * | 2020-07-22 | 2022-12-09 | 比亚迪股份有限公司 | 一种列车运行曲线优化方法、装置及电子设备 |
IT202100007961A1 (it) * | 2021-03-31 | 2022-10-01 | Hitachi Rail S P A | Procedimento e sistema di controllo per generare un profilo di guida ottimale per veicoli provvisti di propulsione elettrica |
CN113635945B (zh) * | 2021-07-19 | 2022-09-09 | 交控科技股份有限公司 | 列车运行图调整方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1511744A (zh) * | 2002-01-31 | 2004-07-14 | 株式会社东芝 | 自动列车运行装置以及列车运行辅助装置 |
CN1772545A (zh) * | 2004-11-09 | 2006-05-17 | 株式会社日立制作所 | 列车行车时刻表评价装置及列车行车时刻表评价方法 |
CN101607559A (zh) * | 2009-06-18 | 2009-12-23 | 株洲南车时代电气股份有限公司 | 一种搜索列车牵引计算运行曲线的方法及装置 |
CN102887150A (zh) * | 2011-07-20 | 2013-01-23 | 株式会社日立制作所 | 列车控制系统 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69126644T2 (de) * | 1990-07-18 | 1997-12-18 | Hitachi Ltd | Verfahren zum Erzeugen eines Zuglaufplanes |
EP0539885B1 (en) * | 1991-10-25 | 1997-04-23 | Kabushiki Kaisha Toshiba | Optimal train running-pattern calculating apparatus and system including the same |
JPH07165080A (ja) | 1993-12-13 | 1995-06-27 | Hitachi Ltd | 移動体運行計画作成方法 |
JPH08156794A (ja) | 1994-10-05 | 1996-06-18 | Hitachi Ltd | 運転曲線作成方法及びその装置 |
AU2008201906B9 (en) * | 2002-05-20 | 2011-06-23 | TTG (Holdings) Pty Ltd | Method for improving timekeeping and saving energy on long-haul trains |
US9733625B2 (en) * | 2006-03-20 | 2017-08-15 | General Electric Company | Trip optimization system and method for a train |
JP4011539B2 (ja) * | 2003-12-15 | 2007-11-21 | 株式会社東芝 | 運転計画作成装置 |
US8370006B2 (en) * | 2006-03-20 | 2013-02-05 | General Electric Company | Method and apparatus for optimizing a train trip using signal information |
US9037323B2 (en) * | 2006-12-01 | 2015-05-19 | General Electric Company | Method and apparatus for limiting in-train forces of a railroad train |
US9233622B2 (en) * | 2008-03-11 | 2016-01-12 | General Electric Company | System and method for managing an amount of stored energy in a powered system |
JP5586308B2 (ja) * | 2010-04-01 | 2014-09-10 | 株式会社東芝 | 目標速度算出機能を備えた列車制御装置 |
DE102011106440B4 (de) * | 2011-07-04 | 2018-05-30 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Energieeffizienter Bahnbetrieb |
IN2014CN02754A (zh) * | 2011-10-19 | 2015-07-03 | Mitsubishi Electric Corp | |
US8838304B2 (en) * | 2012-06-29 | 2014-09-16 | Mitsubishi Electric Research Laboratories, Inc | Method for determining run-curves for vehicles based on travel time |
US8660723B2 (en) * | 2012-06-29 | 2014-02-25 | Mitsubishi Electric Research Laboratories, Inc. | Method for determining run-curves for vehicles in real-time subject to dynamic travel time and speed limit constraint |
US8983690B2 (en) * | 2013-03-14 | 2015-03-17 | Mitsubishi Electric Research Laboratories, Inc. | System and method for optimizing energy consumption in railway systems with energy storage devices |
US8831801B1 (en) * | 2013-03-14 | 2014-09-09 | Mitsubishi Electric Research Laboratories, Inc. | System and method for optimizing energy consumption in railway systems |
JP6334282B2 (ja) * | 2014-06-11 | 2018-05-30 | 株式会社東芝 | 情報処理装置および運転曲線作成方法 |
-
2013
- 2013-07-19 JP JP2013151104A patent/JP6296716B2/ja active Active
- 2013-09-06 WO PCT/JP2013/005289 patent/WO2015008318A1/ja active Application Filing
- 2013-09-06 CN CN201380078342.9A patent/CN105452085B/zh active Active
- 2013-09-06 EP EP13889480.3A patent/EP3023314A4/en not_active Withdrawn
-
2016
- 2016-01-13 US US14/994,904 patent/US10137911B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1511744A (zh) * | 2002-01-31 | 2004-07-14 | 株式会社东芝 | 自动列车运行装置以及列车运行辅助装置 |
CN1772545A (zh) * | 2004-11-09 | 2006-05-17 | 株式会社日立制作所 | 列车行车时刻表评价装置及列车行车时刻表评价方法 |
CN101607559A (zh) * | 2009-06-18 | 2009-12-23 | 株洲南车时代电气股份有限公司 | 一种搜索列车牵引计算运行曲线的方法及装置 |
CN102887150A (zh) * | 2011-07-20 | 2013-01-23 | 株式会社日立制作所 | 列车控制系统 |
Non-Patent Citations (1)
Title |
---|
Numerical Study on Dynamic Programming Applied to Optimization of Running Profile of a Train;Hideyoshi KO et al.;《Industry Applications Society》;20051231;1084-1092 * |
Also Published As
Publication number | Publication date |
---|---|
CN105452085A (zh) | 2016-03-30 |
JP6296716B2 (ja) | 2018-03-20 |
EP3023314A1 (en) | 2016-05-25 |
JP2015020637A (ja) | 2015-02-02 |
US10137911B2 (en) | 2018-11-27 |
US20160129926A1 (en) | 2016-05-12 |
EP3023314A4 (en) | 2017-04-12 |
WO2015008318A1 (ja) | 2015-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105452085B (zh) | 运转曲线制作装置、运转曲线制作方法以及运转曲线控制程序 | |
KR102611938B1 (ko) | 신경망을 사용한 통합 회로 플로어 플랜 생성 | |
JP6759673B2 (ja) | 衝突リスク算出プログラム、衝突リスク算出方法および衝突リスク算出装置 | |
Huang et al. | Fuzzy interpolation and extrapolation: A practical approach | |
Szu et al. | Fast simulated annealing | |
GB2592293A (en) | Operation of a vehicle using motion planning with machine learning | |
CN103630142B (zh) | 对地图数据进行处理以确定节能路线的方法和装置 | |
CN108153867A (zh) | 基于时间规律性的用户轨迹预测方法和装置 | |
Duz et al. | Real time estimation of local wave characteristics from ship motions using artificial neural networks | |
Chen et al. | Port yard storage optimization | |
CN113375689B (zh) | 导航方法、装置、终端和存储介质 | |
Yang et al. | Maritime traffic flow clustering analysis by density based trajectory clustering with noise | |
CN105808582A (zh) | 基于分层策略的决策树并行生成方法和装置 | |
Wang et al. | Personalized route recommendation with neural network enhanced search algorithm | |
CN113415288A (zh) | 分段式纵向车速规划方法、装置、设备及存储介质 | |
CN111967692A (zh) | 一种基于搜寻理论的搜救资源方案生成方法及相关设备 | |
Neumann | Method of path selection in the graph-case study | |
CN106503386A (zh) | 评估光功率预测算法性能优劣的方法及装置 | |
CN116036604B (zh) | 数据处理方法、装置、计算机及可读存储介质 | |
Prpić-Oršić et al. | Uncertainties of ship speed loss evaluation under real weather conditions | |
Sandidzadeh et al. | Genetic algorithm and particle swarm optimization algorithm for speed error reduction in railway signaling systems | |
Huang et al. | Processing continuous K-nearest skyline query with uncertainty in spatio-temporal databases | |
CN114399119A (zh) | 基于条件卷积生成式对抗网络的mmp预测方法及装置 | |
Dramski et al. | Selected shortest path in the graph algorithms with a use of trapezoidal grid | |
Li et al. | Safety-aware adversarial inverse reinforcement learning for highway autonomous driving |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |