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

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

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

Введение.

Обыкновенными дифференциальными уравнениями называют уравнения, связывающие переменную x, неизвестную функцию этой переменной y(x) и ее производные: y0(x); y00(x); y000(x) и так далее. Ðåøå-

íèåì обыкновенного дифференциального уравнения является функция, которая, будучи подставленной, в уравнение превращает его в тождество.

Таким образом, дифференциальные уравнения это уравнения нового типа. Если в алгебраических уравнениях неизвестным являлось некоторое число (для системы уравнений набор чисел), то в дифференциальном уравнении неизвестна уже функция.

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

В качестве примеров приведем несколько дифференциальных уравнений, выражающие законы.

Пример 1. Второй закон Ньютона. Этот закон математически выражается в виде дифференциального уравнения:

m

d2x(t)

= F t; x(t);

dx(t)

;

dt2

 

dt

 

d2x(t)

здесь учтено, что сила F , вызывающая ускорение a = dt2 тела массы m, в общем случае зависит от положения тела (если тело движется в

каком-либо поле), его скорости (учет сопротивления движению), а также может меняться во времени.

Пример 2. Уравнение радиоактивного распада. Закон радиоактивного распада хорошо известен: скорость распада (количество ато-

dx(t)

мов радиоактивного элемента, распавшихся за единицу времени) dt

1

пропорциональна количеству радиоактивного вещества x(t). Этот закон формулируется математически с помощью дифференциального уравне-

dx(t)

íèÿ: dt = kx(t).

Решение этого дифференциального уравнения находится достаточно просто x(t) = e kt. В этом легко убедиться непосредственной подстанов-

кой функции в уравнение.

Пример 3. Описание демографического процесса . Из статисти- ческих данных известно, что в исследуемом регионе число новорожденных и число умерших за единицу времени пропорционально численности жителей региона с коэффициентами k1 è k2 соответственно.

Пусть y(t) число жителей региона в момент времени t. Тогда при-

рост населения за промежуток времени равен разности между числом

родившихся и умерших за это время, т.е. y = ky t, где k = k1 k2 èëèyt = ky. Переходя к пределу при t ! 0, получим уравнение dydt = ky. Его решение имеет вид y = Cekt, где C постоянная, определяемая

начальными условиями (числом жителей региона в начальный момент времени). В этом легко убедиться непосредственной подстановкой функции в уравнение.

Пример 4. Зависимость выручки от реализации товара U(p) от це-

ны p товара задается законом d2U(2p) = (4p2 3)U(p). Решением этого dp

уравнения является функция U(p) = pe 2p2 .

Пример 5. Найдите уравнение кривой, зная, что отрезок, который касательная к кривой в произвольной точке отсекает на оси ординат, равен удвоенной ординате точки касания.

Пусть уравнение кривой y = f(x). Тогда уравнение касательной к этой кривой в точке M с координатами (x; y) имеет вид Y y = y0(x)(X x). Отрезок, отсекаемый кривой на оси Y равен Y = y y0(x) (так как X = 0). По условию Y = 2y. Получаем уравнение y xy0 = 2y èëè xy0 + y = 0.

Умножив обе части этого уравнения на dx, получим xdy + ydx = 0 или d(xy) = 0. Откуда xy = C или y = Cx .

2

Дифференциальное уравнение задает семейство гипербол, асимптотами которых являются оси координат. Для выделения одной кривой необходимо задать на кривой точку. Пусть точка (2; 3) лежит на кривой.

Значит xy = 6, т.е. C = 6.

y

6

3

-

2

x

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

1. Дифференциальные уравнения первого порядка.

Обыкновенным дифференциальным уравнением 1-го порядка

называется уравнение вида

F (x; y; y0) = 0

(1:1)

Разрешив его относительно производной, получим обыкновенное дифференциальное уравнение 1-го порядка, разрешенное относительно производной

y0 = f(x; y)

(1:2)

Уравнение первого порядка можно записать в дифференциальной форме

M(x; y)dx + N(x; y)dy = 0

(1:3)

3

В случае, если f(x; y) = f(x), то есть функция зависит только от одной переменной x, получим простейшее дифференциальное уравнения 1-го порядка y0 = f(x), решение которого y = R f(x)dx + C легко найти.

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

