Σύστημα Λόρεντζ

μικρογραφία|300x300εσ|Ο ελκυστής του Λόρεντζ
Το **Σύστημα Λόρεντζ** είναι μια μορφοκλασματική δομή που αντιστοιχεί στη μακροχρόνια συμπεριφορά του ταλαντωτή Λόρεντζ. Ο ελκυστής δείχνει πώς οι διάφορες μεταβλητές του δυναμικού συστήματος εξελίσσονται με την πάροδο του χρόνου σε μια μη περιοδική τροχιά.

Το 1963, ο μετεωρολόγος Έντουαρντ Λόρεντζ ήταν ο πρώτος που τόνισε την πιθανή χαοτική φύση της μετεωρολογίας. Το μοντέλο Λόρεντζ, επίσης γνωστό ως δυναμικό σύστημα Λόρεντζ ή ταλαντωτής Λόρεντζ, είναι ένα απλουστευμένο πρότυπο των μετεωρολογικών φαινομένων που βασίζεται στη ρευστομηχανική. Το μοντέλο αυτό είναι ένα τρισδιάστατο δυναμικό σύστημα που παράγει χαοτική συμπεριφορά υπό ορισμένες συνθήκες.

Το πρότυπο του Λόρεντζ είχε σημαντικό αντίκτυπο στην ανάδειξη των πιθανών ορίων της ικανότητας πρόβλεψης της μακροπρόθεσμης κλιματικής και μετεωρολογικής εξέλιξης. Αποτελεί σημαντικό στοιχείο της θεωρίας ότι η ατμόσφαιρα των πλανητών και των αστέρων μπορούν να περιλαμβάνουν μια μεγάλη ποικιλία οιονεί περιοδικών καθεστώτων και υπόκεινται σε απότομες και φαινομενικά τυχαίες μεταβολές. Αποτελεί επίσης ένα χρήσιμο παράδειγμα στη θεωρία των δυναμικών συστημάτων, χρησιμεύοντας ως πηγή νέων μαθηματικών εννοιών.

## Πρότυπο Λόρεντζ

Ο ελκυστής και οι σχετικές εξισώσεις δημοσιοποιήθηκαν το 1963 από τον Έντουαρντ Λόρεντζ.

Από μαθηματική άποψη, η σύζευξη της γήινης ατμόσφαιρας με τον ωκεανό περιγράφεται από το σύστημα συζευγμένων μερικών διαφορικών εξισώσεων Ναβιέρ-Στόκες της ρευστομηχανικής. Αυτό το σύστημα εξισώσεων ήταν πολύ περίπλοκο για να επιλυθεί αριθμητικά με τους υπολογιστές που υπήρχαν την εποχή του Λόρεντζ. Έτσι, ο Λόρεντζ είχε την ιδέα να αναζητήσει ένα εξαιρετικά απλουστευμένο πρότυπο αυτών των εξισώσεων για τη μελέτη μιας συγκεκριμένης φυσικής κατάστασης: του φαινομένου της συναγωγής Ρέιλεϊ-Μπενάρ. Το αποτέλεσμα ήταν ένα δυναμικό διαφορικό σύστημα με τρεις μόνο βαθμούς ελευθερίας, το οποίο ήταν πολύ πιο απλό να ολοκληρωθεί αριθμητικά από τις αρχικές εξισώσεις.

### Δυναμικό διαφορικό σύστημα Λόρεντζ
Αυτό το διαφορικό σύστημα γράφεται ως εξής:
:\begin{cases} \dfrac{\mathrm{d}x}{\mathrm{d}t}=\sigma\left[y(t)-x(t)\right]\\
\dfrac{\mathrm{d}y}{\mathrm{d}t}=\rho \, x(t) - y(t) - x(t) \, z(t)\\
\dfrac{\mathrm{d}z}{\mathrm{d}t} =x(t) \, y(t) - \beta \, z(t) \end{cases}

όπου ,  και  είναι τρεις αυστηρά θετικές, σταθερές παράμετροι. Δεδομένου ότι το σύστημα προκύπτει από την απλοποίηση των εξισώσεων που διέπουν τη συναγωγή Rayleigh-Mπενάρ, οι παράμετροι ονομάζονται από τη φυσική τους προέλευση:  είναι ο αριθμός Prandtl, και , που καταχρηστικά αναφέρεται ως "αριθμός Rayleigh", είναι στην πραγματικότητα ο λόγος του αριθμός Rayleigh Ra προς τον κρίσιμο αριθμό Rayleigh Rac (Rac είναι η τιμή του Ra πάνω από την οποία το φυσικό σύστημα δεν μπορεί να παραμείνει αναλλοίωτο με την πάροδο του χρόνου, όπου επομένως υπόκειται σε κινήσεις συναγωγής).

