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

Кудряшов Методы нелинейной математич 2008

.pdf
Скачиваний:
511
Добавлен:
16.08.2013
Размер:
1.5 Mб
Скачать

1.22. Система Лоренца

91

возможны как стационарный, так и хаотический режимы движения.

Когда r → r3 24, 74, седловые предельные циклы стягиваются к стационарным точкам O1 и O2, и они становятся неустойчивыми.

Винтервале r3 < r < r4 = 28 все стационарные точки O, O1 и O2 являются неустойчивыми. Устойчивым притягивающим предельным множеством (аттрактором) остается лишь аттрактор Лоренца. В динамической системе, описываемой системой уравнений Лоренца, вне зависимости от выбранных начальных условий реализуется хаотический режим движения.

Хаотический режим движения, описываемый системой уравнений Лоренца, продолжает существовать при значениях r из интервала: r4 < r < r5 148, 4. При r5 < r в фазовом пространстве вместо странного аттрактора образуется предельный цикл, и движение в системе становится периодическим.

Вработе Д. Рюэля и Ф. Такенса аттрактор Лоренца был назван странным аттрактором, поскольку он является новым ти-

пом стационарного режима движения, отличным от стационарной точки и предельного цикла. Этот аттрактор не является двумерной поверхностью. Поверхность, которая проходит через аттрактор, пересекается фазовыми траекториями в точках, образующих канторово множество [53, 55].

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

Глава 2

Элементы группового анализа дифференциальных уравнений

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

2.1.Однопараметрическая группа преобразований Ли

Взаимно однозначное отображение множества M самого на себя называется преобразованием множества M . Совокупность

2.1. Однопараметрическая группа преобразований Ли

93

всех преобразований множества M образует группу преобразований. Далее в качестве множества M возьмем N -мерное ев-

клидово пространство RN , в котором точку x RN c коор-

динатами

x1, x2, . . . ,

 

xN

 

можно преобразовать в точку x

 

RN

с

координатами x

1,

x

2, . . . , x N

x

 

 

 

 

 

 

 

 

 

 

 

. Пусть отображение

 

в

 

 

задается

 

 

 

 

 

 

 

i

=

x

 

x

 

 

 

 

 

 

 

 

формулой

 

 

 

= f (x, a) или покомпонентно x

 

f i

 

 

 

 

 

 

 

 

 

 

 

 

 

x1, x2, . . . , xN , a , где — a вещественное число: a R (R

множество всех вещественных чисел).

Относительно функций f i (x, a) предположим, что они вместе со всеми своими производными по всем компонентам вектора x и по параметру a принадлежат классу непрерывных функ-

ций. Как правило, отображения f (x, a), которые используются в групповом анализе дифференциальных уравнений, определены не на всем пространстве RN и не при всех вещественных числах R. Поэтому обычно говорят о локальных преобразованиях про-

странства V RN и R. Дадим определение.

Определение 2.1. Локальной однопараметрической группой преобразований Ли называется семейство преобразований

 

 

 

 

 

 

 

 

f (x, a), имеющее свойства:

 

 

 

 

 

1) существует a0

 

 

 

 

 

 

, такое, что f (x, a0) = x;

 

 

2) если x = f (x, a),

то

найдется a˜

 

, такое,

что

x = f (x , a˜);

 

 

 

 

 

= f (x, a) и

3) для любых a

,

b

, таких, что x

x = f (x , b) найдется

c

 

,

такое, что

x =

f (x, c),

где

c = ϕ (a, b)

4)f (x, a) C .

Каждое отображение f (x, a), которое имеет свойства 1) — 4) задает однопараметрическую непрерывную группу преобразований, которую далее обозначаем G.

Приведем примеры однопараметрических групп преобразований.

1. Группа преобразований переноса (сдвига) по оси x:

x = x + a.

(2.1.1)

94

Глава 2. Элементы группового анализа

 

Поскольку для данного преобразования a0 = 0, a˜

= −a и

x

= x + c, где c = a + b, то вышеперечисленные свойства для

однопараметрической группы преобразований выполняются.

 

2. Группа преобразований растяжений на плоскости R2 (x, y):

 

x = xea, y = ye2a.

(2.1.2)

Очевидно, что свойства 1) — 4) также выполняются.

3. Группа преобразований вращения в плоскости R2 (x, y):

x = x cos α + y sin α,

y = y cos α

 

(2.1.3)

x sin α

 

 

является также по определению однопараметрической группой преобразований Ли.

Пусть преобразования x = f (x, a) образуют группу. Можно

выбрать a0 = 0. Разложим f (x, a) в ряд Тейлора по параметру a в окрестности a = 0. Тогда поскольку f (x, 0) = x, то

 

 

 

 

 

 

x = x + a

∂f

 

|a=0 + o (a) .

(2.1.4)

∂a

Обозначим

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f

 

 

 

