小编Ign*_*iña的帖子

第一次在这里提问。我以前使用一个简单的 MATLAB 脚本对 90 个 Hopf 振荡器进行建模，这些振荡器通过矩阵耦合，具有 randn 噪声，并具有简单的欧拉阶跃积分。我想升级这个，所以我进入了 Julia，似乎有很多令人兴奋的特性。

这是我正在求解的方程组

我有点迷失了。我开始使用 Differentialequations.jl (stochastic solver) ，找到了一个解决方案，并发现自己有一个基准，告诉我解决 200 秒占用 4 Gb ！（2.5 Gb with alg_hints=[:stiff]）（我没有修复 dt，以前我使用 dt=0.1）

function Shopf(du,u,p,t)

  du[1:90,1]=(p[1:90,1]-u[1:90,1].^2.0-u[1:90,2].^2.0).*u[1:90,1]-p[1:90,2].*u[1:90,2] + 0.5*(-p[: , end].*u[:,1]+p[:,4:end-1] *u[:,1])
  du[1:90,2]=(p[1:90,1]-u[1:90,1].^2.0-u[1:90,2].^2.0).*u[1:90,1]+p[1:90,2].*u[1:90,1] + 0.5*(-p[: , end].*u[:,2]+p[:,4:end-1] *u[:,2])

end


function ?_Shopf(du,u,p,t)

du[1:90,1]=0.04*ones(90,1)
du[1:90,2]=0.04*ones(90,1)


end


#initial condition
u0=-0.1*ones(90,2);
#initial time
t0=0.0;
#final time
tend=200.0;
#setting parameter matrix
p0=[0.1 , 2*pi*0.04]
push!(p0,-p0[2])
p=p0'.*ones(90,3);
SC=SC;
p=[p SC]
p=[p sum(SC,dims=2)]
#


#col 1 :alpha
#col 2-3 : [w0 -w0]

#col 3-93 : …