Кудряшов Методы нелинейной математич 2008
.pdf2.6. Интегрирование дифференциальных уравнений |
111 |
(2.6.3) будет также интеграл XF = C2. Однако поскольку уравнение (2.6.2) является уравнением 1-го порядка, то эти решения функционально зависимы, и поэтому
ξ |
∂F |
+ η |
∂F |
= f (F ) . |
(2.6.5) |
|
∂x |
∂y |
|||||
|
|
|
|
Учитывая (2.6.3) и (2.6.5) имеем
∂F |
= |
Qf |
|
, |
∂F |
= − |
P f |
|
. |
||
|
|
|
|
|
|
|
|
||||
∂x |
ξQ |
− |
ηP |
∂y |
ξQ |
− |
ηP |
||||
|
|
|
|
|
|
|
|
|
|
Откуда приходим к выражению
Qdx − P dy |
= |
dF |
, |
(2.6.6) |
|
ξQ − ηP |
f (F ) |
||||
|
|
|
являющемуся полным дифференциалом.
Соотношение (2.6.6) приводит к результату Софуса Ли: оператор (2.6.4) допускается уравнением (2.6.2) тогда и только тогда, когда функция
μ = |
1 |
(2.6.7) |
ξQ − ηP |
является его интегрирующим множителем. Поэтому, зная какойлибо оператор, допускаемый уравнением 1-го порядка по формуле (2.6.7), можно найти интегрирующий множитель.
Используя этот подход, проинтегрируем уравнение Рикатти, имеющее вид [32]:
dy |
+ y2 = |
2 |
. |
(2.6.8) |
dx |
|
|||
|
x2 |
|
Подставляя x = ax, y = by, находим, что уравнение (2.6.8) будет инвариантным, если b = 1/a, и поэтому оно допускает группу преобразований растяжения
x = xeα, y = ye−α. |
(2.6.9) |
112 |
Глава 2. Элементы группового анализа |
|
||||
Инфинитезимальный оператор, соответствующий преобразо- |
||||||
ваниям (2.6.9) имеет вид |
|
|
|
|
||
|
X = x |
∂ |
− y |
∂ |
(2.6.10) |
|
|
|
|
. |
|||
|
∂x |
∂y |
Уравнение (2.6.8) можно представить в виде
2 |
|
dx = 0. |
|
|||||
dy + y2 − |
|
|
(2.6.11) |
|||||
x2 |
||||||||
Интегрирующий множитель для (2.6.11) находим по формуле |
||||||||
(2.6.7) |
|
|
x |
|
|
|
|
|
μ = |
|
|
|
. |
(2.6.12) |
|||
|
|
|
||||||
x2y2 − xy − 2 |
||||||||
Он приводит к решению уравнения (2.6.11) в виде: |
|
|||||||
|
xy − 2 |
= |
c |
. |
(2.6.13) |
|||
|
|
|
||||||
|
xy + 1 |
x3 |
|
Метод интегрирующего множителя применяется только для интегрирования уравнений первого порядка. Однако группы преобразований могут быть использованы и для интегрирования уравнений более высокого порядка, поскольку они часто подсказывают упрощающую замену переменных. Применение этого под-
хода основано на следующей теореме [30, 32].
Теорема 2.5. Всякая однопараметрическая локальная группа преобразований x i = f i (x, a) в RN с помощью некоторой невырожденной замены переменных zi = zi (x) может быть
приведена к группе переноса. |
|
|
|
|
|
|
Доказательство. |
Пусть |
|
группе |
преобразований |
||
x = f i (x, a) (i = 1, . . . , |
N ) соответствует оператор |
|||||
|
X = ξi (x) |
|
∂ |
. |
(2.6.14) |
|
|
∂xi |
|||||
|
|
|
|
|
||
Этот оператор может быть представлен в виде |
||||||
X˜ = X zi |
|
∂ |
|
, |
(2.6.15) |
|
|
|
|||||
|
∂zi |
2.6. Интегрирование дифференциальных уравнений |
113 |
||||||||
если использовать соотношение |
|
|
|
|
|
|
|||
ξi |
∂ |
= |
ξi |
∂zj |
|
∂ |
. |
(2.6.16) |
|
∂xi |
∂xi |
∂zj |
|||||||
|
|
|
|
|
|||||
Выбирая в качестве новых переменных любой набор N − 1 |
|||||||||
функционально независимых инвариантов |
|
||||||||
Z1 = I1 (x) , Z2 = I2 (x) , . . . , ZN −1 = IN −1 (x) , |
|
||||||||
а в качестве ZN переменную, которая находится из уравнения |
|||||||||
|
X (z) = 1, |
|
|
(2.6.17) |
получим набор переменных, определяющих группу преобразований вдоль оси ZN .
Рассмотрим интегрирование уравнения (2.6.8) с помощью соответствующей замены переменных. Уравнение (2.6.8) допускает оператор (2.6.10), который в соответствии с вышедоказанной теоремой приводит к группе преобразований переноса после замены
t = t (x, y) , |
u = u (x, y) , |
(2.6.18) |
|||||
которая находится из системы уравнений |
|
||||||
x |
∂u |
− y |
∂u |
= 0, |
(2.6.19) |
||
|
|
|
|
||||
∂x |
∂y |
||||||
x |
∂t |
− y |
∂t |
= 1. |
(2.6.20) |
||
|
|
|
|||||
∂x |
∂y |
||||||
В качестве одной из переменных можно взять инвариант |
|||||||
|
|
u = xy, |
|
(2.6.21) |
|||
а в качестве второй переменной получаем |
|
||||||
|
t = ln x. |
|
(2.6.22) |
114 |
Глава 2. Элементы группового анализа |
Эта замена приводит группу преобразований растяжений к группе преобразований переноса вдоль оси x, поскольку
t = ln x = ln x + a = t + a , |
|
||||||
|
u |
= |
y |
= |
y |
= u . |
|
x 2 |
|
|
|||||
|
|
|
|
x2 |
|
||
Подставляя (2.6.21) в уравнение (2.6.8), приходим к уравне- |
|||||||
нию |
|
|
|
|
|
|
|
|
du |
+ u2 − u − 2 = 0, |
(2.6.23) |
||||
|
dt |
которое в результате интегрирования дает решение
|
u + 1 |
|
|||
ln |
|
|
= 3 t + c, |
|
|
u − 2 |
|
||||
имеющее в старых переменных вид |
|
||||
y = |
2x3 + c |
(2.6.24) |
|||
|
. |
||||
x (x3 − c) |
Рассмотренный способ, указывающий упрощающую замену переменных, обладает универсальностью и может использоваться, как для интегрирования уравнений первого порядка, так и для понижения порядка обыкновенных дифференциальных уравнений второго и выше порядков.
2.7.Группы преобразований для линейного уравнения теплопроводности
Найдем группы преобразований, которые допускаются линейным уравнением теплопроводности
ut = a2uxx. |
(2.7.1) |
2.7. Группы для линейного уравнения теплопроводности 115
Здесь u — температура cреды, которая зависит от x и t, a2 — коэффициент температуропроводности.
Общий вид инфинитезимального оператора можно представить формулой
1 ∂ |
2 ∂ |
|
∂ |
|
|||
X = ξ |
|
+ ξ |
|
+ η |
|
, |
(2.7.2) |
∂x |
∂t |
|
|||||
|
|
|
∂u |
|
где ξ1, ξ2 и η — координаты касательного векторного поля, которые требуется найти. Оператор, продолженный на вторые производные, можно записать в виде:
X = X + ξ1 |
∂ |
+ ξ2 |
∂ |
+ ξ11 |
∂ |
. |
(2.7.3) |
∂ux |
|
|
|||||
2 |
|
∂ut |
∂uxx |
|
Действуя оператором (2.7.3) на уравнение (2.6.1), получаем соотношение
ξ2 − a2ξ11 = 0, |
(2.7.4) |
в котором производные ut и uxx связаны уравнением (2.6.1). Находим выражения для ξ1, ξ2 и ξ11 через ξ1, ξ2и η.
Поскольку
ξ1 = Dx (η) − uxDx ξ1 − utDx ξ2 ,
ξ2 = Dt (η) − uxDt ξ1 − utDt ξ2 ,
ξ11 = Dx (ξ1) − uxxDx ξ1 − uxtDx ξ2 ,
то
ξ1 = ηx + uxηu − uxξ1x − (ux)2 ξ1u − a2uxxξ2x − a2uxuxxξ2u, (2.7.5)
ξ2 = ηt + a2uxxηu − uxξ1t − a2uxuxxξ1u−
(2.7.6)
−a2uxxξ2t − a4 (uxx)2 ξ2u ,
116 |
Глава 2. Элементы группового анализа |
ξ11 = ξ1x + uxξ1u − uxxξ1x − uxuxxξ1u − a2uxxxξ2x − a2uxuxxxξ2u = = ηxx + uxxηu + uxηux − uxxξ1x − uxξ1xx − 2uxuxxξ1u−
−2 (ux)2 ξux1 − a2uxxxξx2 − a2uxxξxx2 − a2 (uxx)2 ξu2 − |
(2.7.7) |
|||||||||||||||||||||||
|
−a2uxuxxxξu2 − a2uxuxxξux2 + uxηxu + (ux)2 ηuu− |
|
||||||||||||||||||||||
− (ux)2 ξxu1 − (ux)3 ξuu1 − a2uxuxxξxu2 − a2 (ux)2 uxxξuu2 − |
||||||||||||||||||||||||
|
−uxxξx1 − uxuxxξu1 − a2uxxxξx2 − a2uxuxxxξu2 . |
|
||||||||||||||||||||||
Подставляя (2.7.6) и (2.7.7) в определяющее уравнение (2.7.4) |
||||||||||||||||||||||||
и приводя подобные члены, получим уравнение: |
|
|
|
|||||||||||||||||||||
|
ηt − a2ηxx + ux −ξt1 + a2ξxx1 − 2a2ηux + |
|
||||||||||||||||||||||
|
|
)2 |
|
|
2ξ1 |
− |
2 |
η |
|
|
3 2ξ1 + |
|
|
|
||||||||||
|
+ (ux |
|
|
|
|
3a |
|
|
|
a |
|
uu |
+ (ux) a |
|
uu |
|
|
|
||||||
|
|
|
|
|
|
ux |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
a2ξ2 |
|
− |
|
u |
+ a4ξ2 |
|
+ |
|
|
|
|
(2.7.8) |
||||
|
+uxx |
− |
+ 2a2ξ1 |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
x |
xx |
|
|
|
|
|
|
|
||
|
+u |
x |
u |
xx |
3a2ξ1 |
|
a2ξ2 |
+ 2a4ξ2 |
|
+ |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
u |
|
|
|
|
ux |
|
|
|
|
|
|
||||
|
4ξ2 |
+ u |
|
u |
xxx |
2a4 |
ξ2 |
2 |
|
|
a4 |
ξ2 |
= 0. |
|
||||||||||
+uxxx2a |
|
x |
|
|
x |
|
|
|
u + (ux) uxx |
uu |
|
|
|
|||||||||||
Откуда находим уравнения для определения координат ξ1, ξ2 |
||||||||||||||||||||||||
и η: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ξ2 |
= ξ2 = 0, |
|
|
|
ξ1 |
= 0, |
|
ξ1 = 0, |
|
|
η |
uu |
= 0, |
(2.7.9) |
||||||||||
u |
x |
|
|
|
|
|
|
|
u |
|
|
|
|
|
xx |
|
|
|
|
|
|
|||
|
ηt = a2ηxx, |
|
|
ξt2 = 2ξx1 , |
ξt1 = −2a2ηux. |
(2.7.10) |
||||||||||||||||||
Из уравнений (2.7.9) находим |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
ξ2 = Ψ (t) , |
|
|
|
|
|
|
|
(2.7.11) |
|||
|
|
|
|
|
|
|
|
ξ1 = ϕ1 (t) + xϕ2 (t) , |
|
|
|
|
|
(2.7.12) |
2.7. Группы для линейного уравнения теплопроводности 117
η = b (x, t) + ub1 (x, t) . |
(2.7.13) |
|||||||||||||||
Подставив полученные соотношения во второе и третье урав- |
||||||||||||||||
нение (2.7.10), имеем |
1 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
Ψt, |
|
|
|
|
||||||||
ϕ2 (t) = |
|
|
|
|
|
|
|
(2.7.14) |
||||||||
2 |
|
|
|
|
||||||||||||
|
|
|
x2 |
|
|
|
|
|
|
x2 |
|
|||||
b1 (x, t) = ϕ (t) − |
|
|
|
|
· ϕ1t − |
|
|
Ψtt. |
||||||||
2a2 |
8a2 |
|||||||||||||||
Принимая во внимание эти выражения, из первого уравнения |
||||||||||||||||
(2.7.10) приходим к соотношениям |
|
|
|
|
|
|
|
|
|
|||||||
bt = a2bxx, |
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||
ϕt = − |
|
Ψtt, |
|
|
|
(2.7.15) |
||||||||||
4 |
|
|
|
|||||||||||||
ϕ1tt = 0, |
|
|
|
Ψttt = 0. |
|
|||||||||||
Откуда находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Ψ (t) = B2 + 2B3t + B6t2, |
|
|||||||||||||||
ϕ1 (t) = B1 + B4t, |
|
|
|
(2.7.16) |
||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
ϕ (t) = B5 − |
|
B6t. |
|
|
|
|
||||||||||
2 |
|
|
|
|
||||||||||||
Таким образом, координаты касательного векторного поля ξ1, |
||||||||||||||||
ξ2и η принимают вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ξ1 = B1 + B4t + (B3 + B6t) x, |
|
|||||||||||||||
ξ2 = B2 + 2B3t + B6t2, |
|
|
|
(2.7.17) |
||||||||||||
1 |
|
|
|
|
|
|
B4x |
B6x2 |
||||||||
|
|
|
|
|
|
|
||||||||||
η = b (x, t) + B5u − |
|
B6t + |
|
|
+ |
|
|
u, |
||||||||
2 |
|
2a2 |
|
4a2 |
где B1, . . . , B6 — произвольные постоянные, b (x, t) — произвольная функция, которая удовлетворяет уравнению теплопроводности (2.7.1).
118 |
Глава 2. |
Элементы группового анализа |
|
||||
Полагая Bi = 0 |
(i = 1, . . . , 6) |
из (2.7.17) находим беско- |
|||||
нечное число операторов |
|
|
|
|
|
||
|
|
X |
|
= b (x, t) |
∂ |
, |
(2.7.18) |
|
|
∞ |
|
||||
|
|
|
|
∂u |
|
которые допускаются линейным уравнением теплопроводности (2.7.1). Из оператора (2.7.18) следует семейство преобразований
u (x, t) = a b (x, t) + u (x, t) , x = x, t = t ,
которое выражает свойство линейности для решений уравнения (2.7.1). Полагая b (x, t) = 0 и последовательно одну из постоянных Bi = 1, а остальные равными нулю, найдем еще шесть
операторов:
X1 |
= |
∂ |
, |
|
|
X2 = |
|
|
∂ |
, |
|
||||||
|
|
|
∂t |
|
|||||||||||||
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
X4 = t |
|
∂ |
xu |
|
∂ |
|
|||||||||
|
|
|
|
|
− |
|
|
|
|
, |
|||||||
|
|
|
∂x |
2a2 |
∂u |
||||||||||||
|
|
|
|
|
∂ |
|
|
∂ |
|
− |
|||||||
|
X6 = xt |
|
+ t2 |
|
|
|
|||||||||||
|
∂x |
∂t |
|
|
X3 |
= x |
∂ |
|
+ 2t |
∂ |
, |
|||
|
∂x |
∂t |
||||||||
|
|
|
|
|
|
|
|
|||
|
|
X5 = u |
∂ |
, |
|
(2.7.19) |
||||
|
|
|
|
|
||||||
|
|
|
|
|
|
∂u |
|
|
||
t |
|
x2 |
∂ |
|
|
|||||
|
+ |
|
u |
|
. |
|
|
|||
2 |
4a2 |
∂u |
|
|
Эти операторы содержат важную информацию о симметрии уравнения теплопроводности.
Используя координаты касательных векторных полей в операторах (2.7.19), можно найти группы преобразований и инварианты, соответствующие этим операторам.
Например, решая уравнения Ли
dx = x , x (a = 0) = x , da
dt = 2t , t (a = 0) = t da
с координатами, соответствующими оператору X3, находим группу преобразований растяжения, которая допускается уравнением (2.6.1):
x = eax, t = e2at, u = u. |
(2.7.20) |
2.8. Группы для нелинейного уравнения теплопроводности 119
Откуда следуют два инварианта
I1 |
= |
x 2 |
= |
x2 |
, I2 = u |
= u. |
(2.7.21) |
|
t |
t |
|||||||
|
|
|
|
|
|
Эти выражения приводят к хорошо известному факту, что линейное уравнение теплопроводности (2.7.1) имеет инвариантное решение в виде
u = C1 w(z), |
z = |
C2 x |
, |
(2.7.22) |
1 |
||||
|
|
t 2 |
|
где C1 и C2 — произвольные постоянные.
Аналогично, используя операторы (2.7.19), можно найти и другие группы преобразований, которые соответствуют уравнению (2.7.1).
2.8.Группы преобразований для нелинейного уравнения теплопроводности
Для нелинейного уравнения теплопроводности
∂u |
= |
∂ |
un |
∂u |
(2.8.1) |
|
∂t |
∂x |
∂x |
||||
|
|
|
группа преобразований находится аналогично, как это было сделано для уравнения (2.7.1).
Пусть оператор X представляется в виде
X = ξ1 |
∂ |
+ ξ2 |
∂ |
+ η |
∂ |
, |
(2.8.2) |
|
∂x |
∂t |
∂u |
||||||
|
|
|
|
|
где ξ1, ξ2и η — координаты касательного векторного поля, которые следует найти. Предполагая, что
X = X + ξ1 |
∂ |
+ ξ2 |
∂ |
+ ξ11 |
∂ |
, |
(2.8.3) |
∂ux |
|
|
|||||
2 |
|
∂ut |
∂uxx |
|
120 |
Глава 2. Элементы группового анализа |
и действуя этим оператором на уравнение (2.8.1), получаем соотношение
ξ2 − n (n − 1) (ux)2 un−2η − 2nun−1uxξ1−
(2.8.4)
−nun−1uxxη − unξ11 = 0.
Вычисляя по формулам продолжения компоненты ξ1, ξ2и ξ11 и подставляя выражения для них в определяющее уравнение (2.8.4), получаем четыре независимых оператора
X1 = |
∂ |
, X2 = |
|
∂ |
, |
X3 = 2t |
∂ |
+ x |
∂ |
, |
||||
∂t |
|
|
|
∂t |
∂x |
|||||||||
|
|
|
∂x |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.8.5) |
|
|
X4 = |
nx ∂ |
+ u |
∂ |
|
|
|
|
|||||
|
|
|
|
|
|
|
. |
|
|
|
|
|||
|
|
|
2 |
∂x |
∂u |
|
|
|
|
В случае n = −4/3 уравнение (2.8.1) допускает дополнительный к (2.8.5) оператор в виде
2 ∂ |
∂ |
(2.8.6) |
||||
|
|
|
|
|
||
X5 = −x ∂x + 3xu |
∂u . |
|||||
|
Уравнение (2.8.1) при n = 0 переходит в линейное уравнение (2.7.1), имеющее группу преобразований, представленную в предыдущем параграфе.
Заметим, что поскольку операторы группы преобразований линейные, то, кроме “основных” операторов (2.8.5), уравнением (2.8.1) допускается и произвольная сумма их, умноженная на постоянные коэффициенты. Например, уравнение (2.8.1) допускает оператор
X = |
1 ∂ |
− |
∂ |
, |
(2.8.7) |
||
c0 |
|
∂t |
∂x |
который ведет к инвариантам
ξ = x − c0t, u (x, t) = v,