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

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

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

 

5.4. Метод простейших уравнений

281

Подставляя (5.4.60) в (5.4.53), получаем B2(1) = 12,

B2(2) = 0.

Рассмотрим первый случай B2 = B2(1) = 12. Имеем

 

 

 

c = 0,

 

 

B1 = 0,

A0 = 8 b − α,

(5.4.62)

A1 = 9a,

A0 = 11b − α,

C0 = β α,

(5.4.63)

 

 

 

 

A(1)

= 12,

A(2) = 0.

(5.4.64)

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Полагая A2 = A(1), находим

 

 

 

 

 

 

a = 0,

C1 =

1

β α2 8

β b2 + 24 β d.

(5.4.65)

 

 

2

Решение в виде перидических волн определяются формулой

 

 

 

 

y(z) = α 4 b −

12 d

(5.4.66)

 

 

 

 

 

 

,

 

 

 

 

Q(z)2

где Q(z) — общее решение уравнения

 

 

 

 

 

 

 

 

Qz2 − Q4 − b Q2 − d = 0.

(5.4.67)

Полагая B2 = 0 в преобразовании (5.4.60), получаем

B1 = 0,

 

A2 = 6,

D1 = ±6

A1 = 3 a,

(5.4.68)

 

 

 

A0 = α − b,

C0 = β α,

(5.4.69)

C1 =

1

 

β α2 +

3

β ac − 6 β d −

1

β b2.

(5.4.70)

 

2

 

 

2

2

Решение в этом случае имееет вид

 

 

 

 

y(z) = α − b − 3 a Q − 6 Q2 + 6 Qz ,

(5.4.71)

где Q(z) удовлетворяет уравнению

 

 

 

 

Qz2 − Q4 − a Q3 − b Q2 − c Q − d = 0.

(5.4.72)

282 Глава 5. Методы построения точных решений

Уравнение (5.4.59) при этом запишется в виде

yzzz + β yzz + α yz +

1

β y2 + y yz + β α y (z) +

2

 

1

 

3

 

 

1

(5.4.73)

 

 

 

 

 

+

β α2 +

β ac − 6 β d −

β b2 = 0.

 

 

 

2

2

2

Решения (5.4.66) и (5.4.71) уравнения (5.4.53) являются периодическими волнами. Более детальное обсуждение применения метода простейших уравнений можно найти в работе [124].

5.5.Применение многоугольников при построении точных решений

Одной из трудностей применения метода простейших уравнений для построения точных решений нелинейных дифференциальных уравнений является выбор простейшего уравнения. В ряде случаев в качестве простейших уравнений можно взять хорошо известные точно решаемые нелинейные обыкновенные дифференциальные уравнения: уравнение Риккати (5.4.5), уравнение для эллиптической функции Якоби (5.4.11) или уравнение для эллиптической функции Вейерштрасса (5.4.12). Однако во многих случаях эта трудность становится препятствием применения метода простейших уравнений, в преодолении которой может существенно помочь применение многоугольников, которые часто называют многоугольниками Ньютона [11, 12, 130].

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

5.5. Применение многоугольников

283

Целиком выпуклая ломаная была рассмотрена Пюизе. Работа самого Пюизе не содержала ссылки на Ньютона, и остается неясным ответ на вопрос, был ли знаком Пюизе с работой Ньютона. Возможно, поэтому многие исследователи многоугольник Ньютона часто называют многоугольником Ньютона—Пюизе [3].

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

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

Mn[y(z), yz (z), yzz (z), . . . , z] = 0.

(5.5.1)

Каждому члену (моному) дифференциального уравнения (5.5.1) ставится в соответствие точка на плоскости по следующему правилу [11, 12, 130]

C1 zq1 yq2 (q1, q2), C2

dk y

(−k, 1),

(5.5.2)

dzk

где C1, C2 — произвольные постоянные. При умножении мономов их координаты складываются. Множество точек на плоскости, сооветствующее дифференциальному уравнению, образует так называемый носитель. Соединяя все внешние точки носителя, получаем многоугольник Ньютона на плоскости. Таким образом, каждому обыкновенному дифференциальному уравнению (5.5.1) ставится в соответствие многоугольник L1 на плоскости [11, 12, 127, 128, 129].

