我想解决的问题是这样的,对于这样一个欧拉方程(章动)(见下图),做出其相平面的轨迹图。
我可以通过plotly作为后端实现这样的效果,
但我希望能做出其运动轨迹,类似matlab的quiver3实现的功能。
我的代码如下:
using DifferentialEquations, Plots
plotly();
function sol_test(du,u,p,t)
x,y,z=u;
du[1] = dx = y*z;
du[2] = dy = -2*x*z;
du[3] = dz = x*y;
end
plotly()
plt = plot()
for i in 1:50
u0 = .-(rand(3)*2,1.0);
u0 = u0./sqrt(sum(u0.^2))
tspan = (0.0,10.0)
prob = ODEProblem(sol_test,u0,tspan)
# sol = solve(prob,Euler();dt=1e-2);
sol = solve(prob,Tsit5());
plot!(plt, sol,vars=(1,2,3), lw=1,xlims=(-1, 1), ylims=(-1,1), zlims=(-1,1))
end
plot(plt,label="")
谢谢!