Tartalomjegyzék

Lorentz system

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

file:   d2017_09_20_Lorentz_system.m
author: Peter Polcz <ppolcz@gmail.com>

Created on 2017. September 20.

Output:

┌d2017_09_20_Lorentz_system
│   - Persistence for `d2017_09_20_Lorentz_system` initialized [run ID: 1754, 2017.09.20. Wednesday, 21:36:37]
│   - Script `d2017_09_20_Lorentz_system` backuped

sigma = 10;
beta = 8/3;
rho = 28;

f = @(t,x) [
    -sigma*x(1) + sigma*x(2)
    rho*x(1) - x(2) - x(1)*x(3)
    -beta*x(3) + x(1)*x(2)
    ];

[t,x] = ode45(f,[0 50],[1 1 1]);
figure, plot3(x(:,1),x(:,2),x(:,3))

figure,
subplot(311), plot(t, x(:,1));
subplot(312), plot(t, x(:,2));
subplot(313), plot(t, x(:,3));

End of the script.

Output:

└ 0.79986 [sec]

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