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

ДУ РФ

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

6.14. y − 4y + 4y = cos 2x.

6.15. y + 9y = cos 3x.

6.16. y + y = x2 + 2x.

6.17. y − 6y + 9y = e3x.

6.18. y + y = (4x − 3) ex.

 

Ответы:

 

 

 

1

 

 

 

 

9

 

 

 

6.1 y = C1 e−x + C2 e2x +

 

(2x2

2x + 1). 6.2. y = C1 + C2 e4x +

 

(2x2 − x).6.3.

4

16

 

 

 

 

 

 

 

 

1

 

1

 

y = e−x(C1 cos x + C2 sin x) +

 

 

e3x. 6.4. y =

C1 ex + C2 e2x

 

x e−x. 6.5.

 

17

3

 

 

 

 

x2

 

 

 

1

 

 

 

 

 

y = (C1 + C2x) ex +

 

 

ex. 6.6.

y = C1 + C2 e2x +

 

(x2 − x) e2x. 6.7. y = C1 ex +

 

2

4

1

 

sin 3x). 6.8. y = C1 cos x + C2 sin x + 2x(sin x − cos x). 6.9.

C2 e6x +

 

(7 cos 3x

150

y = e−x(C1 cos x + C2 sin x) +

1

ex(cos x + sin x). 6.10. y = C1 cos x + C2 sin x +

1

ex.

8

 

 

 

 

 

 

 

 

5

 

 

 

 

 

1

 

2

 

6.11. y = C1 ex + C2 e2x +

e−x. 6.12. y = C1 e2x + C2 e3x +

(18x2 + 30x + 19). 6.13.

3

34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

7+

33 x

733 x

 

 

 

 

 

 

 

 

 

 

y = C1 e

2

 

+ C2 e

2

 

+

 

 

(3 sin x + 7 cos x). 6.14. y = (C1 + C2x) e2x

 

 

sin 2x.

 

 

 

 

 

 

58

8

6.15. y

=

C1 cos 3x + C2 sin 3x +

x

sin 3x. 6.16. y =

C1 + C2 e−x +

x3

. 6.17.

 

 

 

 

 

 

x2

 

 

 

6

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = (C1 + C2x) e3x +

 

e3x. 6.18. y = C1 cos x + C2 sin x +

 

 

(4x − 7) ex.

 

 

 

 

 

 

2

2

 

 

 

 

 

 

Занятие 7. Операторный метод решения линейных уравнений с постоянными коэффициентами

Обозначим через D оператор дифференцирования, D = d . Тогда dky = Dky, и уравнение dx dxk

a0y(n) + a1Y (n−1) + · · · + an−1y + any = f (x)

запишем в операторном виде

(a0Dn + a1Dn−1 + · · · + an−1D + an)y = f (x)

или

F (D)y = f (x),

где F (D) = a0Dn + a1Dn−1 + · · · + an−1D + an — операторный многочлен. Характеристическое уравнение однородного уравнения можно записать в виде

(роль λ выполняет также D)

a0Dn + a1Dn−1 + · · · + an−1D + an = 0.

31

Частное решение неоднородного уравнения ищется с помощью обратного опе-

ратора

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D)

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

 

 

 

f (x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D)

 

 

 

 

 

 

 

Для обратного оператора справедливы следующие свойства

 

 

 

 

1)

1

 

 

 

1

 

 

 

 

 

 

— постоянная,

 

 

 

 

 

 

 

 

 

 

cf (x) = c

 

 

f (x), c

 

 

 

 

 

 

 

 

F (D)

F (D)

 

 

 

 

 

 

 

 

2)

1

(f1(x) + f2(x)) =

 

1

 

 

f1(x) +

1

f2

(x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D)

F (D)

F (D)

 

 

 

 

3)

1

 

·

1

 

 

f (x) =

 

 

1

 

·

 

 

1

 

f (x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1(D)

F2(D)

F2(D)

F1(D)

 

(D) F2

(D) f (x)

,

 

4) F1(D) · F2(D) f (x) = F1

(D)

· F2

(D) f (x) = F1

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

1

 

 

 

5)

Φ(D)

f (x) = Φ(D)

1

f (x) =

 

