python实现小世界网络生成
作者:Alvassss 时间:2022-07-22 07:54:46
没有使用igraph库哦 因为我还没学
小世界网络简介:
1998年, Watts和Strogatz 提出了小世界网络这一概念,并建立了WS模型。实证结果表明,大多数的真实网络都具有小世界特性(较小的最短路径)和聚类特性(较大的聚类系数)。传统的规则最近邻耦合网络具有高聚类的特性,但并不具有小世界特性;而随机网络具有小世界特性但却没有高聚类特性。因此这两种传统的网络模型都不能很好的来表示实际的真实网络。Watts和Strogatz建立的小世界网络模型就介于这两种网络之间,同时具有小世界特性和聚类特性,可以很好的来表示真实网络。
小世界模型构造算法
1、从规则图开始:考虑一个含有N个点的最近邻耦合网络,它们围成一个环,其中每个节点都与它左右相邻的各K/2节点相连,K是偶数。
2、随机化重连:以概率p随机地从新连接网络中的每个边,即将边的一个端点保持不变,而另一个端点取为网络中随机选择的一个节点。其中规定,任意两个不同的节点之间至多只能有一条边,并且每一个节点都不能有边与自身相连。
在上述模型中,p=0对应于完全规则网络,p=1则对应于完全随机网络,通过调节p的值就可以控制从完全规则网络到完全随机网络的过渡。
效果如下:
代码如下:
import matplotlib.pyplot as plt
import random as rd
import numpy as np
plt.rcParams['font.sans-serif'] = ['SimHei'] # 中文字体设置
plt.rcParams['axes.unicode_minus']=False
#小世界项目
def dian(N,K,P):
global ls
tim=[]
for i in range(N):
for j in range(1,K+1):
ls[i]=ls.get(i,set())
ls[i].add((i+j)%N)
ls[i].add((i-j)%N)
ls[(i-j)%N]=ls.get((i-j)%N,set())
ls[(i-j)%N].add(i)
ls[(i+j)%N]=ls.get((i+j)%N,set())
ls[(i+j)%N].add(i)
for i in range(N):
for j in list(ls[i]):
if rd.random()<=P:
aa=ls[i].pop()
a=set(range(N))
a.discard(i)
a=a^ls[i]
for i in range(rd.randint(1,len(a)-1)):
aa=a.pop()
ls[aa].discard(i)
b=a.pop()
ls[i].add(b)
ls[b].add(i)
for i in range(N):
tim.append(len(ls[i])*40-N)
new=[]
for i in range(len(ls)):
l=[]
l.append(i)
l+=list(ls[i])
new.append(l)
return new,tim
def hua(L,S):
x=np.linspace(0,100,len(L))
y=np.sqrt(np.abs(10000-(x-50)**2))
plt.scatter(x,y,s=S,edgecolor='k',alpha=0.7)
for i in range(len(L)):
plt.text(x[i]-0.13,y[i]-0.015,str(S[i]//40+1))
for j in L[i]:
plt.plot(list((x[i],x[j])),list((y[i],y[j]))\
,color='gray',linewidth=1,alpha=0.7)
plt.title('小世界网络初步')
plt.xticks([])
plt.yticks([])
plt.axis('off')
plt.savefig('niu.png')
ls={}
l,k=dian(20,3,0.5) #不要超过40哦~
hua(l,k)
来源:https://blog.csdn.net/Alvassss/article/details/86617594
标签:python,小世界网络,生成
0
投稿
猜你喜欢
SQL Server 2005改进后的几个实用新特性
2008-05-07 19:16:00
mysql通过文档读取并执行命令之快速为mysql添加多用户和数据库技巧
2024-01-24 01:45:02
《JavaScript语言精粹》
2009-04-03 11:27:00
一文带你掌握Python中多线程和线程池的使用方法
2022-10-20 21:53:09
Python序列化模块之pickle与json详解
2023-07-08 05:48:23
Python进程间通信之共享内存详解
2022-01-06 05:47:10
Django中使用Whoosh进行全文检索的方法
2023-03-20 05:35:43
sql server中批量插入与更新两种解决方案分享(存储过程)
2012-05-22 18:29:59
Python正则表达式和元字符详解
2021-02-26 02:54:19
Java基于redis和mysql实现简单的秒杀(附demo)
2024-01-16 16:55:13
python随机生成大小写字母数字混合密码(仅20行代码)
2021-04-03 23:28:47
keras获得某一层或者某层权重的输出实例
2023-04-06 10:54:27
基于pycharm的beautifulsoup4库使用方法教程
2022-12-05 22:47:05
Sql学习第三天——SQL 关于CTE(公用表达式)的递归查询使用
2024-01-23 01:00:10
python+OpenCV人脸识别考勤系统实现的详细代码
2021-11-01 01:28:44
Python面向对象特殊成员
2021-07-03 12:05:48
MySQL备份与恢复之保证数据一致性(5)
2024-01-16 19:16:20
一篇文章介绍redux、react-redux、redux-saga总结
2023-08-22 16:56:32
python绘制汉诺塔
2023-07-09 23:13:19
python开发之str.format()用法实例分析
2022-01-15 22:10:23