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

9869

.pdf
Скачиваний:
0
Добавлен:
25.11.2023
Размер:
3.38 Mб
Скачать

ГЛАВА 3

§1. Дифференциальные уравнения первого порядка

Вкурсе математики средней школы изучались алгебраические уравнения, где неизвестными были числа. Сейчас мы переходим к рассмотрению так называемых дифференциальных уравнений, при решении которых находят неизвестные функции, удовлетворяющие заданным соотношениям, включающим операцию дифференцирования.

Рассмотрим для начала задачу о законе изменения скорости свободного падающего тела. Пусть тело массы m падает с некоторой высоты. Учтем, что кроме силы тяжести, на него действует сила сопротивления воздуха. Запишем второй закон Ньютона

 

 

 

m

dv

mg kv ,

(1.1)

 

 

 

 

 

 

 

 

dt

 

 

предполагая, что сила сопротивления пропорциональна скорости

v

в

каждый момент

времени t

с коэффициентом пропорциональности

k .

Уравнение (1.1),

кроме неизвестной функции v v t , содержит еще и ее

 

dv

 

 

 

 

 

 

производную dt

 

 

 

 

v t . Это и есть дифференциальное уравнение.

 

 

Дадим общие определения. Дифференциальным уравнением

первого порядка называется уравнение

 

 

 

 

 

 

F x, y, y 0 ,

(1.2)

связывающее независимую переменную x и искомую функцию y

с ее

первой производной y . Если

y можно явно выразить через оставшиеся

переменные уравнения (1.2), то оно приобретает вид

 

 

 

 

 

 

y f x, y .

(1.3)

Решением дифференциального уравнения (1.2) называется функция y x , которая при подстановке в уравнение (1.2) обращает его в тождество.

Можно убедиться, в частности, что функция

140

 

 

 

 

 

 

 

v t C e

k

mg

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

m

(1.4)

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

при

любом значении постоянной C

удовлетворяет уравнению (1.1).

Действительно,

подставляя

функцию (1.4) и

ее производную

 

 

k

 

k

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

m Ce

 

 

 

 

 

 

 

v t

 

 

 

в (1.1), получим тождество. Это означает, что функция

вида (1.4) является решением уравнения (1.1).

 

 

Заметим,

что мы

нашли

бесконечно

много функций,

удовлетворяющих дифференциальному уравнению (1.1) – каждому

значению

постоянной

C соответствует свое

решение вида (1.4).

Множество

функций

y x,C , обращающих

уравнение (1.3) в

тождество, называют общим решением дифференциального уравнения

(1.3). Запись общего решения содержит произвольную постоянную C .

Заметим, что решение дифференциального уравнения может быть

записано и в неявном виде x, y,C 0 .

 

 

Допустим, что в рассматриваемой задаче известна скорость

тела в

начальный момент времени t 0. Обозначим её

v0 v 0 .

Чтобы

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

из найденного множества решений (1.4) только одно - то, которое

соответствует начальному

условию

v0 v 0 . При

t 0

и v v0 из

множества решений (1.4)

получим

v

 

C

mg

,

откуда

C v

 

 

mg

.

0

k

0

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

Подставляя найденное значение постоянной в (1.4), получим закон

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

v 0 v0 :

 

k

t

 

mg

 

k

t

 

 

 

 

 

v t v0 e m

 

 

1

e m .

(1.5)

k

 

 

 

 

 

 

 

 

 

141

Согласно последнему равенству, скорость v падающего тела при

t будет стремиться к величине

mg

. Отсюда,

в частности,

можно

k

 

 

 

 

найти нужный коэффициент сопротивления k

(парашют),

чтобы

обеспечить приземление с допустимой скоростью. Функция (1.5)

представляет собой так называемое частное решение уравнения (1.1),

соответствующее начальному условию v0 v 0 .

Частным решением уравнения (1.3) называется одна функция,