1

 

 

 

Φ(D)f (x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6)

F (D)

 

 

 

 

 

F (D)

 

 

 

F (D)

 

 

 

 

 

 

 

 

D f (x) =

f (x) dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При отыскании частного решения помимо этих свойств используются формулы:

 

 

1

 

 

 

 

 

ekx

 

 

 

 

 

 

 

 

 

 

7)

 

 

ekx

=

 

 

 

, если F (k) = 0,

 

 

 

 

 

 

F (D)

F (k)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

sin ax

1

 

 

 

cos ax

, F (−a2) = 0,

8)

 

 

sin ax =

 

 

,

 

 

cos ax =

 

 

 

F (D2)

F (−a2)

F (D2)

F (−a2)

9)

 

1

1 =

 

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D)

 

an

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

10)

 

 

ekxv(x) = ekx

 

v(x).

 

 

 

 

 

 

F (D)

F (D + k)

 

 

 

 

 

Последняя формула носит название формулы смещения.

 

 

П р и м е р 1. Решить уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5y + 7y − 12y = e3x.

 

 

 

 

 

Р е ш е н и е. Запишем уравнение в операторном виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5D2 + 7D − 12)y = e3x.

 

 

 

Находим общее решение однородного уравнения:

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

x

125 x

5D

 

+ 7D − 12 = 0, D1 = 1, D2 =

 

, y0(x) = C1 e

 

+ C2 e .

 

5

 

32

При отыскании частного решения воспользуемся формулой 7):

 

 

 

 

 

 

1

 

e3x =

 

 

 

 

e3x

1

e3x.

 

 

 

 

 

yч =

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

5D2 + 7D − 12

5 · 32 + 7 · 3 12

54

 

 

 

 

 

Ответ: y = C1 e−x + C2 e125 x +

1

 

e3x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и м е р 2. Решить уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y − y = sin x.

 

 

 

 

 

 

 

 

 

 

 

Р е ш е н и е. Запишем уравнение в операторном виде

 

 

 

 

 

 

 

 

 

 

 

 

(D3 1)y = sin x.

 

 

 

 

 

 

 

 

 

 

 

Находим общее решение однородного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D3 1 = 0, (D − 1)(D2 + D + 1) = 0,

 

D1 = 1, D2,3 =

1

± i

3

,

 

 

 

 

2

2

 

 

 

x

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0(x) = C1 ex + e2 C2 cos

3

+ C3 sin

3

 

 

 

 

 

 

2

2

 

 

 

 

 

 

Ищем частное решение неоднородного уравнения.

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч =

 

 

 

sin x.

 

 

 

 

 

 

 

 

 

 

 

 

D3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как обратный оператор содержит нечетные степени, то воспользоваться одной из формул 8) нельзя. Поэтому “дополняем” знаменатель обратного оператора до разности квадратов, а затем воспользуемся свойством 5) и первой формулой 8).

yч =

 

1

 

sin x =

D3

+ 1

sin x = (D3 + 1)

 

1

sin x =

D3

1

D6

 

 

 

D6

 

 

1

 

 

1

1

 

 

 

 

 

 

 

 

1

 

 

=

 

(D3 + 1) sin x =

 

(cos x + sin x).

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: y = C1 ex + e2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

C2 cos

 

 

 

 

+

(cos x − sin x).

 

3

+ C3 sin

 

3

 

 

2

2

2

П р и м е р 3. Решить уравнение

y − 4y + 4y = e2x cos 3x.

Р е ш е н и е. Запишем уравнение в операторном виде

(D2 4D + 4)y = e2x cos 3x.

Находим общее решение однородного уравнения:

D2 4D + 4 = 0, (D − 2)2 = 0, D1,2 = 2, y0(x) = (C1 + C2x) e2x.

33

(x2 − x + 2).

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

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1 sin 3x

cos 3x

yч =

 

 

 

e2x cos 3x = e2x

 

 

cos 3x = e2x

 

 

 

 

 

 

 

= e2x

 

.

(D − 2)2

D2

D

 

 

3

 

9

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: y = (C1 + C2x) e2x

 

e2x cos 3x.

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

