Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Динамический хаос (ИПИС, ФКС)

.pdf
Скачиваний:
64
Добавлен:
13.02.2015
Размер:
5.81 Mб
Скачать

То же самое можно сказать и про аттрактор, соответствующий системе Лоренца: он либо должен переходить сам в себя при замене x x , y y , либо должен состоять из нескольких симметричных частей, переходящих друг в друга при замене x x , y y . Из рис. 63 можно видеть, что аттрактор Лоренца обладает указанной симметрией.

2. Ограниченность области, где может располагаться аттрактор

В фазовом пространстве системы Лоренца можно указать такую ограниченную замкнутую область, в которую фазовые траектории могут только входить, а покидать её они не могут. Именно внутри указанной ограниченной замкнутой области и должен располагаться аттрактор системы Лоренца.

Для доказательства этого рассмотрим систему уравнений Лоренца (110) и введём в них новую переменную w r z , тогда

w z bz x y b w br x y .

В результате система Лоренца примет вид

w

x

( y x)

x

 

y

y x w

y .

(129)

b w br x y

w

 

Умножая первое уравнение (129) на x , второе – на y , а третье – на w и складывая, получим

x x y y z z xy x2 y2 xy w b w2 br w xy w ,

1

d

(x2

 

1

2

 

3

 

1

2

 

1

br2

. (130)

2 dt

y2 z2 ) x

2

y

4

y2 b w

2

r

4

 

 

 

 

 

 

 

 

 

Обозначим правую часть (130) за функцию E :

 

1

2

 

3

 

1

2

 

1

br2

. (131)

E(x, y, w) x

2

y

4

y2 b w

2

r

4

 

 

 

 

 

 

 

 

Функции E(x, y, w) можно дать наглядную геометрическую интерпретацию. Уравнение E 0 , т.е.

 

1

2

 

3

 

1

2

 

1

br2

0

x

2

y

4

y2 b w

2

r

4

 

 

 

 

 

 

 

 

задаёт в пространстве переменных (x, y, w) эллипсоид с центром, смещённым относительно начала координат

(см. рис. 65).

Рис. 65. Поверхность E(x, y, w) 0 при значении параметров: b 83 , r 28 .

Поверхность этого эллипсоида делит

всё

пространство

(x, y, w) на две области: внутри эллипсоида

E 0 , вне

эллипсоида E 0 .

 

 

Введём ещё одну функцию S :

 

 

S(x, y, w) x2 y2

z2 ,

(132)

фигурирующую в левой части уравнения (130). Тогда

уравнение S Const , т.е. уравнение

 

x2 y2 z2 Const

(133)

также будет задавать в пространстве (x, y, w) некоторый эллипсоид. Выберем значение константы таким большим, чтобы эллипсоид S (133) полностью содержал внутри себя предыдущий эллипсоид E 0 (см. рис. 66).

Рис. 66. Взаимное расположение

эллипсоидов

E(x, y, w) 0 и S(x, y, w) 1000

при значении

параметров: 10 , b 83 , r 28 .

Тогда в любой точке поверхности эллипсоида S функция E 0 , а значит,

 

 

d S

 

 

 

 

 

dt 0

.

 

 

 

 

 

 

Другими

словами функция

S x2 y2 z2

будет

убывать с течением времени в любой точке поверхности

эллипсоида

S

. Это и означает, что все траектории,

пересекающие

поверхность

S , ведут только

внутрь

ограниченной ею области. Следовательно, каким бы ни был аттрактор системы уравнений Лоренца, он обязан располагаться внутри эллипсоида S .

3. Диссипативность системы Лоренца

Систему уравнений Лоренца (110) можно записать в другой форме

 

 

x F(x) ,

(134)

где x {x, y, z}, а

F(x) { ( y x), r x y xz, bz x y}

– векторное поле, составленное из правых частей уравнений Лоренца (110), которое имеет смысл поля

скоростей в

фазовом пространстве, т.е. вектор

F определяет

мгновенную

скорость

движения

изображающей

точки с координатами x в

фазовом

пространстве.