用Mathematica绘制洛伦兹模型的状态方程

偶然看到的，原文：你看见蝴蝶翅膀上的数学公式了吗？——哆嗒数学网

看到matlab代码代码就像用麦酱改写一下试试。。

原matlab代码：

f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+28*x(2)-x(3)];
t_final=100;x0=[0;0;1e-10];[t,x]=ode45(f,[0,t_final],x0);
plot3(x(:,1),x(:,2),x(:,3));

其实就几个微分方程和初始条件，改起来还是很容易的

sol=NDSolve[{
x[1]'[t]==-8/3x[1][t]+x[2][t]x[3][t],
x[2]'[t]==-10x[2][t]+10x[3][t],
x[3]'[t]==-x[1][t]x[2][t]+28x[2][t]-x[3][t],
x[1][0]==0,x[2][0]==0,x[3][0]==1*^-10},Array[x,3],{t,0,100}]
ParametricPlot3D[{x[1][t],x[2][t],x[3][t]}/.First@sol,{t,0,100},PlotStyle->Thin]

为了保持与原程序一致，这里使用的x[i]而不是{x,y,z}
结果：