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

ДУ РФ

.pdf
Скачиваний:
27
Добавлен:
27.03.2015
Размер:
385.61 Кб
Скачать

Ре ш е н и е. Матрица системы равна

 

 

 

 

 

4

 

2

5

 

 

 

 

 

 

 

 

6

 

1

6

.

 

 

 

 

 

 

A = 8

3

9

 

 

Характеристическое уравнение

 

 

= 0

 

 

 

A λE =

4

6

1 λ

 

6

или λ3

 

4λ2 + 5λ 2 = 0

| − |

 

 

λ

2

 

5

 

 

 

 

8

− − 3

9 λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет корни λ1 = 2, λ2 = λ3 = 1. Ищем решения, соответствующие этим корням.

λ1 = 2

Система для нахождения собственного вектора имеет вид

 

6

2

5

 

γ11

 

 

 

6

3

6

 

γ12

 

= 0.

8

3

7

γ13

 

Для нахождения вектора γ1 вычислим алгебраические дополнения элементов первой строки определителя этой системы

6 2 5

6 3 6 .

8 3 7

Получим

3

7

 

=

 

7

 

γ11 =

= 3, γ12

8

= 6,

 

3

6

 

 

 

6

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 3

 

 

 

 

 

8 3

 

 

 

γ13 =

= 6

или (сократив на 3) γ11 = 1, γ12 =

2, γ13 = 2.

 

Итак,

 

 

 

 

 

 

 

 

 

 

λ1 = 2,

 

γ1 = (1, −2, 2),

 

x1 = e2t, y1 = 2 e2t, z1 = 2 e2t.

λ2 = λ3 = 1

 

 

 

 

 

 

Имеем

6 2

6

γ22

= 0.

 

5

2

 

5

γ21

 

 

8

3

 

8 γ23

61

Легко видеть, что в этом случае n = 3, m = 2, r = 2 и n − m = r. В отличие от предыдущего примера решение, соответствующее λ2 = λ3 = 1, ищем методом неопределенных коэффициентов.

x = (A + Bt) et, y = (C + Dt) et, z = (E + F t) et.

Подставляя предполагаемые решения в систему дифференциальных уравнений и сокращая на et, имеем

A + B + Bt = 4A − 4Bt + 2C + 2Dt + 5E + 5F t,

C + D + Dt = 6A + 6Bt − C − Dt − 6E − 6F t,

E + F + F t = 8A − 8Bt + 3C + 3Dt + 9E + 9F t.

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

 

 

 

 

6B

2D 6F = 0,

 

 

 

 

5B + 2D + 5F = 0,

 

 

 

 

 

 

 

 

 

 

 

 

8B + 3D + 8F = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5A B + 2C + 5E = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6A − 2C − D − 6E = 0,

 

 

 

 

 

8A + 3C + 8E F = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

,

 

 

 

 

 

Положим E = C2

 

F = C3

 

где C2 и C3 — произвольные постоянные. Тогда B = C3,

D = 0, A = C2 + C3, C = 3C3

и решение соответствующее кратным корням

λ2 = λ3 = 1, имеет вид

 

 

 

 

 

 

x = (C2 + C3 + C3t) et, y = 3C3 et, z = (C2 + C3t) et.

Теперь нетрудно составить общее решение исходной системы дифференциальных

уравнений

x = C1 e2t + (C2 + C3 + C3t) et,

y = 2C1 e2t + 3C3 et,

z = 2C1 e2t + (C2 + C3t) et.

Задачи для практических занятий

13.1.

y˙ = x + y

z,

13.2.

y˙ = x + y,

− z,

 

 

 

 

x˙ = x + z y,

 

x˙ = x − y

 

 

 

 

 

 

 

 

 

 

z˙ = 2x y.

 

z˙ = 3x + z.

 

 

 

 

 

 

 

 

±

 

 

 

 

 

 

 

 

(λ1 = 1, λ2 = 2, λ3 = 1).

 

(λ1 = 1, λ2,3 = 1

 

2i).

62

13.3. y˙ = x + 2y z,

13.4.

y˙ = x + y

 

z,

 

 

x˙ = 4x − y − z,

 

 

x˙ = x − y + z,

 

 

 

 

 

 

 

z˙ = x y + 2z.

).

z˙ = 2z y.

 

).

 

,

 

 

 

 

,

 

 

 

 

 

 

 

 

 

