代码如下
clc
clear
N=10; %设置网络节点数量
K=4; %设置邻域节点
A=zeros(N); %邻接矩阵初始化
for i=1:N %该层循环构造最近邻K耦合网络的邻接矩阵
for j=i+1:i+K/2
jj=(j<=N)*j+(j>N)*mod(j,N); %如果j超过了N,要取除以N的余数
A(i,jj)=1; A(jj,i)=1;
end
end
for i=1:N
B(i,1)=cos((2*pi/N)*(i-1));
B(i,2)=sin((2*pi/N)*(i-1)); %生成一个N*2的矩阵用以存放N个节点在图中的位置
end
gplot(A,B,"-o"); %用gplot函数画出对应邻接矩阵的图
结果如下