Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 1. Учебное пособие
.pdfФЕДЕРАЛЬНОЕ АГЕНТСТВО СВЯЗИ
Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования
Московский технический университет связи и информатики
А.Г. Кюркчан, Н.И. Смирнова
КОНСПЕКТ ЛЕКЦИЙ ПО ОБЫКНОВЕННЫМ ДИФФЕРЕНЦИАЛЬНЫМ УРАВНЕНИЯМ. Ч.1.
Учебное пособие для студентов 2 курса
направления 01.03.04 – «Прикладная математика»
Москва 2015
1
УДК 517.9
Кюркчан А.Г., Смирнова Н.И. Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 1.: Учебное пособие/МТУСИ. – М., 2015 – 48 с.
Предлагаемое учебное пособие предназначено, главным образом, для студентов,
обучающихся по направлению 01.03.04 «Прикладная математика», но может быть рекомендовано и студентам других направлений. Оно составлено на основе существующей литературы по обыкновенным дифференциальным уравнениям, а
также - методических разработок авторов. Для изучения этой дисциплины необходимы знания математического анализа и линейной алгебры, необходимо также уметь выполнять действия с комплексными числами.
Издание 2-е переработанное и дополненное
Библиогр. 7 назв.
Рецензенты: д.ф.-м.н., профессор В.Г.Данилов, МТУСИ,
к.т.н., профессор А.П.Анютин, РосНОУ.
2
Содержание
1.Обыкновенные дифференциальные уравнения. …………………………….4
1.1Общие понятия. Примеры……………………………………………....4
1.2Уравнения первого порядка, разрешенные относительно производной……………………………………………..……………….5
1.2.1 Уравнения 1-ого порядка, интегрируемые в квадратурах…..…12
2.Линейные дифференциальные уравнения 1-го порядка…………………...14
2.1Метод вариации произвольной постоянной………………………….15
2.2Метод подстановки (Даламбера)………………………………...……17
2.3Уравнения, сводящиеся к линейным…………………………….……18
2.4Уравнения в полных дифференциалах………………………….……..21
2.5Интегрирующий множитель……………………………………….…..24
3.Нормальные системы дифференциальных уравнений………………….….25
3.1Условие Липшица……………………………………………………....26
3.2Теорема существования и единственности для нормальных систем дифференциальных уравнений……………………….……………..…27
4.Теорема о дифференцируемости решений дифференциальных
уравнений………………………………………………………………...…...34
4.1Особые точки…………………………………………………………...35
5.Теорема существования и единственности для дифференциальных уравнений, не разрешенных относительно производной………..………..42
5.1Особые решения………………………………………………………..43
5.2Огибающая семейства кривых………………………………………...46
Список литературы………………………………………………………………48
3
1.Обыкновенные дифференциальные уравнения.
1.1Общие понятия. Примеры.
Обыкновенными дифференциальными уравнениями (о.д.у.) называются
уравнения вида:
|
|
(n) |
(x) 0 , |
(1) |
F x, y(x), y (x),...y |
|
|||
где F - известная функция, x – |
независимая переменная, |
y(x) – неизвестная |
||
функция. |
|
|
|
|
Порядком дифференциального уравнения (д.у.) называется наивысший
порядок производной неизвестной функции y(x) , входящей в уравнение.
Функция y(x) называется решением (или интегралом) д.у. (1), если она n
раз непрерывно дифференцируема на некотором интервале [a,b] и при x [a,b]
удовлетворяет уравнению. Процесс нахождения решения д.у. называется его
интегрированием .
Пример 1. Пусть f (x) – непрерывная на интервале [a,b] функция, |
y(x) – еѐ |
первообразная, тогда |
|
|
(2) |
y (x) f (x) , |
т.е. для отыскания первообразной мы имеем о.д.у. 1-ого порядка. Его решения имеют вид:
x
y(x) f (t)dt C , [x, x0 ] [a,b] ,
x0
где C произвольная постоянная.
Д.у. (2) имеет бесчисленное множество решений и это справедливо для всех о.д.у. Чтобы выделить единственное решение уравнения (2), достаточно задать значение первообразной y(x) в какой-либо точке, например, y(x0 ) y0 . Тогда решение единственно и равно
x
y(x) y0 f (t)dt .
x0
4
Все элементарные функции являются решениями о.д.у.
Пример 2. Уравнение:
|
|
y y 0 |
|
(3) |
|
|
имеет решения y cos x, |
y sin x , что легко |
проверить, |
подставив их |
в |
это |
|
уравнение. Функция y sin x, |
удовлетворяет условиям |
|
1, |
а |
||
y(0) 0, y (0) |
||||||
y cosx, удовлетворяет условиям |
|
0. |
|
|
|
|
y(0) 1, y (0) |
|
|
|
|||
К о.д.у. приводят многие задачи естествознания.
Пример 3. Движение материальной точки массы m описывается II законом
Ньютона ma F . Если точка движется по оси x |
и x(t) - еѐ положение в момент |
|||||||||
времени t , то функция x(t) удовлетворяет о.д.у. II-ого порядка: |
||||||||||
|
|
d 2 x |
|
|
|
dx |
|
|||
m |
|
F t, x, |
|
|
. |
(4) |
||||
dt2 |
|
|
||||||||
|
|
|
|
|
|
dt |
|
|||
Чтобы найти единственное решение уравнения (4), необходимо задать |
||||||||||
начальные условия |
|
|
|
|
|
|
|
|
|
|
x(t |
|
) x , |
dx |
(t |
|
|
) v . |
(5) |
||
0 |
|
0 |
||||||||
|
|
0 |
dt |
0 |
|
|||||
|
|
|
|
|
|
|
|
|
||
Если задача об отыскании всех решений о.д.у. сводится к вычислению конечного числа интегралов и производных от известных функций и к алгебраическим операциям, то говорят, что уравнение интегрируется в квадратурах. Класс этих уравнений узок, поэтому для исследования д.у. как правило применяются приближенные и численные методы.
В нашем курсе мы не будем подробно останавливаться на таких методах.
Ниже лишь коротко рассмотрим на конкретных примерах два метода нахождения приближенного решения о.д.у., разрешенных относительно производной.
1.2Уравнения первого порядка, разрешенные относительно производной.
Уравнение, относящееся к этому классу имеет вид:
5
dy |
f (x, y) |
(6) |
|
dx |
|||
|
|
В дальнейшем будет доказано, что при некоторых ограничениях, налагаемых на функцию f (x, y) , уравнение (6) имеет единственное решение, удовлетворяющее условию:
y(x0 ) y0 .
Общее решение уравнения (6), содержащее все без исключения решения, зависит от одной произвольной постоянной.
Дифференциальное уравнение dydx f (x, y) устанавливает зависимость
между координатами точки и угловым коэффициентом касательной dy
dx к
графику решения в этой точке. Следовательно, дифференциальное уравнение рассматриваемого вида определяет поле направлений и задача интегрирования д.у.
заключается в том, чтобы найти кривые, называемые интегральными кривыми1,
направление касательных к которым в каждой точке совпадает с направлением поля. Поле направлений изображают единичными векторами, образующими с
осью Ox угол arctg dydx , а касательные к интегральным кривым в каждой их точке
должны совпадать по направлению с этими единичными векторами.
На этой идее основан первый (графический) метод построения приближенного решения д.у. вида (6), называемый методом изоклин. Он бывает полезен в тех случаях, когда необходимо «прикинуть» вид интегральных кривых д.у., неинтегрируемого в квадратурах.
Изоклиной называется геометрическое место точек, в которых касательные к искомым интегральным кривым имеют одинаковое направление. Т.е. в каждой точке изоклины векторы поля направлений имеют одинаковый наклон. Таким
1) Пусть y (x) – это решение уравнения (6) на отрезке [a,b] оси х. График этой функции, т.е. кривая y (x) , называется интегральной кривой уравнения (6).
6
образом нахождение изоклин позволяет графически изобразить поле |
||||||||
направлений, задаваемое д.у., а следовательно, и интегральные кривые. |
||||||||
Рассмотрим применение метода изоклин на конкретных примерах. |
|
|||||||
|
Пример 4. Методом изоклин найдем решение следующего д.у.: |
|
||||||
|
|
|
dy |
x2 y2 . |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
Определим уравнение изоклин, считая, что |
|
|
|
|
||||
|
|
|
dy |
k , |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
где k – постоянная. Таким образом, получим |
|
|
|
|
||||
|
|
x2 y2 k, |
x2 y2 k 2 , |
|
|
|
||
т.е. в данном случае изоклинами являются окружности с центром в начале |
||||||||
координат, причем угловой коэффициент касательной к искомым интегральным |
||||||||
кривым равен радиусу этих окружностей. Построим изоклины, присвоив k |
||||||||
некоторые определенные значения (например, k=0.5, k=1, k=2,…). Для построения |
||||||||
поля направлений на каждой изоклине, изобразим единичные векторы, |
||||||||
расположенные под углом, равным арктангенсу радиуса соответствующей |
||||||||
окружности. После этого можно приближенно провести искомые интегральные |
||||||||
кривые. |
|
|
|
|
|
|
|
|
|
Пример 5. Методом изоклин найдем решение следующего д.у.: |
|
||||||
|
|
|
y 1 xy . |
|
|
|
|
|
Здесь изоклинами являются гиперболы k 1 xy или |
|
y |
|
|||||
xy k 1, причем при k=1 гипербола распадается на |
y 0 |
|
|
|||||
|
|
y >0 |
||||||
пару прямых x=0 и y=0. При k=0 получаем изоклину |
|
|
|
|||||
1 xy 0 ; эта гипербола разбивает |
плоскость на |
|
0 |
x |
||||
части, в каждой из которых y сохраняет постоянный |
|
|
|
|||||
знак. |
Интегральные |
кривые |
y y(x), |
пересекая |
|
|
|
|
|
|
|
|
7 |
|
|
|
|
гиперболу 1 xy 0, переходят из области возрастания функции y(x) в область еѐ |
|||||||||||||||
убывания или наоборот, и, следовательно, на ветвях этой гиперболы расположены |
|||||||||||||||
точки максимума и минимума интегральных кривых. Определим теперь знаки |
|||||||||||||||
второй производной в различных областях плоскости: |
|
|
|
|
|
|
|||||||||
|
|
y xy y или y x(1 xy) y x (x2 |
1) y . |
|
|
|
|
||||||||
Кривая |
x (x2 1) y 0 или |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
y x2 |
x |
|
|
|
|
|
|
(7) |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|||
разбивает плоскость на две части, |
в одной из которых y <0, и, следовательно, |
||||||||||||||
|
|
|
|
|
|
интегральные кривые выпуклы вверх, а в |
|||||||||
|
|
y |
|
|
|
другой |
y >0, |
и |
значит, |
интегральные |
|||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y >0 |
|
кривые выпуклы вниз. При переходе через |
||||||||||
|
|
|
|
|
|
||||||||||
y 0 |
0 |
|
|
|
x кривую (7) интегральные кривые переходят |
||||||||||
|
|
|
|
||||||||||||
|
|
|
|
|
|
от |
выпуклости |
к |
вогнутости, |
и, |
|||||
следовательно, на этой кривой расположены точки перегиба интегральных |
|||||||||||||||
кривых. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В результате проведенного исследования известны области возрастания и |
|||||||||||||||
убывания интегральных кривых, известно расположение точек максимума и |
|||||||||||||||
минимума, известны области выпуклости и вогнутости и расположение точек |
|||||||||||||||
перегиба. Этих сведений вполне достаточно для того, чтобы сделать набросок |
|||||||||||||||
семейства интегральных кривых. |
|
|
|
|
|
|
|
|
|
|
|
||||
Второй метод нахождения приближенного решения о.д.у основывается на |
|||||||||||||||
теории возмущений и подходит для решения нелинейных д.у., а также систем д.у, |
|||||||||||||||
разрешенных относительно производной. Он состоит в нахождении |
|||||||||||||||
стационарных состояний или точек равновесия неизвестной функции (функций), |
|||||||||||||||
т.е. точек, в которых производная неизвестной функции по времени обращается в |
|||||||||||||||
ноль, и последующем разложении правой части д.у., разрешенного относительно |
|||||||||||||||
производной, |
в |
ряд |
Тейлора |
в |
окрестности |
этих |
точек. |
|
Это |
позволяет |
|||||
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
аппроксимировать исходное нелинейное д.у. линейным д.у., которое уже можно решить аналитически. Рассмотрим этот метод на примере решения часто встречающейся на практике нелинейной системы о.д.у.
Пример 6. Модель типа хищник-жертва. Рассмотрим динамику популяции двух видов, взаимодействующих между собой по типу хищник-жертва. При этом предполагается, что жертва может найти достаточно пищи, но при каждой встрече с хищником последний убивает жертву.
Пусть x(t) и y(t) - количество жертв и хищников соответственно в момент времени t . Пусть также коэффициент рождаемости жертв b и коэффициент естественной смертности (т.е. без учета уничтожения хищником) d являются константами, причем b d . Таким образом, в отсутствие хищников популяция жертв будет расти со скоростью ( b d )x . Кроме того, предположим, что число случаев, когда хищник убивает жертву зависит от вероятности их встречи и,
следовательно, пропорционально произведению xy .
Таким образом, популяция жертв удовлетворяет следующему д.у.: dxdt x xy ,
где b d 0 , а 0.
Относительно хищников предположим, что их число при отсутствии жертв по естественным причинам убывает со скоростью y . В результате встреч с жертвами число хищников увеличивается. Таким образом, имеем уравнение:
dydt y xy ,
где 0, а 0.
Итак, модель хищник-жертва сводится к системе 2-х нелинейных о.д.у.
9
dx |
|
x xy , |
||
|
|
|||
|
(a) |
|||
dt |
|
|||
|
dy |
|
y xy , |
|
|
|
|
||
dt |
|
|
||
где 0, 0, 0, 0. Эти уравнения называют уравнениями Лотки-Вольтерра.
Для однозначной разрешимости систему (а) необходимо дополнить начальными условиями
|
x(0) x0 , |
y(0) y0 . |
(b) |
||||
Явное аналитическое решение задачи (a), (b) неизвестно. Найдем |
|||||||
приближенное решение этой задачи, основываясь на теории возмущений. |
|
||||||
Первый шаг, как уже отмечалось, состоит в нахождении точек равновесия |
|||||||
( xs , ys ). Эти точки находятся из условий: |
|
||||||
dx |
|
|
|
xs ( ys ) 0, |
|
||
|
|
|
|
||||
|
|
|
|
|
|
||
|
|
||||||
dt |
|
|
|
xs , ys |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
dy |
|
|
|
ys |
( xs ) 0. |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
||
dt |
|
xs , ys |
|
|
|||
|
|
|
|||||
откуда |
|
|
|
|
|
|
|
xs / , |
ys / . |
(с) |
|||||
Кроме того, точкой равновесия будет так называемое тривиальное решение или
точка покоя
xs ys 0 .
Разложим правые части уравнений (а) в ряды Тейлора в окрестности точки ( xs , ys ). С учетом (с) получим
x( y) ( ys )(x xs ) xs ( y ys ) ... xs ( y ys ) ...
y( x) ( xs )(y ys ) ys (x xs ) ... ys (x xs ) ....
Таким образом, в окрестности этой точки исходные уравнения (а) можно
аппроксимировать линейными
10
