已知21个城市之间的城市距离矩阵，编写matlab代码使得找到一条遍历所有城市的最短路径使用DJ算法解决上问题

以下是使用DJ算法解决21个城市之间的最短路径问题的Matlab代码：

% 城市距离矩阵
D = [0 17 32 46 32 19 25 18 42 31 29 31 42 25 36 28 27 42 38 22 11;
    17 0 15 30 22 3 8 15 25 16 15 18 27 10 22 13 14 25 21 6 8;
    32 15 0 15 35 18 9 28 12 20 17 23 14 18 9 17 23 13 7 23 23;
    46 30 15 0 50 32 25 42 24 37 34 41 24 32 26 32 38 17 16 35 39;
    32 22 35 50 0 19 30 19 54 47 45 45 53 36 47 40 39 54 49 33 23;
    19 3 18 32 19 0 6 12 22 14 14 16 25 8 19 10 11 22 18 3 5;
    25 8 9 25 30 6 0 21 10 13 12 14 20 5 16 7 9 20 16 5 8;
    18 15 28 42 19 12 21 0 47 38 36 36 45 28 39 31 30 45 41 25 15;
    42 25 12 24 54 22 10 47 0 25 22 29 12 24 17 23 29 9 5 31 34;
    31 16 20 37 47 14 13 38 25 0 4 8 19 13 6 6 12 17 12 16 22;
    29 15 17 34 45 14 12 36 22 4 0 7 20 11 3 5 11 13 9 13 19;
    31 18 23 41 45 16 14 36 29 8 7 0 23 12 10 13 17 15 11 15 19;
    42 27 14 24 53 25 20 45 12 19 20 23 0 18 10 16 22 4 7 33 32;
    25 10 18 32 36 8 5 28 24 13 11 12 18 0 12 4 5 16 12 3 10;
    36 22 9 26 47 19 16 39 17 6 3 10 10 12 0 10 16 10 6 24 27;
    28 13 17 32 40 10 7 31 23 6 5 13 16 4 10 0 2 15 11 12 18;
    27 14 23 38 39 11 9 30 29 12 11 17 22 5 16 2 0 19 15 17 21;
    42 25 13 17 54 22 20 45 9 17 13 15 4 16 10 15 19 0 4 29 32;
    38 21 7 16 49 18 16 41 5 12 9 11 7 12 6 11 15 4 0 26 29;
    22 6 23 35 33 3 5 25 31 16 13 15 33 3 24 12 17 29 26 0 12;
    11 8 23 39 23 5 8 15 34 22 19 19 32 10 27 18 21 32 29 12 0];

% DJ算法求解最短路径
n = size(D, 1); % 城市数目
S = zeros(1, n); % S集合，表示已遍历的城市
S(1) = 1; % 从第一个城市开始遍历
d = Inf(1, n); % d数组，表示起点到各个城市的距离
d(1) = 0; % 起点到自身的距离为0
for i = 2:n % 初始化d数组
    d(i) = D(1, i);
end
for i = 2:n % 遍历n-1个城市
    minDist = Inf; % 最短距离初始化为无穷大
    for j = 2:n % 在未遍历的城市中寻找距离最短的城市
        if S(j) == 0 && d(j) < minDist
            minDist = d(j);
            u = j;
        end
    end
    S(u) = 1; % 标记该城市已遍历
    for v = 2:n % 更新d数组
        if S(v) == 0 && d(v) > d(u) + D(u, v)
            d(v) = d(u) + D(u, v);
        end
    end
end
% 输出最短路径和距离
path = [1]; % 路径数组，初始化为起点
dist = 0; % 距离初始化为0
while length(path) < n % 遍历所有城市
    minDist = Inf; % 最短距离初始化为无穷大
    for i = 2:n % 在未遍历的城市中寻找距离最短的城市
        if S(i) == 0 && d(i) < minDist
            minDist = d(i);
            u = i;
        end
    end
    S(u) = 1; % 标记该城市已遍历
    path = [path u]; % 将该城市加入路径
    dist = dist + D(path(end-1), path(end)); % 更新距离
end
dist = dist + D(path(end), path(1)); % 加上回到起点的距离
disp('最短路径：');
disp(path);
disp('最短距离：');
disp(dist);

这段代码中，我们首先定义了城市距离矩阵 `D`，然后使用DJ算法求解最短路径。在求解过程中，我们使用了两个数组 `S` 和 `d`，分别表示已遍历的城市和起点到各个城市的距离。在每次遍历时，我们选择距离起点最近的未遍历城市作为下一个遍历的城市，并更新 `d` 数组。最后，我们根据 `S` 和 `d` 数组求出最短路径和距离。