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

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

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

§ 5. Собственные значения и собственные функции

Рассмотрим однородное уравнение Фредгольма 2-го рода

y(x) −λ b

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

(5.1)

a

 

 

Отметим, что уравнение (5.1) всегда имеет очевидное решение y(x) 0 , которое называется нулевым (тривиальным) решением.

Значения числового параметра λ, при которых это уравнение имеет ненулевые решения y(x) 0, называются характеристическими

числами (величина μ = λ1 называется собственным значением)

уравнения (5.1) или ядра K (x, t) . Каждое ненулевое решение этого уравнения называется собственной функцией, соответствующей характеристическому числу λ (собственному значению μ).

Подчеркнем, что число λ = 0 не является характеристическим числом, так как при λ = 0 из уравнения (5.1) следует, что y(x) 0 .

Если ядро K (x, t) однородного интегрального уравнения Фред-

гольма 2-го рода (5.1) является вырожденным, то задача о нахождении собственных значений и собственных функций интегрального уравнения сводится к поиску собственных значений некоторой матрицы. В самом деле, как следует из формул (4.2), (4.3), (4.6), (4.7) (здесь f (x) = 0 ), всякое решение однородного интегрального

уравнения (5.1) имеет вид

n

 

y(x) = λ ck pk (x) ,

(5.2)

k =1

где неизвестные числа ck являются решением однородной системы уравнений

n

 

 

 

ck −λ akici

= 0

(k = 1, 2, …, n).

(5.3)

i=1

 

 

 

Система (5.3) может быть записана в матричной форме

 

(I −λA)C = 0

или

(A −μI)C = 0 ,

(5.4)

21

где λ, μ ≠ 0 ; μ = λ1 ; I = (δij )nn – единичная матрица порядка n;

A = (aij )nn – квадратная матрица порядка n; С – матрица-столбец,

состоящая из чисел ci ( i =1, n ); 0 – нулевая матрица-столбец.

Таким образом, собственные значения однородного интегрального уравнения (5.1) совпадают с отличными от нуля собственными значениями матрицы A и могут быть найдены из характеристического уравнения

det (A −μI) = 0 .

(5.5)

Отметим, что если исходное интегральное уравнение Фредгольма 2-го рода является неоднородным

y(x) −λ b

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

a

 

и имеет вырожденное ядро K (x, t) , то его решение можно свести к

решению системы линейных алгебраических уравнений (4.7), которая может быть записана в матричной форме

(I −λA)C = B ,

(5.6)

где B – матрица-столбец, состоящая из чисел bi ( i =1, n ).

Пример 5.1. Найти собственные значения и собственные функции интегрального уравнения

y(x) −λ π

(cos2 x cos 2t +cos3x cos3 t) y(t) dt = 0 .

0

 

 

 

 

Решение. Ядро

 

 

 

 

K (x, t) = cos2 x cos 2t +cos 3x cos3 t

 

 

является вырожденным. Здесь

 

 

 

p (x) = cos2 x ,

q (t) =cos 2t ,

p (x) = cos3x , q

(t) = cos3 t .

1

1

2

2

 

По формулам (4.6) найдем элементы матрицы A:

 

 

a11 = π

p1 (t) q1 (t) dt = π

cos2 t cos 2t dt =

π ,

0

0

 

4

22

 

a12 = π

p2 (t) q1 (t) dt = π

cos3t cos 2t dt = 0,

 

0

 

0

 

 

a21 = π

p1 (t) q2 (t) dt = π

cos2 t cos3 t dt = π

cos5 t dt = 0 ,

0

 

0

 

0

 

a22 = π

p2 (t) q2 (t) dt = π

cos3t cos3 t dt = π .

 

0

 

0

 

8

Характеристическое уравнение (5.5) для нахождения собственных значений имеет вид

 

 

 

 

 

 

 

 

 

π −μ

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det (A −μI) =

4

π

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

0

−μ

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

откуда

 

