CN111144722B - 一种基于伪旅行商问题的船舶引航排班方法与装置 - Google Patents
一种基于伪旅行商问题的船舶引航排班方法与装置 Download PDFInfo
- Publication number
- CN111144722B CN111144722B CN201911300847.0A CN201911300847A CN111144722B CN 111144722 B CN111144722 B CN 111144722B CN 201911300847 A CN201911300847 A CN 201911300847A CN 111144722 B CN111144722 B CN 111144722B
- Authority
- CN
- China
- Prior art keywords
- task
- piloting
- node
- time
- iwd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 72
- 238000013178 mathematical model Methods 0.000 claims abstract description 26
- 239000002689 soil Substances 0.000 claims description 60
- 239000013049 sediment Substances 0.000 claims description 21
- 239000004576 sand Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 238000005314 correlation function Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 9
- 230000007547 defect Effects 0.000 abstract description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000009024 positive feedback mechanism Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06316—Sequencing of tasks or work
-
- 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"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06311—Scheduling, planning or task assignment for a person or group
-
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于伪旅行商问题的船舶引航排班方法与装置,涉及一种船舶引航排班问题模型与基于伪旅行商问题的求解方法。方法依次包括以下步骤:根据问题描述,确定问题约束条件,建立船舶引航排班数学模型;将船舶引航排班问题数学模型转化为伪旅行商问题模型;采用改进智能水滴算法求解该伪旅行商问题模型,找出最优引航排班计划。本发明针对智能水滴算法启发性不足,容易陷入局部最优等缺点,提出了基于最优候选结点子群的选择策略,提高了搜索效率和算法精度,采用改进智能水滴算法,求解引航员指派模型,有效提高了求解结果的精度。
Description
技术领域
本发明属于船舶引航排班中的引航任务排班方法,具体涉及一种船舶引航排班问题模型与基于伪旅行商问题的求解方法。
背景技术
随着我国经济的飞速发展,世界经济日益全球化和一体化,港口作为各国间货物运输的主要通道显得更加重要。借助港口城市的先天优势,航运物流成为沿海港口城市经济发展的引擎,引航质量在一定程度上代表着港口的对外开放程度和市场竞争能力,同时引航排班能力也决定着引航服务效果和港口形象。只有提高调度能力,才能从根本上增强整体的引航质量。我国是一个港口大国,然而,在港口迅猛发展的同时,我国港口在引航排班效率上相对落后,这严重制约了港口的核心竞争力,也造成了较大的资源浪费,迫切需要建立一个高效可行的引航排班方法与系统实现引航资源的优化配置,提高港口竞争力。研究引航排班的建模优化和有效的解决方法,成为亟待解决的关键理论问题。
引航排班问题属于NP难的组合优化问题,非常适合采用智能优化算法进行求解。目前应用于组合优化的算法主要有局部搜索法,模拟退火算法,遗传算法,神经网络等。
智能水滴算法是一种新型群体智能算法,模拟自然界水流和泥沙相互作用形成水道的原理。智能水滴算法具有正反馈机制和较强的鲁棒性。已经成功应用在一些组合优化问题上,像旅行商问题、调度问题、车辆路径规划问题等,但是智能水滴算法存在启发性不足和容易陷入局部最优的缺点。
发明内容
本发明所要解决的技术问题是:
船舶引航排班过程中,现有的智能优化算法求解模型存在大量不可行解的问题。
为解决上述技术问题,本发明提供一种基于伪旅行商问题的船舶引航排班方法,包括以下步骤:
步骤(1a):建立本地数据库,进行数据预处理,具体步骤为:
1a1)读入待引航任务数量;各引航任务开始地点、到达地点;各引航任务开始时间、到达时间;
1a2)置等待时间阈值Twait、工作时间阈值Twork、晚班引航权系数与早班引航权系数;
步骤(1b):根据问题描述,确定问题约束条件,建立船舶引航排班数学模型;
步骤(1c):将船舶引航排班问题数学模型转化为伪旅行商问题模型;具体方法是将各船舶引航排班任务看作旅行商问题的各城市,将船舶引航排班任务的费用看作城市间的距离,将引航排班对各引航任务的访问看作是对旅行商问题城市的遍历,从而把船舶引航排班问题数学模型转化为旅行商问题模型;由于一个可行的引航排班计划由多个可行的引航任务组组成,且其费用不仅与组成引航任务组的各任务有关,且与各引航任务组的附加费用有关,故船舶引航排班问题数学模型转化所得旅行商问题模型为一伪旅行商问题模型;
步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,找出最优引航排班计划;
步骤(1e):根据确定的最优引航排班计划进行引航排班。
对于所述步骤(1b),根据问题描述,确定问题约束条件,建立船舶引航排班数学模型;具体步骤为:
2a)建立数学模型包括以下限制条件:
(2a1)在一个引航排班任务组中,前一个任务的到达地点必须等于下一个任务的出发地点;前一个任务的到达时间必须早于或等于下一个任务的出发时间,下一个任务的出发时间减去前一个任务的到达时间为等待时间,等待过程将会产生等待费用,各任务间的等待时间必须小于或等于等待时间阈值;
(2a2)一个任务组的最后一个任务的到达地点必须与第一个任务的出发地点相同,如任务组的出发地点与到达地点不同,将会产生过夜费用;
(2a3)一个引航排班任务组的总引航时间必须小于或等于工作时间阈值,以保证引航员得到休息,防止疲劳引航;
(2a4)一个引航排班任务只能属于一个引航排班任务组;
(2a5)一个引航任务组在各时间段引航的费用是不一样的,辛苦程度越高的班次费用越高,即早班(0:00-8:00)大于晚班(16:00-24:00),晚班大于白班(8:00-16:00);
(2b)根据以上条件,建立引航排班问题数学模型如下:
目标函数中加号连接五部分,第一部分为引航任务费用;第二部分为引航任务组数量费用,且组数越大,费用越高;第三部分为等待时间费用,表示一个可行的引航任务组组合中各引航任务组的时间等待费用之和;第四部分为过夜费用,表示一个可行的引航任务组组合中各引航任务组的外站过夜费用之和;第五部分为引航任务组工作时间一致性费用,表示一个可行的引航任务组组合中各引航任务组工作时间的标准差;w1,w2,w3,w4分别为权系数一、权系数二、权系数三、权系数四;
式中:costij表示引航完第i个任务后接着进行第j个任务服务所产生的费用,其值为第j个任务的引航时间与权系数的乘积,公式为:
costij=t1+q1t2+q2t3 (2)
其中,i,j为0均代表引航服务中心,t1、t2、t3分别代表第j个任务的引航时间处于白班、晚班、早班的时间段;q1、q2分别代表晚班权系数、早班权系数;
表示第k个引航任务组引航完第i个任务后接着进行第j个任务服务;
表示第i个引航任务属于第k个引航任务组
表示每一个引航任务只属于一个引航任务组;
其中M代表引航任务数,N代表引航任务组数,无需预先知道。
twaitk为第k个引航排班任务组的等待时间,由引航任务组中第二个任务至最后一个任务的等待时间之和求得,各任务等待时间由后一个任务的开始时刻减去前一个任务的到达时刻求得;
postk为第k个引航排班任务组的外站过夜费用,其公式为:
tworkk为第k个引航排班任务组的工作时间,tave为一个任务组平均工作时间,其计算公式如下:
对于所述步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,得到最优引航排班计划;具体步骤为:
3a)读入引航任务数量M,各引航任务的开始地点与到达地点,开始时间与到达时间;
3b)初始化静态参数:设置水滴数量NIWD;结点数量NC;泥沙常量InitSoil,速度变化量系数一av,速度变化量系数二bv,速度变化量系数三cv,泥沙变化量系数一as,泥沙变化量系数二bs,泥沙变化量系数三cs;局部泥沙量更新系数ρs;全局泥沙量更新系数ρn;最大迭代次数Nmax;初始化任意两个结点间的泥沙量soil(i,j);
3c)根据结点数量随机产生全局最优路线TB,并计算最优目标函数值Smin;
3d)初始化动态参数:设置每个水滴的速度velIWD,每个水滴携带的泥沙量soilIWD,每个水滴的访问列表Vc(IWD)={},迭代次数置为1;
3e)设置每个水滴的起始出发点,并将该点添加到访问列表Vc(IWD);
3f)采用基于最优候选结点子群的选择策略,对每个水滴计算各未访问结点的概率,并选择下一个待访问的结点;更新访问列表Vc(IWD);
3g)每个水滴从结点i到结点j后,更新水滴的速度velIWD;
其中,velIWD(t)表示水滴在t时刻的速度(或者理解为:位于结点i时的水滴的速度);
velIWD(t+1)表示水滴在t+1时刻的速度(或者理解为:位于结点j时的水滴的速度),soil(i,j)表示结点i到结点j路径中泥沙量;
3h)每个水滴从结点i到结点j后,计算泥沙变化量Δsoil(i,j);
其中,time(i,j;velIWD)表示水滴从结点i移动到结点j所需的时间:
其中,εv是正常数,一个很小的正常数,用于避免出现分母为零的情况;||c(i)-c(j)|| 为结点i与结点j之间的距离;c(i)表示第i个结点的坐标,c(j)表示第j个结点的坐标,||·||表示欧式距离;
3i)每个水滴从结点i到结点j后,更新水滴的泥沙量soilIWD;
soilIWD(t+1)=soilIWD(t)+Δsoil(i,j) (11)
其中,soilIWD(t)表示水滴在t时刻的携带的泥沙量(或者理解为:位于结点i时的水滴携带泥沙量),soilIWD(t+1)水滴在t+1时刻的携带的泥沙量(或者理解为:位于结点j时的水滴携带泥沙量);
3j)每个水滴从结点i到结点j后,更新结点i到结点j路径中泥沙量soil(i,j);
soil(i,j)=(1-ρs)·soil(i,j)-ρs·Δsoil(i,j) (12)
表示等式右边更新等式左边;
3k)对每个水滴重复步骤3f)~步骤3j)直至所有结点被访问,则本次迭代结束;根据目标函数计算所有水滴所对应的目标函数值,选择目标函数值最小的路线作为本次迭代中总成本最小的路线TIB;
3l)比较全局最优目标函数值Smin与步骤3k)中目标函数值的大小,如果步骤 3k)中的目标函数值更小,则更新Smin和全局最优路线TB;
3m)使用全局最优路径间的泥沙量信息更新路径中的泥沙量soil(i,j),
3n)判断是否达到最大迭代次数Nmax,如果是,输出全局最优路径TB以及全局最优目标函数值Smin;否则,迭代次数加一,返回步骤3e)。
所述步骤3b)初始化任意两个结点间的泥沙量soil(i,j)的具体方法为:采用初始泥沙量随机化机制,产生一个(0~1)间的随机数rn,并把随机数rn与泥沙常量InitSoil的乘积赋给soil(i,j),作为其初始值。
所述步骤(3c)与步骤(3k),根据目标函数对个体进行评估,采用的是基于伪旅行商问题的译码方法,具体步骤如下:
步骤(5a):设gtask=1,Nodtask=1,gtask=1代表当前引航任务组序号,Nodtask=1代表当前个体的结点号,置当前引航任务组的开始和到达时间、开始和到达地点分别为个体当前结点对应引航任务的开始和到达时间、开始和到达地点,置当前引航任务组的等待时间twait=0;
步骤(5b):将当前个体的结点号Nodtask对应引航任务放入gtask任务组;结点号Nodtask加1,如果结点号Nodtask>M,转步骤(5d);
步骤(5c):访问当前结点Nodtask,如果满足当前访问的引航任务的开始地点等于之前访问的引航任务的到达地点,当前访问的引航任务的开始时间不早于前一引航任务的到达时间且其差值小于或等于等待时间阈值Twait,从当前访问的引航任务组的第一个引航任务的开始时间到当前访问的引航任务的到达时间间隔小于工作时间阈值Twork,则把该任务放入gtask任务组,更新任务组等待时间twait,Nodtask加1,转步骤(5d);否则,任务组序号gtask加1,把当前任务放入更新的gtask任务组,结点号Nodtask加1;
步骤(5d):如果结点号Nodtask≤M,转步骤(5c);
步骤(5e):按目标函数计算引航排班问题的总费用,即为个体目标函数值。
所述步骤3f)采用基于最优候选结点子群的选择策略,选择下一个待访问的结点的具体过程为:
6a)读取所有未访问的结点,计算出当前结点i到所有候选结点V(IWD)之间的泥沙量,根据泥沙量由低到高,按照选择比例系数γ,选取一定比例的最优候选结点子群elite(IWD)作为候选节点,其公式如下所示:
elite(IWD)=γ·V(IWD) (14)
6b)根据概率公式Pi IWD(j)计算出每个候选结点的概率,其公式如下所示:
其中,k为候选结点,且为水滴未访问过的结点;f(soil(i,j))为结点i到结点j的路径泥沙量的相关函数,具体表示为如下公式:
其中,εs是正常数,一个很小的正数,用于避免出现分母为零的情况;g(soil(i,j)) 是将结点i到结点j之间的路径含沙量转换成正数的函数,具体表示为如下公
6c)根据各个候选节点的概率Pi IWD(j),使用轮盘赌策略选择下一个结点j;并用所述节点更新访问列表Vc(IWD)。
本发明的有益效果是:本发明在分析引航排班各种规则和习惯的基础上,提出了引航排班问题的数学模型,并提出了把该模型转化为伪旅行商问题模型的方法,较好地解决了采用智能优化算法直接求解模型存在大量不可行解的问题。提出了采用智能水滴算法求解该伪旅行商问题模型的方法。所提模型可行,所提算法搜索速度快,搜索精度高,稳定性强。
附图说明
图1是智能水滴算法的流程图;
图2是船舶引航排班方法译码流程图;
图3是船舶引航排班方法转化为伪旅行商问题模型示意图。
具体实施方式
现在结合附图和优选实施例对本发明作进一步的说明。
本发明的基于伪旅行商问题的船舶引航排班装方法,包括以下步骤:
步骤(1a):建立本地数据库,进行数据预处理;其具体步骤为:
1a1)读入待引航任务数量;各引航任务开始地点,到达地点;各引航任务开始时间、到达时间;
1a2)置等待时间阈值Twait、工作时间阈值Twork、晚班引航权系数与早班引航权系数。
步骤(1b):根据问题描述,确定问题约束条件,建立船舶引航排班数学模型;
步骤(1c):将船舶引航排班问题数学模型转化为伪旅行商问题模型;其具体方法是将各船舶引航排班任务看作旅行商问题的各城市,将船舶引航排班任务的费用看作城市间的距离,将引航排班对各引航任务的访问看作是对旅行商问题城市的遍历,从而把船舶引航排班问题数学模型转化为旅行商问题模型。由于一个可行的引航排班计划由多个可行的引航任务组组成,且其费用不仅与组成引航任务组的各任务有关,且与各引航任务组的附加费用有关,故船舶引航排班问题数学模型转化所得旅行商问题模型为一伪旅行商问题模型。
步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,找出最优引航排班计划;
步骤(1e):根据确定的最优引航排班计划进行引航排班。
对于所述步骤(1b),根据问题描述,确定问题约束条件,建立船舶引航排班数学模型;其具体步骤为:
2a)建立数学模型必须考虑以下限制条件:
(2a1)在一个引航排班任务组中,前一个任务的到达地点必须等于下一个任务的出发地点;前一个任务的到达时间必须早于等于下一个任务的出发时间,下一个任务的出发时间减去前一个任务的到达时间为等待时间,等待过程将会产生等待费用,各任务间的等待时间必须小于等于等待时间阈值;
(2a2)一个任务组的最后一个任务的到达地点必须与第一个任务的出发地点相同,如任务组的出发地点与到达地点不同,将会产生过夜费用;
(2a3)一个引航排班任务组的总引航时间必须小于等于工作时间阈值,以保证引航员得到休息,防止疲劳引航;
(2a4)一个引航排班任务只能属于一个引航排班任务组;
(2a5)一个引航任务组在各时间段引航的费用是不一样的,越辛苦的班次费用越高,即早班(0:00-8:00)大于晚班(16:00-24:00),晚班大于白班(8: 00-16:00)。
(2b)根据以上条件,建立引航排班问题数学模型如下:
目标函数由五部分组成,第一部分为引航任务费用;第二部分为引航任务组数量费用,且组数越大,费用越高;第三部分为等待时间费用,表示一个可行的引航任务组组合中各引航任务组的时间等待费用之和;第四部分为过夜费用,表示一个可行的引航任务组组合中各引航任务组的外站过夜费用之和;第五部分为引航任务组工作时间一致性费用,表示一个可行的引航任务组组合中各引航任务组工作时间的标准差。w1,w2,w3,w4为权系数。
式中:costij表示引航完第i个任务后接着进行第j个任务服务所产生的费用,其值为第j个任务的引航时间与权系数的乘积,公式为:
costij=t1+q1t2+q2t3 (2)
其中,i,j为0代表引航服务中心,t1、t2、t3分别代表第j个任务的引航时间处于白班、晚班、早班的时间段;q1、q2分别代表晚班、早班权系数。
表示第k个引航任务组引航完第i个任务后接着进行第j个任务服务;
表示第i个引航任务属于第k个引航任务组
表示每一个引航任务只属于一个引航任务组;
其中M代表引航任务数,N代表引航任务组数,无需预先知道。
twaitk为第k个引航排班任务组的等待时间,由引航任务组中第二个任务至最后一个任务的等待时间之和求得,各任务等待时间由后一个任务的开始时刻减去前一个任务的到达时刻求得;
postk为第k个引航排班任务组的外站过夜费用,其公式为:
tworkk为第k个引航排班任务组的工作时间,tave为一个任务组平均工作时间,其计算公式如下:
对于所述步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,找出最优引航排班计划;其具体步骤为:
3a)读入引航任务数量M,各任务的开始地点与到达地点,开始时间与到达时间;
3b)初始化静态参数:设置水滴数量NIWD,一般NIWD取40-100;结点数量NC,一般NC取预约老人的数量M;泥沙常量InitSoil=1000,速度变化量系数av=1000, bv=0.01,cv=1;泥沙变化量系数as=1000,bs=0.01,cs=1;局部泥沙量更新系数ρs=0.9;全局泥沙量更新系数ρn=0.9;最大迭代次数Nmax=1000;初始化任意两个结点间的泥沙量soil(i,j);3c)根据结点数量随机产生全局最优路线TB,并计算最优目标函数值Smin;
3d)初始化动态参数:设置每个水滴的速度velIWD=200,每个水滴携带的泥沙量soilIWD=0,每个水滴的访问列表Vc(IWD)={},迭代次数置为一;
3e)设置每个水滴的起始出发点,并将该点添加到访问列表Vc(IWD);
3f)采用基于最优候选结点子群的选择策略,对每个水滴计算各未访问结点的概率,并选择下一个待访问的结点;更新访问列表Vc(IWD);
3g)每个水滴从结点i到结点j后,更新水滴的速度velIWD;
其中,velIWD(t)表示水滴在t时刻的速度(或者理解为:位于结点i时的水滴的速度),velIWD(t+1)表示水滴在t+1时刻的速度(或者理解为:位于结点j时的水滴的速度),soil(i,j)表示结点i到结点j路径中泥沙量;
3h)每个水滴从结点i到结点j后,计算泥沙变化量Δsoil(i,j);
其中,time(i,j;velIWD)表示水滴从结点i移动到结点j所需的时间:
其中,εv是正常数,一个很小的正常数,用于避免出现分母为零的情况;||c(i)-c(j)|| 为结点i与结点j之间的距离;c(i)表示第i个结点的坐标,c(j)表示第j个结点的坐标,||·||表示欧式距离;
3i)每个水滴从结点i到结点j后,更新水滴的泥沙量soilIWD;
soilIWD(t+1)=soilIWD(t)+Δsoil(i,j) (11)
其中,soilIWD(t)表示水滴在t时刻的携带的泥沙量(或者理解为:位于结点i时的水滴携带泥沙量),soilIWD(t+1)水滴在t+1时刻的携带的泥沙量(或者理解为:位于结点j时的水滴携带泥沙量);
3j)每个水滴从结点i到结点j后,更新结点i到结点j路径中泥沙量soil(i,j);
soil(i,j)=(1-ρs)·soil(i,j)-ρs·Δsoil(i,j) (12)
3k)对每个水滴重复步骤3f)~步骤3j)直至所有结点被访问,则本次迭代结束;
根据目标函数计算所有水滴所对应的目标函数值,选择目标函数值最小的路线作为本次迭代中总成本最小的路线TIB;
3l)比较全局最优目标函数值Smin与步骤3k)中目标函数值的大小,如果步骤 3k)中的目标函数值更小,则更新Smin和全局最优路线TB;
3m)使用全局最优路径间的泥沙量信息更新路径中的泥沙量soil(i,j),
3n)判断是否达到最大迭代次数Nmax,如果是,那么输出全局最优路径TB以及全局最优目标函数值Smin;否则,迭代次数加一,返回步骤3e)。
所述步骤3b)初始化任意两个结点间的泥沙量soil(i,j)的具体方法为:采用初始泥沙量随机化机制,产生一个(0~1)间的随机数rn,并把rn与泥沙常量InitSoil 的乘积赋给soil(i,j),作为其初始值。
对于所述步骤(3c)与(3k),根据目标函数对个体进行评估,采用的是基于伪旅行商问题的译码方法,其具体步骤如下:
步骤(5a):设gtask=1,Nodtask=1。gtask=1代表当前引航任务组序号,Nodtask=1代表当前个体的结点号,置当前引航任务组的开始和到达时间、开始和到达地点分别为个体当前结点对应引航任务的开始和到达时间、开始和到达地点,置当前引航任务组的等待时间twait=0;
步骤(5b):将当前个体的结点号Nodtask对应引航任务放入gtask任务组;Nodtask加1,如果Nodtask>M,转步骤(5d);
步骤(5c):访问当前结点Nodtask,如果满足当前访问的引航任务的开始地点等于之前访问的引航任务的到达地点,当前访问的引航任务的开始时间不早于前一引航任务的到达时间且其差值小于等于等待时间阈值Twait,从当前访问的引航任务组的第一个引航任务的开始时间到当前访问的引航任务的到达时间间隔小于工作时间阈值Twork,则把该任务放入gtask任务组,更新任务组等待时间twait,Nodtask加1,转步骤(5d);否则,gtask加1,把当前任务放入更新的gtask任务组,Nodtask加1;
步骤(5d):如果Nodtask≤M,转步骤(5c);
步骤(5e):按目标函数计算引航排班问题的总费用,即为个体目标函数值。
参见图3,设一个个体共有10个结点,1~10是结点编号,其中,箭头代表引航任务组在服务完前一结点对应引航任务后可接着服务下一结点的对应引航任务,虚线代表在服务完前一结点对应引航任务后不可接着服务下一结点的对应引航任务,需要更换引航任务组。根据如上译码方式,可确定本个体10个结点对应引航任务需要组成4个引航任务组,并得到各引航任务组的组成列表。其中第一个引航任务组的组成结点列表为[1],该引航任务组由结点号为1对应的引航任务单独组成;第二个引航任务组的组成结点列表为[2,3,4,5],依次服务结点号为2,3,4,5对应的引航任务;第三个引航任务组的组成结点列表为[6,7],依次服务结点号为6,7对应的引航任务;第四个引航任务组的组成结点列表为 [8,9,10],依次服务结点号为8,9,10对应的引航任务。
所述步骤3f)采用基于最优候选结点子群的选择策略,选择下一个待访问的结点的具体过程为:
6a)读取所有未访问的结点,计算出当前结点i到所有候选结点V(IWD)之间的泥沙量,根据泥沙量由低到高,按照选择比例系数γ,选取一定比例的最优候选结点子群elite(IWD)作为候选节点,其公式如下所示:
elite(IWD)=γ·V(IWD) (14)
6b)根据概率公式Pi IWD(j)计算出每个候选结点的概率,其公式如下所示:
其中,k为候选结点,且为水滴未访问过的结点;f(soil(i,j))为结点i到结点j的路径泥沙量的相关函数,具体表示为如下公式:
其中,εs是正常数,一个很小的正数,用于避免出现分母为零的情况;g(soil(i,j)) 是将结点i到结点j之间的路径含沙量转换成正数的函数,具体表示为如下公式:
6c)根据各个候选节点的概率Pi IWD(j),使用轮盘赌策略选择下一个结点j。并用该节点更新访问列表Vc(IWD)。
上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人士能够了解本发明的内容并加以实施,并不能以此限制本发明的保护范围,凡根据本发明精神实质所作的等效变化或修饰,都应涵盖在本发明的保护范围内。
Claims (4)
1.一种基于伪旅行商问题的船舶引航排班方法,其特征在于:包括以下步骤:
步骤(1a):建立本地数据库,进行数据预处理,具体步骤为:
1a1)读入待引航任务数量;各引航任务开始地点,到达地点;各引航任务开始时间、到达时间;
1a2)置等待时间阈值Twait、工作时间阈值Twork、晚班引航权系数与早班引航权系数;
步骤(1b):根据问题描述,确定问题约束条件,建立船舶引航排班数学模型,具体步骤为:
2a)建立数学模型必须考虑以下限制条件:
(2a1)在一个引航排班任务组中,前一个任务的到达地点必须等于下一个任务的出发地点;前一个任务的到达时间必须早于等于下一个任务的出发时间,下一个任务的出发时间减去前一个任务的到达时间为等待时间,等待过程将会产生等待费用,各任务间的等待时间必须小于等于等待时间阈值;
(2a2)一个任务组的最后一个任务的到达地点必须与第一个任务的出发地点相同,如任务组的出发地点与到达地点不同,将会产生过夜费用;
(2a3)一个引航排班任务组的总引航时间必须小于等于工作时间阈值,以保证引航员得到休息,防止疲劳引航;
(2a4)一个引航排班任务只能属于一个引航排班任务组;
(2a5)一个引航任务组在各时间段引航的费用是不一样的,越辛苦的班次费用越高;
(2b)根据以上条件,建立引航排班问题数学模型如下:
目标函数由五部分组成,第一部分为引航任务费用;第二部分为引航任务组数量费用,且组数越大,费用越高;第三部分为等待时间费用,表示一个可行的引航任务组组合中各引航任务组的时间等待费用之和;第四部分为过夜费用,表示一个可行的引航任务组组合中各引航任务组的外站过夜费用之和;第五部分为引航任务组工作时间一致性费用,表示一个可行的引航任务组组合中各引航任务组工作时间的标准差;w1,w2,w3,w4分别为权系数一、权系数二、权系数三、权系数四;
式中:costij表示引航完第i个任务后接着进行第j个任务服务所产生的费用,其值为第j个任务的引航时间与权系数的乘积,公式为:
costij=t1+q1t2+q2t3 (2)
其中,i,j为0均代表引航服务中心,t1、t2、t3分别代表第j个任务的引航时间处于白班、晚班、早班的时间段;q1、q2分别代表晚班、早班权系数;
表示第k个引航任务组引航完第i个任务后接着进行第j个任务服务;
表示第i个引航任务属于第k个引航任务组
表示每一个引航任务只属于一个引航任务组;
其中M代表引航任务数,N代表引航任务组数;
twaitk为第k个引航排班任务组的等待时间,由引航任务组中第二个任务至最后一个任务的等待时间之和求得,各任务等待时间由后一个任务的开始时刻减去前一个任务的到达时刻求得;
postk为第k个引航排班任务组的外站过夜费用,其公式为:
tworkk为第k个引航排班任务组的工作时间,tave为一个任务组平均工作时间,其计算公式如下:
步骤(1c):将船舶引航排班问题数学模型转化为伪旅行商问题模型;其具体方法是将各船舶引航排班任务看作旅行商问题的各城市,将船舶引航排班任务的费用看作城市间的距离,将引航排班对各引航任务的访问看作是对旅行商问题城市的遍历,从而把船舶引航排班问题数学模型转化为旅行商问题模型,由于一个可行的引航排班计划由多个可行的引航任务组组成,且其费用不仅与组成引航任务组的各任务有关,且与各引航任务组的附加费用有关,故船舶引航排班问题数学模型转化所得旅行商问题模型为一伪旅行商问题模型;
步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,找出最优引航排班计划,具体步骤为:
3a)读入引航任务数量M,各任务的开始地点与到达地点,开始时间与到达时间;
3b)初始化静态参数:设置水滴数量NIWD;结点数量NC;泥沙常量InitSoil,速度变化量系数av,bv,cv,泥沙变化量系数as,bs,cs;局部泥沙量更新系数ρs;全局泥沙量更新系数ρn;最大迭代次数Nmax;初始化任意两个结点间的泥沙量soil(i,j);
3c)根据结点数量随机产生全局最优路线TB,并计算最优目标函数值Smin;
3d)初始化动态参数:设置每个水滴的速度velIWD,每个水滴携带的泥沙量soilIWD,每个水滴的访问列表Vc(IWD)={},迭代次数置为一;
3e)设置每个水滴的起始出发点,并将该点添加到访问列表Vc(IWD);
3f)采用基于最优候选结点子群的选择策略,对每个水滴计算各未访问结点的概率,并选择下一个待访问的结点;更新访问列表Vc(IWD);
3g)每个水滴从结点i到结点j后,更新水滴的速度velIWD;
其中,velIWD(t)表示水滴在t时刻的速度,velIWD(t+1)表示水滴在t+1时刻的速度,soil(i,j)表示结点i到结点j路径中泥沙量;
3h)每个水滴从结点i到结点j后,计算泥沙变化量Δsoil(i,j);
其中,time(i,j;velIWD)表示水滴从结点i移动到结点j所需的时间:
其中,εv是正常数,用于避免出现分母为零的情况;||c(i)-c(j)||为结点i与结点j之间的距离;c(i)表示第i个结点的坐标,c(j)表示第j个结点的坐标,||·||表示欧式距离;
3i)每个水滴从结点i到结点j后,更新水滴的泥沙量soilIWD;
soilIWD(t+1)=soilIWD(t)+Δsoil(i,j) (11)
其中,soilIWD(t)表示水滴在t时刻的携带的泥沙量,soilIWD(t+1)水滴在t+1时刻的携带的泥沙量;
3j)每个水滴从结点i到结点j后,更新结点i到结点j路径中泥沙量soil(i,j);
soil(i,j)=(1-ρs)·soil(i,j)-ρs·Δsoil(i,j) (12)
表示等式右边更新等式左边;
3k)对每个水滴重复步骤3f)~步骤3j)直至所有结点被访问,则本次迭代结束;根据目标函数计算所有水滴所对应的目标函数值,选择目标函数值最小的路线作为本次迭代中总成本最小的路线TIB;
3l)比较全局最优目标函数值Smin与步骤3k)中目标函数值的大小,如果步骤3k)中的目标函数值更小,则更新Smin和全局最优路线TB;
3m)使用全局最优路径间的泥沙量信息更新路径中的泥沙量soil(i,j),
3n)判断是否达到最大迭代次数Nmax,如果是,那么输出全局最优路径TB以及全局最优目标函数值Smin;否则,迭代次数加一,返回步骤3e);
对于所述步骤(3c)与(3k),根据目标函数对个体进行评估,采用基于伪旅行商问题的译码方法,其具体步骤如下:
步骤(5a):设gtask=1,Nodtask=1,gtask=1代表当前引航任务组序号,Nodtask=1代表当前个体的结点号,置当前引航任务组的开始和到达时间、开始和到达地点分别为个体当前结点对应引航任务的开始和到达时间、开始和到达地点,置当前引航任务组的等待时间twait=0;
步骤(5b):将当前个体的结点号Nodtask对应引航任务放入gtask任务组;Nodtask加1,如果Nodtask>M,转步骤(5d);
步骤(5c):访问当前结点Nodtask,如果满足当前访问的引航任务的开始地点等于之前访问的引航任务的到达地点,当前访问的引航任务的开始时间不早于前一引航任务的到达时间且其差值小于等于等待时间阈值Twait,从当前访问的引航任务组的第一个引航任务的开始时间到当前访问的引航任务的到达时间间隔小于工作时间阈值Twork,则把该任务放入gtask任务组,更新任务组等待时间twait,Nodtask加1,转步骤(5d);否则,gtask加1,把当前任务放入更新的gtask任务组,Nodtask加1;
步骤(5d):如果Nodtask≤M,转步骤(5c);
步骤(5e):按目标函数计算引航排班问题的总费用,即为个体目标函数值;
步骤(1e):根据确定的最优引航排班计划进行引航排班。
2.根据权利要求1所述的基于伪旅行商问题的船舶引航排班方法,其特征在于,其特征在于:所述步骤3b)初始化任意两个结点间的泥沙量soil(i,j)的具体方法为:采用初始泥沙量随机化机制,产生一个(0~1)间的随机数rn,并把rn与泥沙常量InitSoil的乘积赋给soil(i,j),作为其初始值。
3.根据权利要求1所述的基于伪旅行商问题的船舶引航排班方法,其特征在于,所述步骤3f)采用基于最优候选结点子群的选择策略,选择下一个待访问的结点的具体过程为:
6a)读取所有未访问的结点,计算出当前结点i到所有候选结点V(IWD)之间的泥沙量,根据泥沙量由低到高,按照选择比例系数γ,选取一定比例的最优候选结点子群elite(IWD)作为候选节点,其公式如下所示:
elite(IWD)=γ·V(IWD) (14)
6b)根据概率公式Pi IWD(j)计算出每个候选结点的概率,其公式如下所示:
其中,k为候选结点,且为水滴未访问过的结点;f(soil(i,j))为结点i到结点j的路径泥沙量的相关函数,具体表示为如下公式:
其中,εs是正常数,用于避免出现分母为零的情况;g(soil(i,j))是将结点i到结点j之间的路径含沙量转换成正数的函数,具体表示为如下公式:
6c)根据各个候选节点的概率Pi IWD(j),使用轮盘赌策略选择下一个结点j,并用所述节点更新访问列表Vc(IWD)。
4.一种基于伪旅行商问题的船舶引航排班装置,包括本地数据库,其特征在于:本地数据库进行数据处理的过程包括以下步骤:
步骤(1a):建立本地数据库,进行数据预处理,具体步骤为:
1a1)读入待引航任务数量;各引航任务开始地点,到达地点;各引航任务开始时间、到达时间;
1a2)置等待时间阈值Twait、工作时间阈值Twork、晚班引航权系数与早班引航权系数;
步骤(1b):根据问题描述,确定问题约束条件,建立船舶引航排班数学模型,
具体步骤为:
2a)建立数学模型必须考虑以下限制条件:
(2a1)在一个引航排班任务组中,前一个任务的到达地点必须等于下一个任务的出发地点;前一个任务的到达时间必须早于等于下一个任务的出发时间,下一个任务的出发时间减去前一个任务的到达时间为等待时间,等待过程将会产生等待费用,各任务间的等待时间必须小于等于等待时间阈值;
(2a2)一个任务组的最后一个任务的到达地点必须与第一个任务的出发地点相同,如任务组的出发地点与到达地点不同,将会产生过夜费用;
(2a3)一个引航排班任务组的总引航时间必须小于等于工作时间阈值,以保证引航员得到休息,防止疲劳引航;
(2a4)一个引航排班任务只能属于一个引航排班任务组;
(2a5)一个引航任务组在各时间段引航的费用是不一样的,越辛苦的班次费用越高;
(2b)根据以上条件,建立引航排班问题数学模型如下:
目标函数由五部分组成,第一部分为引航任务费用;第二部分为引航任务组数量费用,且组数越大,费用越高;第三部分为等待时间费用,表示一个可行的引航任务组组合中各引航任务组的时间等待费用之和;第四部分为过夜费用,表示一个可行的引航任务组组合中各引航任务组的外站过夜费用之和;第五部分为引航任务组工作时间一致性费用,表示一个可行的引航任务组组合中各引航任务组工作时间的标准差;w1,w2,w3,w4分别为权系数一、权系数二、权系数三、权系数四;
式中:costij表示引航完第i个任务后接着进行第j个任务服务所产生的费用,其值为第j个任务的引航时间与权系数的乘积,公式为:
costij=t1+q1t2+q2t3 (2)
其中,i,j为0均代表引航服务中心,t1、t2、t3分别代表第j个任务的引航时间处于白班、晚班、早班的时间段;q1、q2分别代表晚班、早班权系数;
表示第k个引航任务组引航完第i个任务后接着进行第j个任务服务;
表示第i个引航任务属于第k个引航任务组
表示每一个引航任务只属于一个引航任务组;
其中M代表引航任务数,N代表引航任务组数;
twaitk为第k个引航排班任务组的等待时间,由引航任务组中第二个任务至最后一个任务的等待时间之和求得,各任务等待时间由后一个任务的开始时刻减去前一个任务的到达时刻求得;
postk为第k个引航排班任务组的外站过夜费用,其公式为:
tworkk为第k个引航排班任务组的工作时间,tave为一个任务组平均工作时间,其计算公式如下:
步骤(1c):将船舶引航排班问题数学模型转化为伪旅行商问题模型;其具体方法是将各船舶引航排班任务看作旅行商问题的各城市,将船舶引航排班任务的费用看作城市间的距离,将引航排班对各引航任务的访问看作是对旅行商问题城市的遍历,从而把船舶引航排班问题数学模型转化为旅行商问题模型,由于一个可行的引航排班计划由多个可行的引航任务组组成,且其费用不仅与组成引航任务组的各任务有关,且与各引航任务组的附加费用有关,故船舶引航排班问题数学模型转化所得旅行商问题模型为一伪旅行商问题模型;
步骤(1d):基于伪旅行商问题对步骤(1c)所得伪旅行商问题模型进行求解,找出最优引航排班计划,具体步骤为:
3a)读入引航任务数量M,各任务的开始地点与到达地点,开始时间与到达时间;
3b)初始化静态参数:设置水滴数量NIWD;结点数量NC;泥沙常量InitSoil,速度变化量系数av,bv,cv,泥沙变化量系数as,bs,cs;局部泥沙量更新系数ρs;全局泥沙量更新系数ρn;最大迭代次数Nmax;初始化任意两个结点间的泥沙量soil(i,j);
3c)根据结点数量随机产生全局最优路线TB,并计算最优目标函数值Smin;
3d)初始化动态参数:设置每个水滴的速度velIWD,每个水滴携带的泥沙量soilIWD,每个水滴的访问列表Vc(IWD)={},迭代次数置为一;
3e)设置每个水滴的起始出发点,并将该点添加到访问列表Vc(IWD);
3f)采用基于最优候选结点子群的选择策略,对每个水滴计算各未访问结点的概率,并选择下一个待访问的结点;更新访问列表Vc(IWD);
3g)每个水滴从结点i到结点j后,更新水滴的速度velIWD;
其中,velIWD(t)表示水滴在t时刻的速度,velIWD(t+1)表示水滴在t+1时刻的速度,soil(i,j)表示结点i到结点j路径中泥沙量;
3h)每个水滴从结点i到结点j后,计算泥沙变化量Δsoil(i,j);
其中,time(i,j;velIWD)表示水滴从结点i移动到结点j所需的时间:
其中,εv是正常数,用于避免出现分母为零的情况;||c(i)-c(j)||为结点i与结点j之间的距离;c(i)表示第i个结点的坐标,c(j)表示第j个结点的坐标,||·||表示欧式距离;
3i)每个水滴从结点i到结点j后,更新水滴的泥沙量soilIWD;
soilIWD(t+1)=soilIWD(t)+Δsoil(i,j) (11)
其中,soilIWD(t)表示水滴在t时刻的携带的泥沙量,soilIWD(t+1)水滴在t+1时刻的携带的泥沙量;
3j)每个水滴从结点i到结点j后,更新结点i到结点j路径中泥沙量soil(i,j);
soil(i,j)=(1-ρs)·soil(i,j)-ρs·Δsoil(i,j) (12)
表示等式右边更新等式左边;
3k)对每个水滴重复步骤3f)~步骤3j)直至所有结点被访问,则本次迭代结束;根据目标函数计算所有水滴所对应的目标函数值,选择目标函数值最小的路线作为本次迭代中总成本最小的路线TIB;
3l)比较全局最优目标函数值Smin与步骤3k)中目标函数值的大小,如果步骤3k)中的目标函数值更小,则更新Smin和全局最优路线TB;
3m)使用全局最优路径间的泥沙量信息更新路径中的泥沙量soil(i,j),
3n)判断是否达到最大迭代次数Nmax,如果是,那么输出全局最优路径TB以及全局最优目标函数值Smin;否则,迭代次数加一,返回步骤3e);
对于所述步骤(3c)与(3k),根据目标函数对个体进行评估,采用基于伪旅行商问题的译码方法,其具体步骤如下:
步骤(5a):设gtask=1,Nodtask=1,gtask=1代表当前引航任务组序号,Nodtask=1代表当前个体的结点号,置当前引航任务组的开始和到达时间、开始和到达地点分别为个体当前结点对应引航任务的开始和到达时间、开始和到达地点,置当前引航任务组的等待时间twait=0;
步骤(5b):将当前个体的结点号Nodtask对应引航任务放入gtask任务组;Nodtask加1,如果Nodtask>M,转步骤(5d);
步骤(5c):访问当前结点Nodtask,如果满足当前访问的引航任务的开始地点等于之前访问的引航任务的到达地点,当前访问的引航任务的开始时间不早于前一引航任务的到达时间且其差值小于等于等待时间阈值Twait,从当前访问的引航任务组的第一个引航任务的开始时间到当前访问的引航任务的到达时间间隔小于工作时间阈值Twork,则把该任务放入gtask任务组,更新任务组等待时间twait,Nodtask加1,转步骤(5d);否则,gtask加1,把当前任务放入更新的gtask任务组,Nodtask加1;
步骤(5d):如果Nodtask≤M,转步骤(5c);
步骤(5e):按目标函数计算引航排班问题的总费用,即为个体目标函数值;
步骤(1e):根据确定的最优引航排班计划进行引航排班。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911300847.0A CN111144722B (zh) | 2019-12-17 | 2019-12-17 | 一种基于伪旅行商问题的船舶引航排班方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911300847.0A CN111144722B (zh) | 2019-12-17 | 2019-12-17 | 一种基于伪旅行商问题的船舶引航排班方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111144722A CN111144722A (zh) | 2020-05-12 |
CN111144722B true CN111144722B (zh) | 2022-08-16 |
Family
ID=70518537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911300847.0A Active CN111144722B (zh) | 2019-12-17 | 2019-12-17 | 一种基于伪旅行商问题的船舶引航排班方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111144722B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217118A (zh) * | 2014-09-12 | 2014-12-17 | 河海大学常州校区 | 一种船舶引航排班问题模型与求解方法 |
CN104239975A (zh) * | 2014-09-12 | 2014-12-24 | 河海大学常州校区 | 基于改进离散粒子群算法的船舶引航排班方法 |
CN109523151A (zh) * | 2018-11-09 | 2019-03-26 | 上海港引航站 | 一种基于数据驱动的港口船舶引航智能调度方法 |
-
2019
- 2019-12-17 CN CN201911300847.0A patent/CN111144722B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217118A (zh) * | 2014-09-12 | 2014-12-17 | 河海大学常州校区 | 一种船舶引航排班问题模型与求解方法 |
CN104239975A (zh) * | 2014-09-12 | 2014-12-24 | 河海大学常州校区 | 基于改进离散粒子群算法的船舶引航排班方法 |
CN109523151A (zh) * | 2018-11-09 | 2019-03-26 | 上海港引航站 | 一种基于数据驱动的港口船舶引航智能调度方法 |
Non-Patent Citations (1)
Title |
---|
基于改进智能水滴算法的混合时间窗电商物流路径优化研究;王奕璇等;《科技管理研究》;20181231(第11期);第211-218页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111144722A (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020176381A1 (en) | Joint order dispatching and fleet management for online ride-sharing platforms | |
CN106225800B (zh) | 基于实时路况信息的环境友好型车辆导航路径构建方法 | |
CN105841709B (zh) | 一种汽车行驶路径规划方法 | |
CN111722643A (zh) | 仿狼群合作狩猎机制的无人机集群动态任务分配方法 | |
CN113065074B (zh) | 一种基于知识图谱和自注意力机制的轨迹目的地预测的方法 | |
CN107220731A (zh) | 一种物流配送路径规划方法 | |
CN105260785B (zh) | 一种基于改进布谷鸟算法的物流配送车辆路径优化方法 | |
JP2021103566A (ja) | 交通システム内でのサイト及び場所の間のルート又はルート距離を使用してサイトを検索又は比較するための方法及び装置 | |
CN112013829A (zh) | 基于多目标优化的多uav/ugv协同长时作业路径规划方法 | |
CN115713856A (zh) | 一种基于交通流预测与实际路况的车辆路径规划方法 | |
CN112836974B (zh) | 一种基于dqn和mcts的箱区间多场桥动态调度方法 | |
CN106991558A (zh) | 一种港口口岸间主航道的自动生成方法及系统 | |
Li et al. | Dynamic route and departure time choice model based on self-adaptive reference point and reinforcement learning | |
Starkey et al. | iPatch: A many-objective type-2 fuzzy logic system for field workforce optimization | |
CN113240339B (zh) | 一种面向大规模打车平台的任务匹配公平方法 | |
CN111090935A (zh) | 一种公共自行车预约调度与路径规划方法 | |
CN111144722B (zh) | 一种基于伪旅行商问题的船舶引航排班方法与装置 | |
CN110363358A (zh) | 基于多智能体模拟的公共交通分担率预测方法 | |
CN109974708A (zh) | 一种无人船航迹规划方法、终端设备及存储介质 | |
Zhang et al. | APFD: an effective approach to taxi route recommendation with mobile trajectory big data | |
CN116659535A (zh) | 一种基于深度强化学习的风险感知路径规划方法 | |
CN115729106A (zh) | 一种多时间约束下行驶里程优化的轻量化拼车调度方法 | |
CN115841286A (zh) | 一种基于深度强化学习的外卖配送路径规划方法 | |
CN111008788B (zh) | 一种基于改进离散头脑风暴优化算法的引航员指派方法 | |
CN113326980A (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 |