WS小世界网络的Matlab代码

复杂网络学习的经典网络之一，WS小世界网络的概念就不再赘述，直接代码开整：

function WS(N,K,p1)
%N个节点，每个节点与相邻节点共有2K条连边，断边重连概率为p1
%------------------------------------------生成规则图

A = zeros(N,N);
for i=1:N
    for j=1:K
        x=i;
        y=i+j;
        if y>N
            y=mod(y,N);
        end
        A(x,y)=1;A(y,x)=1;
    end
end
for i=1:N
    B(i,1)=cos((2*pi/N)*(i-1));
    B(i,2)=sin((2*pi/N)*(i-1));
end
figure
gplot(A,B,'-o');
title('规则图');
axis square
%-----------------------------------------开始断边重连


C=A;
for i=1:N
    for j=1:K
        x=i;
        y=i+K;
        if y>N
            y=mod(y,N);
        end
        p2=rand();
        if p2<=p1
            C(x,y)=0;C(y,x)=0;
            k=0;
            while k==0
                y=unidrnd(N);
                if C(x,y)==0 && x~=y
                    C(x,y)=1; C(y,x)=1;
                    k=1;
                end
            end
        end      
    end
end
figure
gplot(C,B,'-o')
title('WS小世界网络');
axis square



%-----------------------------------------开始计算特征系数
%Clustering_coefficient=['聚类系数：',num2str(clustering(C))];
%Average_path_length=['平均路径长度：',num2str(Alength(C))];
%disp(Clustering_coefficient);
%disp(Average_path_length);
%Degree(C);
end

最后的一些系数可以选择是否输出，具体代码后续呈现