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

3курс,VIсем Задачи инт ур

..pdf
Скачиваний:
33
Добавлен:
18.03.2016
Размер:
1.24 Mб
Скачать

что

y = Ay ,

т.е. неподвижная точка оператора Р является также и неподвижной точкой

оператора А.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для доказательства единственности неподвижной точки оператора А предположим

противное, т.е. пусть существуют две точки

y y такие что

Ay = y

и Ay = y , тогда

An y Py = y

и

An y Py = y , откуда

следует, что

y = y

в силу

единственности

неподвижной точки оператора Р.

 

 

 

 

 

 

Пример 4.3. Пусть в банаховом пространстве

В заданы два сжимающих оператора А и

D,

т.е.

для

любых

 

элементов

x,

y B

имеют

место

соотношения

|| Ax Ay || αA || x y ||,

0 αA <1 и

|| Dx Dy || αD || x y ||,

0 αD <1.

 

Зададим некоторое число ε > 0 .

Назовем операторы А и D ε -близкими, если для

z B выполнено

|| Az Dz || ε .

 

 

 

 

 

 

 

Доказать, что неподвижные точки x и y

этих операторов находятся на расстоянии

ρ(x , y ) || x y ||

 

 

ε

 

 

,

где α = max{αA ,αD } <1.

 

 

 

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Пусть x - неподвижная точка оператора А, т.е. x = Ax , а y - неподвижная точка оператора D: y = Dy .

Возьмем

x в

качестве

начального

приближения

итерационного

процесса

для

определения

y ,

т.е.

y

0

= x ,

y = Dy Dx ,...,

 

y

k

= Dk x ;

 

 

y = lim Dk x .

 

Тогда

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

k →∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|| x y

k

||

|| x y || +|| y y

2

|| +... +|| y

k 1

y

k

||

|| x y || (1+

α

D

+α2 +... +αk 1 )

|| x Dx || .

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

D

D

 

1

αD

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходя

в

последнем

неравенстве

к

 

 

пределу

 

 

при

k → ∞ ,

получим

|| x y ||

||

x Dx || =

|| Ax Dx ||

|| Ax Dx ||

 

ε

 

,

что

и

требовалось

 

 

1αD

 

1α

доказать.

 

1αD

 

 

 

 

 

 

 

 

1α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4.4. Привести пример оператора, переводящего банахово пространство

В

в

себя, удовлетворяющего

 

условию || Ay Az ||B < ||

y z ||B

 

 

для

любых

y, z B ,

и

не

имеющего неподвижной точки.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ax = x + π2 arctg x . Для

Решение. Рассмотрим следующий оператор в пространстве R1 :

любых

 

y, z B ,

применяя

формулу

Лагранжа

(ξ

 

между y

и

 

z ),

получим

|| Ay Az || 1 = | y arctg y (z arctg z) |

= |1

1

 

 

| | y z |

< |

y z | ||

y z || 1 .

Однако,

1+ξ2

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

для всех

y имеет место неравенство Ay = y + π

arctg y > y , что и означает отсутствие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

неподвижной точки рассматриваемого оператора.

Замечание. Указанный оператор не удовлетворяет теореме о неподвижной точке!!!

31

Пример 4.5. При каких λ оператор Фредгольма Ay = λ1 (x s) y(s) ds является

0

сжимающим

а) при действии C[0,1] C[0,1]; б) при действии h[0,1] h[0,1] .

Решение.

а)

Рассмотрим

две

произвольные непрерывные функции

y(x), z(x) C[0,1] , тогда

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|| Ay Az ||

= | λ | max

 

(x s)[ y(s) z(s)]ds

| λ | max | y(s) z(s) | max | x s | =

 

 

 

 

C[0,1]

 

x [0,1]

 

 

 

s [0,1]

x,s [0,1]

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

= | λ | || y z ||C[0,1] ,

т.е. оператор А: C[0,1] C[0,1]

