положителен (здесь и в аналогичных ситуациях cos α, cos β, cos β суть компоненты единичного вектора нормали n, называемые его направляющими косинусами). Для выяснения знака определителя можно рассмотреть его в какой-то одной точке поверхности, т. е. можно взять какие-то конкретные значения параметров u, v и посчитать определитель для этих значений, в других точках знак будет тот же.
Пересадим форму w в пространство параметров, т. е. сделаем замену переменной в форме ω и перейдем к форме ϕ ω. Делается это так: на места переменных x, y, z ставятся значения x(u, v), y(u, v), z(u, v) соответствующих координатных функций параметризации ϕ, на места dx, dy, dz ставятся дифференциальные формы первого порядка dx(u, v), dy(u, v), dz(u, v), представляющие собой дифференциалы координатных функций, после чего раскрываются скобки во внешних произведениях. Получается форма второго порядка вида g(u, v) du dv на области D в R2, интеграл от которой есть интеграл от функции g по области D. В случае, если параметризация не согласована с заданной ориентацией, то перед интегралом ставится знак минус. Проделав намеченные действия для согласованной с ориентацией параметризацией ϕ можно получить выражение для интеграла:
ZZ Z
ωa2 = P (x, y, z) dydz + Q(x, y, z) dzdx + R(x, y, z) dxdy
S |
S |
Z Z |
|
|
|
|
|
|
|
Z |
|
|
u |
v |
|
||||
2 |
|
|
|
|
|
y |
(u, v) y |
(u, v) |
|
= ϕ ωa |
= |
D |
|
P (x(u, v), y(u, v), z(u, v)) |
zu′′ |
(u, v) zv′′ |
(u, v) |
||
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
|
v |
|
|
+Q(x(u, v), y(u, v), z(u, v)) |
|
z (u, v) |
z |
(u, v) |
|
|
xu′′ (u, v) xv′′ |
(u, v) |
|||||
|
|
|
|
|
|
|
|
|
u |
|
v |
|
|
|
|
xu′ |
(u, v) |
xv′ (u, v) |
||
|
|
|
(u, v) |
|
|
|
+ R(x(u, v), y(u, v), z(u, v)) y′ |
y′ (u, v) dudv. |
|||||
|
|
|
|
|
|
|
Укажем связь между поверхностными интегралами первого и второго рода. Пусть S поверхность, ориентированная вектором нормали n, и пусть ωa2 = P (x, y, z) dydz+Q(x, y, z) dzdx+R(x, y, z) dxdy 2-форма, порожденная полем a. Тогда
Z |
Z |
|
ωa2 = |
ha|ni dS, |
(25.5) |
SS
11
или, подробнее,
Z Z
P (x, y, z) dydz + Q(x, y, z) dzdx + R(x, y, z) dxdy
S
Z
=(P (x, y, z) cos α + Q(x, y, z) cos β + R(x, y, z) cos γ) dS,
S |
(25.6) |
|
где n = (cos α, cos β, cos γ).
25.6. Задачи. Вычислить интегралы
Z Z
(1) (2z −x) dydz +(x+2z) dzdx+3z dxdy, S верхняя сторона
S
треугольника x + 4y + z = 4, x > 0, y > 0, z > 0;
Z Z
(2) yz dydz + zx dzdx + xy dxdy, S внутренняя сторона по-
S
верхности тетраэдра x + y + z 6 1, x > 0, y > 0, z > 0;
Z Z
(3) y dzdx, S внешняя сторона сферы x2 + y2 + z2 = R2;
S
Z Z
(4) (x5 + z) dydz, S внутренняя сторона полусферы x2 +
S
y2 + z2 = R2, z 6 0;
Ответы. (1) 128/3; (2) 0; (3) 4πR3/3; (4) −2πR7/105.
§26. Формулы Грина, Стокса, Гаусса Остроградского
26.1.Будем отталкиваться от общей формулы Стокса для гладких многообразий с краем. Начнем с определения используемой в этой формуле операции дифференцирования форм, подробно определив дифференцирование для 1-форм в R3 и сопроводив определение словами, вполне пригодными для проведения дифференцирования форм любого порядка.
Пусть ω1 = f(x, y, z)dx одна из базисных 1-форм в R3. Определим форму dω1 следующим образом: возьмем дифференциал df(x, y, z), рассмотрим его как 1-форму, разложенную по базисным формам, т. е. запишем
df(x, y, z) = fx′ (x, y, z) dx + fy′ (x, y, z) dy + fz′ (x, y, z) dz,
12
и положим
df(x, y, z) = df(x, y, z) dx
= (fx′ (x, y, z) dx + fy′ (x, y, z) dy + fz′ (x, y, z) dz) dx. (26.1)
Если раскрыть скобки в последнем выражении формулы (26.1), то получаем разложение формы dω1 по базисным 2-формам:
df(x, y, z) dx = fz′ (x, y, z) dz dx + fy′ (x, y, z) dx dy.
Аналогично определяется дифференцирование других базисных 1-форм, дифференциал произвольной 1-формы задается как сумма дифференциалов базисных 1-форм в ее разложении по базису. Запишем результат дифференцирования подробнее:
dωF1 = d(P (x, y, z) dx + Q(x, y, z) dy + R(x, y, z) dz) |
|
|
|
|
|||||||||
= |
∂R |
− |
∂Q |
dydz + |
∂P |
− |
∂R |
dzdx + |
|
∂Q |
− |
∂P |
(26.2)dxdy. |
|
|
|
|
|
|
||||||||
∂y |
∂z |
∂z |
∂x |
∂x |
∂y |
||||||||
Поле, координаты которого получились в разложении формы dωF1 по базисным, называют ротором поля F и обозначают через rot F . В терминах ротора равенство (26.2) можно записать так:
dωF1 = ω2 |
F . |
(26.3) |
rot |
|
|
Общая формула Стокса. Пусть M гладкое ориентированное k-мерное многообразие в Rn, k < n, с краем ∂M, ориентация которого индуцирована ориентацией многообразия, и пусть ω гладкая форма, заданная на всем многообразии M. Тогда имеет место равенство
ZZ
ω = dω. |
(26.4) |
∂M M
Самый простой случай n = 2, k = 1. В этом случае из общей формулы Стокса получается
Формула Грина. Пусть D область в R2 с (кусочно) гладкой границей γ, при этом обход границы осуществляется так, что для
каждой точки |
(x, y) |
|
γ близлежащая часть области остается слева |
||||||
1 |
|
|
|
|
|
|
|||
от γ. |
Пусть ω |
|
= P (x, y) dx + Q(x, y) dy гладкая 1-форма на D. |
||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂Q |
∂P |
||
Z |
P (x, y) dx + Q(x, y) dy = Z Z |
|
(x, y) − |
|
(x, y) dxdy. |
||||
∂x |
∂y |
||||||||
γ |
|
|
|
|
D |
|
|
|
|
13
Следующий случай n = 3, k = 1.
Классическая формула Стокса. Пусть S ориентированная поверхность в R3 с краем γ, ориентация которого индуцирована ориентацией поверхности S. Пусть ω1 = P (x, y, z) dx + Q(x, y, z) dy +
R(x, y, z) dz 1-форма на S. Тогда
Z
P (x, y, z) dx + Q(x, y, z) dy + R(x, y, z) dz
γ
= Z Z |
∂R |
− |
∂Q |
dydz + |
∂P |
− |
∂R |
dzdx + |
|
∂Q |
− |
∂P |
dxdy. |
|
|
|
|
|
|
||||||||
∂y |
∂z |
∂z |
∂x |
∂x |
∂y |
||||||||
S |
|
|
|
|
|
|
|
|
|
|
|
|
(26.5) |
|
|
|
|
|
|
|
|
|
|
|
|
|
Формулу (26.5) называют (классической) формулой Стокса. Используя связь между интегралами первого и второго рода, формулу Сток-
са можно записать так: |
Z Z |
|
|
P |
|
Q |
|
|
R |
|
|
|
|
|
Z |
|
|
β |
|
γ |
|
|
|||||||
γ |
S |
|
cos α |
cos |
cos |
|
|
|
||||||
P dx + Q dy + R dz = |
|
|
|
∂ |
|
|
∂ |
|
|
∂ |
|
|
dS |
(26.6) |
|
|
|
|
|
|
|
||||||||
|
|
|
|
∂x |
|
∂y |
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
(в формуле (26.6) определитель |
в последнем интеграле понимается |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
так: надо раскрывать его по первой строке, на пустые места в операции взятия производных ставить соответствующие функции из третьей строки и совершать нахождения производных)
Можно записать формулу Стокса с использованием операции ро-
тора векторного поля: |
Z |
ωF1 = Z Z |
|
|
|
|
ωrot2 |
F . |
(26.7) |
γS
Этот вид формулы иногда сопровождается такими словами: циркуляция векторного поля по ограничивающему поверхность контуру равен потоку ротора этого поля через поверхность.
Последний частный случай общей формулы Стокса в R3 представляется значениями n = 3, k = 2.
Формула Гаусса Остроградского. Рассмотрим трехмерное многообразие в R3 с краем, или, иначе говоря, замкнутую область D, ограниченную (кусочно) гладкой поверхностью S, снабженной индуцированной ориентацией (выбирается внешняя нормаль к поверхности). Пусть на D задана 2-форма ωa2 = P (x, y, z) dydz + Q(x, y, z) dzdx + R(x, y, z) dxdy. Ее дифференциал равен
dωa2 = dP (x, y, z) dy dz + dQ(x, y, z) dz dx + dR(x, y, z) dx dy
= |
∂P |
+ |
∂Q |
+ |
∂R |
dxdydz. (26.8) |
|
|
|
||||
∂x |
∂y |
∂z |
14
Общая формула Стокса в рассматриваемом случае примет вид
Z Z
P (x, y, z) dydz + Q(x, y, z) dzdx + R(x, y, z) dxdy
S
= Z Z Z |
∂P |
+ |
∂Q |
+ |
∂R |
dxdydz. |
(26.9) |
|
|
|
|||||
∂x |
∂y |
∂z |
D
Равенство (26.9) называют формулой Гаусса Остроградского. Для векторного поля a = (P (x, y, z), Q(x, y, z), R(x, y, z)) скаляр-
ное поле (т. е. функцию) ∂P∂x + ∂Q∂y + ∂R∂z , называют дивергенцией поля
a = (P (x, y, z), Q(x, y, z), R(x, y, z)) и обозначают символом div a. В терминах дивергенции формула Гаусса Остроградского выглядит так: Z Z Z Z
ωa2 = div a dxdydz. (26.10)
SD
Формула (26.9) (или (26.10)) с физической точки зрения говорит о том, что поток векторного поля через ограничивающую тело поверхность равен объемному интегралу от дивергенции этого поля.
26.2.Пункт с примером применения формулы Грина.
26.3.Задачи. Применяя формулу Грина, вычислить криволинейный интеграл по замкнутой кривой γ, пробегаемой так, что ее
внутренность остается слева: |
|
|
|
|
|
|
Z (xy+x+y) dx+(xy+x−y) dy, (a) γ эллипс |
x2 |
|
y2 |
|
(1) |
|
+ |
|
= 1, |
|
a2 |
b2 |
||||
γ
(b) γ окружность x2 + y2 = ax;
Z
(2)(x + y)2 dx − (x2 + y2) dy, γ граница треугольника с вер-
γ
шинами (1, 1), (3, 2), (2, 5);
Z
(3) (y − x2) dx + (x + y2) dy, γ граница кругового сектора
γ |
|
0 < r < R, 0 < ϕ < α 6 π/2; |
p |
Z p |
(4)x2 + y2 dx + y(xy + ln(x + x2 + y2)) dx, γ окружность
γ
x2 + y2 = R2;
15