книги из ГПНТБ / Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие
.pdfядро (1), получим |
|
ф (X) = я X с ‘а‘ (*) + /(*)» |
(3) |
i= 1
ь
где С,- = ^ bi (s) cp (s) ofs— неизвестные числа.
а
Таким образом, решение уравнения (2) с вырожденным ядром надо искать в виде (3). Подставляя эту функцию в уравнение (2 ) и сравнивая коэффициенты при одних и тех же функциях а, (х) справа и слева, получим систему линейных алгебраических уравнений относительно С,-:
Ct — X |
П |
(i — 1,2 , ..., n), |
|
Cj<Xij + Рг |
(4) |
||
a ij = $ «< (5) bj (s) ds, |
Рг = $ / (s) bi (s) ds. |
|
|
a |
|
a |
|
Решив эту систему, мы найдем Ch а следовательно, и решение уравнения (2 ) ф (х).
§ 3. Существование решений
1 . Если ядро вырожденное, то вопрос о существован решения интегрального уравнения Фредгольма сводится к вопросу о существовании решения соответствующей системы алгебраических уравнений (4). В более общем случае мы докажем существование решения уравнения (2)
(при достаточно малых значениях |Я|) |
м е т о д о м |
п о с л е |
||
д о в а т е л ь н ы х п р и б л и ж е н и й . |
|
|
ядро |
|
Для простоты выкладок будем предполагать, что: 1) |
||||
К (х, s) непрерывно |
в квадрате а ^ х , s ^ b; |
тогда |
оно |
|
ограничено некоторой константой А, |
2) функция |
|||
f(x) непрерывна на |
отрезке [а, Ь], |
следовательно, |
она |
|
ограничена на этом отрезке некоторой константой В,
Построим последовательность функций ф1 (х), ф2 (х), ...
.... ф„(х), ... по следующему правилу:
ь
Ф1 {x) = f { x ) + ^ <\ K {X, s) Фо (s) ds, (5)
О.
где фо(s) — произвольная фиксированная интегрируемая
200
функция, |
ь |
|
|
|
|
|
Ч>2 (x) = f(x) + X\>K (х, s) фх (s) ds, |
(6) |
|
а |
|
|
Ъ |
|
|
фП(x) = f (х) + Х \ К (X, s) фя_х (s) ds, |
(7) |
|
а |
|
Т е о р е м а 1 . При значениях j Я j < ~д^ ~ ) |
последова |
|
тельность |
(5)—(7) функций ц>п(х) равномерно |
сходится |
на отрезке |
[а, Ь] к функции ф (х), являющейся решением |
|
уравнения (2 ). |
|
|
Д о к а з а т е л ь с т в о . Преобразуем формулы для полу чения функций фя(лг). Подставляя функцию фх (х) в фор мулу ДЛЯ ф2 (х), получим
ь
ф2 (х) =f(x) + $ К (х, s) f (s) ds +
a
b |
b |
-f №$ К (x, s) $ К (s, t) ф0 (t) dt ds. |
|
a |
a |
Меняя в последнем слагаемом порядок интегрирования,
получим |
ь |
ь |
|
Фа (*) =f(x) + X\Ki (X, s) / (s) ds + |
К25 К2 (х, t) ф0 (О dt, |
а |
а |
где |
|
Ki (х, s) = K(x, |
s), |
|
|
ь |
К2 |
(х, 0 = |
$ Кг (х, s) Ki (s, о ds. |
Аналогично находим |
а |
|
|
||
ь |
|
ь |
фп (X) = f { x ) + X <\jKi (X, s) f (s) ds + №5 Ki (х, s) f(s)ds + ...
a |
a |
b |
b |
... + Я » -1 5 Kn-x (x, s) / ( 0 |
ds + l n 5 Kn (X, t) Фо ( 0 dt, |
a |
a |
201
где Кп(х, t) = J Кг(х, s) K*-i(s, t) ds. Предел функции ц>п(х),
ь
если он существует, равен сумме ряда
ь
|
Ь |
|
|
... + Хп\ К п(х, s)f(s)ds + ... |
(8) |
|
а |
|
|
ь |
|
ф (х) = |
lim Хп $ Кп (х, t) ф0 (t) dt. |
|
п-*со а |
|
|
Ф (*)==/(х) + A.$Ki(x, |
s)f(s)ds + ... |
|
а
и функции
Функции Кп (х, s) называются итерированными ядрами.
Докажем равномерную сходимость этого ряда. Для этого оценим интегралы
ь
5 Кп (х, s) f(s) ds.
а
Очевидно,
ь
I Кг (X, t) I ^ $1^ 1 (*• s) ^ 1 (s>О I ds sC л 2 (b — а),
а
Ъ
I Кз (х, t) I < J ! Ki (х, s) К2 (s, t)\ds < A3 (b - а)2,
а
поэтому
Ь |
|
b |
|
5 К п {X, s) f (s) ds |
sS An(b — a)n 1 ^ i f (s) | ds sg; AnB (b — a)n. |
||
a |
|
a |
|
Следовательно, числовой ряд |
|
||
|
|
СО |
|
|
2] BAn\X\n{b- a)n |
(9) |
|
|
п= 1 |
|
|
является мажорантным для ряда (8). Если | Я | < |
— - |
||
то ряд (9) сходится. |
Следовательно, при таких X ряд (8 ), |
||
а вместе с ним |
и |
последовательность функций |
фп(х), |
202
равномерно сходится к функции ф (*) *). Эта функция явля ется решением уравнения (2). В самом деле, переходя в формуле (7) к пределу при п-> оо, получим
_ ь
Ф(*) s А, $ /С (х, s) ф (s) ds + f(x).
а
Переход к пределу под знаком интеграла здесь законен, так как последовательность сходится равномерно.
Заметим, что предел lim (p„ (х) = Ф (х) не зависит от
п —►ОО
выбора функции <p0 (s) (нулевого приближения). Из этого легко следует единственность решения уравнения (2 ). В самом деле, если существует еще одно решение ф ( х ) уравнения (2 ), то, полагая в процедуре построения функ ций (5)—(7) ф0(х)=ф(х), получим
ф1 (*)=Ф(*). ф*(лг) = |
(jc), .... |
фя (*)=ф(лг), ... |
|
||
Эта последовательность имеет пределом функцию |
ф(х). |
||||
Но вместе с тем очевидно, что |
|
|
|
||
Пт фя (х)=ф(д). |
|
|
|||
ГС—*00 |
|
|
|
|
|
Таким образом, |
|
|
|
|
|
|
Ф(х) = ф(х). |
|
|
|
|
2. Поскольку ряд |
(8) |
сходится |
при |
| X] < А ф _ ау то |
|
|
|
ОО |
|
|
|
при таких же X сходится и ряд 2 |
Ап \Х\п~1{Ь— а)^1. Но |
||||
|
|
п= 1 |
|
ОО |
|
|
|
|
|
|
|
этот ряд является мажорантным для |
ряда |
Хп~гКп |
(х, s). |
||
|
СО |
|
|
П=1 |
|
|
|
|
|
|
|
Следовательно, ряд |
2 №~гКп (х, s) сходится равномерно. |
||||
п= 1 |
|
|
|
|
|
Поэтому ряд (8 ) можно записать в виде |
|
|
|||
|
b С |
со |
|
|
|
Ф (*)=/(*) + М 1 |
2 |
(X, s)| / (s) ds |
|
||
|
|
b |
|
|
|
Ф (x)—f(x)-\-X^R(x,s,X)f(s)ds, |
(10) |
||||
a
:) Так как ф = 0.
203
0°
где функция R(x, s, X) — ХплКп(.х> s) называется резоль-
П—I
вентой уравнения (2 ).
Таким образом, если нам известна резольвента уравне ния (2 ), то по формуле (10) мы получим его решение. Мы определили резольвенту лишь для малых значений |Я|. Однако ее можно определить на любой конечной области комплексной плоскости переменной X путем аналитического продолжения (кроме, может быть, конечного числа особых точек этой области). Если это сделано, то по формуле (10) мы получим решение уравнения (2 ) для любых значений X,
кроме |
упомянутых |
особых точек. Мы не будем подробно |
||||||
на этом останавливаться *). |
описанную выше процедуру |
|||||||
3. |
|
Если мы |
применим |
|||||
уравнению Вольтерра |
|
|
|
|
|
|||
|
|
* |
|
|
|
|
|
|
|
Ф (х) = X $/<" (х, s) у {s) ds + f (х) |
{ a s ^ x s ^ b ) , |
(11) |
|||||
|
|
а |
|
|
|
|
|
|
то получим последовательность функций |
|
|
||||||
|
|
|
|
X |
|
|
|
|
|
|
Ч>1 (x) = f{x) + X \ K (х, s) ф0(s)ds, |
|
|||||
|
|
|
|
а |
|
|
|
|
|
|
|
|
х |
|
|
|
|
|
|
ф2 (х) = |
/ (х) + |
X5 К (х, s) ф! (s) ds, |
|
|||
|
|
|
|
а |
|
|
|
|
|
|
|
|
х |
|
|
|
|
|
|
ФПМ = |
/ (*) + |
X5 К (х, s) фл_! (s) ds, |
|
|||
|
|
|
|
а |
|
|
|
|
Эта последовательность |
равномерно сходится на [а, b] |
|||||||
при |
л ю б ы х з н а ч е н и я х |
п а р а м е т р а |
X. В |
самом |
||||
деле, очевидно, справедливы неравенства |
|
|
||||||
14>i (*) I < |
I / W I + 1Я i ^ | К (х, s) 11 ф0 (s) | ds = |
|
|
|||||
|
|
|
|
|
|
;5 + |
]Я,| АВ0(х — а), |
|
*) |
См. |
М и х л и н |
С. Г., |
Интегральные уравнения и их прило |
||||
жения к некоторым проблемам механики, математической физики и техники, ч. I, гл. I, Гостехиздат, 1949.
204
где |
| ф0(s) | «S B0, |
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
IФ2 |
(x) I < I / (*) I + ! |
^ ! $! к (X, s) i |
I <J>1 (s) i |
ds =s£ |
|
||
|
|
|
|
a |
|
|
|
|
^ 5 |
+ |
|^!Л ^{В + |Х| AB0 (s — a)} ds = |
|
|||
|
|
|
|
a |
|
|
|
|
|
|
|
= B + \ k\ A B (x - a ) + \k\2А2Вй^ - ~ - |
|||
Вообще, |
|
|
|
|
|
|
|
| (fn (x) | «£ В + |
1A | AВ (x — a) + . . . |
|
|
|
|||
|
... + |
1Г " 1 Ап~гВ (x — a)n 1 |
\ \ nAnB0 |
(x — a)n |
|||
|
|
|
CO |
|
( я - 1)1 |
|
n\ |
|
|
|
|
|
|
|
|
Поскольку ряд |
2 |
В I ^ I" An-x |
■ равномерно сходится |
||||
n = l
на отрезке [a, b] и его частичные суммы являются мажо рантными для функций ф„ (х), то последовательность {ф„ (*)}
также сходится равномерно; ср (х) = Нш <р„ (х), очевидно,
П-*оо является решением уравнения (1 1), и притом единственным.
§ 4. Понятие о приближенных методах решения интегральных уравнений Фредгольма второго рода
Описанный в § 3 метод последовательных приближений построения решения может служить приближенным мето дом решения интегральных уравнений Фредгольма второго рода. В качестве приближенного решения надо брать функции ф„(лг), определяемые формулами (5)—(7).
Второй метод нахождения приближенного решения интегральных уравнений состоит в том, что ядро уравне ния К (х, s) аппроксимируют с надлежащей точностью вырожденным ядром
К{х, s)= 2 м * ) Ms). (•= 1
Решение уравнения с ядром К (х, s) и будет приближен ным решением исходного уравнения *).
Третий метод, его называют ме т о д о м сеток, состо ит в следующем.
*) М и х л и н |
С. Г., Лекции по линейным интегральным урав |
нениям, Физматгиз, |
1959. |
205
Отрезки [а, b] |
изменения переменных х и s разбивают |
на п одинаковых |
частей точками деления xit s,-. Интеграл |
ь |
|
^ К {х, s) ф (s) ds в интегральном уравнении заменяют интег-
а
ральной суммой.
Получают соотношение |
|
|
|
|
|
|
|
П |
|
|
|
|
|
ф (х) |
А 2 к |
(х, |
Sj) ф7Asj + f(x). |
|
|
|
|
/'= i |
|
|
|
|
|
Полагая здесь х равным Xi |
(/ =1, 2, |
..., п), |
рассмотрим |
|||
систему уравнений |
|
|
|
|
|
|
П |
|
|
|
|
|
|
Ф; = k |
^ ‘7Ф; As/ + A- |
(/= |
1, 2, ..., |
п), |
(12) |
|
/=) |
|
|
|
|
|
|
где |
|
S/), |
fi f(Xi), ASy = Sy-+1 |
|
||
ф; = ф(хг), Kij~K(,Xit |
S j . |
|||||
Решая эту систему относительно ф,-, получим значения приближенного решения в узловых точках. Мы не будем останавливаться на подробном изложении этих и других методов приближенного решения, отсылая читателя к спе циальной литературе *).
Особого рассмотрения требуют приближенные методы решения интегральных уравнений первого рода. Это сде лано в гл. XII.
§5. Теоремы Фредгольма
Вэтом параграфе мы будем рассматривать лишь инте гральные уравнения Фредгольма второго рода
ь |
|
Ф (х) — А $ К (х, s) ф (s) ds — f (х). |
(13) |
а |
|
1. Однородное уравнение |
|
ь |
|
Ф (х) = А ^ К (х, s) ф (s)ds |
(14) |
а
при любых значениях параметра А, очевидно, имеет три
*) К а н т о р о в и ч Л. В. |
и |
К р ы л о в В. И., |
Приближенные |
методы высшего анализа, изд. |
5-е, |
гл. II, Физматгиз, |
1962. |
206
виальное решение ср (х) = 0. Однако при |
некоторых зна |
||||||||||
чениях X оно может иметь и нетривиальные решения. |
|
||||||||||
О п р е д е л е н и е . Значения параметра X, |
при которых |
||||||||||
уравнение (14) имеет нетривиальные решения |
(т. е. не рав |
||||||||||
ные тождественно |
нулю), |
называются собственными зна |
|||||||||
чениями (с. з.) |
уравнения (14) |
(ядра К (х, s)), а соответ |
|||||||||
ствующие |
им |
решения |
ср (х) — собственными |
функциями |
|||||||
(с. ф.) уравнения (ядра). |
|
|
|
|
|
|
|
||||
Справедлива следующая |
|
|
|
X не равно соб |
|||||||
Т е о р е м а |
1. |
Если в |
уравнении (13) |
||||||||
ственному |
значению соответствующего однородного урав |
||||||||||
нения (14), то уравнение (13) |
может иметь лишь един |
||||||||||
ственное решение. |
|
Пусть |
срДх) |
и |
ф2 (х) — два |
ре |
|||||
Д о к а з а т е л ь с т в о . |
|||||||||||
шения |
уравнения (13). Тогда справедливы тождества |
|
|||||||||
|
|
Фх (*) — X ^ К (х, |
s) срг (s) ds = |
/ |
(х), |
|
|||||
|
|
|
|
а |
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
ср2‘ (х) — X \ к (х, |
s) ф2 (s) ds = |
f (х), |
|
||||||
откуда |
|
|
|
а |
|
|
|
|
|
|
|
|
|
|
ь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ф1 - |
фо) - X 5 К (X, |
s) (фх - ф2) ds = |
0 . |
|
|||||
|
|
|
|
а |
|
|
|
|
|
|
|
Следовательно, |
разность |
ф (х) = фх (х) — ф2 (х) |
является ре |
||||||||
шением однородного уравнения. |
Поскольку X не является |
||||||||||
собственным значением, |
то ф (х) = ф* (х) — ф2 (х) == 0. |
Тео |
|||||||||
рема доказана. |
дальнейшего |
напомним |
|
некоторые |
теоремы |
||||||
2. |
Для |
|
|||||||||
о системах линейных алгебраических уравнений. |
|
||||||||||
Т е о р е м а |
А. |
Для того чтобы однородная система |
|||||||||
уравнений |
|
|
|
|
|
|
|
|
|
|
|
|
|
Y i aijXJ = 0 |
(/=1,2, .... |
п) |
|
(15) |
|||||
|
|
/ = i |
|
|
|
|
|
|
|
|
|
имела лишь тривиальное решение (т. е. решение, состоя щее только из нулей), необходимо и достаточно, чтобы определитель этой системы был отличным от нуля.
Те о р е м а Б. Если определитель однородной системы
(15)равен нулю, то эта система имеет р — п — г линейно независимых решений, где г — ранг матрицы системы.
207
Т е о р е м а В. Если однородная система уравнений (15) имеет лишь тривиальное решение, то соответствующая ей неоднородная система уравнений
П |
|
Y , a ijxj = bi (i= 1 , 2 , .... п) |
(16) |
имеет единственное решение при любых значениях правых частей bi.
Матрица В, полученная из матрицы А = {а;у} системы (16) путем присоединения к ней столбца элементов, стоя щих в правых частях этой системы, называется расши ренной матрицей системы (16).
Т е о р е м а Г. Для того чтобы система (16) была раз решима, необходимо и достаточно, чтобы ранг расширен ной матрицы В системы (16) был равен рангу матрицы А этой системы.
3. Как указывалось в § 4, приближенное решени уравнения (13) можно получить, заменяя это уравнение соответствующей системой линейных алгебраических урав
нений
П
фт ^ |
Дftntyj A*Sy ~ fт |
{Щ=:1 , 2 , • • • , П) (17) |
/= 1
ирешая затем эту систему.
Таким же путем известные теоремы о системах линей ных алгебраических уравнений переносятся на интеграль ные уравнения Фредгольма второго рода. Для интеграль ных уравнений эти теоремы называются теоремами Фред гольма. Ниже мы укажем один из способов получения теорем Фредгольма, не вдаваясь в подробные доказатель ства .
Функцию фд(х), равную решению системы (17) в соот ветствующих узловых точках и линейную между ними,
будем называть полигональной функцией, соответствую щей решению системы (17). Справедлива
Т е о р е м а 2. Полигональная функция <рп (х), соот ветствующая решению системы (17), равномерно стре мится при /г->оо к решению интегрального уравнения (13).
Мы опускаем доказательство этой теоремы. Теперь
опишем способ получения теорем Фредгольма. |
ядра |
|||||
Пусть |
Я |
не |
является |
собственным |
значением |
|
К (х, s). |
Тогда |
однородное |
уравнение |
(14) имеет |
лишь |
|
тривиальное |
решение. Поэтому, имея |
в виду теорему А |
||||
208
п. 2 , можно утверждать, что соответствующая система алгебраических уравнений, которой заменяется интеграль ное уравнение (14), т. е. система
|
|
т |
|
|
фт - |
А, 2 |
Яутфу kSj = 0 |
(т = 1 , 2 , . . . , л), |
|
|
|
/■= 1 |
|
|
имеет не равный нулю определитель. Следовательно, |
||||
система уравнений |
|
|||
|
|
П |
|
|
4>т - |
"X 2 |
^'«Ф/ A sJ = f m |
( т = 1,2,..., л), |
|
|
|
1= 1 |
|
|
которой заменяется неоднородное интегральное уравне |
||||
ние (13), |
имеет |
единственное |
решение. Соответствующая |
|
этому решению полигональная функция <р„(л:) при п-*-оо, |
||||
по теореме |
3, |
равномерно стремится к решению уравне |
||
ния (13). Таким образом, справедлива |
||||
1- |
я т е о р е м а Ф р е д г о л ь м а . Для всякого X, нерав |
|||
ного собственному значению, уравнение (13) имеет реше |
||||
ние, и оно единственное. |
|
|||
З а м е ч а н и е . Поскольку |
определители системы (17) |
|||
и транспонированной системы |
|
|||
|
|
П |
|
|
Фт |
|
X |
ДтуФу Asj = fm |
ipl = 1>2 , . . . , Л) |
|
|
/ = 1 |
|
|
совпадают, |
то для всякого X, не равного с. з. ядра К (х, s), |
|||
сопряженное интегральное уравнение
ь
ф (х) — X§ /( (s, х) ф (s) ds — f (х)
а
также имеет единственное решение.
Теперь обратимся к рассмотрению случая, когда X совпадает с одним из собственных значений. Справедлива
2- |
я т е о р е м а |
Ф р е д г о л ь м а . Если X является соб |
|
ственным значением ядра К (х, s), то как однородное инте |
|||
гральное |
уравнение |
(14), так и сопряженное ему уравне |
|
ние имеют конечное |
число линейно независимых решений. |
||
Эта теорема следует |
из того, что однородная система |
||
алгебраических уравнений, соответствующая уравнению (14), имеет, согласно теореме Б, конечное число линейно
независимых решений. |
Пусть X является соб |
|
3- |
я т е о р е м а Ф р е д г о л ь м а . |
|
ственным значением ядра К(х, s). Тогда, |
для того чтобы |
|
209
