
Московский_государственный_технический_университет
.pdf
Пусть в области G определено векторное поле |
|
P; Q; R , L – замкнутый |
a |
контур, расположенный в области G, σ – поверхность, ограниченная гладким
или кусочно-гладким контуром L, n0 M cos ,cos ,cos - единичный вектор
нормали на выбранной стороне поверхности σ. Пусть функции P; Q; R
непрерывны вместе со своими частными производными. Тогда справедлива
формула Стокса:
|
|
|
R |
|
Q |
P |
|
R |
|
Q |
|
P |
|
|||
Pdx Qdy Rdz |
y |
|
z |
cos |
z |
|
cos |
x |
|
y |
cos d . (15.9) |
|||||
|
|
|
|
|
|
|
x |
|
|
|
|
|||||
L |
|
|
|
|
|
|
|
|
|
|
Ориентация контура L согласована с ориентацией поверхности σ по правилу
правого винта. (См.рис. 15.4) Или:
|
|
R |
|
Q |
P |
|
R |
||
Pdx Qdy Rdz |
y |
|
z |
dydz |
z |
|
x |
||
|
|
|
|
|
|
||||
L |
|
|
|
|
|
|
|
|
|
Q |
|
P |
(15.10) |
dxdz |
|
dxdy . |
|||
|
|
x |
|
y |
|
Левая часть формулы Стокса – это циркуляция векторного поля вдоль a M
контура L, а правая представляет собой поток через поверхность σ
векторного поля . Формулу Стокса можно записать в векторной форме:
rot a
a, dr
L |
|
|
|
|
|
|
|
rot a, n0 |
d rot a, d . |
||||
|
|
|
|
|
|
|
|
|
|
|
Рисунок 15.4
(15.11)
Физический смысл формулы Стокса:
циркуляция векторного поля a M
вдоль замкнутого контура равна потоку ротора этого поля через произвольную поверхность,
натянутую на этот контур.
Замечание. Формула Стокса остаѐтся справедливой и в случае,
когда поверхность σ является плоской областью, параллельной
51

какой-нибудь координатной плоскости. Тогда формула Стокса превращается
в формулу Грина:
|
|
Q |
|
P |
Pdx Qdy |
x |
|
dxdy . |
|
L |
|
|
y |
Пример 15.3. Проверим, является ли потенциальным векторное поле
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a 3x2 y3 z3i 3z3 x3 y2 |
2y j |
3x3 y3z2 2z k . |
|
|
|
|
|
||||||||||
Решение. По формуле (15.8) имеем: |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
i |
|
j |
|
|
k |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
rot a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
x |
|
y |
|
|
z |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
3x2 y3 z3 |
3x3 y2 z3 2 y |
3x3 y3 z2 2z |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
9x2 y2 z3 9x2 y2 z3 |
|
|
|
|
|
|
||||
9x3 y2 z2 9x3 y2 z2 i |
j 9x2 y2 z3 |
9x2 y2 z |
3 |
k |
0 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ответ: данное поле является потенциальным, т.к. rot a 0. . |
|
||||||||||||||||
Пример 15.4. Вычислим циркуляцию векторного поля |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a 3x |
y z i |
2x 4y z j |
5z 3y k |
по контуру треугольника ABC , где A(1,1,0),
B(0,0,2), C(3,0,1), двумя способами: 1) с
помощью криволинейного интеграла, 2) по формуле Стокса.
Выясним, как зависит циркуляция от расположения контура в данном поле.
Рисунок 15.5
Решение.
1) Вычислим циркуляцию по формуле (15.3):
|
|
|
|
|
|
|
|
|
|
Ц |
|
a, dr |
|
3x y z dx 2x 4y z dy 5z 3y dz |
|
|
|
(см. |
|
|
ABCA |
|
|
|
L |
AB |
|
BC |
CA |
рис.15.5)
2) Составим уравнения прямых, содержащих стороны треугольника.
52