удовлетворяющая самому уравнению и начальному условию. Задачу нахождения частного решения дифференциального уравнения (1.3),

удовлетворяющего данному начальному

условию y(x0 ) y0 ,

называют

задачей Коши. Если правая часть

f x, y

уравнения

(1.3) непрерывна в

некоторой области, содержащей

начальную точку

x0 , y0 ,

и имеет

f

непрерывную в этой области частную производную y , то задача Коши

имеет единственное решение. При этих условиях частное решение получается из общего решения при конкретном значении произвольной постоянной C .

Процесс отыскания решения дифференциального уравнения называется его интегрированием, а график решения – интегральной кривой. Рассмотрим геометрическую интерпретацию решений уравнения (1.3) на конкретном примере. Пусть требуется найти частное

решение дифференциального уравнения

 

y 2x ,

(1.6)

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

 

y 2 4 .

(1.7)

Непосредственной подстановкой убеждаемся, что функция вида

 

y x2 C

(1.8)

142

обращает уравнение (1.6) в тождество. Она содержит произвольную постоянную C и является общим решением уравнения (1.6). Построив в плоскости xOy графики этих функций при различных значениях C . мы получим семейство парабол (см. рис.1).

Чтобы выделить из этого семейства интегральных кривых конкретную параболу, соответствующую условию (1.7), рассмотрим точку с координатами 2;4 . Через нее проходит парабола семейства (1.8), для

которой C 0 . Соответствующее решение

y x2 является искомым

частным решением.

 

 

 

 

 

y 2x

 

4

 

 

1

 

 

0

1

2

x

 

 

 

-2 -3

Рис. 1

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

Если правая часть f x, y дифференциального уравнения (1.3) может быть записана в виде произведения функций двух функций Q x и P y ,

зависящих от переменных x и y соответственно, то есть y Q x P y ,

то уравнение называют дифференциальным уравнением с

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

143

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

dx

, перепишем последнее уравнение в виде

 

 

dy

Q x P y

 

или dy Q x P y dx .

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

Умножая обе части последнего уравнения на

1

 

P y 0 , получим

 

 

P y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

Q x dx ,

 

 

(1.9)

 

 

 

 

 

 

 

P y

 

 

в котором каждая из переменных x и y находится в той части уравнения,

где ее дифференциал. Считая y известной функцией от x , равенство (1.9)

можно рассматривать как равенство двух дифференциалов и интегрировать обе части уравнения (1.9). Полученные при этом функции

 

 

 

 

dy

 

 

 

 

x

 

Q x dx

 

 

 

 

 

P( y)

P y

и

Q

будут

отличаться

 

постоянным

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

слагаемым:

P( y) Q x C

. Мы записали

соотношение,

связывающее

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение

y , независимую переменную x и произвольную постоянную C ,

это

 

соотношение

 

и

 

представляет

собой

общее

решение

дифференциального уравнения (1.3).

 

 

 

 

 

 

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

дифференциальной форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P x Q

y dx P

x Q y dy 0 ,

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

2

2

 

 

 

 

решается аналогично.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим для примера дифференциальное уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

x .

 

 

 

 

(1.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функцию

f

x, y

y

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

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

произведения f x, y Q x P y

1

 

y и переписать уравнение (1.10):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

144

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

y

 

 

или

 

dy

y

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

Умножая

 

обе

части последнего

 

уравнения

 

на функцию

 

 

1

 

y 0 ,

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

dx

 

 

 

 

 

dy

 

 

dx

 

 

 

 

 

 

 

 

 

ln

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

. Интегрируя

 

 

 

 

 

, находим

ln

y

x

c

,

 

y

x

y

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или ln

y

ln

Cx

, откуда

y C x – общее решение уравнения (1.10), где

C – произвольная постоянная.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим далее задачу Коши: найдем решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2 dy y2

dx 0 ,

 

 

 

 

 

 

 

 

 

 

 

(1.11)

при условии, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 0 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.12)

