生物医学网络科学Biomedical Network Science
小世界与无标度网络模型06 / 17

Barabási–Albert 优先连接模型

阅读
约 3 分钟

两个要素:增长与优先连接

Barabási 和 Albert 注意到 ER 和 WS 模型都假设节点数固定,而真实网络大多在不断增长:网页不断被创建,论文不断被发表,基因在演化中不断复制。他们提出的模型包含两个要素:

  1. 增长(growth):从 m0m_0 个节点开始,每个时间步加入一个新节点,新节点带 mm 条边(m≤m0m \le m_0)连到已有节点;
  2. 优先连接(preferential attachment):新节点连向已有节点 ii 的概率与 ii 当前的度成正比:

Π(ki)=ki∑jkj.\Pi(k_i) = \frac{k_i}{\sum_j k_j}.

直观上这就是“富者愈富”:被引用多的论文更容易被看到、被继续引用。这一思想可以追溯到 Yule、Simon 以及 Price(1965、1976)对引文网络的研究。Price 把它称为“累积优势”(cumulative advantage),并提出了有向的引文生长模型:新论文引用旧论文的概率正比于旧论文已有引用数加一个常数(否则没被引用过的论文永远不会被引用)。BA 模型可以看作 Price 模型的无向版本,常数项取零。

主要结论

平均场分析(推导见数学推导与算法)给出两个结论:

  • 每个节点的度按 ki(t)=m(t/ti)1/2k_i(t) = m\,(t/t_i)^{1/2} 增长,越早加入的节点(tit_i 越小)度越大——先到者优势;
  • 度分布为 P(k)≈2m2/k3P(k) \approx 2m^2/k^3,即 γ=3\gamma = 3 的幂律,而且指数与 mm 无关。用速率方程可得精确解 P(k)=2m(m+1)k(k+1)(k+2)P(k) = \dfrac{2m(m+1)}{k(k+1)(k+2)},大 kk 时与平均场结果一致(Dorogovtsev、Mendes 与 Samukhin 2000;Krapivsky、Redner 与 Leyvraz 2000;严格证明见 Bollobás 等 2001)。

由 ki(t)=m(t/ti)1/2k_i(t)=m(t/t_i)^{1/2} 还能看出:最大的枢纽是最早的节点,其度 kmax∼mt1/2=mN1/2k_{\max}\sim m\,t^{1/2}=mN^{1/2},正好等于 γ=3\gamma=3 时的自然截断 N1/(γ−1)N^{1/(\gamma-1)}。

本专题页面的 BA 演示可以逐步加点,在双对数图上观察度分布的形成。

两个要素缺一不可

  • 只有增长、没有优先连接(Barabási 书中的“模型 A”:新节点均匀随机地连接已有节点):dki/dt=m/tdk_i/dt=m/t,度只按对数增长 ki=m+mln⁡(t/ti)k_i=m+m\ln(t/t_i),度分布变为指数分布 pk=eme−k/mp_k=\dfrac{e}{m}e^{-k/m},没有幂律尾巴;
  • 只有优先连接、没有增长(“模型 B”:节点数 NN 固定,每步按优先连接加一条边):所有节点的度都近似线性增长 ki≈2t/Nk_i\approx 2t/N,度分布先短暂呈长尾,随后变成以平均度为中心的窄峰,最后趋于完全图——幂律只是暂态。

距离与聚集

  • 平均距离:m≥2m\ge2 时 BA 网络的平均距离按 ln⁡N/ln⁡ln⁡N\ln N/\ln\ln N 增长(Bollobás 与 Riordan 2004;Cohen 与 Havlin 2003),比同规模随机网络还短;m=1m=1 时网络是一棵树,距离按 ln⁡N\ln N 增长。
  • 聚集系数:⟨C⟩∝(ln⁡N)2/N\langle C\rangle\propto (\ln N)^2/N(Fronczak 等 2003),比随机网络的 1/N1/N 下降得慢,但随规模增大仍趋于 0,远低于多数真实网络。

BA 模型的局限与扩展

BA 模型只能产生 γ=3\gamma=3,而真实网络的指数差异很大;它的聚集系数随 NN 增大而趋于 0,与真实网络的高聚类不符;它预测“老节点度大”,但现实中后来者常能反超(Google 晚于许多搜索引擎出现);它也没有说明优先连接本身从何而来。后续的工作沿几个方向展开,本专题依次介绍: