Сандаков Методы решения задач по теме Интегралные уравнен 2012
.pdfция задачи (56), (57), отвечающая собственному значению λK) ,
образуют максимальную ортонормированную систему из собственных функций интегрального уравнения (63), если функции uK (x) образуют максимальную ортонормированную систему с
весом r (x) из собственных функций задачи Штурма−Лиувилля
(56), (57).
Теорема 20. Если λ = 0 не есть собственное значение задачи Штурма−Лиувилля (56), (57), то для любой функции f (x) C [a, b] решение u (x) задачи (58), (59) раскладывается в
ряд Фурье по максимальной ОНС {uK (x)}r , K =1, 2, ..., с весом
r (x) из собственных функций задачи Штурма−Лиувилля (56), (57), т.е.
∞ |
|
∞ |
( f , uK ) |
|
|||
u (x)= ∑(u, uK )r uK = |
∑ |
uK (x), |
|||||
|
|
||||||
K=1 |
K=1 |
|
λ |
K |
|
||
|
|
|
|||||
b |
|
1 |
|
b |
|
|
|
где (u, uK )r = ∫r (y)u (y)uK (y)dy = |
|
|
∫ f (y)uK (y)dy и этот ряд |
||||
|
λ |
|
|||||
a |
|
K a |
|
|
сходится абсолютно и равномерно на [a, b].
Теорема 21 (Cтеклова). Если λ = 0 не есть собственное значение
задачи Штурма−Лиувилля (56), (57), то для |
любой функции |
Φ(x) C 2 [a, b], удовлетворяющей краевым |
условиям R Φ = |
= R2Φ = 0, ее ряд Фурье |
1 |
|
∞
Φ(x)= ∑(Φ, uK )r uK (x)
K=1
по максимальной ОНС {uK (x)}r , K =1, 2, ..., с весом r (x) из соб-
ственных функций задачи Штурма−Лиувилля (56), (57) сходится к функции Φ(x) абсолютно и равномерно на [a, b].
Следствие 5. Если λ = 0 не есть собственное значение задачи Штурма−Лиувилля (56), (57), то максимальная ОНС {uK (x)}r ,
K =1, 2, ..., с весом r (x) из собственных функций задачи Штурма− Лиувилля есть ОНБ и ПОНС в L2 (r; (a, b)), т.е. для любой
51
b
функции g (x) L2 (r; (a, b)) с нормой ∫r (x) g (x)2 dx <+∞ следу-
ет, что |
|
a |
|||
|
|
||||
|
n |
|
|
|
|
|
|
|
|
||
|
g (x)− ∑(g, uK )r uK |
|
|
|
→ 0 при n → ∞. |
|
K=1 |
|
|
|
L2(r; (a, b)) |
|
|
|
Следствие 6. Если λ = 0 есть собственное значение задачи Штурма−Лиувилля (56), (57) (ему соответствует не более двух ли-
нейно независимых собственных функций u0(1) (x), u0(2) (x) ортонормированных с весом r (x)), то ортонормированная система функций {u0(i), uK}r , i =1, 2; K =1, 2, ..., с весом r (x) из собствен-
ных функций задачи |
Штурма−Лиувилля (56), (57) есть ОНБ и |
||||||
ПОНС в L2 (r; (a, b)), |
т.е. для любой функции g (x) L2 (r; (a, b)) |
||||||
следует, что |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
g (x)−(g,u0(1) )r u0(1) −(g, u0(2) )r u0(2) − ∑(g, uK )r uK |
|
|
|
→ 0 |
||
|
|
|
K=1 |
|
|
|
L2(r; (a, b)) |
|
|
|
|
|
|||
|
|
при n → ∞. |
|
|
|
|
|
|
Решение неоднородной задачи |
|
|
|
|
|
|
|
с неоднородными краевыми условиями |
|
|
|
|
||
|
Рассмотрим краевую задачу |
|
|
|
|
|
|
|
|
Lz (x)= f (x), |
(65) |
||||
|
|
R z = d ; |
(66) |
||||
|
|
1 |
1 |
||||
|
|
R2z = d2, |
|
|
|
|
|
где f (x) C [a, b], граничные операторы R1, R2 |
|
удовлетворяют |
условиям (53), (54), а d1 и d2 − постоянные.
Утверждение 11. Если λ = 0 не есть собственное значение задачи Штурма−Лиувилля (56), (57), то краевая задача (65), (66) имеет единственное решение
52
b
z (x)= ∫G (x, y) f (y)dy +C1υ1 (x)+C2υ2 (x),
|
|
a |
|
|
|
|
|
|
|
|
где C1, |
C2 |
− решение системы. |
|
|
|
|
|
|||
|
|
C R υ +C |
|
R υ |
|
= d ; |
(67) |
|||
|
|
|
1 |
1 |
1 |
2 |
1 |
2 |
1 |
|
|
|
C1R2υ1 +C2R2υ2 = d2, |
|
|||||||
где υ1, |
υ2 |
− любые линейно независимые решения |
уравнения |
|||||||
Lυ = 0; G (x, y) − функция Грина задачи (58), (59). |
|
Замечание 22. Система (67) всегда имеет и притом единственное решение.
Замечание 23. Решение краевой задачи (65), (66) можно представить в виде z (x) =u (x) + υ(x), где u (x) есть решение краевой
задачи
Lu = f ,
R1u = R2u = 0,
а υ(x) − решение краевой задачи
Lυ = 0,
R1υ = d1, R2υ = d2.
Решение неоднородной краевой задачи с параметром в случае нулевых граничных условий
Рассмотрим краевую задачу |
|
|
Lu =λr (x)u (x) + g (x), x (a, b), |
(68) |
|
|
R1u = R2u = 0, |
(69) |
здесь операторы L, |
R1 и R2 удовлетворяют условиям (52), (53), |
|
(54) p (x) C1 [a, b]; |
g (x) C [a, b]; r (x) C [a, b] |
и r (x) > 0 на |
[a, b], λ − комплексный параметр.
Из теоремы Гильберта следует, что краевая задача (68), (69) эквивалентна интегральному уравнению
53
b |
|
|
|
ω(x)= λ∫K (x, y)ω(y)dy + F (x) |
(70) |
||
a |
|
|
|
относительно функции ω(x)= |
r (x)u (x) |
с непрерывным самосо- |
|
пряженным ядром K (x, y)= |
r (x)r (y)G (x, y), |
где G (x, y) − |
|
|
|
b |
|
функция Грина задачи (56), (57), а F (x)= |
r (x) ∫G (x, y)g (y)dy . |
a
Из (70) получаем утверждение.
Утверждение 12. Если λ = 0 не есть собственное значение задачи Штурма−Лиувилля (56), (57), то:
1) при λ ≠ λK, K =1, 2, ... , где λK − собственное значение зада-
чи Штурма−Лиувилля (56), (57), существует единственное решение краевой задачи (68), (69), которое может быть найдено по 1-й формуле Шмидта для интегрального уравнения (70) (см. формулу
(47));
2) при λ = λK (λK − собственное значение кратности 1 с собственной функцией uK (x), либо кратности 2, т.е. λK = λK+1 с собственными функциями, отвечающими им uK (x), uK+1 (x)) для раз-
решимости задачи (68), (69) необходимо и достаточно, чтобы функция g (x) была ортогональна всем собственным функциям
задачи Штурма−Лиувилля (56), (57), отвечающим собственному значению λK, т.е. чтобы (g, ui )= 0, где i = K, или i = K, K +1.
При этом решение краевой задачи (68), (69) может быть найдено по 2-й формуле Шмидта для интегрального уравнения (70) (см. фор-
мулу (48)).
Замечание 24. Решение краевой задачи
Lz = λr (x) z (x) + g (x);
R1z = d1; R2z = d2
можно представить в виде z (x) = u (x) + υ(x), где u (x) есть решение краевой задачи (68), (69), а υ(x) − решение краевой задачи
54
|
Lυ = 0; |
|
|
|
(71) |
|
|
R1υ= d1, |
R2υ = d2. |
|
(72) |
||
Если λ = 0 не |
есть собственное |
значение |
задачи |
Штурма− |
||
Лиувилля (56), (57), то решение υ(x) |
задачи (71), (72) существует |
|||||
и единственно и имеет вид υ(x)= C1υ1 (x)+C2υ2 (x), |
где C1 и |
|||||
C2 − решение системы (67), а υ1 (x), |
υ2 (x) − любые линейно не- |
|||||
зависимые решения уравнения |
Lυ = 0. |
Условия существования |
||||
решения u (x) задачи (68), (69) дает утверждение 12. |
|
|||||
Условия неотрицательности и положительности |
|
|||||
спектра задачи Штурма−Лиувилля |
|
|
|
|
||
Рассмотрим задачу Штурма−Лиувилля |
|
|
|
|||
|
Lu =λr (x)u (x), |
x [a, b]; |
|
(73) |
||
|
h1u (a)−h2u′(a)= 0; |
|
(74) |
|||
|
H1u (b)+ H2u′(b)= 0, |
|
||||
где r (x) C [a, b]; |
r (x) > 0 на |
[a, b], |
hi ≥ 0, |
Hi ≥ 0, |
i =1, 2 и |
h1 + h2 > 0, H1 + H2 > 0, оператор L определяется равенством (52) и коэффициенты p (x) и q (x) в (52) удовлетворяют условиям: p (x) C1 [a, b], q (x) C [a, b].
Если λK, K =1, 2, ..., − собственные значения задачи Штурма−
Лиувилля (56), (57), то справедливы следующие теоремы. Теорема 22. Пусть в задаче Штурма−Лиувилля (56), (57)
q (x)≥ 0, p (x) > 0 на [a, b], тогда λK ≥ 0 для любого K =1, 2, ... .
Теорема 23. Если в задаче Штурма−Лиувилля (56), (57) p (x) > 0, q (x)≥ 0 на [a, b]. Тогда для того, чтобы λ = 0 было собственным
значением задачи Штурма−Лиувилля (56), (57), необходимо и достаточно, чтобы h1 = H1 = 0, q (x) = 0 на [a, b], при этом собствен-
ному значению λ = 0 соответствует собственная функция u0 ≡ const ≠ 0.
55
Пример 2.1. Найти решение краевой задачи
x |
2 |
′′ |
′ |
− y = 0 , |
|
y |
+ xy |
||
y (1) =1, y (x) |
ограничено при x →0. |
Решение. Нетрудно заметить, что заданное уравнение есть уравнение Эйлера [5]. Общее решение этого уравнения имеет вид
|
|
|
y |
|
|
= C x + |
C2 |
. |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
общ. |
|
1 |
x |
|
||
Из краевых условий находим C2 = 0, |
C1 =1. |
Следовательно, реше- |
||||||||
ние краевой задачи имеет вид y (x)= x. |
|
|||||||||
Пример 2.2. При каких a краевая задача |
|
|||||||||
|
|
|
y′′+ ay =1; |
|
|
|
||||
|
|
|
y (0)= y (1)= 0 |
|
||||||
не имеет решений. |
|
|
|
|
|
|
|
|
|
|
Решение. При |
|
a = 0 |
краевая |
задача решений не имеет. При |
||||||
a > 0, очевидно, |
общее решение уравнения |
′′ |
||||||||
y + ay =1 имеет вид |
||||||||||
y = C cos a x +C |
|
sin |
a x + |
1 . |
|
|
|
|
||
1 |
2 |
|
|
|
a |
|
|
|
|
|
Заданная краевая задача не будет иметь решения тогда и только тогда, когда система
C = − |
1 |
; |
|
|
|
|
|
1 |
a |
|
|
|
(75) |
|
|
|
|
|
sin a = − 1 |
|
C cos |
|
a +C |
2 |
|
||
|
1 |
|
|
a |
|
|
|
|
|
|
|
|
относительно C1 и C2 несовместна (эта система получена из заданных граничных условий). Система (75) несовместна при усло-
вии, если sin a = 0, cos a = −1 (если sin a = 0, то cos a = ±1).
В случае cos a =1 уравнения системы пропорциональны и, следовательно, система (75) совместна. Из последних условий находим,
что заданная краевая задача не имеет решений при a = (2n −1)2 π2.
Кроме того, как отмечено выше, при a = 0 также не существует решения краевой задачи. При a < 0 решение, очевидно, существует.
56
Пример 2.3. Построить функцию Грина следующей краевой задачи
x |
2 |
y |
′′ |
+ 2xy |
′ |
− 2y = 0 ; |
(76) |
|
|
|
|
||||||
y (1) = 0, y (x) |
ограничено при x →0. |
(77) |
Решение. Сначала покажем, что для данной краевой задачи функция Грина существует и она единственная. Для этого достаточно показать, что краевая задача (76), (77) имеет только тривиальное решение.
Заданное уравнение есть уравнение Эйлера. Стандартным методом [5] находим фундаментальную систему решений этого уравне-
ния y1 = x, y2 = x12 , т.е. общее решение уравнения (76) имеет вид
y (x)=C1x + Cx22 , где C1 и C2 − произвольные постоянные. Из
краевых условий (77) находим, что C1 = C2 = 0. Итак, задача (76),
(77) имеет только нулевое решение, а значит, для нее можно построить и притом единственную функцию Грина.
Будем искать функцию Грина G (x, y) в виде
a x + |
a2 |
, 0 ≤ x ≤ y; |
|||||
|
|
||||||
|
1 |
|
x |
2 |
|
|
|
|
|
|
|
|
|
||
G (x, y)= |
|
b2 |
|
|
|||
b x + |
|
, y ≤ x ≤1. |
|||||
|
|||||||
|
1 |
|
x |
2 |
|
|
|
|
|
|
|
|
|
Из непрерывности функции Грина
|
a |
2 |
|
|
b |
|
|
a1y + |
|
|
− b1y + |
2 |
|
= 0. |
|
y |
2 |
2 |
|||||
|
|
|
|
y |
|
|
Из условия скачка первой производной функции Грина имеем
|
|
2a |
|
|
|
2b |
|
|
1 |
|
|||
a1 |
− |
|
2 |
|
− b1 |
− |
|
2 |
|
= |
|
|
. |
y |
3 |
y |
3 |
y |
2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
Введем обозначения
CK = aK −bK,
где K =1, 2, ... .
(78)
(79)
(80)
(81)
57
Тогда из (79) и (80) для определения CK (y) получим систему
C y + |
C2 |
= 0; |
|
||||||||
|
|
|
|||||||||
|
1 |
|
y |
2 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
2C2 |
|
|
1 |
|
|||||
C − |
|
= |
, |
||||||||
|
|
2 |
|||||||||
|
1 |
y |
3 |
|
|
y |
|
||||
|
|
|
|
|
|
|
|
решение которой имеет вид C1 = 31y2 , C2 = − 3y . Далее функция
Грина должна удовлетворять заданным краевым условиям, которые в силу (78), (81) можно переписать в виде
|
|
|
|
a x + |
a2 |
|
|
ограничено при x → 0; |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
x |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
a |
− |
|
|
|
|
|
|
+ a |
2 |
+ |
|
= 0. |
|
|
|
|||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
3y |
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Откуда |
имеем |
a |
|
|
= |
0; a = |
1− y3 |
. |
|
Из |
(81) находим b |
= − |
y |
, |
||||||||||||
|
|
|
3y2 |
|
|
|
3 |
|||||||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
||||
b = |
y |
. |
Таким образом, искомая функция Грина имеет вид |
|
|
|
||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||
2 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1− y3 )x |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, 0 ≤ x ≤ y ; |
|
|
|
|||||||||
|
|
|
|
G (x, y)= 3y2 |
) |
y |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1− x |
3 |
, |
y ≤ x ≤1. |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3x2 |
|
|
|
|
|
|
|
|||
Пример 2.4. Для краевой задачи |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
y′′− y = f (x), |
|
(82) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
y′(0)= 0; |
|
|
|
(83) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
y′(2)+ y (2)= 0 |
|||||||||||||||
построить функцию Грина. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Решение. Рассмотрим однородное уравнение y′′− y = 0. |
Функ- |
|||||||||||||||||||||||||
ции |
y |
(x)= ch x |
и |
|
y |
2 |
(x)= e−x являются решениями этого одно- |
|||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
родного уравнения и удовлетворяют, соответственно, 1-му и 2-му краевым условиям (83). Следовательно, согласно замечанию 20
58
можно функцию Грина G (x, s) краевой задачи (82), (83) искать в виде
G (x, s)= a ch x, |
0 ≤ x ≤ s ; |
(84) |
b e−x, |
s ≤ x ≤ 2. |
|
Постоянные a и b будем находить из условий непрерывности функции Грина и из условия скачка производной функции Грина при x = s. Из этих условий для нахождения a и b получаем систему
a ch s −b e−s = 0;
−b e−s −a sh s =1.
Решая эту систему, получим a = −e−s, и = −ch s. Следовательно,
согласно (84) искомая функция Грина имеет вид |
|
||
|
−s |
ch x, 0 ≤ x ≤ s ; |
|
G (x, s)= −e |
|
|
|
−e−x ch s, s ≤ x ≤ 2. |
|
||
Пример 2.5. Построить функцию Грина краевой задачи |
|||
y′′− y = 0; |
(85) |
||
y′(a) = y′(b) = 0. |
(86) |
||
Решение. Функции y1 =ch (x − a) и y2 =ch (x −b) |
являются ре- |
шениями уравнения (85) и удовлетворяют, соответственно, первому и второму из условий (86). Ни одна из этих функций не удовлетворяет сразу двум условиям (86). Поэтому функция Грина существует и единственна и ее можно искать в виде
A ch (x − a), 0 ≤ x ≤ s ; G (x, s)=
B ch (x −b), s ≤ x ≤ b.
Для нахождения постоянных A и B из условия непрерывности G (x, s) при x = s и условия скачка производной функции Грина при x = s получаем систему:
A ch (s −a) − B ch (s −b) = 0;
B sh (s −b)− A sh (s − a)=1.
59
Из этой системы без труда находим |
|
|
|
|
|||||||
A = − ch (s −b) |
; |
|
B = − ch (s −b). |
|
|||||||
sh (b −a) |
|
|
|
sh (b −a) |
|
||||||
Следовательно, искомая функция Грина имеет вид |
|
||||||||||
|
ch (x −a)ch (s −b) |
, |
a ≤ x |
≤ s; |
|||||||
− |
|
|
sh (b −a) |
|
|||||||
|
|
|
|
|
|
|
|||||
G (x, s)= − ch (s −a)ch (x −b), |
s ≤ x ≤ b. |
||||||||||
|
|
|
sh (b −a) |
|
|
|
|
||||
|
|
|
|
|
|
|
|||||
Пример 2.6. Построить функцию Грина краевой задачи |
|||||||||||
′ |
′ |
− K |
2 |
x |
−1 |
y |
= 0, K N; |
(87) |
|||
(xy ) |
|
|
|
||||||||
|
|
′ |
|
|
|
|
|
|
|
|
|
y (1)+ y (1)= 0; |
|
|
|
|
(88) |
||||||
y (x) |
|
ограничено при x → 0. |
Решение. Уравнение (87) можно переписать в виде x2 y′′+ xy′−K2 y = 0.
Это уравнение Эйлера. Согласно общей методике [5] нетрудно найти линейно независимые решения y1 = xK и y2 = x−K уравне-
ния (87). Так как задача (87), (88) имеет только тривиальное решение, то функция Грина задачи (87), (88) существует и она единственная. Будем искать ее в виде
a xK + a |
2 |
x−K, 0 ≤ x ≤ s; |
||
|
1 |
|
(89) |
|
G (x, s)= |
|
|
|
|
b xK +b x−K, s ≤ x ≤1. |
||||
|
1 |
2 |
|
Из условия непрерывности, условия скачка производной функции Грина и граничных условий (88) получаем для нахождения коэф-
фициентов a1, |
a2, b1, b2 систему: |
|
|
|
|
|
|||||
a sK + a |
2 |
s−K |
−b sK −b s−K = |
0; |
|
|
|
|
|||
1 |
|
|
|
1 |
2 |
|
|
|
1 |
|
|
b KsK−1 −b Ks−K−1 −a KsK−1 |
+ a |
Ks−K−1 |
= |
; |
|||||||
1 |
+b |
|
|
2 |
|
1 |
2 |
|
|
s |
|
b |
+b K −b K = 0; |
|
|
|
|
|
|
||||
1 |
2 |
|
|
1 |
2 |
|
|
|
|
|
|
a2 |
= 0. |
|
|
|
|
|
|
|
|
|
|
60