π

 

π

 

Получаем

собственные значения

 

 

 

−μ

−μ = 0 .

 

 

 

4

 

8

 

 

 

 

 

 

 

μ =

π

,

μ

 

=

π

(соответственно характеристические числа λ =

4

,

4

 

8

π

1

 

 

2

 

 

 

 

 

1

 

λ2 = π8 ).

1.При μ1 = π4 система уравнений (5.4) принимает вид

 

 

 

0

0

 

 

с

 

,

 

 

 

 

 

 

 

 

π

 

 

 

1

= 0

 

 

 

 

 

 

0

 

 

с2

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Собственная функция, соответст-

вующая собственному значению

μ = π ,

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

 

 

 

 

 

 

 

 

 

1

4

 

 

 

(5.2). Именно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) c p (x) + λ c p

(x) = λ c cos2

x =

1

1

1

 

1

2

 

2

 

 

1

1

 

=

4

c cos2

x = C cos2

x ,

C 0 .

 

 

 

 

π

1

 

 

 

1

 

 

 

 

 

1

 

23

2. При μ2 = π8 система уравнений (5.4) принимает вид

 

 

 

π

0

 

 

с

 

,

 

 

 

 

 

8

 

 

 

1

= 0

 

 

 

 

 

 

 

 

с2

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

Собственная функция, соответст-

вующая собственному значению

 

μ2

= π

, находится по формуле

(5.2). Именно,

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) 2c1 p1 (x) 2c2 p2 (x) = λ2c2 cos3x =

=

8

c cos3x = C

 

cos3x ,

C

 

0 .

 

 

 

 

π 2

 

 

 

2

 

 

 

 

2

 

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

y(x) π

cos (x +t) y(t) dt +1.

0

 

Решение. Данное интегральное уравнение Фредгольма является неоднородным. Запишем его в следующем виде

y(x) π

(cos x cos t sin x sin t) y(t) dt +1.

0

 

Ядро K (x, t) = cos x cost sin x sin t является вырожденным. Здесь

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

(x) = −sin x , q2 (t) =sin t . По формулам

(4.6) найдем элементы матрицы A и матрицы-столбца B:

 

a11 = π

p1 (t) q1 (t) dt = π

cos2 t dt = π

1+ cos 2t

dt =

π ,

 

 

0

 

 

0

 

0

 

2

 

2

a12 = π

p2

(t) q1 (t) dt = −π

sin t cost dt = −π

sin t d (sin t) = 0 ,

0

 

 

0

 

 

0

 

 

 

 

a21 = π

p1 (t) q2 (t) dt = π

cos t sin t dt = 0 ,

 

 

 

0

 

 

0

 

 

 

 

 

24

a22 = π

p2 (t) q2 (t) dt = −π

sin2 t dt = −π

1 cos 2t

dt = −

π

,

 

2

0

 

0

 

0

 

2

 

 

b1 = π

q1 (t) f (t) dt = π

cos t 1 dt = π

cost dt = 0 ,

 

 

 

0

 

0

0

 

 

 

 

b2 = π

q2 (t) f (t) dt = π

sin t 1 dt = π

sin t dt = 2 .

 

 

 

0

 

0

0

 

 

 

 

Характеристическое уравнение для определения характеристических чисел вытекает из (5.4) и имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −λ

π

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det (I −λA) =

 

 

2

 

 

 

 

π

= 0 ,

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1+ λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

откуда

 

−λ

π

 

 

π

=0 .

Получаем характеристические числа

1

2

1 + λ

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

λ = −

2

, λ

 

=

2

. Система уравнений (5.6) принимает вид

π

 

 

 

1

 

 

2

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−λ

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

0

 

 

 

с

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

1

=

 

.

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

с2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. При

λ ≠ ±

2

 

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

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

=0 ,

c =

 

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

В этом случае решение исходного интегрального уравнения согласно формуле (4.4) имеет вид

y(x) = λc p (x) + λc p

 

