假设某山区中有100个村庄现在要在村庄中建立几个医疗点方便村民看病。图1中给出这100个村庄的位置及可选道路连接示意图。附件数据的位置表单给出了这100个村庄的坐标单位：米附件数据的连接道路表单给出了可供选择的道路。现在要在100个村庄中建立3个医疗点并在可选道路中根据需要进行部分道路维修假定村民看病都选择维修后的道路。 如果各村庄村民到医疗点的距离太远不便于看病因此站在村民角度出发希望各村庄村民

代码实现如下：

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

% 构建邻接矩阵 n = height(pos); adj = inf(n); for i = 1:height(roads) adj(roads{i,2}, roads{i,3}) = pdist2(pos{i,2:3}, pos{roads{i,3},2:3}); adj(roads{i,3}, roads{i,2}) = pdist2(pos{i,2:3}, pos{roads{i,3},2:3}); end

% 求解k中心斯坦纳树 k = 3; [st, S] = kspantree(adj, k);

% 输出结果 disp('医疗点位置：'); disp(st); disp(['总距离：', num2str(S)]);

% 绘制维修道路图形 G = graph(adj); highlight(G, st, 'NodeColor', 'r', 'MarkerSize', 10); highlight(G, find(st), 'NodeColor', 'g', 'MarkerSize', 10); highlight(G, G.Edges(S), 'LineWidth', 2); plot(G, 'XData', pos.X坐标, 'YData', pos.Y坐标); title(['k中心斯坦纳树，k=', num2str(k), ', S1+S2=', num2str(S)])