П р и м е р 4. Решить уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y + 4y = sin 2x.

 

 

 

 

 

 

 

 

 

Р е ш е н и е. Запишем уравнение в операторном виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(D2 + 4)y = sin 2x.

 

 

 

 

 

 

 

 

 

Находим общее решение однородного уравнения:

 

 

 

 

 

 

 

 

 

 

 

D2 + 4 = 0,

D1,2 = ±2i,

 

y0(x) = C1 cos 2x + C2 sin 2x.

Ищем частное решение неоднородного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч =

 

1

 

sin 2x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D2 + 4

 

 

 

 

 

 

 

 

 

8) неприменима, так как F (

a2) = 22 + 4 = 0

. Воспользуемся формулой

Формула

i2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эйлера: e

 

= cos 2x + i sin 2x. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

sin 2x = Im

 

 

1

 

 

ei2x = Im ei2x

 

 

 

1

 

 

 

1 =

 

 

yч =

 

 

 

 

 

 

 

 

 

 

 

D2 + 4

D2 + 4

(D + 2i)2 + 4

 

 

 

1

 

 

1

1 = Im

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

= Im ei2x

 

·

 

(cos 2x + i sin 2x)

 

 

=

 

cos 2x.

 

D

D + 4i

4i

4

Ответ: y = C1 cos 2x + C2 sin 2x − x4 cos 2x. П р и м е р 5. Решить уравнение

y − 4y + y = x2 − x + 2.

Р е ш е н и е. Запишем уравнение в операторном виде

(D2 4D + 1)y = x2 − x + 2.

Находим общее решение однородного уравнения:

D2 4D + 1 = 0, D1,2 = 2 ± 3, y0(x) = e2x(C1 ch 3x + C2 sh 3x).

Ищем частное решение неоднородного уравнения:

1

yч = D2 4D + 1

34

Для того, чтобы подействовать обратным оператором на многочлен, нужно поделить числитель (в данном случае 1) на знаменатель, записанный по возрастающем степеням (в данном случае на 1 4D + D2). Деление производится до тех пор, пока не получим слагаемое со степенью D, равной степени многочлена. Итак, имеем

 

 

 

 

1

 

2

 

1 4D + D2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 4D + D2

 

1 + 4D + 15D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4D − D

 

D3

 

 

 

 

 

 

 

 

 

D

D2

 

 

 

 

 

4

 

16 2

 

+ 4 3

 

 

 

 

 

 

 

 

 

 

15D

4D

 

Используя результат от деления, получаем

 

yч =

 

 

1

 

 

 

(x2 − x + 2) = (1 + 4D + 15D2)(x2 − x + 2) =

 

 

 

D2

4D + 1

 

 

2

 

 

 

 

 

 

 

2

 

 

 

= x

− x + 2 + 8x − 4 + 15 · 2 = x + 7x + 28.

Ответ: y = e2x(C1 ch 3x + C2 sh 3x) + x2 + 7x + 28.

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

7.1. y − 2y − 3y = e4x.

7.3. y + y − 2y = 3x ex.

7.5. y + y = cos 4x.

7.7. y − 5y + 4y = 4x2 e2x.

7.9. y + y = x sin x.

7.2. y − y = 2 ex − x2.

7.4. y − 3y + 2y = sin x.

7.6. y + y = 4 cos x.

7.8.y − 9y = e3x cos x.

7.10.y − 3y + 2y = 2x.

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

7.11. y + y = 4x ex.

7.13. y + 2y − 3y = x2 ex.

7.15. y − 2y + y = 6x ex.

7.17. y − 5y = 3x2 + sin 5x.

7.12. y + 3y − 4y = e4x + x e−x.

7.14. y − 4y + 8y = e2x + sin 2x.

7.16. y + 4y + 4y = x e2x.

7.18. y + 2y = x2x.

35

Ответы:

7.1. y = C1 e−x + C2 e3x +

 

1

e4x. 7.2. y = C1 ex + C2 e−x + x ex + x2 + 2. 7.3. y =

5

C1 ex + C2 e2x + + ex

2

 

3

. 7.4. y = C1 ex + C2 e2x +

10 (sin x + 3 cos x). 7.5.

 

 

x2

 

 

 