Пусть решение y(z) исходного уравнения выражается через решения Y (z) простейшего уравнения

y(z) = F (Y (z), Yz (z), . . . , z).

(5.5.3)

284 Глава 5. Методы построения точных решений

Подстановка (5.5.3) в уравнение (5.5.1) ведет к новому дифференциальному уравнению, которому также можно поставить в соответствие многоугольник Ньютона на плоскости. Обозначим этот многоугольник L2. Анализируя многоугольник L2, можно построить многоугольник Ньютона L3, соответствующий про-

 

(1)

 

˜(1)

стейшему уравнению. Пусть ребро Γ1

относится к L2, ребро Γ1

(1)

˜

(1)

имеют равные внешние

к L3. Можно доказать, что если Γ1

и Γ1

нормальные векторы, то соответствующие асимптотики совпадают вплоть до значения коэффициентов [11, 12]. Под внешним нормальным вектором мы имеем в виду вектор, который направлен вне многоугольника. Следовательно, пригодный многоугольник L3 имеет все или несколько ребер, параллельных многоугольнику L2. Кроме того, в случае, когда многоугольник L3 передвигается по плоскости, его ребра должны перекрывать носитель преобразованного уравнения.

Предположим, что многоугольник L3 построен. Тогда анали-

зируя его, можно написать простейшее уравнение

 

Em(Y (z), Yz (z), . . . , z) = 0.

(5.5.4)

Если удается найти простейшее уравнение, такое, что выпол-

няется соотношение

 

ˆ

(5.5.5)

Mn(F (Y, Yz , . . . , z)) = R Em(Y, Yz , . . . , z),

ˆ

где R — дифференциальный оператор, то это означает, что для любого решения Y (z) простейшего уравнения (5.5.4) существует решение (5.5.3) исходного нелинейного дифференциального уравнения (5.5.1).

При выборе простейшего уравнения естественным является требование о том, что порядок простейшего уравнения выбирается меньшим, чем порядок исходного уравнения, соответствующего многоугольнику L2. Самыми важными простейшими уравнениями являются такие уравнения, общие решения которых не имеют подвижных критических особых точек. Если общее (или

5.5. Применение многоугольников

285

частное) решение простейшего решения находится, то для изучаемого уравнения по формуле (5.5.3) находятся точные решения.

Применение многоугольников Ньютона для построения точных решений нелинейных дифференциальных уравнений можно сформулировать в виде следующих шести шагов.

Первый шаг. Построение многоугольника L1, соответствующего исходному уравнению.

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

Третий шаг. Построение многоугольника L2 соответствующего преобразованному уравнению.

Четветый шаг. Построение многоугольника L3, характеризующего простейшее уравнение.

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

Шестой шаг. Определение коэффициентов преобразования (5.5.3) и параметров простейшего уравнения.

В случае, когда полюс общего решения простейшего уравнения совпадает с полюсом исходного уравнения, простейшее уравнение часто выбирается без преобразования исходного уравнения (5.5.3). В этом случае формула (5.5.3) становится тождественной подстановкой: y(z) ≡ Y (z).

5.5.1.Точные решения обобщенного уравнения Курамото—Сивашинского

Продемонстрируем применение метода многоугольников для построения точных решений обобщенного уравнения Курамото— Сивашинского [130]

ut + α um ux + δ uxx + β uxxx + γ uxxxx = 0.

(5.5.6)

286 Глава 5. Методы построения точных решений

Вводя новые переменные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

x = x

 

 

 

δ2

 

 

 

 

 

δ

 

 

 

β

 

 

γ δ

 

 

, t = t

 

, u = u,

σ =

 

,

α =

 

 

,

γ

γ

δ

γ δ

получаем уравнение (штрихи далее опускаются)

 

(5.5.7)

 

 

 

 

ut + α um ux + uxx + σ uxxx + uxxxx = 0.

(5.5.8)