называют интегралом дифференциального уравнения . Дифференциальное уравнение имеет бесконечно много решений (ин-

тегралов). Для обозначения этого множества вводится понятие общего решения (общего интеграла) дифференциального уравнения. Общее решение (общий интеграл) дифференциального уравнения первого порядка содержит независимый параметр C. Наличие в общем решении

констант связано с тем, что при решении дифференциальных уравнений выполняется интегрирование операция, обратная дифференцированию. При интегрировании возникает константа. Между понятиями общее решение и общий интеграл принципиальных различий нет. Просто под термином решение принято понимать обычную функцию, а под термином интеграл функцию, заданную неявно. Так сложилось исторически, и такая терминология принята во всем мире.

Определить значение константы интегрирования можно, если мы потребуем, чтобы решение удовлетворяло дополнительному условию. Задача нахождения решения уравнения, удовлетворяющего условию y(x0) = y0, ãäå x0 è y0 заданные числа называется задачей Коши, а условие y(x0) = y0 начальными условиями.

График какого-либо частного решения дифференциального уравнения первого порядка (при фиксированном значении параметра C = C0) задает на плоскости (x; y) линию, которую называют интегральной кривой. При изменении значения параметра C получится другая интегральная кривая. Перебрав все возможные значения параметра C, мы

получим некоторый набор интегральных кривых, который называют семейством интегральных кривых. Таким образом, общее решение диффе-

4

ренциального уравнения первого порядка (точно так же как и его общий интеграл) задает на плоскости (x; y) семейство интегральных кривых. Из

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

При построении решений дифференциальных уравнений возникает очень важный вопрос: всегда ли можно найти частное решение уравнения, удовлетворяющее заданному начальному условию (решить задачу Коши). Если нам удалось решить задачу Коши, то есть построить какоелибо частное решение, нельзя ли найти еще одно или несколько таких решений. В тех случаях, когда решение построено (или подобрано), актуальной остается лишь вторая часть вопроса. А если сразу не удается найти решение, то важно знать, существует ли искомое решение в принципе (ведь искать то, чего нет, занятие весьма неблагодарное). Ответ на этот очень важный вопрос дает следующая теорема о существовании и единственности решения задачи Коши для уравнения первого порядка.

Теорема 1.1. Если в дифференциальном уравнении y0 = f(x; y)

1. функция f(x; y) непрерывна и ограничена в некоторой области D

плоскости (x; y),

2. имеет в области D ограниченную частную производную: @f@y , èëè, что то же самое, удовлетворяет в области D условию Липшица:

jf(x; y1) f(x; y2)j 6 Njy1 y2j, ãäå N = const > 0, то для любой точки (x0; y0) 2 D в некотором интервале x0 6 x 6 x + 0 + существует и притом единственное решение этого уравнения y = y(x),

удовлетворяющее условию y(x0) = y0.

Геометрический смысл теоремы достаточно прост: через каждую точ- ку области D, в которой задано дифференциальное уравнение, обяза-

тельно проходит, причем только одна, интегральная кривая.

Если в одной или некотором подмножестве точек той области плоскости (x; y), в которой нас интересует решение дифференциального уравнения, нарушается хотя бы одно из условий теоремы 1., то нельзя гаран-

5

тировать наличие частных решений этого уравнения, проходящих через эти точки, а если такие решения существуют, то нельзя гарантировать их единственность.

Те точки плоскости (x; y), через которые в результате нарушения

условий теоремы проходит несколько интегральных кривых, называются особыми точками. А те решения уравнения y0 = f(x; y) (èíòå-

гральные кривые), каждая точка которых является особой называются

особыми решениями (особыми интегральными кривыми ). Отметим, что особое решение дифференциального уравнения не мо-

жет быть получено из общего решения ни при каких значениях константы интегрирования C. Это особое решение (если оно существует для данного уравнения) является особой интегральной кривой.

2. Уравнения с разделяющимися переменными.

Пусть в уравнении

y0 = f(x; y)

(2:1)

