Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 2. Учебное пособие
.pdfФЕДЕРАЛЬНОЕ АГЕНТСТВО СВЯЗИ
Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования
Московский технический университет связи и информатики
А.Г. Кюркчан, Н.И. Смирнова
КОНСПЕКТ ЛЕКЦИЙ ПО ОБЫКНОВЕННЫМ ДИФФЕРЕНЦИАЛЬНЫМ УРАВНЕНИЯМ. Ч.2.
Учебное пособие для студентов 3 курса
направления 01.03.04 – «Прикладная математика»
Москва 2015
УДК 517.9
Кюркчан А.Г., Смирнова Н.И. Конспект лекций по обыкновенным дифференциальным уравнениям. Ч. 2.: Учебное пособие/МТУСИ. – М., 2015
– 48 с.
Предлагаемое учебное пособие предназначено, главным образом, для студентов, обучающихся по направлению 01.03.04 «Прикладная математика», но может быть рекомендовано и студентам других направлений. Оно составлено на основе существующей литературы по обыкновенным дифференциальным уравнениям, а также - методических разработок авторов. Для изучения этой дисциплины необходимы знания математического анализа и линейной алгебры, необходимо также уметь выполнять действия с комплексными числами.
Издание 2-е переработанное и дополненное
Библиогр. 7 назв.
Рецензенты: д.ф.-м.н., профессор П.К.Суетин, МТУСИ,
д.ф.-м.н., профессор А.П.Анютин, РосНОУ.
2
Содержание |
|
|
1. |
Линейные дифференциальные уравнения n-ого порядка |
4 |
2. |
Линейные неоднородные дифференциальные уравнения |
9 |
2.1. Метод вариации произвольных постоянных |
12 |
|
2.2. Метод Коши решения линейного неоднородного уравнения |
16 |
|
3. |
Линейные однородные дифференциальные уравнения с |
|
|
постоянными коэффициентами |
17 |
4. |
Линейные неоднородные уравнения с постоянными |
|
|
коэффициентами и специальным видом правой части |
22 |
5. |
Системы линейных дифференциальных уравнений |
23 |
5.1. Метод вариации постоянных |
30 |
|
5.2.Системы линейных дифференциальных уравнений с
постоянными коэффициентами |
31 |
6. Матричная экспонента |
39 |
6.1.Матрицы и дифференциальные уравнения. Дифференцирование
и интегрирование векторных и матричных функций |
42 |
6.2.Нормы векторов и матриц. Бесконечные ряды векторов и
матриц |
44 |
6.3.Существование и единственность решения векторного
|
дифференциального уравнения |
45 |
6.4. |
Матричная экспонента |
48 |
6.5. Неоднородное векторное дифференциальное уравнение |
49 |
|
6.6. |
Неоднородное уравнение. Сопряженная система |
50 |
Список литературы………………………………………………………………52
3
Линейные дифференциальные уравнения n-ого порядка.
Уравнение вида
a0 (x) y(n) a1 (x) y(n 1) .. an 1 (x) y an (x) y (x), |
(1) |
линейное относительно неизвестной функции и еѐ производных, называется
линейным дифференциальным уравнением n – ого порядка.
|
Если правая часть (x) 0, |
при a x b, то |
уравнение |
называется |
||||||||||||
линейным однородным. |
|
|
|
|
|
|
|
|
|
|||||||
|
Если a0 (x) 0 при a x b, то на этом отрезке однородное уравнение (1) |
|||||||||||||||
эквивалентно следующему |
|
|
|
|
|
|
|
|||||||||
|
y(n) p1(x) y(n 1) .. pn 1(x) y pn (x) y 0, |
(2) |
|
|
||||||||||||
где |
pi (x) |
ai (x) |
. Уравнение (2) запишем также в виде |
|
|
|
||||||||||
ao |
(x) |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
y(n) (x) pi (x)y(n i) (x). |
|
|
(2a) |
|
|
||||||||||
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
||
|
Если коэффициенты |
pi (x) |
непрерывны на отрезке [a,b], то в |
|||||||||||||
окрестности любых начальных значений |
|
|
|
|
||||||||||||
|
|
|
|
0 |
|
|
|
|
0 |
,...,y |
(n 1) |
0 |
|
|
|
|
|
y(x0 ) y1 |
, y (x0 ) |
y2 |
|
(x0 ) yn , |
|
|
|
|
|||||||
где |
x0 – любая точка интервала a<x<b, удовлетворяется условие теоремы |
|||||||||||||||
существования и единственности (см. лекцию 3). |
|
|
|
|
||||||||||||
|
Действительно, правая часть (2а) непрерывна по совокупности |
|||||||||||||||
аргументов |
|
|
и |
|
|
существуют |
|
ограниченные |
по модулю |
частные |
||||||
производные |
|
f |
|
|
p |
|
(x), k 0,1,...,(n 1) , т.к. |
функции |
p |
(x) |
- |
|||||
|
|
|
|
n k |
||||||||||||
|
|
|
y(k ) |
|
|
|
|
|
|
n k |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
непрерывны на отрезке [a,b] и, следовательно, ограничены по модулю. Таким образом, функция f (x, y, y ,.., y(n 1) ) удовлетворяет условию Липшица по
всем переменным, начиная со второй.
Запишем линейное однородное уравнение
y(n) p (x) y(n 1) .. p (x) y 0, |
|||
1 |
n |
|
|
коротко в виде L[ y] 0 , где |
|
|
|
L[ y] y(n) p (x) y(n 1) |
.. p |
n |
(x) y, - |
1 |
|
|
|
линейный дифференциальный оператор (линейное дифференциальное выражение).
Линейный дифференциальный оператор обладает свойствами:
1. L[cy] cL[ y],
т.к. (cy)(n) p1 (x)(cy)(n 1) .. pn (x)(cy)
c y(n) p1 (x) y(n 1) .. pn (x) y 2. L[ y1 y2 ] L[ y1] L[ y2 ],
4
т.к. ( y y |
2 |
)(n) p (x)(y y |
2 |
)(n 1) |
.. p (x)(y y |
2 |
) |
|
|
|
|
||||
1 |
1 |
1 |
|
n |
1 |
|
|
|
|
|
|
||||
[ y(n) |
p (x) y(n 1) |
.. p (x) y ] [ y(n) p (x) y(n 1) |
.. p |
n |
(x) y |
2 |
]. |
||||||||
1 |
1 |
1 |
n |
|
1 |
2 |
1 |
2 |
|
|
|
|
|||
Из 1) и 2) следует, что |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
m |
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L[ Ci yi ] Ci L[ yi ], |
|
|
|
|
|
|
|
|
|
|||
|
|
|
i 1 |
i 1 |
|
|
|
|
|
|
|
|
|
|
|
где Ci |
– постоянные. |
|
|
|
|
|
|
|
|
|
|
|
|||
Опираясь на свойства линейного оператора L, установим ряд теорем о решениях линейных однородных дифференциальных уравнений (л. о. д. у.).
Теорема 1. |
|
|
|
|
|
|
|
|
Если y1 |
является решением л. о. д. |
у. |
L[ y1 ] 0 , |
то и |
Cy1 , |
где |
С – |
|
произвольная постоянная, является решением того же уравнения. |
|
|
||||||
Доказательство. В соответствии со свойством 1): |
|
|
|
|
||||
|
L[Cy1 ] CL[ y1 ] 0 , ч.т.д. |
|
|
|
|
|
|
|
Теорема 2. |
|
|
|
|
|
|
|
|
Сумма |
y1 y2 решений |
л. о. д. у. |
L[y]=0 является решением |
того |
же |
|||
уравнения. |
|
|
|
|
|
|
|
|
Доказательство. Воспользуемся свойством 2) |
|
|
|
|
||||
|
L[ y1 y2 ] L[ y1] L[ y2 ] 0, ч. т. д. |
|
|
|
|
|||
Следствие. Линейная |
комбинация |
с |
произвольными |
постоянными |
||||
|
m |
|
|
|
|
|
|
|
коэффициентами Ci yi |
решений y1, y2 ,...,ym л. о. |
д. у. |
L[y]=0 |
является |
||||
i 1
решением того же уравнения.
Теорема 3.
Если л. о. д. у. L[y]=0 с действительными коэффициентами pi (x) имеет комплексное решение y(x) U (x) iV (x) , то его действительная U(x) и
мнимая V(x) части в отдельности являются решениями того же уравнения. Доказательство. Используя свойства 1) и 2) имеем:
L[U iV ] L[U ] iL[V ] 0,
откуда L[U] 0, L[V] 0, т.к. комплексная величина тождественно равна нулю, тогда и только тогда, когда ее действительная и мнимая части тожественно равны нулю.
Определение. Функции y1(x), y2 (x),...,yn (x) называются линейно зависимыми (предполагается, что ни одна из функций yi (x) тождественно не
равна нулю на [a,b]) на некотором отрезке a x b если существуют постоянные числа 1, 2 ,.., n такие, что на этом отрезке
1 y1(x) 2 y2 (x) ... n yn (x) 0 , |
(3) |
причем хотя бы одно i 0.
5
Если тождество (3) справедливо лишь при 1 2 ... n 0, то функции y1(x), y2 (x),...,yn (x) называются линейно независимыми на отрезке a x b.
Пример 1.
Функции 1, x, x2 ,...,xn линейно независимы на любом отрезке [a,b], т.к. тождество
1 2 x 3 x2 ... n 1xn 0
возможно лишь, если все i 0 . Если хотя бы одно i 0 , то слева стоял бы
многочлен степени не выше n, который может иметь не более n различных корней и, следовательно, обращаться в нуль не более чем в n точках отрезка.
Пример 2.
Функции eq1x ,eq2 x ,...,eqn x , где qi q j , если i j, линейно независимы на любом
отрезке a x b.
Допустим, что эти функции линейно зависимы, т.е.
|
|
|
eq1x |
eq2 x ... |
eqn x 0, |
x [a,b], |
||||
|
|
|
1 |
2 |
|
|
n |
|
|
|
причем не все i |
равны нулю. Пусть, например, n 0 . Разделим тождество |
|||||||||
на eq1 x , а результат продифференцируем, получим |
||||||||||
|
2 |
(q |
q )e(q2 q1 ) x ... |
n |
(q |
q )e(qn q1 ) x 0 |
||||
|
2 |
1 |
|
|
|
n |
1 |
|
||
Продолжая эту процедуру (т.е. деля на e(q2 q1 ) x и дифференцируя и т.д.) (n-1) раз, получим
|
n |
(q |
q )(q |
q )...(q |
q |
)e(qn qn 1 ) x 0, |
|
|
||||
|
n |
1 |
n |
2 |
n |
|
n 1 |
|
|
|
|
|
что невозможно, |
т.к. |
qi q j |
|
при |
i j, а |
n 0 |
по предположению, т.к. в |
|||||
качестве n |
можно выбрать любой коэффициент, то свойство доказано. |
|||||||||||
Теорема 4. |
|
|
|
|
|
|
|
|
|
|||
Если функции y1(x), y2 (x),...,yn (x) |
линейно зависимы на отрезке a x b , то |
|||||||||||
на этом отрезке определитель Вронского |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
y1 |
y2 . . |
yn |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 |
y2 . . |
yn |
|
|
W (x) W[ y1, y2 ,...,yn ] |
|
|
|
|
|
|
||||||
|
y1 |
y2 . . |
yn |
|
||||||||
|
|
|
|
|
|
|
|
. |
. |
. . |
. |
|
|
|
|
|
|
|
|
|
y(n 1) |
y(n 1) . . |
y(n 1) |
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
n |
|
тождественно равен нулю. |
|
|
|
|
|
|
||||||
Доказательство. Дано, что |
|
|
|
|
||||||||
|
|
|
1 y1 |
2 y2 ... n yn 0 |
|
(4) |
||||||
На отрезке [a,b], примем не все i равны нулю. Дифференцируя тождество
(4) (n-1) раз, получим
6
y |
|
|
2 |
y |
2 |
... |
n |
y |
n |
0, |
|
|
||||
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|||||
|
y |
|
2 |
y |
... |
n |
y |
0, |
|
|
||||||
|
1 |
1 |
|
|
2 |
|
|
|
|
n |
|
|
|
|||
.......................................... |
(5) |
|||||||||||||||
|
|
|
|
|
y |
|
... y |
|
0. |
|||||||
y |
(n 1) |
(n 1) |
(n 1) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
1 |
1 |
|
|
|
|
2 |
2 |
|
|
|
|
n n |
|
|||
Эта линейная алгебраическая система имеет нетривиальное решение (по предположению) при любом x [a,b]. следовательно, еѐ определитель, являющийся определителем Вронского, тождественно равен нулю при x [a,b], ч. т. д.
Теорема 5.
Если линейно независимые функции y1(x), y2 (x),...,yn (x) являются решениями л. о. д. у.
y(n) p (x) y(n 1) |
.. p (x) y 0, |
(6) |
1 |
n |
|
c непрерывными на |
отрезке a x b коэффициентами |
pi (x) , то |
определитель Вронского этой системы функций не может обратиться в нуль ни в одной точке отрезка.
|
y1 |
|
|
|
y1 |
W (x) |
|
y1 |
|
|
. |
|
y(n 1) |
|
1 |
y2
y
2
y
2
.
y2(n 1)
. .
. .
. .
. .
. .
yn
y
n
y
n
.
yn(n 1)
Доказательство. Пусть в некоторой точке x x0 отрезка [a,b] определитель Вронского W (x0 ) 0 . Выберем постоянные i , i 1,2,...,n так,
чтобы удовлетворялась система уравнений. Рассмотрим систему алгебраических уравнений:
y (x |
0 |
) |
2 |
y |
2 |
(x |
0 |
) ... |
n |
y |
n |
(x |
0 |
) 0, |
|
|
|
|||||||||||||
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
1 |
y |
(x |
0 |
) |
2 |
y |
(x |
0 |
) ... |
n |
y |
(x |
0 |
) 0, |
|
|
|
||||||||||||
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
||||||||
.......................................... |
|
|
|
|
|
|
|
|
|
|
(7) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(n 1) |
(x |
|
) |
|
y |
(n 1) |
(x |
|
) ... |
|
|
y |
(n 1) |
(x |
|
) 0 |
|||||||||||
y |
|
|
|
0 |
2 |
2 |
|
0 |
n |
n |
0 |
|||||||||||||||||||
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
и чтобы не все i равнялись нулю. Это всегда возможно, т.к. опре-делитель
системы (7) равен нулю, т.е. нетривиальное решение системы (7). При таком выборе i линейная комбинация
y 1 y1(x) 2 y2 (x) ... n yn (x)
Будет решением л. о. д. у. (6), удовлетворяющим, в силу уравнений системы (7), нулевым начальным условиям
|
(n 1) |
(x0 ) 0, |
(8) |
y(x0 ) 0, y (x0 ) 0,...,y |
|
7
Таким условиям, очевидно, удовлетворяет тривиальное решение y=0 уравнения (6) и по теореме единственности начальным условиям (8)
удовлетворяет |
только |
это |
решение. |
Следовательно, |
1 y1(x) 2 y2 (x) ... n yn (x) 0 |
и решения y1(x), y2 (x),...,yn (x) , вопреки |
|||
условию теоремы, линейно зависимы.
Замечание 1. Из теорем 4 и 5 следует, что линейно независимые на отрезке [a,b] решения y1(x), y2 (x),...,yn (x) уравнения (6) линейно независимы
и на любом отрезке [a1, b1] [a, b].
Замечание 2. В теореме 5, в отличие от теоремы 4 предполагалось, что функции y1(x), y2 (x),...,yn (x) являются решениями л. о. д. у. (6) с
непрерывными коэффициентами. Отказаться от этого требования и считать функции y1(x), y2 (x),...,yn (x) произвольными (n-1) раз непрерывно
дифференцируемыми нельзя. Легко привести примеры линейно независимых функций, не являющихся решениями уравнения (6) с непрерывными коэффициентами, для которых определитель Вронского не только обращается в нуль в некоторых точках, но даже тождественно равен нулю.
Теорема 6.
Общим решением при a x b л. о. д. у.
y(n) p (x) y(n 1) |
.. p (x) y 0, |
|
|
(6) |
|
||||
|
1 |
|
|
n |
|
|
|
|
|
с непрерывными на отрезке [a,b] коэффициентами |
pi (x) , i=1,2,..,n, является |
||||||||
|
|
|
n |
|
|
|
|
|
|
линейная комбинация |
y Ci yi |
из |
n линейно |
независимых на том же |
|||||
|
|
|
i 1 |
|
|
|
|
|
|
отрезке частных |
решений |
yi , i 1,2,...,n |
с произвольными постоянными |
||||||
коэффициентами Ci . |
|
|
|
|
|
|
|
|
|
Доказательство. |
Уравнение |
(6) |
при |
x [a,b] |
удовлетворяет условиям |
||||
|
|
|
|
|
|
|
|
n |
|
теоремы существования и единственности. Поэтому решение y Ci yi |
при |
||||||||
|
|
|
|
|
|
|
|
i 1 |
|
a x b будет общим, т.е. будет содержать |
все |
частные решения, |
если |
||||||
удастся подобрать произвольные постоянные Ci |
так, чтобы удовлетворялись |
||||||||
произвольно заданные начальные условия |
|
|
|
|
|||||
0 |
|
|
0 |
(n 1) |
|
0 |
|
(9) |
|
y(x0 ) y1 |
, y (x0 ) y2 ,...,y |
|
(x0 ) yn , |
|
|
||||
где x0 - любая точка (a,b).
Из условия (9) получим систему n линейных относительно Ci , i=1,..,n уравнений
8
n C y ( x ) y0 ,
i i 0 1
i 1
n C y ( x ) y0 ,
i i 0 2
i 1
.............................
n C y(n 1) ( x ) y0 .
i i 0 n
i 1
с n неизвестными Ci , определитель которой отличен от нуля, т.к. это определитель Вронского W (x0 ) для n линейно независимых решений
уравнения (6). Следовательно, эта система разрешима (и однозначно) относительно Ci при любом выборе xo (a, b) и при любых правых частях.
Следствие. Максимальное число линейно независимых решений однородного линейного дифференциального уравнения равно его порядку.
Замечание. Любые n линейно независимые частные решения л. о. д. у. n – ого порядка называются его фундаментальной системой решений.
Для построения фундаментальной системы решений произ-вольно зададим n2 чисел
yik (x0 ); i 1,..,n; k 0,1,...,n 1,
Подчинив их лишь условию
y1 (x0 ) |
|
y2 (x0 ) |
|
||||
|
(x0 ) |
|
|
(x0 ) |
|
||
y1 |
|
y2 |
|
||||
|
|
|
|
(x0 ) |
|
||
y1 |
(x0 ) |
|
y2 |
|
|||
|
. |
|
|
|
. |
|
|
y(n 1) (x |
0 |
) |
y(n 1) (x |
0 |
) |
||
1 |
|
|
2 |
|
|
||
. .
. .
. .
. .
. .
yn (x0 ) |
|
|
|
(x0 ) |
|
yn |
|
|
|
(x0 ) |
0, |
yn |
||
|
. |
|
yn(n 1) (x0 ) |
|
|
где x0 любая точка (a,b). Тогда решения yi (x) , определяемые начальными
значениями |
yk (x );i 1,..,n;k 0,1,...,n 1, |
образуют фундаментальную |
|
i 0 |
|
систему, т.к. их определитель Вронского W(x) в точке x x0 отличен от нуля и, следовательно, в силу теорем 4 и 5 решения y1(x), y2 (x),...,yn (x) линейно независимы.
Линейные неоднородные дифференциальные уравнения.
Будем рассматривать линейные неоднородные уравнения вида
y(n) p (x) y(n 1) |
.. p |
n |
(x) y f (x). |
(1) |
1 |
|
|
|
Это уравнение, сохраняя прежние обозначения, запишем в виде
L[ y] f (x).
9
Если при a x b в уравнении (1) все коэффициенты pi (x) и правая
часть f(x) непрерывны, то оно имеет единственное решение, удовлетворяющее условиям
y(m) (x0 ) y0(m) ,m 0,1,..,n 1,
где y0(m) – любые действительные числа, |
а x0 – любая точка интервала |
||
a x b . |
|
|
|
Действительно, правая часть уравнения |
|
||
y(n) p (x) y(n 1) |
p (x) y(n 2) ... p (x) y f (x) |
(1a) |
|
1 |
2 |
n |
|
В окрестности рассматриваемых начальных значений удовлетворяет условиям теоремы существования и единственности:
1)правая часть непрерывна по аргументам x, y, y ,.., y(n)
2)имеет ограниченные частные производные по всем
y(m) ; m 0,1,..,n 1, т.к. эти производные равны непрерывным по предположению на отрезке [a, b] коэффициентам pn m (x)
На начальные значения не налагается никаких ограничений. Из 2-х основных свойств линейного оператора
L[Cy] CL[ y],
L[ y1 y2 ] L[ y1 ] L[ y2 ],
где C=const, следует |
|
~ |
|
||
1. |
~ |
|
|
|
|
Сумма y |
y1 решения y неоднородного уравнения |
|
|||
|
|
L[ y] f (x) |
(1) |
||
и решения |
y1 соответствующего однородного уравнения |
L[y]=0 является |
|||
решением неоднородного уравнения (1). |
|
||||
|
Доказательство. |
|
|
||
|
|
~ |
|
~ |
|
|
~ |
L[ y y1] L[ y] L[ y1], |
|
||
|
|
] |
0 , поэтому |
|
|
но L[ y] f (x) , а L[ y1 |
|
||||
|
|
~ |
|
|
|
|
|
L[ y y1] f (x). |
|
||
2. |
Если yi является решением уравнения L[ y] fi (x)(i 1,2,..,m) , |
||||
|
m |
|
|
m |
|
то y i yi является решением уравнения L[ y] i fi |
|
||||
|
i 1 |
|
|
i 1 |
|
где i – постоянные. |
|
|
|
||
|
Доказательство. |
|
|
||
|
|
m |
|
m |
|
|
|
L i yi |
|
i L[ yi ], |
(2) |
|
|
i 1 |
|
i 1 |
|
m |
|
m |
но L[ yi ] fi [x], поэтому L i yi |
i fi (x). |
|
i 1 |
|
i 1 |
10
