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

attachments_26-06-2014_11-15-27 / Лекц 28 Пар 20 Мет Функ ГРИНА

.pdf
Скачиваний:
13
Добавлен:
30.05.2015
Размер:
388.32 Кб
Скачать
u
n

§20. Метод функций Грина

Мы можем переписать формулы (161), (162) в общем виде

 

 

 

 

 

 

u

 

u x En x y u y dy

u y

 

 

En x y En x y

dS y

(163)

n

 

G

G

 

 

 

n

 

где En x − фундаментальное решение оператора Лапласа, действующего на функции с n независимыми переменными (см. формулы (133) и (134)).

Рассмотрим теперь такую краевую задачу (называемую задачей Дирихле):

u x f x ,

u G u0 x .

С помощью формулы (163) получаем

u x En x y f y dy

 

y

 

 

En x y En x y

u

u0

 

 

dS y

n

 

G

G

 

 

 

 

n

 

 

 

 

 

 

(164)

(165)

Можем сказать, что функция, заданная формулой (165) является решением задачи Дирихле (164).

Формулой (165) пользоваться, однако, затруднительно, так как в поста-

новке задачи (164) не указано явно поведение производной на границе G .

Идея метода функций Грина состоит в том, чтобы заменить функцию En x y на некоторую другую функцию x, y En x y g x, y , называемую функцией Грина. При этом g x, y подбирается таким образом, чтобы в аналоге

формулы (165) исчез интеграл, содержащий производную nu .

Определение 20.1. Функцией Грина для оператора Лапласа называется функция x, y двух переменных x, y G , где G n − ограниченная область, если выполнены условия:

1.x, y является гармонической функцией в области G \ y по перемен-

ной x и x x, y x y ;

2.x, y x G 0;

3.Γ x, y En x y g x, y , где En x − фундаментальное решение оператора Лапласа и g x, y − регулярная и гармоническая функция по пе-

ременной x , т.е. x g x, y 0 .

Предложение 20.2. Кроме перечисленных свойств, для функции Грина автоматически выполняется также следующее свойство:

4. x, y y, x .

Доказательство. (Б/Д) Пусть G n − область. Возьмем точки y и z в

этой области и окружим эти точки шарами радиуса 0 . Выбираем настоль-

ко малым, что U y,

U z, G .

Обозначим G G \

 

 

 

 

 

 

U z,

U y,

и воспользуемся формулой Грина:

 

v u u v dx

 

 

u

u

v

 

 

v

n

 

dS .

(Г3)

 

 

G

 

 

 

G

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

Вэтой формуле положим u x x, z и v x x, y . Заметим, что

x x, y 0 и x x, z 0 согласно свойству 1 функции Грина. Следова-

 

тельно, левая часть формулы (Г3) равна нулю. По свойству 2 x, y

 

x G 0 и

 

 

 

 

 

x, z

 

x G 0 . Поэтому в правой части формулы (Г3) остается только инте-

 

 

 

 

 

 

 

 

грал по границам окрестностей U y,

и U z, в точках y и z , т.е. по сферам

 

S y и S z . Суммируя сказанное, мы приходим к формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ x, z

 

Γ x, y

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ x, y

 

 

 

 

Γ x, z

 

 

 

dSx 0

,

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

 

 

 

 

 

 

 

 

S y S z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

из которой следует равенство (*):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ x, z

 

 

Γ x, y

 

 

 

Γ x, y

 

 

 

 

Γ x, z

 

Γ x, y

 

 

Γ x, z

 

 

 

 

 

dSx Γ x, z

 

 

 

Γ

x, y

 

dSx

n

n

 

 

 

 

n

n

S y

 

 

 

 

 

 

S z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем предел левой части при 0 . По свойству 3

 

 

 

 

 

 

 

 

 

Γ x, y En x y g x, y , подставим вместо x, y это выражение. Заметим,

 

что lim

 

g x, y

Γ x, z

dSx 0 , так как по условию g x, y является гар-

 

 

 

0

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n
Γ x, z

монической и потому непрерывной функцией, а функция непрерывна

во всех точка шара U y, . Кроме того, т.к. непрерывная функция на компакте

ограничена, то существует константа C такая, что

 

x, z

 

C .

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

Следовательно (ограничиваемся случаем n 3 , иные размерности рассматриваются аналогично),

 

En x y

Γ x, z

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

dS

x

 

 

 

 

1

 

 

2

 

0

 

 

 

 

 

 

 

 

dSx

 

C

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

4

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

4

x y

 

 

4

 

 

 

 

 

 

 

 

 

S y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В левой части формулы (*) первый интеграл стремится к нулю:

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ x, y

Γ x, z

dSx 0