Дифференциальное уравнение (1.11) с разделяющимися

 

 

 

 

 

 

 

 

переменными запишем в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2 dy y2 dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножая

 

обе

части последнего

 

уравнения

 

на

 

 

 

1

 

 

 

 

 

y 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разделим переменные:

 

 

 

 

 

 

 

 

dy

 

 

dx

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

Интегрируя

 

 

 

dy

 

dx

,

 

находим

 

 

1

arctg x C , или

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

y

1 x

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

1y arctg x C , где C – произвольная постоянная.

Итак, общее решение уравнения (1.11) имеет вид

1y arctg x C .

Учет начального условия (1.12) дает

1

arctg 0 C ,

откуда C 1.

 

 

 

1

 

 

 

 

Следовательно, решение задачи Коши записывается в виде

1

arctg x 1

y

 

 

 

 

 

 

145

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

или y

 

.

 

 

 

 

 

 

 

 

 

 

 

1 arctgx

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим далее линейные дифференциальные уравнения

первого порядка, которые, по определению, имеют вид

 

 

 

 

 

 

y p x y q x .

 

 

 

 

(1.13)

Решение уравнения (1.13) будем искать в виде произведения

 

 

 

 

 

 

 

 

 

y u v

 

 

 

 

(1.14)

двух неизвестных функций

u u x и v v x , тогда

 

 

 

 

 

 

 

 

 

 

 

 

v .

 

 

 

(1.15)

 

 

 

y u v u v u

 

 

 

Подставив в уравнение

(1.13)

 

вместо y и

y равенства

(1.14)

и

(1.15)

соответственно, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p x u v q x ,

 

 

 

 

 

 

 

u v uv

 

 

 

 

 

или

 

 

 

 

 

p x v q x .

 

 

 

(1.16)

u v u v

 

 

 

 

Рассмотрение вместо

одной неизвестной

функции

y y(x)

двух

функций u x и v x

дает

возможность ввести

для одной из

них, в

частности v x , дополнительное условие, которое упростит уравнение.

Оно состоит в требовании обращения выражения

v p x v

в нуль, то

есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v p x v 0 .

 

 

 

 

 

(1.17)

 

 

Уравнение (1.17) является дифференциальным уравнением с

разделяющимися

переменными

v

и

x . Его

запишем в виде

 

dv

p x v 0

или

dv

p x v . Умножая

обе

части

последнего

 

 

 

 

dx

 

dx

 

 

 

 

 

 

 

уравнения на

dx

, разделяем переменные:

 

dv

p x dx . Интегрируем

 

 

 

 

 

 

v

 

 

 

 

 

v

 

 

 

 

dvv p x dx и находим одно из решений уравнения (1.17), например,

при постоянной C 0 . Это решение обозначим v v0 x . Для второй

146

неизвестной

функции u x из

(1.16) получим

уравнение u v0 q x .

Снова разделяем переменные

du

q x dx

и,

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

v0 (x)

 

 

 

 

 

 

 

u

q x dx

 

C , где C – произвольная постоянная.

v (x)

0

 

 

 

 

 

 

 

Подставляя найденные u x

и v v0 x

в функцию (1.14), получаем

 

 

 

q x dx

решение уравнения (1.13) в виде

y

 

v0 (x)

 

 

C v0 (x) .

 

 

Найдем для примера общее решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2 y x e2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.18)

В нем по условию p x 2 , q x xe2 x . Подставив в уравнение

 

y u v и

 

 

 

 

 

 

 

v u v

 

2u v

xe

2 x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

u v u v , получим u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

u v u v 2v xe2 x .

 

 

 

 

 

 

 

 

 

(1.19)

В

качестве

функции

v

 

возьмем одно решение

v v0 x уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

2v 0 при значении C

0 . Перепишем его в виде

 

dx 2v ,

 

разделим

 

 

 

переменные

 

dv

2dx

и, интегрируя

dv