x

 

1

 

y = C1 cos x+C2 sin x−

1

cos 4x. 7.6. y = C1 cos x+C2 sin x+2x sin x. 7.7. y = C1 ex +

 

 

 

15

C2 e4x + (2x2 + 2x − 3) e2x. 7.8. y = C1 e3x + C2 e3x +

 

 

1

 

(6 sin x

cos x) e3x. 7.9 y =

 

37

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

 

 

 

 

 

 

 

C1 cos x+C2 sin x+

 

 

sin x−

 

 

cos x. 7.10. y = C1 ex +C2 e2x +

 

 

 

 

 

 

 

 

. 7.11.

4

 

4

ln2 2 3 ln 2 + 2

y = C

cos x+ C

sin x+ 2(x

 

 

1) x

7.12. y = C1

 

x + C2

 

 

4x

x

 

4x

x

+

1

e

−x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

3 e .

2

 

x

e

 

 

 

 

 

e

 

 

 

5 e

 

 

6 36

 

.

7.13. y = C1 ex + C2 e3x +

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

ex. 7.14. y = e2x(C1 cos 2x + C2 sin 2x) +

12

16

32

 

1

 

e2x

+

1

(2 cos 2x+sin 2x). 7.15. y = (C1 +C2x) ex +x3 ex. 7.16. y = (C1 +C2x) e2x +

4

20

 

x

1

 

e2x. 7.17.

y = C1 + C2 e5x

1

x3

 

 

3

x2

 

 

 

6

 

 

 

 

 

1

 

(cos 5x − sin 5x).

 

 

 

 

 

 

 

 

x +

 

 

16

 

32

5

25

125

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

2 ln 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.18. y = C1 cos 2x + C2 sin 2x + 2x

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 + ln2 2

(2 + ln2 2)2

 

 

 

 

 

 

 

 

 

Занятие 8. Операторный метод решения линейных уравнений степени выше второй

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

8.1.

y − 3y + 2y = e−x(4x2 + 4x − 10).

8.2.

y − 3y + 3y − y = ex.

8.3. y(V )

8.4.

y(IV ) + 8y + 16y = cos x.

8.5. y(IV

8.6.y(IV ) + 5y + 4y = 2 sin x cos 4x.

8.7.y + y = sin x + x cos x.

8.8.y − y − y + y = 3 ex + 5x cos x.

+ y = x2 1.

) + 2y + y = cos x.

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

8.9.y − 4y + 5y − 2y = 2x + 3.

8.10.y(V ) + y = x3.

8.11.y(IV ) + 18y + 81y = sin 2x sin 4x.

36

8.12.y(IV ) + 5y + 4y = 3 sin x.

8.13.y(IV ) − y = x ex + cos x.

8.14.y − 2y + 4y − 8y = e2x sin 2x.

8.15.y(IV ) + y = 7x − 3 cos x.

8.16.y − 4y + 3y = x2 + x e2x.

Ответы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3

8.1. y = (C1 + C2x) ex + C3 e2x + (x2 + x−1) e−x. 8.2. y = (C1 + C2x+ C3x2) ex +

 

ex.

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x5

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

8.3. y = C1 + C2x + C3x2 + C4 cos x + C5 sin x +

 

 

 

 

 

. 8.4. y = (C1 + C2x) cos 2x +

60

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

(C3 + C4x) sin 2x +

 

 

 

cos x.

8.5. y = (C1 + C2x) cos x + (C3 + C4x) sin x −

 

 

 

 

 

cos x.

9

8

 

8.6. y = C1 cos x + C2 sin x + C3 cos 2x + C4 sin 2x +

1

sin 5x −

1

sin 3x. 8.7. y =

 

504

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

C1 + C2 cos x + C3 sin x +

 

 

sin x −

 

cos x. 8.8. y = (C1 + C2x) ex

+ C3 e−x +

 

 

x2 ex +

4

 

4

4

5

[(x − 1) cos x − (x + 2) sin x]. 8.9. y = (C1

+ C2x) ex + C3 e2x − x − 4. 8.10. y =

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x6

 

x4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 + C2x + C3x2 + C4 cos x + C5 sin x +

 

 

 

 

