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

Сумин Интегралные уравнения Фредголма и Волтерра 2016

.pdf
Скачиваний:
7
Добавлен:
12.11.2022
Размер:
1.12 Mб
Скачать

Все эти решения определяются формулой

y(x) = Cλ cos2 x = 2πC cos2 x = C1 cos2 x .

Пример 7.3. Исследовать на разрешимость при различных значениях параметра λ интегральное уравнение

 

y(x) −λ 1

(2xt 4x2 ) y(t) dt =1 2x .

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Решение.

Ядро K (x, t) = 2xt 4x2

является вырожденным, где

p (x) = 2x , q (t) =t , p

2

(x) = −4x2 , q (t) =1.

 

1

1

 

 

 

 

 

 

2

 

 

 

Запишем интегральное уравнение в следующем виде:

y(x) 2x 1

ty(t) dt −λ 4x2 1

y(t) dt +1 2x .

 

 

 

 

 

0

 

 

 

 

0

 

 

Обозначим C1 = 1

ty(t) dt , C2

= 1

y(t) dt . Тогда

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

 

0

 

 

 

 

 

 

0

 

 

 

уравнение запишется следующим образом:

 

 

y(x) = C λ 2x C

2

λ 4x2 +1 2x .

(7.11)

 

 

 

 

 

 

1

 

 

 

 

 

Подставим (7.11) в выражения для C1 и C2 и получим систему

уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

C1 =

t(C1 2λt C2λ 4t2 +1 2t) dt;

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

C2 =

 

(C1λ 2t C2λ 4t2 +12t) dt,

 

 

 

0

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

C1

1

 

λ +C2

λ = −

 

;

 

 

 

3

6

 

 

 

 

 

 

 

 

 

(7.12)

 

 

 

 

 

 

 

1

+ 4

λ = 0.

 

 

 

C λ + C

 

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим главный определитель

 

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

тели C и

C

2

системы уравнений (7.12):

 

 

1

 

 

 

 

 

 

 

 

 

 

41

 

1

2

λ

λ

 

 

 

 

 

 

 

2

 

 

 

4

 

 

 

 

 

 

 

λ

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

4

 

 

 

= 1

 

 

λ 1

+

 

 

 

λ

 

=

 

+1

,

 

 

 

 

 

 

 

3

3

 

3

 

−λ 1

+

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

λ

 

 

 

 

 

 

 

1

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 =

 

 

 

4

 

 

 

= −

 

1 +

 

 

λ

,

 

 

 

 

 

 

 

 

 

 

 

 

6

3

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

+

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

1

2

λ

1

 

 

= −

1

λ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2

 

 

3

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−λ

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если 0 , то по формулам Крамера

 

 

 

 

 

 

 

1

 

 

+

4

 

 

 

 

 

 

 

 

 

 

6

1

3

λ

 

C

=

 

C1

=

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

3

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

C2 =

 

 

C2

 

=

 

 

 

 

6

 

 

.

 

 

 

 

 

 

 

 

λ

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Если λ ≠ −3 , то данное неоднородное интегральное уравнение имеет единственное решение

 

2λx

1

 

+

4

 

 

 

λ

 

 

 

1

 

λ 4x2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

y(x) =

 

 

6

 

3

 

+

 

 

 

6

+12x =

 

λ

 

2

 

 

 

λ

+

2

 

 

 

 

 

 

 

 

 

 

+1

 

 

 

 

3

1

 

 

3

 

 

 

 

 

 

 

 

 

= 3x(2λ2 x 2λ2 5λ −6) +(λ +3)2

(λ +3)2 ,

а соответствующее однородное уравнение

42

y(x) −λ 1 (2xt 4x2 ) y(t) dt = 0

0

имеет только нулевое решение (λ не является характеристическим числом, 1-й случай альтернативы Фредгольма).

2. Если λ = −3 , то система уравнений (7.12) принимает вид

 

 

 

