All lectures pdfs / Диф.уравнения / Лекция_9
.pdf§9. Линейные дифференциальные уравнения с переменными коэффициентами.
В этом параграфе рассмотрим Задачу Коши для линейной
системы с переменными коэффициентами: |
|
y0 = A(t)y + f(t); t [¡T; T ]; 0 < T < 1; |
|
(y(0) = y0 2 CN (или2RN ): |
(1) |
Здесь
A = (aij(t)); i; j = 1; N - матрица с непрерывными коэффициен-
тами aij0(t); t 21[¡T; T ]; f1(t)
f(t) = Bf. |
N (t)C |
- вектор-функция правых частей с непрерывны- |
||||
|
|
@ |
A |
|
|
|
ми компонентами fi(t); i = |
1; N |
на отрезке [¡T; T ]; |
||||
|
|
ByN (t)C |
- вектор-функция искомых переменных, |
|||
y(t) = 0.y1(t) 1 |
||||||
|
|
@ |
A |
|
|
|
y0 |
|
ByN0C |
|
|
|
|
= |
0.y10 |
1 - вектор начальных данных. |
||||
|
|
@ |
A |
|
|
|
В связи с вышеизложенным существуют постоянные M; N > 0, такие, что
jjA(t)jj · M; jjf(t)jj · N на отрезке [¡T; T ]:
Предположим, что Задача Коши (1) имеет непрерывное и непрерывно-
1
Лекция №9, НГУ, ММФ, 2009 |
2 |
дифференцируемое решение y = y(t). Тогда
dtd [(y(t); y(t))] = dtd £jjy(t)jj2¤ = (y0(t); y(t)) + (y(t); y0(t)) =
=(A(t)y(t) + f(t); y(t)) + (y(t); A(t)y(t) + f(t)) =
=(A(t)y(t); y(t)) + (y(t); A(t)y(t)) + (f(t); y(t)) + (y(t); f(t)) =
= (B(t)y(t); y(t)) + 2Re(f(t); y(t)); где B(t) = A(t) + A¤(t) = B¤(t),
2Re(f(t); y(t)) = (f(t); y(t)) + (f(t); y(t)). p
Поскольку jjBjj = ¸max(B¤B) = maxfj¸min(B)j; j¸max(B)jg (см. Упражнение 1 к этому параграфу),
jjBjj = jjA + A¤jj · jjAjj + jjA¤jj = 2jjAjj
(см. Упражнение 2 к этому параграфу),
2Re(f; y) · 2j(f; y)j · 2jjfjj ¢ jjyjj · jjf(t)jj2 + jjy(t)jj;
то
dtd Z2 · jjBjjZ2 + L2 + Z2 · C1(M)Z2 + N2;
где Z = jjy(t)jj, L = jjf(t)jj, C1(M) = 2M + 1.
Умножив обе части полученного неравенства на eC1t и интегрируя полученное выражение (полагая, что 0 · t · T ), в итоге получим:
Z2 |
· |
eC1tZ2 |
+ |
N |
2 |
eC1t ¡ 1 |
; Z |
|
= y |
|
0 |
|
|
C1 |
0 |
jj 0jj |
или более грубое неравенство:
Z2 |
· |
eC1T Z2 |
+ |
N |
2 |
eC1T ¡ 1 |
; 0 |
· |
t |
· |
T: |
(2) |
|
0 |
|
|
C1 |
|
|
|
Используя прием из §2 с заменой ¿ = ¡t, получим аналогичную оценку и при ¡T · t · 0. Итак, при всех t 2 [¡T; T ]
Z2 |
· |
eC1jtjZ2 |
+ |
N |
2 |
eC1jtj ¡ 1 |
· |
eC1T Z2 |
+ |
N |
2 |
eC1T ¡ 1 |
; |
t |
j · |
T: (3) |
C1 |
|
|||||||||||||||
|
0 |
|
|
0 |
|
|
C1 |
j |
|
Лекция №9, НГУ, ММФ, 2009 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
||||
Легко видеть, что в случае Задачи Коши |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
y0 = A(t)y + f(t); t · T; |
|
|
|
|
|
|
(10) |
|||||||||||||||||||||||||||
|
|
(y(t0) = y0; t0 |
|
|
|
|
( T;j j |
|
T ) |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 ¡ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
оценка (3) перепишется так: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Z2 |
· |
eC1jt¡t0jZ2 |
+ |
N |
2 eC1jt¡t0j ¡ 1 |
· |
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
C1 |
|
|
|
|
|
|
|
|
||||||||||||
· |
e2C1T Z2 |
+ |
N |
2 |
e2C1T ¡ 1 |
; |
|
|
|
t |
j · |
T: |
|
|
|
(30) |
|||||||||||||||||||||
|
0 |
|
|
|
|
C1 |
|
|
|
|
|
j |
|
|
|
|
|
|
|
||||||||||||||||||
Поскольку p |
|
· a + b при a; b > 0, то из (3); (30) следует: |
|||||||||||||||||||||||||||||||||||
a2 + b2 |
|||||||||||||||||||||||||||||||||||||
Z = jjy(t)jj · e21 C1jtjZ0 + Ns |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
e |
C1 |
jCj1¡ 1 · |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
· e21 C1T Z0 + Ns |
e |
C1T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
¡ 1 |
|
|
|
|
|
|
|
(4) |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
и |
|
|
|
C1 |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
+ Ns |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Z =· e2 C1jt¡t0jZ0 |
|
eC1 |
j |
|
|
C1 |
¡ |
|
|
|
· |
|
|||||||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t¡t0j |
1 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
· eC1T Z0 + Ns |
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
e2C1 ¡ 1 |
: |
|
|
|
|
|
|
(40) |
|||||||||||||||||||||||||||
|
|
C1 |
|
|
|
|
|
|
Как и в §2, покажем, как с помощью полученных оценок доказать единственность решения Задачи Коши (1) (если оно существует). Предположим противное: одним и тем же векторам y0; f(t) отвечают два решения yI(t); yII(t). Их разность ¢(t) = yI(t) ¡ yII(t) есть решение следующей Задачи Коши:
(¢0 = A(t)¢; jtj · T;
(5)
¢(0) = 0:
Лекция №9, НГУ, ММФ, 2009 |
4 |
Из оценки (4) для Задачи Коши (5) имеем: |
|
|||
jj¢(t)jj · e21 C1T ¢ 0 + 0 ¢ s |
|
|
|
= 0; |
|
e |
C1¡ 1 |
||
|
|
|
C1T |
|
т.е. yI(t) ´ yII(t) на отрезке [¡T; T ].
Рассмотрим теперь Задачу Коши для однородной системы (т.е.
f(t) ´ 0 при jtj · T ): |
|
|
|
y0 = A(t)y; |
t |
T; |
|
(y(t0) = y0: |
j j · |
|
(10) |
Предположим, что Задача Коши (10) однозначно разрешима при любом векторе y0. Тогда, как и в случае линейных систем с постоянными коэффициентами, можно утверждать, что:
а) Все решения линейной системы y0 = A(t)y образуют линейное пространство.
б) Размерность этого пространства равна N, т.е. размерности пространства начальных векторов y0.
в) N линейно-независимых решений системы y0 = A(t)y можно построить, решив N раз Задачу Коши (10) и взяв в качестве начальных данных N линейно-независимых векторов y0.
Пусть y[k](t); k = 1; N - какая-либо совокупность линейно-независимых
решений для системы y0 = A(t)y, при этом: |
|
|||||||||||
|
|
8 |
³[k] |
|
|
´ |
0 |
|
|
|
|
|
|
|
|
y[k] |
t |
|
= A(t)y[k](t); t |
T; |
|||||
|
|
|
|
( ) |
|
|
[k] |
; |
j j · |
|
||
|
|
<y |
|
(0) = y0 |
|
|
||||||
где |
0 |
: |
|
|
|
|
0 |
|
|
1. |
|
|
y[k](t) = |
.y1k(t) 1, y0[k] = |
.y1k0 |
|
|||||||||
|
ByNk(t)C |
|
|
|
ByNk0C |
|
||||||
|
@ |
A |
|
|
[k] |
@ |
|
|
A |
|
||
Составим из векторов y |
|
(t) матрицу Y (t): Y (t) = (yij(t)); i; j = |
1; N. Ясно, что как и в случае линейных систем с постоянными коэффициентами, если det Y (0) 6= 0, то и det Y (t) 6= 0, jtj · T . Точно
Лекция №9, НГУ, ММФ, 2009 |
5 |
так же можно утверждать, что если det Y (0) = 0, то det Y (t) = 0 при всех t 2 [¡T; T ]. Очевидно, что матрицу Y (t); det Y (t) =6 0 можно назвать фундаментальной матрицей решений для системы y0 = A(t)y. Любое решение Задачи Коши (10) выражается так:
XN
y(t) = Y (t)C = Cky[k](t);
k=1
0C1 1
где C = B@. CA, при этом вектор C следующим образом выража-
CN
ется через вектор начальных данных y0: C = Y ¡1(0)y0, т.е. решение Задачи Коши (10) записывается так:
y(t) = Y (t)Y ¡1(0)y0: |
(6) |
Сама же фундаментальная матрица решений системы y0 = A(t)y
есть решение следующей Задачи Коши для матричной системы
(
(7)
Как и в случае линейной системы с постоянными коэффициентами справедлива формула (типа формулы (8) из §3), связывающая det Y (t) и det Y0. В самом деле, если мы снова вычислим производную от ¢(t) = det Y (t) (см. §3), то получим:
|
d¢(t) |
= Ã |
N |
akk(t)!¢(t) = T rA(t)¢(t): |
|||||
|
|
||||||||
|
|
|
|
||||||
|
dt |
|
|||||||
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
|
|
|
Xk |
|
|
|
|
|
|
|
|
|
|
|
t |
T rA(s)ds |
|
|
|
|
|
|
|
|
¡ |
|
Умножая это соотношение на e |
0 |
, мы получим: |
|||||||
R |
|||||||||
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
d |
|
¡ T rA(s)ds |
|||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
dt (e |
R |
|
¢(t)) = 0; |
Лекция №9, НГУ, ММФ, 2009 |
6 |
т.е.
t |
|
|
¢(t) = eR |
T rA(s)ds |
(8) |
¢(0); |
||
0 |
|
|
где ¢(0) = det Y0. Формула (8) носит название формулы Лиувилля. Она обобщает аналогичную формулу, полученную для систем с постоянными коэффициентами.
Вернемся к формулам (6). Поскольку (см. неравенство (4)):
Z = jjy(t)jj · e12 C1(M)jtjZ0;
то
и |
jjY (t)Y ¡1(0)y0jj · e21 C1jtjjjy0jj |
|
|
|
||
|
|
jjY (t)Y ¡1(0)y0jj |
|
|
|
|
Y (t)Y ¡1 |
(0) = |
sup |
· |
e21 C1jtj: |
(9) |
|
jj |
jj |
Z0=jjy0jj6=0 |
Z0 |
|
|
|
|
|
|
|
|
|
Используя фундаментальную матрицу Y (t), получим формулу для решения Задачи Коши (1). Умножим систему y0 = A(t)y + f(t) слева на Y ¡1(t). В итоге получим (см. Упражнение 3 к этому параграфу):
т.е. |
©Y ¡1(t)y(t)ª0 = Y ¡1(t)f(t); |
|
|
|
|||||||||||
y(t) = Y (t)Y ¡1(0)y0 + Z0 t |
|
Y (t)Y ¡1(s)f(s)ds: |
(10) |
||||||||||||
Поскольку (см. формулу (9)) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
jjY (t)Y ¡1(s)jj · e21 C1jt¡sj; |
|
|
|
|
|
|||||||||
то из (10) следует: |
|
|
e21 C1jtjZ0 + ¯ |
|
|
e21 C1jt¡sj |
|
|
f(s) ds¯ |
|
|||||
Z = y(t) |
|
|
t |
|
|
|
|||||||||
|
|
|
|
|
¯0 |
|
|
1 |
|
|
|
|
|
¯ |
|
jj |
jj · |
|
|
¯ |
|
|
|
jj |
|
jj |
¯ |
|
|||
|
|
¯Z |
|
|
|
|
|
¯ · |
|
||||||
|
|
e2 C1jtjZ + 2¯ |
e |
|
j j ¡ 1 |
: |
|
¯ |
|
||||||
|
|
|
1 |
|
¯ |
|
|
2 C1 |
t |
|
|
|
|
¯ |
|
|
· |
|
|
0 |
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C1 |
|
|
|
|
|
|
Лекция №9, НГУ, ММФ, 2009 |
7 |
Это еще один вариант априорной оценки для решения Задачи Коши (1).
Кратко остановимся на случае одного линейного неоднородного уравнения высокого порядка с переменными коэф-
фициентами: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
8 |
Lx |
|
x(N) |
|
a |
|
t x(N |
¡ |
1) |
|
::: |
|
|
a |
t x |
0 |
x |
|
F t ; |
|
|||
|
= |
|
+ |
|
1 |
( ) |
|
+ |
|
|
+ |
|
|
N¡1( ) |
|
+ aN (t)t |
= |
T;( ) |
(11) |
||||
> |
|
|
|
|
|
|
|
|
|
|
|
(N |
¡ |
1) |
(0) = ® |
|
|
j j · |
|
|
|||
<x(0) = ® ; x0(0) = ® ; :::; x |
|
|
|
|
|
; |
|
|
|
||||||||||||||
: |
i |
|
1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
||
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где a (t); i = 1; N; F (t) - непрерывные функции от t на отрезке [¡T; T ]. Как и в случае линейного уравнения высокого порядка с постоянными коэффициентами, Задача Коши (11) может быть сведена к Задаче Коши вида (1):
|
8 |
d |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
y(t) = A(t)y(t) + f(t); jtj · T; |
|
|
|
|
||||||||||
|
dt |
|
|
|
(100) |
|||||||||||
|
> |
|
|
|
|
|
|
®1 |
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
= |
|
|
|
|
; |
|
|
|
|
|
|
>y(0) = y0 |
0 |
|
1 |
|
|
|
|
|
|||||||
|
< |
|
|
|
|
|
. |
|
|
|
|
|
|
|||
|
> |
|
|
|
|
|
B |
®N C |
|
|
|
|
|
|
||
|
> |
|
|
|
|
|
@ |
|
A |
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|||
где |
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
|
1 |
|
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
x t |
|
|
|
|
|
|
0 |
1 0 : : : |
|
0 |
|
||||
y(t) = |
0x0((t)) |
|
; A(t) = |
|
|
0 |
0 1 : : : |
|
0 |
; |
||||||
B. |
|
|
C |
B: : :0: : : : : :0: :0: : ::::::: : : : :1: : :C |
||||||||||||
|
B |
|
|
C |
|
|
|
B |
|
|
|
|
|
|
C |
|
|
Bx(N¡1) |
(t)C |
|
|
|
B |
|
aN (t) : : : : : : : : |
|
a1 |
C |
|
||||
|
B |
|
|
C |
|
|
|
B |
|
¡ |
(t)C |
|
||||
|
@ |
|
|
A |
|
|
|
@¡ |
|
|
|
|
A |
|
f(t) = |
0 |
.0 |
1 |
: |
|
BF (t)C |
|
||
|
@ |
|
A |
|
Понятно, что если Задача Коши (100) однозначно разрешима, то однозначно разрешима Задача Коши (11) (в силу эквивалентности этих задач). При F (t) ´ 0 (т.е. f(t) ´ 0) фундаментальная
Лекция №9, НГУ, ММФ, 2009 |
8 |
матрица решений для системы y0 |
= A(t)y такова: |
|
|
|
||||||||
©(t) = (y |
|
(t); :::; y |
|
(t)) = |
B |
x1 |
1) |
: : : |
xN |
1)C |
; |
|
|
|
(N. |
|
(N. |
||||||||
|
[1] |
|
[N] |
|
|
0 |
x10 |
|
: : : |
xN0 |
1 |
|
|
|
|
|
|
|
Bx |
¡ |
|
: : : x |
¡ |
C |
|
|
|
|
|
|
|
B |
1 |
|
|
N |
C |
|
|
|
|
|
|
|
@ |
|
|
|
|
A |
|
при этом любое решение уравнения Lx = 0 представимо так:
XN
x(t) = Ckxk(t):
k=1
Заметим, что вронскиан W (t) = det ©(t) =6 0, причем по формуле Лиувилля 9
=
a1(s)ds;; (12)
Лекция №9, НГУ, ММФ, 2009 |
9 |
Упражнения к §9
1. Доказать, что для эрмитовой матрицы A:
jjAjj = maxfj¸min(A)j; j¸max(A)jg:
2. Дана система y0 = A(t)y; t 2 [0; T ] с непрерывными коэффици-
ентами на отрезке [0; T ]. Пусть ¸max(B) · 0 на отрезке [0; T ], где B = A + A¤. Доказать, что для решения Задачи Коши (10) (если
оно существует): (y0 = A(t)y; t 2 (0; T ];
y(0) = y0
справедлива априорная оценка:
Z = jjy(t)jj · Z0 = jjy0jj; t 2 (0; T ]:
3. Рассмотрим следующие задачи:
(X0(t) = A(t)X(t); jtj · T; X(0) = IN
и |
(Y 0(t) = ¡Y (t)A(t); jtj · T; |
|
|
|
Y (0) = IN : |
Показать, что Y (t) = X¡1(t).
4. Пусть любое решение Задачи Коши
(y0 = By; t 2 R1; y(0) = y0
где B - матрица с постоянными коэффициентами, таково, что jjy(t)jj = jjy0jj. Показать, что это справедливо тогда и только тогда, когда
B¤ = ¡B. |
|
B, то |
eB |
¤ |
= e¡B = |
eB |
¡ |
1 |
, т.е. eB - |
Замечание. Если B¤ = |
|
|
|||||||
унитарная матрица. ¤ |
¡ |
|
¡ ¢ |
|
|
¡ ¢ |
|
|
|