является сжимающим при | λ |<1.

 

б)

Рассмотрим две произвольные непрерывные функции y(x), z(x) h[0,1] , тогда

 

 

 

 

 

 

1

1

2

 

(неравенство Коши-Буняковского)

 

|| Ay Az ||h2[0,1] = λ2 dx

(x s) [ y(s) z(s)]ds

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

1

1

 

1

 

 

 

 

 

1

1

 

1

 

λ2 {(x s)2 ds [ y(s) z(s)]2ds}dx = λ2 || y z ||2h[0,1] ∫∫(x s)2 ds dx = λ2 || y z ||h2[0,1]

,

 

0

0

 

 

 

0

0

 

6

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

= || yz||2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h[ 0,1]

 

 

 

 

 

 

следовательно, оператор А:

h[0,1] h[0,1] является сжимающим при условии | λ |<

6 .

 

Пример 4.6. При каких λ оператор Вольтерра By = λx (x s) y(s) ds является

0

сжимающим

а) при действии C[0,1] C[0,1]; б) при действии h[0,1] h[0,1] .

Решение.

а) Рассмотрим любые две непрерывные функции y(x), z(x) C[0,1] , тогда

 

 

 

 

 

 

x

 

 

 

 

x

 

|| Ay Az ||C[0,1] = | λ | max

 

(x s)[ y(s) z(s)]ds

| λ | max | y(s) z(s) | max

(x s) ds =

 

 

 

 

x [0,1]

 

 

 

s [0,1]

x [0,1]

 

 

 

 

 

 

 

0 0

 

 

 

 

0

 

=

| λ | || y z ||

 

, т.е. оператор В: C[0,1] C[0,1]

является сжимающим при | λ |< 2 .

 

2

C[0,1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Рассмотрим две произвольные непрерывные функции y(x), z(x) h[0,1] , тогда

 

 

 

 

1

x

 

2

(неравенство Коши-Буняковского)

|| Ay Az ||h2[0,1] = λ2 dx

(x s) [ y(s) z(s)]ds

 

 

 

 

0

0

 

 

 

 

 

 

λ2 1

[x (x s)2 ds x [ y(s) z(s)]2ds] dx = λ2 || y z ||h2[0,1] 1

dxx (x s)2 ds =

λ2 || y z ||h2[0,1] ,

 

0

0

 

 

 

 

0

0

12

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

|| yz||2

 

 

 

 

 

 

 

 

 

 

 

 

h[ 0,1]

 

 

 

 

 

 

следовательно,

оператор В: h[0,1] h[0,1] будет сжимающим при | λ | < 12 .

32

Пример 4.7. Считая параметр λ "малым", методом последовательных приближений построить резольвенту и получить решение уравнения Фредгольма

y(x) = λπsin(x + s) y(s) ds +cos x .

0

Сформулировать критерий "малого" λ .

Замечание. Приведенный ниже способ построения резольвенты служит иллюстрацией к доказанной в курсе лекций теореме о возможности представления резольвенты R(x, s, λ) в

виде ряда R(x, s, λ) = λm1Km (x, s) при "малых" λ . Так как ядро оператора Фредгольма в

m=1

данной задаче является вырожденным, то решение уравнения может быть получено более эффективно методами, рассмотренными в примерах 6.1-6.2 (тема 6). Для построения резольвенты также могут быть использованы и другие приемы (сравните с результатом примера 6.3!!!).

Решение.

Начнем с построения резольвенты R(x, s, λ) , выписав итерированные ядра:

 

K1 (x, s) K (x, s) = sin(x + s) ,

 

 

 

 

 

 

 

 

 

 

K2 (x, s) = πK (x,t)K1 (t, s) dt = πsin(x +t) sin(t + s) dt = π

cos(x s) ,

 

 

 

 

 

0

 

 

0

 

 

2

 

 

 

 

 

 

K3 (x, s) = πK (x,t)K2 (t, s) dt =

π πsin(x +t) cos(t s) dt = π2 sin(x + s) ,

 

 

 

 

 

0

 

 

2 0

 

 

 

4

 

 

 

 

 

……………….

 

 

 

 

 

 

 

 

 

 

 

π

π

2m1

 

m =1, 2,3... ,

 

 

 

 

K2m (x, s) = K (x,t)K2m1 (t, s) dt =

 

cos(x s) ,

 

 

 

 

 

0

 

 

2

 

 

 

 

 

 

 

 

 

 

π

π

2m

sin(x + s) , m

=1, 2,3... .

 

 

 

 

K2m+1 (x, s) = K (x,t)K2m (t, s) dt =

 

 

 

 

 

 

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя формулу для резольвенты

R(x, s, λ) = λm1Km (x, s)

и суммируя ряд

 

 

 

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

sin(x + s) +

πλ

cos(x s)

 

отдельно по четным и нечетным номерам,

найдем R(x, s, λ) =

 

 

2

 

.

 

 

πλ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий "малого" λ

в данном случае совпадает с условием сходимости ряда

Неймана

| λ | <

2

, и как

легко

видеть, выполняется на более

широком множестве

 

 

 

π

 

 

 

 

 

 

 

 

 

 

значений λ , чем достаточное условие, сформулированное в теореме о разрешимости

уравнения Фредгольма 2-го рода: | λ | <

 

 

1

 

=

1

.

 

 

 

 

 

 

 

 

 

 

K0

(b a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

Далее, так как

решение представимо

в виде

y(x) = f (x) +λ b R(x, s, λ) f (s) ds ,

то

полагая f (x) = cos x ,

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λπ

 

 

 

 

πλ

2

 

cos x +

λπ

sin x

 

π

 

 

 

2

sin x +

2

 

cos x

 

 

y(x) = cos x +λR(x, s, λ) cos s ds = cos x +

 

 

 

 

 

 

 

 

=

 

 

2

 

 

.

 

 

 

 

 

πλ

2

 

 

 

πλ

2

0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

33

Тот же результат дает и непосредственное применение метода последовательных

приближений.

Положим

y0 0,

 

yn+1 = λ A yn + f ,

n = 0,1, 2,... ,

тогда

y1 = λ πsin(x + s) 0 ds +cos x = cos x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2 = λ πsin(x + s) cos s ds +cos x = λ

π sin x +cos x ;

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

π

 

 

 

 

 

 

 

 

 

 

λπ 2

cos x +λ

π

sin x +cos x ;

y3 = λ sin(x + s)[λ

2

sin s +cos s]ds +cos x =

2

 

 

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

……………….

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2n+1

 

λπ

2

+...

 

λπ 2n

 

 

 

λπ

 

 

λπ

2

 

 

λπ 2n2

= cos x 1+

 

 

+

2

 

 

+sin x

2

 

1

+

 

+... +

 

 

;

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2n+2

= cos x

 

λπ

 

2

 

 

λπ 2n

 

 

 

λπ

 

 

λπ

2

 

 

λπ 2n

1

+

 

 

 

+... +

2

 

 

+sin x

2

 

1+

 

 

+... +

 

 

.

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Легко

видеть,

 

что

если

| λ | <

2

,

то

при

n → ∞ обе

подпоследовательности

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos x +

sin x

 

 

 

сходятся к

одной и

той

же функции

y(x)

=

 

 

 

2

 

 

 

,

которая уже была

 

 

πλ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

получена выше в качестве решения задачи другим способом.

Пример 4.8. Методом последовательных приближений решить уравнение Фредгольма

y(x) = λ1

y(s) ds +1.

Показать, что ряд Неймана сходится лишь в области | λ |<1, однако

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение, полученное при этом условии, существует при всех λ 1.

 

 

Решение.

Решение уравнения Фредгольма y = λAy + f может быть записано в виде ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неймана

y(x) = f (x) +λn An f (x) .

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

В нашем

 

случае

Af (x) = 1

f (x) dx , f (x) =1 ,

поэтому

Af (x) = 1

1 dx =1 и,

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

аналогично, An f (x) =1

(n = 2,3,...) . Следовательно, при условии | λ | <1 ,

ряд Неймана

 

 

 

 

 

 

 

 

 

 

 

 

 

λn An f (x) = λn

 

 

сходится,

 

и

решением

уравнения

является функция

n=1

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

y(x) =1+λn

=

 

 

 

,

однако

если

| λ | 1 , то

ряд расходится, и справедливость

1

λ

 

n=1

 

 

 

 

 

 

 

 

 

полученной формулы оказывается под сомнением.

34

 

С другой стороны, обозначив

1

y(s) ds = C ,

будем иметь

y(x) = λC +1, откуда

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C = 1

(λC +1) ds = λC +1.

Поэтому,

при λ =1 уравнение неразрешимо, а для остальных

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ 1

находим

C =

 

 

1

и

y(x) =1+λC =1+

 

 

λ

=

 

1

.

 

 

 

 

 

 

1

λ

1

λ

1λ

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, решение задачи существует при всех λ

1

и дается формулой

y(x) =

 

 

,

 

1

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что совпадает с выражением, полученным ранее лишь на множестве

| λ | <1 . Заметим, что

при условии

| λ | >1

решение

также

существует,

несмотря на

то, что

ряд Неймана

расходится.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задачи для самостоятельного решения

4.1Доказать, что интегральный оператор Фредгольма с вещественным непрерывным ядром,

действующий

в

пространстве

C[a,b] ,

 

не имеет характеристических чисел на интервале

 

 

 

 

1

 

 

K0

= max

 

K (x, s)

 

λ

0,

 

 

 

 

, где

 

.

 

 

 

 

 

 

K0

(b a)

 

x,s [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.2Доказать, что интегральный оператор Фредгольма с вещественным непрерывным ядром,

 

действующий

в пространстве

h[a,b] ,

 

не имеет характеристических чисел на интервале

 

λ

 

 

1

 

 

где K0

= max

 

K (x, s)

 

 

 

 

 

 

0,

 

 

 

 

,

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K0

(b

a)

 

x,s [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3

Доказать, что задача решения уравнения Фредгольма 2-рода y(x) = λb K (x, s) y(s) ds + f (x)

с

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

вещественным непрерывным ядром корректно поставлена при условии | λ | <

1

, где

 

K0 (b a)

 

K0 = max

 

K (x, s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x,s [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

в пространстве C[a,b] ;

 

 

 

 

 

 

 

 

 

б)

в пространстве h[a,b] .

 

 

 

 

 

 

 

 

4.4

Доказать,

что

оператор Фредгольма

y(x) = 1

1

x s2 y(s) ds +1 является

сжимающим

в

 

пространстве C[0,1] ,

 

 

2

0

 

 

 

 

 

и найти его неподвижную точку.

 

 

 

4.5Методом последовательных приближений построить резольвенту и получить решение уравнения Фредгольма 2-го рода:

а) y(x) = λ 1 x s y(s) ds + x ;

0

1

x

 

б) y(x) = λ 0

y(s) ds +1+ x2 ;

1+ s2

в) y(x) = λ1

y(s) ds +sin π x ;

0

 

 

35

г)

y(x) = λπx sin 2s y(s) ds +cos 2x ;

 

 

 

 

0

 

 

д) y(x) =λ

2π sin(x + s) y(s) ds + 2 ;

 

 

 

 

 

0

 

 

 

π

 

 

 

 

е)

y(x) =2 sin x cos s y(s) ds +1

(λ =1) ;

 

 

0

 

 

 

 

ж) y(x) = 1

1 x es y(s) ds +ex

(λ = 1 )

;

 

2

0

2

 

 

1

 

 

 

 

з) y(x) = 2 y(s) ds + x

(λ =1) .

 

 

0

 

 

 

 

В случаях а) - д) сформулировать критерий "малого" λ ; в примерах е) - з) проверить выполнение соответствующего условия.

Ответы к задачам

4.4

y =

4

x +1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

3xs

 

 

 

 

 

 

 

 

 

 

 

3x

 

 

 

 

 

 

 

 

 

 

 

4.5

а)

Резольвента

R(x, s, λ) =

 

,

 

решение

 

y(x) =

 

;

 

условие

 

сходимости

3 λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 λ

 

 

 

 

 

 

 

 

 

последовательных приближений

| λ | < 3 ,

достаточное условие теоремы выполнено при

 

| λ | <1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Резольвента

R(x, s, λ) =

 

 

 

1

 

 

 

 

x

 

,

решение

 

y(x) =1+

 

λx

 

+ x2 ;

 

 

 

 

 

 

 

