1.一种基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，包括以下步骤：S1确认配送中心、配送车辆数量、各个车辆的容量约束、总配送距离以及选择多个任务点；

S2采用GWOPSO算法，即Grey Wolf Optimizer & Particle Swarm Optimization算法求解多目标带时间窗车辆路径规划，进而得到最终的规划路径以及车辆的分配情况；

所述GWOPSO算法包括以下步骤：

S21设置算法参数，主要包括种群大小N、最大迭代次数Maxt、外部存档数目a、划分网格数量b；

S22采用最短距离规则和最大满意度规则来实施种群的初始化，并对种群个体进行编码操作和计算目标函数值，然后将这些个体加入到外部存档中，最后刷新外部存档；

S23采用多策略方式从外部存档中选出3个最优个体 、 和 ；

S24采用灰狼算法或粒子群算法对种群中的所有个体进行更新，将更新后的个体并入外部存档中，并刷新外部存档；

S25选择外部存档中最优的一些个体进行局部搜索，并将得到的新个体并入外部存档，然后刷新外部存档；

S26如果算法的迭代次数小于最大迭代次数，则转到步骤S23进行循环运算；否则，转到步骤S27；

S27随机从外部存档中选择一定比例的个体进行全局搜索，并将全局搜索得到的新个体并入外部存档，然后刷新外部存档；

S28输出外部存档中三个最优个体对应的车辆路径规划方案。

2.根据权利要求1所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述最短距离规则为：当前车辆从配送中心出发，将满足当前车辆容量约束并能在规定时间窗口内服务的任务点记录下来，且最短距离服务原则为：选择距离最小的任务点作为下一个被服务的任务点，当当前车辆的剩余容量不满足下一个任务点的需求时，当前车辆将返回配送中心，配送中心将派遣新车辆，所述新车辆继续按照上述的最短距离服务原则对任务点进行服务，直到所有任务点被服务；

所述最大满意度规则为：车辆从配送中心出发，将满足车辆容量约束并能在规定时间窗口内服务的任务点记录下来，且满意度最大服务原则为：选择满意度最大的任务点作为下一个被服务的任务点，当车辆的剩余容量不满足下一个任务点的需求时，车辆将返回配送中心，配送中心将派遣新车辆，继续按照上述的满意度最大服务原则对任务点进行服务，直到所有任务点被服务。

3.根据权利要求1所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述目标函数为：；

其中，CT表示车辆违反时间窗口的成本，CU表示车辆启用成本，CD表示车辆的运输成本，I表示任务点的数量， 为第i个任务点用户的满意度， 为总成本， 为用户的总满意度。

4.根据权利要求1所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述刷新外部存档为取出外部存档中所有的个体，组成一个解集合，然后从这个解集合中按照解对应的目标函数值选取所有的非劣解组成非劣解集合，并根据非劣解集合中各个解的密度进行排序，密度小的解排在前面，最后选择排在非劣解集合前面的a个解保存到外部存档中。

5.根据权利要求1所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述步骤S23具体包括：S231利用网格技术对外部存档中的个体进行区域划分，将双目标形成的二维空间划分为多个网格，并计算每个个体的网格坐标；

S232计算个体的拐点距离和加权分布距离；

其中，拐点距离计算方法为：根据外部存档中的2个极端个体，确定一条直线，后计算外部存档每个个体距离该直线的距离；

个体 的加权分布距离计算公式如下：

；

其中，M为目标函数的个数， 为个体 在第m维度上的目标函数值， 为同一网格中拐点距离最大的个体；

S233采用距离最大策略或加权分布距离最大策略来选择最优个体，即：则从两个极端个体所在网格，分别选取网格中拐点距离值或加权分布距离值中最大的个体为最优个体和 ，并选择外部存档中拐点距离值或加权分布距离值最大的个体为最优个体 。

6.根据权利要求5所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，步骤S24中，对种群中的所有个体进行更新包括：选择灰狼算法或粒子群算法更新种群中所有个体；

采用的灰狼算法更新个体过程为：

若当前迭代次数为t，则种群中的每个个体对应一头狼，设狼群的个数为种群数目N，每头狼的位置为 ，并指定 、 和 狼为步骤S23选择的3个最优个体 、 和 ，它们的位置分别为 ；

计算第i头狼与 狼的距离 、第i头狼与 狼的距离 和第i头狼与 狼的距离，表示为：；

然后计算第i头狼向 、 和 狼移动的趋势项，分别记为 ：；

式中， ， 和 为

[0,1]上的随机数，j=1,2,3；

则第i头狼的位置更新为：

；

采用的粒子群算法更新个体公式如下：

；

；

；

其中， 表示第i个个体第t次迭代时的速度， 和 为学习因子；t为当前迭代次数，Maxt为算法最大迭代次数， 为惯性系数， 选0.4， 选0.9；

、 和 为步骤S23选择的3个最优个体 、 和 的位置， 为第i个个体更新后的位置。

7.根据权利要求6所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述步骤S25具体包括：S251选择外部存档中最优的3个个体，对这些个体的编码进行映射，得到对应的车辆路径；

S252对于最优个体对应的路径，先进行路径内部局部搜索操作，搜索次数为任务点数目，再进行路径间局部搜索操作，搜索次数为两条路径中最多的任务点数目，若局部搜索得到的路径优于原路径，则代替原路径；

S252将局部搜索得到的个体并入外部存档，然后刷新外部存档。

8.根据权利要求7所述的基于混合灰狼粒子群算法的车辆路径规划方法，其特征在于，所述路径内部局部搜索操作包括：

1）选取1条路径R1，随机交换路径内2个任务点的位置；

2）选取1条路径R1，随机选择2个任务点，将2个任务点间的路径进行逆转；

3）选取1条路径R1，随机选择一个任务点，将其插入到其它任务点前；

所述路径间局部搜索操作包括：

4）选取2条路径R1和R2，分别在R1和R2中随机选择任务点 和 ，交换 和 ；

5）选取2条路径R1和R2，在R1中随机选择2个连续任务点 和 ，在R2中随机选择任务点 ，将 和 插入到 前；

6）选取2条路径R1和R2，在R1中随机选择2个连续任务点 和 ，在R2中随机选择2个连续任务点 和 ，交换 、 和 、 ；

7）选取2条路径R1和R2，在R1中随机选择2个连续任务点 和 ，在R2中随机选择任务点 ，交换 、 和 。

9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现上述权利要求1‑8任一项所述的方法。

10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现上述权利要求1‑8任一项所述的方法。