draw topology in matlab

figi = figure(1);
seed = 42;  % select a seed
rng(seed, 'twister');router_num = 5;
host_num = 0;     % zero means using the default value
edge_spds_r = [1e6 2e6 4e6];  % edge rates for routers
% dropprobs_r = 0.1;
dropprobs_r = 0.1;    % 原来是0edge_spds_rh = Inf;
dropprobs_rh = 0.2;     % 原来是0attmt_assumpt = 'evenly'; % attachment assumption
topotype = 'grid25';
% topotype = 'twohop';
[A, Aspds, Adrps, ur, vh, topoinfo] = genA(topotype, edge_spds_r, dropprobs_r, router_num, ...host_num, attmt_assumpt, edge_spds_rh, dropprobs_rh);router_num = topoinfo.router_num;
A = triu(A(1:router_num, 1:router_num));
[u, v] = find(sparse(A));
edges = [u, v];Aspds = triu(Aspds(1:router_num, 1:router_num));xy = [topoinfo.x_r(:) topoinfo.y_r(:)];x = xy(:, 1);y = xy(:, 2);G = graph(edges(:, 1), edges(:, 2));idxOut = findedge(G, edges(:, 1), edges(:, 2));Ghdl = plot(G, 'XData', x, 'YData', y, 'NodeLabel', {}, 'EdgeColor', 'k', ...'MarkerSize', 12, 'NodeColor', 'yellow');
% scatter(x, y, 20, 'MarkerEdgeColor', 'k',...
%               'MarkerFaceColor', 'yellow',...
%               'LineWidth', .5);
% labeledge(Ghdl, [3 12 7 16], [3 12 7 16]);
% Ghdl.EdgeFontSize = 14;[u, v] = find(sparse(Aspds == edge_spds_r(1)));
edges_lspds = [u, v];idxOut_lspds = findedge(G, edges_lspds(:, 1), edges_lspds(:, 2));[u, v] = find(sparse(Aspds == edge_spds_r(2)));
edges_mspds  = [u, v];idxOut_mspds = findedge(G, edges_mspds(:, 1), edges_mspds(:, 2));idxOut_hspds = setdiff(idxOut, [idxOut_lspds; idxOut_mspds]);highlight(Ghdl, edges_lspds(:, 1), edges_lspds(:, 2), 'EdgeColor', 'g');
highlight(Ghdl, edges_mspds(:, 1), edges_mspds(:, 2), 'EdgeColor', 'b', 'LineWidth', 1.5);[sOut,tOut] = findedge(G, idxOut_hspds);highlight(Ghdl, sOut, tOut, 'LineWidth', 2);circles(x, y, 3, 'edgecolor','k','facecolor','none', 'facealpha', 1);axis equal% highlight(Ghdl, 1:numel(x), 'MarkerSize', 10);
text(x, y, num2cell(1:numnodes(G)), 'FontSize', 8, 'HorizontalAlignment', 'center', 'VerticalAlignment', 'middle');
axis equal
axis tight
axis('off');fname = topotype;
hold onL(1) = plot(nan, nan, 'g-');
L(2) = plot(nan, nan, 'b-', 'LineWidth', 1.5);
L(3) = plot(nan, nan, 'k-', 'LineWidth', 2);legend(L, {'1 Mbps', '2 Mbps', '4 Mbps'}, 'location', 'eastoutside', 'FontSize', 8)% plot(rand(3))
% lgd = legend('a','b','c');
% c = lgd.TextColor;
% lgd.TextColor = 'red';exportgraphics(figi, [fname '.png'],'Resolution', 600);