[发明专利]一种用于多配送中心抗干扰车辆调度方法

权利要求书

1.一种用于多配送中心抗干扰车辆调度方法，其特征在于，包括如下操作步骤：

步骤1、确定配送营运参数；

步骤2、设置算法参数；

步骤3、初始化种群；

步骤4、对pop进行局部搜索；

步骤5、对pop′进行复制算子；

步骤6、对Rpop进行非支配排序；

步骤7、令新群体P_t+1为空集；

步骤8、计算F_i中每个个体的拥挤距离，并建立F_i的偏序关系；

步骤9：将第i个分类子集F_i中选取个体填入P_t+1，直至P_t+1中个体数目为M；

步骤10：对P_t+1进行局部搜索操作，得到种群P′_t+1；

步骤11：对P′_t+1进行复制算子操作，得到种群Q_t+1；

步骤12：将P_t+1和Q_t+1对Rpop更新；

步骤13：如果当前代数t不超过最大迭代次数iter，重复步骤6-12；否则，输出Pareto前沿解集。

2.根据权利要求1所述的用于多配送中心抗干扰车辆调度方法，其特征在于：所述步骤1确定配送营运参数：根据配送营运的实际情况，给出仓库所在位置、客户所在位置，配送经过路段的最大行驶时间，配送经过路径的旅行时间可能发生异常波动的数量，客户规定的配送时间窗等信息；

旅行时间的干扰信息表示如下：

式(1)中的符号含义如下：

t_i，j表示配送车辆从第i个客户到第j个客户的旅行时间；t_i，j在之间任意取值；表示第i个客户到第j个客户的名义旅行时间；θ_i，j表示第i个客户到第j个客户的名义旅行时间相对于的最大偏移量；ρ_i，j表示旅行时间发生波动的不确定水平；Γ表示配送经过路径的旅行时间可能发生异常波动的数量；

明确优化目标如下：

Minimize ∑_i∈V∑_j∈V∑_k∈Kt_i，jx_i，j，k (3)

式(2)和(3)中，x_i，j，k表示如果第k辆车从第i个客户到第j个客户，则x_i，j，k取值为1，反之，取值为0；s_i，k表示第k辆车到达第i个客户的时刻，LT_i表示第i个客户允许的最晚到达时间，ET_i表示第i个客户允许的最晚到达时间；w₁表示单位正惩罚成本，表示配送车辆晚于客户允许的到达时间提供配送服务引起的正惩罚成本的总和；w₂表示单位负惩罚成本；令表示总惩罚成本；TTC＝∑_i∈V∑_j∈V∑_k∈Kt_i，jx_i，j，k表示总运输成本。

3.根据权利要求1所述的用于多配送中心抗干扰车辆调度方法，其特征在于：所述步骤2设置算法参数：种群大小(M)、初始化因子(α)、比例因子(d)；旅行时间在U_T的端点处取值。

4.根据权利要求1所述的用于多配送中心抗干扰车辆调度方法，其特征在于：所述步骤3种群初始化：为了解决配送中心车辆调度普遍存在的二律背反现象，分别设计了及时点优先的启发式规则，记为PTVH，以及最小生成树优先的启发式规则，记为PMSTH；采用2个启发式规则产生M个可行解，形成初始种群；首先给定一个常数α∈[0，1]，随机生成一个实数β∈[0，1]，如果α不小于β，利用PTVH产生一个可行解；反之，利用PMSTH产生一个可行解；重复上述过程，直到产生M个可行解；因此，采用这种初始解生成策略可以有效保证种群的多样性和初始解的质量；该过程记为init_pop。