- •Канонические и нормальные системы ОДУ. Порядок системы. Сведение системы к одному уравнению и наоборот.
- •Лемма Арцелы (критерий компактности).
- •Ломаные Эйлера и теорема Пеано.
- •Теорема о единственности решения задачи Коши для систем ОДУ. Следствие для ОДУ n-го порядка. Случай линейного уравнения и линейной системы.
- •Лемма о равномерной непрерывности.
- •Непрерывность решения системы ОДУ по начальным данным и параметру.
- •Линейная зависимость и независимость вектор-функций. Определитель Вронского.
- •Фундаментальная система решений (ФСР) для линейной однородной системы ОДУ. Существование ФСР и их взаимосвязь. Общее решение линейной однородной неоднородной системы.
- •Резольвента линейной системы ОДУ и ее свойства.
- •Построение линейной, однородной системы по известной ФСР. Формула Лиувилля.
- •Нахождение частного решения линейной неоднородной системы ОДУ методом вариации постоянных. Формула Коши.
- •Линейные системы с постоянными коэффициентами и методы их решения Случай нормализуемой системы.
- •Линейные системы с постоянными коэффициентами и методы их решения Случай не нормализуемой системы.
- •Нормальные линейные системы с постоянными коэффициентами и методы их решения
- •Устойчивость решения по Ляпунову и асимптотическая устойчивость. Лемма Ляпунова об устойчивости.
- •Лемма Ляпунова об асимптотической устойчивости и ее усиленный вариант.
- •Теорема Ляпунова об асимптотической устойчивости по линейному приближению.
- •Лемма Адамара.
- •Дифференцируемости решения системы ОДУ по начальным данным и параметру.
- •Первые интегралы систем дифференциальных уравнений. Задание общего решения системы с помощью полной системы первых интегралов.
- •Существование полной системы первых интегралов
- •Линейные однородные УрЧП первого порядка. Связь с первыми интегралами соответствующей системы ОДУ. Замечание о квазилинейных уравнениях.
- •Квазилинейные УрЧП первого порядка. Две леммы о характеристиках.
- •Теорема о существовании единственного решения задачи Коши для квазилинейных УрЧП первого порядка в случае пространственных переменных.
- •Оператор Штурма-Лиувилля. Его собственные функции. Лемма о нулевом собственном значении.
Фундаментальная система решений (ФСР) для линейной однородной системы ОДУ. Существование ФСР и их взаимосвязь. Общее решение линейной однородной неоднородной системы.
def: Линейно независимая система решений y1,…, yn системы (2) называется фундаментальной системой решений (ФСР) системы (2), а матрица
y11(x)
Y (x)=
y1n (x)
y1n (x)
ynn (x)
называется фундаментальной матрицей решений (ФМР) системы (2).
Теорема:
ФСР системы (2) существуют. Их бесконечно много. Все они могут быть получены из одной по формуле Z (x)=Y (x) S , где Y (x) – фиксированная ФМР, а S – матрица
n ×n , с условием det S ≠ 0 . |
|
|
||
Доказательство: |
|
|
||
◄ Фиксируем произвольную точку x0 = [a,b]. Пусть B – матрица n ×n |
с условием |
|||
det B ≠ 0 , тогда по теореме о ! решения задачи Коши j = |
|
|
! решение |
y j (x) систе- |
1,n |
||||
мы (2) ( y′ = A(x)y (2)), удовлетворяющее начальным условиям |
ykj (x0 )= bkj . В этом слу- |
чае y1,…, yn – ФСР, так как W (x)= detY (x0 )= det B ≠ 0 , поскольку матриц B с условием det B ≠ 0 бесконечно много, то и ФСР бесконечно много. Пусть теперь Y (x) – фиксиро-
ванная ФМР, а Z (x) – произвольная ФМР системы (2). Положим S =Y |
−1(x |
) Z(x ), тогда |
|
d(Y (x)S ) = |
dY (x) S = (A(x) Y (x)) S = A(x)(Y (x) S ), |
0 |
0 |
|
|
||
dx |
dx |
|
|
то есть Y (x)S – матрица из решений системы (2), причем |
|
|
|
Y (x0 )S =Y (x0 )(Y −1(x0 ) Z(x0 ))= (Y (x0 )Y −1(x0 )) Z(x0 )= Z(x0 ). |
|
По теореме о единственности решения задачи Коши (для каждого столбца) получим |
||||||||||||||||
Y (x)S ≡ Z (x) на [a,b] ► |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теорема: |
|
|
|
|
|
|
|
|
y(x)= C1 y1(x)+…+Cn yn (x) (3), где |
|||||||
Общее решение |
системы |
(2) |
имеет вид |
|||||||||||||
y1,…, yn – ФСР системы (2), а C1,…,Cn |
– произвольные постоянные (из C). |
|||||||||||||||
Доказательство: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
◄ ( ) y(x) из формулы (3) является решением системы (2) в силу ее линейности и од- |
||||||||||||||||
нородности. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
1 |
,…, y |
n |
– ФСР системы (2). |
||
( ) Пусть y(x) – произвольное решение системы (2), а |
|
y |
|
|||||||||||||
Фиксируем точку x0 [a,b] тогда система относительно C1,…,Cn : |
|
|
|
|
||||||||||||
|
|
|
|
|
1 |
|
n |
|
~ |
|
|
|
|
|
|
|
|
|
|
C1 y1 (x0 )+…+Cn y1 |
(x0 )= y1(x0 ) |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
n |
|
~ |
|
|
|
|
|
|
|
|
|
|
|
C1 yn |
(x0 )+…+Cn yn |
(x0 )= yn (x0 ) |
|
|
|
|
|
||||||
имеет единственное решение |
~ |
~ |
, то есть ее определитель равен W (x0 )≠ 0 . Поло- |
|||||||||||||
C1,…,Cn |
||||||||||||||||
~ 1 |
|
~ |
n |
(x), тогда |
y(x) – решение системы (2) в силу ее линейности |
|||||||||||
жим y(x)= C1 y |
(x)+…+ Cn y |
|
||||||||||||||
и однородности, причем |
|
|
|
~ |
(x0 ) по построению. По теореме о единственности ре- |
|||||||||||
y(x0 )= y |
||||||||||||||||
|
|
~ |
|
|
|
|
~ |
~ |
1 |
(x)+… |
|
~ |
n |
(x). ► |
||
шения задачи Коши y(x)≡ y(x) на [a,b], то есть |
y |
(x)= C1 y |
|
+ Cn y |
|
15
Теорема:
Общее решение неоднородной линейной системы (1) y′ = A(x)y + f (x) имеет вид: y(x)= yчастное (x)+C1 y1(x)+…+Cn yn (x),
где yчастное (x) – частное решение системы (1), y1,…, yn – ФСР, соответствующая однородной системе (2), а C1,…,Cn – произвольные постоянные. Это следует из теории линейных пространств и предыдущей теоремы.
Резольвента линейной системы ОДУ и ее свойства.
def: Резольвентой системы (2) называют отображение R : [a,b]×[a,b]→ GL(n) – множество невырожденных матриц порядка n, удовлетворяющих следующим условиям:
1) x0 [a,b]: R(x, x0 ) – ФМР системы (2),
2) R(x0 , x0 )= E = |
1 |
0 |
. |
|
|
||
|
0 |
1 |
|
Свойства:
1° dR(x, x0 ) = A(x)R(x, x0 ) – следует из определения. dx
2° Пусть |
y(x) – |
|
решение системы |
(2) с начальными условием |
y(x0 )= y0 , |
тогда |
||||||||||||||||||
y(x)= R(x, x0 )y0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Доказательство: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
◄ Имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d(R(x, x )y |
|
) |
= |
dR(x, x |
) |
y |
|
â ñèëó 1o |
)y |
|
– решение системы (2), |
|
|
||||||||||
|
|
|
0 |
|
0 |
|
|
0 |
|
0 |
= |
A(x)R(x, x |
0 |
|
|
|||||||||
|
|
dx |
|
|
|
|
|
|
dx |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
причем |
|
R(x0 , x0 )y0 |
≡ Ey0 = y0 . |
|
По |
теореме о |
единственности |
|
задачи |
Коши |
||||||||||||||
R(x, x0 )y0 ≡ y(x) на [a,b]. ► |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
′ |
|
′ ′′ |
|
|
|
|
|
′′ |
) (полугрупповое свойство) |
|
|
|
|
|
|
|
|
|
||||||
3° R(x, x |
)R(x , x |
)= R(x, x |
|
|
|
|
|
|
|
|
|
|||||||||||||
Доказательство: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
◄ Пусть |
|
y0 – произвольный вектор из Cn , тогда R(x′, x′′)y0 – значение в точке |
x′ |
того |
||||||||||||||||||||
решения системы (2), которое в точке |
x′′ равно y0 |
(свойство 2°). R(x, x′)R((x′, x′′)y0 ) – |
||||||||||||||||||||||
значения точки x того решения, которое в x′ равно R(x′, x′′)y0 (свойство 2°), то есть того |
||||||||||||||||||||||||
решения |
|
(в силу |
|
существования |
единственности), |
которое |
в точке |
′′ |
равно |
y0 . |
||||||||||||||
|
|
x |
||||||||||||||||||||||
R((x, x′′)y0 ) – значение |
в точке |
x |
того решения, которое в |
точке |
x′′ |
равно |
y0 |
|
||||||||||||||||
R(x, x′)R(x′, x′′)y0 = |
R(x, x′′)y0 , в силу произвольного y0 получаем требуемое. ► |
|
|
4° [R(x, x0 )]−1 = R(x0 , x). Это следует из 3°, так как R(x, x0 ) R(x0 , x)= R(x, x)= E (свойство 3°). Слева аналогично.
5° dR(x1, x)= −R(x1, x)A(x) dx
Доказательство:
16
◄ |
|
R(x1, x + x)− R(x1, x) |
= R(x , x + x) |
R(x, x1 )− R(x + x, x1 ) |
R(x |
, x) |
= |
[R(x + x, x )]−1 × |
||||||||||
|
|
|
|
|||||||||||||||
|
|
x |
|
|
|
|
|
1 |
x |
|
|
|
|
1 |
|
|
|
1 |
|
R(x + x, x )− R(x, x ) |
|
|
|
−1 |
|
|
dR(x, x ) |
|
|
|
|
||||||
× |
1 |
1 |
R(x |
, x) |
→ |
[R(x, x )] |
|
− |
1 |
R(x |
, x) |
= |
R(x , x)× |
|||||
|
|
|||||||||||||||||
|
|
x |
|
|
1 |
|
x→0 |
1 |
|
|
dx |
|
1 |
|
|
1 |
||
|
|
|
|
|
|
|
|
|
свойство 3î |
|||||||||
×(− A(x)R(x, x1 )) R(x1, x)= −R(x1, x)A(x)E ► |
|
|
|
|
|
|
|
|
|
|||||||||
|
свойство 1î |
свойство 3î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Построение линейной, однородной системы по известной ФСР. Формула Лиувилля.
Теорема:
Пусть y1(x),…, yn (x) – произвольные n-мерные векторы-функции класса C1([a,b]), причем W (x)≠ 0 x [a,b], тогда ! система вида (2), для которой y1,…, yn являются
ФСР. Доказательство:
◄ Рассмотрим систему
|
|
|
dyk |
dy1k |
dykn |
|
|
|
1 |
|
dx |
dx |
dx |
|
|
|
|
|
|
1 |
n |
|
|
|
||
|
|
|
|
|
||||
|
|
|
y1 |
y1 |
y1 |
= 0 , k =1, n . |
||
W (x) |
||||||||
|
|
|
yn |
y1n |
ynn |
|
|
|
Это система вида (2) (надо разложить по первому столбцу)
dyk = ∑n akj (x)y j ,
dx j =1
где akj (x) – известные непрерывные функции на [a,b]. Заметим, что y1,…, yn – решение
этой системы (так как при подстановке y j вместо y получаем два одинаковых столбца).►
Единственность. |
|
|
~ |
Пусть Y (x) – ФМР системы |
y′ = |
|
|
A(x)y и y′ = A(x)y . Тогда |
|||
|
Y |
(x)≡ |
~ |
|
A(x)Y (x)≡ A(x)Y (x). |
||
|
|
′ |
|
Умножаем справа на Y −1(x) |
(Y −1(x) – существует, так как detY (x)=W (x)≠ 0 |
||
~ |
|
|
|
Получим, что A(x)≡ A(x) на [a,b]. |
|
Формула Лиувиля.
Пусть y1,…, yn – ФСР системы (2), тогда W ′(x)= J1 +…+ Jn , где
|
|
y11 |
y1n |
|
y11 |
y1n |
|
|
y11 |
|
|
|
|
|
|||||
|
|
… |
n |
|
n |
n |
линейность |
n |
|
|
|
1 |
|
определите ля |
|
||||
Jk |
= |
dyk |
dyk |
= |
∑aks (x)y1s |
∑aks (x)ysn |
= |
∑aks (x) |
y1s |
|
|
dx k −я строка |
dx |
|
s=1 |
s=1 |
по строкам |
s=1 |
|
|
|
… |
|
|
|
|
|
|
y1n |
|
|
y1n |
ynn |
|
y1n |
ynn |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
k −ya |
|
|
|
|
= akk (x) W (x), stroka
так как при s ≠ k получаем две одинаковые строки
x [a,b]).
y1n
ysn =
ynn
17
|
n |
|
|
|
ÎÄÓ ñ |
ìèñÿ |
|
|
|
(x)= |
разделяющи |
||||
W ′(x)= ∑akk (x)W |
= |
Tr A(x)W (x) |
|||||
|
k=1 |
|
|
|
переменным |
è ñëåä A(x ) |
|
|
|
|
x |
|
|
|
|
|
|
|
∫ |
|
|
– Формула Лиувиля. |
|
|
W (x)=W (x )exp |
|
Tr A(t)dt |
||||
|
0 |
x0 |
|
|
|
|
|
|
|
|
|
|
Нахождение частного решения линейной неоднородной системы ОДУ методом вариации постоянных. Формула Коши.
y′ = A(x)y + f (x) (1)
y′ = A(x)y (2) – собственная однородная система.
Будем искать частное решение системы (1) в виде: |
|
|||
y = C1 |
(x)y′+…+Cn (x)yn =Y C(x) (3) |
|||
|
|
C |
(x) |
|
|
|
1 |
|
|
|
C(x)= |
|
, |
|
|
|
Cn |
(x) |
|
|
|
|
где y1,…, yn – ФСР собственной однородной системы (2), а C1(x),…Cn (x) – функции класса C1([a,b]). Подставим y из (3) в (1):
(Y C(x))′ =Y′C(x)+YC′(x) = A(x)(Y C(x))+ f (x) YC′(x)= f (x)
C′(x)=Y −1(x)f (x)
(Y −1(x) – существует, так как detY (x)=W (x)≠ 0 , x [a,b])
C(x)= ∫x Y −1(t)f (t)dt = [Y (x)= R(x, x0 )S, Y −1(t)= (R(t, x0 )S )−1 = S −1 R(x0 ,t)]=
x0
= ∫x R(x, x0 )S S −1R(x0 , x)f (t)dt = ∫x R(x,t)f (t)dt .
x0 |
|
|
|
|
x0 |
|
|
y′ = A(x)y + |
f (x) |
имеет вид: |
|||
Итак, решение задачи Коши y(x |
)= y |
|
|
|||
|
|
0 |
|
0 |
|
|
|
y(x)= R(x, x0 )y0 + ∫x R(x,t)f (t)dt |
|||||
|
|
|
|
|
|
x0 |
Линейные системы с постоянными коэффициентами и методы их решения Случай нормализуемой системы.
где
qi =
l11(p)y1 +…+ l1n (p)yn = f1 |
|
|
(1), |
ln1 |
(p)y1 +…+ lnn (p)yn = fn |
|
|
p = dxd , lij (p) – многочлен от p с постоянными коэффициентами из C. Положим
max deglij , N = q1 +…+ qn – порядок системы (формальный).
1≤i≤n
18
|
l11(p) |
… l1n (p) |
|
L(p)= |
, |
|
ln1(p) |
… lnn (p) |
Lij (p) |
– алгебраическое дополнение элемента lij (p) в матрице L(p). (p)= det L(p) – |
(дифференциальный оператор, а не число). “Умножим” i-ое уравнение из (1) на Lis (p)
слева и просуммируем по i, от 1 до n. (Если умножать элементы столбца матрицы на соответствующие алгебраические дополнения элементов того же столбца, то получим определитель матрицы, а если умножать на алгебраические дополнения другого столбца, то получим 0, на этом принципе построена обратная матрица.)
n
В итоге получим: (2) (p)ys = fs (x)= ∑Lis (p)fi (x), s =1,n . (В дальнейшем будет
i =1
обосновано, что такое “умножение” возможно при достаточно гладких fi (x)).
Система (2) распадается на n независимых ОДУ, отличающихся только правыми частями. Заметим, что каждое (достаточно гладкое) решение системы (1) (или (1одн.)) является решением системы (2) (или (2одн.)), но не наоборот. Из прошлого семестра заключа-
ем, что общее решение системы (2одн.) имеет вид: λ1,…λs – корни характеристического уравнения
y(x)= p1(x)eλ1 x +…+ ps (x)eλs x (3), где
(λ)= 0 , кратности r1,…rs , а p j (x) –
векторный многочлен, каждая компонента которого является многочленом степени ≤ rj −1. Для нахождения общего решения системы (1одн.) можно было бы подставить
y(x) из (3) в (1одн.), найти соотношение между коэффициентами многочленов p1(x),…, ps (x), выразить базисные неизвестные через свободные и обозначить эти свободные через C1,C2 ,…. Однако это слишком сложно. Для упрощения докажем лемму.
Лемма:
Пусть y(x) – решение системы (1одн.) (как показано выше, y(x) охватывается формулой (3)). Тогда каждое слагаемое в записи (3) также является решением системы (1одн.) (то есть можно привести предыдущие рассмотрения для каждой λj в отдельно-
сти, что значительно проще). Доказательство:
|
|
|
|
|
|
|
|
|
по теореме |
◄ Подставим y(x) в (1одн.): L(p)y(x)= L(p)p1(x)eλ1 x +…+ L(p)ps (x)eλs x |
= |
||||||||
= eλ1 x L(p + λ1 ) |
|
(x)+…+ eλs x L(p + λs ) |
|
(x)= 0 . Поскольку все функции xk eλi x |
о сдвиге |
||||
|
|
|
|||||||
p1 |
ps |
линейно |
независимы при различных k и различных λj (“ФСР в случае кратных корней” – из про-
шлого семестра), то все коэффициенты всех многочленов в этом соотношении равны нулю. То есть
|
|
|
|
|
по теореме |
|
|
|
|
|
L(p + λj )p j (x)≡ 0 L(p)(p j (x)eλj x ) |
= |
eλj x L(p + λj )p j (x)≡ 0 p j (x)eλj x – реше- |
||||||||
|
|
|
|
|
о сдвиге |
|
|
|
|
|
ние системы (1одн.). ►
Таким образом, общая схема решения системы (1одн.) изложена. Осталось ее обосно-
вать. Для этого положим: |
l |
(p)= a pq j + |
младшие степени p, некоторые a |
могут рав- |
||||||||||
|
|
|
|
|
|
ij |
|
|
ij |
|
|
|
ij |
|
няться нулю. |
|
|
|
|
|
|
l11(p) |
… l1n (p) |
|
|
|
|||
ω = |
|
a11 |
… a1n |
|
, |
(p)= |
|
|
=ω pN + младшие степени |
p |
||||
|
|
|
|
|||||||||||
|
|
an1 |
… ann |
|
|
|
|
|
ln1(p) |
… |
lnn (p) |
|
(N =q1 +…+qn ) |
|
|
|
|
|
(p)= |
|
|
|
|
||||||
|
|
|
|
|
|
|
ω pN + + младшие |
|
степени p . |
|
19