1

 

 

 

 

 

1

 

3C1

3C2

= −

 

;

или

C1

C2

= −

 

;

6

18

 

 

 

 

 

 

 

 

 

3C2

= 0,

 

 

 

 

C2

= 0.

 

3C1

 

 

 

C1

 

Данная система уравнений несовместна, т.е. решений не имеет. Отметим, что исходное неоднородное интегральное уравнение также не имеет решений, так как правая часть f (x) =12x не ор-

тогональна к p (x) = 2x

и p

2

(x) = −4x2 , поскольку

 

1

 

 

 

1

(12x) 2x dx 0 и 1

(12x) (4x2 ) dx 0 .

0

 

 

 

0

 

Соответствующее однородное уравнение имеет бесконечное множество решений (λ – характеристическое число, 2-й случай альтернативы Фредгольма), так как однородная система уравнений для определения C1 и C2 , соответствующая (7.12), принимает вид

3C1 3C2 = 0;

3C1 3C2 = 0,

откуда C1 =C2 =C , где C произвольно.

Все эти решения согласно (7.11) определяются формулой y(x) = −6Cx +12Cx2 +12x = C3 x 2C3 x2 +12x .

Пример 7.4. Исследовать на разрешимость при различных значениях параметра λ интегральное уравнение

y(x) −λ π cos (x +t) y(t) dt = cos 3x .

0

Решение. Отметим, что ядро

K (x, t) = cos (x +t) =cos x cost sin xsin t

43

является одновременно и

симметричным и вырожденным, где

p1 (x) =cos x , q1 (t) =cost ,

p2 (x) = −sin x ,

q2 (t) =sin t .

Запишем интегральное уравнение в следующем виде:

y(x) = λcos x π

cost y(t) dt −λsin x π sin t y(t) dt +cos 3x .

0

 

0

 

Обозначим

 

 

 

C1 = π

cost y(t) dt , C2 = π

sin t y(t) dt .

0

 

0

 

Тогда интегральное уравнение запишется следующим образом:

y(x) =C1λcos x C2λsin x +cos3x . (7.13)

Подставим (7.13) в выражения для C1 и C2 и получим систему уравнений:

 

 

π

 

 

 

 

 

 

 

 

 

 

C1

= cost (C1λcost C2λsin t + cos3t) dt;

 

0

 

 

 

 

 

 

 

 

 

 

 

= π

 

 

 

 

 

 

 

 

 

 

C2

sin t (C1λcost C2λsin t +cos 3t) dt,

 

0

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

= 0;

 

 

 

 

 

 

 

C1 1 −λ

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

(7.14)

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

C2 1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим определитель

системы уравнений (7.14):

 

 

1 −λ

π

0

 

 

 

 

π

π

 

 

 

 

 

=

 

2

 

 

 

 

 

 

 

 

π

 

= 1−λ

 

1+ λ

.

 

 

 

0

 

1+ λ

 

 

 

2

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

1. Если λ ≠ ± π2 , то система уравнений (7.14) имеет единствен-

ное решение: C1 =0 , C2 = 0 . Поэтому исходное неоднородное интегральное уравнение имеет единственное решение y(x) =cos3x , а соответствующее однородное уравнение

44

y(x) −λ π cos (x +t) y(t) dt = 0

0

имеет только нулевое решение (λ не является характеристическим числом, 1-й случай альтернативы Фредгольма).

2. Если λ = λ1 = π2 , то система уравнений (7.14) принимает вид

C1 0 = 0;

C2 0 = 0,

откуда C1 произвольно, C2 = 0 .

Исходное неоднородное интегральное уравнение имеет бесконечное множество решений, которые определяются формулой:

y(x) = C1λ cos x + cos3x = C1 π2 cos x +cos3x = C cos x +cos3x .

Отметим, что

правая часть

f (x) = cos3x ортогональна к

