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

Дифференциальные уравнения 1

.pdf
Скачиваний:
47
Добавлен:
18.02.2016
Размер:
537.74 Кб
Скачать

дифференциального уравнения (12.2). На практике его можно найти, используя какие-либо из особенностей задания функций P и Q. Иногда его ищут в виде µ = µ(x) или в виде µ = µ(y). Тогда дифференциальное уравнение (12.2) для искомой функции µ упрощается.

Пример 12.1. Решите дифференциальное уравнение

2(x − y4)dy = ydx.

Решение. В данном случае P = y, Q = 2(y4 − x), ∂P∂y = 1 ≠ ∂Q∂x = 1. Нетрудно видеть, что y = 0 является решением нашего дифференциального уравнения. Пусть теперь y ≠ 0. Будем искать интегрирующий множитель в виде µ = µ(y). Из (12.1) имеем, что

−y∂µ∂y = 3µ.

Решением данного дифференциального уравнения служит, например, функция

1

µ = y3 .

После умножения исходного дифференциального уравнения

1

на y3 получаем уравнение в полных дифференциалах

()

1

 

dx + 2 y

x

dy = 0.

 

 

 

 

 

 

 

y2

y3

Отсюда

 

 

 

 

 

 

 

 

 

 

∂u

1

 

∂u

 

 

 

x

 

 

 

=

 

,

 

= 2(y −

 

).

 

∂x

y2

∂y

y3

Из первого уравнения имеем, что

x

u(x, y) = y2 + φ(y).

41

Подставляя u(x, y) во второе уравнение, получаем, что

φ(y) = 2y,

откуда

φ(y) = y2 − C.

А поэтому все решения исходного дифференциального уравнения задаются формулами

yx2 + y2 = C, y = 0.

В некоторых случаях для решения дифференциальных уравнений (11.1) можно применять метод выделения полных диф-

ференциалов, используя известные формулы:

( )

d(xy) = ydx + xdy, d(yn) = nyn−1dy, d

y

=

xdy − ydx

,

x

x2

 

 

 

и т.д.

Пример 12.2. Решите дифференциальное уравнение

ydx − (4x2y + x)dy = 0.

Решение. Нетрудно видеть, что x = 0 является решением данного дифференциального уравнения. Пусть x ≠ 0. Деля исходное дифференциальное уравнение на −x2, имеем, что

 

y

 

y

 

y

 

d(

 

)

+ 4ydy = 0, d(

 

)

+ d(2y2) = 0, d(

 

+ 2y2) = 0.

x

x

x

Последнее дифференциальное уравнение есть уравнение в полных дифференциалах. Оно имеет общее решение

xy + 2y2 = C.

42

Принимая во внимание случай x = 0, приходим к выводу, что все решения исходного дифференциального уравнения задаются формулами

xy + 2y2 = C, x = 0.

Если в дифференциальном уравнении (11.1) можно выделить полный дифференциал некоторой функции φ(x, y), то иногда данное дифференциальное уравнение упрощается, если от переменных (x, y) перейти к переменным (x, z) или

(z, y), где z = φ(x, y).

Пример 12.3. Решите дифференциальное уравнение

ydx − (x3y + x)dy = 0.

Решение. Нетрудно видеть, что x = 0 является решением данного дифференциального уравнения. Пусть x ≠ 0. Деля

исходное дифференциальное уравнение на −x2, имеем, что

( )

d

y

+ xydy = 0.

x

 

 

Переходя к переменным z = xy и y, получаем дифференциальное уравнение (У–10)

dz + y2 dy = 0, z

которое нетрудно решить (как являющееся дифференциальным уравнением с разделяющимися переменными) (У–11).

43

§ 13. Дифференциальные уравнения первого порядка, не разрешенные относительно производной.

Общий вид дифференциального уравнения первого порядка, не разрешенного относительно производной, имеет вид

F (x, y, y) = 0,

(13.1)

где F есть заданная непрерывная функция в некоторой области G R3с прямоугольными декартовыми координатами (x, y, p), x есть аргумент, y = y(x) есть неизвестная функция.

