网络模体的同步和仿真--复杂网络的动态形成和仿真无标度网络,指数网络,ER模型,随机网络
1.2.2 数值模拟clear N=1000; A=zeros(N,1); A(1)=1;A(2)=1;A(3)=1;A(4)=1; for i=5:N; R=randperm(i);k=R(1); A(k)=A(k)+1; A(i,1)=1; end averge=sum(A)/N B=zeros(N,1); for j=1:N for i=1:N if(A(i)==j) B(j)=B(j)+1; end end end X=1:30;X=X';C=1:30; for i=1:30 C(i)=B(i); end plot(X,C) hold on
数值模拟思路: (1)运用随机排序命令 randperm(i)表示第二步中 的随机选择过程； (2) 为使图像准确，应该取 k<20%*N.左边程序运行结果
m0=4
N=1000, m=1,的数值模拟程序