小世界网络matlab程序,小世界网络简介及及matlab建模.doc

小世界网络简介及及matlab建模.doc

小世界网络MATLAB建模1简介小世界网络存在于数学、物理学和社会学中，是一种数学图的模型。在这种图中大部份的结点不与彼此邻接，但大部份结点可以通过任一其它节点经少数几步就可以产生联系。若将一个小世界网络中的点代表一个人，而联机代表人与人之间是相互认识的，则这小世界网络可以反映陌生人通过彼此共同认识的人而起来产生联系关系的小世界现象。在日常生活中，有时你会发现，某些你觉得与你隔得很“遥远”的人，其实与你“很近”。小世界网络就是对这种现象的数学描述。用数学中图论的语言来说，小世界网络就是一个由大量顶点构成的图，其中任意两点之间的平均路径长度比顶点数量小得多。除了社会人际网络以外，小世界网络的例子在生物学、物理学、计算机科学等领域也有出现。许多经验中的图可以用小世界网络来作为模型。因特网、公路交通网、神经网络都呈现小世界网络的特征。小世界网络最早是由邓肯瓦茨(DUNCANWATTS)和斯蒂文斯特罗加茨(STEVENSTROGATZ)在1998年引进的，将高聚合系数和低平均路径长度作为特征，提出了一种新的网络模型，一般就称作瓦茨斯特罗加茨模型(WS模型)，这也是最典型的小世界网络的模型。由于WS小世界模型构造算法中的随机化过程有可能破坏网络的连通性，纽曼NEWMAN和瓦茨WATTS提出了NW小世界网络模型，该模型是通过用“随机化加边”模式来取代WS小世界网络模型构造中的“随机化重连”。在考虑网络特征的时候，使用两个特征来衡量网络特征路径长度和聚合系数。特征路径长度(CHARACTERISTICPATHLENGTH)在网络中，任选两个节点，连同这两个节点的最少边数，定义为这两个节点的路径长度，网络中所有节点对的路径长度的平均值，定义为网络的特征路径长度。这是网络的全局特征。聚合系数CLUSTERINGCOEFFICIENT假设某个节点有K个边，则这K条边连接的节点之间最多可能存在的边的个数为KK1/2，用实际存在的边数除以最多可能存在的边数得到的分数值，定义为这个节点的聚合系数。所有节点的聚合系数的均值定义为网络的聚合系数。聚合系数是网络的局部特征，反映了相邻两个人之间朋友圈子的重合度，即该节点的朋友之间也是朋友的程度。我们可以发现规则网络具有很高的聚合系数，大世界(LARGEWORLD，意思是特征路径长度很大)，其特征路径长度随着N网络中节点的数量线性增长，而随机网络聚合系数很小，小世界(SMALLWORLD，意思是特征路径长度小)，其特征路径长度随着LOGN增长中说明，在从规则网络向随机网络转换的过程中，实际上特征路径长度和聚合系数都会下降，到变成随机网络的时候，减少到最少。但这并不是说大的聚合系数一定伴随着大的路径长度，而小的路径长度伴随着小的聚合系数，小世界网络就具有大的聚合系数，而特征路径长度很小。试验表明，少量的SHORTCUT的建立能够迅速减少特征路径长度，而聚合系数变化却不大，因为某一个SHORTCUT的建立，不仅影响到所连接的节点的特征路径长度，而且影响到他们邻居的路径长度，而对整个网络的聚合系数影响不大。这样，少量的SHORTCUT的建立就能使整个网络不知不觉地变成小世界网络。实际的社会、生态、等网络都是小世界网络，在这样的系统里，信息传递速度快，并且少量改变几个连接，就可以剧烈地改变网络的性能，如对已存在的网络进行调整，如蜂窝电话网，改动很少几条线路，就可以显著提高性能。2小世界网络构成原则WS小世界网络的构成原则为从一个环状的规则网络开始，网络含有N个结点，每个结点向与它最近邻的K个结点连出K条边，并满足NKINN1。随后进行随机化重连，以概率P随机地重新连接网络中的每个边，即将边的一个端点保持不变，而另一个端点取为网络中随机选择的一个节点。其中规定，任意两个不同的节点之间至多只能有一条边，并且每一个节点都不能有边与自身相连。这样就会产生PNK/2条长程的边把一个结点和远处的结点联系起来。改变P值可以实现从规则网络P0向随机网络P1转变。NW小世界网络的构成原则为从一个环状的规则网络开始，网络含有N个结点，每个结点向与它最近邻的K个结点连出K条边，并满足NKINN1。随后进行随机化加边，以概率P在随机选取的一对节点之间加上一条边。其中，任意两个不同的节点之间至多只能有一条边，并且每一个节点都不能有边与自身相连。改变P值可以实现从最近邻耦合网络P0向全局耦合网络P1转变。在P足够小和N足够大时，NW小世界模型本质上等同于WS小世界模型。3MATLAB建模建立一个初始节点数为20的NW网络。MATLAB程序如下FUNCTIONMATRIXSWBY201121250314TICN20M4初始化网络数据P01以概率P01在随机选取的一对结点之间加上一条边MATRIXSPARSE,,,20,20,0创建一个2020的全0稀疏矩阵建立初始的环状的规则网络结点网络有N个节点每个结点向与它最近邻的M个结点连出边求出邻接矩阵FORIM1NMFORJIMIMMATRIXI,J1ENDENDFORI1MFORJ1IMMATRIXI,J1ENDENDFORINM1NFORJIMNMATRIXI,J1ENDENDFORI1MFORJNMINMATRIXI,J1MATRIXJ,I1ENDEND逆时针的边重连,从节点到NM1FORI1NM1FORJI1IMRRAND1随机选取一个数IFRPUNCONECTFINDMATRIXI,0取出邻接矩阵中的非0元素位置MLENGTHUNCONECT求出非0元素个数R1CEILMRAND1正向取整MATRIXI,UNCONECTR11MATRIXUNCONECTR1,I1连接这一对点MATRIXI,J0MATRIXJ,I0加上这个是SW小世界网络ENDENDEND逆时针的边重新连接,从节点NM到N1FORINM1N1FORJI1N1INMRRAND1IFRPUNCONECTFINDMATRIXI,0R1CEILLENGTHUNCONECTRAND1MATRIXI,UNCONECTR11MATRIXUNCONECTR1,I1MATRIXI,J0MATRIXJ,I0ENDENDEND逆时针的边重新连接,节点NFORINFORJ1MRRAND1IFRPUNCONECTFINDMATRIXI,0R1CEILLENGTHUNCONECTRAND1MATRIXI,UNCONECTR11MATRIXUNCONECTR1,I1MATRIXI,J0MATRIXJ,I0ENDENDEND恢复小世界网络的邻接矩阵1NMATRIXM,M0去掉自身节点形成的环END存储邻接矩阵SAVEDATAMATRIXTOC计算程序耗时END上述程序建立了一个NW小世界网络，求出了其邻接矩阵，用TU_PLOT函数画出邻接矩阵的图，就得出了该小世界网络的图形。FUNCTIONTU_PLOTREL,CONTROL由邻接矩阵画连接图，输入为邻接矩阵REL，必须为方阵；CONTROL为控制量，0表示画出的图为无向图，1表示有向图。默认值为0R_SIZESIZERELASIZEX返回的是一个行向量，该行向量第一个元素是X的行数，第2个元素是X的列数IFNARGIN2NARGIN是用来判断输入变量个数的函数