attachments_26-06-2014_11-15-27 / Лекц 28 Пар 20 Мет Функ ГРИНА
.pdf§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 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
монической и потому непрерывной функцией, а функция непрерывна
во всех точка шара 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 . Устремляя к нулю, приходим к нужному заключению. ■