. 8.11 y = (C1 + C2x) cos 3x + (C3 +

120

4

 

 

 

 

 

cos 2x

 

 

cos 6x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

C4x) sin 3x+

 

 

 

 

 

. 8.12. y = C1 cos x+C2 sin x+C3 cos 2x+C4 sin 2x−

 

cos x.

50

1458

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

8.13. y = C1 ex +C2 e−x +C3 cos x+C4 sin x+

 

 

(x2 3x) ex

 

sin x. 8.14. y = C1 e2x +

8

 

4

 

 

 

 

 

 

 

1

e2x(sin 2x−2 cos 2x). 8.15. y = C1 +C2x+C3 cos x+C4 sin x+

C2 cos 2x+C3 sin 2x+

 

40

7

 

3x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

x3 +

 

sin x. 8.16. y = C1

+ C2 ex + C3 e3x +

 

(3x3 + 12x2 + 26x) +

 

(1 2x) e2x.

 

6

2

27

4

Занятие 9. Метод разложения оператора на дроби

Этот метод применяется к нахождению частного решения уравнения

a0y(n) + a1y(n−1) + · · · + an−1y + any = f (x)

с постоянными коэффициентами, когда f (x) — произвольная функция и формулы, приведенные ранее, не применимы. Частное решение ищем по формуле

1

yч = F (D) f (x).

Разложим операторный многочлен F (D) на множители

 

n

F (D) = a0Dn + a1Dn−1

k

+ · · · + an−1D + an = a0 (D − λk ),

 

=1

37

где Dk = λk — корни уравнения

a0Dn + a1Dn−1 + · · · + an−1D + an = 0

— действительные или комплексные. Используя метод неопределенных коэффициентов разложим обратный оператор на простые дроби

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

n

 

Ak

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D) =

a0(D

λ1)(D

λ2)

· · ·

(D

λn)

D

λk .

 

 

 

 

 

 

 

 

=

k=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

n

 

 

Ak

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч =

 

 

f (x) =

 

 

 

f (x) =

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

F (D)

 

 

k=1 D

 

λk

 

 

 

 

 

 

 

=

 

 

Ak

 

eλkxf (x) e−λkx =

 

Ak eλkx

1

f (x) e−λkx.

 

 

 

 

 

D

λk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, нахождение частного решения сводится к вычислению интегралов

 

 

 

 

 

 

 

D f (x) e−λkx =

f (x) e−λkx dx.

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и м е р ы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) Решить уравнение

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y − 3y + 2y =

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2x + 1

 

 

 

 

 

 

Р е ш е н и е. Запишем уравнение в операторном виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(D2 3D + 2)y =

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2x + 1

 

 

 

 

 

 

Находим общее решение однородного уравнения:

 

 

 

 

 

 

 

D2 3D + 2 = 0, D1 = 1, D2 = 2, yo(x) = C1 ex + C2 e2x.

Ищем частное решение неоднородного уравнения:

 

 

 

 

 

 

 

yч =

1

1

 

 

=

1

 

 

 

 

 

 

 

1

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D2 3D + 2 e2x

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(D − 1)(D − 2) e2x + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

D − 2

D − 1

 

e2x + 1

=

 

 

 

 

 

 

1

e2x

 

 

e2x

1

 

 

ex

 

e−x

 

= e2x

1 e2x

ex

1 e−x

=

 

 

 

 

 

 

 

 

 

 

 

.

D − 2

 

e2x + 1

D − 1

e2x + 1

D

e2x + 1

D

e2x + 1

38

Вычисляем интегралы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2x

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

e2x + 1

 

e2x( e2x + 1)

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

2x

 

 

 

 

 

 

 

 

 

 

 

e2x

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

dx =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx =

 

e2x

e2x + 1

 

 

2

 

 

e2x( e2x + 1)

e2x

1

 

 

 

 

 

 

 

d e2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2x

1

 

 

 

 

 

1

 

 

1

 

d e2x =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

e2x( e2x + 1)

 

2

 

 

 

2

e2x

e2x + 1

 

 

 

 

e2x

 

 

 

1

ln e2x +

1

ln( e2x + 1) =

 

