nw_net.m

复杂网络中有一些常见的基本模型

function A=NW_net()
%%% 从有N个节点，每个节点有2K个邻居节点的最近邻耦合网络图通过随机化重连生成NW小世界网路
%% A ——————返回生成网络的邻接矩阵
disp('该程序生成NW小世界网路：');
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(x,y,'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(j,N);
        end
      A(i,jj)=1; A(jj,i)=1;     %%生成最近邻耦合网络的邻接矩阵
    end
end
[m,n]=find(A==0);          %找出所有无边相连的节点对，注意A(i,i)=0,要排除。
for i=1:length(m)
    if m(i)~=n(i)
        p1=rand(1,1);
        if p>p1
            A(m(i),n(i))=1;  %加边策略：对未加边的两节点，生成的随机数小于p时，则加边
            A(n(i),m(i))=1;
        end
    end
end

% for i=1:N
%     for j=i+1:i+K
%         jj=j;
%         if j>N
%             jj=mod(j,N);         %%%有问题
%         end
%         p1=rand(1,1);
%         if p1<p              %% 生成的随机数小于p，则边进行随机化重连,否则，不进行加边
%           %  A(i,jj)=0;A(jj,i)=0;  
%             A(i,i)=inf; a=find(A(i,:)==0);     %重连策略：不断开原来的边，直接在未连的边中随机选择另一个节点，
%             m=length(a);                       %      并与原节点连接。
%             if m==0
%                 disp(['节点',int2str(i),'已与其它所有节点相连，不可再对该节点加边']);
%                 A(i,i)=0;
%                 break;              
%             else                
%             rand_data=randint(1,1,[1,length(a)]);  
%             jjj=a(rand_data);
%             A(i,jjj)=1;A(jjj,i)=1;
%             A(i,i)=0;
%             end
%         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;      %% 画出NW小世界网络图
        end
    end
end
axis equal;
hold off  

[C,aver_C]=Clustering_Coefficient(A);
[DeD,aver_DeD]=Degree_Distribution(A);
[D,aver_D]=Aver_Path_Length(A);   
 disp(['该随机图的平均路径长度为：',num2str(aver_D)]);  %%输出该网络的特征参数
 disp(['该随机图的聚类系数为：',num2str(aver_C)]);
 disp(['该随机图的平均度为：',num2str(aver_DeD)]);