(λ1 = 2 λ2 = λ3 = 3

 

(λ1 = λ2 = 1 λ3 = 2

 

Домашнее задание

13.5.

 

x˙ = x 2y z,

 

 

 

 

13.6.

x˙ = 2x + y,

 

z,

 

 

y˙ = y

x +z,

 

 

 

 

 

y˙ = x + 3y

 

 

 

 

 

 

 

z.

 

 

 

 

 

 

 

 

 

 

 

 

x.

 

 

 

z˙ = x

 

 

 

 

 

 

 

z˙ = 2y + 3z

 

 

 

(λ1

= 0

,

 

= 2

,

λ3

= 1

).

 

(λ1

= 2

,

λ2,3

 

 

 

 

 

 

λ2

 

= 3 i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.7.

 

x˙ = y 2x − 2z,

 

 

 

 

13.8.

x˙ = y − 2z − x,

 

 

y˙ = x

2y + 2z,

 

 

 

 

y˙ = 4x + y,

 

 

 

 

 

 

 

 

 

 

 

3y + 5z.

 

 

 

 

 

 

 

 

 

 

z.

 

 

 

z˙ = 3x

 

 

 

 

z˙ = 2x + y

 

 

 

(λ1

= 3

,

 

 

 

= 1

).

 

 

(λ1

= 1

,

 

 

= 1

).

 

 

λ2 = λ3

 

 

λ2 = λ3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответы:

 

t

e

e

 

e

,

 

 

e t

t

 

 

t

13.1. x = C1 t + C2

 

2t + C3

−t

 

y = C1

t

3C3 e−t, z = C1 et + C2 e2t

 

5C3 e−t.

13.2.

x = e (

2C2 sin 2t +

2C3 cos 2t)

 

y = C1 e

+ e (C2 cos 2t + C3 sin 2t), z = C1 e +

e

(3C2 cos 2t + 3C3 sin 2t). 13.3. x = C1 e

+ (C2 + C3) e

, y = C1 e

+ C2 e

, z =

t

 

 

 

 

 

,

 

2t

3t

 

2t

3t

 

C1 e2t + C3 e3t. 13.4. x = (C1 + C2t) et + C3 e2t, y = (C1 2C2 + C2t) et, z = (C1 −C2 + C2t) et + C3 e2t. 13.5. x = C1 + 3C2 e2t, y = 2C2 e2t + C3 e−t, z = C1 + C2 e2t 2C3 e−t.

13.6. x = C1 e2t + e3t(C2 cos t + C3 sin t), y = e3t((C2 + C3) cos t + (−C2 + C3) sin t), z = C1 e2t + e3t((2C2 − C3) cos t + (C2 + 2C3) sin t). 13.7. x = C1 e3t + (C2 2C3) e−t, y = −C1 e3t + C2 e−t, z = 3C1 e3t + C3 e−t. 13.8. x = (C2 + C3t) e−t, y = 2C1 et

(2C2 + C3 + 2C3t) e−t, z = C1 et (C2 + C3 + C3t) e−t.

Занятие 14. Решение дифференциальных уравнений с помощью степенных рядов

Не вдаваясь в подробности исследований сходимости рядов и существования решения, покажем этот метод на примерах.

П р и м е р ы

1) Для уравнения

(1 − x2)y − xy − y = 0

найти линейно-независимые решения (фундаментальную систему) y1(x) и y2(x) в виде степенных рядов по степеням x и построить общее решение.

63

Р е ш е н и е. Ищем решение уравнения в виде ряда

 

 

 

 

 

 

 

 

y =

ck xk ,

 

 

 

k=0

 

Тогда

 

 

 

 

 

 

 

y =

kck xk−1,

y = k(k − 1)ckxk−2.

 

k=1

k=2

 

Подставляем выражения для y, y и y через ряды в уравнение.

 

 

 

 

k(k − 1)ckxk−2 − k(k − 1)ckxk − kck xk − ck xk = 0.

k=2

k=2

k=1

k=0

Для удобства ”подравняем“ показатели степеней у переменной x. С этой целью заменим индекс суммирования в первой сумме, полагая k − 2 = m, а в остальных суммах положим k = m.

 

 

(m + 2)(m + 1)cm+2xm − m(m − 1)cmxm

m=0

m=2

 

 

− mcmxm − cmxm = 0.

m=1

m=0

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

x0 :

1 · 2c2 − c0 = 0 c2 =

 

c0

 

 

 

 

,

 

 

2!

2c1

 

x1 :

2 · 3c3 − c1 − c1 = 0

c3 =

,

 

3!

xm :

(m + 2)(m + 1)Cm+2 − m(m − 1)cm − mcm − cm = 0 (m ≥ 2).

Из последнего равенства выводим рекуррентное соотношение

 

1 + m2

cm+2 =

(m + 1)(m + 2) cm (m 2).

1. Для нахождения решения y1(x) положим c0 = 1, c1 = 0. Тогда

c3 = c5 = · · · = c2m+1 = · · · = 0

в силу рекуррентного соотношения. Далее имеем (снова используя рекуррентное соотношение)

c2 =

1

 

,

 

c4 =

1 + 22

c2 =

1 + 22

,

 

2!

 

3 · 4

4!

 

 

1 + 42

(1 + 22)(1 + 24)

 

· · ·

c6 =

 

 

 

 

 

c4 =

 

 

 

,

 

5

·

6

6!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

64

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

1 + 22

 

4

(1 + 22)(1 + 24)

 

6

y1(x) = 1 +

 

x

 

+

 

x

+

 

x

+ · · ·

2!

 

4!

6!

2.

Для нахождения решения y2(x) положим c0 = 0, c1 = 1. Тогда

 

c0 = c2 = c4 = · · · = c2m = · · · = 0,

 

 

 

 

 

 

2

 

 

 

1 + 32

 

 

2(1 + 32)

 

 

 

 

c3 =

 

 

,

c5 =

 

 

 

c3 =

 

 

,

 

 

 

 

3!

4 · 5

5!

 

 

 

 

 

c7 =

1 + 52

 

2(1 + 32)(1 + 52)

 

 

 

 

 

 

 

 

 

 

c5 =

 

 

 

 

 

,

 

 

 

 

 

 

 

6 · 7

 

 

 

 

7!

 

 

 

 

 

2

 

3

 

 

 

2(1 + 32)

 

5

 

 

2(1 + 32)(1 + 52)

7

+ · · ·

 

y2(x) = x +

 

x

 

+

 

 

 

x

 

+

 

 

 

 

x

 

 

3!

 

5!

 

 

 

 

7!

 

 

 

3.

Общим решением уравнения будет

 

 

 

 

 

 

 

 

 

 

y= C1y1(x) + C2y2(x).

2)С помощью степенных рядов найти фундаментальную систему решений уравнения

x2y − x2y + (x − 2)y = 0.

Р е ш е н и е. Данное уравнение имеет особенность в точке x = 0. Поэтому его решение будем искать в виде обобщенного степенного ряда

y = xσ ckxk ,

k=0

где показатель σ подлежит определению. Имеем

 

 

y = ck xk+σ , y = (k + σ)ckxk+σ−1,

k=0

k=0

 

 

(k + σ)(k + σ − 1)ckxk+σ−2.

y =

k=0

 

Подставляем эти выражения для y, y и y в уравнение:

 

 

(k + σ)(k + σ − 1)ckxk+σ (k + σ)ckxk+σ+1+

k=0

k=0

 

 

+ ckxk+σ+1 2ckxk+σ = 0,

k=0

k=0

 

 

[(k + σ)(k + σ − 1) 2]ck xk+σ (k + σ − 1)ckxk+σ+1 = 0.

k=0

k=0

65

Заменим в первой сумме k на m, а во второй положим k + 1 = m:

 

 

 

[(m + σ)(m + σ − 1) 2]cmxm+σ (m + σ − 2)cm−1xm+σ = 0.

m=0

m=1

Откуда следует

(σ(σ − 1) 2)c0 = 0,

[(m + σ)(m + σ − 1) 2]cm = (m + σ − 2)cm−1.

Полагая в первом равенстве c0 = 1, находим показатели σ = 2 и σ = 1. Второе

равенство будем использовать как рекуррентное.

 

 

 

 

При σ = 2 из рекуррентного равенства получаем

 

 

 

 

 

 

cm =

cm−1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

m + 3

 

 

 

 

 

 

 

Откуда последовательно находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c0 = 1, c1 =

1

, c2

=

 

 

1

 

,

 

c3 =

 

1

 

и т.д.

 

 

 

·

5

 

 

· 5 · 6

 

4

 

 

4

 

 

 

4

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1(x) = x2 1 +

x

 

2

 

 

 

 

x3

 

 

 

.

 

+

x

+

 

 

+ · · ·

4

4 · 5

4 · 5 · 6

При σ = 1 рекуррентное равенство принимает вид

 

m(m − 3)cm = (m − 3)cm−1.

 

Если m = 3, то оно превращается в рекуррентное соотношение

cm =

cm−1

.

 

 

 

m

 

 

Откуда следует

1

 

c0 = c1 = 1, c2 =

.

 

2

При m = 3 выполняется тождество m(m −3)cm = (m −3)cm−1. Поэтому в качестве коэффициента c3 можно взять любое значение. Полагая c3 = 0, из рекуррентного соотношения получаем, что cm = 0 для всех m ≥ 3. В этом случае

y2

(x) = x 1 + x +

2 .

 

1

 

x

Подстановкой в уравнение убеждаемся, что это будет второе фундаментальное решение.

Общее решение исходного уравнения имеет вид y = C1y1(x) + C2y2(x).

66

Задачи для практических занятий

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

14.1.(1 − x2)y − 4xy − 2y = 0.

14.2.(1 − x)y − 2y + y = 0.

14.3.y + y sin x = 0.

Найти линейно независимые решения следующих уравнений в виде обобщенных степенных рядов.

14.4.xy + 2y + xy = 0.

14.5.2x2y + (3x − 2x2)y − (x + 1)y = 0.

Домашнее задание

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

14.6.y − xy − 2y = 0.

14.7.(x2 + 1)y + 5xy + 3y = 0.

14.8.xy + y ln(1 − x) = 0.

Найти линейно независимые решения в виде обобщенных степенных рядов.

14.9. 9x2y − (x2 2)y = 0.

Ответы:

14.1. y1

= 1 + x2

+ x4 + x6 + · · ·

=

 

 

 

1

 

,

y2 = x + x3 + x5 + · · ·

=

 

 

 

x

 

.

1

x2

 

1

x2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14.2. y1 = 1

 

 

x2

 

 

x3

 

 

 

 

x4 − · · ·, y2 = x + x2 +

 

x3 +

 

 

 

x4 + · · · . 14.3. y1 =

2

2

24

6

4

 

 

 

1

x3 +

1

 

x5 + · · · , y2 = x−

1

 

x4

 

 

1

 

x6 + · · ·

. 14.4. y1 = 1

 

x2

 

 

 

x4

 

 

 

 

 

 

 

sin x

,

1

 

 

 

 

 

 

 

+

 

 

 

 

 

+

 

−· · ·

=

 

 

 

6

120

12

180

3!

5!

 

 

x

y2

= x1

11 2! + 4! − · · ·

 

=

 

 

x

 

 

. 14.5. y1

= x1

1 + x + 2!

+ 3! + · · · =

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

x4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

x

, y2 = x 2

1 + 5 x +

 

 

 

2

7

+ 5 7

9 + · · · . 14.6. y1 = 1+

1

 

+ 1

 

 

3 + 1

 

 

3

 

5 +· · ·,

 

ex

 

 

 

 

 

 

 

2

 

 

 

 

(2x)2

 

 

 

 

 

 

(2x)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

x4

 

 

x6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

5

 

 

 

 

 

 

 

 

 

·7

 

 

 

 

 

 

 

 

 

·

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

·

 

 

·

 

 

 

 

 

 

y

 

= x +

x

+

x

 

 

+

 

 

 

 

x

 

 

+

· · ·

. 14.7.

y

 

= 1

 

 

 

x2

+

 

3 · 5

x4

3 · 5 · 7

x6 +

· · ·

,

 

 

 

 

 

 

 

2

· 4 · 6

 

 

 

2 · 4

2 · 4 · 6

 

 

2

 

 

 

 

2

2 · 4

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

y

 

= x

 

4

x3 +

4

·

6

x5

4

·

6

· 8

x7

+

· · · . 14.8.

y

 

= 1 +

1

x2

+

 

 

1

x4 +

5

x6 +

· · ·

,

2

3

 

 

 

5

3

 

5

 

1

2

12

72

 

 

 

 

 

 

3

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

y2

= x + 6 x3 + 24 x4 + · · · . 14.9. y1

= x3

 

x2 + 5 6

+ 5

 

 

 

 

11 12 + · · ·

, y2 =

 

 

 

 

 

 

 

 

1

 

 

1

·

 

 

 

 

 

 

 

 

·

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

x4

 

 

 

 

 

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

6 7 12 13 + · · · .

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

 

·

 

·

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3 x2 + 6x 7 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

 

 

·

 

·

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

67