Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие

.pdf
Скачиваний:
145
Добавлен:
22.10.2023
Размер:
13.99 Mб
Скачать

ядро (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

 

 

 

 

 

 

 

 

 

 

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 ~

{Щ=: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

Соседние файлы в папке книги из ГПНТБ