生物医学网络科学Biomedical Network Science
网络结构指标与随机图12 / 19

ER 随机网络

阅读
约 2 分钟

模型定义

Erdős–Rényi 随机网络:G(N,p)G(N,p) 模型中,NN 个节点的每一对以概率 pp 独立连边;G(N,M)G(N,M) 模型中,从所有可能的边中均匀随机地选出 MM 条。NN 很大时两者性质相同。(G(N,p)G(N,p) 也常归于 Gilbert 1959。)

注意 NN 和 pp 并不唯一决定一张图:同样的参数,每次生成的网络都不同。我们研究的是这一类网络的统计性质,以及当 N→∞N\to\infty 时“几乎必然”成立的性质。

边数:(N2)\binom N2 个节点对各以概率 pp 连边,故边数服从二项分布,期望 ⟨M⟩=pN(N−1)/2\langle M\rangle=p\,N(N-1)/2,平均度 ⟨k⟩=2⟨M⟩/N=p(N−1)\langle k\rangle=2\langle M\rangle/N=p(N-1)。

度分布

每个节点与其余 N−1N-1 个节点各以概率 pp 相连,因此度服从二项分布,在 N→∞N\to\infty、⟨k⟩\langle k\rangle 固定时趋于泊松分布

P(k)=(N−1k)pk(1−p)N−1−k→N→∞,⟨k⟩固定e−⟨k⟩⟨k⟩kk!,P(k) = \binom{N-1}{k} p^k (1-p)^{N-1-k} \;\xrightarrow{N\to\infty,\ \langle k\rangle \text{固定}}\; e^{-\langle k\rangle}\frac{\langle k\rangle^k}{k!},

其中 ⟨k⟩=p(N−1)\langle k\rangle = p(N-1)。泊松分布的方差等于均值,标准差为 ⟨k⟩\sqrt{\langle k\rangle},相对宽度 σk/⟨k⟩=1/⟨k⟩\sigma_k/\langle k\rangle=1/\sqrt{\langle k\rangle}:平均度越大,节点的度越集中在平均值附近,几乎不存在枢纽。

最大度与最小度

要求“度不小于 kmaxk_{\max} 的节点期望只有一个”:NP(K≥kmax)≈1N\,P(K\ge k_{\max})\approx1;同理 NP(K≤kmin)≈1N\,P(K\le k_{\min})\approx1。Barabási 书中以 N=109N=10^9、⟨k⟩=1000\langle k\rangle=1000(粗略对应全球熟人网络)为例:数值求解得 kmaxk_{\max} 约 1.2×1031.2\times10^3(书中用级数近似得 1185)、kmin≈816k_{\min}\approx816。也就是说,即使在十亿规模的随机社会中,“最受欢迎的人”的朋友也只比平均多约 20%。真实社交网络中存在朋友数比平均多几个数量级的人,这本身就否定了随机网络模型。

聚集系数

节点 ii 的任意两个邻居之间是否连边,与其他边无关,概率仍为 pp,所以 Crand≈p=⟨k⟩/(N−1)C_{\text{rand}} \approx p = \langle k\rangle/(N-1)。在平均度固定时,CrandC_{\text{rand}} 随 NN 增大按 1/N1/N 趋于 0,并且 C(k)C(k) 与 kk 无关。真实网络的平均聚集系数则大体与规模无关(Barabási 书图 3.13)。

平均路径长度

Lrand≈ln⁡N/ln⁡⟨k⟩L_{\text{rand}} \approx \ln N/\ln \langle k\rangle,只随 NN 对数增长:即便有几十亿个节点,平均路径也只有几步到十几步。推导与修正见小世界性与平均距离。

巨连通分支

逐渐增大 pp,ER 网络会经历一个相变(phase transition):⟨k⟩<1\langle k\rangle<1 时只有小碎片,⟨k⟩=1\langle k\rangle=1 时巨连通分支出现,⟨k⟩>ln⁡N\langle k\rangle>\ln N 时整个网络几乎必然连通。巨分支所占比例 SS 满足 S=1−e−⟨k⟩SS = 1 - e^{-\langle k\rangle S}。四个区间的细节见ER 网络的演化与相变,各结果的推导见数学推导与算法。本专题页面的 ER 演示可以拖动 pp 观察这一过程。