注意
Go to the end 下载完整示例代码。
埃尔德什-雷尼图
这个示例演示了如何使用igraph.GraphBase.Erdos_Renyi()生成Erdős–Rényi图。图有两种变体:
Erdos_Renyi(n, p)将从所谓的 \(G(n,p)\) 模型中生成一个图,其中任意两个节点之间的每条边都有独立的概率p存在。Erdos_Renyi(n, m)将从所有具有n个节点和m条边的图中随机均匀地选择一个图。这被称为 \(G(n,m)\) 模型。
我们为每个生成两个图,以便确认我们的图生成器是真正随机的。
import igraph as ig
import matplotlib.pyplot as plt
import random
首先,我们设置一个随机种子以确保可重复性
random.seed(0)
然后,我们生成两个具有相同参数的\(G(n,p)\) Erdős–Rényi图:
g1 = ig.Graph.Erdos_Renyi(n=15, p=0.2, directed=False, loops=False)
g2 = ig.Graph.Erdos_Renyi(n=15, p=0.2, directed=False, loops=False)
为了比较,我们还生成了两个具有固定边数的\(G(n,m)\) Erdős–Rényi图:
g3 = ig.Graph.Erdos_Renyi(n=20, m=35, directed=False, loops=False)
g4 = ig.Graph.Erdos_Renyi(n=20, m=35, directed=False, loops=False)
我们可以打印出每个图的摘要以验证它们的随机性
ig.summary(g1)
ig.summary(g2)
ig.summary(g3)
ig.summary(g4)
# IGRAPH U--- 15 18 --
# IGRAPH U--- 15 21 --
# IGRAPH U--- 20 35 --
# IGRAPH U--- 20 35 --
IGRAPH U--- 15 23 --
IGRAPH U--- 15 28 --
IGRAPH U--- 20 35 --
IGRAPH U--- 20 35 --
最后,我们可以绘制图表来说明它们的结构和差异:
fig, axs = plt.subplots(2, 2)
# Probability
ig.plot(
g1,
target=axs[0, 0],
layout="circle",
vertex_color="lightblue"
)
ig.plot(
g2,
target=axs[0, 1],
layout="circle",
vertex_color="lightblue"
)
axs[0, 0].set_ylabel('Probability')
# N edges
ig.plot(
g3,
target=axs[1, 0],
layout="circle",
vertex_color="lightblue",
vertex_size=15
)
ig.plot(
g4,
target=axs[1, 1],
layout="circle",
vertex_color="lightblue",
vertex_size=15
)
axs[1, 0].set_ylabel('N. edges')
plt.show()

脚本的总运行时间: (0 分钟 2.056 秒)