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

UMF-BOOK

.pdf
Скачиваний:
47
Добавлен:
21.03.2016
Размер:
948.39 Кб
Скачать

М.А. Греков Уравнения математической физики

Будем считать (без ограничения общности), что каждая из систем {pj(x)},

{qj (y)}, j = 1, N линейно независима.

Рассмотрим интегральное уравнение Фредгольма II рода с вырожденным ядром:

ϕ(x) − λ j=1 Zab pj (x)qj (y)ϕ(y)dy = f(x).

(2.2)

 

N

 

 

 

 

 

X

 

 

 

 

Его можно записать в виде разложения ϕ(x) по системе функций pk(x)

 

 

N

 

 

 

 

 

X

 

 

 

(2.3)

ϕ(x) = f(x) + λ

ckpk(x),

 

k=1

 

 

 

 

где

Zab qk(y)ϕ(y)dy,

 

 

 

 

 

 

 

 

(2.4)

ck =

k = 1, N

неизвестные постоянные. Найдем ck так, чтобы функция (2.3) являлась бы решением

уравнения (2.2). Подставим (2.3) в (2.2). Тогда имеем

λ

N

ckpk(x) − λ

N N

abpj(x)qj (y)λckpk(y)dy+ abpj(x)qj (y)f(y)dy!!= 0,

 

X

 

 

X X

 

Z

 

k=1

 

 

j=1 k=1 Z

 

 

 

N

pj (x)

"cj − λ

N

 

abqj (y)ckpk(y)dy− abqj (y)f(y)dy# = 0,

 

 

X

 

 

X

Z

Z

 

 

j=1

 

 

k=1

откуда, вследствие линейной независимости функций pk(x), k = 1, N, следует, что

cj − λ k=1 ck Zabqj (y)pk(y)dy−Zabqj (y)f(y)dy = 0

 

 

N

 

 

 

 

X

 

 

 

 

или

N

 

 

 

 

 

 

 

 

 

 

X

 

 

 

(2.5)

 

cj − λ αjkck = γj,

j = 1, N,

 

k=1

 

 

 

 

где

Zabqj (y)pk(y)dy, γj =

Zabqj (y)f(y)dy.

 

αjk =

(2.6)

Таким образом, задача отыскания решения ϕ(x) интегрального уравнения (2.2)

редуцирована (сведена) к решению системы (неоднородных) алгебраических уравнений (2.5).

Точно так же соответствующее однородное уравнение

Z b XN

ϕ0(x) − λ

pj (x)qj (y)ϕ0(y)dy = 0

(2.7)

a

j=1

 

151

М.А. Греков

Уравнения математической физики

 

 

 

 

 

 

 

 

N

 

 

 

 

сводится равенством ϕ0(x) = λ

 

ck0pk(x) к однородной линейной алгебраической си-

стеме

k=0

 

 

 

 

X

N

 

 

X

 

cj0 − λ

αjkck0 = 0, j =

1, N

.

(2.8)

k=1

Как известно из линейной алгебры, система (2.5) имеет единственное решение при любой правой части, если выполнено

det M(λ) =6 0,

где

 

1 − λα11

−λα12

 

 

 

...

...

M =

 

−λα21

1 − λα11

 

 

 

 

λαN1

λαN2

. . .

−λα1N

. . .

−λα2N

.

.

 

.

 

.

.

 

 

.

. . .

1 − λαNN

(2.9)

.

Так как det M(λ) - полином степени N, то существует конечное число m ≤ N значений λ1, . . . , λm нулей этого полинома, и, следовательно, для этих λ условие (2.9)

нарушается. Эти значения λi, i = 1, m называются собственными значениями ядра K(x, y), все остальные λ - правильные.

Таким образом, для любого λ < ∞, λ =6 λi, i = 1, m система (2.5) имеет единственное решение {c1, . . . , cN }. Подставив это решение в (2.3), получим решение инте-

грального уравнения (2.2). Итак, доказана теорема:

Первая теорема Фредгольма. Если λ не является собственным числом K(x, y), то интегральное уравнение (2.2) при любой непрерывной правой части f(x) имеет, и

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

Союзное с (2.7) интегральное уравнение в силу (0.3) имеет вид:

N

 

b

 

 

 

 

X a

 

 

ψ0(x) − λ j=1

Z

pj(y)qj(x)ψ0(y)dy = 0, x [a, b].