ln 2

 

1+ s2

 

 

 

 

ln 2

 

 

 

 

 

 

 

1λ

 

 

 

 

 

 

 

 

 

 

 

1λ

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

условие сходимости последовательных

приближений

| λ | <

 

,

достаточное

условие

 

 

 

 

теоремы выполнено лишь при | λ | <1 .

 

 

 

 

 

 

 

 

 

ln 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в)

Резольвента

R(x, s, λ) =

 

 

 

1

 

 

,

 

решение

y(x) = sin πx +

 

2λ

 

;

 

условие

 

1

λ

 

 

π(1λ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сходимости последовательных приближений | λ | <1

совпадает с достаточным условием

 

теоремы.

 

 

 

 

x sin 2s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г)

 

Резольвента

R(x, s, λ) =

,

решение

y(x) = cos 2x ; последовательные

 

 

 

 

 

 

 

 

 

 

1+

λπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

приближения сходятся при | λ | < π2 , достаточное условие теоремы выполнено только для

| λ | < π12 .

36

д)

Резольвента

R(x, s, λ) = sin(x + s) +λπ cos(x s) ,

решение

y(x) = 2 ;

 

 

 

 

 

 

 

 

 

 

 

1λ2π2

 

1

 

 

 

 

 

 

 

последовательные приближения сходятся

при | λ | <

 

,

достаточное условие теоремы

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

выполнено только при | λ | <

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е)

Резольвента

R(x, s, λ) = 2sin x cos s

,

 

 

 

решение

y(x) =1+ 2sin x ;

 

 

 

 

 

 

 

 

 

 

2 λ

| λ | < 2 ,

 

 

 

 

 

 

последовательные приближения

 

сходятся при

однако достаточное

условие

"малого" λ : | λ | =

1 <

2

, сформулированное в теореме, -

не выполнено.

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ж)

Резольвента

R(x, s, λ) =

 

xes

 

,

решение

y(x) = x +ex ;

последовательные

1λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| λ | = 1

< 1

приближения сходятся при | λ | <1 ,

однако достаточное условие теоремы

не выполняется.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

e

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1 ;

 

 

 

 

з)

Резольвента

R(x, s, λ) =

 

 

,

решение

 

y(x) = x +

последовательные

 

2 λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

| λ | =1 < 2

 

 

приближения сходятся при | λ | < 2 ,

 

достаточное условие теоремы

также

выполняется.

37

ТЕМА 5

Линейное уравнение Вольтерра 2-го рода.

Основные определения и теоремы.

Уравнение

y = λ x

K (x, s) y(s) ds + f (x) ,

x, s [a, b] , или в операторной

форме

 

 

a

 

 

 

 

y = λ B y + f , называется уравнением Вольтерра 2-го рода.

 

 

Пусть ядро K (x, s)

непрерывно по совокупности переменных на своей треугольной

области определения

∆ ={x, s : a s x b}

и не равно нулю тождественно,

f (x) –

непрерывная на

[a,b] функция.

 

 

 

