my_matrix2.m

对一个50个结点（更多的节点的网络只需要修改模块中的标量维数就行）的复杂非线性耦合网络进行同步化仿真。首先生成K矩阵

function k=my_matrix2(n)
k=zeros(n,n);
for i=2:n-2
    k(i,i)=1;
    k(i,i+1)=-2;
    k(i,i+2)=1;
end
k(1,1)=1;
k(1,2)=-2;
k(1,3)=1;
k(n-1,1)=1;
k(n-1,n-1)=1;
k(n-1,n)=-2;
k(n,n)=1;
k(n,2)=1;
k(n,1)=-2;