(2.10)

 

N

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

Так как ψ0(x) = λ

dk0qk(x), то уравнение (2.10) равносильно союзной с (2.8)

однородной системе

k=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

b

 

 

 

a

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

dj0 = −λ k=1 βjkαk = 0,

βjk = Z

pj(y)qk(y)dy= αkj,

 

 

 

 

 

 

 

или

 

 

b

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

dj0 − λ k=1 αkjdk0 = 0, αkj = Z

qk(y)pj(y)dy, j = 1, N,

(2.11)

152

М.А. Греков

Уравнения математической физики

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk0 = Za b pk(y)ψ0(y)dy,

 

 

 

 

 

 

 

 

 

j = 1, N.

 

Из линейной алгебры

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det(I − λA) = det(I − λA)T

 

rank(I − λA) = rank(I − λA)T ,

 

A = {αij}.

 

 

 

 

 

 

 

 

Если λ = λi, i = 1, m и ранг матрицы M(λ) равен r, то однородная система (2.8)

и союзная с ней (2.11) имеют по N − r линейно независимых решений

 

{cl01i, . . . , clN0i } и {dli1, . . . , dlNi

},

 

 

 

 

(2.12)

l = 1, N − r

Подставив эти решения систем (2.8) и (2.11) в правые части формул

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

ϕl0i(x) = λi