Теорема. Уравнение Вольтерра 2-го рода при любом значении λ

имеет

единственное решение для любой непрерывной функции

f (x) . Это решение может быть

найдено методом

последовательных приближений

yn+1 = λ B yn + f , y0 C[a,b] ,

f (x) C[a, b] .

 

 

 

 

 

 

Следствие 1. При любом λ однородное уравнение имеет только тривиальное решение.

Следствие 2. Оператор Вольтерра, действующий C[a,b] C[a,b] , не имеет

характеристических чисел. Таким образом, оператор Вольтерра является примером вполне непрерывного оператора, не имеющего ни одного характеристического числа.

Метод последовательных приближений для уравнения Вольтерра 2-го рода называется методом Пикара и выглядит так: для любого начального приближения

y0

C[a,b]

 

определим

 

 

yn+1 = λ x K (x, s) yn (s) ds + f (x), n = 0,1, 2,... ,

или

 

 

 

 

 

 

 

 

 

a

 

y

n

+1

= λ B y

n

+ f ,

причем

y

n

(x) y(x) .

 

 

 

 

 

 

n→∞

 

 

 

 

Полагая y0

= 0 , получаем ряд Неймана y = f +λ B f +λ2 B2 f +... +λn Bn f +... .

 

 

 

 

Если записать уравнение Вольтерра 2-го рода в операторной форме y = λAy + f

или

(I λA) y = f ,

то так как решение существует и единственно при любой непрерывной

функции f (x) и

любом λ , его

 

(решение)

можно

 

представить в виде

y = (I λA)1 f = (I +λRλ ) f = f +λRλ f , где

Rλ - интегральный оператор с непрерывным

по x, s ядром R(x, s, λ) , т.е. y(x) = f (x) +λx

R(x, s, λ) f (s)ds .

 

 

 

 

 

a

 

 

 

 

 

Ядро R(x, s, λ)

оператора Rλ

называется резольвентой.

 

 

Функции

K1 (x, s) K (x, s) ,

Kn (x, s) = x K (x,t) Kn1 (t, s) dt ,

n = 2,3,... называются

 

 

 

 

 

s

 

 

 

повторными

(итерированными)

ядрами.

 

Ряд K (x, s) +λK (x, s) +... +λn1 K (x, s) +…

 

 

 

 

 

 

2

n

 

 

 

 

 

1

 

=K ( x,s)

сходится равномерно по x, s [a, b] при любых λ, в отличие от аналогичного ряда для

38

уравнения Фредгольма, сходимость которого гарантировалась лишь при

 

λ

 

<

1

, и

 

 

 

 

K0 (b a)

 

 

 

 

 

 

резольвента R(x, s, λ) может быть получена по формуле R(x, s, λ) = λm1Km (x, s) .

m=1

Уравнения Вольтерра с ядрами специального вида могут также решаться путем сведения к дифференциальному уравнению, либо с использованием преобразования Лапласа (примеры 5.4 и 5.5).

Примеры решения задач.

Пример

5.1.

Методом последовательных приближений

построить резольвенту

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x2

 

интегрального

уравнения

Вольтерра

2-го

рода

 

y(x) = λe( xs) y(s) ds + xe

 

 

и найти

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

решение этого уравнения при λ =1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

Вычислим повторные ядра этого уравнения: K (x, s) = K (x, s) = e( xs) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

K2 (x, s) = x

K (x,t) K1 (t, s) dt = x e( xt ) e(ts) dt = (x s) e( xs) ,

 

 

 

 

 

 

 

 

s

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K3 (x, s) = x

K (x,t) K2 (t, s) dt = x e( xt ) (t s) e(ts) dt =

(x s)2

e( xs) ,

 

2

 

 

s

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

………….

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Km+1 (x, s) = x K (x,t) Km (t, s) dt = (x s)m e( xs) .

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

m!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

эти

выражения

в

формулу

для

