Обыкновенные дифференциальные уравнения в задачах. Учебное пособие
.pdfãäå c1, c2 произвольные константы.
Построим теперь решение краевой задачи. Учитывая условие ограниченности решения при t > 0, очевидно, в формуле (4.5)
нужно положить c1 = 0. Следовательно, ограниченное при t > 0 общее решение уравнения имеет вид
y(t) = c2e2it + |
|
1 |
e t; |
|
i |
2 |
|||
|
|
ãäå c2 произвольная константа. Используя начальное условие y(0) = 1, однозначно находим
c2 = |
i |
3 |
: |
|
|
||
i |
2 |
Следовательно, краевая задача имеет единственное решение, которое можно записать в виде
y(t) = |
|
1 |
|
((i 3)e2it + e t) : |
|||
i |
2 |
||||||
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|||
Вариант 1. При каких вещественных параметрах , , |
âñå |
|||||||||
решения системы |
|
|
|
|
|
|
||||
|
|
dy |
0 |
9 |
|
cos2(3t) + |
|
|
||
|
|
|
|
= |
( 4 |
0 ) y + |
( |
sin(6t) |
) |
|
|
|
|
dt |
|
||||||
ограничены на числовой оси |
< t < 1? |
|
|
|||||||
Вариант 2. При каких вещественных параметрах , , |
âñå |
|||||||||
решения системы |
|
|
|
|
|
|
||||
|
dy |
|
2 |
|
sin2(2t) |
|
|
|||
|
|
= ( |
08 |
0 ) y + ( |
sin(3t) cos(t) + |
) |
|
|||
|
dt |
|
||||||||
91
ограничены на числовой оси |
< t < 1? |
|
|
||||||||
Вариант 3. При каких вещественных параметрах , , |
âñå |
||||||||||
решения системы |
|
|
|
|
|
|
|||||
|
|
dy |
0 |
|
sin(6t) |
|
|
||||
|
|
|
|
= ( |
12 |
03 |
) y + ( sin2(3t) + |
) |
|
||
|
|
dt |
|
||||||||
ограничены на числовой оси |
< t < 1? |
|
|
||||||||
Вариант 4. При каких вещественных параметрах , , |
âñå |
||||||||||
решения системы |
|
|
|
|
|
|
|||||
|
dy |
0 |
4 |
sin2(4t) + |
|
|
|||||
|
|
= ( |
16 |
0 ) y + ( sin(5t) cos(3t) |
) |
|
|||||
|
dt |
|
|||||||||
ограничены на числовой оси |
< t < 1? |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Решение (Вариант 1). Приведем два способа решения этой задачи.
Способ 1. Первый способ основан на общей формуле решения
системы |
|
|
dy |
|
|
|
|
|
|
|
|
= Ay + f(t) |
|
|
|||
|
|
|
dt |
|
|
|||
|
|
|
|
|
|
|
||
следующего вида: |
|
|
|
|
||||
|
|
y(t) = etAC + ∫0 t e(t |
s)Af(s) ds: |
|
|
|||
В нашем случае |
|
|
|
|
|
|
|
|
0 |
9 |
) ; C = ( |
C1 |
cos2 3t + |
) : |
|||
A = ( 4 |
0 |
C2 ) ; |
f(t) = ( |
sin 6t |
||||
92
Поскольку собственными значениями матрицы A
1 = 6i, 2 = 6i, то матричную экспоненту etA
âèòü â âèäå
etA = 1(t)I + 2(t)(A 1I);
являются числа можно предста-
ãäå
|
|
|
t |
|
|
|
|
|
|
|
1(t) = e6it; |
2(t) = ∫0 |
e 2(t s) 1(s) ds = |
sin 6t |
: |
||||||
|
||||||||||
6 |
||||||||||
Тогда после простых вычислений получаем |
|
|
|
|
||||||
|
0 |
2 cos 6t |
3 |
sin 6t |
1 |
|
|
|
||
etA = |
|
|
: |
|
|
|||||
|
2 |
|
|
|||||||
|
@ |
|
|
|
|
|
A |
|
|
|
|
B |
3 |
sin 6t |
|
|
cos 6t |
C |
|
|
|
Следовательно, решение системы можно переписать в виде |
||
y(t) = etA |
2C + ∫t e |
sAf(s) ds3 |
|
4 |
5 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
cos 6t |
3 |
|
sin 6t |
[ |
( |
C1 |
) |
|||||||||||||
|
|
|
|
|
= |
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
2 |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
@ |
|
2 |
|
|
|
|
|
|
|
|
|
A |
C2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
B |
3 |
sin 6t |
|
|
cos 6t |
C |
|
|
|
||||||||||||||
|
|
|
t |
0 |
|
|
2 cos 6s |
|
3 |
sin 6s |
1 |
( |
|
|
|
|
|
|
|
|
|
||||||||||||
+ |
|
|
|
|
|
|
|
2 |
+ |
|
|
|
+ |
|
cos 6s |
||||||||||||||||||
0 |
|
|
|
|
|
|
2 |
2 |
|||||||||||||||||||||||||
|
|
B |
3 |
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
∫ |
|
@ |
|
|
|
|
sin 6s |
|
|
|
cos 6s |
A |
|
|
|
|
|
|
sin 6s |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
2 cos 6t |
3 |
|
sin 6t |
1 |
|
|
|
|
|
C1 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
= |
2 |
|
|
3 |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
B |
|
|
sin 6t |
|
|
cos 6t |
C |
[( C2 |
) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
3 |
|
||||||||||||||||||||||
+ |
∫ |
t |
0 |
( + |
|
) cos 6s + |
|
cos2 6s + |
|
|
sin2 6s |
||||||||||||||||||||||
2 |
2 |
2 |
|||||||||||||||||||||||||||||||
|
|
|
@ |
2 |
|
( |
|
|
|
|
) |
|
( |
|
|
) |
|
|
|
|
|
|
|
||||||||||
|
0 |
|
|
B |
|
|
|
|
+ |
|
|
|
sin 6s + |
|
|
|
|
|
sin 6s cos 6s |
||||||||||||||
|
|
|
3 |
|
|
2 |
|
|
|
|
3 |
|
|
||||||||||||||||||||
)]
ds
1]
C
A ds :
93
Используя формулы
2 cos2 6s = 1+cos 12s; 2 sin2 6s = 1 cos 12s; 2 sin 6s cos 6s = sin 12s;
получаем
|
|
|
|
|
y(t) = 0 |
|
|
cos 6t |
3 |
|
sin 6t |
1 |
[ |
( |
|
C1 |
) |
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
@ |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
C2 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
B |
3 |
sin 6t |
|
|
|
cos 6t |
C |
|
|
|
|
|
|
|
|||||||||||||||||||||
∫ |
|
0 ( + |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
] |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
2 |
|
|
|
|
|
|
4( |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
+ |
t |
2 ) cos 6s + |
3 ) cos 12s + 4( + 3 ) |
1 ds : |
|
||||||||||||||||||||||||||||||||||||||||||
|
|
@ |
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
A |
|
|
|||
0 |
|
B |
|
|
|
|
|
|
+ |
|
sin 6s + |
|
|
|
|
|
|
|
|
|
|
|
sin 12s |
|
C |
|
|
||||||||||||||||||||
|
|
|
|
|
|
3 |
2 |
2 |
|
|
|
|
|
|
3 |
|
|
|
|
|
|||||||||||||||||||||||||||
Следовательно, общее решение системы имеет вид |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
y(t) = 0 |
|
2 cos 6t |
|
|
3 |
|
sin 6t |
1 |
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
3 |
sin 6t |
|
|
|
|
|
cos 6t |
C |
|
|
|
|
|
|
||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||||||||
0 |
C1 + |
|
|
( + |
|
|
) sin 6t |
|
|
( |
|
|
|
|
|
|
) sin 12t + |
|
|
( + 3 )t |
1 |
|
|||||||||||||||||||||||||
6 |
2 |
16 |
|
|
|
|
3 |
4 |
|
: |
|||||||||||||||||||||||||||||||||||||
B |
|
C2 + 9 |
( |
+ 2 |
) |
(cos 6t 1) + |
24 |
|
|
|
|
|
|
3 |
(1 |
|
|
cos 12t) |
C |
|
|||||||||||||||||||||||||||
@ |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
( |
|
|
|
|
) |
|
|
|
|
A |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Из этого выражения, очевидно, вытекает, что необходимым и достаточным условием ограниченности решений на всей числовой оси является равенство
+ 3 = 0:
Параметр может быть произвольным.
Способ 2. Второй способ решения задачи основан на сведении системы к двум дифференциальным уравнениям второго порядка
94
на компоненты решения y1(t) è y2(t). Перепишем исходную систему уравнений в виде
8 |
dy |
= 9y2 |
+ cos2(3t) |
dt1 |
|||
> |
|
|
|
> |
dt |
|
|
> dy2 |
|
|
|
> |
|
|
|
: |
|
= 4y1 + |
sin(6t): |
< |
|
||
+ ;
(4:6)
Все решения этой системы являются аналитическими функциями (почему?). Тогда, очевидно, компоненты решений будут также удовлетворять следующему уравнению
d2y2 |
= 4 |
dy1 |
+ 6 cos 6t; |
2 |
dt |
||
dt |
|
|
которое получается после дифференцирования второго уравнения системы (4.6). Подставляя теперь вместо производной dydt1 правую часть первого уравнения системы, а также учитывая, что
2 cos2 3t = 1 + cos 6t, получаем уравнение, которому должна удовлетворять вторая компонента решения системы
d2y2 |
+ 36y2 = 2 + 4 + 2( + 3 ) cos 6t: |
(4:7) |
2 |
||
dt |
|
|
Для нахождения общего решения вначале рассмотрим характеристическое уравнение
2 + 36 = 0:
Его корни 1 = 6i, 2 = 6i. Поэтому общее решение однородного
уравнения
d2y2
dt2 + 36y2 = 0
имеет вид
y2;îäí(t) = c1 cos 6t + c2 sin 6t:
95
Найдем теперь частное решение уравнения (4.7). Поскольку все параметры , , вещественные, то при нахождении частного решения можно рассмотреть вспомогательное уравнение
d2v |
+ 36v = 2 + 4 + 2( + 3 )e6it; |
(4:8) |
|
dt2 |
|||
|
|
найти его частное решение v÷àñò(t), а затем получить частное решение уравнения (4.7) в виде y2;÷àñò(t) = Re v÷àñò(t). Учитывая,
что корни характеристического уравнения имеют вид 1;2 = 6i, а в правой части уравнения (4.8) присутствует функция e6it, òî
частное решение следует искать в виде v÷àñò(t) = a + bte6it:
Для нахождения констант a, b нужно подставить эту функция в
уравнение (4.8). Если она является решением, то должно выполняться тождество
36a + 12ibe6it 2 + 4 + 2( + 3 )e6it:
Отсюда получаем
|
a = |
+ 2 |
; |
b = |
+ 3 |
; |
||||
|
|
6i |
||||||||
|
|
|
18 |
|
|
|
|
|
||
ò. å. |
|
+ 2 |
|
|
+ 3 |
|
|
|
||
v÷àñò |
(t) = |
+ |
t(sin 6t |
i cos 6t) |
||||||
|
6 |
|||||||||
|
|
|
18 |
|
|
|
|
|
||
частное решение вспомогательного уравнения (4.8). Следовательно,
y2;÷àñò(t) = |
+ 2 |
+ |
+ 3 |
t sin 6t |
|
18 |
6 |
||||
|
|
|
частное решение уравнения (4.7).
Из проведенных рассуждений вытекает формула общего решения уравнения (4.7)
y2(t) = c1 cos 6t + c2 sin 6t + |
+ 2 |
+ |
+ 3 |
t sin 6t: (4:9) |
|
18 |
6 |
||||
|
|
|
96
Аналогичным образом получается уравнение, которому должна удовлетворять первая компонента решения системы
d2y1 |
+ 36y1 = 3( + 3 ) sin 6t: |
(4:10) |
dt2 |
Точно так же можно показать, что функция
y1;÷àñò(t) = |
+ 3 |
t cos 6t |
|
4 |
|||
|
|
является частным решением этого уравнения. Тогда формула общего решения уравнения (4.10) имеет вид
y1(t) = b1 cos 6t + b2 sin 6t + |
+ 3 |
t cos 6t: |
(4:11) |
|
4 |
||||
|
|
|
Следовательно, если вектор-функция
y(t) = (y1(t); y2(t))T
является решением системы дифференциальных уравнений (4.6), то первая компонента y1(t) имеет вид (4.11), а вторая y2(t) вид (4.9), где константы b1, b2, c1, c2 связаны соотношениями
|
3 |
|
|
|
3 |
3 |
|
|
||
b1 = |
|
c2 |
; |
b2 = |
|
|
|
|
c1 |
: |
2 |
|
24 |
2 |
|||||||
Тогда для того, чтобы все решения системы (4.6) были ограниче- ны на числовой прямой < t < 1, необходимо и достаточно выполнение условия
+ 3 = 0:
2001
Вариант 1. Найти общее решение уравнения
p p
y′′ 4iy = e 2t cos( 2t):
97
Вариант 2. Найти общее решение уравнения
p p
y′′ 16iy = e2 2t cos(2 2t):
Вариант 3. Найти общее решение уравнения
p |
|
t |
p |
|
|
|
2 |
||||||
|
||||||
y′′ + 4iy = e |
sin( 2t): |
|||||
Вариант 4. Найти общее решение уравнения
p p
y′′ + 16iy = e2 2t sin(2 2t):
Решение (Вариант 1). Рассмотрим характеристическое урав-
нение
2 4i = 0:
Его корни
p p p p
1 = 2 + i 2; 2 = 2 i 2:
Тогда общее решение однородного уравнения y′′ 4iy = 0
имеет вид p p p p
yîäí(t) = c1e( 2+i 2)t + c2e ( 2+i 2)t:
Найдем частное решение неоднородного уравнения. Перепишем правую часть следующим образом:
p |
|
|
|
|
1 |
p |
|
p |
|
|
1 |
|
p |
|
p |
|
|
e 2t cos(p |
2 |
t) = |
|
e( |
2+i 2)t + |
|
e( |
2 |
i 2)t: |
||||||||
2 |
2 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
98
Тогда, учитывая формулы корней характеристического уравнения, частное решение можно искать в виде
p p p p y÷àñò(t) = ate( 2+i 2)t + be( 2 i 2)t:
После подстановки в уравнение, очевидно, получаем
|
1 i |
|
|
i |
|
|||
a = |
8p |
|
|
; |
b = |
|
: |
|
16 |
||||||||
2 |
||||||||
Следовательно, общее решение уравнения имеет вид
p p p p y(t) = c1e( 2+i 2)t + c2e ( 2+i 2)t
|
1 i |
(p |
2+ip |
|
|
|
i |
(p |
|
ip |
|
|
|
|||
|
2)t |
|
2 |
2)t |
|
|||||||||||
+ |
8p |
|
te |
|
|
|
|
+ |
|
e |
|
|
|
|
|
: |
|
|
|
|
16 |
|
|
|
|
|
|||||||
2 |
|
|
|
|
|
|
|
|
|
|||||||
Отметим, что формулу общего решения можно было бы полу- чить, применяя метод вариации произвольной постоянной, а так-
же используя -функции:
∫t p p
y(t) = c1 1(t) + c2 2(t) + 2(t s)e 2s cos( 2s) ds;
0
ãäå |
|
|
|
|
|
|
e 1t |
e 2t |
||
1(t) = e 1t; |
|
|
|
|||||||
2 |
(t) = |
|
|
|
: |
|||||
|
2 |
|||||||||
|
|
|
|
|
1 |
|
||||
|
|
|
2002 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
Вариант 1. Решить задачу Коши |
|
|
|
|||||||
>8 y′′ + 2y′ |
j yj = t2 |
|
4t 2; |
|
t > 0; |
|||||
> |
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
y(0) = 0; |
1: |
|
|
|
|
|
|
|
|
|
> y′(0) = |
|
|
|
|
|
|
|
|
||
> |
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
: |
|
|
|
|
|
|
|
|
|
|
99
Сколько существует решений? Вариант 2. Решить задачу Коши
8 y′′ + 4y′ |
4jyj = 8t2 16t 4; t > 0; |
|
> |
|
|
> |
|
|
< |
|
|
> y(0) = 0; |
|
|
> |
|
|
> |
решений? |
|
Сколько существует: |
||
> y′(0) = |
|
1: |
Вариант 3. Решить задачу Коши |
||
8 y′′ 2y′ |
j yj = t2 + 4t 2; t > 0; |
|
> |
|
|
> |
|
|
< |
|
|
> y(0) = 0; |
|
|
> |
|
|
> |
решений? |
|
Сколько существует: |
||
> y′(0) = |
2: |
|
Вариант 4. Решить задачу Коши |
||
8 y′′ 4y′ |
4jyj = 8t2 + 16t 4; t > 0; |
|
> |
|
|
> |
|
|
< |
|
|
> y(0) = 0; |
|
|
> |
|
|
> |
решений? |
|
Сколько существует: |
||
> y′(0) = |
|
2: |
Решение (Вариант 1). Очевидно, рассматриваемая задача Коши имеет единственное непродолжаемое решение y(t), опреде-
ленное на интервале (t ; t+). Учитывая вид начальных условий,
100
