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

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

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

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

101

(2.3.1), то x = f (x, a) будет также решением той же системы уравнений

Fk x = 0, (k = 1, . . . , s) . (2.3.2)

При этом говорят, что система уравнений (2.3.1) инвариантна относительно группы преобразований G или что эта система уравнений допускает группу преобразований G.

Справедливы две важные теоремы, которые приведем без доказательства.

Теорема 2.3. Система уравнений (2.3.1) инвариантна относительно группы преобразований G тогда и только тогда, когда

XFk |M = 0 (k = 1, . . . , s) .

(2.3.3)

Доказательство необходимого условия этой теоремы проводится аналогично доказательству теоремы 2.2., а доказательство достаточного условия можно посмотреть в [30, 31].

Теорема 2.4. Поверхность M , заданная системой уравнений (2.3.1), инвариантная относительно некоторой группы преобразований G, может быть задана системой уравнений выраженной через независимые инварианты группы, если инфинитезимальный оператор не обращается в ноль на поверхности

M .

Суть этой теоремы выражает тот факт, что поверхность M может быть записана в виде

Φk (I1, . . . , IN −1) = 0,

(k = 1, . . . , s) ,

(2.3.4)

где I1 (x) , . . . , IN −1 (x) — инварианты группы преобразований

G.

102

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

2.4.Групповой анализ дифференциальных уравнений

Дифференциальные уравнения содержат независимые, зависимые переменные и производные. Предположим, что в En+m число независимых переменных n: x = x1, . . . , xn , а число зависимых переменных m: u = u1, . . . , um . Причем они преобразуются в соответствии с формулами

x = f (x, u, a) ,

(2.4.1)

u = g (x, u, a) .

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

Поскольку в дифференциальные уравнения входят производные от координат вектора u по компонентам вектора x, то введем дополнительные переменные

 

 

 

 

 

u = uik

(i = 1, . . . , n; k = 1, . . . , m) .

(2.4.2)

1

 

 

 

 

Эти координаты являются производными 1-го порядка

 

 

uik =

∂uk

.

 

 

∂xi

 

 

 

 

 

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

u = uijk

(i, j = 1, . . . , n; k = 1, . . . , m) ,

(2.4.3)

2

 

 

и т.д.

Когда точка в En+m преобразуется в соответствии с (2.4.1), то при этом преобразуются и производные. Предположим, что компоненты uki преобразуются по формуле:

uk = hk

x, u, u, a .

(2.4.4)

i

i

1

 

 

 

 

2.4. Групповой анализ дифференциальных уравнений

103

Введем пространство EN , где N = n + m + nm, в котором точка имеет координаты xi, uk , uki . Семейство преобразований (2.4.1) и (2.4.4) образует продолженную группу G1. Если инфинитезимальный оператор группы G имеет вид

X = ξi

+ ηk

(i = 1, . . . , n; k = 1, . . . , m) , (2.4.5)

∂xi

∂uk

 

 

 

то оператор группы G1, продолженный на пространство EN , представляется в виде

 

 

 

 

k

 

 

 

i

 

 

k

 

k

 

 

 

 

 

X = X + ξi

 

 

 

= ξ

 

 

 

+

η

 

 

+ ξi

 

 

 

.

(2.4.6)

 

 

∂uk

∂xi

 

∂uk

∂uk

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

Координаты касательного векторного поля в этой формуле

определяются, как и прежде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

=

∂f i

 

k

=

∂gk

 

 

 

k

=

∂hk

 

 

 

 

|a=0 ,

 

 

 

 

|a=0

,

 

i

 

|a=0 .

(2.4.7)

ξ

 

 

η

 

 

 

ξi

 

 

 

∂a

 

 

∂a

∂a

 

Поскольку решение дифференциального уравнения является со-

отношением между u = ϕ (x), то компоненты ξki можно выразить через ξi и ηk .

Решение дифференциального уравнения в преобразованных координатах имеет вид соотношения

u = ϕ x ,

 

которое можно записать в виде

f (x, u, a) .

 

g (x, u, a) = ϕ

(2.4.8)

После введения переменных

u,

u, . . . , дифференциальное

 

1

2

 

уравнение может рассматриваться как алгебраическое уравнение от переменных x, u, u, u и т. д., некоторые переменные в кото-

12

ром связаны дифференциальными соотношениями.

104

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

При переходе от переменных xi к новым переменным x i операторы дифференцирования связаны соотношениями

 

 

Di = Dif j Dj ,

 

 