|
|
|
x 1 |
|
|
|
|
|
y 1 |
|
|
|
|
z |
y x |
|
|
|
dy dx |
|
|||||||||
AB 1,1, 2 AB : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
2 |
z 2x 2 |
|
dz |
2dx |
|
||||||
|
|
|
|
x |
|
|
y |
|
|
|
|
z 2 |
|
|
|
y 0 |
|
dy 0 |
|
|
|||||||||
BC 3,0, 1 AB : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
3 0 |
|
|
|
|
|
|
|
|
|
|
x 3z 6 |
|
dz 3dz |
|
|||||||||||
|
|
|
|
x 1 |
|
|
|
|
y 1 |
|
|
|
z |
x 2t 1 |
|
dx 2dt |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
AC 2, |
1,1 AC : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t y t 1; |
|
dy dt . |
|
||||||||
|
2 |
|
|
|
|
|
1 |
|
1 |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z t |
|
|
|
dz dt |
|
Вычислим криволинейные интегралы, преобразовав их в определѐнные.
|
0 |
|
|
|
|
|
|
|
|
|
|
22x 20 dx 9 |
|
|
|
|
|||||
AB |
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
29 |
|
|
|
|
||
|
29z 54 dz |
|
|
|
Ц 9. |
|||||
|
|
4 |
Просуммировав, получим ответ: |
|||||||
2 |
||||||||||
BC |
2 |
|
|
|
|
|
AB BC CA |
|||
|
0 |
19 |
|
|
|
|
|
|||
19t 1 dt |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
||||
|
|
2 1 |
|
|
||||||
|
|
|
|
|
||||||
CA |
1 |
|
|
|
|
|
|
|
|
2) Проверим полученный ответ, вычислив циркуляцию по формуле Стокса:
|
|
|
|
ds . (см. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Ц rot a, n0 |
|
(15.11)). В |
|
качестве |
поверхности |
σ, на которой |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
располагается заданный контур, |
выберем плоскость |
треугольника ABC. |
||||||||||||||||||||||||||||||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
Единичный вектор |
нормали, |
составляющий |
|||||||||||||||||||||||||
n |
BA BC 1, 5, 3 . |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
1, 5, 3 . |
|
|
|
|
|
|
|
|
|
|
|
||||||||
острые углы с осями координат: n0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
35 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
по формуле (15.8) и скалярное произведение: |
|
|
4,1,3 , |
||||||||||||||||||||||||||||||
Найдѐм rot a |
rot a |
|||||||||||||||||||||||||||||||||||||
|
|
|
1 |
|
|
18 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
rota, n0 |
|
|
|
|
|
4 5 9 |
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
35 |
35 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
18 |
|
|
|
18 |
|
|
|
|
|
18 |
|
|
1 |
|
|
|
|
9 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Получим |
Ц rot a, n0 |
ds |
|
|
|
ds |
|
|
|
S ABC |
|
|
|
|
|
|
BA BC |
|
|
|
|
35 9 , |
||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||||||||||||||
|
35 |
|
35 |
|
35 |
35 |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. тот же ответ.
3) Исследуем поведение циркуляции при перемещении нашего контура в пространстве. Для этого используем определение скалярного произведения.
Ц rot a, n0 ds |
rota |
n0 cos ds |
27 cos ds , где – угол между |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
векторами |
|
|
|
|
|
|
|
|
rot a |
и n0 . |
|
|
|
|
|
53