p1 (x) =cos x и p2 (x) = −sin x , так как

π

cos3t cost dt = 0 ,

π

cos 3t sin t dt = 0 .

0

 

0

 

Соответствующее однородное интегральное уравнение имеет бесконечное множество решений, которые определяются формулой:

y(x) = C1λ cos x = C1 π2 cos x = C3 cos x

( λ1 является характеристическим числом, 2-й случай альтернативы Фредгольма).

3. Если λ = λ2 = − π2 , то система уравнений (7.14) принимает вид

C1 2 = 0;

C2 0 = 0,

откуда C1 =0 , C2 произвольно.

45

Исходное неоднородное интегральное уравнение имеет бесконечное множество решений, которые определяются по формуле:

y(x) = −C2λ sin x + cos3x = C2 π2 sin x +cos 3x = C sin x +cos 3x .

Соответствующее однородное интегральное уравнение имеет бесконечное множество решений, которые определяются формулой:

y(x) = C2λ sin x = C2 π2 sin x = C4 sin x

( λ2 является характеристическим числом, 2-й случай альтернативы Фредгольма).

46

Глава 2. Интегральные уравнения Вольтерра

§ 8. Уравнения Вольтерра. Основные понятия и определения

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

зывается уравнение вида

y(x) −λ x

K (x, t) y(t) dt = f (x) .

(8.1)

a

 

 

Интегральное уравнение вида

x

K(x, t) y(t) dt = f (x)

(8.2)

a

 

 

называется интегральным уравнением Вольтерра 1-го рода. Все обозначения в уравнениях (8.1), (8.2) аналогичны обозначениям в уравнениях Фредгольма (1.1), (1.2).

Формально уравнения Вольтерра (8.1), (8.2) можно рассматривать как частный случай уравнений Фредгольма (1.1), (1.2), полагая в уравнениях Фредгольма K (x, t) 0 при t > x или рассматривая

ядро интегрального уравнения Вольтерра в треугольнике a x b , a t x . Однако физические задачи, приводящие к уравнениям Вольтерра и Фредгольма, а также свойства решений этих уравнений существенно различны. Поэтому уравнения Вольтерра выделяют в отдельный тип интегральных уравнений.

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

Уравнения Вольтерра в отличие от уравнений Фредгольма всегда имеют единственное решение. В частности, уравнение Вольтерра 2-го рода (8.1) имеет единственное решение при любом значении числового параметра λ. Это решение может быть найдено

47

методом последовательных приближений. Выбирая произвольным образом нулевое приближение y0 (x) , можно с помощью рекуррентной формулы

yn (x) = f (x) + λ x

K (x, t) yn1 (t) dt

(9.1)

a

 

 

построить последовательность функций {yn (x)}, которая всегда сходится к единственному решению интегрального уравнения при

n → ∞ , т.е. y(x) = lim yn (x) .

n→∞

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

 

 

 

 

 

 

 

 

y(x) = x

y(t) dt +

x2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Отметим, что в данном уравнении λ =1,

 

 

K (x, t) =1 .

В качестве нулевого приближения возьмем

 

y0 (x) =1 . Строим

последовательность функций {yn (x)} по формуле (9.1):

 

y1 (x) = f (x) + x

 

K (x, t) y0 (t) dt =

 

x2

 

+

x

1 dt =

x2

+ x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

y2 (x) = f (x) + x

 

K (x, t) y1 (t) dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

x

 

2

 

 

 

 

 

 

 

 

x

2

 

 

x

3

 

 

 

 

x

2

 

 

 

 

 

 

 

 

 

x

3

 

 

 

 

 

 

=

 

+

 

t

 

 

 

+ t dt =

 

 

 

+

 

 

 

 

+

 

 

 

= x2 +

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

0

2

 

 

 

 

 

 

 

2

 

 

 

6

 

 

 

 

 

2

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

x

 

 

2

 

 

 

t3

 

 

 

 