clj0ipj (x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

(2.13)

 

 

 

jP

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

ψl0i(x) = λi

=1 dlj0iqj (x),

i = 1, m, l = 1, N − r,

 

получаем по N − r линейно-независимых решений однородного интегрального уравнения (2.7) и (2.10) для каждого собственного числа λi.

Вторая теорема Фредгольма. Однородное интегральное уравнение (2.7) и союзное с ним (2.10) имеют ровно N − r линейно-независимых решений при любом λ = λi

(det M(λ) = 0), где r ранг матрицы M(λi). Функции ϕ0l i(x), l = 1, N − r называются собственными функциями ядра K(x, y), соответствующими собственном

числу λi.

Вернемся к решению неоднородного уравнения (2.2). Соответствующая система (2.5) относительно коэффициентов разложения функции ϕ(x) имеет решение не для любых правых частей, то есть не для любого вектора γ = (γ1, . . . , γN ), если

λ = λi, i = 1, m, m ≤ N.

Для разрешимости системы (2.5) необходимо и достаточно, чтобы вектор γ был ортогонален решению {d01il, . . . , d0Nil} системы уравнений (2.11), союзной системе уравне-

ний (2.8), т. е.

XN

γkdkil = 0, l =

 

 

 

 

(2.14)

1, N − r, i = 1, m,

k=1

 

где

 

γk = Za b f(y)qk(y)dy.

(2.6)

153

М.А. Греков Уравнения математической физики

На основании (2.6), (2.13) и (2.14) получим, что функции f(x) и ψ0i(x) ортогональ-

ны

 

 

 

 

 

 

 

 

 

 

 

l

j=1 djl0iγj = λi

j=1 djl0i

Z

f(x)qj (x)dx = Z

b

"λi

j=1 djl0iqj (x)# f(x)dx =

λi

 

N

N

 

b

 

 

 

 

 

N

 

X

X

a

 

 

 

a

 

 

 

X

 

 

 

= Za b

f(x)ψl0i(x)dx = 0,

(2.15)

 

 

 

 

 

 

 

 

 

 

l = 1, N − r, i = 1, m.

 

Таким образом пришли к формулировке теоремы

Третья теорема Фредгольма. При значениях λ = λi, i = 1, m для разрешимости

интегрального уравнения (2.2) необходимо и достаточно, чтобы правая его часть f(x) была ортогональна ко всем решениям ψl0i(x), l = 1, N − r однородного урав-

нения (2.10), союзного с однородным уравнением (2.7).

§ 3. Понятие итерированного ядра и резольвенты

Вернемся к методу последовательных приближений, рассмотренному ранее при услови-

ях

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

|λ| <

 

,

 

 

(3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

где M :

a

|K(x, y)| dy ≤ M.

 

 

x, y

 

и f

 

x

непрерывны. Если подставить после-

Предполагалось, что функции K

(

)

 

 

R

 

 

 

 

 

(

)

 

 

довательно каждую предыдущую итерацию ϕn−1, ϕn−2, . . . в выражение

 

 

 

ϕn(x) = λ Za b K(x, y)ϕn−1(y)dy + f(x),

 

(3.2)

то в итоге получим

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕn(x) = λ Zb n

λj−1Kj (x, y)f(y)dy + f(x),

 

(3.3)

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

a

j=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

Zb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.4)

 

 

K1 = K(x, y), Kj (x, y) =

 

K(x, y1)Kj−1(y1, y)dy1,

j = 2, 3, . . .

 

 

 

a

 

 

 

 

 

 

 

 

 

 

Определение. Ядро Kj (x, y) называется итерированным или повторным.

154

М.А. Греков Уравнения математической физики

Kj ≤ Mj−1 max |K|.

(3.5)

Повторяя рассуждения, примененные при доказательстве сходимости последова-

тельности (3.2), можно получить, что ряд

 

 

X

 

λj−1Kj(x, y)

 

j=1

 

при условии (3.1) сходится равномерно в квадрате a ≤ x ≤ b, a ≤ y ≤ b.

 

Определение. Его сумма

 

 

X

 

R(x, y, λ) = λj−1Kj (x, y)

(3.6)

j=1

 

называется резольвентой или разрешающим ядром ядра K(x, y)

или инте-

грального уравнения Фредгольма II рода

 

ϕ(x) − λ Za b K(x, y)ϕ(y)dy = f(x).

(3.7)

Так как ϕ(x) = lim ϕn(x), то переходя в (3.3) к пределу, получим

 

n→∞

 

ϕ(x) = λ Za b R(x, y, λ)f(y)dy + f(x).

(3.8)

Заметим, что R(x, y, λ) непрерывная функция (в следствие равномерной сходимости ряда (3.6)) относительно x, y в квадрате x [a, b], y [a, b] и аналитична (то есть разложима в степенной ряд (3.6)) по λ как для вещественных, так и для комплексных значений λ в круге |λ| < M1 .

Поэтому из (3.8) следует, что, в силу непрерывности f(y), непрерывным является решение ϕ(x).

§ 4. Общий случай интегрального уравнения Фредгольма II рода с непрерывным ядром.

Альтернатива Фредгольма

Изучим теперь решение уравнения (3.7), если параметр λ не обязательно удовлетворяет условию (3.1). Известно, что для функции K(x, y), непрерывной в квадрате [a, b] × [a, b],

155

М.А. Греков Уравнения математической физики

для любого ε > 0 существуют такие линейно-независимые системы непрерывных функ-

ций

 

 

 

 

 

 

 

 

 

 

 

 

(4.1)

{pj(x)}, {qj(y)},

 

x, y [a, b], j = 1, N,

что

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

(4.2)

K(x, y) =

 

 

pj (x)qj (y) + Kε(x, y),

 

j=1

 

 

 

 

 

 

 

 

 

где непрерывная функция Kε(x, y) удовлетворяет условию

 

(b − a) |Kε(x, y)| < ε,

 

x, y [a, b].

(4.3)

В качестве систем (4.1) могут служить, в частности, полиномы (на основании тео-

ремы Вейерштрасса). Подставив (4.2) в (3.7), получим:

 

ϕ(x) − λ Zab Kε(x, y)ϕ(y)dy = F (x),

(4.4)

где

 

 

 

 

 

 

 

 

 

 

 

 

F (x) = f(x) + λ Z b XN

pj(x)qj (y)ϕ(y)dy.

(4.5)

 

 

 

 

a j=1

 

 

 

 

 

для любого фиксированного λ, |λ| < ∞ возьмем такое ε > 0, что

 

|

λ

|

<

 

1

<

b − a

.

(4.6)

 

 

 

 

 

 

ε

ε

 

Тогда для уравнения (4.4), в силу (4.3), выполнено условие (3.1) и, следовательно, уравнение (4.4) однозначно обращается.

Обозначив через Kε(x, y, λ) резольвенту ядра Kε(x, y), уравнение (4.4) запишем в

виде

 

 

 

ϕ(x) = F (x) + λ Zab Rε(x, y, λ)F (y)dy.

(4.7)

Заменяя F (x) выражением (4.5), получим

 

 

ϕ(x) = f(x) + λ j=1 Zab pj (x)qj (y)ϕ(y)dy + λ Zab

Rε(x, y, λ)F (y)dy+

 

N

 

 

 

X

 

 

 

 

N

 

 

 

X

 

 

Zab Rε(x, t, λ) · λ j=1 Zab pj (t)qj (y)ϕ(y)dydt.

 

Отсюда

 

 

 

ϕ(x) − λ Z b XN

rj(x)qj (y)ϕ(y)dy = g(x).

(4.8)

a j=1

где Z b

rj(x) = pj(x) + λ Rε(x, t, λ)pj(t)dt,

a

156

М.А. Греков

Уравнения математической физики

 

 

 

g(x) = f(x) + λ Zab Rε(x, y, λ)f(y)dy.

Таким образом, для любого конечного фиксированного λ интегральное уравне-

ние (3.7) эквивалентно интегральному уравнению Фредгольма II рода (4.8) с вырожденным ядром.

Используя три теоремы Фредгольма для интегрального уравнения с вырожденным ядром, приходим к следующей альтернативе:

Альтернатива Фредгольма. Для любого фиксированного конечного λ

1)либо соответствующее уравнению (3.7) однородное уравнение не имеет нетривиального (не равного нулю) решения и тогда уравнение (3.7) всегда имеет и при том единственное решение при любой правой части f(x),

