Lorenz Attractor

Question

James Hendren 2012-3-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/32673-lorenz-attractor

回答： ANTONINO BIANCUZZO 2022-3-19

Does anyone have a script written to solve lorenz attractors and them graph them? If so could you send the script to me please?

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

请先登录，再进行评论。

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

Answer 1

ANTONINO BIANCUZZO 2022-3-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/32673-lorenz-attractor#answer_921864

在 MATLAB Online 中打开

clear all, close all, clc
Beta = [10; 28; 8/3]; % chaotic values
x0 = [0; 1; 20]; % initial condition
dt = 0.001;
tspan = dt:dt:50; % 50.000 time span
%% ode options
options = odeset('RelTol',1e-12,'AbsTol',1e-12*ones(1,3));
[t,x] = ode45(@(t,x)lorenz(t,x,Beta),tspan,x0,options);
%% lorenz plot
plot3(x(:,1),x(:,2),x(:,3));
function dx = lorenz(t,x,Beta)
%% Lorenz equations
dx = [                                         
    Beta(1)*(x(2)-x(1));
    x(1)*(Beta(2)-x(3)) - x(2);
    x(1)*x(2) - Beta(3)*x(3);
];