ξ (x) =

∂a

|a=0

(2.1.5)

за вектор касательного векторного поля группы преобразований G имеющего компоненты:

= ( 1, 2, . . . , N ) .

ξ ξ ξ ξ

Таким образом, преобразование x = f (x, a) можно представить в виде

x = x + a ξ (x) .

(2.1.6)

Справедлива следующая теорема Софуса Ли [30, 31, 32].

Теорема 2.1. Пусть функции f i (x, a) , (i = 1, . . . , N )

удовлетворяют групповым свойствам и имеют разложения

2.1. Однопараметрическая группа преобразований Ли

95

(2.1.6). Тогда они являются решениями обыкновенных дифференциальных уравнений первого порядка (называемыми уравнениями Ли) с начальными условиями

∂f i

f ,

f i |a=0 = x.

 

da = ξi

(2.1.7)

Справедливо и обратное утверждение, что при любом глад-

 

 

 

решение задачи Коши (2.1.7) образует

ком векторном поле ξ

f

группу

преобразований.

 

 

 

 

Доказательство.

Пусть f (x, a) удовлетворяет групповым свойствам, тогда, придавая параметру a приращение b, имеем

f (x, a + a) = f f (x, a) a .

 

Откуда получаем два соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f

 

+ o a) ,

 

f (x, a +

a) = f (x, a) +

 

a

∂ a

(2.1.8)

f f (x, a) a

= f (x, a) +

a

∂f

|

a=0 + o a) ,

(2.1.9)

∂ a

приравнивая которые, приходим к уравнению

 

 

 

 

 

∂f

f (x, a) , a

 

f (x, a) ,

 

 

 

= ξ

(2.1.10)

 

∂ a

совпадающему с уравнением Ли.

Пусть f (x, a) — решение задачи Коши (2.1.7). Рассмотрим две функции

u(b) = f i x , b = f i f i (x, a) , b ,

v (b) = f i (x, a + b) ,

тогда в силу уравнений Ли получаем две задачи:

dudb = ξi (u) , u |b=0 = f i (x, a) ;

96

Глава 2. Элементы группового анализа

 

 

dv

= ξi (v) ,

v |b=0 = f i (x, a) .

 

 

 

 

 

db

Поскольку функции u (b) и v (b) удовлетворяют одному и тому же обыкновенному дифференциальному уравнению и одному начальному условию, то в силу единственности решения задачи Коши приходим к равенству u (b) = v (b), что равносильно групповому свойству преобразований xi = f i (x, a).

Для группы преобразований растяжений на плоскости (2.1.2)

имеем

 

 

 

 

 

 

 

 

 

 

 

∂x

 

 

 

∂y

 

 

 

ξ1 =

 

|a=0 = x,

ξ2 =

 

|a=0 = 2y,

(2.1.11)

∂a

∂a

поэтому задача Коши для системы уравнений Ли имеет вид

 

∂x

 

 

 

∂y

 

 

 

= x ,

x |a=0 = 0,

 

 

= 2y , y |a=0 = 0.

(2.1.12)

 

∂a

∂a

Решая задачу (2.1.12), можно найти однопараметрическую группу преобразований (2.1.2). Таким образом, зная группу преобразований в силу определения (2.1.5), можно найти компоненты касательного векторного поля и, наоборот, задаваясь компонентами касательного векторного поля и затем решая задачу Коши для уравнений Ли, можно восстановить группу преобразований.

2.2.Инварианты. Инфинитезимальный оператор группы преобразований

Возьмем функцию F (x), определенную в некоторой области

RN .

Определение 2.2. Функция F (x) является инвариантом

группы преобразований , если для всех допустимых

x = f (x, a)

2.2. Инварианты. Оператор группы преобразований

97

x и a выполняется равенство

 

F x = F (x) .

(2.2.1)

Справедлива теорема, которая является критерием инвари-

антности для функции F (x).

Теорема 2.2. Функция F (x) является инвариантом отно-

сительно группы преобразований тогда и только

x = f (x, a)

тогда, когда она удовлетворяет уравнению

ξi (x)

∂F (x)

= 0.

(2.2.2)

∂xi

 

 

 

Доказательство. Пусть F (x) удовлетворяет условию (2.2.1),

тогда

f (x, a) = F x + ξ a + o (a) =

 

F

(2.2.3)

 

= F (x) + a ξi (x)

∂F (x)

+ o (a) .

 

 

∂xi

 

 

 

 

 

Откуда следует, что при условии (2.2.1) выполняется (2.2.2). Пусть теперь F (x) — решение уравнения (2.2.2), тогда, по-

скольку оно справедливо для любой точки, запишем его для

x = f (x, a):

 

ξi x

 

∂F (x )

 

 

 

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

∂x i

 

 

 

Это условие можно представить в виде

 

 

 

dF (x )

