DRUN
.pdfгде u(x) новая функция, уравнение (1.20) преобразуется в дифференциальное уравнение с разделяющимися переменными.
|
|
du |
|
|
|
du |
= Á(u) ¡ u; |
||||
|
u + x |
|
|
= Á(u); x |
|
|
|||||
|
dx |
dx |
|||||||||
|
|
|
|
|
du |
|
dx |
||||
|
|
|
|
|
|
|
= |
|
|
: |
|
Интегрируя, получаем |
|
|
|
|
Á(u) ¡ u |
|
x |
||||
Z |
|
|
du |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||||
|
|
= ln jxj + C: |
|||||||||
|
Á(u) ¡ u |
Подставляя после интегрирования вместо u его значение u = xy , получим общее решение уравнения.
Отдельно нужно рассмотреть случай u = u0 , если Á(u0) ¡ u0 = 0 . Уравнение вида
P (x; y)dx + Q(x; y)dy = 0
является однородным, если функции P (x; y) и Q(x; y) - однородные функции одного и того же измерения.
Пример 2. Найти общее решение уравнения
(x2 ¡ y2)dx + 2xydy = 0:
P (x; y) = x2 ¡ y2
и
Q(x; y) = 2xy
однородные функции второго измерения.
Делаем подстановку y = xu; dy = xdu + udx . Заданное уравнение перепишется в
виде
x2(1 ¡ u2)dx + 2x2u(xdu + udx) = 0;
x2((1 + u2)dx + 2xudu) = 0:
Отсюда x = 0; или
(1 + u2)dx + 2xudu = 0: |
(1.22) |
x = 0 является решением исходного уравнения. Уравнение (1.23) - уравнение с разделяющимися переменными
Так как u = xy , то xy22
|
|
|
2udu |
|
+ |
dx |
= 0; |
|
|
|
|
|
|
||||
|
|
1 + u2 |
|
x |
||||
Z |
udu |
+ Z |
dx |
|
|
|||
2 |
|
|
|
= ln jCj; C 6= 0; |
||||
1 + u2 |
|
x |
ln(1 + u2) + ln jxj = ln jCj;
1 + u2 = Cx :
+ 1 = Cx . Окончательно
y2 + x2 = C; x
11
или |
|
|
|
(1.23) |
y2 + x2 = Cx; ; C 6= 0 |
||||
общий интеграл уравнения. Решение |
x = 0 не входит в (1.23) ни при каких значениях |
|||
константы C: |
|
|
|
|
Ответ. |
|
|
|
|
y2 + x2 = Cx; C 6= 0; x = 0: |
|
|||
Уравнение вида |
µkx + my + n¶ |
|
||
y0 = f |
|
|||
|
|
ax + by + c |
|
|
приводится к однородному с помощью переноса начала координат в точку пересечения прямых
ax + by + c = 0;
kx + my + n = 0:
Если x0; y0 решение системы, то в переменных x1; y1 ,
x1 = x ¡ x0; y1 = y ¡ y0
уравнение однородное. Если же прямые параллельны, то делается замена
ax + by + c = z:
В некоторых случаях уравнение можно привести к однородному заменой y = zm . Число m подбирается из условия, чтобы полученное уравнение было однородным.
Пример 3. Решить уравнение x2yy0 + y4 = 2x2:
После замены y = zm получаем уравнение x2mz2m¡1z0+z4m = 2x2 , у которого степени однородности всех его слагаемых должны быть равны, то есть
2 + (2m ¡ 1) = 4m = 2
. Отсюда m = 12 . Замена y = z1=2 приводит уравнение к однородному.
z0 = ¡2 z2 + 4: x2
Задачи
Решить уравнения.
1. |
2x3y0 = y(2x2 ¡ y2) . |
|||||||
3. |
(x2 + y2)y0 = 2xy . |
|
|
|||||
5. |
xy0 = y ¡ xey=x . |
2 |
|
|||||
7. |
y0 |
= 2 µx + y ¡ 1¶ |
. |
|
||||
|
|
|
|
y + 2 |
|
|
||
9. |
y0 |
= |
y + 2x |
+ tg |
y ¡ 2x |
. |
||
x + 1 |
|
|||||||
|
|
|
|
x + 1 |
.
2. (x ¡ y) dx + (x + y) dy = 0 .
4. xy0 ¡ yp= (x + y) ln x +x y . 6. xy0 = x2 ¡ y2 + y .
8. (y + 2) dx = (2x + y ¡ 4) dy .
10. (x ¡ y ¡ 1) + (y ¡ x + 2)y0 = 0 .
11.Найти кривую, у которой точка пересечения любой касательной с осью абсцисс одинаково удалена от точки касания и от начала координат.
12.Найти кривую, у которой расстояние любой касательной от начала координат равно абсциссе точки касания.
12
Ответы
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
1. |
x = §y |
ln Cx , |
y = 0 . 2. |
ln(x |
|
+ y |
) = C ¡ 2arctg |
|
|
. 3. y |
|
¡ x |
|
= Cy , |
y = 0 . |
||||||||||||||||||||||||||
|
x |
|
|
|
|||||||||||||||||||||||||||||||||||||
4. |
ln |
|
x + y |
= Cx . 5. |
y = ¡x ln ln Cx . |
6. |
arcsin |
y |
= ln Cx sign x , |
y = §x . 7. |
y + 2 = |
||||||||||||||||||||||||||||||
|
x |
|
|
x |
|
||||||||||||||||||||||||||||||||||||
C exp |
µ¡ |
2arctg |
y + 2 |
¶ |
. 8. |
(y + 2)2 = C(x + y |
¡ |
1) |
, |
y = 1 |
¡ |
x |
. 9. |
sin |
y ¡ 2x |
= C(x + 1) |
. |
||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
x |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x + 1 |
|
||||||||||||||||
10. |
(y ¡ x + 2) |
2 ¡ |
|
|
|
11. |
y = C(x |
2 |
|
2 |
) . 12. |
|
x |
2 |
+ y |
2 |
= Cx . |
|
|
|
|
|
|
|
|
||||||||||||||||
|
+ 2x = C . |
|
+ y |
|
|
|
|
|
|
|
|
|
|
|
1.5Линейные дифференциальные уравнения первого порядка
Линейным дифференциальным уравнением первого порядка называется уравнение, линейное относительно неизвестной функции и ее производной. Оно имеет вид
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
+ p(x)y = q(x): |
(1.24) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|||||||
|
|
Уравнение |
|
|
|
|
|
|
dy |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ p(x)y = 0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|||||
называется линейным однородным дифференциальным уравнением. |
|
||||||||||||||||||||
|
|
Методы решения |
1. Метод Бернулли. |
|
|||||||||||||||||
|
|
Ищем решение уравнения (1.24) в виде произведения двух функций, зависящих от x : |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
y(x) = u(x)v(x): |
(1.25) |
|||||||
|
|
Продифференцируем |
обе части |
(1.25). Найденное значение |
производной |
||||||||||||||||
|
dy |
= u |
dv |
+ v |
du |
и функцию y(x) = u(x)v(x) подставляем в уравнение (1.24): |
|||||||||||||||
dx |
|
|
|||||||||||||||||||
dx |
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
u |
dv |
+ v |
du |
|
+ p(x)uv = q(x); |
|
|||||||
|
|
|
|
|
|
|
|
dx |
dx |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
или |
|
|
|
|
du |
|
|
|
|
|
|
dv |
|
|
|||||||
|
|
|
|
|
|
|
|
v |
|
+ u[ |
+ p(x)v] = q(x): |
(1.26) |
|||||||||
|
|
|
|
|
|
|
|
dx |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
||||||
Функцию v |
выберем так, чтобы она являлась частным решением уравнения |
dxdv + p(x)v = 0:
Тогда из (1.26)
v dudx = q(x):
Получили два уравнения с разделяющимися переменными для деления функций v(x); u(x):
Разделяя переменные и интегрируя уравнение |
dvv = ¡p(x)dx |
||||||
ln jvj = ¡ Z |
p(x)dx; |
v(x) = e¡ R p(x)dx: |
|||||
Найденное значение v подставляем в уравнение (1.28). |
|||||||
e¡ R p(x)dx |
du |
|
du |
= q(x)eR p(x)dx |
|||
|
= q(x); или |
|
|
||||
dx |
dx |
(1.27)
(1.28)
последовательного опре-
, получим
(1.29)
:
13
Отсюда |
Z |
du = Z ¡q(x)eR p(x)dx¢dx + C; |
|
|
|
|
|||
|
|
u(x) = Z |
q(x)eR p(x)dxdx + C: |
(1.30) |
Подставляя найденное значение u(x); |
v(x) в формулу (1.25), находим общее решение |
|||
уравнения (1.24) |
|
|
|
|
µZ ¶
R R
y(x) = u(x)v(x) = e¡ p(x)dx q(x)e p(x)dxdx + C :
2. Метод вариации постоянной (метод Лагранжа).
Вначале находим общее решение однородного уравнения (уравнение (уравнение с разделяющимися переменными)
y0 + p(x)y = 0; |
|
y = Ce¡ R a(x)dx: |
|
Общее решение исходного уравнения (1.24) ищем в виде |
|
y = C(x)e¡ R p(x) dx; |
(1.31) |
где C(x) произвольная дифференцируемая функция g(x):
После подстановки в (1.24), получаем уравнение для определения функции C(x)
C0(x) = q(x)eR p(x)dx;
Отсюда |
q(x)eR p(x) dx dx; |
C(x) = C + Z |
где C –произвольная постоянная. Подставляя найденнную функцию C(x) в (1.31), получаем общее решение исходного уравнения (1.24)
y = ·C + Z q(x)eR p(x) dx dx¸e¡ R p(x)dx:
Пример 1. Решить уравнение
y0 + y cos x = 2xe¡ sin x:
Вначале находим общее решение однородного уравнения y0 + y cos x = 0;
y = Ce¡ sin x:
Решение исходного уравнения ищем в виде:
y = C(x)e¡ sin x:
Имеем:
C0(x)e¡ sin x ¡ cos xC(x)e¡ sin x + C(x)e¡ sin x cos x = 2xe¡ sin x; C0(x) = 2x;
14
C(x) = x2 + C:
Общее решение:
y = (x2 + C)e¡ sin x:
Некоторые уравнения становятся линейными, если переменную x считать функцией, а y аргументом. Например,
dy |
1 |
|
|
|
= |
|
: |
dx |
x + y2 |
Отсюда
dxdy = x + y2
- линейное уравнение относительно неизвестной функции x(y):
1.6Уравнение Бернулли
Уравнение Бернулли имеет вид |
|
|
dy |
+ p(x)y = q(x)yn; n 6= 0; 1: |
(1.32) |
dx |
При n = 1 уравнение (1.32) линейное, а при n = 0 - c разделяющимися переменными. В остальных случаях уравнение (1.32) преобразуется в линейное относительно новой
функции z(x)
z = y1¡n
Замечание. Уравнение Бернулли можно решать также как линейное уравнение с помощью подстановки y = u(x)v(x) (методом Бернулли), не вводя вспомогательную функцию z .
Пример 1. Найти общее решение уравнения Бернулли
xy0 ¡ y = ¡3x2y2
Решение. Преобразуем уравнение к виду y0 ¡ xy = ¡3xy2 . Воспользуемся подстановкой y = uv; y0 = u0v + v0u . Тогда
u0v + v0u ¡ uvx = ¡3xu2v2;
далее
u0v + u(v0 ¡ xv ) = ¡3xu2v2:
Приравниваем нулю выражение, стоящее в скобках
|
|
|
|
|
|
dv |
¡ |
v |
= 0; или |
dv |
= |
v |
: |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
dx |
x |
dx |
x |
||||||||
Разделяем переменные |
dv |
= dx |
и интегрируем: |
|
|
|
|
|
|||||||||
|
dv |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z |
dv |
= Z |
dx |
; ! ln jvj = ln jxj; ! v = x: |
||||||||||||
|
|
|
|
|
|
|
|||||||||||
|
|
v |
|
x |
Найденное v = x подставим в уравнения для определения функции u(x) :
dudxx = ¡3xu2x2:
15
Сокращаем на x , разделяем переменные и интегрируем
|
|
|
|
|
|
|
u2 |
= ¡3x2dx; |
Z |
|
|
u2 = ¡3 Z |
x2dx; (u 6= 0): |
|||||||
|
|
|
|
|
|
|
du |
|
|
|
|
|
du |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
3x3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
1 |
|
|||||||
|
|
|
|
|
¡ |
|
= ¡ |
|
+ C; |
|
|
= C ¡ x3; ! u = |
|
|
||||||
|
x |
|
|
|
u |
3 |
u |
C ¡ x3 |
||||||||||||
y = uv = |
|
|
общее решение уравнения. При |
u = 0 , y = 0 является решением |
||||||||||||||||
|
|
|
||||||||||||||||||
C ¡ x3 |
|
|||||||||||||||||||
уравнения и не следует из общего решения. |
|
|
|
|||||||||||||||||
Ответ. y = |
|
x |
|
|
; y = 0: |
|
|
|
|
|
|
|
|
|
||||||
C¡x |
3 |
|
|
|
|
|
|
|
|
|
Замечание. Окончательное нахождение общего решения рассмотренных типов уравнений (уравнений с разделяющимися переменными, однородные, линейные, Бернулли) сводится к вычислению интегралов. В таком случае говорят, что решение определяется (записывается) в квадратурах.
Рассмотрим уравнение общее решение которого в общем случае не записывается в квадратурах.
1.7 |
Уравнение Риккати |
|
|
|
||
Уравнение |
|
|
|
|||
|
y0 = p(x)y2 + q(x)y + r(x) |
(1.33) |
||||
называется уравнением Риккати. |
|
|
|
|||
При |
p(x) = 0 уравнение линейное, при r = 0 - уравнение Бернулли. |
|
||||
Свойства уравнения Риккати. |
|
|
|
|||
1) Заменой |
|
|
|
|||
|
1 |
|
|
|
|
|
|
y = § |
|
z(x) |
|
|
|
|
p(x) |
|
|
|
||
уравнение принимает вид |
|
|
|
|||
|
z0 = §z2 § pr: |
|
|
|
||
2) При замене |
q |
|
||||
|
y = u(x) + ®(x); ®(x) = ¡ |
|
||||
|
|
|
|
|||
|
2p |
|
у полученного уравнения коэффициент при u будет равен нулю. Коэффициент при u2 не изменяется.
u0 = pu2 + r + p®2 + q® ¡ ®0:
Комбинируя обе подстановки, уравнение Риккати можно привести к виду
y0 = §y2 + R(x):
3) Если известно частное решение y = y1(x) , то заменой y = z + y1 уравнение Риккати сводится к уравнению Бернулли. В этом случае общее решение уравнения Риккати выражается в квадратурах.
Пример 1. Найти общее решение уравнения
y0 ¡ 2xy + y2 = 5 ¡ x2:
Имеем уравнение Риккати. Частное решение можно пробовать искать например, в виде y1 = anxn + an¡1xn¡1 + : : : ; y1 = aebx; y1 = xan :
16
Ищем частное решение в виде
y = ax + b:
После подстановки в уравнение и приравнивания коэффициентов при одинаковых степенях x в левой и правой частях полученного уравнения, получаем систему для определения
коэффициентов a и b
a + b2 = 5;
¡2a + a2 = 5:
Отсюда a = 1; b = §2 . Возьмем
y1 = x + 2:
Далее
y = z + x + 2:
После подстановки в исходное уравнение получаем относительно z уравнение Бернулли
z0 + 4z = ¡z2:
4
z = x + 2 + ce4x ¡ 1
его общее решение. Следует обратить внимание, что решение исходного уравнения y1 = x + 2 не содержится в полученном общем решении уравнения Бернулли.
Ответ:
4
z = x + 2 + ce4x ¡ 1; y = x + 2:
В общем случае общее решение уравнения Риккати не может быть выражено в квадратурах. Так для специального уравнения Риккати
y0 = ay2 + bx®
имеется схема определения тех значений ® при которых решение определяется в квадратурах. Для остальных значений ® доказано, что решение не может быть выражено в квадратурах от элементарных функций.
|
|
|
|
|
|
|
|
|
Задачи |
|
|
|
Решить уравнения. |
|
|
|
|
|
|
|
|
||||
1. |
y0 + ytgx = cos¡1 x |
. |
2. |
y0 |
¡ |
|
2xy = 3x2 |
2x4 |
. |
|||
2x(x |
2 |
|
|
|
|
3¡ |
|
|||||
3. |
|
+ y) dx = dy . |
4. |
y dx ¡ (2x + y ) dy = 0 . |
||||||||
5. |
(xy0 ¡ 1) ln x = 2y . |
|
6. |
(1 ¡2 |
2xy)y0 = y(y ¡ 1) . |
|||||||
7. |
y = x(y0 ¡ x cos x) . |
|
8. |
(sin |
|
y + xctgy)y0 |
= 1 . |
|||||
Решить задачу Коши. |
|
|
|
|
|
|
|
|||||
9. |
y0 ¡ ytgx = cos¡3 x , |
y(0) = 0 . |
|
|
|
10. y0 + y cos x = cos x , y(0) = 1 . |
||||||
Решить уравнения Бернулли. |
|
|
|
|
|
|
||||||
11. |
(x3 + ey)y0 = 3x2 . |
|
|
12. |
|
y0 + 2xy = y2ex2 . |
||||||
13. |
y0 ¡ 2yex = 2p |
yex . |
|
14. |
|
y0 ¡ y cos x = y2 cos x . |
||||||
15. |
y0 + 2xy = 2x3y3 . |
|
|
16. |
2y2y0 + y3 + x = 0 . |
17
Решить уравнения Риккати.
17. xy0 ¡ (2x + 1)y + y2 = ¡x2 . 18. y0 + 2yex ¡ y2 = e2x + ex .
19.Найти решение уравнения y0 sin x¡y cos x = ¡x¡2 sin2 x; которое стремится к нулю при x ! 1 .
20.Найти решение уравнения y0 sin 2x = 2(y + cos x) , которое остается ограниченным
¼.
2
21. Найти кривые, у которых площадь треугольника, ограниченного касательной, осью абсцисс и отрезком от начала координат до точки касания, есть величина постоянная,
равная |
a2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
22. Даны два различных решения y1 и |
y2 |
линейного уравнения первого порядка. |
|||||||||||||||||||||||||||||||||||||||
Выразить через них общее решение этого уравнения. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
23. С помощью дифференцирования решить уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(x) = Z0x(x ¡ t)y(t) dt ¡ Z0x y(t) dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ответы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1. |
y |
|
|
C |
|
|
|
x |
|
x |
. 2. |
y |
|
|
|
|
C exp(x2) + x3 |
. |
3. |
y = C exp(x2) |
¡ |
x2 |
¡ |
1 |
. 4. |
x = (y + C)y2 |
. |
||||||||||||||||
|
= |
|
cos2 |
|
+ sin |
|
|
= |
2 |
|
|
|
y = 0 ; |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
5. |
y = C ln |
|
|
x ¡ ln x . 6. |
|
(y ¡ 1) x = y ¡ ln Cy ; |
y = 1 . 7. |
y = x(C + sin x) . |
|||||||||||||||||||||||||||||||||||
8. |
x = (C ¡ cos y) sin y . |
9. |
|
|
y = sin x cos¡2 x . 10. |
y = 1 . |
11. x3e¡y |
= C + y . 12. |
y = |
||||||||||||||||||||||||||||||||||
exp(¡x2)(C ¡ x)¡1 . |
13. |
p |
|
|
+ 1 = C exp exp x . |
14. y = |
|
|
|
|
|
|
1 |
|
|
|
|
. 15. |
y¡2 = |
||||||||||||||||||||||||
y |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
C exp( |
¡ |
sin x) |
¡ |
1 |
|||||||||||||||||||||||||||||||||||||||
|
|
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
. 16. y3 = Ce¡x ¡x+1 . 17. y = x+x(x+C)¡1 ; |
y = x . 18. y = ex ¡(x+C)¡1 ; |
|||||||||||||||||||||||||||||||||||
Ce2x |
+x2 + |
|
|
||||||||||||||||||||||||||||||||||||||||
2 |
|
||||||||||||||||||||||||||||||||||||||||||
y |
¡ |
ex |
. 19. |
y +x¡1 sin x |
. 20. |
|
y = tgx |
¡ |
cos¡1 x |
. 21. |
xy = a2 +Cy3 |
. 22. |
y = y |
+C(y |
2 ¡ |
y |
) |
. |
|||||||||||||||||||||||||
|
|
|
|
x |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
||||||||||||||||||||
23. |
|
y = ¡2e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.8Уравнения в полных дифференциалах
Уравнение
P (x; y)dx + Q(x; y)dy = 0 |
(1.34) |
называется уравнением в полных дифференциалах, если его левая часть является полным дифференциалом некоторой функции U(x; y) , то есть
P dx + Q dy = Ux0 dx + Uy0 dy = dU(x; y): |
(1.35) |
В предположении, что функции P (x; y) , Q(x; y) непрерывно-дифференцируемые в односвязной области, условие (1.35) справедливо тогда и только тогда, когда
|
|
Py0 = Qx0 : |
(1.36) |
|
В этом случае общий интеграл уравнения (1.34) имеет вид |
|
|||
|
U(x; y) = C; |
C = const: |
|
|
В соответствии с (1.35) |
|
|
|
|
U0 |
|
= P (x; y); |
U0 = Q(x; y): |
(1.37) |
x |
|
y |
|
18
Интегрируя первое равенство в (1.37) по x , считая y параметром, получаем |
|
|
U(x; y) = Z |
P (x; y)dx + '(y) = ©(x; y) + '(y): |
(1.38) |
'(y) – произвольная дифференцируемая функция. Подставляя U(x; y) из (1.38) во второе равенство из (1.37), получим уравнение для определения функции '(y) :
@©(x; y) + d' = Q(x; y): @y dy
Пример 1. Найти общий интеграл уравнения:
(x2 + y2 + x)dx + (2xy + ey + y)dy = 0;
P (x; y) = x2 + y2 + x; Q(x; y) = 2xy + ey + y;
@P (x; y) |
= 2y; |
@Q(x; y) |
= 2y: |
||
@y |
|
@x |
|||
|
|
Условие (1.36) выполняется - уравнение в полных дифференциалах. Из
Ux0 = P (x; y) = x2 + y2 + x;
U(x; y) = Z (x2 + y2 + x)dx + '(y) = x3 + y2x + 3
Далее
Uy0 = Q(x; y) = 2xy + ey + y; 2xy + '0(y) = 2xy + ey + y:
Отсюда
'(y) = (ey + y)dy = ey + y2 :
2
Теперь можно записать общее решение
x3 |
x2 |
y2 |
|
|||
|
+ y2x + |
|
+ ey + |
|
|
= C: |
3 |
2 |
2 |
|
x2 + '(y):
2
1.9Интегрирующий множитель
Если уравнение (??) не является уравнением в полных дифференциалах (условие (1.36) не выполняется), то его можно с использованием интегрирующего множителя преобразовать в уравнение в полных дифференциалах.
Интегрирующим множителем для уравнения (??) называется функция m(x; y) 6= 0 , после умножения на которую оно превращается в уравнение в полных дифференциалах.
То есть для уравнения
mP dx + mQ dy = 0
уже выполняется условие
(mP )y0 = (mQ)x0 : |
(1.39) |
Если функции P (x; y); Q(x; y) имеют непрерывные частные производные и не обращаются в ноль одновременно, то интегрирующий множитель существует.
19
После преобразования (1.39) получается уравнение для нахождения интегрирующего множителя m(x; y)
m(Py0 ¡ Qx0 ) = Qmx0 ¡ P my0 : |
(1.40) |
Общего метода нахождения решения уравнения (1.40) не имеется. В некоторых случаях его решение можно искать путем подбора. Например,
m = m(x); m = m(y); m = m(xy); m = m(x + y); m = m |
µy ¶ |
: |
|
||||
|
|
|
|
x |
|
|
|
В частности, если выражение |
|
|
|
|
|
||
|
Py0 ¡ Qx0 |
|
|
|
|
|
(1.41) |
Q |
|
|
|
|
|||
|
|
|
|
|
|
||
зависит только от x , то интегрирующий множитель можно искать в виде |
m = m(x); |
||||||
если выражение |
|
|
|
|
|
||
|
Py0 ¡ Qx0 |
|
|
|
|
|
(1.42) |
P |
|
|
|
|
|||
|
|
|
|
|
|
зависит только от y , то интегрирующий множитель можно искать в виде m = m(y): Замечание. В некоторых случаях общее решение уравнения (??) можно получить
выделяя в нем полные дифференциалы с помощью формул
d(xy) = y dx + x dy;
µy ¶ |
|
y2 |
|||
d |
x |
|
= |
y dx ¡ x dy |
; |
|
|
Пример 1. Решить уравнение
d(yn) = nyn¡1 dy; d(ln y) = |
dy |
; |
|
y |
|||
|
|
x dx + y dy = 12 d(x2 + y2) и т.д.
y dx ¡ (4x2y + x) dy = 0:
Здесь P = y , Q = ¡(4x2y + x) ,
Py0 = 1; Q0x = ¡(8xy + 1)
условие (1.36) не выполнено. |
|
|
|
|
|
|
Выражение (1.41) |
|
|
|
|
|
|
Py0 ¡ Qx0 |
|
1 + (8xy + 1) |
2 |
|||
|
Q |
= |
|
¡(4xy + 1)x |
= ¡ |
x |
не зависит от y . Поэтому, интегрирующий множитель ищем в виде m = m(x) и определяем его из уравнения (1.40):
|
1 dm |
2 |
|
|
1 |
|
||||||
|
|
|
|
|
= ¡ |
|
; m(x) = |
|
: |
|||
|
m dx |
x |
x2 |
|||||||||
Умножая рассматриваемое уравнение на |
1 |
, получаем уравнение в полных диффе- |
||||||||||
x2 |
||||||||||||
ренциалах |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
xy2 dx ¡ (4y + x1 ) dy = 0:
Его решение.
xy + 2y2 = C:
20