- •СОДЕРЖАНИЕ
- •От автора
- •§1. Основные определения
- •§2. Интегрирование в квадратурах дифференциальных уравнений
- •§3. Теоремы существования и единственности решения задачи Коши
- •§1. Существование и единственность решения
- •§2. Интегрирование в квадратурах
- •§1. Дифференциальные уравнения высших порядков
- •§2. Интегрирование уравнений с помощью степенных рядов. Асимптотика
- •§3. Системы дифференциальных уравнений в нормальной форме
- •§1. Непрерывная зависимость решений от параметров
- •§2. Дифференцирование по параметрам и начальным значениям
- •§3. Уравнения в вариациях. Линеаризация
- •§1. Линейные однородные уравнения n-го порядка
- •§2. Линейные неоднородные уравнения
- •§1. Линейные однородные системы
- •§2. Линейные неоднородные системы
- •§3. Функции от матриц. Матричная экспонента
- •ГЛАВА 7. ЛИНЕЙНЫЕ УРАВНЕНИЯ ВТОРОГО ПОРЯДКА
- •§1. Линейное уравнение второго порядка
- •§2. Краевая задача. Функция Грина
- •§3. Собственные значения и интегральные уравнения
- •§4. Интегрирование уравнений с помощью рядов Фурье
- •§5. Уравнение Бесселя
- •ГЛАВА 8. ЭЛЕМЕНТЫ ТЕОРИИ УСТОЙЧИВОСТИ
- •§1. Устойчивость по Ляпунову
- •§2. Устойчивость линейных систем
- •§3. Устойчивость по первому приближению
- •§4. Нелинейные системы. Функция Ляпунова. Теоремы Ляпунова
- •§5. Теоремы Барбашина и Красовского
Предположим, что в задаче (1.11) коэффициент Q(t, ) 0, а P от не зависит: P(t, ) = P(t). Причем P(t) непрерывен при t 0; параметром считается y0 . Докажем, что если P(t) 0, то y(t, y0) непрерывна по y0 равномерно относительно t при t 0. Действительно, из формулы (1.12) имеем
t |
~ |
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
P(t )dt |
|
|
|
|
|
|
|
|
|
y(t,y0) y0 e0 |
. |
Отсюда |
сразу |
|
получаем |
следующее |
представление: |
||
|
|
|
t |
~ ~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P(t )dt |
|
|
|
|
||
y y(t,y0 y0) y(t,y0) y0 e0 |
. |
Поэтому, если |
P(t) 0, |
то |
|||||
| y | | y0 | и, таким образом, |
| y | |
при | y0 |
| ( ) для всех t 0, |
что |
|||||
и требовалось доказать. |
|
то утверждение о непрерывности y(t, y0 ) |
|||||||
Докажем, что если P(t) 0, |
|||||||||
по y0 , равномерной |
относительно |
t |
при |
t 0, неверно. |
Рассуждаем |
от |
|||
противного. Допустим, что имеет место утверждение, доказанное ранее, т.е. для
0 |
( ) такое, что |
| y | для всех |
t 0, если |
| y0 | . |
Возьмем |
|||||||||||||
| y |
0 |
|
1 |
|
|
|
t |
~ |
~ |
|
|
|
|
|
||||
|
| |
|
|
( ). В силу условия P(t) 0 имеем exp P(t )dt exp( t) |
и тем |
|||||||||||||
|
2 |
|||||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
~ ~ |
при достаточно большом t становится больше |
||||||||||
самым экспонента exp P(t )dt |
||||||||||||||||||
|
|
|
|
|
0 |
|
t |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
~ ~ |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
P(t )dt |
|
|
|
|
|
|
|
||
4 |
. Но тогда из представления |
y y0 e0 |
получаем |
| y | |
4 |
2 , |
||||||||||||
|
|
|
||||||||||||||||
|
||||||||||||||||||
|
|
|
|
|
|
|
y | . |
|
|
2 |
|
|
|
|||||
что противоречит неравенству | |
|
|
|
|
|
|
|
|
||||||||||
Замечание. Данный пример показывает, что теорема 2 не распространяется на бесконечный промежуток изменения t и для ее справедливости нужны еще некоторые дополнительные требования, например, условие на знак P(t).
§2. Дифференцирование по параметрам и начальным значениям
Кратко исследуем возможность дифференцирования решений по параметрам и начальным значениям. Не ограничивая общности, достаточно рассмотреть этот вопрос для какой-либо одной переменной . Запишем систему в виде
dyi |
= fi(y(t, ),t, ) (i = 1, …, m). |
(2.1) |
|
||
dt |
|
|
Будем считать, что выполнены условия теоремы 1. Тем самым решение начальной задачи для системы (2.1) существует и является непрерывной
функцией параметра при |
| 0 | |
c . Построим конечно-разностные |
отношения – функции zi (t, ): |
|
|
55
|
|
|
|
zi (t, ) = |
yi (t, 0 ) yi (t, 0) |
, |
(2.2) |
|
|
|
|
||||
|
|
|
|
|
|
|
|
которые являются решениями следующей системы: |
|
||||||
|
dzi |
= |
1 |
{ fi(y(t, 0 |
),t, 0 ) fi(y(t, 0),t, 0)}. |
(2.3) |
|
|
|
|
|||||
|
dt |
|
|
|
|
||
Из (2.2) получаем следующее представление: |
|
||||||
|
yi(t, 0 ) yi(t, ) zi(t, ) . |
(2.4) |
Предположим, что дополнительно к условиям теоремы 1 функции |
fi (y,t, ) в |
|
~ |
обладают непрерывными частными производными по |
. Тогда |
области D |
||
пользуясь представлением (2.4) и тождеством Адамара, запишем (2.4) в виде
|
|
|
|
|
|
dzi |
|
= |
a |
|
(t, )z |
|
(t, ) (t, ), |
(2.5) |
||||
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
ik |
|
k |
i |
|
|
|
где |
|
|
|
|
|
dt |
|
|
k 1 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
aik (t, ) = |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
1 f |
i |
(y1(t, 0), …, |
|
|
|
|
|
|
|
|
|
|
ym (t, 0 ), 0 ) d , |
|
||||
= |
|
yk (t, 0) zk , …, |
(2.6) |
|||||||||||||||
yk |
||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
1 f |
i |
|
|
|
|
|
|
|
|||
|
|
|
|
i (t, ) |
|
|
(y1(t, 0),...,ym(t, 0), 0 |
)d . |
(2.7) |
|||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
|
|
fi |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Здесь |
|
равны нулю, если |
|
|
является одним из |
yi0 и отличны, вообще |
||||||||||||
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
говоря, от нуля, если является одним из k . Начальные условия для функций zi (t, ) также имеют различный вид в зависимости от того, является ли
параметр 0 каким-либо из начальных значений задачи или нет. Если 0 |
|
yj0 |
|||||||||||
и поскольку |
yi (t0, 0) yi0 , |
yi (t0, 0 ) yi0 , то в силу (2.4) |
zi (t, ) 0. |
||||||||||
Если же 0 = |
yj0, то yi (t0, 0 |
) ( 0 |
) ij |
и zi (t, ) ij ( ij |
– символ |
||||||||
Кронекера). |
|
|
|
|
|
|
|
|
|
|
|
|
|
Из (2.6), (2.7) видно, что ai,k (t, ) |
и i (t, ) |
– непрерывные функции t |
|||||||||||
и при | t t0 |
| H , | | c . Действительно, aik и i зависят от t и как |
||||||||||||
сложные функции: и непосредственно, и |
через |
y(t, 0 ). |
Но в |
силу |
|||||||||
непрерывности частных производных от функций |
fi и |
доказанной |
выше |
||||||||||
непрерывности |
y(t, 0 ) по t и |
(непрерывность |
y(t, ) |
по |
t и |
||||||||
0 |
эквивалентна непрерывности |
y(t, 0 |
) |
по |
t |
и |
) |
эти |
|||||
сложные функции будут также непрерывными по t и . |
|
|
|
|
|
|
|||||||
Поэтому правые части (2.5) удовлетворяют условиям теоремы 1, из которой |
|||||||||||||
следует, что |
zi (t, ) – непрерывная функция |
t |
и |
при |
| t t0 |
| H , |
|||||||
| | c . Это означает, в частности, что существуют предельные значения
56
limzi (t, ) |
zi (t,0) lim |
|
yi (t, 0 ) yi (t, 0 ) |
|
yi |
|
, |
(2.8) |
||||||
|
|
|
||||||||||||
0 |
|
|
|
|
|
0 |
|
|
|
0 |
|
|
||
|
|
yi |
|
|
|
|
|
|
|
|
|
|
|
|
т. е. производные |
|
|
при 0. Эти производные удовлетворяют (2.5), где |
|||||||||||
|
|
|||||||||||||
aik aik (t,0), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i i (t,0). Из |
(2.6), (2.7) видно, |
что |
|
эти |
величины |
|||||||||
представляют собой |
|
|
fi |
(y(t, 1,..., k ,y10,...,ym0),t, 1,..., s) и |
fi |
() при |
||||||||
|
|
|
||||||||||||
|
|
|
|
|||||||||||
|
|
|
|
|
yk |
|
|
|
|
|
|
|
||
0. Поэтому, пользуясь той же теоремой 1, можно сделать заключение о
непрерывной зависимости производных yi от t параметров и начальных
значений.
Теорема |
|
3. |
Если |
функции |
fi (y,t, ) |
непрерывны |
вместе |
с |
частными |
|||||||||||||
производными |
по |
y1,..., ym |
1,..., s |
в |
~ |
то |
существуют |
|
производные от |
|||||||||||||
D, |
|
|||||||||||||||||||||
решения задачи (1.1), (1.2) по начальным значениям |
y10,..., ym0 |
и параметрам |
||||||||||||||||||||
,..., |
s |
, |
непрерывные при |
|t t |
c |
| H , |
| y |
i0 |
y0 |
| |
i |
|, |
| |
k |
|
k0 |
| c |
|||||
1 |
|
|
|
|
|
|
|
|
|
i0 |
|
|
|
|
|
k |
||||||
(i 1,...,m; |
k 1,...,s). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Замечание. В теореме сформулированы условия существования первых непрерывных производных по параметрам решения задачи Коши. Вопрос о существовании и непрерывности производных высших порядков по параметрам исследуется аналогично. Можно показать, что существование непрерывных частных производных до порядка k функций fi (y,t, ) является достаточным условием существования непрерывных частных производных
k -го порядка по параметрам y10,..., ym0 , |
1,..., s |
решения задачи (1.1), (1.2). |
|||||||||||
|
|
|
|
|
|
|
|
|
k y |
(i 1,...,m; |
|||
Существование непрерывных частных производных |
|
i |
|
||||||||||
k |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||
k 1,...,r 1) позволяет при достаточно |
малых 0 |
искать решение |
|||||||||||
задачи (1.1), (1.2) в виде асимптотического разложения |
|
|
|
||||||||||
r |
k |
yi |
|
|
( ) |
k |
|
|
|
||||
yi(t, ) yi(t, 0) |
|
|
(t, 0) |
|
O[( )r 1]. |
(2.9) |
|||||||
|
|
|
k! |
|
|||||||||
k 1 k |
|
|
|
|
|
|
|
||||||
Идея разложения решения по малому параметру , от которого правые части зависят регулярно, принадлежит выдающемуся французскому математику Анри Пуанкаре. На формуле (2.9) основаны многие методы асимптотического исследования и численного интегрирования дифференциальных уравнений, а также некоторые методы теории устойчивости.
57
§3. Уравнения в вариациях. Линеаризация |
|
|
|||||||||||
|
|
|
Предельные при 0 |
|
значения функций zi (t, ) – функции zi (t,0) = |
||||||||
|
yi |
|
|
– являются решениями задачи Коши (для определенности рассмотрим |
|||||||||
|
|
||||||||||||
|
|
0 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
fi и равно yj0 ): |
|
|
случай, когда не входит явно в правые части |
|
||||||||||||
|
|
|
|
|
dz |
i |
|
m |
f |
i |
|
|
|
|
|
|
|
|
|
|
|
|
(y,t)zk , |
zi (t0,0) ij , |
(3.1) |
||
|
|
|
|
|
dt |
|
|
|
|||||
|
|
|
|
|
|
k 1 |
yk |
перехода при 0 |
|
||||
которую |
получим в результате предельного |
в (2.5). |
|||||||||||
Нетрудно видеть, эта же система может быть получена путем формального дифференцирования исходной задачи (1.1), (1.2) по параметру . Уравнения (3.1) часто называются системой уравнений в вариациях относительно рассматриваемого решения y(t, ).
Переход к уравнениям в вариациях связан с идеей линеаризации уравнений (1.1) в окрестности некоторого выбранного решения. Пусть, например, свойства
какого-то частного решения |
yi |
(i 1,...,m) нам известны и нужно исследовать |
||||||||||||||||||||
другое решение |
|
|
|
|
|
этого же уравнения, близкое к yi |
. Так, например, |
|||||||||||||||
yi |
|
(i 1,...,m) |
||||||||||||||||||||
yi |
и |
|
|
|
|
могут различаться по начальным значениям и т.п. |
|
Введем величины |
||||||||||||||
|
yi |
|||||||||||||||||||||
yi |
|
|
|
yi (i 1,...,m). Предполагая, что правые части уравнений достаточно |
||||||||||||||||||
yi |
||||||||||||||||||||||
гладки, подставим |
|
|
yi yi |
в (2.1) и разложим fi |
по степеням yi . Тогда |
|||||||||||||||||
|
yi |
|||||||||||||||||||||
получается уравнение вида |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
d |
|
|
|
m |
|
|
|
|
fi |
|
|
|
|
||
|
|
|
|
|
|
|
( yi) aik (t) yk i( y1 , …, ym,t), |
aik (t) |
(y,t). |
(3.2) |
||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
dt |
|
|
k 1 |
|
|
|
|
yk |
|
||||||
|
Если параметр, значениями которого различаются |
yi |
и |
|
, |
явно не |
||||||||||||||||
|
yi |
|||||||||||||||||||||
входит в |
fi , то i |
o( | y| ), |
получаем линейную систему уравнений |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( yi) aik (t) yk , |
|
|
|
|
|
|
|
(3.3) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|
|
которая представляет собой систему в вариациях, так как совпадает с точностью до обозначений с (3.1). Если рассмотрение ведется на конечном
промежутке [t0,T ], то в силу доказанных выше теорем yi и yi близки на всем [t0,T ] (так как их начальные значения или правые части уравнений отличаются
мало) и, зная yi , можно получить приближенное выражение для yi в виде yi yi , где yi определяется из линейных уравнений (3.4).
Если же, зная yi , мы хотим исследовать yi на неограниченном промежутке, то задача становится сложнее и относится к теории устойчивости. Линейные уравнения изучены значительно полнее, чем системы общего вида. Поэтому процесс линеаризации – замена (2.1) уравнениями вида (3.4) – приводит в тех случаях, когда эта операция законна, к серьезному упрощению задачи.
58
Замечание. В теореме 3 рассматривался вопрос о существовании производных решения по параметрам y10,...,ym0 , 1,..., s, входящем в начальные условия и правые части системы (1.1), (1.2). Аналогично может быть исследован вопрос
и о существовании производных |
|
dyi |
. Эти производные также являются |
||||||||
|
|
||||||||||
|
|
|
|
|
dt0 |
|
|
|
|
|
|
решениями начальной задачи для системы в вариациях |
|
||||||||||
|
d |
|
yi |
|
m |
fi |
|
yk |
|
|
|
|
|
|
|
(y,t) |
, |
(3.4) |
|||||
|
|
|
|
|
|||||||
|
|
|
t0 |
|
yk |
t0 |
|||||
|
dt |
|
k 1 |
|
|
|
|||||
получающейся дифференцированием системы (1.1) по t0 и является линейной
системой |
относительно |
yi |
. |
Чтобы получить |
начальные |
условия задачи, |
|
||||||
|
|
t0 |
|
|
|
|
заменим |
исходную систему |
дифференциальных |
уравнений |
(1.1) системой |
||
|
|
|
|
t |
|
|
интегральных уравнений yi(t) yi0 fi(y, )d . Если подставить решение y(t) |
||||||
|
|
|
|
t0 |
|
|
исходной системы, то мы получим тождества. Дифференцируя тождества по t0,
|
y |
|
|
t |
|
m f |
i |
|
y |
k |
|
|
|
|
|
||||||
получим |
i |
|
|
|
|
|
|
|
|
|
|
|
, откуда при |
t t |
|
в силу (1.2) |
|||||
|
t0 |
f (y(t0),t0) |
|
|
|
|
|
|
|
|
d |
|
|
0 |
|
||||||
|
|
|
t |
k 1 yk |
|
t0 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
будем иметь |
yi |
|
|
fi ( |
y10, …, ym0 , t0). Эти условия и представляют собой |
||||||||||||||||
|
|
||||||||||||||||||||
t0 |
|||||||||||||||||||||
|
|
|
|
t t0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
начальные данные для системы (3.4).
ИСТОРИЧЕСКАЯ СПРАВКА
Красовский Николай Николаевич. Выдающийся математик.
Ученик Е. А. Барбашина, И. Г. Малкина. Родился 7 сентября 1924 г. в Свердловске. В 1949 г. окончил Уральский политехнический институт (УПИ) (Свердловск). Специалист в области математики и механики. Избран членомкорреспондентом по Отделению механики и процессов управления (механика) в 1964 г., академиком по Отделению механики и процессов управления (теория устойчивости и регулирования) в 1968 г. Награжден золотой медалью им. А. М. Ляпунова за цикл статей «Работы по теории устойчивости и по теории оптимального управления» (1992 г.); золотой медалью им. М. В. Ломоносова за выдающиеся достижения в области математической теории управления и теории
дифференциальных игр (1996 г.). Основатель крупной научной школы. В настоящее время является главным научным сотрудником отдела динамических систем Уральского филиала института математики и механики РАН, членом президиума Российской Академии Наук. Среди его учеников академики: президент РАН Ю. С. Осипов, А. Б. Куржанский, А. И. Субботин, члены-корреспонденты РАН А. Г. Ченцов, В. Е. Третьяков и многие другие члены-корреспонденты РАН, доктора и кандидаты наук, инженеры и преподаватели.
59