резольвенты,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R(x, s, λ) = λm1Km (x, s) = e( xs) eλ( xs) = e(λ1)( xs) . Заметим, что ряд сходится при любых

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ, что обеспечивается не

малостью λ, как было

в случае уравнения Фредгольма, а

наличием множителя

m!

в знаменателях повторных ядер.

 

 

 

 

 

 

 

Далее,

положив λ =1 , получим R(x, s, λ) = R(x, s,1) =1 и запишем

решение

уравнения при

f (x) = xe

x2

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x2

x

 

s2

 

 

x2

 

y(x) = f (x) +λ R(x, s, λ) f (s) ds = xe

 

+1 s e

 

ds = (x +1) e

 

1 .

 

2

2

2

 

 

 

 

 

a

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Пример 5.2. Методом последовательных приближений решить уравнение Вольтерра y(x) = x (s x) y(s) ds +1 .

0

Решение. Итерационный

процесс для данного уравнения выглядит так:

yn+1 (x) = x (s x) yn (s) ds +1 .

Выберем в качестве начального приближения y0 (x) 0 , тогда

0

 

последовательно найдем:

 

39

y1 (x) = x (s x) 0 ds +1 =1 ,

0

y2 (x) = x (s x) y1 (s) ds +1 = x (s x) 1 ds +1 =1

x2

,

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

y3 (x) = x (s x) y2

(s) ds +1 =1

x2

+

x4

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2!

 

4!

 

 

 

 

 

 

 

 

 

 

……………..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn+1 (x) = x (s x) yn (s) ds +1 =1

x2

+

x4

... + (1)n x2n .

 

 

 

 

 

 

 

 

0

2!

4!

 

 

(2n)!

 

 

 

 

Переходя к пределу при n → ∞, получим функцию y(x) = lim yn (x) = cos x , которая и

является решением рассматриваемого уравнения.

 

n→∞

 

 

 

 

 

 

Пример 5.3. Доказать, что если ядро уравнения Вольтерра

y(x) = λx

K (x, s) y(s) ds + f (x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

зависит только от разности аргументов,

т.е.

K (x, s) = K (x s) , то все повторные ядра, а

следовательно и резольвента, также являются функциями лишь от разности

(x s) .

Решение. Пусть K (x, s) = K (x s) , тогда

K2 (x, s) = x

K (x t) K (t s) dt . Произведя замену

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

переменной

интегрирования

по

 

 

формуле

t s =ξ ,

получим

K2 (x, s) = xs K (x s ξ) K (ξ) dξ = xs F2 (x s,ξ) dξ = Φ2 (x s) .

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

Действуя

далее

 

 

 

 

 

аналогичным

путем,

найдем,

Km (x, s) = x K (x,t) Km1 (t, s) dt = x K (x t) Km1 (t s) dt = xs K (x s ξ) Km1 (ξ) dξ = Φm (x s) ,

s s 0

что и требовалось доказать.

Пример 5.4. С помощью преобразования Лапласа найти резольвенту и записать решение

уравнения Вольтерра y(x) = x sin(x s) y(s) ds + f (x) .

0

Решение. Поставим в соответствие функциям, входящим в уравнение, их изображения:

y(x) ÷Y ( p) , f (x) ÷ F( p) ,

sin x ÷

 

1

 

K ( p) . Учитывая, что изображение свертки двух

 

p2 +1

 

 

 

 

 

 

 

 

 

 

 

функций есть произведение их изображений, получим

Y ( p) = K ( p) Y ( p) + F( p) , откуда

Y ( p) =

 

 

F( p)

= F( p) +

 

 

K ( p)

F( p) .

 

1

K ( p)

1

K ( p)

 

 

 

 

 

 

 

Используя результат предыдущей задачи, можно сделать вывод, что резольвента

зависит

лишь

от разности (x s) ,

следовательно,

решение исходного уравнения

представимо в виде y(x) = f (x) +x R(x s,1) f (s) ds . Сравнивая последние две формулы,

a

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]