2 dx , находим

ln

 

v

 

2x C .

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При C 0 получим

 

 

v v

 

e2 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.20)

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим функцию (1.20) в (1.19), получим u e

 

 

xe

 

или

du

 

 

 

 

2 x

2 x

dx x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Снова разделяя переменные

du x dx и интегрируя du x dx ,

 

 

 

 

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

x2

C ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(1.21)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где C – произвольная постоянная.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя найденные функции (1.20) и (1.21) в

равенство

 

y u v ,

получим общее решение данного уравнения (1.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

147

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

2 x

 

y

 

C

e

 

.

 

 

 

2

 

 

 

 

§ 2. Линейные однородные дифференциальные уравнения

второго порядка

Начнём с задачи из механики. Рассмотрим прямолинейное движение

материальной точки массы

m по оси Ox . Отклонение точки от

положения равновесия будем

определять функцией x x(t) . Пусть

движение происходит под действием трёх сил: силы, притягивающей

точку к началу координат и имеющей

проекцию на ось Ox , равную

ax,a 0 ,

силы

сопротивления

среды,

которую

считаем

пропорциональной первой степени скорости bx,b 0 и возмущающей силы, направленной по оси Ox и равной F (t) в момент времени t .

Применяя второй закон Ньютона к движущейся массе, получим mx bx ax F (t) .

Разделим обе части уравнения на m и после введения новых

обозначений b / m 2h 0, a / m k2 и F (t) / m f (t) приведем его к виду

x 2hx k2 x f (t) .

(2.1)

Полученное уравнение относится к классу так называемых

линейных дифференциальных уравнений второго порядка, имеющих

вид

 

 

(x) y f ( x) .

 

(2.2)

y

a1 (x) y

a0

 

В них неизвестная функция

y(x)

 

 

 

и ее производные y (x),

 

y (x)

входят линейно. В качестве коэффициентов уравнения a0 ( x),a1 ( x)

и

f ( x)

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

(a, b) .

При этих условиях существует единственное решение уравнения

(2.2),

удовлетворяющее заданным начальным условиям

y(x0 ) y0 ,

 

,

x0 (a, b) .

y (x0 ) y1

 

 

 

148

Если правая часть уравнения (2.2) равна нулю:

 

a1

 

a0 ( x) y 0 ,

(2.3)

y

(x) y

то оно называется однородным,

в противном случае (если

f ( x) 0 ) –

неоднородным.

Уравнение вида (2.2) служит математической моделью разнообразных колебательных физических процессов, то есть процессов,

которые характеризуются определенной повторяемостью во времени.

Колебательные процессы широко распространены в природе и технике,

например, качания маятника часов, переменный электрический ток и т.д.

При колебательном движении маятника или груза на пружине изменяется координата центра масс, в случае переменного тока колеблются напряжение и сила тока. Физическая природа колебаний может быть разной, однако они описываются сходными характеристиками и уравнениями одинакового типа. Математические модели, сводящиеся к линейному дифференциальному уравнению второго порядка с постоянными коэффициентами, называют линейным осциллятором.

Методика решений рассматриваемых уравнений базируется на

следующем утверждении.

Если

y1 (x) и

y2 ( x) – два каких-либо

непропорциональных

друг

другу

решения

уравнения (2.3), т.е.

y2 ( x) y1 ( x) ,

то

общее

решение

yoo ( x)

однородного

дифференциального уравнения второго порядка имеет вид yoo (x) C1 y1 (x) C2 y2 (x) ,

где C1 ,С2 – произвольные постоянные. Следовательно, два любых непропорциональных решения однородного линейного

дифференциального уравнения второго порядка формируют его общее

решение. Однако нет общего метода отыскания функций y1 (x) и y2 ( x) . Их легко найти в случае, когда коэффициенты уравнения (2.3) являются числами. Обозначим их a0 и a1 :

y a1 y a0 y 0 .

(2.4)

149

 

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