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

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

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

Вариант 4. Для решения задачи Коши

8 x2y′′ + 5xy

2 ctg y = 2p

 

; x > 0;

3

>

 

 

 

 

 

<

 

 

 

 

 

> yjx=1 =

2

+ ;

 

 

>

 

 

 

 

 

:

= 0

 

 

 

 

> yjx=1

 

 

 

 

найти @@y

= 23 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение (Вариант 1). Перепишем рассматриваемую задачу. По аналогии с уравнением Эйлера сделаем замену

x = et; y(et; ) = ye(t; ):

В новых переменных уравнение принимает вид

(dt

I) dty + 3dty + 2 sin y = p2:

 

d

 

d

 

d

 

 

 

 

 

 

 

 

e

 

e

e

Тогда задачу Коши можно переписать следующим образом

 

 

 

+ 2 sin y = p

 

;

 

y

+ 2y

2

8 y′′

 

=

 

;

(4:43)

 

2

> ejt=0

e

 

e

>

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

: e

 

=

0

 

 

 

 

> ejt=0

 

 

 

 

 

y

 

 

 

 

 

 

 

Поскольку

@@ y(et; ) @@ ye(t; );

то в силу тождества

y(x; ) ye(ln x; );

x > 0;

141

нахождение производной @@y

= 4

сводится к вычислению @@y = 4 .

 

 

 

 

 

e

-

Для этого воспользуемся следствием

1 к теореме 3.3.4 о диффе

 

 

 

 

 

 

 

 

ренцируемости решения задачи Коши по компонентам начального вектора.

Перепишем задачу Коши (4.43) в эквивалентном виде (см. лем-

ìó 3.1.6)

8

 

 

(

 

 

 

 

 

 

) ;

 

dw

 

 

w2

 

>

 

=

 

 

2 sin w1 + p

 

 

 

dt

2w2

2

 

 

 

 

 

 

 

 

 

 

 

 

>

j

(

0

)

 

 

 

(4:44)

 

>

 

 

 

 

>

 

 

 

 

 

 

 

 

 

ãäå

>

 

 

 

 

 

 

 

 

 

 

>

 

 

=

2

 

;

 

 

 

 

< w t=0

 

 

 

 

 

:

 

 

 

 

 

 

(t; ) = y

 

 

w

(t; ) = y(t; );

w

(t; ):

 

1

 

 

 

e

 

 

2

et

 

Тогда вектор-функция

(

2

 

)

 

@

 

 

 

v(t) =

v1

(t)

=

0

@

w1(t; )

1

=

 

 

 

 

B

@

 

C =

4

 

v (t)

 

@

@

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

w2(t; )

C

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

0

 

@

y(t; )

1

 

@

B

 

 

 

t

 

C

@

 

 

B

ye

 

C =

@

 

 

 

A

 

 

 

(t; )

 

 

 

 

 

 

 

 

 

e

 

 

4

является решением задачи Коши для системы уравнений в вариациях вида (3.3.8)

>

 

dv

 

 

 

0

1

 

>

 

 

@

 

 

 

A

>

 

 

=

1

 

 

 

 

 

8 dt

0

 

 

1 v;

<

 

 

 

 

 

 

 

 

>

 

 

 

 

2 cos(w1(t; 4 ))

2

(4:45)

>

 

 

 

= (

 

) :

 

 

>

 

 

 

 

 

 

>

 

 

 

 

 

 

>

vjt=0

 

 

 

:

0

 

 

 

 

 

Найдем функцию w1(t; 4 ). Для этого нужно решить задачу Ко-

142

øè (4.44) ïðè =

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

4 , т. е. задачу Коши

 

 

 

 

 

dw

 

(

 

 

 

 

 

w2

 

 

) ;

 

>

 

 

 

 

 

 

 

 

dt =

 

 

 

 

2 sin w1 + p2

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

(

 

 

 

)

 

 

 

 

 

 

>

 

j

 

 

0

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

=

 

4

 

 

;

 

 

 

 

 

< w t=0

 

 

 

 

 

 

 

ãäå

:

 

 

 

 

 

 

 

 

 

w1(t; 4 )

1

 

 

 

 

 

 

 

w (t)

 

 

 

w(t) =

 

w1

(t)

)

= 0

 

:

 

 

 

 

