NW小世界网络求特征根程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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