Lorenz Attractor
3 views (last 30 days)
Show older comments
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 Comments
Answers (1)
ANTONINO BIANCUZZO
on 19 Mar 2022
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);
];
0 Comments
See Also
Categories
Find more on Graphics Object Identification in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!