Файл: Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 195

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

функция,

ь

 

 

 

 

Ч>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