Используя переменные бегущей волны

 

 

 

 

 

 

 

 

 

u(x, t) = y(z),

z = x − C0 t,

 

 

(5.5.9)

и интегрируя полученное уравнение по переменной z, приходим к уравнению

 

α

 

M3[y] = yzzz + σ yzz + yz − C0 y +

 

ym+1 = 0,

(5.5.10)

m + 1

где постоянная интегрирования полагается равной нулю. Точные решения уравнения (5.5.10) находятся в соответствии со следую-

щей теоремой.

Теорема 5.1. Пусть y(z) — решение уравнения

yz =

 

 

 

 

9 α

 

1

 

 

 

 

 

 

 

3

m+3

 

 

 

 

 

 

 

y 3

 

2 m3 + 11 m2 + 18 m + 9

 

(5.5.11)

 

 

 

 

3

y.

 

 

 

 

 

 

 

 

 

 

 

 

 

2 m

2

 

 

 

 

 

 

 

 

 

+ 18 m + 27

 

 

 

 

 

Тогда y(z) является также решением уравнения (5.5.10) при

m = 0 m = 1 m =

 

 

 

m = 3 m = 9 ± 3

 

 

3

 

3

,

,

− ,

2 ,

 

− ,

 

 

2

 

 

 

 

C0 =

 

 

 

3 (2 m2 + 9 m + 9)

,

 

 

 

 

 

 

(2 m2 + 18 m + 27)3/2

 

 

(5.5.12)

 

 

 

 

 

 

 

3 (3 + m)

 

 

 

 

 

 

 

σ = ±

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 m

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 18 m + 27

 

 

 

 

5.5. Применение многоугольников

287

Доказательство. Мономам уравнения (5.5.10) соответствуют следующие точки на плоскости: M1 = (3, 1),

M2 = (2, 1), M3 = (1, 1), M4 = (0, 1) и M5 = (0, m + 1). Многоугольник L1, соответствующий уравнению (5.5.10), является большим треугольником на рис. 5.7 (здесь левый треугольник соответствует m > 0 и, правый треугольник соответствует m < 0). Подходящий многоугольник L3 может быть найден без преобразования исходного уравнения это малый треугольник на рис.

4.10, который содержит точки: Q1 = (1, 1), Q2 = M4 = (0, 1), и

Q3 = (0, 1 + m3 ).

Рис. 5.7. Многоугольники Ньютона, соответствующие дифференциальным уравнениям (5.5.10) и (5.5.13) при m > 0 (слева) и при m < 0 (справа)

Простейшее уравнение, соответствующее многоугольнику L3, можно представить в виде

m+3

− B y = 0,

(5.5.13)

E1[y] = yz − A y 3

288 Глава 5. Методы построения точных решений

где A и B — параметры уравнения, которые следует найти. Общее решение уравнения (5.5.13) определяется формулой

 

B

m z

 

A

3

 

 

 

 

m

 

y(z) = C1 exp

 

 

.

(5.5.14)

 

3

B

Подставляя

 

 

 

 

 

 

 

 

 

 

m+3

+ B y

(5.5.15)

yz = E1[y] + A y 3

в уравнение (5.5.10) и приравнивая нулю коэффициенты при одинаковых степенях y(z), приходим к алгебраическим уравнениям относительно параметров A, B, C0 и σ

 

A3 (2 m + 3) (m + 3) (m + 1) + 9 α = 0,

(5.5.16)

 

2

(m + 1)

B3 + σ B2 + B − C0

= 0,

(5.5.17)

 

A (3 + m) (m + 1)

(B m + 3 B + σ) = 0,

(5.5.18)

(m + 1)

2m2 + 18σ2m +

27σ2 9m2 54m − 81

= 0. (5.5.19)

Решая уравнения (5.5.16), (5.5.17), (5.5.18) и (5.5.19), получаем

A1 = (9 α)1/3 ( m3 + 11 m2 + 18 m + 9)1/3 ,

1

 

 