(x) + f (x) = −

 

2λ

 

sin x +1 =

2

 

 

π

1

1

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

25

=1 2 +4λλπ sin x .

2.При λ = π2 решения системы уравнений

0

0

c

 

 

0

 

0

2

1

 

=

2

 

 

c2

 

образуют пары (c1 , c2 ) , где c1 произвольно, c2 =1. Соответственно, решения интегрального уравнения имеют вид

y(x) = λc1 p1 (x) + λc2 p2 (x) + f (x) = π2 c1 cos x π2 1 sin x +1 =

=C cos x π2 sin x +1.

3.При λ = − π2 система уравнений

 

2

0

c

 

 

0

 

0

0

1

 

=

2

 

 

c2

 

решений не имеет. Значит, при данном λ не имеет решений и исходное интегральное уравнение.

Отметим, что однородное интегральное уравнение Фредгольма 2-го рода может вообще не иметь характеристических чисел и собственных функций, или не иметь действительных характеристических чисел и собственных функций.

Пример 5.3. Показать, что интегральное уравнение

y(x) = λ 1

( x t t x)y(t) dt

0

 

не имеет действительных характеристических чисел и собственных функций.

Решение. Ядро K (x, t) =

x t t x

является

вырожденным.

Здесь p1 (x) = x , q1 (t) =t ,

p2 (x) = −x ,

q2 (t) = t .

По формулам

(4.6) найдем элементы матрицы A:

 

 

26

a11 = 1

p1 (t) q1 (t) dt = 1

t t dt = 1

t3/ 2 dt =

2

,

 

 

5

 

 

 

0

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

a12 = 1

p2 (t) q1 (t) dt = −1

t2 dt = −

1

,

 

 

 

 

 

3

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

a21 = 1

p1 (t) q2 (t) dt = 1

t dt =

1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

2

 

 

 

 

 

 

 

a22 = 1

p2 (t) q2 (t) dt = −1

t t dt = −1

t3/ 2 dt = −

2

.

5

0

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

Характеристическое уравнение для определения характеристических чисел вытекает из (5.4) и имеет вид

 

 

 

 

 

 

1 −λ

2

 

 

 

1

λ

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

det (I −λA) =

 

 

 

1

 

 

 

 

 

3

 

2

 

=

 

 

 

 

λ

 

 

 

1

+ λ

 

 

 

 

 

 

 

 

 

2

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

1

 

2

 

 

 

 

λ2

.

= 1

 

λ 1

+

 

 

λ

+

 

 

 

λ

 

=1 +

 

 

 

 

5

5

6

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

При действительных λ выполнено соотношение det (I −λA) 0 ,

поэтому при всех действительных λ исходное интегральное уравнение имеет только тривиальное решение y(x) 0 . Итак, данное

уравнение не имеет действительных характеристических чисел и собственных функций.

Пример 5.4. Показать, что интегральное уравнение

y(x) = λ 1

sin πx cos πt y(t) dt

0

 

не имеет характеристических чисел и собственных функций. Решение. Запишем исходное уравнение в виде

y(x) = λsin πx 1

cos πt y(t) dt

(5.7)

0

 

 

27

и обозначим

C = 1 cos πt y(t) dt .

0

Тогда

y(x) =Cλsin πx .

Подставив полученное выражение для y(x) в обе части уравнения (5.7), получим

Cλsin πx = Cλ2 sin πx 1 cos πt sin πt dt .

 

 

 

 

0

Но

 

 

 

 

1

cos πt sin πt dt =

1

π

cos z sin z dz = 0 ,

π

0

 

0

 

 

 

 

поэтому Cλsin πx = 0 , откуда C = 0 , и значит y(x) 0 . Таким об-

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

§ 6. Интегральные уравнения с симметричным ядром

Ядро интегрального уравнения называется симметричным, если K (x, t) = K (t, x) , где аргументы x и t определены на основном

