UMF-BOOK
.pdfМ.А. Греков Уравнения математической физики
то есть |
|
|
|
|
|
|
|
|
|
|
|
|
|
Z |
∂E(x, ξ) |
|
|
Z |
|
∂E |
|
|
|
|
|||
|
|
d − |
|
|
|
d 2 |
= 0, |
x Ω, |
(2.5) |
||||
|
∂ν |
∂r |
|||||||||||
|
|
|
|
|
|ξ−x|=ε |
|
|
|
|
||||
Z |
|
∂E(x, ξ) |
− Z |
∂E |
|
|
|
|
|||||
|
|
|
d 1 |
|
d 2 |
= 0, |
x , |
(2.6) |
|||||
|
∂ν |
|
∂r |
||||||||||
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
(минус перед вторыми интегралами в (2.5) и (2.6) из-за изменения направления внешней нормали 2 на противоположное. ) Соотношения (2.5), (2.6) равносильны следующим:
Z |
∂E(x, ξ) |
|
Z |
|
2 |
|
|
|
|||||
|
|
|
|
d = − |
|
|
|
= −|S1 |
|, x Ω, |
(2.7) |
|||
|
∂ν |
|
εm−1 |
||||||||||
|
|
|
|
|
|
|ξ−x|=ε |
|
|
|
|
|
|
|
|
Z |
|
∂E(x, ξ) |
d 1 |
= − Z |
|
d 2 |
= 0, |
x . |
(2.8) |
|||
|
|
|
|
|
|
||||||||
|
|
|
∂ν |
|
εm−1 |
||||||||
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
Равенство (2.7) подтверждает первую формулу в (2.2). Для подтверждения третьей формулы в (2.2) необходимо в (2.8) устремить ε → 0.
Следует отметить, что равенства (2.7) и (2.8) получены с учетом того, что свойство 3 гармонических функций верно для кусочно гладких поверхностей .
Однако, чтобы получить третью формулу в (2.2) при ε → 0 в (2.8) необходимо ужесточить требование к гладкости поверхности .
А именно, потребуем, чтобы была поверхностью Ляпунова. Тогда x существует касательная плоскость к , то есть существует нормаль к в точке x и 2 стремится к полусфере при ε → 0. Тогда
|
|
|
|
|
|
|
|
|
1 |
|S1|εm−1 |
|
|
|
|
|
|
|
|||
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
2 |
|
|
|
|
(2.9) |
||||||||
lim |
|
|
d |
|
= lim |
|
|
|
|
|
|
= |
|
| |
S |
1| |
. |
|||
εm−1 |
|
|
|
|
εm−1 |
|
||||||||||||||
ε→0 |
Z |
|
2 |
|
ε→0 |
|
|
2 |
|
|
|
|||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из теоремы, или из (2.2) для объемного потенциала |
|
|
|
|
|
|||||||||||||||
|
|
|
u(x) = ZΩ |
E(x, ξ)ρ(ξ)dξ |
|
|
|
|
|
|
|
|||||||||
вытекает формула Гаусса |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z |
∂u(x) |
dγ = −|S1| |
Z |
ρ(ξ)dξ, |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
(2.10) |
|||||||||||||
|
|
∂n |
|
|
|
|
||||||||||||||
|
γ |
|
|
|
|
|
|
|
Ω∩D |
|
|
|
|
|
|
|
|
|
где n нормаль в точке x γ, D Rm любая область с достаточно гладкой границей γ (поверхностью Ляпунова).
141
М.А. Греков |
|
|
|
|
|
|
|
Уравнения математической физики |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Действительно, имеем |
∂nx |
= ZΩ |
|
|
∂nx |
|
|
x γ |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
ρ(ξ) |
dξ, |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
∂u(x) |
|
|
|
|
∂E(x, ξ) |
|
|
|
|
|
|
|
|
|
|
|||||||
Тогда |
Zγ |
∂nx |
dγ = Zγ |
dγ ZΩ |
ρ(ξ) |
|
∂nx |
dξ = ZΩ |
ρ(ξ)dξ Zγ |
|
∂nx |
|
|
dγ = |
|||||||||||||
|
|
∂u(x) |
|
|
|
|
|
|
|
∂E(x, ξ) |
|
|
|
|
|
|
|
|
|
∂E(x, ξ) |
|
||||||
Z |
|
|
|
|
|
dγx + Z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
ρ(ξ)dξ Zγ |
∂nx |
|
|
ρ(ξ)dξ Zγ |
∂nx |
|
dγx = −|S1| |
Z |
ρ(ξ)dξ. |
||||||||||||||||||
|
|
|
|
∂E(x, ξ) |
|
|
|
|
|
|
|
|
∂E(x, ξ) |
|
|
|
|
|
|||||||||
Ω∩D |
|
|
|
|
|
|
Ω1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Ω∩D |
|
|||||
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
{z |
|
|
|
|
} |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
=0 в |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
силу (2.2) |
|
|
|
|
|
|
|
|
|
§ 3. Потенциал двойного и простого слоя
Пусть Ω Rm конечная область с достаточно гладкой границей , а σ(ξ) заданная на вещественная непрерывная функция (достаточно быть суммируемой).
Определение 1. Потенциал двойного слоя масс (зарядов) распределенных по
с плотностью σ называется функция |
|
|
Z σ(ξ) |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
u(x) = |
|
1 |
|
∂E(x, ξ) |
d . |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.1) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|S1| |
∂ν |
|
|
||||||||||||
|
|
Так как E(x, ξ) гармоническая функция при x = ξ и |
∂E |
|
0, то u(x) гармо- |
||||||||||||||||||||
|
|
∂ν −→ |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
ническая функция. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
→∞ |
|
||||||||
|
|
|
|
|
|
|
0. |
|
|
|
|
|
|
|
|
||||||||||
|
|
x |
|
Rm, x , причем u(x) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
−→ |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Обозначим |
|
|
|
|
|
|
|
→∞ |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
Ω+ = Ω, |
|
Ω− = C(Ω) = Rn\(Ω ). |
|
||||||||||||
|
|
При m = 3 выражение под интегралом в (3.1) имеет простой физический смысл. |
|||||||||||||||||||||||
|
|
Пусть ξ0 Ω+, ξ00 Ω− – две симметричные точки относительно точки ξ , |
|||||||||||||||||||||||
лежащие на нормали ν к в точке ξ. |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
Пусть в точках ξ0, ξ00 сосредоточены электрические заряды соответственно −σ0 и |
|||||||||||||||||||||||
+σ0, такие, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
σ0|ξ00 |
− ξ0| = σ(ξ), |
(σ0 > 0) |
|
|
(3.2) |
|||||||||
при |
| |
ξ00 |
− |
ξ0 |
| → |
0. Потенциал поля, созданного этими зарядами в любой точке x = ξ |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|||
равен |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
σ0 |
− |
|
|
σ0 |
|
|
= σ0 [E(x, ξ00) − E(x, ξ0)] . |
(3.3) |
|||||||
|
|
|
|
|
|
| |
ξ00 |
x |
| |
ξ0 |
− |
x |
| |
||||||||||||
|
|
|
|
|
|
|
− | |
|
|
|
|
|
|
|
|
|
|
|
142
М.А. Греков |
Уравнения математической физики |
|
|
Определение 2. Предельное расположение зарядов при |ξ00 − ξ0| → 0 называется диполем, а величины σ(ξ) и ~ν – его моментом и осью соответственно.
По определению производной в данном направлении, очевидно, что
|
lim |
σ |
[E(x, ξ00) |
− |
E(x, ξ0)] = |
|
|
||||
|
|ξ00−ξ0|→0 |
0 |
|
|
|
|
|
|
|||
σ(ξ) lim |
|
E(x, ξ00) − E(x, ξ0) |
|
= σ(ξ) |
∂ |
E(x, ξ) |
− |
(3.4) |
|||
|
|
∂ν |
|||||||||
|ξ00−ξ0|→0 |
|ξ00 − ξ0| |
|
|
|
|
|
потенциал диполя.
3.1Прямое значение потенциала двойного слоя
Будучи гармонической функцией вне потенциал двойного слоя u(x) имеет смысл и на самой . Имеет место теорема
Теорема 9. Если замкнутая ляпуновская поверхность и плотность σ(ξ) непрерывна на , то потенциал двойного слоя (3.1) имеет вполне определенное значе-
ние x0 и является непрерывной функцией точки x0 . |
|
|||||||
|
1 |
Z |
|
∂ |
|
|
|
|
u(x0) = |
|
σ(ξ) |
|
E(x0 |
, ξ)d , x0 |
. |
(3.5) |
|
|S1| |
∂ν |
Доказательство теоремы можно найти в книге [С. Г. Михлин. Курс математической физики. 2002, стр. 365]
Определение 3. Значение потенциала двойного слоя на называется прямым
значением этого потенциала.
Частный случай потенциала двойного слоя – с точностью до множителя интеграл Гаусса
u(x) = |S1|w0(x)
w0(x) = Z |
∂ |
∂ν E(x, ξ)d , σ(ξ) ≡ 1, |
который, как было показано, равен
w0(x) = − |
1 |
|S1|, x . |
2 |
3.2Предельные значения потенциала двойного слоя
На примере интеграла Гаусса видно, что, вообще говоря, потенциал двойного слоя терпит разрыв при переходе через поверхность .
143
М.А. Греков Уравнения математической физики
Обозначим |
Ω+ конечную область, ограниченную поверхностью , а |
|||
Ω− = Rm\{Ω+ S |
)}, а так же |
|
|
|
|
u±(x0) = lim u(x) |
|||
|
x |
→ |
x0 |
|
|
|
|
± |
|
|
x Ω |
|
есть предельные значения u(x) при x → x0 .
Имеет место теорема
Теорема 10. Если замкнутая ляпуновская поверхность и σ(ξ) – плотность, непрерывная на , то для потенциала двойного слоя справедливы предельные соотно-
шения
u+(x0) = − |
1 |
σ(x0) + u(x0), |
|
|||
|
|
|||||
2 |
(3.6) |
|||||
|
1 |
|
|
|
||
u−(x0) = |
σ(x0) + u(x0), x0 |
. |
||||
|
||||||
2 |
Из (3.6) следует, что при переходе через потенциал двойного слоя испытывает скачок
u+(x0) − u−(x0) = −σ(x0), x0 . |
(3.7) |
Теорема 11. При условиях предыдущей теоремы выполнено u(x) u±(x0)
.
Замечание 1. Если σ(ξ) не является непрерывной, но суммируема на , то u±(x0) существует почти всюду и выражаются по формулам (3.6).
Замечание 2. Для нормальных производных потенциала двойного слоя u(x)
верна теорема Ляпунова:
Теорема 12. (Ляпунова.) При тех же условиях, если существует |
|
|||||||||||
x |
|
Ω |
∂u |
|
|
|
± |
|
||||
|
|
|
|
|
∂u |
|
|
|||||
lim |
|
|
= |
|
|
|
|
, |
|
|||
|
→ ± |
∂n |
|
∂n |
|
|
|
|||||
x |
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂u |
+ |
|
∂u |
|
− |
|
|||
|
|
|
= |
|
|
|
(3.8) |
|||||
∂n |
∂n |
в точке x0 .
144
М.А. Греков |
Уравнения математической физики |
|
|
3.3Потенциал простого слоя
Рассмотрим потенциал простого слоя
u(x) = |
1 |
Z |
E(x, ξ)µ(ξ)d |
|
|S1| |
(3.9) |
и его свойства. Имеют место две теоремы.
Теорема 13. Если замкнутая ляпуновская поверхность, а µ(ξ) измерима и ограничена, то потенциал простого слоя (3.9) непрерывен в Rm, включая .
Теорема 14. Если замкнутая ляпуновская поверхность, а µ(ξ) непрерывна на , то на существуют равномерные пределы нормальных производных потенциала простого слоя (3.9) как изнутри, так и извне .
Справедливы формулы |
|
|
|
|
|
|
|
|
|
||
|
∂u±(x0) |
= ± |
1 |
µ(x0) + |
∂u(x0) |
, x0 , |
|||||
|
∂n |
|
2 |
∂n |
|||||||
где |
|
|
|
|
|
|
|
|
|
||
|
|
∂u(x0) |
= |
1 |
Z |
∂E(x0, ξ) |
µ(ξ)d , |
||||
|
|
|
|
|
|
|
|||||
|
|
∂n |
|
|S1| |
|
∂n |
прямое значение нормальной производной потенциала простого слоя
∂u±(x0) = lim ∂u(x). ∂n x→x0 ∂n
x Ω±
(3.10)
(3.11)
(3.12)
Скачок нормальной производной потенциал простого слоя определяется согласно (3.10) по формуле
∂u+(x0) |
− |
∂u−(x0) |
= µ(x0), x0 . |
(3.13) |
|
∂ν |
∂ν |
||||
|
§4. Интегральные уравнения теории потенциала
4.1Сведение задач Дирихле и Неймана к интегральным уравнениям
Рассмотрим замкнутую ляпуновскую поверхность , ограничивающую две области: Ω+внутреннюю (конечную) и Ω− внешнюю.
Поставим одновременно 4 краевые задачи для однородного уравнения Лапласа.
145
М.А. Греков |
Уравнения математической физики |
|||||
|
|
|
||||
Найти функцию u(x), гармоническую в Ω+ или Ω− и удовлетворяющую либо усло- |
||||||
вию задачи Дирихле |
|
|
|
|
|
|
|
|
u| = ϕ(x), |
|
(4.1) |
||
либо условию задачи Неймана |
|
|
|
|
|
|
|
|
∂u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂n |
= ψ(x). |
|
(4.2) |
|
Функции ϕ(x) и ψ(x) будем |
считать непрерывными на |
|
. |
|||
|
|
|
|
|
Обозначим: D+(D−) внутреннюю (внешнюю) задачу Дирихле, N+(N−) внут-
реннюю (внешнюю) задачу Неймана.
Решение задачи Дирихле будем искать в виде потенциала двойного слоя |
|
|||||||
u(x) = |
1 |
Z |
σ(ξ) |
∂E(x, ξ) |
d , |
|
||
|
|
|
(4.3) |
|||||
|S1| |
∂ν |
|||||||
решение задачи Неймана в виде потенциала простого слоя |
|
|||||||
u(x) = |
1 |
|
Z |
µ(ξ)E(x, ξ)d , |
|
|||
|
|
(4.4) |
||||||
|
|S1| |
потребуем, чтобы плотности σ(ξ) и µ(ξ) были непрерывны на .
Представления (4.3), (4.4) являются гармоническими функциями как в Ω+, так и в Ω− и, следовательно, осталось только удовлетворить краевым условиям (4.1), (4.2).
Особый случай задача D−, решение которой на бесконечности может иметь порядок O(|x|2−m), в то время как потенциал (4.3) убывает быстрее и имеет порядок O(|x|1−m). Отсюда следует, что не всякую гармоническую в Ω− функцию можно пред-
ставить в виде (4.3).
Рассмотрим, например, задачу D+. Условие (4.1) следует понимать как |
|
|||||||||||||||||||||||
|
|
|
lim u x |
ϕ x |
|
, |
|
x |
Ω |
+, x |
|
. |
|
|
|
|
(4.5) |
|||||||
|
|
x→x0 |
( ) = |
( 0) |
|
|
|
|
|
0 |
|
|
|
|
|
|||||||||
Так как u(x) потенциал двойного слоя, то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
lim u x |
1 |
σ x |
1 |
|
|
σ |
|
ξ |
|
∂E(x0 |
, ξ) |
d , |
x |
|
|
+, x |
. |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
x→x0 ( |
) = −2 |
( 0) + |S1| Z |
( |
|
) |
|
∂ν |
|
|
|
|
Ω |
0 |
|
(4.6) |
|||||||||
Подставив (4.6) в (4.5) и заменив x0 на x получим интегральное уравнение отно- |
||||||||||||||||||||||||
сительно неизвестной функции σ(ξ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
D+ |
: σ(x) − |
2 |
Z |
σ(ξ) |
∂E(x, ξ) |
d = −2ϕ(x), |
x , |
|
|
|||||||||||||||
|
|
|
|
|
|
(4.7) |
||||||||||||||||||
|S1| |
|
|
∂ν |
|
|
146
М.А. Греков |
|
|
|
|
Уравнения математической физики |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Аналогично получаются интегральные уравнения для остальных трех задач |
|
||||||||||||||||||||
D− : |
σ(x) + |S1| Z |
|
σ(ξ) |
|
∂ν |
d = 2ϕ(x), |
x , |
(4.8) |
|||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
∂E(x, ξ) |
|
|
|
|
||||
N+ : |
µ(x) + |S1| Z |
|
µ(ξ) |
|
∂n |
d = 2ψ(x), |
x , |
(4.9) |
|||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
∂E(x, ξ) |
|
|
|
|
||||
N− : |
µ(x) − |
|
|
2 |
|
Z |
µ(ξ) |
∂E(x, ξ) |
d = −2ψ(x), |
x . |
|
||||||||||
|
|
|
|
|
|
(4.10) |
|||||||||||||||
|S1| |
|
|
∂n |
||||||||||||||||||
Свойства уравнений (4.7) (4.10). |
|
|
|
|
|
|
|
|
|
||||||||||||
Свойство 1. |
|
∂E |
|
|
∂E |
|
|
|
cos(n, r) |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
= |
|
|
|
cos(n, r) = − |
|
|
|
, |
|
|
|||||||
|
|
∂n |
∂r |
|
|
|
rn−1 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
∂E |
|
cos(ν, r) |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
= − |
|
|
. |
|
|
|
||||||
|
|
|
|
|
|
|
|
∂ν |
rn−1 |
|
|
|
Можно показать, что имеет место оценка для поверхности Ляпунова
| cos(n, r)| ≤ Crα,
C = const, α показатель Липшица (Гельдера), 0 < α ≤ 1. В пределе cos(n, r) → 0, r = |x−ξ| → 0. То же для | cos(ν, r)|. Таким образом уравнения (4.7) и (4.10) интегральные
|
|
|
|
|
1 |
|
|
уравнения со слабой особенностью, то есть особенность имеет порядок |
O |
|
, |
||||
rm−1−α |
|||||||
меньший размерности поверхности , то есть интегрируемая особенность. |
|
|
|
||||
Свойство 2. Ядра |
∂ |
E(x, ξ) и |
∂ |
E(x, ξ) получаются одно из другого перестановкой |
|||
|
|
||||||
|
∂ν |
∂n |
|
|
|
точек x и ξ, так как эти ядра вещественны, то они сопряженные. Отсюда следует, что
уравнения (4.7) и (4.10), а также (4.8) и (4.9) попарно сопряженные.
Такие ядра называются союзными, или транспонированными. По определению:
1.Ядра K(x, ξ) и K (x, ξ) = K(ξ, x) называются сопряженными.
2.K(x, ξ) и K(ξ, x) союзные, или транспонированные.
Если K(x, ξ) вещественная функция, то K (x, ξ) = K(ξ, x) = K(ξ, x).
В области вещественных функций (ядер) понятия союзные ядра (уравнения) и сопряженные ядра (уравнения) совпадают.
147
Глава 7
Интегральные уравнения Фредгольма II рода
Речь пойдет об интегральных уравнениях Фредгольма II рода вида:
Z
ϕ(x) − λ K(x, y)ϕ(y)dy = f(x), x Ω( ), (0.1)
Ω( )
где Ω Rm - конечная область, - ее граница, ядро K(x, y) и функция f(x) - заданные вещественные функции точек x, y; ϕ(y) - искомая функция, λ - параметр.
Интегральное уравнение |
|
ϕ0(x) − λ ZΩ( ) K(x, y)ϕ0(y)dy = 0, x Ω( ) |
(0.2) |
называется соответствующим (0.1) однородным уравнением, а однородное уравнение
Z
ψ(x) − λ K(y, x)ψ(y)dy = 0, x Ω( )− (0.3)
Ω( )
союзным (транспонированным). Ядра K(x, y) и K(y, x) по аналогии называются со-
юзными или транспонированными.
Далее будем рассматривать только одномерные интегральные уравнения, т.е. урав-
нения вида: |
|
ϕ(x) − λ Zab K(x, y)ϕ(y)dy = f(x). |
(0.4) |
Отметим, что уравнение с переменным верхним пределом вида: |
|
ϕ(x) − λ Zax K(x, y)ϕ(y)dy = f(x), x > a |
(0.5) |
называется уравнением Вольтерра II рода. В силу линейности интегрального оператора в (0.4) легко видеть, что, если существует общее решение этого уравнения, то оно имеет вид:
Φ(x) = ϕ0(x) + ϕ(x), |
(0.6) |
где ϕ0(x) - общее решение соответствующего однородного уравнения |
|
ϕ0(x) − λ Zab K(x, y)ϕ0(y)dy = 0, |
(0.7) |
а ϕ(x) - частное решение неоднородного уравнения (0.4).
148
М.А. Греков |
Уравнения математической физики |
|
|
§ 1. Метод последовательных приближений решения интегральных уравнений Фредгольма II рода
Если параметр λ удовлетворяет условию |
|
|
|
|
|
|λ| < |
1 |
, |
(1.1) |
|
|
|||
|
M |
|||
где M > 0 и такое, что |
Zab |K(x, y)|dy ≤ M, |
|
|
|
|
a ≤ x ≤ b, |
(1.2) |
то решение ϕ(x) уравнения (0.4) существует и его можно построить методом последо-
вательных приближений.
|
Суть метода. Решение ϕ(x) ищем в виде предела последовательности |
|
||
|
|
ϕ(x) = lim ϕn(x), |
(1.3) |
|
|
|
|
n→∞ |
|
|
ϕ0(x) = f(x), ϕn(x) = f(x) + λ Zab K(x, y)ϕn−1(y)dy, n = 1, 2, . . . |
(1.4) |
||
|
Известно, что сходимость последовательности {ϕn(x)} равносильна сходимости ря- |
|||
да |
∞ |
|
||
|
|
|
|
|
|
|
|
X |
(1.5) |
|
|
ϕ0(x) + |
[ϕn(x) − ϕn−1(x)], |
n=1
так как частичная сумма Sn = ϕn, а значит, Sn → S.
В силу (1.2) имеем оценки, полученные методом математической индукции:
где n = 1, 2, . . . ,
Пусть
Тогда
|ϕn+1
|ϕ0(x)| ≤ m, |ϕn(x) − ϕn−1(x)| ≤ m|λ|nMn,
m = maxx [a,b] f(x). Докажем это.
|ϕn(x) − ϕn−1(x)| ≤ m|λ|nMn.
(x) − ϕn(x)| = λ |
|
ab |
K(x, y)ϕn(x)dx − |
ab |
K(x, y)ϕn−1 |
|
Z |
|
|
Z |
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.6)
(x)dx ≤
≤ λm|λ|nMn |
Z |
ab |
|
= m|λ|n+1Mn+1. |
|
K(x, y)dx| |
|||
|
|
|
|
|
|
|
|
|
|
Таким образом, любой член ряда (1.5) не превышает соответствующего члена положительного числового ряда (мажорируется):
X∞
m|λ|nMn,
n=0
149
М.А. Греков |
Уравнения математической физики |
|
|
который в силу (1.1) сходится. Следовательно, ряд (1.5) и, значит, последовательность непрерывных функций ϕn(x), сходится абсолютно и равномерно к непрерывной функции ϕ(x):
|
∞ |
|
X |
ϕ(x) = nlim ϕn(x) = ϕ0(x) + |
[ϕn(x) − ϕn−1(x)]. |
→∞ |
n=1 |
Переходя к пределу в равенстве (1.4)
Z b
ϕn(x) = f(x) + λ K(x, y)ϕn−1(y)dy,
a
получим, что имеет место (0.4) и, значит, ϕ(x) является решением этого уравнения. Покажем, что это решение единственно. Действительно, пусть кроме ϕ(x) решени-
ем уравнения (0.4) является еще и функция ψ(x). Тогда
θ(x) = ϕ(x) − ψ(x)
есть решение однородного уравнения (0.7), т.е.
Z b
θ(x) = λ K(x, y)θ(y)dy,
a
откуда находим, что
θ0 ≤ |λ|Mθ0, θ0 = max |θ(x)|,
x [a,b]
а значит,
|λ| ≥ M1 ,
что противоречит неравенству (1.1), если θ0 =6 0. Тогда пусть θ0 = 0, значит, θ(x) ≡ 0, т.е. ϕ(x) ≡ ψ(x), что и требовалось доказать.
Аналогичные рассуждения, касающиеся интегрального уравнения Вольтерра II ро-
да, позволяют построить решение для него в виде |
|
ϕ0(x) = f(x), ϕn(x) = f(x) + λ Zax K(x, y)ϕn−1(y)dy, |
(1.7) |
ϕ(x) = lim ϕn(x).
n→∞
§ 2. Теоремы Фредгольма для интегральных уравнений с вырожденным ядром
Определение. Ядро K(x, y) уравнения (0.4) называется вырожденным, если оно
имеет вид
XN
K(x, y) = |
pj (x)qj (y), |
(2.1) |
|
j=1 |
|
где pj (x), qj(y), x [a, b], y [a, b] - заданные вещественные непрерывные функции.
150