Дифференциальные уравнения первого порядка (96
..pdf
Выполним подстановку z = zu , тогда
xu′ + u = 1 |
+ |
5 |
u |
2 |
, |
|
|
2du |
|
= |
dx |
, |
2 |
|
du |
|
= |
dx |
. |
|||
|
|
|
|
|
|
|
|
|
|
|
1 2 |
3 2 |
|
|
||||||||
2 |
|
5u2 |
− 2u + 2 |
|
x |
5 |
|
x |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u − |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
5 |
|
|
|
|
После интегрирования получим
23 arctg 5u3−1 = ln x +C ,
откуда
u= 15 1 + 3 tg 32 ln x +C .
Сучетом y = 1 / xu, найдем окончательно общее решение:
y = |
|
|
5 |
|
|
|
|
. |
|
|
3 |
ln |
|
x |
|
||
|
x 1 |
+ 3 tg |
|
|
+C |
|||
|
|
|
2 |
|
|
|
|
|
Основные типы уравнений и способы их решения представлены в табл. 8.1.
31
Тип уравнения
С разделяющимися переменными
Однородные
уравнения
Линейные
уравнения
Таблица 8.1
|
|
|
Вид уравнения |
Способ решения |
|||||
y′ = f (x, y) |
P (x, y)dx +Q (x, y)dx = 0 |
||||||||
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||
y′ = f1 (x) f2 (y) |
p1 (x) p2 (y)dx + q1 (x)q2 (y)dy = 0 |
С учетом того, что y′ = |
dy |
||||||
|
|
|
|
|
, преобра- |
||||
|
|
|
|
dx |
|||||
|
|
|
|
зовать уравнение к виду M (y)dy = |
|||||
|
|
|
|
= N (x)dx и проинтегрировать обе |
|||||
|
|
|
|
части |
|||||
|
|
|
|
|
|||||
y′= f ( x, y) , где f (x, y) |
Можно привести к виду P(x, y)dx + |
Использовать подстановку u(x) = |
y |
, |
|||||
|
|
|
+Q (x, y)dy = 0 , где P(x, y) и Q (x, y) |
|
|||||
– однородная функция |
|
|
|
x |
|||||
нулевого измерения. |
должны быть однородными функциями |
что приведет уравнение к уравнению с |
|||||||
Можно привести к виду |
разделяющимися переменными. Здесь |
||||||||
y |
одинакового k-го измерения: |
y = xu( x) , а y′ = u′x +u |
|||||||
P(tx,ty) = tk P(x, y) Q(tx,ty) = tkQ(x, y) |
|||||||||
y′ = f1 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
||||
x |
|
|
|
|
|
|
|||
|
|
y′+ p(x) y |
= q(x), или x′+ p(y)x = q(y) |
Использовать: 1) метод вариации про- |
|||||
|
|
|
|
извольной постоянной; 2) подстановку |
|||||
|
|
|
|
Бернулли y = uv , причем y′ = u′v + |
|||||
|
|
|
|
+ uv′ , или соответственно x = uv и |
|||||
|
|
|
|
x′ = u′v +uv′ |
|||||
|
|
|
|
|
|
|
|
|
|
Тип уравнения
Уравнения
Бернулли
Уравнения в полных дифференциалах
Уравнения
Лагранжа
Уравнения
Клеро
|
|
|
|
|
|
|
Окончание таблицы 8.1 |
|
|
|
|
|
|
|
|
||
|
Вид уравнения |
|
|
|
|
Способ решения |
||
|
|
|
|
|
|
|
||
y′ = f (x, y) |
P (x, y)dx +Q (x, y)dx = 0 |
|||||||
|
||||||||
|
|
|
|
|
||||
y′+ p(x) y = q(x) yα |
или x′+ p( y)x = q( y)xα , где α≠0;1 |
Использовать: 1) введение новой пе- |
||||||
|
|
|
|
|
|
|
ременной по формуле z (x) = y1−α ; |
|
|
|
|
|
|
|
|
2) метод вариации произвольной по- |
|
|
|
|
|
|
|
|
стоянной; 3) подстановку Бернулли |
|
|
|
|
|
|
|
|
y = uv или соответственно x = uv |
|
|
|
|
|
|
|
|||
P(x, y)dx +Q (x, y)dy = 0 , |
где |
∂P |
= |
∂Q |
|
Восстановить функцию U (x, y) по ее |
||
|
∂x |
полному дифференциалу dU ( x, y) |
||||||
|
|
|
∂y |
|||||
|
|
|
||||||
a (y′) x +b(y′) y +c(y′) = 0 , или |
y = x f (y′)+ g (y′) |
Ввести параметр p = y′ и привести |
||||||
|
|
|
|
|
|
|
это уравнение к линейному |
|
|
|
|
|
|
|
|
||
|
y = x y′ + g (y′) |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
33
9. МЕТОД ИЗОКЛИН |
|
Уравнение, разрешенное относительно производной, |
|
y′ = f ( x, y) , |
(9.1) |
определяет в каждой точке плоскости ( x, y) , где существует функция f (x, y) , значение y′ , т. е. угловой коэффициент касательной к
интегральной кривой в этой точке.
Если в каждой точке области D задано значение некоторой величины, то говорят, что в D задано поле этой величины. Следовательно, дифференциальное уравнение (9.1) задает в каждой точке области определения функции f (x, y) направление касательной к интегральной кривой, т. е. определяет поле направлений. Совокупность отрезков этих касательных дает геометрическую картину поля направлений.
Задача интегрирования дифференциального уравнения (9.1) может быть сформулирована следующим образом: найти такие кривые, касательные к которым в каждой точке имели бы направление, совпадающее с направлением поля в этой точке.
Имея поле направлений, можно приближенно построить интегральные кривые. Построение поля направлений путем перебора координат большого количества точек весьма трудоемкая задача. Для ее облегчения удобно использовать изоклины.
Определение. Изоклиной называется геометрическое место точек одинакового угла наклона касательных к интегральным кривым.
Семейство изоклин дифференциального уравнения (9.1) определяется уравнением
f ( x, y) = k , |
(9.2) |
где k – параметр.
В каждой точке изоклины (9.2) касательные к интегральным кривым имеют один и тот же угловой коэффициент k , т. е. касательные параллельны между собой. Придавая k различные число-
34
вые значения, получаем сеть изоклин, с помощью которых можно приближенно построить поле направлений.
Пример 9.1. Методом изоклин построить семейство интеграль-
ных кривых уравнения |
|
|
dy = x . |
(9.3) |
|
dx |
|
|
Решение. Используя (9.2), найдем, что семейство изоклин урав- |
||
нения (9.3) есть множество прямых x = k |
на плоскости. |
|
Построим в плоскости (x, y) |
поле направлений: через точки |
|
изоклины, соответствующей числу |
k , |
проведем отрезки прямых |
таким образом, чтобы тангенс угла tg α между проводимым отрез-
ком и осью Ox равнялся k .
Построим приближенно кривые в плоскости ( x, y) , которые в
каждой своей точке касаются поля направлений (рис. 9.1). Следует отметить, что построенные кривые приблизительно иллюстрируют поведение интегральных кривых уравнения (9.3).
Рис. 9.1
Пример 9.2. Методом изоклин построить семейство интегральных кривых дифференциального уравнения
dy |
= ln(2x − y) . |
(9.4) |
dx |
|
|
|
|
35 |
Решение. Областью определения функции является множество точек, координаты которых удовлетворяют неравенству 2x − y > 0 ,
т. е. y < 2x . |
|
|
|
Функция f (x, y) = ln(2x − y) |
определена и непрерывна при |
||
y < 2x , а ее частная производная |
|
|
|
f ′ = − |
|
1 |
|
|
|
|
|
y |
2x − y |
||
|
|||
определена в этой же области. По теореме Коши через каждую точку области y < 2x проходит интегральная кривая и при том только
одна. |
|
|
|
|
|
k = ln(2x − y) , |
т. е. у = |
Изоклины |
задаются уравнениями |
||||||
= 2x − ek |
(рис. 9.2). Определим области возрастания и убывания |
||||||
частных |
решений, |
а также |
существование точек экстремума: |
||||
y′ = 0 ln(2x − y) = 0 , поэтому |
y = 2x −1 – геометрическое место |
||||||
точек возможных |
экстремумов; |
y′ > 0 y < 2x −1, а |
y′ < 0 |
||||
2x −1 < y < 2x , |
поэтому |
y < 2x −1 |
– область возрастания, а |
||||
2x −1 < y < 2x |
– область убывания интегральных кривых. На пря- |
||||||
мой y = 2x −1 находятся точки минимумов интегральных кривых. Исследуем вторую производную:
y′′ = (ln(2x − y))′ = 2 −ln(2x − y) . 2x − y
Интегральные кривые имеют выпуклость вниз при 2x −e2 < < y < 2x, и вверх – при y < 2x −e2 . Найдем асимптоты интеграль-
ных кривых. Проверим, существуют ли интегральные кривые, совпадающие с изоклинами (такие кривые являются асимптотами). Определим, при каком k график изоклины задается функцией
y = 2x −ek , являющейся решением дифференциального уравнения:
d (2x −ek ) = ln(2x −(2x −ek )) . dx
36
Рис. 9.2
Выполнив дифференцирование, найдем k = 2 , т. е. y = 2x −e2
является асимптотой и интегральной кривой. На границе области определения лежит прямая y = 2x , которая не является интеграль-
ной кривой. При y → 2x ln(2x − y) → −∞, т. е. y′ → −∞ и угол наклона интегральных кривых стремится к − π2 .
Замечание. Несмотря на то, что по построенному полю направлений возможно провести гладкую кривую, пересекающую ли-
нию y = 2x −e2 с угловым коэффициентом касательной k = 2 в точке пересечения, интегральные кривые не пересекают прямую y = 2x −e2 , поскольку сама эта прямая является интегральной кри-
вой, и в случае пересечения ее с другой интегральной кривой через точку пересечения проходили бы как минимум две интегральные кривые, а это противоречило бы условию единственности решения задачи Коши, которое должно соблюдаться всюду в области
2x − y > 0 .
Таким образом, метод изоклин позволяет по заданному полю направлений для уравнения, разрешенного относительно производной, определить приблизительное поведение семейства интегральных кривых данного уравнения.
37
10. РЕШЕНИЕ ТИПОВЫХ ЗАДАЧ
Задача 1. Определить тип дифференциального уравнения 2x2 yy′+ y2 = 2 и найти его общий интеграл.
Решение. Это уравнение с разделяющимися переменными. Разделим обе части уравнения на x2 (2 − y2 )≠ 0 и перепишем его в виде
|
|
|
|
|
|
2 ydy |
|
= |
dx |
. |
|
|
|
|
|
|
2 − y2 |
|
|
x2 |
|
Проинтегрировав, получим |
|
|
|
|
||||||
ln |
|
y2 −2 |
|
= 1 |
+C , |
или y2 −2 = C e1 x . |
||||
|
|
|||||||||
|
|
|
|
x |
|
|
1 |
|||
|
|
|
|
|
|
|
|
|
|
|
Задача 2. Определить тип |
|
дифференциального уравнения |
||||||||
xy′+2 y + x5 y3ex = 0 и найти его общий интеграл.
Решение. Это уравнение Бернулли. Решим его, выполнив замену z = y12 , y ≠ 0 . С учетом того, что z′ = −2 yy3′ , получим −x2 z′+ + 2z = −x5ex . Это линейное неоднородное уравнение. Решим соответствующее однородное уравнение
xdz +2z = 0 , −2dx
или
dzz = 4 dxx ln z = 4ln x +ln C z = Cx4 .
Методом вариации постоянной найдем C (x):
−x2 (C′x4 +4x3C) +2Cx4 = −x5ex C′ = 2ex C = 2ex +C1 ,
тогда z = 2x4ex +C1x4 .
38
Общее решение имеет вид y−2 = 2x4ex +C x4 . |
|
||
|
|
1 |
|
Заметим, |
что y = 0 является |
частным решением |
исходного |
уравнения. |
|
|
|
Задача 3. |
Определить тип |
дифференциального |
уравнения |
y = 3xy′−7 y′3 |
и найти его общий интеграл. |
|
|
Решение. |
Это уравнение Лагранжа. Введя параметр |
p = y′ , по- |
|
лучим y = 3xp −7 p3. Возьмем полный дифференциал, учитывая, что dy = pdx : pdx = 3pdx +3xdp −21p2dp или −2 pdx = (3x −21p2 )dp .
Перепишем последнее уравнение в виде
dpdx + 23xp = 212p .
Это линейное неоднородное уравнение. Решим соответствующее однородное уравнение
dx |
+ |
3x |
= 0 x = |
C |
. |
|
dp |
2 p |
|
||||
|
|
p |
3/ 2 |
|
||
Методом вариации постоянной найдем C (x):
|
C′ |
− |
3 |
|
|
C |
+ |
3 |
|
|
C |
= |
21 |
p C′ = |
21 |
p5/ 2 C = 3 p7 / 2 +C |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
p |
|
3/ 2 2 |
|
p |
|
3/ 2 2 |
|
p |
|
3/ 2 |
|
2 |
|
|
|
|
|
|
2 |
|
1 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
C1 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = |
3 p |
|
+ |
|
|
|
; |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
3/ 2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3xp − |
7 p |
3 |
. |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = |
|
|
|
||||||||
|
|
|
Задача 4. |
Определить тип |
|
дифференциального уравнения |
|||||||||||||||||||||||||
y − xy′ = 4 y′ |
и найти его общий интеграл. |
||||||||||||||||||||||||||||||
|
|
|
Решение. |
Это уравнение |
Клеро. |
|
Перепишем его в виде |
||||||||||||||||||||||||
y = x y′ − 4 y′ . Положив y′ = p , получим |
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = xp −4 |
p . |
|
|
(10.1) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
39 |
Продифференцировав уравнение (10.1) и заменив dy на pdx , найдем
|
|
pdx = pdx + xdp − 2dp , |
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
dp x − |
|
|
|
= 0 . |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Приравняв нулю первый множитель, |
получим dp = 0 , |
откуда |
||||||||||||||||||||||||
p = C . Общее решение исходного уравнения есть y = Cx −4 |
C . |
|||||||||||||||||||||||||
Приравняв нулю второй множитель, получим x = |
2 |
. |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
Исключив р из этого уравнения и из уравнения (10.1), получим |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
y = − |
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
– решение исходного уравнения. |
|
|
|
|
|
|
|
|
|
( |
) |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Задача 5. Решить задачу Коши: x( y′− y) = ex , y 1 = e . |
|
|
||||||||||||||||||||||||
Решение. Это линейное неоднородное уравнение. Разделим обе |
||||||||||||||||||||||||||
части уравнения на |
x ≠ 0 и решим соответствующее однородное |
|||||||||||||||||||||||||
уравнение y′ − y = 0, |
или |
|
dy |
= dx ln |
|
y |
|
= x +C , или y = C ex . |
||||||||||||||||||
|
|
|
||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
||
Методом вариации постоянной найдем C1 (x): |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
C′ex +C ex −C ex |
= |
ex |
, или |
C′ = 1 |
C = ln |
|
x |
|
+C |
2 |
. |
|
||||||||||||||
|
|
|
||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||
1 |
1 |
1 |
|
|
x |
|
|
|
1 |
|
|
|
x |
1 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Тогда общее решение будет иметь вид y = ex ln x +C2ex .
Частное решение найдем, подставив начальные данные в общее решение:
e = C2e C2 =1 y = ex (ln x +1) .
40