=

∂F (x ) df i (x, a)

= ξi x

∂F (x )

= 0

 

da

 

∂x i

 

 

da

∂x i

и, следовательно, F (x ) как функция от a является решением задачи Коши:

dF

= 0, F |a=0 = F (x) ,

da

откуда получаем условие (2.2.1).

98Глава 2. Элементы группового анализа

Вкритерии инвариантности (2.2.2) появился оператор

N

X = ξi ∂xi , (2.2.4)

i=1

который действует на функцию F (x).

Определение 2.3. Оператор (2.2.4) называется инфинитезимальным оператором группы G.

Координаты касательного векторного поля ξi называются также координатами инфинитезимального оператора. В определении оператора X слово “инфинитезимальный” часто опускается и тогда говорят об операторе группы G.

Критерий инвариантности (2.2.2) представляет собой линейное однородное уравнение первого порядка в частных производных, решением которого является функция от N − 1 функционально независимых инвариантов. Любой другой инвариант является функцией от независимых инвариантов.

В качестве независимых инвариантов можно взять первые интегралы

I1 (x) = c1, . . . , IN −1 (x) = cN −1

системы уравнений

 

 

 

 

 

 

 

dxi

dx2

 

dxN

 

 

 

=

 

=

. . . =

 

,

(2.2.5)

 

ξ1 (x)

ξ2 (x)

ξN (x)

полученной как решения уравнения в частных производных

1 ∂F

+ ξ

∂F 2

+ . . . + ξ

N

∂F

= 0.

(2.2.6)

ξ

 

 

 

 

∂x1

∂x2

 

∂xN

 

 

 

 

 

 

Выше мы говорили о том, что задание компонент касательного векторного поля эквивалентно заданной однопараметрической группе преобразований. При известном касательном векторном поле можно построить соответствующий инфинитезимальный оператор и определить из решения системы уравнений (2.2.5)

2.2. Инварианты. Оператор группы преобразований

99

базисные инварианты и функции, которые являются инвариантными относительно соответствующей группы преобразований.

В качестве примера рассмотрим группу преобразований рас-

тяжения в RN :

λia ,

 

x i = xi exp

(2.2.7)

где (λi — постоянные, i = 1, . . . , N ), для которой компоненты касательного векторного поля выражаются формулами

ξi = λixi.

(2.2.8)

При этом оператор группы преобразований растяжения имеет

вид

N

X = λixi ∂xi λixi ∂xi , (2.2.9)

i=1

Инварианты группы преобразований растяжения находятся как решения системы уравнений

dx1

dx2

dxN

 

 

 

 

=

 

 

 

= . . . =

 

 

 

(2.2.10)

1

x

1

2

x

2

N

x

N

λ

 

 

λ

 

 

λ

 

 

и выражаются формулами:

I1 =

x1

I2 =

x1

. . . , IN −1 =

x1

 

,

 

,

 

,

(x2)a2

(x3)a3

(xN )aN

(2.2.11)

ak = λ1/λk .

Функции, которые зависят от инвариантов (2.2.11), являются инвариантами в RN относительно группы преобразований растяжения (2.2.7).

В качестве второго примера возьмем группу преобразований вращения на плоскости R2 (x, y). Компоненты касательного векторного поля для этого случая имеют вид

ξ1 =

∂x

ξ2 =

∂y

 

 

|α=0 = −y,

 

|α=0 = x.

(2.2.12)

α

α

100

Глава 2. Элементы группового анализа

 

Оператор

группы преобразований

вращения на

плоскости

R2 (x, y) выражается формулой

 

 

 

 

 

 

 

X = −y

 

+ x

(2.2.13)

 

 

 

 

 

.

 

∂x

∂y

Решение уравнения

 

dy

 

 

 

 

 

dx

 

 

 

 

 

 

=

 

 

 

 

 

(2.2.14)

 

y

 

x

 

 

дает инвариант

 

 

 

 

 

 

 

I = x2 + y2,

 

 

(2.2.15)

который показывает, что любая функция от x2 + y2 будет являться инвариантом относительно группы преобразований вращения на плоскости.

2.3.Инвариантные уравнения

Пусть задана (N − s)-мерная поверхность M Rn следующей системой уравнений

 

 

F1 (x) = 0, . . . , Fs (x) = 0, s N.

(2.3.1)

Предположим, что ранг матрицы

 

 

"

∂Fk

"

 

 

s

"

∂xi

"

(k = 1, . . . , s; i = 1, . . . , N )

 

 

"

"M

 

равен .

"

 

"

 

 

Определение 2.4. Поверхность M, заданная системой уравнений (2.3.1), называется инвариантной относительно группы

преобразований , если в результате отображения

x = f (x, a)

каждая точка остается на этой поверхности.

Если поверхность инвариантна относительно соответствующей группы преобразований и x — решение системы уравнений