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

Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 1. Учебное пособие

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

 

 

 

 

.

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

xs ( y ys ),

y ys (x xs )

 

 

 

 

(d)

Положим

 

 

x a2 ,

y

b2 ,

(x x ) U ,

y y

s

V

 

 

 

 

 

 

s

 

 

s

 

 

 

 

 

 

 

s

 

 

 

 

 

Тогда система (d) примет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U a2V ,

V

b2U

 

 

 

 

 

(e)

Продифференцируем по t, например, первое уравнение и подставим в него второе

 

 

 

 

 

 

 

..

a2b2U 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

Легко проверить, что общим решением этого уравнения является

 

 

 

 

 

 

 

U c1 cos(abt) c2 sin(abt) ,

 

 

 

 

(f)

где с1 и с2 - постоянные, определяемые из начальных условий (b). Из (e)

получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

1

.

b

c sin(abt) c cos(abt) .

 

 

 

 

 

 

 

 

U

 

 

 

(g)

 

 

 

 

 

a2

 

a

 

1

 

 

 

 

 

2

 

 

 

 

 

 

Исключая из (f) и (g) время, имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U 2

 

V 2

 

c2

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

 

1

a2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возвращаясь к исходным обозначениям, получим

 

 

 

 

 

 

 

 

 

 

 

(x x )2

 

( y y

 

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

s

 

1,

 

 

 

 

(h)

 

 

 

 

 

 

c2 x

 

c2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

где

с2 с2

.

Соотношение

 

(h)

это

 

 

 

 

 

 

c2 1

2

 

y

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнение эллипсов (в зависимости от с2, т.е.

 

 

 

 

 

 

от

начальных

условий)

с

центром

 

в точке

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(xs , ys ) ). Стрелки указывают направление,

 

 

 

 

 

 

соответствующее возрастанию времени. Из (d)

 

 

 

 

 

 

видно, что при

x xs , растет y

(т.к. >0)

и

 

 

 

x

s

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

наоборот.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

Видно, что изменение популяций носит циклический характер: через

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

1.2.1 Уравнения 1-ого порядка, интегрируемые в квадратурах

Рассмотрим теперь некоторые уравнения 1-ого порядка, которые

интегрируются в квадратурах.

1) Уравнения с разделяющимися переменными.

Уравнением с разделяющимися переменными называется уравнение вида:

 

dy

f (x)g( y)

(13)

 

dx

 

 

 

будем считать, что функции f (x) , g( y) непрерывны при a x b

,c y d , и

g( y) 0 , т.к. если, например, в какой-то точке y1 , функция g( y1 ) 0 , то из (13)

следует, что y y1 есть решение. Запишем уравнение (13) в виде:

dy f (x)dx g( y)

и проинтегрируем

 

dy

f (x)dx C .

(14)

g( y)

 

 

 

Это и есть общее решение уравнения (13).

Пример 7. Найдем решение следующего д.у. с разделяющимися переменными:

 

y y2 .

Перепишем уравнение в виде:

 

dy

 

dx .

 

y2

 

 

 

После интегрирования получим

 

1

x C .

y

 

 

 

 

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

y 1 (x C) ,

где С - произвольная постоянная. Кроме того, имеется еще решение y 0.

12

Найдем теперь решение следующей задачи (такую задачу называют задачей

Коши): y y2 , y(x0 ) y0

при y0 0 .

Определив значение C из начального

условия, получим:

 

 

 

 

 

 

 

y(x)

 

 

 

y0

 

 

 

 

.

 

1 y

0

(x x )

 

 

 

 

0

 

2)Уравнения, приводящиеся к уравнениям с разделяющимися переменными.

Втаких уравнениях переменные можно разделить, выполнив некоторую

подходящую замену.

а) Рассмотрим уравнение вида:

dydx f (ax by) ,

где a,b const .

Перейдем к новым переменным z ax by , в результате будем иметь

 

dz

a b

dy

, т.е.

dz

a bf (z)

 

 

 

dx

 

dx

 

 

dx

 

или

 

 

 

 

 

 

 

 

 

 

 

dz

dx ,

 

 

 

 

 

 

 

 

a bf (z)

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

dz

x a bf (x) C .

б) Рассмотрим теперь следующее уравнение

dy

 

y

 

 

f

 

 

(x, y) ,

 

 

dx

x

 

где (x, y) однородная функция нулевого порядка. Такое уравнение называют

однородным д.у. 1-ого порядка.

Функция (x, y) называется однородной функцией нулевого порядка, если

(tx,ty) t 0 (x, y) (x, y)

13

Выполним замену y x z , т.е. y xz , тогда

dy

x

dz

z, x

dz

z f (z),

dz

 

dx

.

 

 

 

