无尺度网络

title: 源起.
content: 无尺度网络的概念是随着对复杂网络的研究而出现的。“网络”其实就是数学中图论研究的图,由一群顶点以及它们之间所连的边构成。在网络理论中则换一套说法,用“节点”代替“顶点”,用“连结”代替“边”。复杂网络的概念,是用来描述由大量节点以及这些节点之间错综复杂的联系所构成的网络。这样的网络会出现简单网络中没有的特殊拓扑特性。  自二十世纪60年代开始,对复杂网络的研究主要集中在随机网络上。随机网络,又称随机图,是指通过随机过程制造出的复杂网络。最典型的随机网络是保罗·埃尔德什和阿尔弗雷德·雷尼提出的ER模型。ER模型是基于一种“自然”的构造方法:假设有formula_1个节点,并假设每对节点之间相连的可能性都是常数formula_2是自然数formula_3的概率: formula_4  也就是说formula_5 的概率正比于formula_3 的某个幂次(一般是负的,记为formula_7)。因此formula_3越大,formula_5 的概率就越低。然而这个概率随formula_3增大而下降的“速度”是比较缓慢的:在一般的随机网络中,下降的速度是指数性的,而在无尺度网络中只是以多项式类的速度下降。  在现实中许多大规模的无尺度网络中,度分布的formula_11值介于2与3之间。在对数坐标系中,度分布将会是一条斜率介于-2至-3之间的直线。如左下图中,横坐标为节点的度,从formula_12一直到formula_13;纵坐标为找到这样的节点的概率从formula_14一直到formula_12。最高度数的节点有882条连接。所有的蓝点大致成一条直线分布(绿色的直线)。
title: 度-度相关性.
content: 仅仅是将度分布的幂律分布作为无尺度网络的定义有其不够完善之处。由于幂律分布是方差可能无穷的高可变分布,对于度分布是同一个幂律分布的不同网络,其拓扑结构和特性可能存在巨大的差异。2005年,Lun Li和大卫·阿尔德森(David Alderson)等人在论文《迈向无标度图的理论》("Towards a Theory of Scale Free Graphs")中提出了一种补充性的标度性测度。设formula_16为所有具有(依照幂律分布的)度分布formula_17的网络formula_18的集合,对于其中每一个网络formula_19,定义度-度相关数: formula_20 其中formula_21表示formula_18中所有连接的集合。根据排序原理,如果度数大的点之间相互连接的话,那么formula_23会比较大。设formula_24为最大的,那么定义度-度相关系数: formula_25 度-度相关系数formula_26介于0与1之间。formula_26越靠近1,则称此网络formula_18“无尺度”的,formula_26靠近0,则称formula_18是“富尺度”的。在此定义下,无尺度网络中的节点度数分布特征体现了自相似的性质,而凸显了“无尺度”的特征,与富尺度网络之间有相当的差异。
title: 例子.
content: 不少现实中的网络结构都属于无尺度网络,或者有无尺度的特性。以下是一些无尺度网络的例子:
title: BA模型.
content: Albert-László Barabási与Réka Albert在1999年的论文中提出了一个模型来解释复杂网络的无尺度特性,称为BA模型。这个模型基于两个假设: 在这种假设下,BA模型的具体构造为: 这样,在经过formula_44次之后,得到的新网络有formula_45个节点,一共有formula_46条边。  分析BA模型网络的渐近度分布(当节点数量很大时的度分布)主要有连续场理论、主方程法和速率方程法。这三种方法得到的渐近结果都是相同的。2001年,Béla Bollobás证明了在节点数量很大时,BA模型网络的度分布遵从formula_47的幂律分布。之后,其它的无尺度网络模型也开始被提出。
title: 其它无尺度模型.
content: BA模型成功的为无尺度网络找到了一个简单而合理的形成机制。然而,BA模型也有其自身的局限。例如,它只能描述formula_47的无尺度网络,对于真实网络的一些非幂律特征如指数截断(exponential cutoff)、小变量饱和(saturation of small variables)等无法描述。因此,各种BA模型的推广、变化版本开始出现。Bollobás在2001年提出了线性弦图模型(LCD模型),允许节点自己与自己相连。而后又出现了只允许重复连线而不允许自连线的模型与不允许重复连线、自连线而是在选中的旧节点的邻域随机联线的模型。
title: 适应度模型.
content: 在BA模型的制造过程中,人们发现,存在越久的节点具有越高的度数。然而,现实生活的网络中并非存在越久的元素就能有更多的联系。BA模型并没有包括“后起之秀”的现象。于是,出现了基于BA模型的适应度模型。适应度模型主要是修正了优先连接的机制,对每个节点加上一个吸引因子formula_49,这样新节点的相连概率改正为:
title: 局域世界演化模型.
content: 另一种基于BA模型的推广版本是局域世界演化模型。这个模型假设每个新节点在引入时并不能在全域进行优先连接。比如说一家新上市的公司可能只会与同地区或同国家的公司展开贸易联系,居民搬入新社区时只会与同一幢楼的人开始认识等等。局域世界演化模型将BA模型优先连接的机制改为:新加入的节点时,先选择全部节点的一部分(随机选取的formula_51个节点)作为局域世界,然后再在局域世界中进行优先连接。模拟结果指出,当formula_51从formula_53变化到formula_54时,产生的网络从服从指数分布逐渐过渡到服从幂律分布。
title: 复制模型.
content: 在BA模型中,度分布实际上是和增长的时间formula_55(或说增长次数)相关的,只是在formula_55十分大时近似于度分布。复制模型是一个与增长时间无关的模型。复制模型的做法是每次随机地“复制”一个原有的节点:即随机选定一个节点formula_57,再加入一个新节点,然后新节点按照formula_57与其它旧节点连接的方式与旧节点相连,最后与formula_57也相连。
title: 分层模型.
content: 2001年,Barabási提出了第一个确定性的分层网络模型。这个模型是为了解释生物学中formula_60的代谢网络。分层模型的想法是从模体(motif)出发,通过自相似的层次叠加而得到复杂网络。这种思想类似于分数维图形。Barabási的模型是:  这样形成的网络是无尺度网络,Barabási算出它的formula_61。后来有使用5节点4连结作为模体,得到formula_62,而4节点3连结作为模体得到formula_63,近似于代谢网络。需要注意的是此模型中不少度数是0概率的,所以需要使用补分布绕过。类似的确定性模型还有伪分形图(pseudofractal scale-free network)以及阿波罗模型(Apollonian network)。
title: 无尺度网络的阿喀琉斯之踵.
content: 2000年7月27日,《自然》杂志的封面文章标题是《因特网的阿喀琉斯之踵》("Achilles' Heel of the Internet")。阿喀琉斯是古希腊神话中的英雄,他出生后,他的母亲捏着他的脚踝将他浸泡在冥河中,从此他的身体刀枪不入,只有踵部没被浸到,是为其致命弱点。因此如今“阿喀琉斯之踵”常被用来称呼一个系统的致命缺陷。这篇文章中从因特网的无尺度特性出发,探讨它对意外故障的承受能力。  假设在一个网络中移除一个节点,以及与其相关的连接,那么原网络中的其他点也可能受到影响:原本相连的两个节点可能不再相连;即使相连,从其中一处到另一处可能需要经过更多的路途。总的来说,网络的连通性降低了。文章比较了ER随机网络模型与BA模型在移除少量节点时对网络连通性的影响。这个影响主要使用最大连通子图的大小formula_64与平均路径长度formula_65来衡量。在执行“随机攻击策略”,也就是在网络中随机地去除一些节点时,无尺度网络的formula_64比随机网络下降慢得多,formula_65的增长也缓慢得多。但是在执行“蓄意攻击策略”,也就是选择移除连接度最高的节点时,则会得到相反的结果。受到随机攻击的随机图会分裂成几个较小的子图,而无尺度网络则有很大概率保持连通;然而面对蓄意攻击(或称协同攻击)时,只需要移除5-10%的高于5度的节点,就能彻底瘫痪无尺度网络。
title: 流行病临界值.
content: 流行病或网络病毒在复杂网络中的传播也是复杂网络研究的方向之一。在均匀网络如ER模型随机网络或小世界网络中,如果考虑易感(S)→感染(I)→易感(S)的SIS模型,那么存在一个与网络特性相关的临界值,当有效传播率高于这个临界值的时候,传染病会在网络中传播并稳定在某个恒定密度上(激活相态)。而当有效传播率低于这个临界值时,传染病会很快逐渐消亡(吸收相态)。对于无尺度网络,由于度分布不均匀,临界值比较小。对于BA模型,临界值为0。也就是说,只要有效传播率大于0,病毒就能有效传播并达到稳定。而对于有限规模的无标度网络,临界值大于0,但会在均匀网络的十分之一左右。因此,无标度网络对于病毒传播的抵抗性较均匀网络脆弱得多。  由于无尺度网络应对流行病感染的脆弱性,人们提出不同的免疫策略来弥补。主要研究的免疫策略有三种:随机免疫、选择免疫与熟人免疫。
title: 随机免疫.
content: 随机免疫是在网络中随机抽取一部分节点进行免疫。研究表明,采取这种策略的话,需要对网络中几乎所有的节点都进行免疫才能保证最终消灭传染病。
title: 选择免疫.
content: 选择免疫是在网络中抽取度最大的节点进行免疫。就BA模型而言,采取这种策略的话,即使有效传播率变化,也可以只免疫很小一部分节点就保证消灭传染病。
title: 熟人免疫.
content: 由于选择免疫需要知道全局节点的度数情况,才能找到度数最大节点进行免疫,这在面对互联网等庞大的复杂网络时会导致难以操作。熟人免疫采取的是随机抽取一部分节点,然后对每个节点随机选一个与之相连的“邻居”节点来进行免疫。由于在无尺度网络中,度大的节点可以与非常多的节点相连,因此选择“邻居”免疫的话,碰到度大节点的概率会比碰到度小节点的概率大得多。所以熟人免疫要比随机免疫有效得多,只略差于选择免疫。
title: 同步性.
content: 音乐会或歌剧完场时,台下的观众不间断地鼓掌。在很短几次后,鼓掌的频率就会变得同步。这种现象显示出网络的同步性。研究表明,网络动力系统的同步性取决于节点动力系统的特性,节点的耦合方式与网络的结构。对于BA模型网络,节点数的增加不会降低网络同步的稳定性。而面对随机攻击和蓄意攻击,BA模型网络的同步性与连通性表现出相同的特征:对于随机攻击承受性强,而对蓄意攻击则显得脆弱。

Summary: 在网络理论中,无尺度网络(Scale-free network,或称无标度网络)是带有一类特性的复杂网络,其典型特征是在网络中的大部分节点只和很少节点连接,而有极少的节点与非常多的节点连接。这种关键的节点(称为“枢纽”或“集散节点”)的存在使得无尺度网络对意外故障有强大的承受能力,但面对协同性攻击时则显得脆弱。现实中的许多网络都带有无尺度的特性,例如因特网、金融系统网络、社会人际网络等等。