MATLAB 优化医疗点位置：最小化总距离

% 读取数据
pos = readtable('位置.xlsx');
link = readtable('连接道路.xlsx');

% 构建距离矩阵
n = size(pos,1);
d = zeros(n);
for i = 1:n
    for j = i+1:n
        idx = find((link.2==i & link.3==j) | (link.2==j & link.3==i));
        if ~isempty(idx)
            d(i,j) = link.4(idx);
            d(j,i) = d(i,j);
        end
    end
end

% 枚举所有可能的医疗点位置
min_S1 = inf;
for i = 1:n-2
    for j = i+1:n-1
        for k = j+1:n
            % 计算各村庄到这三个医疗点的距离
            dist = [d(:,i) d(:,j) d(:,k)];
            [min_dist,idx] = min(dist,[],2);
            % 计算总距离S1
            S1 = sum(min_dist);
            % 更新最小值
            if S1 < min_S1
                min_S1 = S1;
                ans_idx = [i j k];
                ans_dist = dist;
            end
        end
    end
end

% 输出结果
fprintf('问题1的答案为：\n');
fprintf('医疗点位置分别为：%d %d %d\n',ans_idx);
fprintf('总距离S1为：%.1f\n',min_S1);

% 作图
figure;
hold on;
plot(pos.2,pos.3,'ko','MarkerSize',4);
plot(pos.2(ans_idx),pos.3(ans_idx),'ro','MarkerSize',8,'LineWidth',2);
for i = 1:n
    [~,min_idx] = min(ans_dist(i,:));
    switch min_idx
        case 1
            plot([pos.2(i) pos.2(ans_idx(1))],[pos.3(i) pos.3(ans_idx(1))],'r-');
        case 2
            plot([pos.2(i) pos.2(ans_idx(2))],[pos.3(i) pos.3(ans_idx(2))],'r-');
        case 3
            plot([pos.2(i) pos.2(ans_idx(3))],[pos.3(i) pos.3(ans_idx(3))],'r-');
    end
end
axis equal;
title('问题1的结果');

代码解释：

1. 读取数据：
   • 使用 readtable 函数读取两个 Excel 文件：'位置.xlsx' 和 '连接道路.xlsx'，分别存储村庄坐标和道路信息。
2. 构建距离矩阵：
   • 使用嵌套循环遍历所有村庄对，查找连接它们的道路信息，并将距离存储在距离矩阵 d 中。
3. 枚举优化算法：
   • 使用三层嵌套循环枚举所有可能的三个医疗点组合。
   • 对于每个组合，计算每个村庄到最近医疗点的距离，并累加所有村庄的距离得到总距离 S1。
   • 使用 min_S1 变量记录最小总距离，并保存对应的医疗点位置 ans_idx 和距离矩阵 ans_dist。
4. 输出结果：
   • 使用 fprintf 函数输出最小总距离和对应的医疗点位置。
5. 作图：
   • 使用 plot 函数绘制村庄坐标和医疗点位置。
   • 使用 plot 函数绘制每个村庄到其最近医疗点的连线。

说明：

• 代码假设所有村庄之间都有连接道路。
• 代码可以根据实际情况进行修改，例如增加更多医疗点、考虑其他优化算法等。
• 该代码使用 MATLAB 语言实现，需要在 MATLAB 环境中运行。

改进建议：

• 可以使用更有效的优化算法，例如遗传算法或模拟退火算法。
• 可以考虑其他因素，例如医疗点之间的距离、医疗点附近的人口密度等。
• 可以将代码封装成函数，方便调用。
• 可以使用更清晰的注释和代码格式，提高代码的可读性。