NW小世界网络求特征根程序

clc;
clear all;
N=input('输入节点总数N：');
K=input('输入耦合数值K:');
pp=input('输入随机化加正边的概率p+:');
pr=input('输入随机加入的抑制连接概率p-:');
T=input('输入最大运行次数T:');
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);
plot(x,y,'ro','MarkerEdgeColor','g','MarkerFaceColor','r','markersize',8);
hold on;
for t=1:T;
A=zeros(N);
for i=1:N
for j=i+1:i+K/2
jj=j;
if j>N
jj=mod(j,N);
end
A(i,jj)=1;A(jj,i)=1;
end
end %%生成邻接矩阵
[m,n]=find(A==0);
for i=1:length(m)
if m(i)>n(i) % m>n
p1=rand(1,1);
A(m(i),n(i))=1;A(n(i),m(i))=1;
end
end
r(t)=max(eig(A));
fid=fopen('f:\eigenvalue.txt','wt');
fprintf(fid,'%f\n',r(t));
end
for 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;
end
end
end
axis equal;
hold off