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

以下是用Matlab编程实现RNG算法的示例代码：

% 设置节点数量和通信半径
n = 50;
r = 0.3;

% 生成随机节点
x = rand(1, n);
y = rand(1, n);

% 建立全连通网络
G = complete_graph(n);

% 绘制执行算法前的拓扑对比图
figure;
subplot(1, 2, 1);
plot(G, 'XData', x, 'YData', y, 'MarkerSize', 10);
title('Initial Topology');

% 计算节点之间的欧几里得距离
dist = pdist2([x', y'], [x', y']);

% 确定节点之间是否应该建立链路
adj = dist < r;

% 生成随机的链路
G = graph(adj);
G = simplify(G);

% 绘制执行算法后的拓扑对比图
subplot(1, 2, 2);
plot(G, 'XData', x, 'YData', y, 'MarkerSize', 10);
title('Final Topology');

% 统计链路数量的变化
initial_edges = numedges(complete_graph(n));
final_edges = numedges(G);
fprintf('Initial Number of Edges: %d\n', initial_edges);
fprintf('Final Number of Edges: %d\n', final_edges);

% 统计链路数量的变化并绘制统计图
figure;
bar([initial_edges, final_edges]);
title('Edge Count Comparison');
xlabel('Topology');
ylabel('Edge Count');
set(gca, 'xticklabel', {'Initial', 'Final'});

这段代码生成50个随机节点，并通过计算节点之间的欧几里得距离来确定节点之间是否应该建立链路。然后，它使用igraph库来生成随机的链路。最后，它绘制了执行算法前后的拓扑对比图，并统计了链路数量的变化，并将其显示在一个统计图中。