请问Julia有提供类似matlab的quiver3的函数吗?

我想解决的问题是这样的,对于这样一个欧拉方程(章动)(见下图),做出其相平面的轨迹图。
code_fig

我可以通过plotly作为后端实现这样的效果,
euler_n

但我希望能做出其运动轨迹,类似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="")

谢谢!

试试这个

1赞

多谢!!我再去看看这个issue!

京ICP备17009874号-2