最经典的小世界网络WS的matlab源程序

function A=WS_net（）
%%% 从有N个节点，每个节点有2K个邻居节点的最近邻耦合网络图通过随机化重连生成WS小世界网路
%% A ——————返回生成网络的邻接矩阵
disp（‘该程序生成WS小世界网路：‘）;
N=input（‘请输入最近邻耦合网络中节点的总数N：‘）;
K=input（‘请输入最近邻耦合网络中每个节点的邻居节点的个数的一半K：‘）;
p=input（‘请输入随机化重连的概率p:‘）;
if K>floor（N/2）
    disp（‘输入的K值不合法‘）
    return;
end
angle=0:2*pi/N:2*pi-2*pi/N;  %%生成最近邻耦合网络的各节点坐标
x=100*sin（angle）;
y=100*cos（angle）;
plot（xy‘ro‘‘MarkerEdgeColor‘‘g‘‘MarkerFaceColor‘‘r‘‘markersize‘8）;
hold on; 

A=zeros（N）;
for i=1:N
    for j=i+1:i+K
        jj=j;
        if j>N
            jj=mod（jN）;
        end
      A（ijj）=1; A（jji）=1;     %%生成最近邻耦合网络的邻接矩阵
    end
end

for i=1:N
    for j=i+1:i+K
        jj=j;
        if j>N
            jj=mod（jN）;
        end
        p1=rand（11）;
        if p1