多目标应用：多目标蜣螂优化算法求解多旅行商问题（旅行商个数和起点城市可修改）

IT猿手

已于 2023-06-26 21:26:42 修改

阅读量1.8k

点赞数 2

分类专栏：多目标应用 TSP MATLAB 文章标签：算法 TSP MTSP 人工智能

于 2023-02-07 14:30:04 首次发布

本文链接：https://blog.csdn.net/weixin_46204734/article/details/128917784

版权

MATLAB 同时被 3 个专栏收录

546 篇文章 125 订阅

订阅专栏

多目标应用

112 篇文章 11 订阅

订阅专栏

TSP

64 篇文章 1 订阅

订阅专栏

一、多旅行商问题

多旅行商问题（Multiple Traveling Salesman Problem, MTSP）是著名的旅行商问题（Traveling Salesman Problem, TSP）的延伸，多旅行商问题定义为：给定一个𝑛座城市的城市集合，指定𝑚个推销员，每一位推销员从起点城市出发访问一定数量的城市，最后回到终点城市，要求除起点和终点城市以外，每一座城市都必须至少被一位推销员访问，并且只能访问一次，需要求解出满足上述要求并且代价最小的分配方案，其中的代价通常用总路程长度来代替，当然也可以是时间、费用等。围绕着各推销员的起始点和终止点来划分，多旅行商问题大致可以分为四种，其中单仓库多旅行商问题是其中一种。多旅行商问题
单仓库多旅行商问题（Single-Depot Multiple Travelling Salesman Problem, SD-MTSP）：𝑚个推销员从同一座中心城市出发，访问其中一定数量的城市并且每座城市只能被某一个推销员访问一次，最后返回到中心城市，通常这种问题模型被称之为SD-MTSP。

1.1多旅行商多目标问题描述

对于推销员人数为𝑚，所需要访问的城市总数为𝑛的多旅行商问题而言，记中心城市为0，其双目标的数学模型可以表述为：
$\min f=\left[f_{1}, f_{2}\right]^{T}$
$\begin{array}{c} f_{1}=\sum_{k=1}^{m} \sum_{i=0}^{n} \sum_{j=0}^{n} d_{i j} x_{i j k} \\ f_{2}=\max _{1 \leq k \leq m} \sum_{i=0}^{n} \sum_{j=0}^{n} d_{i j} x_{i j k}-\min _{1 \leq k \leq m} \sum_{i=0}^{n} \sum_{j=0}^{n} d_{i j} x_{i j k} \end{array}$
其中：
$\begin{array}{l} x_{i j k}=\left\{\begin{array}{lc} 1 & \text { 推销员 } k \text { 由城市 } i \text { 到达城市 } j \\ 0 & \text { 否则 } \end{array}\right. \\ \text { s.t. } \\ \sum_{k=1}^{m} \sum_{i=1}^{n} x_{i 0 k}=m \\ \sum_{k=1}^{m} \sum_{j=1}^{n} x_{0 j k}=m \\ \sum_{k=1}^{m} \sum_{i=1}^{n} x_{i j k}=1 \quad \forall j=1, \ldots, n \\ \sum_{k=1}^{m} \sum_{j=1}^{n} x_{i j k}=1 \quad \forall i=1, \ldots, n \\ \end{array}$
其中， $f_{1}$ 表示所有推销员的总路程， $f_{2}$ 代表推销员中最长路线与最短路线的差值（平衡度）。约束条件1和约束条件2表示编号为 0 的中心城市的入度和出度都必须为𝑚，即代表𝑚个推销员均从中心城市出发并且完成行程后都回到了中心城市。约束条件3和约束条件4表示除中心城市外的其他所有城市的出度和入度均为 1，即每个城市能且只能被推销员中的一个人访问。

参考文献：
[1]杨帅. 求解多旅行商问题的进化多目标优化和决策算法研究[D].武汉科技大学,2020.

二、多目标蜣螂优化算法

非支配排序的蜣螂优化算法（Non-Dominated Sorting Dung beetle optimizer，NSDBO）

三、NSDBO求解多旅行商问题

本文选取国际通用的TSP实例库TSPLIB中的测试集bayg29，bayg29中城市分布如下图所示：
在这里插入图片描述

本文采用NSDBO求解bayg29，两个目标函数分别是： $f_{1}$ 表示所有推销员的总路程， $f_{2}$ 代表推销员中最长路线与最短路线的差值（平衡度）。设置城市1作为多个旅行商的起点城市，旅行商和起点可自行修改。部分代码如下：

Tnum=5;%旅行商个数（可以自行更改）
StartPoint=1; %选择起点城市（可以自行更改）
TestProblem=1;%1
MultiObj = GetFunInfo(TestProblem);
MultiObjFnc=MultiObj.name;%问题名
% Parameters
params.Np = 100;        % Population size
params.Nr = 200;        % Repository size
params.maxgen =100;    % Maximum number of generations
numOfObj=MultiObj.numOfObj;%目标函数个数
D=MultiObj.nVar;%维度
f = NSDBO(params,MultiObj);
X=f(:,1:D);%PS
Obtained_Pareto=f(:,D+1:D+numOfObj);%PF
save Obtained_Pareto Obtained_Pareto
save X X