(2.4.9)

где

 

 

 

 

 

 

 

 

Di =

+ up

+ up

+ . . . .

(2.4.10)

∂xi

 

p

 

i ∂up

ij ∂u

j

 

 

 

 

 

 

 

 

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

Дифференцируя соотношение (2.4.8) по переменной xi, имеем

Digk = ujk Dif j .

Дифференцируя последнее выражение по a и полагая a = 0, приходим к выражению

D

ηk = ξkδ

ij

+ ukD ξj

(i = 1, . . . , n; k, j = 1, . . . , m) ,

i

j

j i

(2.4.11)

 

 

 

 

где δij — символ Кронекера. Из (2.4.11) имеем

ξik = Diηk − ujk Diξj .

(2.4.12)

Формулы для второго, третьего и более высших продолжений могут быть получены аналогичным путем. Например, координаты ξkij вычисляются по формулам

ξijk = Dj ξik − uilk Dj ξl .

(2.4.13)

Как следует из предложенных выражений, для построения

продолженного оператора X требуется знать координаты ξi и ηk

1

оператора X.

2.4. Групповой анализ дифференциальных уравнений

105

Соответственно, дифференциальное уравнение p-го порядка можно рассматривать как некоторую поверхность в пространстве

переменных x, u, . . . , u:

p

F x, u, u, . . . , u = 0.

(2.4.14)

1p

При этом симметрия дифференциальных уравнений исследуется, как симметрия алгебраических уравнений от переменных x, u, u, . . . , u.

1p

Пусть X будет инфинитезимальным оператором продолжен-

p

ной группы Gp, тогда обобщение теоремы 2.3 на случай дифференциального многообразия (2.4.14) приводит к критерию инвариантности в виде

X F

x, u, u, . . . ,

u

F =0 = 0.

(2.4.15)

p

1

p

|

 

Определение 2.5. Уравнение (2.4.15), полученное в резуль-

тате действия оператора X на дифференциальное многообразие

p

(2.4.14), называется определяющим уравнением.

Вследствие независимости компонент векторов u, u, . . . , u и

1 p

их произведений уравнение расщепляется на переопределенную систему уравнений, из которой находятся координаты векторного поля ξi и ηk. Определяющее уравнение (2.4.15) всегда линейно относительно координат касательного векторного поля, и поэтому задача нахождения группы преобразований G также является линейной, что особенно важно при анализе нелинейных дифференциальных уравнений. Линейность определяющего уравнения относительно координат касательного векторного поля следует

из линейности оператора X группы преобразований Gp продол-

p

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

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

106

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

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

Для того чтобы найти определяющее уравнение, зададим оператор X с компонентами касательного векторного поля ξi и ηk :

X = ξi ∂xi + ηk ∂uk .

Затем продолжим этот оператор, учитывая соответствующее число производных в исходном дифференциальном уравнении:

X = X + ξi

+ ξi

+ . . .

(2.4.16)

 

 

p

k ∂uik

jk ∂uijk

 

 

В этом операторе компоненты продолженного векторного поля ξik , ξijk , . . . определяются по формулам (2.4.12) и (2.4.13) и по формулам продолжения более высокого порядка.

Действуя оператором (2.4.16) на исходное дифференциальное уравнение, получаем определяющее уравнение. Приравнивая выражения при разных степенях координат xi, uk , uki и т. д., приходим к системе переопределенных дифференциальных уравнений для определения компонент ξi и ηk . Решая эту систему уравнений, находим множество инфинитезимальных операторов, допускаемых исходным дифференциальным уравнением.

2.5. Группы преобразований уравнения 2-го порядка

107

2.5.Группы преобразований, допускаемые обыкновенным дифференциальным уравнением 2-го порядка

Вкачестве примера рассмотрим групповой анализ обыкновенного дифференциального уравнения 2-го порядка

yxx = f (x, y, yx) ,

(2.5.1)

где f — гладкая функция от x, y и yx.

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

В этом уравнении четыре переменные x, y,

yx и yxx, поэтому раз-

мерность продолженного пространства N = 4.

Пусть общий вид оператора группы G, которая допускается

уравнением (2.5.1), имеет вид

 

 

 

 

X = ξ

+ η

.

(2.5.2)

∂x

 

 

 

∂y

 

Обозначая координаты касательного векторного поля за ξ1 и ξ11, имеем оператор X в виде

 

2

 

 

 

 

 

 

 

 

X = ξ

+ η

+ ξ1

+ ξ11

