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

Обыкновенные дифференциальные уравнения в задачах. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
688 Кб
Скачать

ãäå 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]