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

Диф-уравнения-1 порядка

.pdf
Скачиваний:
9
Добавлен:
11.05.2015
Размер:
293.5 Кб
Скачать

Задание 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

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

Задание 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

IV. Уравнения Бернулли.

Для решения неоднородного уравнения применим метод вариа-

öèè

произвольной

постоянной,

 

считая, что

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