
- •Оглавление
- •Литература
- •Многие задачи механики и физики приводят к исследованию дифференциальных уравнений с частными производными второго порядка. Так, например:
- •Упражнения
- •Примеры
- •Рисунок 3.
- •Пусть требуется найти решение уравнения
- •Предположим, что функция f(x) разложима в равномерно и абсолютно сходящийся ряд Фурье по синусам
- •Рассмотрим краевую задачу для уравнения эллиптического типа

40
46. Найти решение уравнения Лапласа в области, заключенной между двумя кон-
центрическими окружностями радиусов R1 |
и R2 с |
центром |
в начале координат, |
|||||||
удовлетворяющее краевым условиям |
¶u |
|
|
= j (q ), |
u |
|
|
= j |
|
(q ). |
|
|
|
|
|
||||||
|
|
|
|
|||||||
|
¶r |
r =R |
1 |
|
|
r = R2 |
|
2 |
|
|
|
|
|
|
|
||||||
1 |
|
|
|
|
|
|
|
|||
Указание. Решение представить в виде ряда Фурье по cos kq и sin kq. |
15. Метод функций Грина.
Рассмотрим краевую задачу для уравнения эллиптического типа
D(u) = f (M ), |
(15.1) |
|
æ |
|
|
¶u ö |
|
|
|
|
|
|
|
|
|
|||
|
ça u +a |
|
|
|
|
÷ |
= j(M ), |
s = ¶V , |
|
|
(15.2) |
|||||
|
|
|
|
|
|
|
||||||||||
|
è |
1 |
2 |
¶n øs |
|
|
|
|
|
|
a2 |
+a2 |
|
|||
где a |
= a |
(M ), |
a |
2 |
= a |
2 |
(M ), |
a |
,a |
2 |
³ 0, |
¹ 0. |
||||
1 |
1 |
|
|
|
|
|
1 |
|
|
1 |
2 |
|
Функцией Грина называют решение задачи (15.1),(15.2) при специальных значениях функций f и j , именно: f (M ) = -d (M , P) (d - дельта функция), j(M ) º 0. Решение этой задачи, т.е. функцию Грина обозначим черезG(M,P). Если функция Грина найдена, то с ее помощью легко найти и решение исходной задачи (15.1), (15.2). Для этого применим формулу Грина к функциямv=G(M,P) и к исходному решению u(M):
|
æ |
¶u |
ò(GD(u) - uD(G))dv =òçG |
|
|
¶ |
||
V |
è |
n |
S |
|
¶G ö
- u ÷ds. ¶n ø
Поскольку в области V D(u) = f (M ), a D(G) = -d (M , P), то
ò f (M )G(M , P)dvM +òu(M )d (M , P)dvM
V V
Второй интеграл левой части по свойству d - следнее соотношение можно записать в виде
|
ò |
æ |
¶u |
|
¶G ö |
|
= |
çG |
|
- u |
|
÷ds M . |
|
¶n |
|
|||||
|
è |
|
¶n ø |
|||
|
S |
|
|
|
|
|
функции равен u(P). Поэтому по-
|
ò |
æ |
|
¶u |
|
|
¶G ö |
|
|
|
ò |
|
|
|
|
|
|
|
|
||
u(P) = |
çG |
|
|
- u |
|
|
÷ds |
M |
- |
G(M , P) f (M )dv |
M |
. |
|
|
(15.3) |
||||||
|
¶n |
|
|
||||||||||||||||||
|
è |
|
|
|
¶n ø |
|
|
|
|
|
|
|
|
|
|||||||
|
S |
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
Здесь интегрирование производится по координатам точки М. |
|
||||||||||||||||||||
Для первой граничной задачи (a1 º1, a2 º 0) |
G |
|
S = 0, u |
|
S |
= j, из формулы |
|||||||||||||||
|
|
||||||||||||||||||||
|
|
|
|||||||||||||||||||
(15.3) получим решение задачи (15.1),(15.2) |
|
|
|
|
|
|
|
||||||||||||||
u(P) = -òj |
(M ) |
¶G |
- òG(M , P) f (M )dvM . Для второй граничной задачи |
||||||||||||||||||
|
dsM |
||||||||||||||||||||
¶n |
|||||||||||||||||||||
|
|
S |
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|

