Tartalomjegyzék

Rössler attractor

TELJES MATLAB SCRIPT KIEGÉSZÍTŐ FÜGGVÉNYEKKEL

file:   d2017_09_20_Rossler_attractor.m
author: Peter Polcz <ppolcz@gmail.com>
Created on 2017. September 20.
Output:
┌d2017_09_20_Rossler_attractor
│   - Persistence for `d2017_09_20_Rossler_attractor` initialized [run ID: 8148, 2017.09.20. Wednesday, 20:59:46]
│   - Script `d2017_09_20_Rossler_attractor` backuped

Rössler attractor on the Wikipedia

syms x y z real

% Chaotic
a = 0.2;
b = 0.2;
c = 5.7;

% Periodic
a = 0.1;
b = 0.1;
c = 12;

f_sym = [
    -y - z
    x + a*y
    b + z*(x - c)
    ];

f_fh = matlabFunction(f_sym, 'vars', {'t' [x;y;z]});

[t,r] = ode45(f_fh, [0,1000], [1;1;0]);

% Plot only the last part of the simulation (do not show transients)
I = round(numel(t)/2):numel(t);

plot3(r(I,1),r(I,2),r(I,3))

End of the script.

Output:
└ 2.2904 [sec]