Дифференциальные уравнения. Устойчивость решений дифференциальные уравнения старшего порядка
.pdf
№ 4245 М И Н И С Т Е РС Т В О Н АУ К И И В ЫС Ш Е ГО О Б РА З О ВА Н И Я РФ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
«НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ «МИСиС»
ИНСТИТУТ БАЗОВОГО ОБРАЗОВАНИЯ
Кафедра математики
Е.В. Твердохлебова
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ. УСТОЙЧИВОСТЬ РЕШЕНИЙ
Дифференциальные уравнения старшего порядка
Учебное пособие
Рекомендовано редакционно-издательским советом университета
Москва 2020
УДК 517 Т26
Р е ц е н з е н т д-р физ.-мат. наук, проф. Б.Ф. Мельников
(Российский государственный социальный университет)
Твердохлебова Е.В.
Т26 Дифференциальные уравнения. Устойчивость решений: Дифференциальные уравнения старшего порядка / Е.В. Твердохлебова. – М. : Изд. Дом НИТУ «МИСиС», 2020. – 128 с.
ISBN 978-5-907226-91-3
Учебное пособие содержит необходимые теоретические материалы, решения основных типов задач и задания для самоподготовки по дифференциальным уравнениям старших порядков, а также системам дифференциальных уравнений. Рассмотрены алгоритмы решения основных типов уравнений, допускающих понижение степени, алгоритм построения решений линейных дифференциальных уравнений n-го порядка с постоянными коэффициентами и приводящихся к ним, а также основные методы построения частных решений неоднородных линейных уравнений. Для решения систем дифференциальных уравнений рассматривается метод исключения, приводящий к решению линейных дифференциальных уравнений n-го порядка. Пособие содержит описание основ применения операционного метода решения линейных уравнений и систем уравнений, основанного на использовании преобразования Лапласа для функции комплексной переменной.
Материал, изложенный разделе «Дифференциальные уравнения старшего порядка», предназначен для студентов всех специальностей. Раздел «Решение линейных дифференциальных уравнений операционным методом» предназначено для студентов специальностей 09.03.01, 09.03.03, 01.03.04, 09.03.02, 11.03.04, 22.03.01.
|
УДК 517 |
|
Е.В. Твердохлебова, 2020 |
ISBN 978-5-907226-91-3 |
НИТУ «МИСиС», 2020 |
ОГЛАВЛЕНИЕ |
|
Введение........................................................................... |
5 |
1. Дифференциальные уравнения старших порядков............... |
6 |
1.1. Уравнения, допускающие понижение порядка.............. |
7 |
1.2. Линейные дифференциальные уравнения |
|
n-го порядка................................................................. |
24 |
1.2.1. Линейные однородные уравнения |
|
с постоянными коэффициентами.................................. |
28 |
1.2.2. Однородные уравнения Эйлера |
|
и приводящиеся к ним................................................ |
37 |
1.2.3. Неоднородные линейные уравнения |
|
с постоянными коэффициентами.................................. |
44 |
1.3. Метод исключения для решения систем |
|
дифференциальных уравнений........................................ |
70 |
2. Решение линейных дифференциальных уравнений |
|
операционным методом..................................................... |
76 |
2.1. Понятие оригинала и изображения функции, |
|
преобразование Лапласа................................................. |
76 |
2.2. Свойства изображений и оригиналов.......................... |
80 |
2.2.1. Изображение для линейной |
|
комбинации оригиналов.............................................. |
80 |
2.2.2. Изображение функции f(at) по теореме подобия....... |
83 |
2.2.3. Изображение функции eatf(t) |
|
по теореме смещения изображения............................... |
84 |
2.2.4. Изображение функции η(t −a) f(t −a) |
|
по теореме смещения оригинала................................... |
86 |
2.2.5. Изображение функции tn f(t) |
|
по теореме дифференцирования изображения................. |
87 |
2.2.6. Изображение функции f(t) по теореме |
|
t |
|
об интегрировании изображения.................................. |
89 |
2.2.7. Умножение изображений.................................... |
93 |
2.3. Восстановление оригинала по изображению................ |
96 |
2.3.1.Восстановление оригинала разложением |
|
на простейшие дроби................................................... |
96 |
2.3.2. Восстановление оригинала |
|
при помощи вычетов.................................................. |
101 |
3
2.3.3. Восстановление оригинала |
|
по изображению F(p)·G(p)............................................ |
103 |
2.4. Решение линейных дифференциальных уравнений..... |
106 |
2.5. Решение систем линейных дифференциальных |
|
уравнений первого порядка............................................ |
112 |
Заключение .................................................................... |
125 |
Библиографический список.............................................. |
126 |
4
ВВЕДЕНИЕ
Учебное пособие предназначено для самостоятельного повторения и усвоения материала, который предлагается на занятиях
ваудитории. Изложение материала в данном пособии соответствует последовательности преподавания его в аудитории.
Каждый раздел сопровождается теоретическими сведениями
вобъеме, необходимом для понимания рассматриваемых алгоритмов. При необходимости изложение дополнено сведениями из других разделов математики, например линейной алгебры или основ операционного исчисления теории функции комплексной переменной.
Изложение материала в данном пособии предполагает, что студентами уже изучена часть курса по дифференциальным уравнениям первого порядка, представленная в пособии 3457. При использовании алгоритмов, описанных в этом пособии,
втексте указываются ссылки на соответствующие разделы. Каждый раздел содержит большое количество задач с подроб-
но изложенным решением и задач для самостоятельного выполнения с ответами, что окажет помощь при самостоятельном изучении предмета.
5
1. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ СТАРШИХ ПОРЯДКОВ
Обыкновенным дифференциальным уравнением (ЛДУ) n-го порядка называется уравнение вида
F(x, y, y′, y′′, ..., y ( n )) = 0,
где F – известная функция (n + 2) переменных, определенная
вобласти G Rn+2 ;
x – независимая переменная из интервала (a, b); y = y(x) – неизвестная функция;
n – порядок уравнения.
Дифференциальные уравнения, разрешенные относительно старшей производной, называют уравнениями в нормальной форме:
y(n)(x) = f(x, y, y′, y′′, ..., y(n – 1)).
Функция y(x) называется решением дифференциального уравнения n-го порядка, если она n раз непрерывно дифференцируема на промежутке (a, b) и удовлетворяет уравнению для всех x (a, b).
Дифференциальное уравнение n-го порядка имеет бесконечное множество решений, зависящих от значений констант. Общее решение имеет вид
y(x,C1,C2, , Cn ) .
Чтобы выделить единственное частное решение уравнения, при котором константы С1, С2, …, Сn принимают определенные значения, достаточно задать условия однозначности, которые в данной задаче называют начальными условиями:
y(x0 ) = y0, y′(x0 ) = y0,1, y′′(x0 ) = y0,2, ,y(n−1) (x0 ) = y0,n−1 ,
где y0, y0,1, y0,2, ,y0,n−1 – константы.
Задача отыскания частного решения, удовлетворяющего заданным начальным условиям, называется задачей Коши. При определенных ограничениях на правую часть уравнения в нормальной форме задача Коши имеет единственное решение.
6
Справедлива следующая теорема существования и един-
ственности решения задачи Коши: если функция f в правой части уравнения y(n) = f(x, y, y′, y′′, ..., y (n – 1)) и ее частные про-
изводные по переменным y, y′, y′′, ..., y(n – 1) непрерывны в обла-
сти G Rn+2 , то для любой точки (x0, y0, y0,1, y0,2, ..., y0,n-1) G на некотором интервале (x0 – h, x0 + h) существует единственное
решение y(x) уравнения, удовлетворяющее начальным условиям
y(x0) = y0, y′(x0) = y0,1, y′′(x0) = y0,2, ..., y(n – 1) (x0) = y0, n – 1.
1.1. Уравнения, допускающие понижение порядка
Рассмотрим несколько типов уравнений, допускающих понижение порядка.
I. Уравнение разрешимо относительно производной, правая часть не содержит неизвестной функции и ее производных:
y(n) = f(x).
Тогда решение находится в результате последовательного n-кратного интегрирования правой части.
Пример 1.1. Найти общее решение y(3) = 2 .
Решение. x3
Обозначим y(2) (x) = g(x) , тогда
y(3) (x) = |
dy(2) (x) |
= |
dg(x) |
, |
|||||
dx |
|
dx |
|||||||
|
|
|
|
|
|
|
|||
dg |
= |
2 |
∫dg = ∫ |
|
2 |
dx , |
|||
dx |
3 |
3 |
|||||||
|
x |
|
|
x |
|
||||
g(x) = −x12 +C1 . Далее аналогично предыдущему
|
|
|
|
g(x) = y(2) (x) = |
dy′ |
, |
|
|
||||
|
|
|
|
dx |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|||
dy′ |
|
1 |
|
|
|
|
|
1 |
|
|
||
|
= − |
|
|
+C1 |
∫dy′= ∫ |
− |
|
|
|
+C1 |
dx , |
|
dx |
x |
2 |
x |
2 |
||||||||
|
|
|
|
|
|
|
|
|
||||
7
|
|
|
|
|
|
|
|
y′= |
1 |
+C1x +C2 , |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
1 |
|
|
|
|
|
|
|
∫ |
|
|
∫ |
|
1 |
|
|
|
|
|
|
|
= |
|
|
+C1x +C2 |
dy |
= |
|
|
+C1x |
+C2 |
dx , |
|||||||||
|
dx |
x |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
C x2 |
|
|
|
|
|
|
|
|
|
|
|
|||
Ответ: y =ln |
x |
+ |
1 |
|
+C x +C |
, C |
, C |
, C |
3 |
– константы. |
|||||||||||
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
2 |
2 |
|
3 |
|
|
1 |
|
2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
II. Уравнение содержит только старшие производные неизвестной функции y(x) и не содержит саму функцию и ее младшие производные:
F(x, y(k), y(k + 1), … y(n) = 0.
Подстановка y(k) (x) =z(x) |
приводит к уравнению n–k-го по- |
||
рядка: |
|
|
|
′ |
|
(n−k) |
) =0 . |
F(x,z,z , ,z |
|
||
Пусть z =ϕ(x,C1,C2,…,Cn−k ) – общее решение этого уравнения, тогда функция y(x) может быть найдена интегрированием урав-
нения |
|
|
|
|
|
|
|
|
|
y(k) =ϕ(x,C ,C ,…,C |
) (см. п. I). |
||||||||
|
1 |
2 |
n−k |
|
|
||||
Пример 1.2. Найти решения y′′′= − |
1 |
(y′′)3 . |
|||||||
2 |
|||||||||
Решение. |
|
|
|
|
|
|
|
||
|
|
|
функцию: z(x) = y′′(x) , тогда |
||||||
Введем новую неизвестную |
|||||||||
z′(x) = y′′′(x) . |
|
|
|
|
|
|
|
|
|
Отсюда |
|
|
|
|
|
|
|
|
|
z′= − |
1 |
z3 |
dz + |
1 |
z3dx =0 . |
||||
2 |
|
||||||||
|
|
|
2 |
|
|
|
|||
Получено уравнение с разделяющимися переменными (см.
раздел 2.1 пособия 3457): |
|
|
|
|
|
|
|
|
||||||
|
dz |
|
1 |
|
|
|
dz |
+ |
1 |
|
dx =0, |
(1) |
||
3 |
+ |
=0 |
|
|
|
|
|
|||||||
3 |
2 |
|
||||||||||||
z |
|
|
|
dx |
|
z |
|
|
|
|
||||
3 |
2 |
|
|
|
||||||||||
|
z |
|
|
|
|
3 |
=0. |
|
(2) |
|||||
|
|
|
|
|
|
|
|
z |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8
Общее решение получается интегрированием уравнения (1):
− |
1 |
+ |
|
1 |
x = − |
C1 |
|
, |
|||
2 |
2 |
|
|||||||||
|
2z |
|
|
|
2 |
|
|
||||
|
z2 |
|
1 |
|
|
|
|
|
|||
|
= |
|
|
. |
|
||||||
|
x +C |
|
|||||||||
|
|
|
1 |
|
|
|
|
||||
Из соотношения |
|
|
|
|
|
|
|
|
|
|
|
z = ± |
|
1 |
|
= y′′ |
|
||||||
|
|
|
|
|
|
|
|||||
|
x +C |
|
|||||||||
|
|
|
|
|
|
|
|||||
|
|
|
1 |
|
|
|
|
|
|||
после двукратного интегрирования находим общее решение:
3
y(x) = ±34 (x +C1)2 +C2x +C3 , где C1, C2, C3 – константы.
Рассмотрим решение алгебраического уравнения (2): z = y′′=0 ,
откуда после двукратного интегрирования найдем
y(x) = C4 x + C5,
где C4, C5 – константы.
Очевидно, что это решение невозможно получить из общего ни при каких значениях констант и поэтому решение является самостоятельным.
|
4 |
3 |
|
|
|
|
|
||
Ответ: y(x) = ± |
(x +C ) |
2 |
|
+C x +C , |
y(x) =C x +C . |
||||
|
|||||||||
|
3 |
1 |
|
|
2 |
3 |
4 |
5 |
|
|
|
|
|
|
|
|
|
||
Пример 1.3. |
|
Найти |
|
общее |
решение |
уравнения |
|||
(1+x2 )y′′+2xy′=12x3 .
Решение.
Дифференциальное уравнение не содержит y(x). Замена функции y′=z(x) , y′′=z′(x) .
Получим дифференциальное уравнение (1+x2 )z′+2xz =12x3 ,
которое можно записать также в виде
9
z′+ |
2x |
12x3 |
||
|
z = |
|
, |
|
1+x2 |
1+x2 |
|||
где 1+x2 ≠0, x R .
Теперь видно, что уравнение является линейным относительно z(x) (см. разд. 2.3 пособия 3457).
Линейное уравнение может быть решено, например, методом Бернулли: будем искать решение уравнения в виде произведения двух функций z(x) = u(x) · υ(x). Для определения вида функций u(x) и υ(x) выполним действия в соответствии с алгоритмом Бернулли.
1. Решаем однородное линейное уравнение:
υ′+1+2xx2 υ=0 , dxdυ = −12+xx2 υ ,
dυυ = −12+xx2 dx , ln υ = −ln 1+x2 ,
ненулевое частное решение этого уравнения υ(x) =1+1x2 .
2. Полученное частное решение подставляем в исходное урав-
нение, чтобы найти u(x): |
|
|
|
|
|
|
|
|
|||
|
|
u′υ= |
|
12x3 |
, |
|
|||||
|
|
|
|
|
|
||||||
|
|
|
|
|
1+x2 |
|
|||||
|
|
u′ |
|
1 |
|
|
= |
|
12x3 |
, |
|
|
|
|
+x2 |
1+x2 |
|||||||
|
|
1 |
|
|
|||||||
|
du |
=12x3 |
u =3x4 +C . |
||||||||
|
dx |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
10
