ba.m

一种构建BA无标度复杂网络模型的matlab代码

/一种随机生成BA的模型/
NO=20;
WO=[];
N=95;
W=zeros(N,N);
%for i=1:N
%   W(i,i)=0;
%end
for i=1:N-1
    for j=(i+1):N
        if(rand>0.899)
           W(i,j)=floor(rand*10);
           if(W(i,j)==9)
               if(rand<0.99)
                   W(i,j)=3;
               end
           end
           if(W(i,j)==8)
               if(rand<0.93)
                   W(i,j)=4;
               end
           end
           if(W(i,j)==7)
               if(rand<0.91)
                   W(i,j)=2;
               end
           end
            if(W(i,j)==6)
               if(rand<0.90)
                   W(i,j)=1;
               end
           end
       end
    end
end
for i=2:N
    for j=1:i-1
        W(i,j)=W(j,i);
    end
end
WO=W;

%以上是生成50节点的随机网络,下面是增加节点NO-N
m=6;

for i=1:NO-N
   m_flag=1; 
    sum=0;
    for j=1:length(WO)-1
       for k=j+1:length(WO)
           if WO(j,k)~=0
               sum=sum+1;
          end
       end
    end
    for j=1:length(WO)
        sumj=0;
        for k=1:length(WO)
            sumj=sumj+WO(j,k);
        end
        temp(j)=sumj/sum;
    end
           
     s=zeros(1,length(WO));
     l=1;
     while l<=length(WO) & m_flag<length(WO)/15
        if(rand<=temp(l))
          Q=floor(rand*10);
           if(Q~=0)
               m_flag=m_flag+1;  
           end
         s(l)=Q; 
        end
         l=l+1;
     end
     while m_flag<=m
         T=floor(rand*length(WO));
         if T>0    
            Q=floor(rand*10);
            if(Q~=0)
               m_flag=m_flag+1;
            end
           s(T)=Q; 
         end
     end
     WO=[WO;s];
     WO=[WO,[s';0]];
         
end
       
Dogx=fopen('D:\MATLAB7\work\network.doc','w');
for i=1:NO
    for j=1:NO
       Dogy=fprintf(Dogx,'%d',WO(i,j));
       Dogy=fprintf(Dogx,',');
      end
        Dogy=fprintf(Dogx,'\n');
end
fclose(Dogx);