%生成小世界網(wǎng)絡(luò)的鄰接矩陣
tic %程序運(yùn)行時(shí)間控制
%設(shè)置網(wǎng)絡(luò)可調(diào)控制參數(shù)
N=input('輸入網(wǎng)絡(luò)節(jié)點(diǎn)個(gè)數(shù)： ');
m=input('輸入節(jié)點(diǎn)的最小度m： ');
p=input('輸入概率p,其中（1>=p>=0）： ');
%N=10;m=3;
%p=0.1;%平均度:2m,m<<N
%初始的環(huán)形規(guī)則格點(diǎn)網(wǎng)絡(luò)
matrix=sparse([]);
for i=m+1:N- m
   for j=i- m:i+m
        matrix(i,j)=1;
   end
end
for i=1:m
     for j=1:i+m
            matrix(i,j)=1;
     end
end
for i=N- m+1:N
      for j=i- m:N
           matrix(i,j)=1;
      end
end
for i=1:m
      for j=N- m+i:N
           matrix(i,j)=1;matrix(j,i)=1;
      end
end
%逆時(shí)針的邊重連,從節(jié)點(diǎn)1 到N- m- 1
for i=1:N- m- 1
     for j=i+1:i+m
         r=rand(1);
         if r<=p
            unconect=find(matrix(i,:)==0);
            r1=ceil(length(unconect)*rand(1));
            matrix(i,unconect(r1))=1;
            matrix(unconect(r1),i)=1;
            matrix(i,j)=0; matrix(j,i)=0;
         end
     end
end
%逆時(shí)針的邊重新連接,從節(jié)點(diǎn)N- m到N- 1
for i=N- m+1:N- 1
       for j=[i+1:N 1:i- N+m]
           r=rand(1);
           if r<=p
               unconect=find(matrix(i,:)==0);
               r1=ceil(length(unconect)*rand(1));
               matrix(i,unconect(r1))=1;
               matrix(unconect(r1),i)=1;
               matrix(i,j)=0;matrix(j,i)=0;
           end
       end
end
%逆時(shí)針的邊重新連接,節(jié)點(diǎn)N
for i=N
      for j=1:m
          r=rand(1);
          if r<=p
             unconect=find(matrix(i,:)==0);
             r1=ceil(length(unconect)*rand(1));
             matrix(i,unconect(r1))=1;
             matrix(unconect(r1),i)=1;
             matrix(i,j)=0;matrix(j,i)=0;
          end
     end
end
%恢復(fù)小世界網(wǎng)絡(luò)的鄰接矩陣
for m=1:N
     matrix(m,m)=0;
end
%存儲(chǔ)鄰接矩陣
save data matrix;
toc %程序結(jié)束
r=0:2*pi/N:2*pi;
x=cos(r);
y=sin(r);
plot(x,y,'o'),
hold on,
for i=1:N
    for j=i:N
       if matrix(i,j)==1
            n=[cos(i*2*pi/N) cos(j*2*pi/N)];
             m=[sin(i*2*pi/N) sin(j*2*pi/N)];
             plot(n,m,'b-');
            
        end
    end
end