|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41 |
|
|
|
|
|
|
|
|
|
|||
(a º 0, a |
2 |
º1) |
¶G |
|
|
|
|
|
= 0, |
|
|
¶u |
|
|
= j(M ) и из формулы (15.3) получаем решение |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
1 |
|
|
¶n |
|
|
S |
|
|
|
|
|
¶n |
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
задачи (15.1), (15.2) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
u(P) = òG(M , P)j(M )dsM - òG(M , P) f (M )dvM . |
|
|
|
|
|
|||||||||||||||||||||||||||||
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для третьей граничной задачи (a1 ¹ 0 u a2 ¹ 0) |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
¶G |
|
|
= - |
a1 |
|
G |
|
S |
, |
¶u |
|
|
= - |
a1 |
u |
|
S |
+ |
j(M ) |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
¶n |
|
|
¶n |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
S |
|
|
a |
2 |
|
|
|
|
|
|
S a2 |
|
|
a2 |
|
S |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
j(M ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
В этом случае u(P) = ò |
|
|
|
G(M , P)ds M -òG(M , P) f (M )dvM . |
||||||||||||||||||||||||||||||
a |
2 |
(M ) |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
Таким образом, исходная задача сводится к задаче о нахождении функции Грина.
16. Нахождение функции Грина методом электростатических изображений.
Функция Грина находится явно лишь для областей частного вида. При построении функции Грина полезно воспользоваться следующей ее физической - ин терпретацией. Из курса физики известно, что электрический заряд величиныq, помещенный в точку Р, создает в свободном неограниченном пространстве электростатическое поле, потенциал которого (при определенном выборе системы единиц)
равен: |
u0 (M ) = |
|
q |
|
. |
|
(16.1) |
4prMP |
|
||||||
|
|
|
|
|
|||
Поскольку в замкнутой области V функция Грина G(M,P) имеет вид |
|||||||
|
G(M , P) = |
1 |
|
+ g(M , P), |
(16.2) |
||
|
4prMP |
где g(M,P), как функция точки М, гармоническая в области V и непрерывна вместе с первыми производными в замыкании [V] области V, то первое слагаемое в правой части (16.2) является потенциалом точечного единичного заряда, помещенного в точке Р области V. Второе слагаемое g(M,P) можно также интерпретировать, как потенциал электростатического поля, созданного одним или несколькими зарядами, но расположенными обязательно вне областиV. Это возможно сделать потому, что потенциал электростатического поля, т.е. функция вида (16.1), является гармонической функцией в любой области, свободной от зарядов (т.е. при М¹Р). Заряды вне V надо выбрать так, чтобы они уничтожили на поверхности S действие заряда в точке P, т.е. чтобы выполнялось соотношение G S = 0;эти заряды называют электроста-
тическими изображениями единичного заряда в Р, а сам метод нахождения функции Грина - методом электростатических изображений.
Пример 1. Функция Грина для полупространства. Пусть S есть плоскость z=0, а

42
V - полупространство z>0. Если в точке PÎ V поместить единичный положительный заряд, то его действие на S уничтожится, очевидно, единичным отрицательным
зарядом, помещенным в точке P , которая является зеркальным изображением точки
1
P относительно S. Потенциал поля, созданного зарядом в P , есть
1
g(M , P) = - |
|
1 |
|
|
и, следовательно, функция Грина - |
||||
4prMP |
|
||||||||
|
|
|
|
|
1 |
|
|
|
|
|
1 |
æ |
|
1 |
|
|
1 |
ö |
|
G(M , P) = |
ç |
|
|
- |
÷ |
||||
|
|
|
|
|
|
||||
|
4p |
ç |
|
|
|
|
rMP1 |
÷ |
|
|
è rMP |
|
ø |
||||||
Пример 2. Функция Грина для шара. Пусть V - шар, ограниченный сферой S: |
|||||||||
x2 + y2 + z2 = R2 |
с |
центром в начале координат. Поместим единичный заряд в |
точку P, расположенную внутри сферы S. Покажем, что действие этого заряда наS |
|
может быть уничтожено некоторым зарядом, помещенным в точке P , являющейся |
|
|
1 |
инверсией точки P относительно сферы S; точка P , лежит на прямой ОР вне шара, |
|
причем |
1 |
|
|
rOP × rOP = R2 . |
(16.3) |
1 |
|
Пусть P - произвольно зафиксированная точка сферы S. Рассмотрим два треуголь-
ника OPP и OP P. Эти треугольники подобны, так как они имеют общий угол при
1
вершине О и стороны, образующие этот угол, пропорциональны в силу (16.3). Из
|
|
|
|
|
r |
|
|
|
|
= |
r0 |
, где r0 = rOP ; откуда |
|
|||||
подобия треугольников следует |
PP |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
r |
|
|
|
|
R |
|
|
|||||
|
|
|
|
|
PP |
|
|
|||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||
1 |
|
|
- |
R |
× |
1 |
|
|
|
= 0 |
(16.4) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
4pr |
|
|
|
|
r0 |
4pr |
|
|
|
|
|||||||
PP |
|
PP1 |
|
|
при любом положении точки P на сфере.
Рисунок 7.
Из (16.4) следует, что действие заряда q =1 в Р уничтожается на S зарядом

43
q = - |
R |
, |
помещенном в P . Следовательно, функция Грина - потенциал поля, со- |
|||||||||||
|
||||||||||||||
|
r0 |
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
1 |
æ |
1 |
|
R |
|
1 |
ö |
|
|
зданного |
этими зарядами- |
есть G(M , P) = |
ç |
- |
× |
÷ |
Метод элек- |
|||||||
|
|
|
|
|||||||||||
|
|
|
|
|
4p ç rMP |
|
r0 |
|
rMP |
÷ |
|
|||
|
|
|
|
|
|
è |
|
|
|
|
1 |
ø |
|
тростатических изображений можно применить и в плоском случае, однако физическая интерпретация здесь будет несколько иной. Именно, если на прямой, проходящей через точку Р ортогонально плоскости (x,y), разместить положительные электрические заряды с единичной плотностью, то они создадут плоское поле (т.е. поле, не зависящее от координаты z), потенциал которого (при соответствующем выборе
системы единиц) будет равен u0 |
(M ) = |
1 |
ln |
1 |
. |
2p |
|
||||
|
|
|
rMP |
17. Решение задачи Дирихле для шара.
Зная функцию Грина, можно построить решение задачи Дирихле для уравнения Лапласа. Поскольку в этом случае f (M ) = 0, то искомое решение примет вид
|
ò |
|
|
|
|
|
|
|
|
¶n |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
u(P) = - |
|
|
j |
(M ) |
¶G |
ds M . Поскольку в этом случае |
|
||||||||||||||||||
S |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
1 |
æ |
|
1 |
|
|
R |
|
|
|
1 |
ö |
|
|
|
|
|||||
|
|
|
|
|
|
ç |
|
|
|
|
|
|
÷ |
|
|
|
|
||||||||
G(M , P) = |
|
|
|
|
|
|
|
- |
|
× |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
, то |
|
|
|
||||||||||
|
|
|
|
4p ç rMP |
r0 rMP |
÷ |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
è |
|
|
|
|
|
|
|
|
|
|
|
1 |
ø |
|
|
|
|
|
|
|
|
1 |
|
|
ò |
¶ |
æ |
|
1 |
|
|
R |
|
1 |
ö |
|
|
|
||||||
|
|
|
|
|
ç |
|
|
|
|
÷ |
|
|
|
||||||||||||
u(P) = - |
|
|
|
|
|
|
|
|
|
ç |
|
|
|
- |
|
|
|
× |
|
|
÷j(M )dSM |
, где j(M ) = u |
S |
. Таким образом, |
|
|
4p |
|
|
|
|
|
|
|
|
r0 |
|
|
|||||||||||||
|
|
|
S |
¶n è rMP |
|
|
rMP1 ø |
|
|
|
остается произвести дифференцирование в подынтегральном выражении. Пусть Т - переменная точка, расположенная внутри шара.
|
Рисунок |
8. |
|
|
|
|
|
|
|
|
|
|
|
Из треугольников ОТР и OTP получим r |
|
= (r 2 + r 2 |
|
|
1 |
|
|||||||
PT |
- 2rr |
0 |
cosg ) |
2 |
, |
||||||||
|
|
|
|
1 |
|
0 |
|
|
|
|
|||
|
|
= (r2 |
+ r2 |
1 |
|
|
|
|
|
|
|
|
|
r |
TP |
- 2rr cosg ) |
2 |
, |
|
|
|
|
|
|
|
||
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|

44
где r = rOT , r 0 = rOP , r1 = rOP1 , g - угол при вершине О. Учитывая, что направление внешней нормали к сфере совпадает с направлением радиуса (т.е. с направлением роста r ), получим
|
|
¶ |
æ |
|
|
1 |
|
|
|
R |
|
|
1 |
ö |
|
¶ |
|
æ |
|
1 |
|
|
|
R |
|
1 |
ö |
|
|
|
|||||||
|
|
|
|
ç |
|
|
|
- |
|
|
|
× |
|
|
÷ |
= |
|
|
ç |
|
|
|
|
- |
|
|
|
|
|
× |
|
|
÷ |
|
|
. Учитывая выражения для |
|
|
¶n |
|
|
|
|
|
r |
|
r |
|
¶P |
|
|
|
|
|
r |
|
r |
|
|||||||||||||||||
|
|
ç r |
P M |
|
|
|
0 |
|
|
÷ |
|
ç r |
TP |
|
|
0 |
|
|
÷ |
|
|
|
|||||||||||||||
|
|
M è |
|
|
|
|
|
|
|
MP1 |
ø |
|
|
|
è |
|
|
|
|
|
|
|
|
|
TP1 |
ø |
|
r =R |
|
||||||||
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
R2 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
TP |
и для r |
TP |
, а также, что r = |
подсчитаем правую часть последнего выра- |
||||||||||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
r0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
жения и получим, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
u(P) = |
|
1 |
|
ò |
|
|
|
|
|
R2 - r02 |
|
|
|
|
|
|
|
|
j(M )dSM . Эта формула называется фор- |
||||||||||||||||||
4pR |
|
|
|
|
|
|
- 2Rr |
|
|
|
|
|
|
3 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
cosg )2 |
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
S (R2 + r2 |
0 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мулой Пуассона, а интеграл, стоящий справа - интегралом Пуассона.
Можно проверить, что функция u(P) действительно является решением задачи Дирихле для шара при любой непрерывной функции j(M ). Для этого следует записать подынтегральное выражение в декартовых координатах
R2 - x2 - y2 - z 2
0 0 0 3 j(M ),
((x0 -x )2 + ( y0 -h)2 + (z0 -V )2 )2
где (x0 , y0 , z0 ), (x,h,V ) -соответственно координаты точек P(x0 , y0 , z0 ) и
M (x,h,V).
Непосредственным дифференцированием убеждаемся, что это выражение как функция точки P(x0 , y0 , z0 ) удовлетворяет уравнению Лапласа.