帮我生成matlab代码，要求如下：选取一种拓扑控制算法，编程实现该算法并用图形显示效果。要 求有执行算法前后拓扑对比图、链路数量统计对比（拓扑发现算法要显示拓扑最 终状态并对节点角色进行标注）。节点数目最少 50 个，随机分布，其他如通信半 径、部署区域大小等自行选取。考虑到连通性问题，要求应用拓扑控制的初始网 络为全连通网络，即每个节点与其他任意节点至少存在一条链路（可在随机生成 节点位置时多试几次或增加节点密度）

以下是一个基于最小生成树拓扑控制算法的Matlab代码，可以生成50个节点的随机网络，并使用拓扑控制算法进行优化，最终显示效果。

% 生成50个节点的随机网络
num_nodes = 50;
min_distance = 0.1;
max_distance = 0.5;
positions = rand(num_nodes, 2);
distances = pdist2(positions, positions);
distances(distances < min_distance) = inf;
distances(distances > max_distance) = inf;

% 确保网络是连通的
while ~all(sum(isfinite(distances)) >= num_nodes)
    positions = rand(num_nodes, 2);
    distances = pdist2(positions, positions);
    distances(distances < min_distance) = inf;
    distances(distances > max_distance) = inf;
end

% 生成初始拓扑图
G = graph(distances);
subplot(2,2,1);
plot(G, 'XData', positions(:,1), 'YData', positions(:,2));
title('初始拓扑图');

% 执行拓扑控制算法
MST = minspantree(G);
distances = distances .* ismember(G.Edges.EndNodes, MST.Edges.EndNodes);
G = graph(distances);
subplot(2,2,2);
plot(G, 'XData', positions(:,1), 'YData', positions(:,2));
title('拓扑控制后的拓扑图');

% 显示链路数量统计对比
num_links_before = numedges(G);
num_links_after = numedges(MST);
subplot(2,2,3);
bar([num_links_before, num_links_after]);
set(gca, 'xticklabel', {'拓扑控制前', '拓扑控制后'});
ylabel('链路数量');
title('链路数量统计对比');

% 显示拓扑最终状态并对节点角色进行标注
bins = conncomp(MST);
subplot(2,2,4);
gscatter(positions(:,1), positions(:,2), bins);
title('拓扑最终状态并对节点角色进行标注');

在上述代码中，我们首先生成50个随机节点，并确保网络是连通的。然后，我们使用Matlab自带的图形库来绘制初始拓扑图。接着，我们执行最小生成树算法，得到一个基本的拓扑控制方案，并用该方案优化网络拓扑，最后再用图形显示效果。我们还绘制了一个条形图来显示拓扑控制前后的链路数量对比，并使用`conncomp`函数对节点进行分组，以便更好地显示拓扑最终状态并对节点角色进行标注。

你可以根据需要选择其他的拓扑控制算法来编程实现，但需要注意的是，Matlab的图形库可能不够灵活，因此你可能需要另外的工具来实现更复杂的可视化效果。