.

(2.5.3)

∂x

∂y

∂yx

 

2

 

 

 

∂yxx

 

Действуя X на уравнение (2.5.1), получаем соотношение

2

ξ11 ξ

∂f

η

∂f

ξ1f

∂f

= 0.

(2.5.4)

 

 

 

∂x

∂y

∂yx

Пусть для определенности f = 0, тогда из (2.5.3) получаем

ξ11 = 0.

(2.5.5)

108

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

По формулам первого и второго продолжения (2.4.12) и (2.4.13) находим выражения координат ξ1 и ξ11 через координаты ξ и η:

ξ1 = ηx + yxηy − yx x + yxξy ) = ηx + yx y ξx) − yx2 ξy , (2.5.6)

ξ11 = Dx 1) − yxxDx (ξ) = ηxx + yxx y ξx) +

+yx

yx ξxx) 2yxyxxξy (yx)2 ξxy +

 

+yx

ηxy + yx yy ξxy ) (yx)2 ξyy

(2.5.7)

−yxx x + yxξy ) .

Из выражения (2.5.7) следует, что дифференциальному уравнению второго порядка

yxx = 0

(2.5.8)

соответствует определяющее уравнение в виде:

 

ηxx + yx (2ηxy ξxx) + yx2 yy xy ) (yx)3 ξyy = 0,

(2.5.9)

из которого получаем расщепленные уравнения относительно координат ξ и η в виде

ηxx = 0, ξyy = 0, yx ξxx = 0, ηyy yx = 0.

(2.5.10)

Из первого и второго уравнений (2.5.10) имеем

η= ϕ1 (y) x + ϕ2 (y) ,

ξ= Ψ1 (x) y + Ψ2 (x) .

Подставляя эти выражения в третье и четвертое уравнения, получаем равенства

1y = Ψ1xxy + Ψ2xx,

(2.5.11)

ϕ1yy x + ϕ2yy = 2Ψ1x ,

2.5. Группы преобразований уравнения 2-го порядка

109

откуда находим

Ψ1 (x) = c1x + c2,

Ψ2 (x) = c3x2 + c4x + c5,

(2.5.12)

ϕ1 (y) = c3y + c6, ϕ2 (y) = c1y2 + c7y + c8.

Координаты касательного векторного поля ξ и η с учетом (2.5.12) выражаются формулами

ξ (x, y) = (c1x + c2) y + c3x2 + c4x + c5,

(2.5.13)

η (x, y) = (c3y + c6) x + c7y + c8 + c1y2.

(2.5.14)

В выражениях ci (i = 1, . . . , 8) — произвольные постоянные. Полагая постоянные ci в (2.5.13) и в (2.5.14) последовательно

равными единице, а остальные нулю, получаем операторы, которые допускаются уравнением (2.5.8):

 

 

X1 = xy

+ y2

, X2

= y

 

 

,

 

 

 

 

∂x

∂y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x

 

 

X3

= x2

 

+ xy

 

,

X4 = x

 

,

X5 =

 

 

,

(2.5.15)

∂x

 

 

 

 

 

 

∂x

 

 

 

 

∂y

 

 

 

∂x

 

 

 

 

 

 

 

 

X6

= x

 

, X7 = y

 

, X8 =

 

.

 

 

 

 

∂y

∂y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂y

 

 

Операторы (2.5.15) соответствуют полной симметрии обыкновенного дифференциального уравнения второго порядка.

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

110

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

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

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

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

dy

=

Q (x, y)

,

(2.6.1)

dx

P (x, y)

(где P (x, y) и Q (x, y) — некоторые функции от x и y) может быть проинтегрировано с помощью интегрирующего множителя.

Уравнение (2.6.1) можно записать в виде

Q (x, y) dx − P (x, y) dy = 0,

(2.6.2)

которое равносильно уравнению в частных производных первого порядка

P

∂F

+ Q

∂F

= 0

(2.6.3)

∂x

∂y

 

 

 

 

в том смысле, что интеграл F (x, y) = C (C — постоянная) является решением уравнения (2.6.3) и, наоборот, всякое решение уравнения (2.6.3) — решение уравнения (2.6.2).

Предположим, что уравнение (2.6.2) допускает некоторую группу преобразований с инфинитезимальным оператором в виде

X = ξ

+ η

.

(2.6.4)

∂x

 

 

 

∂y

 

Под действием этого оператора всякое решение F (x, y) = C1 уравнения (2.6.3) переходит снова в решение, так что решением