Οι δυναμικές μεταβλητές ,  και  αντιπροσωπεύουν την κατάσταση του συστήματος σε κάθε χρονική στιγμή. Η φυσική ερµηνεία είναι η εξής:  είναι ανάλογη της έντασης της κίνησης συναγωγής,  είναι ανάλογη της διαφοράς θερµοκρασίας µεταξύ των ανοδικών και καθοδικών ρευµάτων και  είναι ανάλογη της απόκλισης του κατακόρυφου προφίλ θερµοκρασίας από ένα γραµµικό προφίλ..

Συχνά υποθέτουμε  = 10  και  = 8/3, με το ρ να παραμένει μεταβλητό. Το σύστημα παρουσιάζει χαοτική συμπεριφορά για  = 28, μετά από κάποιο χρονικό διάστημα στο ψευδο-περιοδικό καθεστώς (της τάξης των 30 δευτερολέπτων). Ο χρόνος αυτός εξαρτάται από το  : μπορεί να φτάσει τα 180 δευτερόλεπτα για  = 24,3. Το κατώφλι για την εγκαθίδρυση του χαοτικού καθεστώτος είναι μεταξύ 24,29 και 24,30 για αυτές τις τιμές  και .

### Σημεία ισορροπίας
Τα σημεία ισορροπίας ή σταθερά σημεία του συστήματος είναι οι σταθερές (*x*,*y*,*z*) λύσεις του διαφορικού συστήματος. Υπάρχουν ένα ή τρία τέτοια σημεία:
* το σταθερό σημείο (0, 0, 0) υπάρχει ανεξάρτητα από τις τιμές των πραγματικών παραμέτρων ,  και . Είναι σταθερό όταν \rho (σύγκλιση στην κατάσταση ηρεμίας, χωρίς συναγωγή), ασταθές διαφορετικά (συναγωγή) ,
* τα δύο συμμετρικά σταθερά σημεία \left( - \sqrt{\beta(\rho - 1)},-\sqrt{\beta( \rho - 1)},  \rho - 1\right) και \left( \sqrt{\beta( \rho - 1)},\sqrt{\beta( \rho - 1)}, \rho - 1 \right) υπάρχουν μόνο για την  \rho\ge1. Ανάλογα με τις τιμές των  και , είναι σταθερές (σταθερή συναγωγή) ή ασταθείς (περιοδικά κυμαινόμενη συναγωγή ή χαοτική).

## Περιγραφή

Ο ελκυστής Λόρεντζ ορίζεται ως το σύνολο των μακροχρόνιων τροχιών του παραπάνω δυναμικού συστήματος Λόρεντζ.

Πράγματι, όταν οι παράμετροι ,  και  παίρνουν τις ακόλουθες τιμές:  = 10,  = 28 και  = 8/3, το διαφορικό δυναμικό σύστημα Λόρεντζ παρουσιάζει έναν παράξενο ελκυστή με τη μορφή φτερών πεταλούδας, που απεικονίζεται στο διπλανό σχήμα.

Για σχεδόν όλες τις αρχικές συνθήκες (εκτός από εκείνες για τα σταθερά σημεία), η τροχιά του συστήματος πλησιάζει γρήγορα τον ελκυστή, με την περιφορά της τροχιάς αρχικά να τυλίγεται γύρω από τη μία πτέρυγα, στη συνέχεια να πηδά από τη μία πτέρυγα στην άλλη για να αρχίσει να τυλίγεται γύρω από την άλλη πτέρυγα, και ούτω καθεξής, με έναν προφανώς ακανόνιστο τρόπο.

Η ύπαρξη ενός παράξενου ελκυστή για ορισμένες τιμές παραμέτρων είχε πιθανολογηθεί από τον Έντουαρντ Λόρεντζ ήδη από το 1963 βάσει αριθμητικών προσομοιώσεων. Ωστόσο, έπρεπε να περιμένουμε μέχρι το 2001 για μια αυστηρή απόδειξη του γεγονότος αυτού από τον Ουόρικ Τάκερ.

Σε αυτές τις περιπτώσεις, ο ελκυστής είναι ένα φράκταλ με διάσταση Χάουστορφ μεταξύ 2 και 3.