при 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S y

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рассмотрим второй интеграл. Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ x, z

Γ x, y

 

 

 

 

 

 

 

 

 

 

 

Γ x, z

1

 

 

 

 

1

 

 

 

 

 

g x, y

 

 

dSx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dSx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

4

n

 

x y

 

 

n

 

 

S y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Γ x, z

 

 

1

 

 

dSx

1

 

 

 

 

 

Γ x, z g dSx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

x y

 

 

 

 

 

 

n

x y

 

 

 

 

 

 

4

 

x y

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как функция непрерывна, а функция g − гармоническая (следовательно,

имеет непрерывную производную

 

g ), то второй интеграл в последнем выра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

жении будут стремиться к нулю при 0 . Далее заметим, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, где вектор n

− нормаль, направлен-

 

n

 

 

x y

 

 

 

 

 

 

 

 

 

r

 

 

 

 

r2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ная внутрь шара по радиусу r

 

x y

 

 

к его центру. Следовательно,

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Γ x, z

 

 

 

 

1

 

 

 

 

dSx

 

1

 

 

 

Γ x, z dSx .

 

 

 

 

 

 

 

 

4

 

 

 

n

 

x y

 

 

4 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как функция x, z на сфере

 

x y

 

непрерывна, то можно воспользо-

 

 

ваться теоремой о среднем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Γ x, z dSx

1

 

Γ

ˆ

 

 

 

2

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 2

 

 

 

 

 

 

4 2

x, z 4

 

Γ x, z ,

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где xˆ − некоторая точка на сфере x y . Тогда

Γxˆ, z Γ y, z при 0 .

Витоге левая часть формулы (*) при 0 стремится к Γ y, z . Аналогично доказывается, что правая часть этой формулы при 0 стремится к Γ z, y .

Теорема 20.3. Пусть G n − ограниченная область, f x C G и

u0 x C G . Тогда решение задачи Дирихле (164) может быть найдено по формуле:

 

 

 

u x Γ x, y f y dy

u0 y

Γ x, y

dS y ,

 

(166)

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Γ x, y

− функция Грина для оператора Лапласа.

 

 

 

 

 

 

 

 

Доказательство. Из пунктов 2 и 3 определения 20.1 следует, что

 

 

 

 

 

 

 

g x, y

 

x G En x y

 

x G

 

 

 

 

 

 

 

(**)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По формуле Грина (Г3) при v g , u f , получаем:

 

 

 

 

 

 

 

 

 

 

 

0 g x, y f x dx

 

 

x, y

u

u

g

 

 

 

 

 

g

n

dSx .

(***)

 

 

 

 

 

G

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В формулах (163) и (***) отдельно сложим правые части и левые. Тогда

u x

 

g x, y E

 

x y

f y

dy

 

u

En x y

 

 

g x, y

dS

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

g x, y

E x y dS.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но последний интеграл в этой формуле равен 0 ввиду формулы (**).

Замечание 20.4. Построение функции Грина сводится к решению задачи

 

 

 

 

 

 

 

y

g x, y 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(167)

 

 

 

 

 

g x, y |

 

 

 

E x y |

 

 

 

 

, x G.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y G n

 

 

 

 

y G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 20.5. Если функция Грина существует, то она единственна.

Доказательство. Допустим, что в области G существуют две функции Грина:

Γ1 x, y En x y g1 x, y , Γ2 x, y En x y g2 x, y .

Тогда их разность g1 x, y g2 x, y будет функцией, гармонической внутри

 

 

 

 

 

 

 

 

 

 

 

области G и непрерывной на еѐ замыкании G

, причем на границе G

эта

функция тождественно равна 0. По следствию 19.5,

g1 x, y g2 x, y 0 всю-

ду в области G . Отсюда сразу получаем, что Γ1 x, y Γ2 x, y .

 

 

Теорема 20.6. Функция Грина для области G является неположительной

во всех еѐ точках.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Нам известно, что функция En x y

стремится к

при y x (например, при n 3 E

x y

1

 

1

 

). Функция g x, y

 

 

 

 

 

 

 

n

 

 

4

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гармоническая, в частности, непрерывная функция. Следовательно,

 

 

 

 

 

 

при y x . Для точки x G найдется сфера

 

 

 

 

Γ

 

x, y

 

 

y G;

y x

до-

статочно малого радиуса, на которой Γ x, y C , где C 0 . По свойству 2

функции Грина Γ x, y 0 на границе G . Функция Γ x, y является гармонической в области G \ y; y x как разность двух гармонических функций,

значит для нее выполнен принцип максимума (минимума). Следовательно, функция Γ x, y принимает неположительные значения во всех точках области

G \ y; y x . Устремляя к нулю, приходим к нужному заключению. ■