2)либо однородное уравнение имеет нетривиальные решения и тогда как однородное уравнение Фредгольма II рода, так и союзное с ним уравнение имеют одинаковое конечное число линейно-независимых решений. В этом случае для разрешимости неоднородного уравнения Фредгольма II рода (3.7) необходимо и достаточно, чтобы его правая часть f(x) была ортогональна ко всем решениям

союзного с однородным уравнением Фредгольма II рода уравнения, т. е.

 

 

Za b f(x)ψl0(x)dx = 0,

 

 

 

 

 

 

l = 1, L,

 

где ψ0(x), l =

 

все линейно-независимые решения однородного уравне-

1, L

l

 

 

 

 

 

 

 

ния, союзного с однородным уравнением Фредгольма II рода.

 

Определение. Значение λ, для которого однородное уравнение Фредгольма II рода с

любым непрерывным ядром имеет решения φ0(x),

 

 

 

l = 1, L, как и в случае вырож-

 

 

 

l

 

 

 

 

 

 

денного ядра, называется собственным числом ядра K(x, y), а функции φ0(x)

 

 

 

 

 

 

 

 

 

l

- собственными функциями этого же ядра, соответствующими собственному

числу λ.

 

 

 

 

 

 

 

Замечание 1. Заметим, что записав уравнение Фредгольма II рода в виде

 

 

 

ϕ(y) − λ Zab K(y, t)ϕ(t)dt = f(y)

(4.9)

и умножив его на λK(x, y), после интегрирования получим:

 

 

 

ϕ(x) − λ2

Zab K2(x, t)ϕ(t)dt = f2(x),

(4.10)

 

 

f2(x) = f(x) + λ Zab K(x, y)f(y)dy

 

157

М.А. Греков

Уравнения математической физики

 

 

 

 

 

 

Продолжая процесс, получаем

 

 

 

 

 

ϕ(x) − λm Zab Km(x, t)ϕ(t)dt = fm(x),

(4.11)

 

fm(x) = fm−1(x) + λ Zab K(x, t)fm−1(t)dt,

f1(x) = f(x),

Km(x, y) = Zab K(x, t)Km−1(t, y)dt,

K1(x, y) = K(x, y),

m = 2, 3, . . .

Из формулы (4.11) следует:

Утверждение. Если λ и ϕ(x) собственное число и соответствующая ему собственная функция ядра K(x, y), то λm и ϕ(x) собстывенное число и соответствующая ему собственная функция итерированного ядра Km(x, y).

Можно доказать, что верно и обратное утвердение.

Замечание 2. Все изложенное относительно одномерного интегрального уравнения Фредгольма II рода распространяется непосредственно на случай m-мерного уравнения Фредгольма II рода с непрерывным ядром K(x, y) и непрерывной правой частью f(x).

Замечание 3. Более широким классом функций, для которого справедлива теория Фредгольма, является множество функций, суммируемых с квадратом на рассматриваемом измеримом множестве Ω. В случае лебеговой меры для уравнения Фред-

гольма II рода

u(x) − λ ZΩ

 

 

 

K(x, y)u(y)dy = f(x).

(4.12)

Предпологается, что