Image:Lorenz attractor yb.svg| = 28,  = 10,  = 8/3
Image:LorenzAttractor.png
Image:Lorenz Ro28-200px.png| = 28
Image:Lorenz Ro15-200px.png| = 15
File:Lorenz apparition small.gif|Κινούμενη εικόνα του  = 0 ως 28


## Παραλλαγές
Προκειμένου να μελετηθεί μαθηματικά ο ελκυστής Λόρεντζ, εμφανίστηκαν στην επιστημονική βιβλιογραφία αυστηρά αλλά διακριτά πρότυπα . Το κύριο πλεονέκτημά τους ήταν ότι ήταν ευκολότερο να μελετηθούν και, κυρίως, ότι μπορούσαμε να είμαστε σίγουροι για την ύπαρξή τους (αν και δεν γνωρίζαμε αν ο ελκυστής Λόρεντζ περιγράφεται σωστά από αυτά τα πρότυπα, βλ. παραπάνω).

## Ανάλυση

### Προσομοίωση στο Julia

```julia
using Plots
# define the Lorenz attractor
@kwdef mutable struct Lorenz
    dt::Float64 = 0.02
    σ::Float64 = 10
    ρ::Float64 = 28
    β::Float64 = 8/3
    x::Float64 = 2
    y::Float64 = 1
    z::Float64 = 1
end

function step!(l::Lorenz)
    dx = l.σ * (l.y - l.x);         l.x += l.dt * dx
    dy = l.x * (l.ρ - l.z) - l.y;   l.y += l.dt * dy
    dz = l.x * l.y - l.β * l.z;     l.z += l.dt * dz
end

attractor = Lorenz()

# initialize a 3D plot with 1 empty series
plt = plot3d(
    1,
    xlim = (-30, 30),
    ylim = (-30, 30),
    zlim = (0, 60),
    title = "Lorenz Attractor",
    marker = 2,
)

# build an animated gif by pushing new points to the plot, saving every 10th frame
@gif for i=1:1500
    step!(attractor)
    push!(plt, attractor.x, attractor.y, attractor.z)
end every 10
```


### Προσομοίωση με το MATLAB

```matlab
% Resolver para el intervalo de tiempo [0,100] con condiciones iniciales [1,1,1]
% ''f'' es un conjunto de ecuaciones diferenciales
% ''a'' es un arreglo que contiene variables x, y, z
% ''t'' es la variable de tiempo

sigma = 10;
beta = 8/3;
rho = 28;
f = @(t,a) [-sigma*a(1) + sigma*a(2); rho*a(1) - a(2) - a(1)*a(3); -beta*a(3) + a(1)*a(2)];
[t,a] = ode45(f,[0 100],[1 1 1]);     % Solución de EDO de Runge-Kutta de 4.º/5.º orden
plot3(a(:,1),a(:,2),a(:,3))
```


### Προσομοίωση με το Python

```numpy
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
from mpl_toolkits.mplot3d import Axes3D

rho = 28.0
sigma = 10.0
beta = 8.0 / 3.0

def f(state, t):
    x, y, z = state  # Desempaqueta el vector de estado
    return sigma * (y - x), x * (rho - z) - y, x * y - beta * z  # Derivadas

state0 = [1.0, 1.0, 1.0]
t = np.arange(0.0, 40.0, 0.01)

states = odeint(f, state0, t)

fig = plt.figure()
ax = fig.add_subplot(projection='3d')
ax.plot(states[:, 0], states[:, 1], states[:, 2])
plt.show()
```


## Βιβλιογραφία
* 
* 

## Δημοσιεύσεις
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* Shen, B.-W. (2015-12-21). "Nonlinear feedback in a six-dimensional Lorenz model: impact of an additional heating term". *Nonlinear Processes in Geophysics*. **22** (6): 749–764. doi:10.5194/npg-22-749-2015. ISSN 1607-7946.
* 
* 
* 
* 
* 

## Δείτε επίσης
*Αλγόριθμος διαμαντιού τετραγώνου
*Καμπύλη που γεμίζει το χώρο
*Νιφάδα του Κοχ 
*Καμπύλη Χίλμπερτ
*Καμπύλη του δράκου
*Σπόγγος του Μένγκερ

*Σύνολο Μάντελμπροτ
*Σύνολο Julia
*Σύνολο Κάντορ
*Οικοδομήσιμο Σύμπαν

## Παραπομπές

 
 
Κατηγορία:Φράκταλ
Κατηγορία:Χάος
Κατηγορία:Δυναμικά συστήματα