функция f(x; y) представима в виде произведения двух функций, каждая из которых зависит только от одной переменной: f(x; y) = f1(x)f2(y): Тогда уравнение (2.1) принимает вид

y0 = f1(x) f2(y);

(2:2)

и его называют уравнением с разделяющимися переменными . Запишем дифференциал функции y(x)

dy

dy = dx dx = f1(x)f2(y)dx;

Ограничимся рассмотрением таких значений y, при которых f2(y) 6= 0. Тогда, разделив обе части полученного уравнения на f2(y) и введя (для

удобства) функцию g(y) =

1

 

 

f2(y), получим уравнение, которое называют

уравнением с разделенными переменными

 

 

g(y)dy = f1(x)dx:

(2:3)

6

Действительно, если в исходном уравнении в правой части были обе переменные (и x и y), то в последнем уравнении переменная x встречается

только в правой части, а переменная y только в левой.

Òàê êàê dy = y0(x)dx, то обе части уравнения с разделенными пе-

ременными представляют собой дифференциалы некоторых функций. Это значит, что функция g(y(x))y0(x) является производной некоторой

функции (первообразной от g(y(x))y0(x)) è f1(x) также есть производная некоторой функции (первообразной от f1(x)). Поскольку равны дифференциалы функций, то с точностью до постоянной равны и сами функции. Следовательно, можно записать следующее равенство:

R g(y(x))y0(x)dx = R f1(x)dx + C:

Перейдя в интеграле правой части этого выражения к интегрированию по y, то есть сделав замену переменных, получим выражение, которое и

определяет зависимость y от x.

Z

Z

 

g(y)dy =

f1(x)dx + C:

(2:4)

Пример 2.1. Решите уравнение y0 = xy2 + 4y. Найдите частное решение, y 2

удовлетворяющее начальным условиям y(3) = 1.

Решение . Данное уравнение является уравнением с разделяющимися переменными. Вынесем в правой части за скобку общий множитель и запишем производную

как отношение дифференциалов

dy

=

y(x + 4)

. Обе части уравнения умножим сна-

dx

y

2

2

2

 

чала на dx, затем на

y

2

2

 

 

 

 

 

 

 

. В уравнении

 

y

2dy = (x+4)dx правая часть зависит

 

 

 

y

 

 

 

 

 

 

y

 

только от x, а левая только от y. Найдем первообразные для левой и правой частей

уравнения:

 

 

dy =

 

 

 

 

 

 

y2

2

 

 

2

y2

 

 

 

 

R

y

 

dy =

x2 y y

2 2 ln jyj;

 

 

 

(x + 4)dx = R

+ 4x.

 

 

 

 

 

 

R

 

 

 

2

 

 

 

y2

 

= x2

 

 

 

 

 

 

 

 

 

 

Общий интеграл уравнения имеет вид

2 ln y

+ 4x + C.

 

 

 

 

 

 

 

 

2

j j

2

 

Чтобы найти частное решение, подставим начальные условия ( x = 3, y = 1) в

 

1

2 ln 1 =

9

+12+C. Следовательно, C = 16 и частный

общий интеграл. Получим 2

2

 

y2

 

x2

 

 

интеграл имеет вид

2 2 ln jyj =

2 + 4x 16.

Задания для самостоятельного решения.

7

Решите уравнение. Если заданы начальные условия, решите задачу Коши.

Задание 2.1. x(1 + y2)dx y(1 + x2)dy = 0.

 

 

 

 

 

 

 

 

Задание 2.2.

y0

p

 

 

2

 

 

 

 

 

 

 

 

2

è

 

 

 

 

 

p

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1 x

 

= (x + 3)y

 

 

 

 

 

y( 3=2) =

.

 

 

 

Задание 2.3. y0

(x2 3x + 2) ln2 y = 3xy 5y.

 

 

 

 

 

 

 

 

 

y2 sin3 x = cos x

 

 

 

 

 

 

 

 

 

è y( =6) = p3

 

.

 

 

 

 

 

 

 

 

y3

 

 

 

 

 

Задание 2.4. y

 

 

 

1

2

 

 

 

Задание 2.5.

 

0

 

 

 

 

x x2

p

 

 

 

 

 

 

 

2