ZΩ

ZΩ

 

 

 

 

|K(x, y)|2 dxdy < ∞,

(4.13)

 

 

 

ZΩ

|f(x)|2 dx < ∞.

 

Тогда все сказанное относительно непрерывного решения уравнения Фредгольма II рода справедливо в отношении решения уравнения (4.12), принадлежащего классу функций L2(Ω) Z

|u(x)|2 dx < ∞. (4.14)

Ω

Т.е. u(x) квадратично суммируемая функция.

158

М.А. Греков Уравнения математической физики

Замечание 4. Альтернатива Фредгольма верна и в случае ядер со слабой особенно-

стью, т.е. ядер вида

 

 

 

K(x, y) =

K (x, y)

, 0 < α < m0,

(4.15)

|x − y|α

где K (x, y) непрерывная функция по x и y, m0 размерность области Ω (или ее

границы). Справедливость этого вытекает из того факта, что итерированное ядро уравнения (4.11) для ядра (4.15) при достаточно большом m представляет собой непрерывную функцию переменных x и y.

В частности, если в (4.15) α < m0/2, то (4.15) ядро Фредгольма, удовлетворяю-

щее (4.13) и соответствующее уравнение (4.12) с ядром (4.15) является уравнением Фредгольма II рода.

Замечание 5. Частным случаем уравнения Фредгольма II рода, когда K(x, y) = 0 при

y > x, является уравнение Вольтерра II рода.

 

u(x) − λ Zax K(x, y)u(y)dy = 0.

(4.16)

§ 5. Понятие спектра ядра Фредгольма

Определение. Множество всех собственных чисел ядра K(x, y) называется спектром

этого ядра.

Изучение спектра ядер занимает значительное место в теории интегральных уравнений. Спектр вырожденного ядра уравнения Фредгольма II рода состоит из конечного числа элементов, а спектр ядра уравнения Вольтерра II рода является пустым.Уравнение Вольтерра II рода всегда разрешимо, а соответствующее однородное уравнение имеет только тривиальное решение.

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

Определение. Ядро K(x, y) называется симметричным, если для любых x, y Ω вы-

полнено

K(x, y) = K(y, x).

Если ϕ1(x), ϕ2(x) собственные функции симметричного ядра, соответствующие

отличным друг от друга собственным числам λ1, λ2, то

 

Za b ϕ1(x)ϕ2(x)dx = 0.

(5.1)

159

М.А. Греков Уравнения математической физики

Действительно

Rb

 

b

ϕ1(x) − λ1 K(x, y)ϕ1(y)dy = 0 | × λ2ϕ2(x)

a

R

ϕ2(x) − λ2 K(x, y)ϕ2(y)dy = 0 | × λ1ϕ1(x),

a

Затем проинтегрируем Za b

dx. Получим

R

R

R

b

b

b

2 − λ1) a ϕ1(x)ϕ2(x)dx = λ1λ2 a

ϕ2(x)dx a K(x, y)ϕ1(y)dy−

 

R

R

 

b

b

 

−λ1λ2 ϕ1(x)dx K(x, y)ϕ2(y)dy =

 

a

a

Z b Zb

= λ1λ2 ϕ2(x)dx [K(x, y) − K(y, x)]ϕ1(y)dy = 0,

a

a

откуда и следует равенство (5.1), если только λ1 =6 λ2. Из (5.1) следует, что собственные

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

Действительно, если собственное число λ и соответствующая ему собственная

функция - комплексные, т.е.

λ = λ1 + iλ2, ϕ(x) = ϕ1(x) + iϕ2(x),

то λ и ϕ(x) тоже собственное число и собственная функция соответственно. Если λ2 =6 0, т.е. λ =6 λ, то

Zb Zb

ϕ(x)ϕ(x)dx = |ϕ(x)|2 dx = 0

a a

ϕ(x) ≡ 0, что невозможно в силу определения собственной функции λ2 = 0.

§ 6. Применение теории Фредгольма к внутренней задаче Дирихле

Альтернатива Фредгольма применяется в теории краевых задач для гармонических функций. Рассмотрим внутреннюю задачу Дирихле.Ее решение ищется в виде потенци-

ала двойного слоя

Z σ(ξ)

∂nξ

d = 0.

(6.1)

u(x) = |S1|

1

 

∂E(x, ξ)

 

 

160

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