python 构建建模无标度网络

python 构建建模无标度网络用简洁的python代码生成无标度网络,并绘制网络度的分布

大家好,欢迎来到IT知识分享网。

用简洁的python代码生成无标度网络,并绘制网络度的分布

假设有如下一个简单网络

在这里插入图片描述

网络的邻接矩阵和度数如下

——- node 1 node 2 node 3 degree
node 1 ——– 1 1 2
node 2 1 —– 0 1
node 3 1 0 0 1

假设一个新的点4要和这三个点进行连接

  • 和点1、2、3连接的概率分别是: 3 3 + 2 + 2 \frac{3}{3+2+2} 3+2+23 2 3 + 2 + 2 \frac{2}{3+2+2} 3+2+22 2 3 + 2 + 2 \frac{2}{3+2+2} 3+2+22
  • 把这三个概率加在一起为1
  • 那可以使用一个随机数,看落在三段的那一段里,来判断和对应的哪个点发生连接

思路是嵌套循环+随机数+if条件判断

  1. 第一层循环,每次循环代表一个新生成的点
  2. 第二层循环(内层),每次循环代表判断是否和这个点产生连接

生成一个平均度为 2 n − 1 2 \frac{2n-1}{2} 22n1的网络,具体代码如下

import numpy as np
nodes_degree = []
N = 1000 # 生成1000个点
link = []
for i in range(N):
            link.append([])

for i in range(N): # i为每一个点,进行连接
       nodes_degree = [int(np.sum(node))+1 for node in link]
       p = np.random.randint(np.sum(nodes_degree))
       degree_control = 0
       for j, degree in enumerate(nodes_degree):
           if p < degree_control:
               link[i].append(j)
               link[j].append(i)
               break
           else:
               degree_control += degree

绘制网络度的分布

import matplotlib.pyplot as plt
mxdeg=0
for i in range(N):
    mxdeg=max(mxdeg,link[i].__len__())
mxdeg+=1

deg=[0]*mxdeg
for i in range(N):
    deg[link[i].__len__()]+=1
for i in range(mxdeg):
    deg[i]/=1.0*N
plt.plot(np.linspace(0,mxdeg-1,mxdeg),deg)
plt.xlabel('degree')
plt.ylabel('P')
plt.title('The probability distribution of degree')
plt.show()

度的分布:

在这里插入图片描述

免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://yundeesoft.com/21427.html

(0)
上一篇 2024-01-20 12:00
下一篇 2024-01-20 17:00

相关推荐

发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注

关注微信