【问题标题】:Plotting ODEs on maple在 maple 上绘制 ODE
【发布时间】:2016-06-22 12:52:19
【问题描述】:

我是 Maple 的初学者,我不知道如何绘制这样的 ODE 系统。

这是代码

g := 9.81;
m := 0.25e-2;
h := 0.73e-3;
xi := -.5;
yi := -.5;
zi := .3;
vxi := 8;
vyi := 2;
vzi := 3;
wx := (0*2)*Pi;
wy := (13*2)*Pi;
wz := (0*2)*Pi;
eqx := m*(diff(x(t), t, t)) = -h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+   (diff(z(t), t))^2)*(diff(x(t), t))+a*(wy*(diff(z(t), t))-wz*(diff(y(t), t)));
eqy := m*(diff(y(t), t, t)) = -h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+(diff(z(t), t))^2)*(diff(y(t), t))+a*(wz*(diff(x(t), t))-wx*(diff(z(t), t)));
eqz := m*(diff(z(t), t, t)) = -m*g-h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+(diff(z(t), t))^2)*(diff(z(t), t))+a*(wx*(diff(y*t, t))-wy*(diff(x(t), t)));
ei := x(0) = xi, y(0) = yi, z(0) = zi, (D(x))(0) = vxi, (D(y))(0) = vyi, (D(z))(0) = vzi;
F := dsolve({ei, eqx, eqy, eqz}, {x(t), y(t), z(t)}, numeric, output = listprocedure);

提前致谢。

【问题讨论】:

    标签: plot ode maple


    【解决方案1】:

    命令是

    plots:-odeplot(F, [x,y,z](t), 0..100);
    

    【讨论】:

    • 警告,无法评估超过初始点的解决方案,问题可能很复杂,最初是单一的或设置不正确
    • 该错误信息与绘图命令无关。该消息是关于系统本身的数学特性的。今晚晚些时候我会处理它。
    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2019-07-17
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2020-12-17
    相关资源
    最近更新 更多