复杂网络幂律分布matlab仿真程序

有关复杂网络无标度网络的幂律分布MATLAB仿真程序，希望对大家有帮助！

function Nodes=Powerlaw_with_Expected_Exponent（Nmgamanode_attribution）
%This program generates a power law degree distribution graph with expected gama
%Ref[Universal behavior of load distribution in Scale free networksPRL V87200112]
%Input:
%N--the number of nodes
%m--the mean degree =2m
%gama--the expected exponentblong to （2infinit）whose corelation with the control
%      parameter alfa is: gama=1+1/alfawhere alfa is [01）
%node_attribution--if node_attribution==1Nodes is adjacent matrix; if
%                   node_attribution=2Nodes is adjacent list.
%Output
%Nodes--return the adjacent matrix if node_attribution=1.the adjacent list
%       if node_attribution=2
%



if node_attribution==1
    Nodes=zeros（NN）;
    p=zeros（N1）;
    pp=zeros（N1）;
    node_1_Len=zeros（N1）;
    node_2_Len=zeros（N1）;
    alfa=1/（gama-1）;
    edge_num=N*m;
    
    %1.get nodes‘ weight
    for i=1:N
        p（i）=1/（i^alfa）;        
    end
    %2.normalized node weights
    sum_p=sum（p）;
    pp（1）=p（1）/sum_p;
    for i=2:N
        pp（i）=pp（i-1）+p（i）/sum_p;        
    end
    
    %3.add edge
    for i=1:edge_num
       %i
       ADD_ONE_EDGE=0;
       while ADD_ONE_EDGE==0
            node_1_Len=find（pp>rand（1））; 
            node_1=node_1_Len（1）;
            node_2_Len=find（pp>rand（1））;
            node_2=node_2_Len（1）;
            while node_2==node_1%avoid self-loop
                node_2_Len=find（pp>rand（1））;
                node_2=node_2_Len（1）;
            end
            if Nodes（node_1node_2）==0
                Nodes（node_1node_2）=1;
                Nodes（node_2node_1）=1;
                ADD_ONE_EDGE=1;
            else
                ADD_ONE_EDGE=0;
            end%Nodes（node_1node_2）
       end%while ADD_ONE_EDGE==0
    end% for i=1:edge_num    
end%node_attribution


if TEST==1
    t=toc
    for i=1:N
        Degree（i）=nnz（Nodes（i:））;
    end
    plot_distribution（DegreeLine_color1）;
    fname=[‘PLRG‘num2str（N）‘k‘num2str（m）‘gama‘num2str（gama）];
    Write_Sparse_Matrix（Nodes[fname‘.adj‘]）;
    Write_into_Pajek（Nodes[fname‘.net‘]）;
end

return

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2180  2008-07-24 23:15  Powerlaw_with_Expected_Exponent.m

----------- ---------  ---------- -----  ----

                 2180                    1