n=size(C,1);%n表示问题的规模（城市个数） D=zeros(n,n);%D表示完全图的赋权邻接矩阵 for i=1:n for j=1:n if i~=j D(i,j)=((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2)^0.5; else D(i,j)=eps; %i=j时不计算，应该为0，但后面的启发因子要取倒数，用eps（浮点相对精度）表示 end D(j,i)=D(i,j); %对称矩阵 end end Eta=1./D; %Eta为启发因子，这里设
简单复杂网络算法，运行效果好，效率高。 N=input('请输入最近邻耦合网络中节点的总数N：'); %%参数输入 K=input('请输入最近邻耦合网络中每个节点的邻居数K：'); if K > floor (N-1) | mod (K,2) ~=0 disp ('参数输入错误：K值必须是小于网络节点总数且为偶数的整数'); return; end angle=0:2*pi/N:2*pi-2*pi/N;x=100*sin(angle);y=100*cos(angle);%%生成各节点坐