f (z) z

 

dx

 

dx

 

dx

 

 

x

Таким образом, мы снова получили д.у. с разделенными переменными. После интегрирования получим:

 

dz

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

ln

x

ln C, x C exp

 

 

.

 

 

 

f (z) z

 

 

 

 

 

 

 

 

 

f (z) z

 

 

 

 

 

 

 

 

в) Уравнение вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M (x, y)dx N(x, y)dy 0

 

 

будет являться однородным, если

 

 

 

 

 

 

 

 

 

 

dy

 

M (x, y)

 

y

 

 

 

 

 

 

 

f

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

N (x, y)

 

x

 

 

т.е.

M (x, y) , N (x, y)

-

однородные

функции

одного

порядка,

( M (tx,ty) t n M (x, y),

N(tx,ty) t n N(x, y) ),

и

будет

уравнением

с

разделяющимися переменными, если

 

 

 

 

 

MM1(x)M2 ( y), N N1(x)N2 ( y) .

2.Линейные дифференциальные уравнения 1-го порядка.

Линейным дифференциальным уравнением 1-го порядка называется уравнение 1-го порядка, линейное относительно неизвестной функции и еѐ

производной.

 

 

 

dy

p(x) y f (x)

(1)

 

 

 

 

 

 

 

 

 

dx

 

 

где p(x), f (x) – непрерывны в области, где ищется решение уравнения (1).

 

Если

f (x) 0 , то

уравнение (1) называется линейным однородным,

в

противном

случае это

уравнение называется

линейным неоднородным.

В

14

линейном однородном уравнении переменные разделяются, и его решение не составляет труда:

 

dy

p(x) y 0, т.е.

dy

 

p(x)dx

 

 

dx

y

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

ln

 

y

 

p(x)dx ln C,

C 0,

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

y C exp{ p(x)dx}, C 0

(2)

При делении на y мы потеряли решение y=0, однако оно может быть включено в (2), если считать, что C может быть равно нулю.

Для интегрирования же неоднородных уравнений применяются метод вариации произвольной постоянной и метод подстановки, которые и будут рассмотрены ниже.

2.1 Метод вариации произвольной постоянной

Для интегрирования неоднородного уравнения (1) можно применять метод вариации произвольной постоянной.

В соответствии с этим методом вначале находится решение соответствующего однородного уравнения, которое задается формулой (2). Затем в полученном решении «варьируют постоянную», т.е. полагают C C(x) , и

подставляют это решение в исходное неоднородное уравнение (1). Тогда,

учитывая что:

dydx dCdx e p( x)dx C(x) p(x)e p(c)dx ,

будем иметь:

dCdx e p( x)dx C(x) p(x)e p( x)dx p(x)C(x)e p( x)dx f (x) ,

или

15

dCdx f (x)e p( x)dx .

Таким образом, мы получили линейное однородное д.у. относительно C(x),

интегрируя которое, находим

C(x) f (x)e p( x)dxdx C1 ,

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

y C(x)e

p( x)dx

C1e

p( x)dx

e

p( x)dx

f (x)e

p( x)dx

dx .

(3)

 

 

 

 

Итак из формулы (3) видно, что общее решение неоднородного линейного дифференциального уравнения 1-го порядка равно сумме общего решения соответствующего однородного уравнения

С1e p( x)dx

и частного решения неоднородного уравнения

e p( x)dx f (x)e p( x)dxdx ,

получающегося из (3) при С1 0 .

Пример 1. Методом вариации произвольной постоянной найдем решение неоднородного уравнения

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

y

x2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x

 

 

 

 

 

 

 

 

 

 

 

 

Вначале интегрируем соответствующее однородное уравнение:

 

dy

 

y

0, т.е.

dy

 

 

dx

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

 

y

 

ln

 

x

 

ln C; y Cx .

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x

 

 

y

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

dC

 

 

 

Затем полагаем

C C(x) ,

тогда y C(x) x,

и

 

x C(x) . Подставляя эти

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

dx

выражения в исходное уравнение, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dC

x C(x) C(x) x2 , т.е.

dC

x,

dC xdx ,

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

и после интегрирования будем иметь:

16

C(x) x2 C1. . 2

Таким образом, общее решение имеет вид:

y C1x x3 . 2

Пример 2. Методом вариации произвольной постоянной найдем решение неоднородного уравнения:

dydx yctgx 2x sin x .

Интегрируя однородное уравнение, находим

dyy cossin xx dx , т.е. dydx yctgx 0, ln y ln sin x ln C, y C sin x .

Варьируем постоянную:

y C(x)sin x, y C (x)sin x C(x)cos x .

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

C (x) sin x C(x) cos x C(x) cos x 2x sin x ,

