LorenzSystem, ode23

Question

Jovos 2016-4-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/278187-lorenzsystem-ode23

编辑： Star Strider 2016-4-10

Hi, I get a question to model the LorenzSystem. It has

dx/dt = 10(y-x)  
dy/dt = -xz+28x-y
dz/dt = xy-8z/3

Setting the initial value x = 1; y = 2; z =3. It should get a graph like this. But mine is so weird.

function chaotic
 clc;clear;
 %LorenzSystem
 y0=[1;2;3];
 soln = ode23(@f,[0 100],y0)
 t = linspace(0,100,300);
 y(:,1)=deval(soln,t,1); 
 y(:,2)=deval(soln,t,2); 
 y(:,3)=deval(soln,t,3); 
   figure
   plot3(y(:,1),y(:,2),y(:,3));
   hold on;grid on;
   end
   function dxdt = f(t,x)
   dxdt=[0;0;0];
   dxdt(1) = 10*(x(2)-x(1));
   dxdt(2) = -x(1)*x(3)+28 * x(1)- x(2);
   dxdt(3) = x(1)*x(2)-8*x(3)/3;
   end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2016-4-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/278187-lorenzsystem-ode23#answer_217228

编辑：Star Strider 2016-4-10

在 MATLAB Online 中打开

It’s always best not to ‘overthink’ problems (and I admit that I’ve don that so often, I’ve lost count).

This works:

% Lorenz System
f = @(t,x) [10*(x(2)-x(1));  -x(1).*x(3)+28 .* x(1)- x(2);  x(1).*x(2)-8*x(3)/3];
y0=[1;2;3];
[t,y] = ode23(f,[0 100],y0);
figure(1)
plot3(y(:,1),y(:,2),y(:,3));
grid on