WS小世界网络模型的程序代码(matlab)

程序仿真实例

例一、请输入最近邻耦合网络中节点的总数N:30 请输入最近邻耦合网络中每个节点的邻居数K:4

请输入随机化重连的概率p:0.9

例二、请输入最近邻耦合网络中节点的总数N:40 请输入最近邻耦合网络中每个节点的邻居数K:2

请输入随机化重连的概率p:0.7

Matlab的m文件代码如下：

N=input('请输入最近邻耦合网络中节点的总数N:');

K=input('请输入最近邻耦合网络中每个节点的邻居数K:');

if K>floor(N-1)|mod(K,2)~=0;

disp('参数输入错误：K值必须是小于网络节点总数且为偶数的整数'); return ;

end

angle=0:2*pi./N:2*pi-2*pi/N;

angle=0:2*pi/N:2*pi-2*pi/N;

x=100*sin(angle);

y=100*cos(angle);

plot(x,y,'ro','MarkerEdgeColor','g','MarkerFaceColor','r','MarkerSize',8);

hold on;

A=zeros(N);

for i=1:N

for j=i+1:i+K/2

jj=j;

if j>N

jj=mod(j,N);

end

A(i,jj)=1;

A(jj,i)=1;

end

end

%WS小世界网络的代码

p=input('请输入随机化重连的概率p:');

for i=1:N

for j=i+1:i+K/2

jj=j;

if j>N

jj=mod(j,N);

end

p1=rand(1,1);

if p1

A(i,jj)=0;A(jj,i)=0;A(i,i)=inf;a=find(A(i,:)==0);

rand_data=randi([1,length(a)],1,1);

jjj=a(rand_data);A(i,jjj)=1;A(jjj,i)=1;A(i,i)=0;

end

end

end

for i=1:N

for j=i+1:N

if A(i,j)~=0

plot([x(i),x(j)],[y(i),y(j)],'linewidth',1.2);

hold on;

end

end

end

axis equal;

hold off