получим:

 

C(x) x

2

C1

,

C (x) 2x,

 

Откуда

y x2 sin x C sin x .

 

 

 

1

 

 

 

2.2 Метод подстановки (Даламбера)

Для решения линейного неоднородного уравнения 1-ого порядка (1) можно применять также метод подстановки (Даламбера). В соответствии с этим методом в уравнении (1) применяют подстановку y(x) U (x)V (x) , тогда это уравнение примет вид:

V ( dUdx pU ) (U dVdx f ) 0.

Выберем функцию U так, чтобы первая скобка обратилась в нуль:

dUdx p(x)U 0, или dUU p(x)dx .

17

Обозначим через U0 (x) одно из частных решений этого уравнения, тогда для нахождения V имеем уравнение:

U0 (x) dVdx f (x) 0 .

Находим общее решение этого уравнения V V (x,C) . Теперь общее решение уравнения (1) будет иметь вид:

y(x) U0 (x)V (x, C) .

Пример 3. Методом подстановки найдем решение следующего

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2xy 3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Перепишем это уравнение в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

2x

 

 

3

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

y

 

 

 

y2

 

 

 

 

Оно линейно относительно x и

 

dx

 

. Положим x U ( y)V ( y) , тогда имеем:

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V (

dU

 

 

2

U ) (U

dV

 

3

) 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

y

 

 

 

 

 

 

 

 

dy

y2

 

 

Теперь

dU

 

2

U 0 , или

dU

 

 

 

2dy

,

т.е. U y2 C .

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

y

 

 

 

 

U

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

Положим U

0

y2

, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

dV

 

 

3

0, или

 

dV

3

 

dy , т.е. V ( y, C)

1

C ,

 

 

 

 

 

 

 

y4

 

y3

 

 

 

 

 

dy

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда окончательно имеем

x Cy 2 1y .

2.3 Уравнения, сводящиеся к линейным

18

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

1) Уравнение Бернулли

Уравнением Бернулли называется нелинейное относительно неизвестной

функции уравнение вида:

 

dy

p(x) y f (x) yn , (n 1) ,

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

y n

dy

p(x) y1 n

f (x) .

(4)

 

dx

 

 

 

 

 

 

 

 

 

 

Выполним замену переменной

y1 n z,

 

 

 

 

тогда

(1 n) y n

dy

 

dz

.

 

 

 

 

 

 

 

 

 

dx

 

dx

 

Подставив эти выражения в (4), получим линейное дифференциальное уравнение:

 

 

 

1 dz

 

p(x)z f (x) .

 

 

 

 

 

 

 

 

 

 

 

 

1 n dx

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Найдем решение следующего уравнения Бернулли:

 

 

 

 

dy

 

 

y

 

 

x2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

2x 2 y

Перепишем его в виде:

2 y

dy

 

y2

 

 

x2 , (n 1) .

 

 

 

 

 

dx

x

 

 

 

 

 

 

 

 

 

 

Выполним замену: y2 z.

Тогда 2 y

dy

 

 

dz

, и мы получим линейное уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

dx

из примера 1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

z

 

x2 ,

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

x

 

 

 

 

 

которое может быть решено как методом вариации произвольной постоянной, так и методом подстановки.

19

Отметим, что уравнение Бернулли можно решать и непосредственно методом подстановки.

2) Уравнение Риккати

Уравнением Риккати называется нелинейное относительно неизвестной

функции уравнение вида:

dydx p(x) y q(x) y2 f (x) .

В общем виде это уравнение не интегрируется в квадратурах, однако,

заменой переменной может быть преобразовано в уравнение Бернулли, если известно одно частное решение y1 (x) этого уравнения.

Полагая y y1

z , получим:

 

 

 

 

 

 

 

 

 

 

 

y1 z p(x)(y1 z) q(x)(y1 z

2

) f (x) ,

 

 

 

но так как

 

p(x) y1 q(x) y1 f (x) , то для z(x) получим уравнение Бернулли:

y1

 

 

 

 

2

 

 

z p(x) 2y1(x)q(x) z q(x)z2 0 .

Пример 5. Найдем решение следующего уравнения Риккати:

dydx y2 x22 .

Здесь легко подобрать частное решение: y1 1 x . Полагая y1 1 x z , получим:

 

 

 

 

 

1

 

 

 

 

 

 

1

 

1 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

z

x2

, т.е.

z

x2

 

(z x)

x2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

z

 

 

 

 

или:

 

 

 

 

 

z

z

 

2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее уравнение является уравнением Бернулли. Перепишем его в виде

z 2 1, z2 xz

Выполним замену:U 1z , тогда dUdx zz2 , а уравнение примет вид:

dUdx 2Ux 1; .

20

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