Определение 13.1. Вектор-функцию x = φ(t), y = ψ(t), t I, где φ и ψ есть непрерывно дифференцируемые на I функции, причем φ(t) ≠ 0, t I, будем называть параметрическим решением дифференциального уравнения

(13.1), если

F

(φ(t), ψ(t), φ(t))

0, t I.

 

 

ψ (t)

 

Из этого решения при x = t вытекает определение явного решения y = ψ(x), x I, дифференциального уравнения (13.1).

Интегральной кривой дифференциального уравнения (13.1) будем называть кривую на плоскости (x, y), описываемую кривой x = φ(t), y = ψ(t), t I. Из определения параметрического решения вытекает, что интегральная кривая дифференциального уравнения (13.1) является гладкой кривой.

Для решения дифференциального уравнения (13.1) можно попытаться разрешить его относительно y. В этом случае получается одно или несколько уравнений вида y= f(x, y), к которым можно попытаться применить исследованные нами методы интегрирования.

44

(13.2)

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

Положим y= p и рассмотрим смешанную систему уравне-

ний {

F (x, y, p) = 0, dy = pdx.

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

и наоборот. В самом деле, если x

= φ(t), y

= ψ(t), t

I, является параметрическим решением дифференциально-

го уравнения (13.1), то

p = p(t)

ψ(t)(φ(t))1(t)

. Отсюда

 

x

 

dy(t) ≡ p(t)dx(t) и F (φ(t), ψ(t), p(t)) 0 на I, т.е. функции φ(t), ψ(t), p(t), удовлетворяют системе (13.2) при всех t I. Обратно, если φ(t), ψ(t), p(t), удовлетворяют при t I системе (13.2), то x = φ(t), y = ψ(t) при t I определяют параметрическое решение дифференциального уравнения (13.1), так как из второго уравнения системы (13.2) находим,

что p(t) = ψ(t), а из первого уравнения этой системы – что

( φ(t) )

ψ(t)

F φ(t), ψ(t), φ(t) 0 на I.

В дальнейшем будем полагать, что уравнение F (x, y, p) = 0 определяет на R3 такую гладкую поверхность S, для которой известно параметрическое представление. Это означает, что существуют такие непрерывно дифференцируемые функции x = φ(u, v), y = ψ(u, v), p = χ(u, v), что на некоторой обла-

45

сти Ω R2 сумма квадратов якобианов положительна:

 

D(x, y)

2

D(y, p)

2

 

D(p, x)

2

 

 

 

 

 

(

 

)

+ (

 

)

 

+ (

 

) > 0

D(u, v)

D(u, v)

 

D(u, v)

и, кроме того,

F (φ(u, v), ψ(u, v), χ(u, v)) 0, (u, v) Ω.

Далее в системе (13.2) перейдем от переменных x, y, p, к переменным u, v по заданным формулам. Первое уравнение системы (13.2) удовлетворяется тождественно на области Ω. Второе уравнение этой системы дает дифференциальное

уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ψ

 

 

∂ψ

 

 

 

∂φ

∂φ

 

 

 

 

 

 

du +

 

 

dv ≡ χ(u, v)(

 

 

du +

 

dv)

 

 

 

 

∂u

∂v

∂u

∂v

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ψ

 

∂φ

 

∂ψ

 

 

∂φ

 

 

 

 

(

 

 

− χ

 

)du +

(

 

− χ

 

)dv = 0.

(13.3)

 

∂u

∂u

∂v

∂v

Переход от дифференциального уравнения (13.1) к дифференциальному уравнению (13.3) называется общим методом введения параметров. В результате мы получили дифференциальное уравнение (13.3) в симметрической форме, некоторые методы решения которого мы рассмотрели ранее. Если, например, множество решений дифференциального уравнения (13.3) задается формулой v = g(u, C), где C есть произвольная постоянная, то функции x = φ(u, g(u, C)), y = ψ(u, g(u, C)), p = χ(u, g(u, C)), удовлетворяют системе (13.2). Поэтому в силу эквивалентности дифференциального уравнения (13.1) и системы (13.2), функции x = φ(u, g(u, C)), y = ψ(u, g(u, C)), задают параметрические решения дифференциального уравнения (13.1).

46

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

Пусть y = f(x, y) есть решение дифференциального уравнения (13.1) относительно y. В этом случае система (13.2)

принимает вид

{

 

y = f(x, p), dy = pdx.

Полагая, что u = x, v = p, получаем, что дифференциальное уравнение (13.3) имеет вид

()

∂f

− p

 

∂f

 

 

dx +

 

dp = 0.

(13.4)

∂x

∂p

Если p = ω(x, C), где C есть произвольная постоянная, задает решения дифференциального уравнения (13.4), то y = f(x, ω(x, C)) есть решения дифференциального уравнения y = f(x, y).

Пусть теперь x = g(y, y) есть решение дифференциального уравнения (13.1) относительно x. В этом случае система (13.2)

принимает вид

{

 

x = g(y, p), dy = pdx.

Полагая, что u = y, v = p, имеем, что дифференциальное уравнение (13.3) имеет вид

()

 

∂g

1

 

∂g

 

p

 

dy + p

 

dp = 0.

(13.5)

∂y

∂p

Если p = ω(y, C), где C есть произвольная постоянная, задает решения дифференциального уравнения (13.5), то формула x = g(y, ω(y, C)) задает решения дифференциального уравнения x = g(y, y).

47

Переход от дифференциального уравнения y = f(x, y) к дифференциальному уравнению (13.4) и переход от дифференциального уравнения x = g(y, y) к дифференциальному уравнению (13.5) называют методом введения параметра.

48

§ 14. Уравнения Лагранжа и Клеро.

Метод введения параметра в общем случае приводит к дифференциальным уравнениям (13.3), не интегрируемым в квадратурах. Приведем два класса дифференциальных уравнений (13.1), для которых метод введения параметра всегда приводит к дифференциальным уравнениям, интегрируемым в квадратурах.

Уравнением Лагранжа будем называть дифференциальное уравнение

y = a(y)x + b(y),

(14.1)

где a(p) ̸≡p и b(p) есть заданные непрерывно дифференцируемые функции по своему аргументу на некотором промежутке.

Положим y= p. От уравнения Лагранжа перейдем к эквивалентной системе{

y = a(p)x + b(p), dy = pdx.

Найдя dy из первого уравнения системы и подставив его во второе уравнение, получаем дифференциальное уравнение xa(p)dp + a(p)dx + b(p)dp = pdx

или

 

(a(p)x + b(p))dp = (p − a(p))dx.

(14.2)

Так как a(p) ̸≡p, то отсюда получаем линейное дифференциальное уравнение первого порядка относительно x:

dx

=

a(p)

x +

b(p)

,

dp

p − a(p)

p − a(p)

 

 

 

которое всегда интегрируется в квадратурах. Если существует p0, такое, что p0 = a(p0), то p = p0 есть решение дифференциального уравнения (14.2). Поэтому

y = a(p0)x + b(p0)

49

8y = 4x2 + 4xp − p2, dy = pdx.

есть частное решение уравнения Лагранжа.

Уравнением Клеро будем называть дифференциальное уравнение

y = xy+ b(y),

(14.3)

где b(p) есть заданная непрерывно дифференцируемая функция по своему аргументу на некотором промежутке.

Полагая, как и в предыдущем случае, y= p, получаем дифференциальное уравнение (У–12)

(x + b(p))dp = 0.

(14.4)

Отсюда приходим к выводу, что либо dp = 0, либо x = −b(p). Если dp = 0, то p = C и решениями уравнения Клеро (14.3)

будут

y = Cx + b(C).

Если же x = −b(p), то равенства

x = −b(p), y = −pb(p) + b(p)

определяют параметрическое решение уравнения Клеро. Замечание 14.1. При решении дифференциального уравнения (14.1) возможно появление составных решений.

Пример 14.1. Решите дифференциальное уравнение

y(4x − y) = 4(2y − x2).

Решение. Разрешив уравнение относительно y, перей–

дем к эквивалентной системе, положив y= p (У–13)

{

Найдя из первого уравнения dy и подставив его во второе уравнение, получаем

2(2x − p)(2dx + p) = 0.

50