y3 (x) = f (x) + K (x, t) y2

(t) dt =

 

 

 

 

+

 

t

 

+

 

 

 

 

 

 

dt =

 

 

2

 

 

 

6

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

=

x2

 

 

+

x3

 

+

 

x4

 

=

 

x2

+

 

 

x3

 

+

 

x

4

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

6

 

24

 

 

2!

 

 

3!

 

4!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

…,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn (x) = f (x) + x

K (x, t) yn1 (t) dt =

x2

 

+

x3

+

x4

 

+... +

xn

 

=

 

 

4!

n!

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

 

3!

 

 

 

 

 

 

 

 

 

 

 

 

 

48

 

 

x

2

 

x

3

 

x

4

 

x

n

 

n

x

k

 

= 1

+ x +

 

+

 

+

 

+... +

 

 

(x +1)

=

 

x 1.

 

 

 

 

 

 

 

 

k!

 

2! 3!

4!

 

n!

 

k =0

 

Отсюда точное решение y(x) определяется как предел

 

 

 

n

xk

 

 

x

 

y(x) = lim

yn

(x) = lim

 

x 1

= e

 

x 1 .

k!

 

n→∞

 

n→∞ k =0

 

 

 

 

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

 

y(x) = x +1 x

y(t) dt .

 

0

 

Отметим, что в данном уравнении λ =1, K (x, t) =1 .

В

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

f (x)

данного уравнения, т.е. y0 (x) = f (x) = x +1 .

Строим по формуле (9.1) последовательность функций {yn (x)}:

y1 (x) = f (x) + x

K (x, t) y0 (t) dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= x +1 x

(t +1) dt = x +1

x2

x =1

x2

 

 

,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

t

2

 

 

 

y2 (x) = f (x) + K (x, t) y1 (t) dt = x +11

 

 

 

 

dt =

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

2

 

 

 

 

= x +1 x +

 

x3

=1 +

 

x3

,

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y3 (x) = f (x) + x

K(x, t) y2 (t) dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

t

3

 

 

 

 

 

 

 

x

4

 

 

 

 

x

4

 

 

= x +1

1+

 

dt = x +1x

 

 

=1

 

 

 

,

 

 

 

 

 

 

 

 

 

0

 

6

 

 

 

 

 

 

 

 

24

 

 

24

 

…,

49

xn+1

x

yn (x) = f (x) + 0 K (x, t) yn1 (t) dt =1+ (1)n (n +1)!.

Отсюда точное решение y(x) определяется как предел

y(x) = lim

y

 

(x) = lim

1+ (1)n

xn+1

=1 .

n

 

 

 

 

 

 

 

 

n→∞

 

 

n→∞

(n +1)!

 

§10. Решение интегрального уравнения путем сведения его

кдифференциальному уравнению

Если в интегральных уравнениях (8.1), (8.2) ядро K (x, t) и свободный член f (x) имеют непрерывные производные по перемен-

ной x, то эти уравнения могут быть почленно продифференцированы несколько раз. Это позволяет свести решение интегрального уравнения к задаче Коши для некоторого обыкновенного дифференциального уравнения. Отметим, что производная от интеграла вычисляется по формуле

d

x

K (x, t) y(t) dt = K (x, x) y(x) + x

K (x, t)

y(t) dt . (10.1)

dx

x

0

0

 

Пример 10.1. Решить интегральное уравнение y(x) = x + x xty(t) dt.

0

Решение. Данное уравнение перепишем в следующем виде

 

x

 

y(x) = x 1

+ ty(t) dt .

 

0

 

Обозначим

z(x) =1+ x ty(t) dt .

0

Продифференцируем последнее уравнение и найдем z(x) = xy(x) .

Так как y(x) = x z(x) , то получаем обыкновенное дифференциальное уравнение z(x) = x2 z(x) . Отсюда z(x) = Cex3 /3 .

50