МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ Интегральные уравнения

advertisement
Федеральное агентство по образованию Российской Федерации
Государственное образовательное учреждение высшего профессионального образования
Ивановский государственный химико-технологический университет
МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Интегральные уравнения
Методические указания
Составитель: Г.А. Зуева
Иваново 2006
Составитель: Г.А. Зуева
УДК 613.19
Методы математической физики. Интегральные уравнения: Методические указания / Иван. гос. хим.-технол. ун-т; Сост. Г.А. Зуева. - Иваново,
2006. – 32 с.
Методические указания содержат изложенный в краткой форме теоретический материал, относящийся к уравнениям математической физики (интегральным уравнениям). Приведены контрольные задачи и упражнения с ответами, представлен разбор решения типовых задач.
Методические указания рекомендуются студентам, изучающим методы математической физики, а также аспирантам высших технических учебных
заведений.
Рецензент
доктор химических наук, профессор В.В. Рыбкин
(ГОУ ВПО «Ивановский государственный химико-технологический университет»)
2
Вопросы и задачи для самопроверки
1.
2.
3.
4.
5.
Что такое неопределенный интеграл?
Напишите формулу интегрирования по частям.
Напишите формулу Ньютона-Лейбница.
В чем состоит геометрический смысл определенного интеграла?
Чему равна производная определенного интеграла по переменному
верхнему пределу?
6. Что такое несобственный интеграл?
7. Как определяются пространства R2 , R3 , Rn ?
11
8. Найдите интеграл
∫∫x
2 2
y dxdy .
00
9. Найдите общее решение для уравнения
1) у″-6у′+9у=е3х;
2) у″+3′+2у=cos2x+2sin2x.
10. Найдите интегралы вида:
n
1) ∫ x dx;
2) ∫ x cos(nx )dx ;
2
4) ∫ cos (nx)dx; 5) ∫ sin 2 (nx)dx;
те же интегралы найдите в пределах
9) от -π до π,
8) от 0 до π,
3) ∫ x sin(nx )dx ;
6) ∫ cos( nx ) sin( mx )dx ;
10) от 0 до l.
l
11) найдите ∫ cos( πn / l )dx , n=1,2,…
0
11. Что такое скалярное произведение? Какими свойствами оно обладает?
Какие векторы называются ортогональными?
12. Какой может быть область сходимости степенного ряда?
13. Найдите область сходимости ряда
∞
∑
xn
2
n =1 n + 1
.
14. Сформулируйте задачу Коши для дифференциального уравнения первого
порядка.
15. Сформулируйте теорему единственности решения задачи Коши для дифференциального уравнения первого порядка.
3
1. Интегральные уравнения. Основные понятия и
определения. Классификация уравнений
Определение: Интегральными уравнениями называются уравнения, содержащие неизвестную функцию под знаком интеграла.
Многие задачи математической физики сводятся к линейным интегральным
уравнениям.
Определение: Интегральное уравнение называется линейным, если в него
неизвестная функция входит линейно.
Пример линейного уравнения
1
x
e
ϕ ( x) − ∫ e ϕ ( y )dy = x −
, 0 ≤ x ≤ 1.
x
xy
0
Пример нелинейного уравнения
1
ϕ ( x) − ∫
xyϕ ( y )
2
0 1 + ϕ ( y)
dy = f ( x), 0 ≤ x ≤ 1.
Решить интегральное уравнение – значит найти такую функцию, которая
обращает данное уравнение в верное тождество.
Интегральные уравнения подразделяются на уравнения первого рода и уравнения второго рода. В уравнения первого рода неизвестная функция входит
только под знаком интеграла. В уравнения второго неизвестная функция входит
как под знаком интеграла, так и вне интеграла.
Уравнения первого и второго рода с постоянными пределами интегрирования называются уравнениями Фредгольма, а уравнения с переменным верхним
пределом называются уравнениями Вольтерра.
Таким образом, линейное интегральное уравнение Фредгольма первого рода
имеет вид:
b
∫ K ( x, y)ϕ ( y)dy = f ( x),
a ≤ x ≤ b.
a
Линейное интегральное уравнение Фредгольма второго рода имеет вид
b
ϕ ( x ) − ∫ K ( x, y )ϕ ( y )dy = f ( x ), a ≤ x ≤ b.
a
Линейное интегральное уравнение Вольтерра первого рода имеет вид
x
∫ K ( x, y)ϕ ( y)dy = f ( x),
a ≤ x.
a
Линейное интегральное уравнение Вольтерра второго рода имеет вид
4
x
ϕ ( x ) − ∫ K ( x, y )ϕ ( y )dy = f ( x ), a ≤ x.
a
Здесь
ϕ ( x) – неизвестная искомая функция;
f(x) – заданная непрерывная функция, называемая свободным членом интегрального уравнения;
K(x,y) – заданная непрерывная функция, называемая ядром интегрального
уравнения.
Уравнения второго рода иногда записывают с параметром λ так:
b
ϕ ( x) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x), a ≤ x ≤ b.
(1)
a
или
x
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x ), a ≤ x.
(2)
a
Тогда уравнения (1) или (2) представляют собой не одно уравнение, а семейство уравнений, зависящее от числового параметра λ .
Если f(x) = 0, то интегральное уравнение называется однородным, в противном
случае оно называется неоднородным.
Уравнения Вольтерра более просты, чем уравнения Фредгольма; уравнения
второго рода более просты, чем уравнения первого рода. Уравнение Вольтерра
можно при некоторых ограничениях рассматривать как частный случай уравнения Фредгольма.
2. Метод последовательных приближений решения
интегральных уравнений
2.1. Построение решения уравнения Фредгольма второго рода при малых
значениях параметра методом последовательных приближений.
Будем рассматривать интегральное уравнение Фредгольма второго рода
b
ϕ ( x) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x), a ≤ x ≤ b,
(3)
a
где ядро K(x,y) и правая часть f(x) – заданные непрерывные функции переменных x, y, ϕ (x ) - искомая функция, λ - параметр.
В случае, когда параметр λ мал, т.е. удовлетворяет условию
1
λ< ,
M
где М – положительное число, такое, что
5
1
⎞2
b
⎛b b
2
⎜
⎟
K ( x, y ) dxdy ≤ M или ∫ K ( x, y ) dy ≤ M , a ≤ x ≤ b ,
⎜∫ ∫
⎟
a
⎝a a
⎠
решение ϕ (x ) уравнения (3) существует и его можно построить методом последовательных приближений.
Суть метода: функция ϕ (x) ищется в виде предела последовательности
ϕ 0 ( x ) = f ( x ),
b
ϕ n ( x ) = f ( x ) + λ ∫ K ( x, y )ϕ n −1 ( y )dy , n = 1,2,...
(4)
a
т.е.
ϕ (x ) = lim ϕ n ( x )
n→∞
и других решений уравнение (3) не имеет.
Методом последовательных приближений найти решение уравнения Фредгольма 2-го рода, предварительно убедившись, что выполнено условие
λ<
1
⎞ 2
⎛bb
1 ⎜
2
= ∫ ∫ K ( x, y ) dxdy ⎟
⎟
M ⎜
⎝a a
⎠
−
:
(5)
b
1. ϕ ( x ) − ∫ xyϕ ( y )dy = 2 x.
a
Р е ш е н и е. Здесь f(x) = 2х, K(x,y) = xy, λ = 1. Имеем
bb
⎛bb
⎞
⎜ K ( x , y ) 2 dxdy ⎟ = M 2 = x 2 y 2 dxdy = 1 / 9.
∫∫
⎜∫∫
⎟
aa
⎝aa
⎠
1
Значит M = 1/3, λ = 1 <
= 3 . Следовательно, условие (5) выполнено.
1/ 3
Построим последовательность приближенных решений (4)
ϕ 0 ( x ) = 2 x,
b
ϕ n ( x ) = 2 x + 1 ⋅ ∫ xyϕ n −1 ( y )dy , n = 1,2,...
a
т.е.
b
1
3
ϕ1 ( x ) = 2 x + 1 ⋅ ∫ xy 2 ydy = 2 x + 2 x ;
a
6
b
2
3
1
3
1
9
ϕ 2 ( x ) = 2 x + 1 ⋅ ∫ xy (2 y + y )dy = 2 x + 2 x ( + );
a
b
2
3
2
9
1
3
ϕ 3 ( x) = 2 x + 1 ⋅ ∫ xy (2 y + y ( + ))dy = 2 x + 2 x( +
a
1 1
+ );
9 27
…
1
3
1
9
ϕ n ( x) = 2 x + 2 x( + +
….
ϕ (x ) = lim ϕ n ( x ) .
1
1
+ ... + n −1 );
27
3
n→∞
Т.к. сумма убывающей геометрической прогрессии
b
1 1 1
1/ 3
1
+ +
+ ... = 1 =
= , то
3 9 27
1 − q 1 − 1/ 3 2
1
ϕ (x ) = 2x + 2x · = 3x.
2
2. ϕ ( x ) +
π
1
cos
π∫
2
yϕ ( y )dy = 1.
0
Ответ: ϕ ( x ) =
1
3. y ( x) −
2π
π
2
.
3
∫ sin x y(t )tdt = 2 sin x.
0
Ответ: y ( x) = 4 sin x .
1
4. y ( x) − π ∫ (1 − x) sin 2πt y (t )dt =
0
1
(1 − x).
2
Ответ: y (x ) = 1 − x.
1
1
5. ϕ ( x ) − ∫ ϕ ( y )dy = sin πx.
2
0
Ответ: ϕ ( x) = sin πx +
2.2. Построение решения уравнения Вольтерра второго рода методом последовательных приближений.
Будем рассматривать интегральное уравнение Вольтерра второго рода
7
2
π
.
x
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x ), a ≤ x.
(6)
a
Методом последовательных приближений функция ϕ ( x ) ищется в виде
предела последовательности
ϕ 0 ( x) = f ( x),
x
ϕ n ( x ) = f ( x ) + λ ∫ K ( x, y )ϕ n −1 ( y )dy , n = 1,2,... ,
(7)
т.е.
ϕ (x ) = lim ϕ n ( x ) .
(8)
a
n→∞
Интегральное уравнение Вольтерра (6) при требовании непрерывности его
ядра K ( x, y ) и правой части f (x ) имеет единственное решение (8) для каждого
конечного значения параметра λ .
Этим существенно отличается интегральное уравнение Вольтерра второго
рода от интегрального уравнения Фредгольма 2-го рода, которое, как будет показано ниже, не для каждого λ может иметь решение, а при некоторых λ , оно
может иметь даже несколько решений.
Обычно полагают ϕ ( x) = f ( x), однако это вовсе не обязательно: удачный
0
выбор ненулевого приближения часто позволяет ускорить сходимость последовательности ϕ (x) к точному решению.
n
Методом последовательных приближений найти решение уравнения Вольтерра
2-го рода:
x
6. ϕ ( x) = 1 − ∫ ( x − y )ϕ ( y )dy.
0
Р е ш е н и е. Положим
ϕ 0 ( x) = 1.
Тогда
x
2
2
x
x
ϕ1 ( x) = 1 − ∫ ( x − y )dy = 1 −
=1−
,
2
2!
0
x
2
2
4
2
4
y
x
x
x
x
ϕ 2 ( x) = 1 − ∫ ( x − y )(1 −
)dy = 1 −
+
=1−
+
,
2!
2
24
2!
4!
0
x
2
4
2
4
6
y
y
x
x
x
ϕ 3 ( x) = 1 − ∫ ( x − y )(1 −
+
)dy = 1 −
+
−
.
2!
4!
2!
4!
6!
0
… .
8
Для n-ого приближения
2n
2k
n
x2 x4 x6
n x
k x
+
−
+ ... + ... + (−1)
= ∑ (−1)
.
ϕ n (x) = 1 −
2!
4!
6!
(2n)! k =0
(2k )!
Тогда решение уравнения
ϕ (x ) = lim ϕ n (x) =
n→∞
∞
∑ (−1)
k x
k =0
2k
( 2k )!
= cos x.
П р о в е р к а: Подставим ϕ ( x ) = cos x в исходное уравнение
x
ϕ ( x) = 1 − ∫ ( x − y )ϕ ( y )dy.
0
Получим
x
x
cos x = 1 − ∫ x cos xdy + ∫ y cos xdy.
0
(9)
0
Находим первый интеграл
x
∫ x cos xdy = x sin y
0
x
0
= x sin x .
(10)
Находим второй интеграл
u=x
x
x
dv = cos ydy
=
y
sin
y
= x sin y + cos x − 1.
y
cos
xdy
=
∫
0
v = sin y
0
du = dy
Подставим (10), (11) в (9)
cos x = 1 − x sin x + x sin x + cos− 1,
cos x = cos x - верно.
(11)
Ответ: ϕ ( x) = cos x.
x
7. ϕ ( x) = 1 + ∫ ϕ ( y )dy.
0
Ответ: ϕ ( x ) = exp x.
2
x
8. y ( x ) = 1 − x + ∫ xy (t )dt.
0
a) y 0 ( x ) = 1 − x 2 ;
2
б) y 0 ( x ) = 1.
x
9. y ( x ) = x + x − ∫ y (t )dt.
0
9
Ответ: y ( x) = 1.
б) y 0 ( x ) =
a) y ( x) = 1;
0
x
10. ϕ ( x ) = 1 + ∫ yϕ ( y )dy ,
x2
+ x.
2
Ответ: y ( x) = x.
ϕ 0 ( x) = 1 .
0
Ответ: ϕ ( x ) = exp( x 2 / 2).
x
11. ϕ ( x ) = x − ∫ ( x − y )ϕ ( y )dy , ϕ 0 ( x ) = 0 .
0
Ответ: ϕ ( x ) = sin x .
x
12. ϕ ( x ) = 1 − ∫ ( x − y )ϕ ( y )dy , ϕ 0 ( x ) = 0 .
0
Ответ: ϕ ( x ) = ch x .
x
13. ϕ ( x ) = 2 + ∫ 2 x − y ϕ ( y )dy , ϕ 0 ( x ) = 0 .
x
0
Ответ: ϕ ( x ) = ( 2e) x .
3. Интегральные уравнения с вырожденными ядрами
3.1 Решение уравнений Фредгольма 2-го рода с вырожденным ядром.
Ядро K(x,y) называется вырожденным, если оно имеет вид
K ( x, y ) =
n
∑ pi ( x )qi ( y ) .
(12)
i =1
Соответствующее интегральное уравнение (1)
b⎛ n
⎞
(13)
ϕ ( x ) − λ ∫ ⎜⎜ ∑ pi ( x ) qi ( y ) ⎟⎟ϕ ( y )dy = f ( x )
⎠
a ⎝ i =1
решается путем сведения к системе линейных алгебраических уравнений следующим образом.
Перепишем уравнение (13) в виде
n
ϕ ( x ) − λ ∑ ci pi ( x ) = f ( x ),
(14)
i =1
где неизвестные ci определяются через искомое решение ϕ (x ) равенствами
b
сi = ∫ qi ( y )ϕ ( y )dy , i = 1,2,..., n
(15)
a
10
Умножая тождество ((14) последовательно на qi ( x ), i = 1, 2..., n, и далее интегрируя обе части на отрезке [a, b] , с учетом (15) получим для неизвестных чисел
ci следующую систему линейных алгебраических уравнений:
n ⎛b
b
⎞
⎜
⎟
ci − λ ∑ ∫ qi ( x ) p j ( x )dx c j = ∫ qi ( x ) f ( x )dx, i = 1, 2,..., n.
⎜
⎟
j =1⎝ a
a
⎠
Введем обозначения
b
(16)
α ij = ∫ qi ( x ) p j ( x )dx ,
(17)
γ i = ∫ qi ( x ) f ( x )dx.
(18)
a
b
a
Тогда система (16) запишется в виде
n
ci − λ ∑α ij c j = γ i , i = 1,2,..., n.
(19)
j =1
Если с1, с2, …, сn – какое-нибудь решение системы (19), то в соответствии с (14) функция
n
ϕ ( x ) = f ( x ) + λ ∑ ci p i ( x )
(20)
i =1
будет решением исходного интегрального уравнения (13). Если же система (19)
несовместна, то и интегральное уравнение не имеет решения.
Этот метод применим, конечно, и в том частном случае, когда уравнение (13) однородное, т.е. f ( x ) = 0 .
Найти все решения или установить неразрешимость заданных уравнений Фредгольма 2-го рода с вырожденным ядром:
π
⎛1
⎞
∫ ⎜⎝ π sin x sin y + y ⎟⎠ϕ ( y )dy = sin 2 x.
−π
Р е ш е н и е:
1
Ядро K ( x, y ) = sin x sin y + y вырожденное, λ = 1. Полагая
14. ϕ ( x ) −
π
p1 ( x ) =
1
π
p 2 ( x ) = 1,
sin x,
q1 ( y ) = sin y ,
q2 ( y ) = y ,
по формулам (17), (18) вычисляем
α11 =
π
1
∫ π sin
−π
2
x dx = 1 ,
α12 =
π
∫ sin x dx = 0 ,
−π
11
α 21 =
γ1 =
π
1
∫ π x sin x dx = 2 ,
−π
α 22 =
∫ sin x ⋅ sin 2 x dx = 0 ,
γ2 =
π
−π
π
∫ x dx = 0 ,
−π
π
∫ x ⋅ sin 2 x dx = −π .
−π
Система (19) принимает вид
c1 ⋅ 0
+ c2 ⋅ 0 = 0 ,
c1 ⋅ ( −2) + c2 ⋅ 1 = −π .
Ее общее решение: c1 = C , c2 = −π + 2C , где С – произвольная постоянная.
Следовательно, любая функция вида
C
1
ϕ ( x ) = sin 2 x + sin x − π + 2C = sin 2 x + C ( sin x + 2) − π
π
π
есть решение заданного уравнения и других решений это уравнение не имеет.
15. Решить уравнение
1
y ( x ) − 2 ∫ xt y (t )dt = x.
0
Р е ш е н и е:
Ядро K ( x, t ) = x t вырожденное, λ = 2. Полагая
p1 ( x ) = p( x ) = x , q1 (t ) = q( x ) = t ,
решение ищем в виде
y ( x ) = f ( x ) + c p или y ( x ) = x + c x .
Система редуцируется к уравнению
c − 2α c = γ ,
где
b
1
1
α = ∫ q( x ) p( x )dx = ∫ x x dx = ,
2
a
0
b
1
a
0
2
5
γ = ∫ q( x ) f ( x )dx = ∫ x xdx = .
Тогда получаем уравнение
2
c−c = .
5
Последнее уравнение не имеет решения относительно с, следовательно, исходное интегральное уравнение также не имеет решения.
16. ϕ ( x ) −
1 2π
∫ cos x sin yϕ ( y )dy = sin x .
π 0
12
Ответ: ϕ ( x ) = cos x + sin x.
17. ϕ ( x ) −
24 1
3
(1 − x 2 )(1 − y )ϕ ( y )dy = x.
∫
7 0
2
Ответ: ϕ ( x ) = x + C (1 − x 2 ).
1
18. ϕ ( x ) − ∫ (1 + x )сos 2πyϕ ( y )dy = x.
0
Ответ: ϕ ( x ) = x.
1
19. ϕ ( x ) − ∫ (2 x − y )ϕ ( y )dy = cos 2πx.
0
Ответ: ϕ ( x ) = cos 2πx.
1
1
6
20. y ( x ) − ∫ (1 + 2 xt ) y (t )dt = − ( x + 3).
0
1
Ответ: y ( x ) = x + .
2
1
3
2
−1
21. y ( x ) − ∫ ( xt + x 2 (t − 1)) y (t )dt = 0.
Ответ: y ( x ) =
22. y ( x ) −
4π
cos
π 0∫
2
5
x + x2.
2
( x − t ) y (t )dt = sin 2 x.
Ответ: Решения нет.
1
5
3
1
6
23. y ( x ) − ∫ ( x − t ) y (t )dt = x + x − .
0
Ответ: y ( x ) = 1 + x .
24. y ( x ) −
1 π
π −∫π
cos( x − t ) y (t )dt = 0.
Ответ: y ( x ) = C1 cos x + C 2 sin x.
1
25. y ( x ) − 3∫ ( x 2 t 2 − 4 xt + 1) y (t )dt = 2π 2 cos 2πx.
0
5
Ответ: y ( x ) = 2π 2 cos 2πx + ( 2 x 2 − 1).
3
13
1
26. y ( x ) − ∫ ( xt + x 2 ) y (t )dt = 0.
−1
Ответ: y ( x ) = 0.
3.2. Решение уравнений Вольтера 2-го рода с вырожденным ядром.
Пусть исходное интегральное уравнение (2) при λ = 1имеет вид:
x⎛ n
⎞
⎟
⎜
⎠
a ⎝ i =1
(уравнение с вырожденным ядром)
ϕ ( x ) − ∫ ⎜ ∑ pi ( x ) qi ( y ) ⎟ϕ ( y )dy = f ( x )
(21)
Запишем его следующим образом:
x⎛ n
⎞
ϕ ( x ) = f ( x ) − ∫ ⎜ ∑ pi ( x ) qi ( y ) ⎟ϕ ( y )dy .
⎟
⎜
⎠
a ⎝ i =1
Вводя функции
(22)
x
u1 = ∫ q1 (t )ϕ (t )dt ,
0
. . . . . . .
(23)
x
un = ∫ qn (t )ϕ (t )dt ,
0
и подставляя их в (22), заключаем, что решение интегрального уравнения (21)
имеет вид
n
ϕ ( x ) = f ( x ) + ∑ pi ( x )ui ( x ).
(24)
i =1
Далее, дифференцируя соотношения (23) и подставляя вместо ϕ (x ) выражение
(24), получаем для неизвестных функций ui ( x ) систему дифференциальных
уравнений
n
u' i = q1 ( x ) f ( x ) + ∑ q1 ( x ) pi ( x )ui ( x ),
i =1
. . . . . . . . . . . .
n
u' n = qn ( x ) f ( x ) + ∑ qn ( x ) pi ( x )ui ( x ).
i =1
Из (23) при x = 0 находим начальные условия: ui (0) = ... = un (0). Определив
функции ui (x ) и подставив их в (24), получим решение ϕ (x ) интегрального
уравнения (21).
14
27. Решить интегральное уравнение
x ch y
ϕ ( x) = 1 + ∫
ϕ ( y )dy.
ch
x
0
Р е ш е н и е:
x
Полагая u( x ) = ∫ ch (t ) ϕ ( y )dy , получим
0
1
u( x ).
ch ( x )
Далее, дифференциальное уравнение для u ( x ) имеет вид
1
⎛
⎞
u ' ( x ) = ch ( x ) y ( x ) = ch ( x )1 + ⎜1 +
u( x ) ⎟
⎝ ch ( x )
⎠
или
ϕ ( x) = 1 +
u ' − u = ch x .
Решая это уравнение с учетом начального условия u (0) = 0 , находим
1
u ( x ) = ( xe x + sh x ) , откуда
2
1 xe x + sh x
ϕ ( x) = 1 +
.
2 ch ( x )
x
x
28. ϕ ( x ) = e + ∫ ϕ ( y )dy.
0
Ответ: ϕ ( x ) = e x ( x + 1).
x
29. ϕ ( x ) = x − 1 + ∫ ( x − y )ϕ ( y )dy.
0
30. y ( x ) =
1
1+ x2
Ответ: ϕ ( x ) = − e − x .
x
+ ∫ sin ( x − t ) y (t ) dt.
0
Ответ: y ( x ) =
1
1
+ x arct x − ln(1 + x 2 ).
2
1+ x2
x
31. y ( x ) = 1 + ∫ t y (t ) dt.
0
2
Ответ: y ( x ) = e x / 2 .
15
4. Понятие итерированного ядра и резольвенты.
Решение интегрального уравнения Фредгольма второго рода с
помощью резольвенты
Часто вместо одного уравнения рассматривают семейство уравнений
b
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x ), a ≤ x ≤ b,
(25)
a
соответствующих различным значениям числового параметра λ . Предполагается, что λ фиксировано. Будем решать уравнение (25) методом последова1
. Взяв в
тельных приближений при условии выполнения неравенства (5) λ <
M
качестве нулевого приближения ϕ 0 ( x ) = f ( x ), получим
b
b
a
a
ϕ1 ( x ) = f ( x ) + λ ∫ K ( x, y ) f ( y )dy = f ( x ) + λ ∫ K1 ( x, y ) f ( y )dy ,
где K1 ( x, y ) = K ( x, y );
(26)
b
b
b
a
a
a
ϕ 2 ( x ) = f ( x ) + λ ∫ K ( x, y )ϕ1 ( y )dy = f ( x ) + λ ∫ K ( x, y )( f ( y ) + λ ∫ K1 ( x, t ) f (t )dt )dy =
b
b
b
a
a
a
= f ( x ) + λ ∫ K1 ( x, y ) f ( y )dy + λ2 ∫ K ( x, y )( ∫ K1 ( x, t ) f (t )dt )dy =
/ поменяем порядок интегрирования /
b
b b
a
b
a a
b
a
b
a
b
= f ( x ) + λ ∫ K1 ( x, y ) f ( y )dy + λ2 ∫ ( ∫ K ( x, y ) K1 ( x, t )dy ) f (t )dt =
= f ( x ) + λ ∫ K1 ( x, y ) f ( y )dy + λ2 ∫ K 2 ( x, t ) f (t )dt =
= f ( x ) + λ ∫ K1 ( x, y ) f ( y )dy + λ2 ∫ K 2 ( x, y ) f ( y )dy ,
a
(27)
a
b
где K 2 ( x, y ) = ∫ K ( x, t ) K1 (t , y )dt;
a
n
b
j =1
a
Вообще ϕ n ( x ) = f ( x ) + ∑ λ j ∫ K j ( x, y ) f ( y )dy =
b⎛ n
⎞
= f ( x ) + λ ∫ ⎜ ∑ λ j −1K j ( x, y ) ⎟ f ( y )dy ,
⎜
⎟
a ⎝ j =1
⎠
(28)
16
b
где K j ( x, y ) = ∫ K ( x, t ) K j −1 (t , y )dt; j = 2,3,...
(29)
a
Ядра K j ( x, y ) называются итерированными (повторными) ядрами. Пользуясь понятием итерированных ядер, последовательным приближениям (26)-(28)
можно придать вид
b⎛ n
⎞
ϕ n ( x ) = f ( x ) + λ ∫ ⎜ ∑ λ j −1K j ( x, y ) ⎟ f ( y )dy.
(30)
⎜ j =1
⎟
a⎝
⎠
1
При λ <
(31)
M
ряд
∞
R ( x , y , λ ) = ∑ λ j −1 K j ( x , y )
(32)
j =1
сходится равномерно при a ≤ x ≤ b, a ≤ y ≤ b к функции R( x, t , λ ) , называемой
резольвентой ядра K ( x, t ). Следовательно, (30) в пределе при n → ∞ переходит
в формулу
b
ϕ ( x ) = f ( x ) + λ ∫ R ( x, y , λ ) f ( y )dy ,
(33)
a
выражающую решение интегрального уравнения через резольвенту.
33. С помощью итерированных ядер найти резольвенту и решение интегрального уравнения
1 1 x
ϕ ( x) −
ϕ ( y )dy = 1 + x 2 .
∫
2
ln 2 0 1 + y
Р е ш е н и е. В данном случае K ( x, y ) =
x
1 + y2
и для итерированных ядер на
основании (29) получаем
x
,
K1 ( x, y ) = K ( x, y ) =
1+ y2
1
1
0
1
2
2
01 + s 1 + y
1
K 2 ( x, y ) = ∫ K ( x, t ) K1 (t , y )dt = ∫
x
s
dy =
ln 2 x
,
2 1+ y2
2
ln 2
x
s
x
⎛ ln 2 ⎞
,
K 3 ( x, y ) = ∫ K ( x, t ) K 2 (t , y )dt =
dy = ⎜
⎟
∫
2 1 + s2 1+ y2
2 ⎠ 1+ y2
⎝
0
0
............
⎛ ln 2 ⎞
K j ( x, y ) = ⎜
⎟
⎝ 2 ⎠
j −1
x
1+ y
2
.
17
Поэтому резольвента ядра равна
∞
R( x, t , λ ) = ∑ λ
j −1
j =1
j −1
∞ ⎛ ln 2
K j ( x, t ) = ∑ ⎜
λ ⎞⎟
⎠
j =1 ⎝ 2
x
1+ y2
причем этот ряд сходится в области
2
λ<
.
ln 2
Заметим, что в рассматриваемом случае
M
2
11
2
x2
11
= ∫ ∫ K ( x, y ) dxdy = ∫ ∫
=
x
1
,
ln 2 1 + y 2
λ
1−
2
(34)
dxdy =
π +2
,
2 2
24
(
1
y
)
+
00
00
т.е. условие (31) приводит к неравенству
6
.
(35)
λ <2
π +2
6
2
<
, то из сравнения (34) и (35) видно, что в рассматриваемом
Так как
π + 2 ln 2
случае область сходимости ряда Неймана для резольвенты шире, чем это гарантируется условием
λ<
⎛bb
1
⎞ 2
1 ⎜
2
= ∫ ∫ K ( x, y ) dxdy ⎟
⎟
M ⎜⎝ a a
⎠
−
.
Далее, для заданного уравнения λ =
R( x, y ,
1
и, следовательно,
ln 2
1
x
)=2
. Решение уравнения на основании (33) равно
2
ln 2
1+ y
ϕ ( x) = 1 + x 2 +
1 1
x
2
2
(1 + y 2 )dy = 1 +
x + x2.
∫
2
ln 2 0 1 + y
ln 2
Методом итерированных ядер найти резольвенту и решение заданных интегральных уравнений:
1 π
34. ϕ ( x ) −
∫ ϕ ( y )dy = sin x.
2π 0
1
,
1 − πλ
ϕ ( x ) = sin x.
Ответ: R ( x, y , λ ) =
18
ln 2 1 x + y
35. ϕ ( x ) −
ϕ ( y )dy = x.
∫2
2 0
1
,
3λ
1−
2 ln 2
2 ln 2 − 1 x +1
ϕ ( x) =
2
+ x.
ln 2
Ответ: R( x, y , λ ) = 2 x + y
1
36. y ( x ) − π ∫ x sin 2πt y (t )dt = cos 2πx.
0
Ответ: R ( x, t , λ ) = x sin 2πt ,
y ( x ) = cos 2πx.
11
37. y ( x ) − ∫ xe t y (t )dt = e − x .
20
xe t
Ответ: R ( x, t , λ ) =
,
1− λ
y ( x ) = e − x + x.
π /2
38. y ( x ) − ∫ sin x cos t y (t )dt = 1.
0
2
sin x cos t ,
2−λ
y ( x ) = 1 + 2 sin x.
Ответ: R ( x, t , λ ) =
5. Характеристические числа и собственные функции. Теоремы
Фредгольма
5.1. Собственное число (значение) и собственный вектор матрицы.
⎛x ⎞
⎛a a ⎞
X = ⎜⎜ 1 ⎟⎟ - собственный вектор матрицы A = ⎜⎜ 11 12 ⎟⎟ ,
⎝ x2 ⎠
⎝ a 21 a 22 ⎠
если AX = λX ,
т.е. если после преобразования вектора X с помощью матрицы A получаем
вектор Y = AX , параллельный вектору A . λ - называется собственным числом
матрицы A .
AX − λX = 0 или ( A − λE ) X = 0.
Получаем однородную систему уравнений
a12
⎞ ⎛ x1 ⎞
⎛ ( a11 − λ )
⎟⎜ ⎟=0
⎜⎜
( a 22 − λ ) ⎟⎠ ⎜⎝ x 2 ⎟⎠
a 21
⎝
19
или
( a11 − λ ) x1 +
a12 x 2 = 0,
(36)
a 21 x1 + ( a 22 − λ ) x 2 = 0.
Для того чтобы однородная система (36) имела ненулевое решение необходимо
и достаточно, чтобы определитель системы был равен нулю.
( a11 − λ )
a12
=0
( a 22 − λ )
a 21
или
det( A − λE ) = 0.
(37)
Это характеристическое уравнение матрицы A . Из этого характеристического
уравнения находятся собственные значения λ матрицы A . Каждому собственному значению λ соответствует собственный вектор, координаты которого определяются из системы (36) при соответствующем λ .
5.2. Характеристические числа и собственные функции интегрального
уравнения.
Значение параметра λ , при которых однородное уравнение Фредгольма
b
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = 0
(38)
a
имеет ненулевые (нетривиальные) решения ϕ ( x ) ≠ 0 , называются характеристическими числами этого уравнения или ядра K ( x, y ) , а каждое ненулевое
решение – собственной функцией, соответствующей характеристическому числу λ . Заметим, что число λ = 0 не является характеристическим, т.к. при λ = 0
уравнение (38) имеет лишь нулевое решение. Если λ - характеристическое
число, то число μ = 1 / λ называется собственным числом интегрального уравнения. При этом μ ≠ 0.
Из результатов п. 3. следует, что в случае уравнения с вырожденным ядром
b⎛ n
⎞
ϕ ( x ) − λ ∫ ⎜⎜ ∑ pi ( x ) qi ( y ) ⎟⎟ϕ ( y )dy = 0
(39)
⎠
a ⎝ i =1
всякое решение имеет вид
n
ϕ ( x ) = λ ∑ ci pi ( x ) ,
(40)
i =1
где C = ( c1 , c2 , ..., cn )T - решение однородной системы
n
ci − λ ∑ α ij c j = 0, i = 1,2,..., n
j =1
или в матричной форме
20
( E − λA) C = 0
с матрицей A = α ij . По формуле (17)
(41)
{ }
b
α ij = ∫ qi ( x ) p j ( x )dx , i, j = 1,2,..., n.
(42)
a
Заметим, что если заменить λ на 1 / μ , то система ( 41) принимает вид
( A − μE ) C = 0, μ ≠ 0.
(43)
Отсюда следует, что собственные числа интегрального уравнения (39) совпадают с собственными числами матрицы А, а собственные функции определяются соотношением (40), где C = ( c1 , c2 , ..., cn )T - соответствующие собственные
векторы этой матрицы.
39. Найти характеристические числа и собственные функции уравнения
1
ϕ ( x ) − λ ∫ ( xy − 2 x 2 )ϕ ( y )dy = 0.
0
Р е ш е н и е. Ядро K ( x, y ) = x y − 2 x 2 - вырожденное. Полагая
p1 ( x ) = x,
p 2 ( x ) = −2 x 2 ,
q1 ( y ) = y ,
q2 ( y ) = 1,
Найдем элементы матрицы А в (43)
1
1
1
α11 = ∫ q1 ( x ) p1 ( x )dx = ∫ x 2 dx = ,
3
0
0
1
1
0
0
1
2
1
1
0
1
0
1
0
0
1
2
α12 = ∫ q1 ( x ) p2 ( x )dx = −2 ∫ x 3dx = − ,
2
3
α 22 = ∫ q2 ( x ) p2 ( x )dx = −2 ∫ x 2 dx = − .
α 21 = ∫ q2 ( x ) p1 ( x )dx = ∫ xdx = ,
Характеристическое уравнение для определения собственных чисел матрицы А
имеет вид
1
1
−μ −
1
1
1
3
2
= μ2 + μ +
= (1 + ) 2 = 0,
det( A − μE ) =
1
2
3
36
6
− −μ
2
3
откуда μ = −1 / 6 - единственное собственное число матрицы А. Соответствующие собственные векторы находим из системы уравнений
1⎞
⎛1
⎜ − ⎟⎛ c ⎞ 0
1
2 ⎟⎜ 1 ⎟ = ⎛⎜ ⎞⎟,
( A + E) C = ⎜ 2
⎜1
1 ⎟⎜⎝ c ⎟⎠ ⎜⎝ 0 ⎟⎠
6
−
⎜
⎟ 2
⎝2
2⎠
общее решение которой c = C , c = C , где С – произвольная постоянная. Сле1
2
довательно, окончательно получаем, что заданное интегральное уравнение име-
21
ет единственное характеристическое число λ =
1
μ
= −6,
а соответствующие
собственные функции имеют вид
ϕ ( x) = −6(c x − 2c x 2 ) = C ( x − 2 x 2 ),
1
2
где С – произвольная постоянная.
Интегральное уравнение может вообще не иметь характеристических чисел
(например, в том случае, когда ядро K ( x, t ) - вырожденное, матрица А в (41)
нулевая) либо не иметь действительных характеристических чисел.
40. Найти характеристические числа и собственные функции уравнения
π
ϕ ( x ) − λ ∫ x cos y ϕ ( y )dy = 0.
−π
Р е ш е н и е. Имеем
p1 ( x ) = x,
q1 ( y ) = cos y ,
ϕ ( x) − λxc = 0,
c=
π
∫ cos yϕ ( y )dy,
−π
откуда
π
c − λc ∫ x cos x dx = 0.
π
Но
−π
∫ x cos x dx = 0,
поэтому при любом λ последнее уравнение имеет только
−π
одно решение: c = 0 . Следовательно, при любом λ интегральное уравнение
имеет только тривиальное решение, т.е. не имеет характеристических чисел.
Найти характеристические числа и собственные функции заданных интегральных уравнений с вырожденными ядром (ограничиться случаем действительных характеристических чисел):
1
41. ϕ ( x ) − λ ∫ (1 + 2 x ) y ϕ ( y )dy = 0.
0
6
Ответ: λ = , ϕ ( x) = C (1 + 2 x).
7
1
42. y ( x ) − λ ∫ ( x + t ) y (t )dt = 0.
0
Ответ: μ =
1, 2
22
1
1
±
, y ( x) = C ( 3x ± 1).
2
3 1, 2
1
43. y ( x ) − λ ∫ (1 + 2 x )t y (t )dt = 0.
0
6
Ответ: λ = , y ( x) = C (1 + 2 x).
7
44. ϕ ( x ) − λ
1
∫ x ϕ ( y )dy = 0.
−1
Ответ: λ = 1, ϕ ( x) = C x .
5.3. Теоремы Фредгольма.
Для уравнения Фредгольма 2-го рода вида
b
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = f ( x ),
(44)
a
где a и b - конечные числа, а ядро K ( x, y ) и свободный член f (x) интегрируемы с квадратом в области a ≤ x, y ≤ b и на отрезке [a, b] (в частности, непрерывны), справедливы следующие теоремы Фредгольма (при формулировке которых мы ограничимся случаем действительного ядра K ( x, y ) .
1. Однородное уравнение
b
ϕ ( x ) − λ ∫ K ( x, y )ϕ ( y )dy = 0
(45)
a
имеет либо конечное, либо счетное множество характеристических чисел; если этих чисел счетное множество, то они стремятся к бесконечности.
2. Если λ - характеристиченсоке число, то уравнение (45) и сопряженное ему однородное уравнение
b
ϕ ( x ) − λ ∫ K * ( x, y )ϕ ( y )dy = 0,
(46)
a
*
где K ( x, y ) = K ( y , x ) , имеют одно и то же, и при том конечное, число независимых решений.
3. Альтернатива Фредгольма: либо неоднородное уравнение (44) имеет
одно и только одно решение для любой функции f ( x) ∈ L [a, b] , либо соответ2
ствующее однородное уравнение (45) имеет по крайней мере одно нетривиальное решение. (Другими словами, если число λ не является характеристическим,
то уравнение (44) имеет, и при том единственное, решение для любой функции
f ( x) ∈ L [a, b] ).
2
4. Если λ - характеристическое число, то для того чтобы уравнение (44)
имело решение, необходимо и достаточно, чтобы свободный член f (x) был
23
ортогонален любому решению ϕ * ( x) однородного сопряженного уравнения
(46), т.е.
b
*
∫ f ( x ))ϕ
( y )dy = 0.
a
Проиллюстрируем теоремы Фредгольма на примере интегрального уравнения с
вырожденным ядром.
45. Исследовать решения интегрального уравнения
π
ϕ ( x ) − λ ∫ ( x 2 cos y + x sin y )ϕ ( y )dy = cos x
−π
в зависимости от значений параметра λ .
Р е ш е н и е. Решение интегрального уравнения сводится к решению неоднородной системы
( E − λA)C = F ,
(47)
где F = ( f1 , f 2 ,..., f n ) T ,
2
α11 =
α 21 =
∫ cos x ⋅ x
a
p 2 ( x ) = x,
q2 ( y ) = sin y.
p1 ( x ) = x ,
q1 ( y ) = cos y ,
π
b
f i = ∫ qi ( x ) f ( x )dx. В рассматриваемом случае имеем
2
π
α12 = ∫ cos x ⋅ x dx = 0,
dx = 4π ,
−π
π
∫ sin x ⋅ x
2
π
π
α 22 = ∫ sin x ⋅ x dx = −2π ,
dx = 0,
−π
π
f1 = ∫ cos 2 x dx = π ,
−π
π
π
f 2 = ∫ sin x cos x dx = 0.
−π
Система (47) имеет вид
⎛1 − 4πλ 0
⎞ ⎛ c1 ⎞ ⎛ π ⎞
⎜⎜
⎟ ⎜ ⎟ = ⎜ ⎟.
1 + 2πλ ⎟⎠ ⎜⎝ c2 ⎟⎠ ⎜⎝ 0 ⎟⎠
⎝ 0
(48)
Характеристическое уравнение
det( E − λA) = (1 + 2πλ )(1 − 4πλ ) = 0
1
1
λ
=
−
, являющиеся характеристическими
имеет корни λ1 =
и 2
2π
4π
числами соответствующего однородного уравнения.
1
1
При любом λ ≠
система (48) имеет единственное решение
,−
4π 2π
24
π
, c2 = 0;
1 − 4πλ
соответствующее решение интегрального уравнения:
λπ
1
1
ϕ ( x ) = cos x +
x2, λ ≠
.
,−
1 − 4πλ
4π 2π
1
При λ = λ1 =
из (48) получаем
4π
⎛0 0 ⎞ c
⎟ ⎛ 1 ⎞ ⎛π ⎞
⎜
3
⎟ ⎜⎜ c ⎟⎟ = ⎜⎜ 0 ⎟⎟.
⎜0
⎟⎝ 2 ⎠ ⎝ ⎠
⎜
⎝
2⎠
Эта система, а вместе с ней и исходное интегральное уравнение решения не
имеют.
1
При λ = λ2 = −
система (37) принимает вид
2π
⎛ 3 0 ⎞ ⎛ c1 ⎞ ⎛ π ⎞
⎜⎜
⎟⎟ ⎜⎜ ⎟⎟ = ⎜⎜ ⎟⎟
0
0
⎝
⎠ ⎝ c2 ⎠ ⎝ 0 ⎠
c1 =
и имеет решения c1 =
π
3
, c2 = C. Соответствующие решения интегрального
уравнения таковы:
1
6
ϕ ( x ) = cos x + λ2 ( c1 x 2 + c2 x ) = cos x − + Cx,
где С – произвольная постоянная.
Исследовать решения заданных уравнений с вырожденным ядром при различных значениях параметра λ .
1
46. ϕ ( x ) − λ ∫ x (1 + y ) ϕ ( y ) dy = x 2 .
0
6
1 7λ
Ответ: при λ ≠ , ϕ ( x ) = x 2 +
x,
5
2 6 − 5λ
6
при λ = решения нет.
5
1
47. ϕ ( x ) − λ ∫ x ϕ ( y ) dy = sin 2πx.
0
1
48. y ( x ) − λ ∫ ( x + t ) y (t ) dt =
−1
Ответ: при λ ≠ 2, ϕ ( x ) = sin 2πx,
при λ = 2 ϕ ( x ) = sin 2πx + Cx.
1 3
+ x.
2 2
25
3
1 3
y( x) = + x +
2
2 2
2
λ
+
((1 + 2λ ) x + 1 + λ ),
4
3
1 − λ2
3
3
решений нет.
при λ = ±
2
Ответ: при λ ≠ ±
π
49. y ( x ) − λ ∫ cos( x + t ) y (t ) dt = 1.
0
Ответ: при λ ≠ ±
при λ =
2
π
2
2λ
1+
λπ
sin x ,
2
y ( x ) = 1 − sin x + C cos x,
π
при λ = −
y( x) = 1 +
2
π
решений нет.
6. Физические примеры.
Многие задачи математической физики приводят к линейным интегральным
уравнениям. Рассмотрим некоторые примеры.
1. Если внешнее воздействие на какую-либо линейную систему описывается
функцией f ( x ) ( a ≤ x ≤ b), то результат этого воздействия описывается функцией
b
f€( x ) = ∫ G ( x;ξ ) f (ξ ) dξ ,
(49)
a
где G ( x;ξ ) - функция влияния, определяемая рассматриваемой системой. Например, f (x ) может означать плотность нагрузки, распределяемой вдоль балки, а f€( x ) - соответствующий прогиб и т.п.
Допустим, что вид воздействия нам неизвестен, но известен отклик системы на это воздействие и требуется по этому отклику восстановить воздействие. Тогда в соотношении (49) функция f€( x ) (как и G ( x;ξ ) ) будет заданной, а
f (x ) - искомой, т.е. мы приходим к интегральному уравнению - линейному
интегральному уравнению Фредгольма первого рода.
2. Бывают случаи, когда известной оказывается некоторая линейная комбинация ϕ ( x ) = α f ( x ) + β f€( x ) функций, описывающих внешнее воздействие и
соответствующий отклик. Тогда для восстановления внешнего воздействия потребуется решить интегральное уравнение
26
b
α f ( x ) + β ∫ G ( x;ξ ) f (ξ ) dξ = ϕ ( x ).
(50)
a
Это линейное интегральное уравнение Фредгольма второго рода с искомой
функцией f (x ) .
3. Приведем еще одну задачу, сводящуюся к уравнению вида (50). Рассмотрим уравнение вынужденных поперечных колебаний струны, закрепленной при
x = 0 и x = l [11]:
'
ρutt' ' = Pu 'xx
+ f ( x, t ) .
(51)
Если внешнее воздействие является гармоническим,
f ( x.t ) = ϕ ( x ) cos ω t ,
и нас интересует вынужденное колебание, происходящее с той же частотой ω ,
т.е. u( x, t ) = v ( x ) cos ω t ,
то после подстановки в (51) мы приходим к краевой задаче
Pv ' ' = − ρω 2 v − ϕ ( x ) (0 ≤ x ≤ l ), ν (0) = 0, v (l ) = 0.
При ω = 0 получилась бы задача на стационарное отклонение струны под действием внешней нагрузки. Эта задача рассмотрена, например, в [11]. Решение
получено в виде
l
⎡ 1
⎤
v ( x ) = ∫ G ( x;ξ ) ⎢ − ϕ (ξ )⎥ dξ ,
⎣ P
⎦
0
(52)
где функция влияния в данной задаче равна
⎧⎪− (l − ξ ) xl −1 (0 ≤ ξ ≤ x ≤ l ),
G ( x,ξ ) = ⎨
⎪⎩− ξ (l − x )l −1 (0 ≤ x ≤ ξ ≤ l ).
Однако при ω ≠ 0 к внешней нагрузке ϕ ( x ) добавляется инерционный член
ρω 2 v ( x ), зависящий от искомого решения. Если на минуту считать его известным и воспользоваться решением (52), то придем к соотношению
ρω 2
l
l
0
0
1
v( x ) = −
G ( x;ξ )v (ξ ) dξ − ∫ G ( x;ξ )v (ξ ) dξ .
∫
P
P
Но так как v ( x ) на самом деле неизвестна, то это соотношение представляет
собой интегральной уравнение, причем того же типа, что и (50).
4. Уравнения Вольтерра второго рода типичны при описании физических
процессов, связанных с явлениями последействия. В этих уравнениях переменная х обычно обозначает время. Тогда состояние системы, характеризуемое
функцией y ( x ) , определяется внешним воздействием f ( x ) и зависит от состояния системы в предшествующие моменты времени. Ядро K ( x, s ) описывает
величину последствия состояния системы в момент s на состояние системы в
момент x > s.
27
В качестве примера рассмотрим электрическую цепь, изображенную на рис.
1. Пусть в катушке индуктивности не проявляется явление гистерезиса. Тогда
поток индукции в катушке Ф связан с током I L соотношением Φ = L ⋅ I L . Согласно известным формулам электродинамики имеем
1
du
dΦ
⋅ IC = ,
= u.
R ⋅ I R = u,
C
dt
dt
IC
v
IL
C
IR
L
R
Рис. 1.
Используя закон Кирхгофа I R + I C + I L = 0, приходим к следующему дифференциальному уравнению относительно u :
u
u&
+ u&&C + = 0.
L
R
Пусть теперь катушка снабжена магнитным сердечником, в котором проявляется гистерезис. Тогда вместо соотношения Φ = L ⋅ I L нужно использовать
более сложное, учитывающее зависимость Ф не только от значения I L в момент t , но и в предшествующие моменты времени (эффект последействия). Это
видоизмененное соотношение таково:
t
Φ (t ) = L ⋅ I L + ∫ M (t − τ ) I L (τ )dτ .
t0
Здесь M (t − τ ) - функция, учитывающая влияние значения I L в момент τ на
величину Ф в момент t и определяемая обычно эмпирическим способом.
Имеем
u
I L = − I R − I C = − − Cu&.
R
Отсюда
t −τ
I (τ )
u (t ) = u (t 0 )e R C L
dτ
−
C
и, следовательно,
t
t −
t0
t0
ξ − t0
Φ (t ) = ∫ u (ξ ) dξ = u(t0 ) ∫ e R C dξ −
ξ − ξ −τ
t
R C I L (τ )
d
ξ
e
dτ =
∫ ∫
C
t0
t0
28
t
~
= f (u(t0 ), t ) + ∫ K (t − τ ) I L (τ ) dτ ,
t0
где
t −
ξ −τ
1
~
K (t − τ ) = − ∫ e R C dξ .
C
τ
Подставляя полученное выражение для Ф в уравнение, связывающее Ф и I L ,
получим
t
t
t0
t0
~
f (u (t0 ), t ) + ∫ K (t − τ ) I L (τ ) dτ = L I L (t ) + ∫ M (t − τ ) I L (τ )dτ
или окончательно для
I L (t ) =
I L (t ) имеем
t
1
∫ K (t − τ ) I L (τ ) dτ + f (u(t0 ), t ) L ,
t0
[
]
1
~
где K (ξ ) = K (ξ ) − M (ξ ) . Таким образом, приходим к интегральному уравнеL
нию Вольтера второго рода.
Ряд других конкретных физических задач, приводящих к интегральным
уравнениям, содержится, например, в [2, 11].
Появление интегральных уравнений при исследовании краевых задач является естественным, т.к. такие уравнения связывают между собой значения известных и неизвестных функций на конечном интервале, а не на бесконечно
малом, как дифференциальные уравнения.
7. Связь интегральных уравнений с дифференциальными
Интегральные уравнения Вольтера 2-го рода используются обычно при описании динамики различных процессов в системах.
В частности, всякая задача Коши для линейного дифференциального уравнения
dny
dx n
+ a1 ( x )
d n −1 y
dx n −1
+ ... + a n ( x ) y = f ( x ),
y ( x0 ) = y 0 , y ' ( x0 ) = y1 , …, y ( n −1) ( x0 ) = y n −1
может быть сведена к решению неоднородного линейного интегрального уравнения Вольтера 2-го рода.
50. Составить интегральное уравнение, соответствующее задаче Коши
u"+2u'+ u = x 2 , u(0) = 1, u' (0) = 0.
Р е ш е н и е: Положим
29
u ' ( x ) = y ( x ).
Интегрируя (53) с учетом начальных условий, последовательно находим
x
x
u' ( x ) = u ' (0) + ∫ y (t )dt = ∫ y (t )dt ,
0
(54)
0
x
s
x
u ( x ) = u (0) + ∫ ds ∫ y (t )dt =1 + ∫ ( x − t ) y (t )dt.
0
(53)
0
(55)
0
Подставляя (53) – (55) в исходное дифференциальное уравнение, получаем
x
⎛ x
⎞
⎜
y ( x ) + 2 ∫ y (t )dt + 1 + ∫ ( x − t ) y (t )dt ⎟ = x 2 ,
⎜
⎟
0
⎝ 0
⎠
или
x
2
y ( x ) = x − 1 − ∫ ( 2 + x − t ) y (t )dt.
(56)
0
Таким образом, показано, что если u (x ) – решение исходной задачи Коши,
то функция y ( x ) = u" ( x ) удовлетворяет интегральному уравнению (56). Обратно, если y (x ) - решение этого уравнения, то функция u (x ) , определяемая соотношением (55), удовлетворяет как исходному дифференциальному уравнению, так и начальным условиям. Следовательно, рассматриваемая задача Коши
эквивалентна интегральному уравнению (56).
Проверить, что данные функции являются решениями соответствующих интегральных уравнений:
51. y ( x ) = e
2x
x
y ( x ) = e + ∫ e x − t y (t )dt.
x
,
0
52. y ( x ) = xe
x2 / 3
x
,
y ( x ) = x + ∫ x t y (t )dt.
0
53. y ( x ) = e
−x x
(
2
2
+ 1),
y( x) = e
−x
x
+ ∫ e − ( x − t ) sin( x − t ) y (t )dt.
0
Составить интегральные уравнения, соответствующие следующим задачам
Коши:
54. u '+2 xu = e x , u(0) = 1.
x
x
Ответ: y ( x ) = e − 2 x − ∫ 2 xy (t )dt.
0
55. u"−2u '+u = 0, u( 2) = 1, u' ( 2) = −2.
30
x
Ответ: y ( x ) = 2 x − 9 + ∫ ( 2 − x + t ) y (t )dt.
2
x
56. u"− sin x ⋅ u'+ e u = x, u(0) = 1, u ' (0) = −1.
x
x
Ответ: y ( x ) = e ( x − 1) − sin x + x + ∫ (sin x − e x ( x − t )) y (t )dt.
0
CПИСОК РЕКОМЕНДУЕМОЙ ЛИТЕРАТУРЫ
1. Голосков Д.П. Уравнения математической физики. Решение задач в системе Maple. Учебник для вузов – СПб.: Питер, 2004. – 539 с.
2. Васильев А.Б., Тихонов Н.А. Интегральные уравнения. М.: Изд-во Моск.
ун-та, 1989. 156 с.
3. Сборник задач по уравнениям математической физики / Владимиров В.С. ,
Михайлов В.П., Вашарин А.А., и др. – М.: Наука, 1982. – 256 с.
4. Сборник задач по математике для втузов. Ч. 4. Методы оптимизации.
Уравнения в частных производных. Интегральные уравнения: Учеб. пособ. /Вуколов Э.А., Ефимов А.В. и др. – М.: Наука, 1980. 304 с.
5. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1981. 542 с.
6. Арсенин В.Я. Математическая физика. Основные уравнения и специальные функции.: М. – Наука, 1965. 367 с.
М.: Наука, 1977.
7. Владимиров В.С., Жаринов В.В. Уравнения математической физики. –
М.: ФИЗМАТЛИТ, 2003. 400 с.
8. Владимиров В.С. Уравнения математической физики. М. – Наука, 1976.
527 с.
9. Забрейко П.П. и др. Интегральные уравнения. – М.: Наука, 1968.
10. Трикоми Ф. Интегральные уравнения. – М.: ИЛ, 1960. 300 с.
11. Мышкис А.Д. Математика для втузов. Специальные курсы. М.: Наука,
1971. 632 с.
31
СОДЕРЖАНИЕ
1.
2.
3.
4.
5.
6.
7.
Вопросы и задачи для самопроверки ………….………………..………
Интегральные уравнения. Основные понятия и определения. Классификация уравнений …………………………………………………..
Метод последовательных приближений решения интегральных
уравнений………………………………………………………………..
2.1. Построение решения уравнения Фредгольма второго рода при
малых значениях параметра методом последовательных приближений…….………………………………………………………………….
2.2 Построение решения уравнения Вольтерра второго рода методом последовательных приближений………………………………….
Интегральные уравнения с вырожденными ядрами…………………
3.1. Решение уравнений Фредгольма 2-го рода с вырожденным ядром ………………………………………………………………………
3.2. Решение уравнений Вольтерра 2-го рода с вырожденным ядром.
Понятие итерированного ядра и резольвенты. Решение интегрального уравнения Фредгольма второго рода с помощью резольвенты...
Характеристические числа и собственные функции. Теоремы Фредгольма…………………………………………………………………….
5.1. Собственное число (значение) и собственный вектор матрицы...
5.2. Характеристические числа и собственные функции……………..
5.3. Теоремы Фредгольма……………………………………………….
Физические примеры …………………………………………………..
Связь интегральных уравнений с дифференциальными ……………..
Список рекомендуемой литературы……………………………………..
32
3
4
5
5
7
10
10
14
16
19
19
20
23
26
29
31
Download