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: 7733, 2017.09.20. Wednesday, 21:34:43]
│   - Script `d2017_09_20_Rossler_attractor` backuped

System model

Rössler attractor on the Wikipedia

syms x y z real
syms a b c real


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

Chaotic

a = 0.2;
b = 0.2;
c = 5.7;
f_fh = matlabFunction(subs(f_sym), 'vars', {'t' [x;y;z]});
[t,r] = ode45(f_fh, [0,1000], [1;1;0]);
plot3(r(:,1),r(:,2),r(:,3))
xlabel 'x'; ylabel 'y'; zlabel 'z'

Periodic

a = 0.1;
b = 0.1;
c = 12; % 4,6,8
f_fh = matlabFunction(subs(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))
xlabel 'x'; ylabel 'y'; zlabel 'z'

End of the script.

Output:

└ 2.9666 [sec]

Published with MATLAB® R2015a, and sanitarized by BeautifulSoup (Python script using BSoup is written by Polcz).