Rownanie Lorenza Przeradzki p3


Równanie Lorenza
ńł
x = (y - x)
ł
y = rx - y - xz (1)
ół
z = xy - bz,
gdzie stale dodatnie , r i b maja interpretacje fizyczna (stale Prandtla i
8
Rayleigha). Lorenz przyja l  = 10, b = .
3
V (x, y, z) := rx2 + y2 + (z - 2r)2
 jej pochodna wzdluż trajektorii:
Ł
V = -2(rx2 + ry2 + b(z - r)2 - br2).
Ł
V < 0 na zewnatrz elipsoidy E := {(x, y, z) : rx2 + ry2 + b(z - r)2 d" br2}.
Wprowadzmy oznaczenia m := supE V i dla  > 0
E := {(x, y, z) : V (x, y, z) d" m + }.
Ł
Wtedy E " int E i sup \E V := - < 0, wiec wartość funkcji V wzdluż
3
trajektorii startujacej na zewnatrz E maleje, aż do chwili, gdy trajektoria
Ł
ta ,,wejdzie do tej elipsoidy. Ponieważ V d" 0 na E \ E, wiec trajektoria
ta nigdy już nie opuści zbioru E. Zatem O+(P) jest zbiorem zwartym dla
3
każdego P " i tym samym zbiór (P ) jest niepusty, zwarty i spójny.
Punkty stale ukladu: dla r d" 1 jest tylko jeden punkt staly p0 = (0, 0, 0),
a dla r > 1 pojawiaja sie jeszcze dwa
pą = ą b(r - 1), ą b(r - 1), r - 1 .
Z analizy stabilności punktów stalych:
dla r > rH, gdzie
 + b + 3
rH := ,
 - b - 1
wszystkie punkty stale sa niestabilne. Dla r = rH mamy jedna rzeczy-
wista wartość wlasna przybliżenia liniowego ukladu w punktach pą i dwie
sprzeżone, czysto urojone: ma miejsce bifurkacja Hopfa.
1
Dywergencja pola wektorowego danego przez prawa strone ukladu Lo-
renza:
f(x, y, z) := [(y - x), rx - y - xz, xy - bz], div f = -( + b + 1) < 0.
Na podstawie tw. o dywergencji, jeśli wezmiemy dowolny zbiór mierzalny
3
D " , przez D(t) oznaczymy zbiór
{(x(t), y(t), z(t)) : (x(0), y(0), z(0)) " D},
a przez O(t) jego objetość, to O (t) = -( + b + 1)O(t), czyli
O(t) = O(0) exp(-( + b + 1)t) - 0
przy t +". Zatem mamy zbiór
A := lim D(t) = D(s),
t+"
te"0 se"t
gdzie D(0) = D = E z wcześniejszych rozważań. Zbiór A jest domkniety,
niezmienniczy i ma objetość 0. Jest on globalnym atraktorem dla ukladu
Lorenza.
K. Mishaikov, M. Mrozek: Chaos in the Lorenz equations, a computer
assisted proof, Bull. Amer. Math. Soc. 32 (1995), 66 72.
K. Mishaikov, M. Mrozek, A. Szymczak: Chaos in the Lorenz equations,
a computer assisted proof, part II  details, Mathematics of Computations
67 (1998), 1023 1046.
Uklad dynamiczny Lorenza zredukowany do atraktora A ma dwie nastepujace
cechy:
1) topologiczna tranzytywność, tzn. dla dowolnych zbiorów otwartych
U, V " A istnieje t > 0 takie, że
Ą(t, U) )" V = ";

2) wrażliwość na warunki poczatkowe, tzn. istnieje dodatnia liczba M
taka, że dla dowolnego p " A istnieja ciagi pn p i (tn) " (0, ") o wlasności
||Ą(tn, pn) - Ą(tn, p)|| e" M.
2
Warunek 1) gwarantuje niepodzielność zbioru A na mniejsze zbiory do-
mkniete i niezmiennicze. W szczególności jest on spelniony, gdy istnieje tra-
jektoria, która jest gestym podzbiorem A. Bylby też spelniony, gdyby caly
zbiór A byl jedna trajektoria okresowa.
Warunek 2) oznacza, że żadne rozwiazanie o trajektorii zawartej w A nie
jest stabilne nawet po redukcji ukladu do zbioru A. Gdyby A byl trajektoria
okresowa, to ten warunek nie bylby spelniony. Niektórzy autorzy do tych
dwóch warunków dokladaja jeszcze
3) gestość zbioru trajektorii okresowych w zbiorze A.
uklad Rsslera
ńł
x =
ł -(y + z)
y = x + ay, a = b = 0.2, c > 6, (2)
ół
z = b + xz - cz
uklad Chua
ńł
d-c
x = a(y - x - g(x)) g(x) = cx + (|x + 1| - |x - 1|)
ł
2
y = x - y + z, a = 15, b = 25.58, (3)
ół
z = -by c = -5, d = -8.
7 7
3


Wyszukiwarka

Podobne podstrony:
Zastosowania do RRCz Przeradzki p3
2002 p3 answers
uklady rownan (1)
Zestaw 1 Funkcja kwadratowa Funkcja homograficzna Równanie liniowe
przerazajacy staruch
modele rownan
Rownanie ruchu pojazdu samochodowego
Równania kwadratowe matematyka
bilans wodny metoda najmniejszych kwadratow rownanie bubendeya
rownania dynamiczne cw

więcej podobnych podstron