(

 

 

3

 

 

A2,3 = (9α)1/3 '

 

i

 

 

 

2

2

 

m3 + 11m2 + 18m + 9 1/3 ,

 

 

 

C0 = B B2 + σ B + 1 ,

B =

 

σ

,

 

 

3 + m

σ = ±

3 (3 + m)

 

.

 

 

 

2 m2 + 18 m + 27

 

(5.5.20)

(5.5.21)

(5.5.22)

5.5. Применение многоугольников

289

Таким образом, мы нашли значения параметров уравнения

(5.5.14), зависимость (5.5.12) и соотношение

 

ˆ

(5.5.23)

M3[y] = R E1[y],

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где R — дифференциальный оператор.

 

 

 

 

 

 

 

 

Общее решение уравнения (5.5.12) можно представить в виде

 

 

 

 

 

 

 

 

 

 

m z

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

y(z) = C + C1 exp ±

 

 

 

 

 

 

, (5.5.24)

2 m2 + 18 m + 27

 

 

 

где C1 — произвольная постоянная, C определяется выражением

 

 

 

 

3

 

 

 

 

.

 

 

 

 

 

C =

±

9 α

2 m2 + 18 m + 27

 

 

(5.5.25)

 

 

 

3 3

 

 

 

 

 

 

 

 

(m + 3) (2 m2 + 5 m + 3)

 

 

 

Скорость уединенной волны (5.5.24) при m

 

 

 

 

 

0 стремится к

 

 

 

 

 

 

 

 

 

 

 

 

C0 = 1/ 27. Однако при m → ∞ скорость C0 0.

Полагая m = 1 в (5.5.24), получим решение в виде уединенной

волны

 

30

 

 

1

 

 

 

 

±

 

47

 

(

 

'

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

y(z) =

 

47

 

225 α

+ C

 

exp

 

 

x − C0 t

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.5.26)

 

 

 

 

C0 =

 

 

60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

47

47

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а также значение

параметра

 

σ: σ

=

±12/

47

. Это решение

уравнения Курамото—Сивашинского представляет собой кинк [41, 45, 113, 120, 124]. Некоторые другие решения обобщенного уравнения Курамото—Сивашинского (5.5.6) получены в работе [131].

5.5.2.Автомодельные решения уравнения Кортевега— де Вриза пятого порядка

Используя многоугольники Ньютона, найдем решения уравнения Кортевега—де Вриза пятого порядка [130]

290 Глава 5. Методы построения точных решений

ut + uxxxxx 10 u uxxx 20 uxuxx + 30 u2 ux = 0. (5.5.27)

Решения уравнения (5.5.27) будем искать, используя автомодельные переменные

1

y(z), z =

x

 

 

u(x, t) =

 

 

,

(5.5.28)

(5t)2/5

(5t)1/5

 

 

 

 

где y(z) удовлетворяет нелинейному обыкновенному дифференциальному уравнению

yzzzzz 10 y yzzz 20 yz yzz + 30 y2 yz − z yz 2 y = 0. (5.5.29)

Результаты можно сформулировать в виде следующей теоре-

мы.

Теорема 5.2. Пусть Y (z) — решение уравнения

Yzzzzz 40 Y Yzz Yz 10 Y 2 Yzzz 10 Yz3+

(5.5.30)

+30 Y 4Yz − Y − z Yz = 0.

 

Тогда

 

y(z) = ± Yz − Y 2

(5.5.31)

является решением уравнения (5.5.29).

Доказательство. Мономам уравнения (5.5.29) соответствуют точки на плоскости: M1 = (5, 1), M2 = (3, 2), M3 = (3, 2), M4 = (1, 3), M5 = (0, 1) и M6 = (0, 1). Носитель уравнения (5.5.29) определяется точками: Q1 = M1 = (5, 1), Q2 = M4 = = (1, 3), Q3 = M5 = M6 = (0, 1) и Q4 = M2 = M3 = (3, 2). Соответствующий многоугольник L1 представлен на рис. 5.8.

Общее решение уравнения (5.5.29) имеет полюс второго по-

рядка. Поэтому решение будем искать в виде

 

y = A1Y + A2Y 2 + A3Yz ,

(5.5.32)