квадрате [a, b]×[a, b] . Одним из методов решения однородного

интегрального уравнения Фредгольма 2-го рода с симметричным ядром является сведение этого уравнения к краевой задаче.

Пример 6.1. Найти характеристические числа и собственные функции однородного уравнения

y(x) = λ π K (x, t) y(t) dt ,

0

где симметричное ядро определяется следующей формулой

cos x sin t, 0 x t;

K (x, t) =

cost sin x, t x ≤ π.

28

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

 

y(x) = λ x

K (x, t) y(t) dt + λ π

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

 

 

0

 

x

 

 

 

или

 

 

 

 

 

 

y(x) = λsin x x

cost y(t) dt + λcos x π

sin t y(x) dt .

(6.1)

 

 

0

 

x

 

 

Дифференцируя уравнение (6.1), получаем

 

 

 

 

x

 

π

 

 

(x) = λcos x cost y(t) dt −λsin x sin t y(t) dt .

(6.2)

y

 

 

0

 

x

 

 

Повторное дифферецирование дает

′′

 

x

y (x) = λy(x) − λsin x cos t y(t) dt

 

0

+ λcos x π

sin t y(t) dt .

x

 

Но выражение в квадратных скобках равно y(x) , поэтому последнее уравнение можно записать следующим образом:

y′′(x) y(x) y(x) .

Из равенств (6.1) и (6.2) получаем, что y(π) = 0 , y(0) =0 .

Итак, исходное интегральное уравнение сведено к следующей краевой задаче

y′′(x) (λ −1) y(x) = 0;

(6.3)

 

 

(6.4)

y(π) = 0,

y (0) = 0.

Здесь возможны три случая.

 

 

 

Случай 1. Пусть λ −1 = 0 , т.е.

λ =1. Уравнение (6.3) принимает

′′

 

+C2 . Используя краевые условия

вид y (x) =0 , откуда y(x) =C1 x

(6.4), получаем систему уравнений

 

C π+ C = 0;

 

 

1

2

 

C1 = 0,

 

которая имеет единственное решение: C1 =0 , C2 = 0 . Значит, интегральное уравнение имеет только тривиальное решение y(x) 0 .

29

Случай 2. Пусть λ −1 > 0 , т.е. λ >1. Общее решение уравнения (6.3) имеет вид

y(x) =C1 ch ( λ −1 x)+C2 sh ( λ−1 x),

откуда

 

 

 

 

sh ( λ −1 x)+ C2 ch ( λ −1 x) .

y (x) = λ −1 C1

Краевые условия (6.4) дают следующую систему уравнений

C1 ch (π λ −1)+C2 sh (π λ −1)= 0;

 

 

 

 

 

λ −1 C2

1

= 0.

 

Откуда следует, что C2 = 0 ,

C1 =0 (так как ch (π λ −1)0 ). Ин-

тегральное уравнение имеет только тривиальное решение y(x) 0 .

Итак, при λ ≥1 интегральное уравнение не имеет характеристических чисел, а значит, и собственных функций.

Случай 3. Пусть λ −1 < 0 , т.е. λ <1. Общее решение уравнения (6.3) имеет вид

y(x) = C1 cos ( 1−λ x)+C2 sin (

1−λ x),

(6.5)

откуда

 

 

 

 

 

 

1 −λ −C1 sin ( 1 −λ x)+C2

cos ( 1 −λ x) .

 

y (x) =

 

Краевые условия (6.4) дают следующую систему уравнений

 

C1 cos (π 1 −λ)+C2 sin (π 1 −λ)= 0;

(6.6)

 

 

 

 

1−λ C2 1 = 0.

 

 

 

 

 

Из-за того, что система уравнений (6.6) – однородная, она будет иметь ненулевое решение, только если ее основной определитель равен нулю, т.е.

cos (π 1−λ)

sin (π 1− λ)

 

= 1 −λ cos (π 1−λ)= 0 .

 

0

1−λ

 

 

30