Диф-уравнения-1 порядка
.pdfЗадание 3.8. xy0 |
= y + |
|
2x2 |
ln y |
|
|
|||||||||
y |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
x. |
|
|
|||
|
|
|
2 |
|
2 |
|
|
|
|
|
|
||||
Задание 3.9. y0 = |
|
y |
|
x |
|
|
|
y(1) = 1. |
|||||||
|
2 |
|
2 |
è |
|||||||||||
|
y |
|
2x |
|
|
|
|
|
|
|
|||||
Задание 3.10. y0 |
= |
|
2x2 |
xy +22y2 |
. |
|
|||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
2xy + x |
|
|
|||||||
Задание 3.11. y0 = |
|
|
|
5xy + 2y2 |
|
|
|||||||||
|
5x2 |
+ xy y2 |
. |
|
|
||||||||||
Задание 3.12. y0 |
= |
2x2 |
+ 2xy + y2 |
è |
y(1) = 4. |
||||||||||
|
|
|
|
|
|
xy + 2x2 |
|
|
|||||||
Задание 3.13. y0 |
= |
3 |
x2 |
|
xy + y2 |
|
y(1) = 4. |
||||||||
|
|
|
|
|
2 |
|
è |
||||||||
|
|
|
|
|
|
xy x |
|
|
|
|
|||||
Задание 3.14. y0 |
= |
y3 + 3x2y |
è y(2) = 1. |
||||||||||||
|
x3 + xy2 |
|
4. Линейные уравнения первого порядка.
Линейными дифференциальными уравнениями первого порядка называют уравнения линейные (т.е. первой степени) относительно неизвестной функции и ее производной:
y0 + p(x)y = q(x); |
(4:1) |
где p(x) и q(x) непрерывные функции.
Если q(x) = 0, то уравнение (4.1) называется линейным однород-
íûì. В противном случае (когда q(x) 6= 0) линейным неоднород-
íûì.
В линейном однородном уравнении переменные разделяются. В уравнении y0 +p(x)y = 0 перенесем слагаемое p(x)y в правую часть. Разделив
переменные и проинтегрировав получившееся уравнение, получим
R
y(x) = Ce p(x)dx; (4:2)
где C константа интегрирования.
Для решения неоднородного уравнения можно применить метод вариации произвольной постоянной. Будем считать, что в решении (4.2) однородного уравнения C является не константой, а функцией от x: C(x).
11
Дифференцируя выражение (4.2), получим, что
dy(x) C(x) |
e R p(x)dx C(x)p(x)e R p(x)dx: |
|
||
dx |
= |
dx |
(4:3) |
Подставляем функцию и ее производную в неоднородное уравнение (4.1):
C(x) |
e R p(x)dx C(x)p(x)e R p(x)dx + p(x)e R p(x)dx = q(x): |
(4:4) |
dx |
R
Сократив слагаемые C(x)p(x)e p(x)dx и умножив обе части уравнения
R
íà e p(x)dx, получим уравнение для определения функции C(x):
dC(x) |
= q(x)eR p(x)dx: |
(4:5) |
dx |
Проинтегрировав это уравнение и подставив результат в выражение (4.2), получим решение неоднородного линейного дифференциального уравнения первого порядка (4.1):
Z
R R
p(x)dx ~ p(x)dx
y(x) = e C + q(x)e dx : (4:6)
Здесь ~
C не зависящая от x константа интегрирования уравнения.
Другой метод решения неоднородного линейного дифференциального уравнения первого порядка заключается в том, что мы представляем искомую функцию y(x) в виде произведения двух функций:
y(x) = u(x)v(x): |
(4:7) |
Подставляя выражение (4.7) в уравнение (4.1), приводим это уравнение к виду:
v(x) |
dx |
+ p(x)u(x) + |
u(x) dx |
q(x) = 0 |
(4:8) |
||
|
|
du(x) |
|
|
dv(x) |
|
|
Приравняв нулю каждое из выражений, стоящее в квадратных скобках (поскольку неизвестных функций у нас две, то для их определения требуется два уравнения), получим линейное однородное уравнение для функции u(x):
du(x) |
+ p(x)u(x) = 0 |
(4:9) |
|
dx |
|||
|
|
12
и уравнение с разделяющимися переменными для функции v(x):
u(x) |
dv(x) |
q(x) = 0: |
(4:10) |
dx |
Решив эти уравнения, получим решение неоднородного линейного дифференциального уравнения первого порядка, в точности совпадающее с выражением (4.6).
Отметим, что и метод вариации постоянной, и представление искомой функции в виде произведения по сути две разные формы одного и того же метода решения уравнения. Как правило, ни один из этих двух способов не дает каких-либо преимуществ, и выбирают тот способ, который больше нравится.
|
|
|
|
dy |
|
y |
|
|
|
|
|
|
Пример 4.1. Решите уравнение dx |
+ |
|
= sin x . |
|
||||||||
x 1 |
|
|||||||||||
Решение . Решим однородное уравнение |
dy + |
y |
|
= 0. Перенесем второе сла- |
||||||||
x 1 |
||||||||||||
|
dy |
y |
dx |
dy |
dx |
|||||||
гаемое в правую часть dx = |
|
|
, разделим переменные y |
= x и проинте- |
||||||||
x |
1 |
|||||||||||
грируем |
|
|
. После преобразований получим решение одно- |
|||||||||
|
ln jyj = ln jx 1j + ln C |
|
|
|
|
|
|
|
родного уравнения y0 = |
C |
|
|
|
|
|
|
|
|
|
|
||||
x 1 |
. Для решения неоднородного уравнения будем счи- |
||||||||||||||
тать, что C функция переменной x и будем искать решение в виде y = |
C(x) |
||||||||||||||
x 1 |
. Âû- |
||||||||||||||
числим производную y0 = |
C0(x)(x 1) 2 |
C(x) |
и подставим y è y0 в неоднородное |
||||||||||||
|
|
||||||||||||||
|
|
|
|
(x 1) |
|
|
|
|
|
|
|
|
|||
уравнение. Получим |
C0(x 1) C |
+ |
|
C |
= sin x, откуда C0 = (x |
|
1) sin x. Èí- |
||||||||
(x |
1)2 |
||||||||||||||
|
(x 1)2 |
|
|
|
|
|
|
теграл для нахождения C(x) берется по частям: u = x 1, dv = sin xdx, du = dx,
та интегрирования). |
|
R |
~ |
|
|
~ ~ |
|
|
|
|
C(x) в решение. Итак, общее решение уравнения |
||||
v = cos x è C(x) = |
|
(x |
1) sin xdx = (x 1) cos x + sin x + C (C констан- |
||||
|
|
Подставим |
|
||||
y = (x 1) cos x + sin x + C |
|
|
|||||
|
x 1 |
|
|
|
. |
|
|
Отметим, что многие дифференциальные уравнения, которые формально не являются линейными, путем замены переменных могут быть сведены к линейным. Одним из распространенных способов замены поменять местами переменные x и y, т.е. искать зависимость не y = y(x),
à x = x(y). |
|
|
|
|
|
|
|
|
|
|
|
|
Пример 4.2. Решите уравнение |
dy |
+ |
1 |
|
|
= 0. |
|
|
|
|
|
|
dx |
2x y |
2 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
dy |
|
dx |
|
1 |
|
|
|
|
|
|
. В результате |
|
|
|||||
|
|
|
|
|
|
|
= |
|
|
|||
Будем считать аргументом y, а функцией x. Тогда с учетом того, что |
dx |
dy |
|
|
||||||||
|
|
dx |
|
|
|
|
|
|
|
|
||
исходное уравнение преобразуется к виду |
+ x = y |
|
|
получили неод- |
||||||||
dy |
|
|
13
нородное линейное дифференциальное уравнение относительно функции x(y), инте-
грирование которого не представляет труда (проделайте это самостоятельно).
Рассмотрим еще один тип уравнений, который с помощью замены переменных сводится к линейным, уравнение Бернулли, имеющее вид
dy |
+ p(x)y = q(x)yn; |
(4:11) |
|
dx |
|||
|
|
которое заменой переменных z = y1 n сводится к линейному. Действи-
тельно, дифференцируя z(x) как сложную функцию от y, находим
dz |
= |
dz |
|
dy |
= |
d |
(y1 n) = (1 |
n)y n |
dy |
(4:12) |
|
|
|
|
|
||||||
dx |
dy |
dx |
dy |
dx |
Разделив обе части уравнения (4.11) на y n и используя выражение
(4.12), получаем для функции z(x) линейное неоднородное дифферен-
циальное уравнение первого порядка.
|
1 dz |
|
|
|||
|
|
|
|
+ p(x)z = f(x) |
|
(4:13) |
1 n |
dx |
|
||||
Интегрирование уравнения (4.13) дает |
|
|
||||
y(x) = e R p(x)dx C~ + Z f(x)eR p(x)dxdx 1=(n 1) |
: |
(4:14) |
Уравнение Бернулли может быть решено и представлением искомой функции в виде произведения двух функций.
Задания для самостоятельного решения.
Решите уравнение. Если заданы начальные условия, решите задачу Коши.
Задание 4.1.
Задание 4.2.
Задание 4.3.
Задание 4.4.
Задание 4.5.
y0 |
|
y |
= x2 è y(2) = 5. |
||||
x 1 |
|
||||||
y0 |
|
|
y |
2 |
. |
||
+ |
|
|
= ex +2x |
|
|||
x + 1 |
|
||||||
y0 |
|
2y |
= x 2 è y(1) = 10. |
||||
x 3 |
|
||||||
y0 2y ctg x = cos x. |
|||||||
y0 |
+ |
|
y |
= ln x. |
|
||
x + 1 |
|
|
14
Задание 4.6. y0 |
|
|
y |
|
= ln(x + 1) è y(0) = 1. |
|||||||||||
x + 1 |
||||||||||||||||
Задание 4.7. y0 |
|
2xy |
|
= arctg x è y(1) = 2 + 1. |
||||||||||||
x2 + 1 |
|
|||||||||||||||
Задание 4.8. y0 |
+ |
2xy |
|
|
= 3x + 1. |
|
||||||||||
2 |
|
|||||||||||||||
|
|
|
x + 3 |
|
|
|
|
|
|
|
||||||
Задание 4.9. y0 |
3x2y = (x2 + 1)ex3 . |
|
||||||||||||||
Задание 4.10. y0 |
+ y tg x = sin3 x. |
|
||||||||||||||
Задание 4.11. y0 |
+ 3y tg x = sin x è y( =3) = 3. |
|||||||||||||||
Задание 4.12. y0 |
|
|
4y |
|
= x 2 è y( 2) = 7. |
|||||||||||
x + 3 |
||||||||||||||||
Задание 4.13. y0 |
|
|
2xy |
|
|
|
= 2x 7 è y(3) = 5. |
|||||||||
x2 |
|
4 |
|
|||||||||||||
Задание 4.14. y0 |
|
|
|
y |
|
|
|
|
|
= tg x è y( =3) = 1. |
||||||
|
|
|
3 |
|
|
|
|
|
||||||||
|
|
|
|
|
|
sin x |
|
|
|
|||||||
|
|
cos |
x |
|
|
|
||||||||||
Задание 4.15. y0 |
|
|
(x + 3)y |
|
||||||||||||
|
|
= 2. |
|
|||||||||||||
x2 + 6x + 8 |
|
|||||||||||||||
Задание 4.16. 2y0 |
|
3y |
|
|
|
= (x2 3x + 5)y3 |
. |
|||||||||
x + 1 |
|
5. Уравнения в полных дифференциалах.
Если мы рассматриваем уравнения первого порядка в дифференциальной форме
M(x; y)dx + N(x; y)dy = 0; |
(5:1) |
то может оказаться, что его левая часть является полным дифференциалом некоторой функции двух переменных U(x; y):
dU(x; y) = @U@x dx + @U@y dy = M(x; y)dx + N(x; y)dy = 0;
то есть функции M(x; y) и N(x; y) являются частными производными функции двух переменных U(x; y):
M(x; y) = |
@U |
; N(x; y) = |
@U |
: |
(5:2) |
|
@x |
@y |
|||||
|
|
|
|
Если выполнены условия теоремы о равенстве смешанных производных, порядок взятия частных производных не влияет на результат, то есть
15
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
@ |
|
@U |
= |
|
@ |
|
@U |
, то функции M(x; y) и N(x; y) должны удовле- |
||||||
|
@y |
@x |
|
@x |
@y |
|
|
|
|
|
|
|
|||
творять условию: |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
@M(x; y) |
= |
@N(x; y) |
(5:3) |
||
|
|
|
|
|
|
|
|
|
|
@y |
|
@x |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
В этом случае говорят, что уравнение (5.1) является уравнением в полных дифференциалах. Общий интеграл этого уравнения записывается в виде
U(x; y) = C; |
(5:4) |
а удовлетворяющий начальному условию y(x0) = y0 частный интеграл в âèäå
U(x; y) = U(x0; y0): |
(5:5) |
Найдем сначала частный интеграл уравнения (5.1), удовлетворяющий условию y(x0) = y0. Проинтегрируем, например, первое из уравнений
(5.2) по x, используя при этом понятие определенного интеграла с пере-
x
R
менным верхним пределом: U(x; y) = M(x; y)dx + '(y).
x0
При вычислении интеграла считается, что y = const, и интегрирование по x выполняется как от обычной функции одной переменной. Поэтому, в общем случае константа интегрирования зависит от y: '(y). Для нахождения функции '(y) воспользуемся вторым уравнением (5.2). Продифференцировав полученное выражение по y (считая при этом x = const), получим:
|
@y |
|
= |
@y |
x0 |
M(x; y)dx + '(y)! = x0 |
@y |
|
dx + |
|
|
dy |
= |
||||||
|
@U(x; y) |
|
@ |
x |
x |
@M(x; y) |
|
|
|
d'(y) |
|
||||||||
|
x @N(x; y) |
|
R d'(y) |
R |
|
d'(y) |
|
|
|
|
|
||||||||
здесьR0 |
|
|
|
|
|
|
|
|
|||||||||||
|
мы воспользовались соотношением (5.3). В итоге уравнение для |
||||||||||||||||||
= |
|
@x |
|
dx + |
dy |
= N(x; y) N(x0 |
; y) + |
dy |
|
= N(x; y), |
|||||||||
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
отыскания функции '(y) будет иметь простой вид |
|
d'(y) |
= N(x0; y). |
||||||||||||||||
|
dy |
|
|
Важно, что левая часть полученного уравнения зависит только от пере- менной y (ведь x0 константа!). Проинтегрировав обе части уравнения
y |
|
по y в пределах от y0 до y, получим '(y) = yR0 |
N(x0; y)dy, что позволяет |
16
записать частный интеграл уравнения (5.1) в виде
x |
y |
ZZ
U(x; y) = M(x; y)dx + N(x0; y)dy = 0: |
(5:4) |
|
x0 |
y0 |
|
Если начать интегрирование со второго уравнения из (5.2), то совершенно аналогично можно получить, что
x |
y |
ZZ
U(x; y) = M(x; y0)dx + N(x; y)dy = 0: |
(5:5) |
|
x0 |
y0 |
|
Отметим, что в выражении (5.4) в первом интеграле переменная y во
время интегрировании по x считается параметром, а потом опять пе-
ременной. Аналогично, в выражении (5.5) во втором интеграле переменная x при интегрировании по y считается параметром, а потом опять переменной.
Пример 5.1. Найдите решение уравнения
(x2 3y + 2x)dx + (y 3x + sin y)dy = 0, удовлетворяющее начальному
условию y(2) = 0.
Решение . Äàíî M(x; y) = x2 3y + 2x, N(x; y) = y 3x + sin y. Это уравнение является уравнением в полных дифференциалах, если выполнено условие
|
@M(x; y) |
= |
@N(x; y) |
|
|
@M(x; y) |
= |
3, |
@N(x; y) |
= 3. Следовательно, за- |
|||||||||||||||||||
|
@y |
|
|
|
@x |
|
. Имеем |
|
|
@y |
|
|
|
@x |
|
||||||||||||||
дано уравнение в полных дифференциалах. Для нахождения функции |
U(x; y) âû- |
||||||||||||||||||||||||||||
Найдем функцию '(y). ДляR |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
x3 |
@U(x; y) |
2 |
||||||||||
этого найдем Rчастную производную |
|
||||||||||||||||||||||||||||
числим интеграл U(x; y) = |
|
M(x; y)dx = |
|
(x |
|
3y + 2x)dx = 3 |
3xy + x + '(y). |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@y |
|
и прирав- |
|
няем ее функции N(x; y). Итак, |
|
|
@U(x; y) |
|
= 3x + '0(y) = y 3x + sin y, откуда |
||||||||||||||||||||||||
|
|
@y |
|
|
|||||||||||||||||||||||||
'0(y) = y + sin y è '(y) = |
|
(y + sin y)dy = |
y2 |
cos y. Общий интеграл уравнения |
|||||||||||||||||||||||||
|
2 |
||||||||||||||||||||||||||||
|
|
|
|
( |
) = 3 3 R |
|
+ x |
|
|
+ 2 cos |
|
= |
|
|
|
|
|
|
|
||||||||||
имеет вид U x; y |
x3 |
|
|
|
2 |
|
y2 |
|
|
|
y |
|
|
C. |
|
|
|
|
|
||||||||||
|
xy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Найдем |
|
частный |
интеграл. |
|
Подставим |
|
â |
|
общий интеграл начальные условия: |
||||||||||||||||||||
|
|
|
23 |
|
|
|
2 |
|
02 |
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
2 |
|
|
||
U(2; 0) = |
3 3 2 0 + 2 + |
|
2 cos 0 = C. Значит, C = |
3 + 4 1 = 5 |
3 и частный |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
x3 |
|
|
|
|
2 |
|
|
y2 |
|
|
|
|
2 |
|
|
|
|
||||
интеграл задан уравнением |
3 3xy + x |
|
+ |
2 cos y = 53. |
|
|
|
|
Задания для самостоятельного решения.
Найдите общий интеграл уравнения
17
Задание 5.1. Задание 5.2.
Задание 5.3.
Задание 5.4.
(2x2 cos y y2 sin x 1)dx + (2 x3 sin y y2 sin x)dy = 0.
(sin y + y cos x)dx + (x cos y + sin x)dy = 0.
(2x ln y + x2)dx + (x2 3)dy = 0. y
(x y + 1)ex+ydx + (x y 1)ex+ydy = 0.
6. Основные типы дифференциальных уравнений первого порядка.
I. Уравнения с разделяющимися переменными. Уравнение y0 = f1(x) f2(y) называется уравнением с разделяющимися переменными . С помощью алгебраических преобразований приводим уравнение к виду
dy
f2(y) = f1(x)dx.
Обе части уравнения с разделенными переменными представляют собой дифференциалы некоторых функций. Поскольку равны дифференциалы функций, то с точностью до постоянной равны и сами функции. Следовательно, общий интеграл уравнением с разделяющимися перемен-
ными имеет вид |
|
|
1 |
dy = f1(x)dx + C. |
|
|
|
|
|
f2 |
(y) |
y0 = f(x; y) |
, ãäå |
|
|||
|
|
|
|
|
||||
II. |
Однородные уравнения. Уравнения |
f(x; y) |
||||||
|
R |
|
R |
|
однородная функция нулевого измерения называются однородными уравнениями.
С помощью замены переменной z(x) = |
y(x) |
||
|
|
||
x или, что тоже самое |
|||
|
y(x) = x z(x) привoдим его к уравнению с разделяющимися перемен-
íûìè.
уравнения первого порядка это уравнения линейные относительно неизвестной функции и ее производной: y0 + p(x)y = q(x); где p(x) и q(x) непре-
рывные функции. Если q(x) = 0, то уравнение называется однородным, åñëè q(x) 6= 0 неоднородным линейным уравнением.
Линейное однородное уравнение является уравнением с разделяю-
R
щимися переменными. Его решение имеет вид y(x) = Ce p(x)dx, ãäå
C = const.
18
Для решения неоднородного уравнения применим метод вариа-
öèè |
произвольной |
постоянной, |
|
считая, что |
C является |
функци- |
|||
åé îò x: C(x). |
Для определения функции |
C(x) получаем урав- |
|||||||
нение |
|
dC(x) |
= |
q(x)eR p(x)dx |
. |
Решение неоднородного |
линейно- |
||
|
dx |
|
|
|
|
|
го дифференциального уравнения первого порядка имеет вид
Решение |
|
R |
|
y(x) = e R p(x)dx |
|
C~ + |
q(x)eR p(x)dxdx , ãäå C~ константа. |
линейного неоднородного уравнения можно предствить в виде суммы общего решения однородного уравнения и частного решения
неоднородного уравнения: yон = yоо + yчн.
Уравнение dxdy + p(x)y = q(x)yn называ- åòñÿ уравнением Бернулли. Заменой переменных z = y1 n оно сводится
к линейному уравнению первого порядка..
V. Уравнения в полных дифференциалах. Уравнение
M(x; y)dx + N(x; y)dy = 0 называется уравнением в полных дифференциалах, если его левая часть является полным дифференциалом некоторой функции U(x; y), т.е. M(x; y)dx + N(x; y)dy = dU(x; y). Тогда
M(x; y) = @U ; N(x; y) |
= @U |
|
|
|
M(x; y) и N(x; y) должны |
||||||
|
|
|
@x |
|
@y и функции |
|
|
||||
удовлетворять условию: |
|
@M(x; y) |
= |
@N(x; y) |
|
||||||
|
|
|
|
||||||||
|
@y |
|
@x . |
||||||||
|
|
|
|
|
|
||||||
где для нахождения функции '(y) нужно |
|
|
x |
||||||||
|
|
R |
|||||||||
|
Решение этого уравнения ищется в виде |
U(x; y) = M(x; y)dx+'(y), |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
проинтегрировать выражение |
||
|
d'(y) |
= N(x0; y). |
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|
19