i want to get adjacency matrix of a network

58 ビュー (過去 30 日間)
mks
mks 2023 年 7 月 30 日
コメント済み: Steven Lord 2023 年 7 月 30 日
i have write this but is not running
clc;
clear all
n=10;
p=0.4;
adj_matrix = generate_adjacency_matrix(n, p);
function adj_matrix = generate_adjacency_matrix(n, p)
% Initialize an n x n matrix with all zeros
adj_matrix = zeros(n, n);
% Loop through all possible node pairs
for i = 1:n
for j = 1:n
% Skip diagonal elements (no self-loops)
if i == j
continue;
end
% Generate a random number between 0 and 1
random_number = rand();
% If the random number is less than p, add a link between nodes i and j
if random_number < p
adj_matrix(i, j) = 1;
adj_matrix(j, i) = 1; % Since it's an undirected network
end
end
end
end
disp(adj_matrix);

回答 (2 件)

Matt J
Matt J 2023 年 7 月 30 日
編集済み: Matt J 2023 年 7 月 30 日
Easier:
function adj_matrix = generate_adjacency_matrix(n, p)
adj_matrix=triu(rand(n)<p,1);
adj_matrix=adj_matrix+adj_matrix'; % Since it's an undirected network
end

Steven Lord
Steven Lord 2023 年 7 月 30 日
If you want to generate both the graph object and its adjacency matrix, tell MATLAB to build the graph using just the upper triangular part of the random matrix.
n = 10;
p = 0.6;
A = rand(n) < p;
G = graph(A, 'upper');
adj = full(adjacency(G))
adj = 10×10
0 0 1 0 0 1 0 1 0 0 0 1 0 1 1 1 0 1 1 1 1 0 1 1 0 1 0 1 0 1 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 0 1 1 1 1 1 1 1 1 1 0 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 1 1 1 0 1 1 0 0 0 1 0 1 1 1 0 0 1 1 0 1 1 1 1 0 1 0 1 1
check = triu(A)+triu(A, 1).'
check = 10×10
0 0 1 0 0 1 0 1 0 0 0 1 0 1 1 1 0 1 1 1 1 0 1 1 0 1 0 1 0 1 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 0 1 1 1 1 1 1 1 1 1 0 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 1 1 1 0 1 1 0 0 0 1 0 1 1 1 0 0 1 1 0 1 1 1 1 0 1 0 1 1
isequal(adj, check)
ans = logical
1
  2 件のコメント
mks
mks 2023 年 7 月 30 日
移動済み: Bruno Luong 2023 年 7 月 30 日
how can draw it ?
Steven Lord
Steven Lord 2023 年 7 月 30 日
Since I created the graph object before creating the adjacency matrix, just plot it.
plot(G)

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeNetworks についてさらに検索

製品


リリース

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by