+ 1)

è

 

 

 

 

.

 

 

y0

(2y + 1)e

= (2x

 

 

1)(y

 

 

 

y(1) = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание 2.6. y0

(2y 3)ex x2

= (2x 1)(y2 3y + 1) è y(1) = 3.

Задание 2.7. y0

(x 4)2 tg y = (x + 1) cos4 y.

 

 

 

 

 

 

 

 

 

 

y0 y2

 

 

y + 1) 4

 

2x

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание 2.8.

 

 

(3

2 x 1 p

 

 

 

 

 

 

 

 

 

 

=

 

 

3 y3 y2 + y 5 è y(1) = 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание 2.9.

y0

 

 

 

 

2

 

 

 

3

 

2

 

 

 

è

 

 

 

 

 

p3 .

 

 

 

 

 

 

 

(2xy + xy

) = y

 

ln

 

 

x

 

 

 

 

 

y(e ) = 1

 

 

 

 

 

 

Задание 2.10. y0(2y + 1)ex2+y = x è y(1) =

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание 2.11.

y0(2y

 

 

 

 

 

 

 

 

 

 

y2

 

3y

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

3) 2x + 7 = e

 

 

 

 

(x + 2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание 2.12. y0(2y + 5)(x2 + 2x

 

 

3) = (x

 

5)(y2 + 5y

 

4) è y(3) = 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Однородные уравнения.

Некоторые дифференциальные уравнения путем замены переменных или введением некоторой специальной комбинации переменных можно привести к виду уравнений с разделяющимися переменными. Одним из таких уравнений являются однородные уравнения.

Функция f(x; y) называется однородной измерения k, если для любого действительного числа t выполняется условие f(tx; ty) = tkf(x; y). Функция f(x; y) называется однородной измерения 0, если при умножении x и y на один и тот же множитель функция не меняется, т.е.

f(tx; ty) = f(x; y). Однородную функцию нулевого измерения можно представить в виде f(x; y) = ' xy .

8

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

2

 

 

 

 

3

 

 

 

Пример 3.1. Проверьте, что функция f(x; y) =

3x

x

 

y2+ 6y2

 

ÿâëÿ-

3

 

 

ется однородной нулевого измерения.

 

 

4x

5xy

+ x

y

 

 

 

 

 

 

 

 

 

 

tx

3

2

ty

ty 3

 

t3

x3

 

x2y

+ 6

y3

 

Решение . Имеем f(tx; ty) =

3(

)

 

(tx) (

2) + 6( 2)

=

(3

 

 

 

 

 

 

2 )

=

4(tx)

3

 

3

 

3

 

 

2

 

 

 

3

2

 

3

 

5(tx)(ty)

+ (tx) (ty)

t (4x

 

5xy

 

+ x y)

 

=

3x

x

y2

+ 6y2

 

= f(x; y).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4x

5xy

+ x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Однородные уравнения первого порядка записываются в виде

y0 = f(x; y)

(3:1)

где f(x; y) однородная функция нулевого измерения. Сделав замену

 

y(x)

z(x) =

x или, что тоже самое y(x) = x z(x), вычислив производную

y0 = z(x) + x z0(x) и обозначив f(x; y) = '(z), уравнение (3.1) можно

записать в виде z(x) + x

 

z0(x) = '(z), è ïðè '(z)

 

z = 0 получим

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

уравнение с разделяющимися переменными

 

 

 

 

 

 

 

 

 

 

 

dz

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

:

 

 

 

 

 

 

 

 

 

 

'(z) z

x

 

 

 

 

 

 

 

Интегрирование обеих частей уравнения дает

 

 

 

 

 

 

 

 

Z

'(z) z = ln jxj ln C:

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(здесь, для удобства, константа интегрирования записана в виде

ln C).

Из этого выражения можно получить функцию x = x(z): x = CeR

 

dz

 

'(z) z .

Пример 3.2. Решите уравнение

 

dy

y

2 y

 

 

 

 

 

 

 

dx = x + sin

x.

 

 

 

 

 

 

Решение . Сделаем замену z(x) =

 

y(x)

, откуда получаем, что y(x) = x z(x).

 

 

x

 

Вычислим производную y0

= z(x) + x

z0(x). Подставим y0 и y в уравнение и получим

2

dz

 

 

2

 

 

 

 

 

 

 

 

dz

dx

 

 

 

z + xz0 = z + sin z èëè xdx

= sin z. Разделим переменные

 

= x и найдем

sin2 z

первообразные для обеих частей уравнения: ctg z = ln jxj C. Возратившись к

y

исходным переменным, получим общий интеграл уравнения ctg x = C ln jxj.

dy

4x3 + 2x2y + y3

 

Пример 3.3. Решите уравнение dx =

x3 + xy2 .

 

Решение . Легко показать, что функция f(x; y) =

4x3 + 2x2y + y3

является одно-

 

 

 

x3 + xy2

 

 

 

y(x)

 

родной нулевого измерения. Сделаем замену z(x) =

 

 

èëè y(x) = x z(x). Âû-

 

x

числим производную y0 = z(x) + x z0(x). Подставим y0 и y в уравнение и получим

9

z(x) + x z0(x) =

4x3

+ 2x3z + x3z3

 

 

 

 

4 + 2z + z3

 

 

 

 

x3 + x3z2

 

 

èëè z(x) + x z0

(x) =

 

 

1 + z2

. Перенеся z â

правую часть и приведя к общему знаменателю, получим x z0(x) =

 

4 + z

 

 

1 + z2

. Разделим

переменные

 

(z2 + 1)dz

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= x и найдем первообразные для обеих частей уравнения:

z + 4

 

 

(z2 + 1)dz

 

z

17

 

 

dz =

z2

 

 

 

 

 

 

R

dx

 

 

 

 

 

 

 

=

4 +

 

 

 

 

2 4z + 17 ln jz + 4j

è

 

x

= ln jxj. Тогда

z2

z + 4

z + 4

 

R

 

 

 

R

 

 

 

 

 

 

 

 

. Возратившись к исходным

 

 

 

 

 

 

 

2 4z + 17 ln(z + 4) = ln jxj + C

 

 

 

 

 

 

 

 

переменным, получим

общий интеграл уравнения

y2

 

 

 

4y + 17 ln

j

y + 4xj

= ln

x

+ C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x2

x

 

jxj

j

j

 

 

 

 

 

 

 

 

 

 

К однородным приводятся уравнения, правая часть которых зави-

сит от линейной комбинации переменных: dy = f

 

a1x + b1y + c1

ïðè

a1b2

a2b1

= 0, если положить z = a1x + b1y + c1; t2

= a2x + b2y + c2,

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

a x + b2y + c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выразить y и x через z и t и считать z функцией от t. Можно также ввести новые переменные u = x ; v = y , причем и подбирают

так, чтобы выражение в правой части стало однородным.

Пример 3.4. Решите уравнение

dy

= x + y 3

 

 

 

dx

y x + 1.

du

u + v

Сделав замену x = u 2; y = v

 

 

1, приведите уравнение к виду

dv

= v u è

решите получившееся однородное уравнение.

Задания для самостоятельного решения.

Решите уравнение. Если заданы начальные условия, решите задачу Коши.

Задание 3.1. 2xy0 = y +

y2

 

10

x 4x è y( 1) =

3 .

Задание 3.2. (y x)ydx + x2dy = 0.

 

Задание 3.3. y0 =

x2

2

xy + y2

 

 

 

2 è y(1) = 4.

 

 

xy 3x

 

 

Задание 3.4.

Задание 3.5.

Задание 3.6.

Задание 3.7.

y0

=

2

x2

 

3

xy + y2

 

 

 

 

 

 

2 .

 

 

 

 

 

 

xy 3x

 

 

 

 

y0

=

y + cos y sin 3 y

y(2) =

2

 

 

x

 

 

 

x

 

x è

 

3 .

y0

=

x2

2

xy

+

4y2

 

 

 

 

 

 

2 .

 

 

 

 

 

 

3xy 4x

 

 

 

y0

=

 

 

2

 

 

 

 

 

y2

 

 

3x2

2xy + 52 .

 

 

 

 

 

x

 

+ 2xy y

 

 

10