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;endangle=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;     %%生成最近鄰耦合網絡的鄰接矩陣    endend[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    endend% 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% endfor 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    endendaxis 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)]);