CN111091286B - 一种公共自行车调度方法 - Google Patents
一种公共自行车调度方法 Download PDFInfo
- Publication number
- CN111091286B CN111091286B CN201911289386.1A CN201911289386A CN111091286B CN 111091286 B CN111091286 B CN 111091286B CN 201911289386 A CN201911289386 A CN 201911289386A CN 111091286 B CN111091286 B CN 111091286B
- Authority
- CN
- China
- Prior art keywords
- service
- transport vehicle
- station
- individual
- time
- 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 36
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 50
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 230000035772 mutation Effects 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 230000000474 nursing effect Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims 1
- 230000007723 transport mechanism Effects 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012549 training Methods 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- 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/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种公共自行车调度方法，包括以下步骤：获取公共自行车本地数据，并进行预处理；根据预处理后的数据预先构建公共自行车调度模型；将公共自行车调度模型转化为伪旅行商模型；采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解；输出最优公共自行车调度模型的目标函数值，调度路径规划结果。优点：本发明有效的消除了常用头脑风暴优化算法产生的不可行解问题；能够合理的分配运输车资源，能够满足城市公共交通机构对运输车资源合理的调度，达到对运输车服务路径的优化。
Description
技术领域
本发明涉及一种公共自行车调度方法，属于城市智能公共交通系统技术领域。
背景技术
随着城市化进程的不断加快，城市交通问题日益严峻，严重阻碍了城市的发展。公共自行车因其低碳环保、方便灵活的特点以及独特的短距离出行优势,在各大城市快速发展起来。由于公共自行车可与公共交通实现“无缝对接”，所以“公交+公共自行车”的绿色交通出行方式受到了广大市民的青睐。由此，随之带来了“租车难、换车难”的新问题，即某些自行车租赁站点在一些时间段内的自行车数量不够多，使得用户不能租到自行车，而另外一些自行车站点在一些时间段内没有停车位，导致用户不能归还自行车。解决对公共自行车进行合理调度的问题，能提高公共自行车使用效率，从而提高用户对公共自行车使用的满意度，在促进公民绿色出行、缓解交通拥堵方面具有广泛的应用价值。
目前，主要的调度算法主要有：神经网络算法、蚁群算法、遗传算法等。神经网络法需要大量的训练数据，算法收敛慢，搜索能力不高；蚁群算法具有搜索时间长，容易陷入局部最优的缺点；遗传算法有计算速度慢，容易过早收敛等缺点。
发明内容
本发明所要解决的技术问题是克服现有技术的缺陷，提供一种公共自行车调度方法。
为解决上述技术问题，本发明提供一种公共自行车调度方法，包括以下步骤：
获取公共自行车本地数据，并进行预处理；
根据预处理后的数据预先构建公共自行车调度模型；
将公共自行车调度模型转化为伪旅行商模型；
采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解；
输出最优公共自行车调度模型的目标函数值，调度路径规划结果。
进一步的，获取公共自行车本地数据，并进行预处理的过程为：
读入公共自行车服务站点数量，服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵；读入服务站点的预约时间段和调度量；
读入公共自行车服务中心运输车数量，运输车装载量和运输车服务不同需求的服务站点的时间；
设定运输车以相同速度出发并设定运输车速度，由所得服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵，求出服务中心与服务站点之间的距离时间向量，服务站点与服务站点之间的距离时间矩阵。
进一步的，所述公共自行车调度模型包括以下限制条件：
在整个区域内公共自行车服务中心只有1个，运输车从服务中心出发，沿着某条路线逐个服务站点，最终回到服务中心，形成一个封闭的行驶路径；
每个站点接受服务的时长与站点的调度量有关，调度量越大，装载时间越长，每辆运输车装载自行车数量固定；
每个站点仅由一辆运输车服务，且一次完成对该站点的服务；
每个站点的调度量不能超过运输车最多可装载的自行车数量，且在调度过程中，运输车上的自行车数量不能超过其载运能力；
运输车服务站点的时间必须满足对应的时间窗要求，不允许站点在时间窗上限之后访问对象，在时间窗下限之前访问对象会产生等待成本，在时间窗之内访问对象不产生等待成本；
根据以上限制条件，建立公共自行车调度模型：
公式(1)表示运输车路程成本、服务时间、等待时间、站点服务数量、运输车使用数的加权总和取得最小值；
表示路程成本、
表示服务时间的惩罚、
表示等待时间惩罚、w3N1表示未服务站点惩罚，w4(K-K1)表示服务站点的运输车车辆数惩罚，运输车车辆数越多，成本越高；式(1)需要满足以下几个约束条件：
公式(2)表示第k辆运输车是否从第i个服务站点到第j个服务站进行服务，i，j都为0时，代表服务中心；k＝1,2…K，i,j＝1,2…Nc；
公式(3)表示第k辆运输车是否为第i个服务站点进行服务；
公式(4)表示每一个站点至多由一辆运输车进行服务，且一次完成对该站点的服务；
0≤Qk≤Q,k＝1,2,…K (5)
公式(5)表示第i辆运输车的车载量不得超过最多可装载的自行车数量Q；
|qi|≤Q，i＝1,2,…,Nc (6)
公式(6)表示第i个站点的调度量不得超过运输车最多可装载的自行车数量Q；
公式(7)表示运输车必须在站点要求的预约时间段上限之前到达预约地点；其中：
K1：未服务的运输车车辆数；
N1：未被服务站点的数量；
dij：第i个站点与第j个站点之间的距离矩阵；
waiti：第i辆运输车等待时间；
[ai,bi]：站点预约接受服务的时间；
N1：未被服务站点的数量；
K1：未服务的运输车车辆数；
dij：第i个站点与第j个站点之间的距离矩阵；
Tki：第i个站点接受第k辆运输车服务所需的服务时长；
waitk：第k个运输车的等待成本；
[ai,bi]：第i个站点接受服务的预约时间段；
Qk：第k辆运输车实际服务的装载量；
qi：第i个服务站点调度量，qi>0表示第i个服务站点从运输车调入qi辆自行车；qi＜0，表示运输车从第i个服务站点带走|qi|辆自行车；
wδ,δ＝1,2,3,4：目标函数各部分的加权系数。
进一步的，所述采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，包括以下步骤：
步骤1、设置概率参数p1、p2、p3、p4、p5、p6，随机生成NBSO个个体，每个个体有Nc个结点，对应Nc个服务站点，聚类数为m；
步骤2、根据老人和护工信息数据，计算所有个体的目标函数值，并采用基于目标空间聚类的策略，将NBSO个个体聚为m个类；
步骤3、对每个类中的个体按照目标函数值的大小进行排序，并从中选出目标函数值最小的个体作为该类的聚类中心；
步骤4、在0和1之间随机生成一个值r1，如果r1小于概率参数p1，则随机选择一个聚类中心；否则，随机产生一个个体代替该聚类中心；
步骤5、采用同一类中个体向类中随机个体学习的方法产生新个体；计算新产生个体的目标函数值，与当前个体相比，目标函数值小的个体作为下次迭代的新的个体；
步骤6、如果产生NBSO个新的个体，转至步骤7，否则转至步骤5；
步骤7、如果达到最大的迭代次数，则结束，否则跳至步骤2。
进一步的，所述步骤2和步骤5中均通过基于伪旅行商问题模型的启发式译码装置对目标函数进行计算，其译码方法包括以下步骤：
5a)置各辆运输车的初始时刻均为tinit，选择一辆运输车，置当前服务运输车号为k＝1，置当前选中结点号i＝1；
5b)设当前运输车的初始时间点为tm，且tm＝tinit，依次顺序选择对应未处理结点；置初始化车载量范围[lower,higher]为[0,Q]，即lower＝0,higher＝Q，设第k运输车的等待时间waitk＝0，临时等待时间ttemp＝0；
5c)如果k＞K，转至步骤5k)；
5d)如果当前结点对应服务站点的调度量大于0，则置lower＝lower+qi；否则，置higher＝higher+qi；
5e)如果lower＞high，即置k＝k+1，转至步骤5b)；
5f)根据运输车初始时间点tm和路程时间，计算运输车到达当前结点对应服务站点的时间点tn，如果该时间点在站点预约时间上限以上，置k＝k+1，转至步骤5b)；
5g)如果该时间点在站点预约时间下限以下，计算临时等待时间ttemp；
5h)更新运输车等待时间waitk＝waitk+ttemp；
5i)由运输车到达当前结点对应服务站点的时间点tn，等待时间，服务时间，求得运输车离开当前结点对应服务站点的时间点tm，并将该服务站点序号添加到运输车的服务列表；
5j)i＝i+1，判断结点号i是否小于Nc；若是，转至步骤5d)；
5k)计算目标函数，包括路程成本，等待成本，未服务站点的惩罚成本和提供服务的运输车车辆数的惩罚成本；
5l)输出目标函数值，调度路径规划结果，即各运输车服务站点的次序，所用运输车车辆数和各运输车初始所带公共自行车车辆数。
进一步的，步骤2中采用基于目标空间聚类的策略将个体聚类的步骤为：
6a)根据个体的目标函数值，对NBSO个个体目标函数值进行升序排列，得到排序后的个体目标函数值；
6b)由后一个个体目标函数值减去前一个个体的目标函数值，求得NBSO-1个相邻个体目标函数值的差值；
6c)找出m-1个最大的相邻个体目标函数值的差值，并记录各差值对应的序号；
6d)根据求得的m-1个序号进行切分，把已排序后个体中序号小于等于该序号的个体划分为一类，大于该序号的个体划分到另一类。
进一步的，步骤5中采用同一类中个体向类中随机个体学习的方法产生新个体，步骤为：
7a)在0和1之间随机生成一个值r2；
7b)如果r2小于概率参数p2：
⑥在m个聚类中随机选择一个类，并在0和1之间随机生成一个值r3；
⑦如果r3比概率参数p3小，选择该类聚类中心加上随机扰动产生一个新个体；
⑧否则，随机产生一个0到1间的数值r4；
⑨如果r4小于概率参数p4，从该聚类中随机选择一个个体，加上随机扰动产生一个新个体；
⑩否则，从该聚类中随机选择2个个体，互相融合加上随机扰动产生两个新个体；
7c)否则，在m个类中随机地挑选两个聚类产生新个体
⑥在0和1之间随机生成一个值r5：
⑦如果r5小于概率参数p5，将两个类的聚类中心互相融合加上随机扰动产生两个新个体；
⑧否则，随机产生一个0到1间的数值r6：
⑨如果r6小于概率参数p6，将第一个类的聚类中心与从第二个类中随机选择的一个个体互相融合加上随机扰动产生两个新个体；
⑩否则，分别从两个类中选择一个个体互相融合加上随机扰动产生两个新个体；
其中，步骤7b)步骤7c)中2个个体的相互融合采用交叉操作实现；个体施加随机扰动采用变异操作实现。
本发明所达到的有益效果：
本发明通过把公共自行车调度模型转化为伪旅行商问题，有效的消除了常用头脑风暴优化算法产生的不可行解问题；
本发明通过基于启发式的译码规则，合理的分配运输车资源，能够满足城市公共交通机构对运输车资源合理的调度，达到对运输车服务路径的优化。
附图说明
图1改进离散头脑风暴优化算法的流程图；
图2公共自行车预约调度与路径规划模型的启发式译码算法流程图；
图3公共自行车预约调度与路径规划模型转化为伪旅行商问题模型示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
一种公共自行车调度方法，包括步骤：
1)获取公共自行车本地数据，并进行预处理；
2)根据预处理后的数据预先构建公共自行车调度模型；
3)将公共自行车调度模型转化为伪旅行商模型；
4)采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解；
5)输出最优公共自行车调度模型的目标函数值，调度路径规划结果。
所述步骤1)数据预处理的具体过程为：
1a)读入公共自行车服务站点数量Nc，服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵；读入服务站点的预约时间段和调度量；
1b)读入公共自行车服务中心运输车数量K，运输车装载量Q和运输车服务不同需求的服务站点的时间；
1c)设定运输车以相同速度出发并设定运输车速度，由1a)所得服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵，求出服务中心与服务站点之间的距离时间向量，服务站点与服务站点之间的距离时间矩阵。
所述步骤2)根据预处理后的数据预先构建公共自行车调度模型，其建立过程包括步骤：
在整个区域内公共自行车服务中心只有1个，运输车从服务中心出发，沿着某条路线逐个服务站点，最终回到服务中心，形成一个封闭的行驶路径；
每个站点接受服务的时长与站点的调度量有关，调度量越大，装载时间越长，每辆运输车装载自行车数量固定；
每个站点仅由一辆运输车服务，且一次完成对该站点的服务；
每个站点的调度量不能超过运输车最多可装载的自行车数量，且在调度过程中，运输车上的自行车数量不能超过其载运能力；
运输车服务站点的时间必须满足对应的时间窗要求，不允许站点在时间窗上限之后访问对象，在时间窗下限之前访问对象会产生等待成本，在时间窗之内访问对象不产生等待成本；
根据以上限制条件，建立公共自行车调度模型：
公式(1)表示运输车路程成本、服务时间、等待时间、站点服务数量、运输车使用数的加权总和取得最小值；
表示路程成本、
表示服务时间的惩罚、
表示等待时间惩罚、w3N1表示未服务站点惩罚，w4(K-K1)表示服务站点的运输车车辆数惩罚，运输车车辆数越多，成本越高；式(1)需要满足以下几个约束条件：
公式(2)表示第k辆运输车是否从第i个服务站点到第j个服务站进行服务，i，j都为0时，代表服务中心；k＝1,2…K，i,j＝1,2…Nc；
公式(3)表示第k辆运输车是否为第i个服务站点进行服务；
公式(4)表示每一个站点至多由一辆运输车进行服务，且一次完成对该站点的服务；
0≤Qk≤Q,k＝1,2,…K (5)
公式(5)表示第i辆运输车的车载量不得超过最多可装载的自行车数量Q；
|qi|≤Q，i＝1,2,…,Nc (6)
公式(6)表示第i个站点的调度量不得超过运输车最多可装载的自行车数量Q；
公式(7)表示运输车必须在站点要求的预约时间段上限之前到达预约地点；其中：
K1：未服务的运输车车辆数；
N1：未被服务站点的数量；
dij：第i个站点与第j个站点之间的距离矩阵；
waiti：第i辆运输车等待时间；
[ai,bi]：站点预约接受服务的时间；
N1：未被服务站点的数量；
K1：未服务的运输车车辆数；
dij：第i个站点与第j个站点之间的距离矩阵；
Tki：第i个站点接受第k辆运输车服务所需的服务时长；
waitk：第k个运输车的等待成本；
[ai,bi]：第i个站点接受服务的预约时间段；
Qk：第k辆运输车实际服务的装载量；
qi：第i个服务站点调度量，qi>0表示第i个服务站点从运输车调入qi辆自行车；qi＜0，表示运输车从第i个服务站点带走|qi|辆自行车；
wδ,δ＝1,2,3,4：目标函数各部分的加权系数。
所述步骤3)将公共自行车调度模型转化为伪旅行商模型，具体过程为：将站点视为城市点，运输车视为旅行商人。问题的实质是寻找一条拜访每个城市的回路，路径的限制是每个城市只能拜访一次，而且最后要回到原来出发的城市。路径的选择目标是这条回路的目标函数值最小；
所述步骤4)采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，其具体做法包括如下步骤：
步骤1、设置概率参数p1、p2、p3、p4、p5、p6，随机生成NBSO个个体，每个个体有Nc个结点，对应Nc个服务站点，聚类数为m；
步骤2、根据老人和护工信息数据，计算所有个体的目标函数值，并采用基于目标空间聚类的策略，将NBSO个个体聚为m个类；
步骤3、对每个类中的个体按照目标函数值的大小进行排序，并从中选出目标函数值最小的个体作为该类的聚类中心；
步骤4、在0和1之间随机生成一个值r1，如果r1小于概率参数p1，则随机选择一个聚类中心；否则，随机产生一个个体代替该聚类中心；
步骤5、采用同一类中个体向类中随机个体学习的方法产生新个体；计算新产生个体的目标函数值，与当前个体相比，目标函数值小的个体作为下次迭代的新的个体；
步骤6、如果产生NBSO个新的个体，转至步骤7，否则转至步骤5；
步骤7、如果达到最大的迭代次数，则结束，否则跳至步骤2。
如图2所示，所述步骤4)采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，其特征在于，所述步骤2和步骤5中均通过基于伪旅行商问题模型的启发式译码装置对目标函数进行计算，其译码方法包括以下步骤：
a1)置各辆运输车的初始时刻均为tinit，选择一辆运输车，置当前服务运输车序号为k＝1，置当前选中结点号i＝1；
b1)设当前运输车的初始时间点为tm，且tm＝tinit，依次顺序选择对应未处理结点；置初始化车载量范围[lower,higher]为[0,Q]，即lower＝0,higher＝Q，设第k辆运输车的等待时间waitk＝0，临时等待时间ttemp＝0；
c1)如果k＞K，转至步骤k1)；
d1)如果当前结点对应服务站点的调度量大于0，则置lower＝lower+qi；否则，置higher＝higher+qi；
e1)如果lower＞high，即置k＝k+1，转至步骤b1)；
f1)根据运输车初始时间点tm和路程时间，计算运输车到达当前结点对应服务站点的时间点tn，如果该时间点在站点预约时间上限以上，置k＝k+1，转至步骤b1)；
g1)如果该时间点在站点预约时间下限以下，计算临时等待时间ttemp；
h1)更新运输车等待时间waitk＝waitk+ttemp；
i1)由运输车到达当前结点对应服务站点的时间点tn，等待时间，服务时间，求得运输车离开当前结点对应服务站点的时间点tm，并将该服务站点序号添加到运输车的服务列表；
j1)i＝i+1，判断结点号i是否小于Nc；若是，转至步骤5d)；
k1)计算目标函数，包括路程成本，等待成本，未服务站点的惩罚成本和提供服务的运输车车辆数的惩罚成本；
l1)输出目标函数值，调度路径规划结果，即各运输车服务站点的次序，所用运输车车辆数和各运输车初始所带公共自行车车辆数。
如图3所示，设一个个体共有10个结点，1～10是结点编号，其中，箭头代表运输车在服务完前一结点对应需服务站点后可接着服务下一结点的对应需服务站点，虚线代表在服务完前一结点对应需服务站点后不可接着服务下一结点的对应需服务站点，需要更换运输车。根据如上译码方式，可确定本个体10个结点对应需服务站点需要3个运输车进行服务，并得到各运输车的服务列表。其中第一个运输车的服务结点列表为[2,3,4],服务结点号为2,3,4对应的站点；第二个运输车的服务结点列表[5,7]，依次服务结点号为5,7对应的站点；第三个运输车的服务结点列表为[8,9,1,6,10]，依次服务结点号为8,9,1,6,10对应的站点。
所述步骤4)采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，其特征在于，步骤2中采用基于目标空间聚类的策略将个体聚类的步骤为：
a2)根据个体的目标函数值，对NBSO个个体目标函数值进行升序排列，得到排序后的个体目标函数值；
b2)由后一个个体目标函数值减去前一个个体的目标函数值，求得NBSO-1个相邻个体目标函数值的差值；
c2)找出m-1个最大的相邻个体目标函数值的差值，并记录各差值对应的序号；
d2)根据求得的m-1个序号进行切分，把已排序后个体中序号小于等于该序号的个体划分为一类，大于该序号的个体划分到另一类。
所述步骤4)采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，其特征在于：步骤5中采用同一类中个体向类中随机个体学习的方法产生新个体，步骤为：
a3)在0和1之间随机生成一个值r2；
b3)如果r2小于概率参数p2：
①在m个聚类中随机选择一个类，并在0和1之间随机生成一个值r3；
②如果r3比概率参数p3小，选择该类聚类中心加上随机扰动产生一个新个体；
③否则，随机产生一个0到1间的数值r4；
④如果r4小于概率参数p4，从该聚类中随机选择一个个体，加上随机扰动产生一个新个体；
⑤否则，从该聚类中随机选择2个个体，互相融合加上随机扰动产生两个新个体；
c3)否则，在m个类中随机地挑选两个聚类产生新个体
①在0和1之间随机生成一个值r5：
②如果r5小于概率参数p5，将两个类的聚类中心互相融合加上随机扰动产生两个新个体；
③否则，随机产生一个0到1间的数值r6：
④如果r6小于概率参数p6，将第一个类的聚类中心与从第二个类中随机选择的一个个体互相融合加上随机扰动产生两个新个体；
⑤否则，分别从两个类中选择一个个体互相融合加上随机扰动产生两个新个体；
其中，步骤b3)、步骤c3)中2个个体的相互融合采用交叉操作实现；个体施加随机扰动采用变异操作实现。
其中对于交叉操作，本发明中采用部分匹配交叉操作，具体步骤为：
产生两个[1,10]区间内的随机整数r1和r2，确定两个位置，对两个位置的中间数据进行交叉，如r1＝4,r2＝7
交叉为：
交叉后，同一个个体中有重复的编号，不重复的数字保留，有冲突的数字(带*位置)采用部分映射的方法消除冲突，即利用中间段的对应关系进行映射。结果为：
其中变异操作，本发明中采用序列倒置，具体操作如下：
产生两个[1,10]区间内的随机整数r1和r2，确定两个位置，将其对换位置，如r1＝4,r2＝7
序列倒置后为：
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.一种公共自行车调度方法，其特征在于，包括以下步骤：
获取公共自行车本地数据，并进行预处理；
根据预处理后的数据预先构建公共自行车调度模型；
将公共自行车调度模型转化为伪旅行商模型；
采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解；
输出最优公共自行车调度模型的目标函数值，调度路径规划结果；
所述公共自行车调度模型包括以下限制条件：
在整个区域内公共自行车服务中心只有1个，运输车从服务中心出发，沿着某条路线逐个服务站点，最终回到服务中心，形成一个封闭的行驶路径；
每个站点接受服务的时长与站点的调度量有关，调度量越大，装载时间越长，每辆运输车装载自行车数量固定；
每个站点仅由一辆运输车服务，且一次完成对该站点的服务；
每个站点的调度量不能超过运输车最多可装载的自行车数量，且在调度过程中，运输车上的自行车数量不能超过其载运能力；
运输车服务站点的时间必须满足对应的时间窗要求，不允许站点在时间窗上限之后访问对象，在时间窗下限之前访问对象会产生等待成本，在时间窗之内访问对象不产生等待成本；
根据以上限制条件，建立公共自行车调度模型：
公式(1)表示运输车路程成本、服务时间、等待时间、站点服务数量、运输车使用数的加权总和取得最小值；
表示路程成本、
表示服务时间的惩罚、
表示等待时间惩罚、w3N1表示未服务站点惩罚，w4(K-K1)表示服务站点的运输车车辆数惩罚，运输车车辆数越多，成本越高；式(1)需要满足以下几个约束条件：
公式(2)表示第k辆运输车是否从第i个服务站点到第j个服务站进行服务，i，j都为0时，代表服务中心；k＝1,2…K，i,j＝1,2…Nc；
公式(3)表示第k辆运输车是否为第i个服务站点进行服务；
公式(4)表示每一个站点至多由一辆运输车进行服务，且一次完成对该站点的服务；
0≤Qk≤Q,k＝1,2,…K (5)
公式(5)表示第i辆运输车的车载量不得超过最多可装载的自行车数量Q；
|qi|≤Q，i＝1,2,…,Nc (6)
公式(6)表示第i个站点的调度量不得超过运输车最多可装载的自行车数量Q；
公式(7)表示运输车必须在站点要求的预约时间段上限之前到达预约地点；其中：
K1：未服务的运输车车辆数；
N1：未被服务站点的数量；
dij：第i个站点与第j个站点之间的距离矩阵；
Tki：第i个站点接受第k辆运输车服务所需的服务时长；
waitk：第k个运输车的等待成本；
[ai,bi]：第i个站点接受服务的预约时间段；
Qk：第k辆运输车实际服务的装载量；
qi：第i个服务站点调度量，qi>0表示第i个服务站点从运输车调入qi辆自行车；qi＜0，表示运输车从第i个服务站点带走|qi|辆自行车；
wδ,δ＝1,2,3,4：目标函数各部分的加权系数。
2.根据权利要求1所述的公共自行车调度方法，其特征在于，获取公共自行车本地数据，并进行预处理的过程为：
读入公共自行车服务站点数量，服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵；读入服务站点的预约时间段和调度量；
读入公共自行车服务中心运输车数量，运输车装载量和运输车服务不同需求的服务站点的时间；
设定运输车以相同速度出发并设定运输车速度，由所得服务中心与服务站点之间的距离向量，服务站点与服务站点之间的距离矩阵，求出服务中心与服务站点之间的距离时间向量，服务站点与服务站点之间的距离时间矩阵。
3.根据权利要求1所述的公共自行车调度方法，其特征在于，所述采用改进离散头脑风暴优化算法对所得伪旅行商模型进行求解，包括以下步骤：
步骤1、设置概率参数p1、p2、p3、p4、p5、p6，随机生成NBSO个个体，每个个体有Nc个结点，对应Nc个服务站点，聚类数为m；
步骤2、根据老人和护工信息数据，计算所有个体的目标函数值，并采用基于目标空间聚类的策略，将NBSO个个体聚为m个类；
步骤3、对每个类中的个体按照目标函数值的大小进行排序，并从中选出目标函数值最小的个体作为该类的聚类中心；
步骤4、在0和1之间随机生成一个值r1，如果r1小于概率参数p1，则随机选择一个聚类中心；否则，随机产生一个个体代替该聚类中心；
步骤5、采用同一类中个体向类中随机个体学习的方法产生新个体；计算新产生个体的目标函数值，与当前个体相比，目标函数值小的个体作为下次迭代的新的个体；
步骤6、如果产生NBSO个新的个体，转至步骤7，否则转至步骤5；
步骤7、如果达到最大的迭代次数，则结束，否则跳至步骤2。
4.根据权利要求3所述的公共自行车调度方法，所述步骤2和步骤5中均通过基于伪旅行商问题模型的启发式译码装置对目标函数进行计算，其译码方法包括以下步骤：
5a)置各辆运输车的初始时刻均为tinit，选择一辆运输车，置当前服务运输车号为k＝1，置当前选中结点号i＝1；
5b)设当前运输车的初始时间点为tm，且tm＝tinit，依次顺序选择对应未处理结点；置初始化车载量范围[lower,higher]为[0,Q]，即lower＝0,higher＝Q，设第k运输车的等待时间waitk＝0，临时等待时间ttemp＝0；
5c)如果k＞K，转至步骤5k)；
5d)如果当前结点对应服务站点的调度量大于0，则置lower＝lower+qi；否则，置higher＝higher+qi；
5e)如果lower＞high，即置k＝k+1，转至步骤5b)；
5f)根据运输车初始时间点tm和路程时间，计算运输车到达当前结点对应服务站点的时间点tn，如果该时间点在站点预约时间上限以上，置k＝k+1，转至步骤5b)；
5g)如果该时间点在站点预约时间下限以下，计算临时等待时间ttemp；
5h)更新运输车等待时间waitk＝waitk+ttemp；
5i)由运输车到达当前结点对应服务站点的时间点tn，等待时间，服务时间，求得运输车离开当前结点对应服务站点的时间点tm，并将该服务站点序号添加到运输车的服务列表；
5j)i＝i+1，判断结点号i是否小于Nc；若是，转至步骤5d)；
5k)计算目标函数，包括路程成本，等待成本，未服务站点的惩罚成本和提供服务的运输车车辆数的惩罚成本；
5l)输出目标函数值，调度路径规划结果，即各运输车服务站点的次序，所用运输车车辆数和各运输车初始所带公共自行车车辆数。
5.根据权利要求3所述的公共自行车调度方法，其特征在于，步骤2中采用基于目标空间聚类的策略将个体聚类的步骤为：
6a)根据个体的目标函数值，对NBSO个个体目标函数值进行升序排列，得到排序后的个体目标函数值；
6b)由后一个个体目标函数值减去前一个个体的目标函数值，求得NBSO-1个相邻个体目标函数值的差值；
6c)找出m-1个最大的相邻个体目标函数值的差值，并记录各差值对应的序号；
6d)根据求得的m-1个序号进行切分，把已排序后个体中序号小于等于该序号的个体划分为一类，大于该序号的个体划分到另一类。
6.根据权利要求3所述的公共自行车调度方法，其特征在于，步骤5中采用同一类中个体向类中随机个体学习的方法产生新个体，步骤为：
7a)在0和1之间随机生成一个值r2；
7b)如果r2小于概率参数p2：
①在m个聚类中随机选择一个类，并在0和1之间随机生成一个值r3；
②如果r3比概率参数p3小，选择该类聚类中心加上随机扰动产生一个新个体；
③否则，随机产生一个0到1间的数值r4；
④如果r4小于概率参数p4，从该聚类中随机选择一个个体，加上随机扰动产生一个新个体；
⑤否则，从该聚类中随机选择2个个体，互相融合加上随机扰动产生两个新个体；
7c)否则，在m个类中随机地挑选两个聚类产生新个体
①在0和1之间随机生成一个值r5：
②如果r5小于概率参数p5，将两个类的聚类中心互相融合加上随机扰动产生两个新个体；
③否则，随机产生一个0到1间的数值r6：
④如果r6小于概率参数p6，将第一个类的聚类中心与从第二个类中随机选择的一个个体互相融合加上随机扰动产生两个新个体；
⑤否则，分别从两个类中选择一个个体互相融合加上随机扰动产生两个新个体；
其中，步骤7b)步骤7c)中2个个体的相互融合采用交叉操作实现；个体施加随机扰动采用变异操作实现。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911289386.1A CN111091286B (zh) | 2019-12-13 | 2019-12-13 | 一种公共自行车调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911289386.1A CN111091286B (zh) | 2019-12-13 | 2019-12-13 | 一种公共自行车调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111091286A CN111091286A (zh) | 2020-05-01 |
| CN111091286B true CN111091286B (zh) | 2022-08-16 |
Family
ID=70395017
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911289386.1A Active CN111091286B (zh) | 2019-12-13 | 2019-12-13 | 一种公共自行车调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111091286B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111090935B (zh) * | 2019-12-13 | 2022-09-23 | 河海大学常州校区 | 一种公共自行车预约调度与路径规划方法 |
| CN113203417B (zh) * | 2021-04-06 | 2022-08-19 | 广东安恒电力科技有限公司 | 基于头脑风暴优化算法的无人机巡检路径规划方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
| CN108805335A (zh) * | 2018-05-17 | 2018-11-13 | 杭州电子科技大学 | 一种公共自行车调度方法 |
-
2019
- 2019-12-13 CN CN201911289386.1A patent/CN111091286B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
| CN108805335A (zh) * | 2018-05-17 | 2018-11-13 | 杭州电子科技大学 | 一种公共自行车调度方法 |
Non-Patent Citations (2)
| Title |
|---|
| 公共自行车交通系统调度算法;刘登涛等;《计算机系统应用》;20111231;第20卷(第9期);第112-116页 * |
| 目标空间聚类的差分头脑风暴优化算法;吴亚丽等;《控制理论与应用》;20171231;第34卷(第12期);第1583-1593页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111091286A (zh) | 2020-05-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111862579B (zh) | 一种基于深度强化学习的出租车调度方法及系统 | |
| CN111785045B (zh) | 基于演员-评论家算法的分布式交通信号灯联合控制方法 | |
| CN105513400B (zh) | 动态规划出行路径的方法 | |
| CN108417031B (zh) | 一种基于Agent仿真的智能停车泊位预约策略优化方法 | |
| Jung et al. | Design and modeling of real-time shared-taxi dispatch algorithms | |
| CN111091286B (zh) | 一种公共自行车调度方法 | |
| CN113780808B (zh) | 基于柔性公交接驳系统线的车辆服务属性决策优化方法 | |
| CN109612488B (zh) | 基于大数据微服务的混合出行方式路径规划系统及方法 | |
| Zhao et al. | Dynamic path planning of emergency vehicles based on travel time prediction | |
| CN110019569B (zh) | 一种获取城市轨道交通运营状态信息的方法 | |
| CN112085249B (zh) | 一种基于强化学习的定制公交线路规划方法 | |
| CN111898793B (zh) | 一种组合出行方式下考虑用户感知差异性的路径选择方法 | |
| CN108960539A (zh) | 一种需求响应式接驳公交路径优化方法 | |
| CN111090935B (zh) | 一种公共自行车预约调度与路径规划方法 | |
| CN111445094B (zh) | 一种结合时间要求的快递车辆路径优化的方法及系统 | |
| CN114066503A (zh) | 一种基于构建虚拟服务片区进行出租车需求预测的方法 | |
| CN115936240A (zh) | 共享单车需求预测与投放调度方法 | |
| CN114742340A (zh) | 一种大规模路网中智能网联共享电动车充电站的优化布设求解方法 | |
| Zhao et al. | Robust Optimization of Mixed-Load School Bus Route Based on Multi-Objective Genetic Algorithm. | |
| Zhang et al. | The research on planning of taxi sharing route and sharing expenses | |
| Li et al. | Flexible Bus Route Setting and Scheduling Optimization Adapted to Spatial-temporal Variation of Passenger Flow. | |
| CN114444400A (zh) | 一种自主代客泊车群车库位分配方法 | |
| CN114037480A (zh) | 面向新城市的新能源车充电桩需求预测及部署优化方法 | |
| Lejdel | A conceptual framework for modeling smart parking | |
| CN113313412A (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 |