Обыкновенные дифференциальные уравнения в задачах. Учебное пособие
.pdf
Ÿ 2.3. Краевые задачи на ограниченном интервале
Рассмотрим на ограниченном интервале (a; b) R краевую за-
дачу для систем линейных дифференциальных уравнений с переменными коэффициентами
|
8 |
dy |
= A(t)y + f(t); t 2 (a; b); |
|
|
||||
|
|
|
|
|
|||||
|
dt |
|
(2:3:1) |
||||||
Будем |
: |
|
y |
|
A(t) |
t=b = φ: |
n n |
||
|
< |
|
t=a + y |
|
|
||||
|
|
L j |
|
R j |
|
|
|
||
|
предполагать, что |
|
матрица размера |
|
, элемен- |
||||
|
|
|
|
|
|
|
|
|
|
ты которой принадлежат классу C[a; b], f(t) 2 C[a; b], L è R постоянные матрицы размера n n, φ постоянный вектор.
Определение. Краевая задача (2.3.1) называется однородной, åñëè f(t) 0 è φ = 0.
Ниже мы сформулируем теорему, которая дает необходимое и достаточное условие существования и единственности решения
y(t) 2 C1(a; b) \ C[a; b]
краевой задачи (2.3.1) при любых данных f(t) 2 C[a; b] è φ.
Теорема 2.3.1. Следующие три утверждения эквивалентны:
1)однородная краевая задача (2.3.1) имеет только тривиальное решение y(t) 0;
2)краевая задача (2.3.1) разрешима при любых f(t) 2 C[a; b] è
φ;
3) определитель ∆ = det(LY (a) + RY (b)) не равен нулю, где
Y (t) фундаментальная матрица решений для системы dydt =
A(t)y.
Следствие. Краевая задача (2.3.1) однозначно разрешима при любых f(t) 2 C[a; b] è φ тогда и только тогда, когда ∆ ≠ 0, ïðè
41
этом решение имеет вид
y(t) = Y (t)(LY (a) + RY (b)) 1(φ R g) |
|
+ ∫a t Y (t)(Y ( )) 1f( ) d ; |
(2:3:2) |
ãäå
∫b
g = Y (b)(Y ( )) 1f( ) d :
a
Замечание 1. Условие однозначной разрешимости:
∆ = det(LY (a) + RY (b)) ̸= 0 |
(2:3:3) |
не зависит от выбора фундаментальной матрицы решений Y (t).
Замечание 2. Åñëè ∆ = 0, то краевая задача (2.3.1) разрешима не при любых f(t) è φ. Но если при заданных f(t) è φ
решение краевой задачи (2.3.1) существует, то оно определено не однозначно.
Замечание 3. Во многих учебниках рассматриваются краевые задачи вида
|
|
dy |
(a; b); |
|
8 dt = A(t)y + f(t); t 2 |
||||
> |
|
|
|
|
> |
|
|
|
|
> |
|
|
|
|
> |
|
|
|
|
< |
Lyjt=a = φ; |
|
||
> |
Ryjt=b = ; |
|
||
> |
|
|
|
|
> |
|
|
|
|
> |
|
|
|
|
: |
|
|
|
|
ãäå L постоянная матрица размера l n, R постоянная матрица размера r n, причем l + r = n. Очевидно, такие задачи входят в класс краевых задач вида (2.3.1).
42
Упражнение 1. Показать, что краевая задача
|
> |
dy |
|
|
0 |
|
1 |
0 |
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
8 |
|
= |
0 |
|
|
1 y + 0 |
1 ; |
t |
2 |
(0; 1); |
|
dt |
|
2 |
||||||||
|
> |
|
|
|
|
0 |
f(t) |
|
|
||
имеет |
> |
|
|
= y (1) = 0 |
A |
|
|
удовлетворяет |
|||
|
> y1(0) |
@ |
1 |
|
A @ |
f(t) |
|
||||
|
решение тогда и только тогда, когда |
|
|||||||||
|
: |
|
|
|
|
|
|
|
|
||
условию
∫1
f(t) sin( t) dt = 0:
0
Рассмотрим теперь на интервале (a; b) R краевую задачу для
линейных дифференциальных уравнений произвольного порядка с переменными коэффициентами
|
8 |
|
|
|
d |
|
|
|
|
|
|
|
|
||||
|
|
|
|
Pn (t; |
|
) u = h(t); |
t 2 (a; b); |
|
|
|
|||||||
ãäå |
: |
dt |
|
|
(2:3:4) |
||||||||||||
L j |
|
R j |
|
|
|
|
|
|
|
|
|||||||
|
< |
|
u t=a + u t=b = φ; |
|
|
|
|
|
|
||||||||
|
|
d |
|
dn |
|
|
dn 1 |
|
|
d |
|
||||||
|
Pn (t; |
|
) = |
|
+ a1(t) |
|
+ + an 1(t) |
|
|
+ an(t)I; |
|||||||
|
dt |
dtn |
dtn 1 |
dt |
|||||||||||||
aj(t) 2 C[a; b], j = 1; : : : ; n, h(t) 2 C[a; b], |
|
|
1 |
|
|||||||||||||
|
|
|
|
0 |
u′(t) |
1 |
|
0 φ2 |
|
||||||||
|
|
u(t) = B |
u(t) |
C ; |
φ1 |
C |
|
||||||||||
|
|
. |
|
φ = B . |
; |
||||||||||||
|
|
|
|
|
B |
|
|
|
C |
B |
n |
C |
|
||||
|
|
|
|
@ |
|
|
|
A |
@ |
|
|
A |
|
||||
|
|
|
|
|
B u(n 1)(t) |
C |
B φ |
|
|
C |
|
||||||
L è R постоянные матрицы размера n n.
Определение. Краевая задача (2.3.4) называется однородной, если h(t) 0, φi = 0, i = 1; : : : ; n.
43
Решение краевой задачи (2.3.4) естественно искать в классе функций Cn(a; b) \ Cn 1[a; b]. Следующая теорема дает необхо-
димое и достаточное условие однозначной разрешимости задачи (2.3.4) при любых данных h(t) 2 C[a; b] è φ.
Теорема 2.3.2. Следующие три утверждения эквивалентны: 1) однородная краевая задача (2.3.4) имеет только тривиаль-
ное решение u(t) 0;
2) краевая задача (2.3.4) разрешима при любых h(t) 2 C[a; b],
φi, i = 1; : : : ; n;
3)определитель e
∆= det(LU(a) + RU(b)) не равен нулю, где
U(t) фундаментальная матрица решений для уравнения
( ) d
Pn t; dt u = 0.
Замечание 4. Условие однозначной разрешимости: e
∆ ≠ 0 не зависит от выбора фундаментальной матрицы решений U(t).
Вернемся теперь к краевой задаче (2.3.1) для системы дифференциальных уравнений. Предположим, что условие однозначной разрешимости (2.3.3) выполнено. Тогда, как отмечалось выше, решение краевой задачи (2.3.1) может быть записано в виде (2.3.2).
В частности, при φ = 0 оно имеет вид
y(t) = Y (t)(LY (a) + RY (b)) 1Rg |
|
+ ∫a t Y (t)(Y ( )) 1f( ) d : |
(2:3:5) |
Поскольку
∫b
g = Y (b)(Y ( )) 1f( ) d ;
a
перепишем формулу (2.3.5) следующим образом: y(t) = Y (t)(LY (a) + RY (b)) 1
44
[ ∫t ∫b ]
(LY (a) + RY (b)) (Y ( )) 1f( ) d R Y (b) (Y ( )) 1f( ) d
a a
= Y (t)(LY (a) + RY (b)) 1
[ ∫t ∫b ]
LY (a) (Y ( )) 1f( ) d R Y (b) (Y ( )) 1f( ) d :
a t
Рассмотрим матрицу
G(t; ) = 8 |
|
|
|
1 |
|
1 |
|
Y (t)(LY (a) + RY (b)) |
1 LY (a)(Y ( )) |
1; |
a < < t; |
||||
Очевидно,:используя эту матрицу, формулу (2.3.5) можно пере- |
|||||||
< |
Y (t)( Y (a) + |
R |
Y (b)) |
R |
Y (b)(Y ( )) ; |
t < < b: |
|
|
L |
|
|
|
|
||
писать в виде
∫b
y(t) = G(t; )f( ) d :
a
Из определения вытекает, что матрица G(t; ) обладает следу-
ющими свойствами:
1) ïðè t ≠ имеет место тождество dtd G(t; ) A(t)G(t; );
2) матрица G(t; ) имеет скачок на диагонали квадрата (a; b)
(a; b):
G( + 0; ) G( 0; ) = I;
3) матрица G(t; ) удовлетворяет граничному условию
LG(a; ) + RG(b; ) = 0; 2 (a; b):
Определение. Матрица G(t; ), удовлетворяющая условиям 1) 3), называется матрицей Грина краевой задачи (2.3.1).
45
Упражнение 2. Показать, что матрица Грина краевой задачи (2.3.1) определяется условиями 1) 3) единственным образом.
Вернемся теперь к рассмотрению краевой задачи (2.3.4) для дифференциального уравнения. Очевидно, функция u(t) является решением краевой задачи (2.3.4) тогда и только тогда, когда вектор-функция u(t) является решением следующей краевой за-
äà÷è: |
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
e |
|
t=be |
|
t 2 (a; b); |
|
|
|
|
||||
|
|
|
|
t=a |
|
|
|
(2:3:6) |
|||||||||
|
|
|
|
|
|
|
|||||||||||
ãäå |
|
8 dt |
u = A(t)u + f(t); |
|
|||||||||||||
|
: L j |
|
+ |
R j |
|
= φ; |
|
|
|
|
|
|
|
||||
|
|
|
u |
|
u |
|
|
|
|
|
|
|
|
||||
|
0 |
0 |
|
|
1 |
|
0 |
: : : |
0 |
1 |
|
0 0 |
1 |
|
|||
|
B |
0 |
|
|
0 |
|
1 |
: : : |
0 |
C |
|
B |
0 |
C |
|
||
A(t) = |
B |
|
|
|
. |
|
|
. |
. |
.. |
. |
C |
; f(t) = |
B |
|
C |
: |
B . |
|
|
|
|
|
C |
B . |
C |
|||||||||
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
B |
|
C |
|
e |
B |
0 |
|
|
0 |
|
0 |
: : : |
1 |
C |
e |
B |
0 |
C |
|
||
B |
|
|
|
C |
B |
C |
|
||||||||||
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
B |
|
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
B |
|
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
B |
|
C |
|
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
B |
|
C |
|
|
@ |
|
|
|
|
|
|
|
|
|
|
A |
|
@ |
|
A |
|
|
B an(t) an 1(t) an 2(t) : : : a1(t)C |
|
Bh(t)C |
|
|||||||||||||
Будем предполагать, что выполнено условие однозначной разрешимости: e
∆ ≠ 0. Тогда из сказанного выше следует, что при φ = 0 решение задачи (2.3.6) может быть записано в виде
u(t) = ∫b G(t; )f( ) d ; |
||
a |
e |
e |
ãäå e
G(t; ) матрица Грина этой краевой задачи. С учетом вида вектор-функции fe(t) получаем
∫b
u(t) = g(t; )h( ) d ;
a
46
ющими свойствами: |
g(t; ) |
e |
|
ãäå g(t; ) последний элемент первой строки матрицы G(t; ). |
|||
Из определения вытекает, что функция |
|
|
обладает следу- |
10) ïðè t ̸= имеет место тождество Pn |
(t; dt) g(t; ) 0; |
||
|
|
d |
|
20) в квадрате (a; b) (a; b) функция g(t; ) и все ее производные
gt(k)(t; ), k = 1; : : : ; n |
|
2, непрерывны, а на диагонали производ- |
||||||
íàÿ (n |
1)-го порядка имеет скачок , равный 1, ò. å.: |
|
||||||
|
g( + 0; ) g( |
0; ) = 0; |
|
|
|
|||
|
g′( + 0; ) g′( |
0; ) = 0; |
|
|
||||
|
: : : : : : : : : : : : : : : : : : : : : : : : : : : |
|
|
|||||
|
g(n 2)( + 0; ) g(n 2)( |
0; ) = 0; |
|
|
||||
|
g(n 1)( + 0; ) g(n 1)( |
0; ) = 1; |
|
|
||||
30) функция g(t; ) удовлетворяет граничным условиям |
||||||||
0 |
g(a; ) |
1 |
|
0 |
g(b; ) |
1 |
|
|
B |
|
C |
|
B |
|
C |
2 |
|
L B |
|
C |
|
R B |
|
C |
|
|
B |
g′(a; ) |
C |
+ |
B |
g′(b; ) |
C = 0; |
|
(a; b): |
B |
. |
C |
|
B |
. |
C |
|
|
B |
|
C |
|
B |
|
C |
|
|
B |
|
C |
|
B |
|
C |
|
|
B |
|
C |
|
B |
|
C |
|
|
@ |
|
A |
|
@ |
|
A |
|
|
B g(n 1)(a; ) |
C B g(n 1)(b; ) |
C |
|
|
||||
Определение. Функция g(t; ), удовлетворяющая условиям 10) 30), называется функцией Грина краевой задачи (2.3.4).
Упражнение 3. Показать, что функция Грина краевой зада- чи (2.3.4) определяется условиями 10) 30) единственным образом.
Доказательство приведенных утверждений можно прочитать, например, в [2, 7].
47
Ÿ 2.4. Задача Штурма Лиувилля
Рассмотрим на ограниченном интервале (a; b) R краевую задачу для линейного дифференциального уравнения второго по-
рядка |
|
d2y |
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
(a; b); |
|
|
|
|
|
|
dy |
|
|
|
|
|
|||
|
8 dt2 |
+ q(t)y + y = 0; t 2 |
|
|
|||||||
|
> |
|
|
dt |
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
= 0; |
|
(2:4:1) |
|
|
> l1y + l2 |
|
|
|
t=a |
|
|
|
|||
|
> |
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
dt |
|
|
|
|
||
|
> |
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
ãäå |
вещественнозначная |
непрерывная функция на |
, |
|
|||||||
> |
|
|
|
|
|
|
|
|
|
||
|
> r1y + r2 dy t=b = 0; |
|
|
|
|||||||
q(t) |
: |
|
|
|
|
|
|
|
[a; b] |
|
|
числовой параметр, lj, rj, j = 1; 2, вещественные числа, причем
jl1j + jl2j ̸= 0; jr1j + jr2j ̸= 0;
т. е. краевая задача (2.4.1) невырождена.
Определение. Значение параметра , при котором краевая
задача (2.4.1) имеет нетривиальное решение, называется собственным значением краевой задачи (2.4.1), а соответствующее ему ре-
шение y(t; ) ̸ 0 собственной функцией.
Определение. Задача о нахождении собственных значений и собственных функций краевой задачи (2.4.1) называется задачей
Штурма Лиувилля.
Приведем основные утверждения о задаче Штурма Лиувилля.
1.Для любого собственного значения краевой задачи (2.4.1) существует только одна линейно независимая собственная функция y(t; ).
2.Åñëè y(t; 1) è y(t; 2) собственные функции, соответству-
48
ющие двум различным собственным значениям 1 è 2, òî
∫b
y(t; 1)y(t; 2) dt = 0;
a
ò.е. собственные функции ортогональны.
3.Собственные значения краевой задачи (2.4.1) могут быть только вещественными.
4.Краевая задача (2.4.1) имеет счетное число собственных зна- чений j, j = 1; 2; : : : , причем они ограничены снизу и неограни-
ченны сверху, т. е.
< 1 < 2 < < n < : : : ; |
jlim j = 1: |
|
!1 |
Отметим, что при доказательстве четвертого утверждения используется следующая теорема, которая представляет самостоятельный интерес.
Теорема Штурма (о сравнении). Пусть функция y(t) ̸ 0
является решением уравнения
|
d2y |
|
+ Q(t)y = 0; |
|
|
|
dt2 |
|
|||
а функция z(t) ̸ 0 решением уравнения |
|
||||
|
|
d2z |
+ q(t)z = 0; |
(2:4:2) |
|
|
|
dt2 |
|||
ãäå Q(t), q(t) 2 C[a; b], причем Q(t) q(t). Предположим, что t1, t2 2 [a; b] соседние нули функции z(t), ò. å. z(t1) = z(t2) = 0
è z(t) ≠ 0, t 2 (t1; t2). Тогда существует по крайней мере одно значение t 2 [t1; t2] такое, что y(t ) = 0.
Следствие 1. Пусть в условии теоремы выполнено строгое неравенство Q(t) > q(t), t 2 [t1; t2]. Тогда на открытом интервале (t1; t2) существует по крайней мере один нуль функции y(t).
49
Следствие 2. Любое решение z(t) ̸ 0 уравнения (2.4.2) при q(t) 0, t 2 [a; b], не может иметь на [a; b] более одного нуля.
В заключение приведем теорему В. А. Стеклова, которая имеет очень важное значение при решении краевых задач для уравнений с частными производными.
Теорема (В. А. Стеклов). Система собственных функций задачи Штурма Лиувилля (2.4.1) полна в L2(a; b).
Доказательство приведенных утверждений можно прочитать, например, в [7].
Замечание. Все перечисленные утверждения о задаче (2.4.1) обобщаются для задачи Штурма Лиувилля
> |
d |
|
|
dy |
|
|
|
|
|
dy |
|
|
|
||
> |
|
|
|
|
|
|
|
8 dt |
(p(t) dt ) + q(t)y + y = 0; t 2 (a; b); |
||||||
> |
|
|
dt |
|
|
|
|
> |
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
> l1y + l2 |
|
|
|
t=a |
= 0; |
||
> |
|
|
|
dt |
|
||
> |
|
|
|
|
|
||
> |
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
> r1y + r2 dy t=b = 0; |
|||||||
: |
|
|
|
|
|
|
|
где вещественнозначные функции
p(t) 2 C1[a; b]; p(t) > 0; q(t) 2 C[a; b]:
Упражнения (А. Кнезер).
1"
1.Пусть q(t) 4t2 + t2 , " > 0, t 2 [t0; 1), t0 > 0. Äî- казать, что любое решение уравнения (2.4.2) имеет бесконечное
число нулей.
1
q(t) 4t2 , t 2 [t0; 1), t0 > 0. Доказать, что любое решение z(t) ̸ 0 уравнения (2.4.2) имеет не более одного нуля.
50