Видим, что циркуляция зависит от этого угла . Отсюда следует, что если треугольник ABC перемещается в пространстве параллельно своему исходному положению, то угол не меняется, и циркуляция по контуру остаѐтся равной 9, а при повороте контура меняется направление вектора что влечѐт изменение циркуляции.
Циркуляция достигнет максимального значения, когда cos 1, т.е. когда
|
|
|
Цmax |
|
|
|
|
|
3 |
|
|
|
|
|
27ds |
27 S ABC |
|
|
|
|
|||||||
rot a |
n0 |
. В этом случае |
|
105 . |
|||||||||
2 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Ответ: Ц=9, Цmax 32 105 .
Занятие 16.
Оператор Гамильтона. Повторные дифференциальные операции.
Ауд.: ОЛ-6 №№ 2382, 2385(б), 2398(б) или ОЛ-5 №№ 10.122, 126, 126, 127, 149 (а, в).
1.Оператор Гамильтона.
Определение. Оператором Гамильтона называют символический вектор
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(набла): |
|
i |
|
|
j |
|
|
k |
или |
|
|
; |
|
; |
|
. |
|
x |
y |
z |
|
|
|||||||||||||
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
x |
|
y |
|
z |
С помощью оператора Гамильтона можно кратко записывать операции над скалярными и векторными полями, доказывать некоторые утверждения векторного анализа. Данный оператор имеет и самостоятельное значение.
2.Дифференциальные операции первого порядка.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1) |
Умножая вектор |
на скалярную функцию U x, y, z , получим grad U : |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
U |
U |
U |
|
|
|||||
|
U |
|
i |
|
|
j |
|
|
|
|
k U |
|
|
i |
j |
k |
grad U . |
|
||||
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
x |
|
y |
|
|
|
z |
|
|
|
|
x |
|
y |
z |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P; Q; R даѐт: |
2) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Скалярное произведение вектора на вектор-функцию a |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
div a |
|
, a |
|
|
P |
|
|
Q |
|
|
R div a . |
|
|
|
|||||||
|
|
|
|
y |
|
z |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
54
|
|
|
|
|
|
|
|
|
|
|
|
|
P; Q; R даѐт |
3) Векторное произведение вектора |
на вектор-функцию a |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
i |
|
j |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rot a |
: |
a |
|
|
|
|
|
|
|
rot a. |
|
|
|
x |
|
y |
|
z |
|
|
|||||||
|
|
|
|
|
P |
|
Q |
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Применять оператор нужно, соблюдая следующие правила:
1) Линейность оператора . Пусть 1, 2 , 3 , . - вещественные или комплексные числа, F,G, H , - скалярные или векторные функции.
Тогда 1F 2G 3H 1 F 2 G 3 H
2) Действие оператора на произведение функций.
|
|
|
|
|
F G H F GH |
F |
|||
|
|
|
|
|
|
|
|
|
G H |
FG H . |
||
|
|
|
|
(Символ «↓» указывает на объект действия оператора).
Докажем, например, что a
div
|
|
|
|
|
|
|
|
|
|
|
|||||||
div a |
b |
a |
b |
a |
b |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
b |
rota |
a |
rotb . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a b |
a ,b b a |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
, b |
b |
, a |
b |
rota |
a |
rotb . |
При доказательстве воспользовались свойством смешанного произведения
|
|
|
|
|
|
|
|
|
|
|
|
|
векторов p a |
, b |
p, |
a |
b . |
|
|
|
|
|
|
||
|
2.Дифференциальные операции второго порядка. |
|
|
|||||||||
Пусть в области G заданы скалярное поле U x, y, z и |
векторное поле |
|||||||||||
|
|
|
|
|
|
|
имеют в области G непрерывные |
|||||
a P; Q; R , такие, что функции U , P, Q, R |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
частные |
производные |
|
второго порядка. |
Тогда |
и |
gradU - |
||||||
|
rot a |
|||||||||||
дифференцируемые векторные поля, |
|
- |
дифференцируемое |
скалярное |
||||||||
div a |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
поле. |
В этих |
случаях к полям |
и |
gradU применимы |
операции |
|||||||
rot a |
55