(

2

 

 

 

@ w2(t;

4 )

A

 

Очевидно, отсюда мы приходим к задаче Коши на нахождение функции w1(t)

>8 w1′′ + 2w1

+ 2 sin w1 = p

 

;

2

>

 

 

 

<

 

 

 

>w1jt=0 = ;

>4 : w1jt=0 = 0:

Нетрудно догадаться, что w1(t) 4 . Подставляя эту функцию в (4.45), получаем задачу Коши для системы уравнений в вариаци-

ях в явном виде

>

 

dv

 

 

0

 

1

>

 

 

@

 

 

A

>

 

 

=

1

 

 

8 dt

0

 

 

1 v;

<

 

 

 

 

 

 

 

>

 

 

 

 

p

2

2

>

 

 

 

= (

 

 

) :

>

 

 

 

 

 

>

 

 

 

 

 

>

vjt=0

 

 

:

0

 

 

 

Тогда нахождение искомой функции

 

@

 

 

 

 

e

 

 

 

 

v1(t) =

@ y(t; ) = 4

143

сводится к решению задачи Коши

 

 

v1′′ + 2v1+ p

 

 

v1 = 0;

 

>

2

 

<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 v

= 1;

 

 

 

 

 

 

 

 

 

>

1jt=0

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

>

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

= 0:

 

 

 

 

 

 

 

 

 

 

> v1jt=0

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

2 + 2 + p

 

= 0:

 

 

 

 

 

 

2

Его корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1;2 = 1

 

 

 

= 1 i

 

 

 

:

 

 

p

 

p

 

1

 

1

2

2

Тогда общее решение дифференциального уравнения имеет вид

(√p ) (√p ) v1(t) = c1e t cos 2 1 t + c2e t sin 2 1 t :

Учитывая начальные условия, получаем

c1 = 1; 1 +

 

 

 

 

 

 

 

:

p

 

1

p

 

+ 1

2

c2 = 0; ò. å. c2 =

2

Отсюда

 

 

 

 

 

( (√p ) √p (√p )) v1(t) = e t cos 2 1 t + 2 + 1 sin 2 1 t :

Следовательно, возвращаясь к исходным переменным, имеем

 

 

 

 

 

 

@

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

y(x; ) = 4

= v1(ln x)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

(cos (√p2

1 ln x) + p2 + 1 sin (√p2 1 ln x)) :

x

144

2016

Вариант 1. Решить задачу Коши

>8 tyy′′

3t(y)2 + 5yy= 0; t > 0;

<

 

> yjt=1 = 4; : yjt=1 = 8:

Вариант 2. Решить задачу Коши

>8 tyy′′

2t(y)2 + 3yy= 0; t > 0;

<

 

> yjt=2 = 5; : yjt=2 = 5:

Вариант 3. Решить задачу Коши

>8 tyy′′ + t(y)2

3yy= 0; t > 0;

<

 

yjt=3 = 3;

 

> yjt=3 = 2:

 

:

 

Вариант 4. Решить задачу Коши

>8 tyy′′ + 2t(y)2

5yy= 0; t > 0;

<

 

yjt=4 = 2;

 

> yjt=4 = 1:

 

:

 

Решение (Вариант 1). Решить задачу Коши можно несколькими способами. Приведем некоторые из них.

145

Способ 1. Рассматриваемое уравнение имеет вид

F (t; y; y; y′′) = 0;

где функция

F (t; u1; u2; u3) = tu1u3 3tu22 + 5u1u2

является однородной относительно u1, u2, u3 степени 2, ò. å.

F (t; cu1; cu2; cu3) = c2F (t; u1; u2; u3):

Тогда, сделав замену

y(t) = y(t)z(t);

(4:46)

порядок уравнения можно понизить (см. Ÿ 3.5). Действительно, поскольку

y′′(t) = y(t)z(t) + y(t)z(t) = y(t)[z2(t) + z(t)];

то уравнение можно переписать в виде

 

Учитывая, что

y2

(

t[z2

+ z] 3tz2

)

не является реше-

 

 

+ 5z = 0:

 

тривиальная функция y(t) 0

нием рассматриваемой задачи, то по крайней мере в достаточно малой окрестности начальной точки t = 1 эта задача сводится к задаче Коши для уравнения Бернулли

tz 2tz2 + 5z = 0;

 

{ z

t=1 = 2:

(4:47)

j

 

 

Функция z(t) 0 не является ее решением, поэтому, делая замену u(t) = z 1(t), приходим к задаче для линейного дифференциаль-

ного уравнения

 

 

5

 

 

 

 

 

<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 u= t u

 

2;

(4:48)

 

> u

t=1 =

1

:

 

 

2

 

 

>

j

 

 

 

 

 

 

:

 

 

 

 

 

 

 

146

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

uîäí(t) = ct5;

а частным решением неоднородного уравнения является функция

t

u÷àñò(t) = 2:

Тогда общее решение неоднородного уравнения имеет вид

u(t) = ct5 + 2t :

Отсюда

u(t) = 2t

решение задачи Коши (4.48), а

z(t) = 2t

решение задачи Коши (4.47). Учитывая теперь (4.46), имеем

y(t) = 2t y(t):

В силу начального условия

yjt=1 = 4

решением исходной задачи является функция

y(t) = 4t2:

Способ 2. Умножая уравнение на t, имеем

y(t2y′′) 3(ty)2 + 5y(ty) = 0; t > 0:

147

p= y3ep 4:

Тогда, делая замену переменных, как в уравнении Эйлера

 

t = es;

y(es) = y(s);

(4:49)

уравнение можно переписать в виде

e

 

 

yy′′

3(y)2 + 4yy= 0:

 

Поскольку функция

ee

e

 

ee

 

F (u1; u2; u3) = u1u3

3(u2)2 + 4u1u2;

 

определяющая уравнение, не зависит от t, то, сделав замену

 

 

p(y) = y;

(4:50)

 

 

e

e

 

 

понизим порядок уравнения

ype p 3p2 + 4ype = 0:

Учитывая начальные условия в исходной задаче, получаем, что случаи p(ye) 0, ye(s) 0 исключены. Тогда приходим к уравне-

íèþ

Нетрудно найти общее решение уравнения

p(ye) = cye3 + 2y:e

С учетом замены (4.50) начальное условие имеет вид

pjye=4 = 8:

Тогда c = 0, ò. å. p(ye) = 2ye. Отсюда в силу замен (4.49), (4.50) приходим к задаче Коши

{ye= 2y;e yejs=0 = 4:

148

Очевидно, решение этой задачи

ye(s) = 4e2s;

и следовательно,

y(t) = 4t2:

Способ 3. В силу начальных условий функция y(t) const

не является решением задачи. Поэтому случаи y(t) 0, y(t) 0 исключены. Тогда исходное уравнение может быть переписано в

âèäå

y′′

 

y

5

 

 

 

 

 

 

3

 

+

 

= 0:

 

y

y

t

Нетрудно заметить, что его можно представить следующим обра-

çîì

 

 

 

 

 

j

 

 

j

 

 

j

j)

Понижая порядок уравнения( j j

 

 

 

 

 

d

ln y

3 ln

 

y

 

 

+ 5 ln t

= 0:

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln jyj

3 ln jyj + 5 ln jtj = const;

отсюда получаем

 

 

 

 

 

 

 

 

 

y3

 

 

 

 

 

 

y= c

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t5

 

 

 

 

 

 

 

 

 

 

 

 

 

В силу начальных условий c =

 

 

1

. Тогда мы приходим к задаче

 

8

 

 

 

 

 

 

 

 

 

 

 

 

Êîøè

 

 

 

8 y=

 

y3

 

 

 

 

 

 

 

;

 

 

 

 

 

8t5

 

Очевидно, общее решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< yjt=1 = 4:

 

 

 

 

 

уравнения определяется из соотноше-

íèÿ

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

=

 

 

 

 

+ ec:

 

 

 

 

 

 

 

 

 

 

 

y2(t)

 

16t4

 

149

Из начальных условий ec = 0 и решение задачи Коши имеет вид

y(t) = 4t2:

2017

Вариант 1. При каких начальных данных y0 решение задачи

Êîøè

t2y + ty y2 = 0;

 

 

 

 

t

 

1;

 

{ y

t=1

= y0

 

 

 

j

 

 

 

 

 

существует на всей полуоси ft 1g? Вычислить предел lim y(t).

t!1

Вариант 2. При каких начальных данных y0 решение задачи

Êîøè

t2y + 2ty y2 = 0;

 

 

 

 

t

 

1;

 

{ y

t=1

= y0

 

 

 

j

 

 

 

 

 

существует на всей полуоси ft 1g? Вычислить предел lim y(t).

t!1

Вариант 3. При каких начальных данных y0 решение задачи

Êîøè

t2y + 3ty y2 = 0;

 

 

 

 

t

 

1;

 

{ y

t=1

= y0

 

 

 

j

 

 

 

 

 

существует на всей полуоси ft 1g? Вычислить предел lim y(t).

t!1

Вариант 4. При каких начальных данных y0 решение задачи

Êîøè

t2y + 4ty y2 = 0;

 

 

 

 

t

 

1;

 

{ y

t=1

= y0

 

 

 

j

 

 

 

 

 

существует на всей полуоси ft 1g? Вычислить предел lim y(t).

t!1

150

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