5

我们有系统:

x'[t] == x[t] - 5 y[t] + z[t]
y'[t] == 3 x[t] - 3 y[t] - 3 z[t]
z'[t] == -2 x[t] + 10 y[t] + 4 z[t]

和初始条件:

x[0] == .01
y[0] == 3
z[0] == 0

我制作了具体的情节:

eqn = {x'[t] == x[t] - 5 y[t] + z[t],  y'[t] == 3 x[t] - 3 y[t] - 3 z[t], 
z'[t] == -2 x[t] + 10 y[t] + 4 z[t]}; 

sol = NDSolve[{eqn, x[0] == .01, y[0] == 3, z[0] == 0}, {x[t], y[t], 
z[t]}, {t, -5, 5}]

{xde[t_], yde[t_], zde[t_]} = {x[t], y[t], z[t]} /. Flatten[sol]

 ParametricPlot3D[{xde[t], yde[t], zde[t]}, {t, 0, 10}, 
 AxesLabel -> {"x", "y", "z"}, 
 PlotRange -> {{-15, 15}, {-15, 15}, {-15, 15}}]

我知道如何选择一个随机点来绘制整个轨迹,但我找不到一种方法来动画沿着绘制的轨迹移动的点。在此特定示例中,该点应位于 t == 0 并一直移动直到 t == 2。

4

1 回答 1

6

这在 Mathematica 中非常简单 - 使用交互式界面:

Animate[Show[ParametricPlot3D[{xde[t], yde[t], zde[t]}, {t, 0, 10}, 
   AxesLabel -> {"x", "y", "z"}, 
   PlotRange -> {{-5, 15}, {-5, 5}, {-5, 15}}],
  Graphics3D[{Red, PointSize[.05], Point[{xde[T], yde[T], zde[T]}]}]], {T, 0, 2}]

在此处输入图像描述

于 2012-12-17T22:08:07.547 回答