вычисления дивергенции и ротора, |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
а к полю div a - операция вычисления |
||||||||||||||||||||||||||||
градиента. Таким образом, получим пять повторных операций: |
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
||||
div gradU , rot grad U |
, div rot a |
rot rot a |
, grad div a . |
|
|
|||||||||||||||||||||||
Определение. |
|
Операция |
|
|
|
|
|
|
|
|
|
|
называется |
оператором |
Лапласа и |
|||||||||||||
|
|
|
div grad U |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U . Оператор Лапласа |
|
обозначается |
|
символом |
(лапласиан): div grad U |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
можно записать с помощью оператора набла : |
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
div grad U |
, U |
2U U . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
2 |
2 |
|
|
|||||||
При этом |
2 |
|
|
|
i |
|
|
j |
|
|
|
k |
|
|
|
|
|
|
|
|
|
. |
|
|
||||
|
y |
|
|
|
x |
2 |
y |
2 |
z |
2 |
|
|
||||||||||||||||
|
|
|
|
|
x |
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
2U |
|
|
2U |
|
|
2U |
|
|
|
|
|
|
|
|
||||||
Тогда |
|
U 2U x2 |
|
y2 |
|
z2 . |
|
|
|
|
|
|
|
|||||||||||||||
Определение . |
|
Уравнение |
|
|
в |
|
|
частных |
|
производных |
U 0 или |
|||||||||||||||||
2U |
2U |
2U |
0 называют уравнением Лапласа. |
|
|
|||||||||||||||||||||||
x2 |
y2 |
z2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дважды дифференцируемую функцию |
U x, y, z , |
удовлетворяющую этому |
||||||||||||||||||||||||||
уравнению, называют |
гармонической в области |
G. Для плоской области |
уравнение Лапласа примет вид: 2U 2U 0 .
x2 y2
Примеры гармонических функций:U ax by c , U 3x 2xy 7zx ,
U r , r x2 y2 z2 , r 0 .
4.Запись дифференциальных операций второго порядка в операторной
|
|
|
|
|
форме. |
|
|
|
|
|
|
rot grad U U 0 - потенциальное поле grad U является безвихревым |
|||||
полем. |
|
|
|
|
|
|
|
|
|
|
|
div rot a , |
|
|
|||
, a |
0 - векторное поле роторов соленоидально. |
||||
|
|
|
|
|
|
rot rot a |
grad div a |
a |
56

|
|
|
|
|
|
grad diva |
rot rota |
a ,где |
|
|
|
|
|
|
|
|
|
a P |
x, y, z i Q x, y, z j |
R x, y, z k P Q R |
- вектор-функция, |
||
полученная в результате применения оператора Лапласа функциям P,Q, R . |
|||||
Определение. Векторное поле |
|
|
|||
q M называется гармоническим в области G, |
если в этой области оно одновременно является потенциальным и соленоидальным, т.е. для любой точки гармонического поля выполняются
|
0 |
|
0 . Для односвязных областей из первого условия |
||
условия rot q |
и div q |
||||
следует, что |
|
|
|
|
|
q |
grad , |
а из второго, что div q |
div grad 0 , т. е. |
функция φ, потенциал поля – решение уравнения Лапласа и гармоническая функция.
Пример 17.1. Проверим, является ли гармонической функция
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
U ln |
|
, r x2 |
y2 , r 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Решение. Вычислим частные производные второго порядка функции |
|
||||||||||||||||||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2U |
|
|
2U |
|
|
|
|
|
|
|||
U |
|
ln x2 y2 и вычислим U x2 |
|
y2 . |
|
|
|
|
|
|
|||||||||||||||||||
2 |
|
|
|
|
|
|
|||||||||||||||||||||||
U |
|
|
|
|
|
|
x |
U |
|
|
y |
|
|
2U |
|
x2 y2 |
2U |
|
y2 x2 |
|
|||||||||
|
|
|
|
, |
|
|
|
|
, |
x2 |
|
|
, |
y2 |
|
|
, |
|
|||||||||||
x |
x2 y2 |
y |
x2 y2 |
x2 y2 2 |
x2 y2 2 |
|
|||||||||||||||||||||||
U 2U2 |
2U2 |
0 функция |
U ln |
1 |
- гармоническая. |
Тогда поле |
|||||||||||||||||||||||
r |
|||||||||||||||||||||||||||||
|
|
|
x |
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
x |
|
|
|
y |
|
|
|
|
|
|
2U 2U |
|
|
|
|||||||||
q |
grad U |
|
|
i |
|
|
j -гармоническое и div q |
x2 |
y2 |
0, rot q |
0 |
||||||||||||||||||
x2 |
y2 |
x2 y2 |
57