e2x

 

 

 

 

 

1

ln( e2x + 1).

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− x +

 

 

 

 

 

 

2

 

 

2

2

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

e−x

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

e2x + 1

 

 

ex( e2x + 1)

 

 

 

 

 

 

 

 

 

 

ex

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

=

 

 

 

 

 

dx =

 

 

 

 

 

 

 

 

d ex

=

 

 

 

 

arctg ex.

 

 

e2x( e2x + 1)

 

e2x

e2x + 1

ex

Итак,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч = e2x

 

 

2x

1

ln( e2x + 1) + ex

 

1

 

 

 

 

 

 

 

 

 

 

 

=

 

e

 

 

 

− x +

 

 

 

+ arctg ex

 

2

 

 

 

2

ex

 

 

 

 

 

 

 

 

= −x e2x +

 

e2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln( e2x + 1) + ex arctg ex +

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C1 ex + C2 e2x − x e2x +

 

e2x

 

 

ln( e2x + 1) + ex arctg ex +

1

.

 

2

 

 

 

 

2

2) Решить уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y + y =

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р е ш е н и е. Ищем частное решение неоднородного уравнения

 

 

 

 

 

 

 

 

 

 

 

yч =

 

1

 

 

 

 

1

 

 

 

=

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D2 + 1 sin2 x

 

 

(D − i)(D + i) sin2 x

 

 

 

 

 

 

1

1

 

1

 

 

 

1

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

=

 

 

 

2i

D − i

D + i

sin2 x

 

= 2i

D − i eix sin2 x

D + i e−ix sin2 x eix

=

1

1

 

 

 

 

e−ix

1

 

 

 

1

 

 

= 2i eix D sin2 x

e−ix D sin2 x eix .

 

 

1

 

 

 

1 e−ix

1

 

1

 

 

 

39

Вычисляем интегралы

 

 

D sin2 x =

 

 

1 e−ix

=

 

 

 

 

 

sin2 x dx − i

 

cos x

 

 

D sin2 x =

 

 

1 eix

=

 

 

 

 

 

sin2 x dx + i

 

cos x

cos x − i sin x

dx =

 

 

 

 

 

 

sin2 x

 

 

,

sin x

= sin x − i ln tg

2

dx

1

 

 

 

x

 

 

cos

x + i

sin x

 

dx =

 

 

 

 

 

 

 

 

 

 

 

 

sin2 x

 

 

.

sin x

= sin x + i ln tg

2

dx

1

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

2

 

 

 

 

yч = 2i (cos x + i sin x) sin x − i ln tg

 

 

 

1

 

 

 

 

 

 

1

 

 

 

x

 

 

 

 

 

1

 

+ i ln

tg

x

 

 

 

 

 

tg

x

.

= (cos x − i sin x) sin x

2

= 1 cos x ln

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что решение общее решение однородного уравнения имеет вид

y = C1 cos x + C2 sin x,

окончательно получаем

y= C1 cos x + C2 sin x − 1 cos x ln tg x2 .

За м е ч а н и е. В этом примере существует другой путь нахождения частного решения.

yч

= D2 + 1 sin2 x =

2i eix D sin2 x e−ix D sin2 x eix =

 

 

1

1

 

 

1

 

1 e−ix

 

 

1

 

 

1

 

 

 

 

 

 

== 2i eix

 

sin2 ξ dξ − e−ix

sin2 ξ =

 

 

 

 

 

 

1

 

 

 

 

 

e−iξ

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

sin2 ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ei(x−ξ) e−i(x−ξ)

 

1

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2 ξ

 

 

 

 

sin2 ξ

 

 

sin2

ξ

 

=

 

sin(x − ξ)

= sin x

 

cos ξ

 

 

cos x

sin ξ

=

 

 

 

 

 

 

 

 

 

= sin x

sin2 x dx − cos x

 

sin2 x dx = 1 cos x ln

 

tg 2 .

 

 

 

cos x

 

 

 

 

 

 

sin x

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ex

9.1. y + 4y + 3y = ex + 1 .

9.2. y − a2y = f (x), y(0) = 0, y (0) = 0.

1 9.3. y + 6y + 11y + 6y = ex + 1 .

40