14 Криволинейные интегралы
Определение 14.1. Кривая в пространстве R3 множество вида
` = {X(t) = (x(t); y(t); z(t)) : t [α; β]}, |
(24) |
ãäå x(t), y(t), z(t) непрерывные функции на [α; β].
Кривая называется гладкой, если существуют непрерывные x˙ = x0(t), y˙ = y0(t), z˙ = z0(t), не обращающиеся в 0 одновременно (ïðè t = α имеются в виду правые производные, при t = β левые).
Кривая называется кусочно-гладкой, если существуют такие tj :
α=t0<t1<. . . <tN =β, что составляющие е кривые {X(t) : t [tj−1; tj]} гладкие.
Определение 14.2. Длина кривой L(`) = супремум длин вписанных в не ломаных
X0X1 . . . Xn, Xi = X(ti), ò.å.
n |
−→ |
|
X |
|
|
L(`) = sup |Xi−1Xi| : α = t0 < t1 < . . . < tn = β, n N .
k=1
Как известно, длина кусочно-гладкой кривой (24) равна
β
Z
p
L(`) = x˙ 2 + y˙2 + z˙2 dt.
α
Будем называть разбиением кривой (24) набор точек
{Pi = (x(ti), y(ti), z(ti)) : α = t0 < t1 < . . . < tm = β};
обозначим кусочки кривой `i = {(x(t), y(t), z(t)) : t [ti−1; ti]}. Диаметр разбиения
m
равен diam{Pi} = max L(`i).
i=1
Определение 14.3. Функция f(x, y, z) интегрируема по кривой ` вида (24), если ` D(f) и существует число J R : ε > 0 δ > 0 : если diam{Pi} < δ, то при любом выборе точек Xi `i имеем
m
S{Xi}(f) = XL(`i)f(Xi) (J − ε; J + ε).
{Pi}
i=1
Такое число J называется криволинейным интегралом I рода функции f по кривой ` и обозначается
Z
J = f(x, y, z)ds.
`
Здесь ds единый символ. Однако, s имеет и самостоятельный смысл. Параметр s, заданный на кривой ` = {X(s) = (x(s); y(s); z(s)) : s [α; β]}, называется натуральным, если t1, t2 [α; β], t1 < t2, выполнено равенство
L{X(s) : t1 ≤ s ≤ t2} = t2 − t1.
54
Теорема 14.1. Пусть ` кусочно-гладкая кривая вида (24), ` D(f) R3; f(x, y, z)
непрерывная скалярная функция. Тогда f интегрируема по `, и е криволинейный интеграл I рода равен
β
ZZ
p
f(x, y, z)ds = f(x(t), y(t), z(t)) x˙ 2 + y˙2 + z˙2 dt. (25)
` α
Доказательство. Интеграл в правой части (25) существует, поскольку функция кусоч- но непрерывна. Сделаем в интеграле замену переменной, введя натуральный пара-
Z t
p
ìåòð s = L({X(τ) : α ≤ τ ≤ t}) = x˙ 2 + y˙2 + z˙2 du. Тогда интеграл примет вид
α
Z L(`)
f x(s(t)), y(s(t)), z(s(t)) ds. (26)
0
Для (26) все интегральные суммы и диаметры разбиений будут такими же, как для левой части (25). 
Свойства криволинейного интегралa I рода
1. |
Линейность: |
Z` |
(cf + kg)ds = c Z` |
f ds + k Z` |
g ds. |
||||
2. |
Аддитивность: если конец кривой `1 совпадает с началом `2, òî |
||||||||
|
|
|
Z |
|
f ds = Z |
f ds + Z |
f ds. |
||
|
|
≤ Z |
`1 `2 |
`1 |
|
`2 |
|
|
|
|
( ) ` |
f ds |
≤ |
` |
|
|
|
|
|
3. |
Оценка: L ` inf f |
|
|
L(`) sup f. |
|
|
|
||
`
4. Независимость от ориентации кривой: если
` = {X(t) : t [α; β]} = {Y (t) : t [a; b]}, причем X(α) = Y (b), X(β) = Y (a),
то значения интеграла получатся одинаковыми.
Пример 14.1. Проинтегрируем функцию по четверти окружности
` = {x2 |
+ y2 = 1; x ≥ 0; y ≥ 0}, параметризовав кривую двумя разными способами. |
|||||||||||||||
1) ` = {x = t, y = √ |
1 − t2 |
, |
z = 0; 0 ≤ t ≤ 1}. Получаем: |
|
|
|
|
|
||||||||
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
√1−t t2 |
|
t · √11 |
t2 dt = −√1 − t2 0 = 1. |
|||||||||
Z x ds = Z t · s12 + |
2 |
dt = Z |
||||||||||||||
` |
0 |
|
|
|
|
− |
|
0 |
|
− |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2) ` = |
{ |
x = cos s, y = sin s, z = 0; 0 |
≤ |
s |
≤ |
π/2 |
} |
. Параметр s натуральный. Тогда |
|||
|
|
0 |
π/2 |
= |
|
|
|
||||
` x ds = |
π/2 cos s ds = − sin s 0 |
1. |
|
|
|
|
|
||||
Z |
|
Z |
|
|
|
|
|
|
|
|
|
55
Очевидно, если s è u два натуральных параметра на `, òî ëèáî u = s+const, ëèáî
u = const − s. В первом случае говорят, что s è u задают одинаковую ориентацию,
во втором случае противоположные ориентации. Если на кусочно-гладкой кривой задана ориентация, то в каждой е точке (кроме нескольких точек излома) задан единичный касательный вектор ~τ(X(s)) = dsd X(s). При cмене ориентации кривой все
вектора ~τ(X) умножаются на (−1). Если на кривой задан параметр t, возрастающий в том же направлении, что и s, то вектор ~τ равен
|
~τ(X(t)) = |
2 |
1 |
2 2 |
y˙˙ |
|
. |
(27) |
|||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
z˙ |
|
|
|
||||
|
|
p |
|
|
|
||||||
|
|
x˙ + y˙ + z˙ |
|
|
|||||||
Определение 14.4. Пусть ` R3 кусочно-гладкая кривая, на которой выбрана |
|||||||||||
ориентация, задаваемая натуральным параметром |
s |
. Пусть ~ |
~ |
||||||||
~ |
~ векторное поле, |
|
~ |
|
|
|
V (x, y, z) = P (x, y, z)i+ |
||||
|
|
. Тогда криволинейным инте- |
|||||||||
Q(x, y, z)j + R(x, y, z)k |
|
|
` D(V ) |
|
|
|
|
|
|||
гралом II рода векторного поля ~ |
по кривой ` называется интеграл |
|
|||||||||
|
V |
|
|||||||||
Z` |
P (x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz = Z` |
(V~ (X), ~τ(X))ds. |
|
||||||||
Физический смысл криволинейного интеграла II рода это работа силы ~
V , дейст-
вующей на точку, движущуюся по кривой ` в направлении возрастания s.
Выведем формулу для вычисления криволинейного интеграла II рода от непрерывного векторного поля по кусочно-гладкой кривой с произвольным параметром
t, таким, что dt/ds > 0. Для этого в формулу (25) подставим выражение (27) для вектора ~τ :
(V~ , ~τ)ds = |
β |
|
2 2 |
2 |
Q |
x˙ 2 |
+ y˙2 + z˙2 |
dt = |
β |
|
|
P x˙ + Qy˙ + Rz˙ dt. |
|||||||
Z |
Z |
|
|
|
P |
|
|
|
Z |
|
x˙ + y˙ + z˙ |
|
R |
|
|
|
|||
|
|
|
(x˙ ; y˙; z˙) |
|
p |
|
|
|
|
` |
α |
p |
|
|
|
|
|
α |
Эта формула мотивирует традиционное обозначение криволинейного интеграла II рода.
56
Свойства криволинейного интегралa II рода
Z Z Z
1. Линейность: ~ ~ ~ ~
(cV + kW , ~τ)ds = c (V , ~τ)ds + k (W , ~τ)ds.
` ` `
2. Аддитивность: если конец кривой `1 совпадает с началом `2 и ориентация кривых `1 è `2 совпадает с ориентацией кривой ` = `1 `2, òî
Z |
(V~ , ~τ)ds = Z |
(V~ , ~τ)ds + Z (V~ , ~τ)ds. |
` |
`1 |
`2 |
Z
3. Оценка:
`
~
(V
) |
≤ |
( ) |
` |
| |
~ |
| |
, ~τ ds |
|
L ` |
sup V . |
|||
|
|
|
|
|
|
|
4. Если поменять ориентацию кривой, то криволинейный интеграл II рода поменяет |
||||||||||||||||||||||||||
çíàê. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ ~ |
|
|
|
|
|
|
||
Пример 14.2. Проинтегрируем векторное поле ~ |
|
по кривой |
` = |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V (x, y) = yi + 2xj |
|
|
|
|
|
||||
{x2 + y2 = 1; x ≥ 0; y ≥ 0}, опять двумя способами. |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
t, y = √ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
1) ` |
= |
|
x |
1 |
|
|
t2, z = 0; 0 |
|
|
t |
|
1 . Соответственно, |
dx = dt, dy = |
|||||||||||||
|
{ |
|
= |
. Получаем: |
− |
|
|
|
|
≤ |
|
≤ |
} |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
−tdt/√1 − t2 |
|
|
|
|
|
|
|
|
|
|
|
√ |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
2 |
|
|
|
π |
|
|||
|
Z |
ydx + 2xdy = Z |
√1 − t2 |
|
|
|
|
|
|
0 |
|
|||||||||||||||
|
+ 2t · |
√ |
− |
|
|
dt = |
3t |
1 − t − arcsin t |
= − |
|
. |
|
||||||||||||||
|
|
|
2 |
4 |
|
|||||||||||||||||||||
|
1 |
|
t2 |
|
||||||||||||||||||||||
|
` |
|
|
0 |
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2) ` = {x = cos s, y = sin s, z = 0; 0 ≤ s ≤ π/2}; dx = − sin s ds, dy = cos s ds.
Ориентация другая.
Z` |
π/2 |
π/2 |
1 + 32 |
ds = 4 . |
||
ydx + 2xdy = Z0 |
(sin s(− sin s) + 2 cos s · cos s)ds = |
Z0 |
||||
|
|
|
|
cos 2s |
|
π |
Теорема 14.2. (Формула Ньютона Лейбница для криволинейного интеграла) :
Пусть Ω R |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
область; Φ скалярная функция на Ω; векторное поле V = rΦ |
||||||||||||||||||
непрерывно. Пусть кусочно-гладкая кривая |
` Ω, на ней задана ориентация, при |
||||||||||||||||||
которой точка A начало `, точка B e конец. Тогда |
|
|
|
|
|||||||||||||||
|
Z` |
( |
z)dx + Q(x,dΦ ) |
dy |
+ |
( |
|
) |
dz |
= Φ( ) − Φ( |
) |
. |
|
|
|||||
|
|
|
P x, y, |
|
y, z |
|
R |
x, y, z |
B |
A |
|
(28) |
|||||||
|
|
|
| |
|
|
{z |
|
|
|
|
|
|
|
} |
|
|
|
|
|
Доказательство. Вычисляем: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
β |
∂x (X(t))x˙ + ∂y |
|
|
|
|
|
|
β |
dtΦ(X(t)) dt = Φ(X(t)) t=α = |
||||||||||
Z |
(X(t))y˙ + ∂z (X(t))z˙ dt = Z |
||||||||||||||||||
α |
∂Φ |
|
|
∂Φ |
∂Φ |
|
|
|
α |
d |
|
|
|
|
t=β |
||||
= Φ(B) − Φ(A).
57
15 Площадь поверхности
Вспомним, что множество Σ R3 называется простой гладкой двумерной поверхностью, если оно является образом связного открытого множества U R2 ïîä действием инъективного отображения ~r : U 7→R3, ò.å.
Σ = ~r(U) = x(u, v); y(u, v); z(u, v) : (u, v) U , |
(29) |
x(u, v)
причeм отображение ~r = y(u, v) имеет непрерывные частные производные, и его
z(u, v)
матрица Якоби
|
|
∂x/∂u |
∂x/∂v |
|
|
∂~r ∂~r |
||||
J~r(u, v) = |
∂y/∂u |
∂y/∂v |
= |
|||||||
|
|
|
||||||||
∂u |
|
∂v |
|
|||||||
|
|
∂z/∂u |
∂z/∂v |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
имеет всюду ранг 2 (ò.å. ∂~r/∂u è ∂~r/∂v линейно независимы в каждой точке (u, v)
U).
Связное множество Σ R3 называется гладкой поверхностью, если A Σ ε >
0 : множество Σ ∩ Uε(A) является простой гладкой поверхностью. |
|
Важнейшим частным случаем гладкой поверхности является |
график функции |
z = h(x, y), имеющей непрерывные частные производные, т.е. множество вида |
|
h = (x, y, z) : (x, y) U, z = h(x, y) . |
(30) |
Любая гладкая поверхность локально представима в виде графика x = f(y, z), èëè
y = g(x, z), èëè z = h(x, y).
Как определить площадь гладкой поверхности? Если действовать по аналогии с кривой, т.е. взять супремум площадей многогранных поверхностей, вписанных в
глакую поверхность, то даже для ограниченного цилиндра получится ∞ (парадокс
½сапог Шварца ). Поэтому определим площадь поверхности не геометрически, а аналитически.
Определение 15.1. Пусть Σ простая гладкая поверхность вида (29), причем множество U квадрируемо, а частные производные ∂x/∂u, . . . , ∂z/∂v непрерывны и ограничены (в этом случае мы будем говорить, что простая поверхность Σ квадрируема). Тогда площадь поверхности Σ равна
S(Σ) = |
|
det J~rT · J~r |
|
dudv = |
s |
|
|
|
|
|
|
|
|
|
dudv. |
(31) |
|
|
|
~rv |
· |
~ru |
~rv |
· |
~rv |
|
|||||||
ZZ |
|
|
|
|
ZZ |
~ru |
· |
~ru |
~ru |
· |
~rv |
|
|
|||
U |
q |
|
U |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Матрица J~rT ·J~r = G(u, v) (матрицa Грама) является симметричной и положительно опрелеленной. Выясним геометрический смысл корня из е определителя.
58
s
|
|
|
|
|
|
|
|
|
|
~ru |
~ru |
~ru |
|
~rv |
|
|
|
|
|
pdetG(u, v) = |
|
·· |
= p|~ru|2 · |~rv|2 − (~ru · ~rv)2 = |
||||||||||||||
|
~rv ·· |
~ru |
~rv |
~rv |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||
= r|~ru| |
2 |
· |~rv| |
2 |
− |
|
|
|
|
|
|
|
|
||||||
|
|
|
|~ru| · |~rv| · cos(~r[u; ~rv) = |~ru| · |~rv| · sin(~r[u; ~rv) = |~ru × ~rv| . |
|||||||||||||||
δ |
2 Если взять на U квадрат {uo ≤ u ≤ uo +δ; vo ≤ v ≤ vo +δ} с малым δ (его площадь |
||||||||
), то его образ на Σ |
|
2 |
|
|
|
|
|
δ ·~ru è |
|
|
|
будет близок к параллелограмму, натянутому на вектора |
|
||||||
увеличения площади при отображении·| |
~r : U 7→Σ. |
p |
|
|
|
||||
δ |
·~rv, площадь которого δ |
|
~ru × ~rv|. Таким образом, |
|
detG(u, v) это коэффициент |
||||
Чтобы можно было пользоваться определением 15.1, необходимо доказать его корректность, т.е. независимость S(Σ) от способа задания квадрируемой гладкой
поверхности Σ через параметры (u, v). Пусть
n o n o
~
Σ = ~r(u, v) : (u, v) U = R(s, t) : (s, t) W .
Тогда u и v дифференцируемы по s и t; обозначим соответствующую матрицу Якоби через T (s, t). Тогда по теореме о дифференцировании сложной векторнозначной фун-
кции имеем J~ (s, t) = J~r(u, v) · T (s, t).
R
S(Σ) = ZZ |
|
|
|
|
det(J~rT (u, v) · J~r(u, v)) dudv = |
||
U |
q |
||
[сделаем замену переменных; теперь u = u(s, t), v = v(s, t)]
ZZ
q
=det(J~rT (u, v) · J~r(u, v)) · |detT (s, t)| dsdt =
W
ZZ
q
=det (T T (s, t) · J~rT (u, v) · J~r(u, v) · T (s, t)) dsdt =
W
= ZZ |
|
|
|
|
det(JRT~ (s, t) · JR~ (s, t)) dsdt, |
||
W |
q |
||
что и требовалось доказать.
Если квадрируемую гладкую поверхность разрезать на несколько частей по кусоч- но-гладким кривым, то е площадь равна сумме площадей е частей (это следует
59
из аддитивности двойного интеграла). Благодаря этому свойству можно корректно определить площадь ограниченной кусочно-гладкой поверхности как сумму площадей е гладких частей.
Посмотрим, как будет выглядеть формула площади для графика (30):
y(u, v) = v |
|
= J~r(u, v) = |
|
|
|
0 |
|
|
1 |
|
; |
|
|
|
|
||||||||||
x(u, v) = u |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
0 |
|
|
|
|
|
|
||||
z(u, v) = h(u, v) |
|
|
|
|
|
∂h/∂u |
∂h/∂v |
|
|
|
|
||||||||||||||
|
|
1 + (∂h/∂u)2 |
∂h/∂u |
∂h/∂v |
|
|
|
|
|
∂h |
|
2 |
|
∂h |
|
2 |
|||||||||
det(G(u, v)) = |
∂h/∂v · ∂h/∂u |
1 + (∂h/∂v· |
)2 |
= 1 + ∂u |
+ |
∂v |
. |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следовательно, переобозначив |
u = x, v = y, мы получаем |
формулу |
|
|
|
|
|||||||||||||||||||
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S( h) = ZZ |
|
|
|
|
|
|
|
|
|
|
dxdy. |
|
|
|
|
|
(32) |
|||||||
|
1 + |
∂x |
|
2 |
+ |
|
∂y |
|
|
2 |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
∂h |
|
|
∂h |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
U |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Подинтегральное выражение в этой формуле имеет такой геометрический смысл: это секанс угла наклона касательной плоскости к плоскости Oxy, или, что то же самое,
угла между нормалью и осью Oz. Действительно,
|
−∂h/∂x |
|
1 |
|
|
|
|
|||
~ |
∂h/∂y |
~ ~ |
|
|
|
|
|
|
|
|
N(x,y,h(x,y)) h k N = |
− 1 |
= cos(Nd, k) = q1 + ∂h/∂x 2 + ∂h/∂y 2 . |
||||||||
цилиндром {(x − 1)2 + y2 |
= 1}. |
|
Σ = {z = p |
|
|
|
} |
|
|
|
|
4 − x |
|
− y |
|
|
|||||
Пример 15.1. Найти площадь части полусферы |
|
|
|
2 |
2 |
|
, вырезаемой |
|||
1) Представим Σ как график:
p
z = 4 − x2 − y2; (x; y) D = {(x − 1)2 + y2 < 1}.
60
Тогда секанс угла наклона равен
s
1 + |
|
x |
|
|
2 |
|
|
y |
|
|
2 |
2 |
|
|
|
|
|
|
+ |
|
|
|
= |
|
|
|
. |
||||
4 − x2 |
− y2 |
|
4 − x2 |
− y2 |
4 − x2 |
− y2 |
|||||||||
|
p |
|
|
|
|
|
p |
|
|
|
|
p |
|
|
|
Вычислим площадь по формуле (32), перейдя к полярным координатам:
S(Σ) = ZZ |
4 x2 |
|
y2 |
= |
π/2 |
dϕ |
2 cos ϕ |
4 |
ρ2 = |
|
|
|
||||||
|
Z |
Z |
|
|
|
|
||||||||||||
|
|
2dxdy |
|
|
|
|
|
|
|
|
|
2ρdρ |
|
|
|
|||
|
p |
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
||
D |
− − π/ |
2 |
|
−π/2 |
|
|
0 |
− |
dϕ = |
|
π/2 |
4(1 − | sin ϕ|)dϕ = 4π − 8. |
||||||
|
= |
|
Z |
|
−2 4 − ρ2 |
|
0 |
|
|
Z |
||||||||
|
|
|
− |
π/2 |
|
p |
|
|
2 cos ϕ |
|
|
− |
π/2 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2) Параметризуем Σ широтой ψ и долготой ϕ :
x = 2 cos ψ cos ϕ; y = 2 cos ψ sin ϕ; z = 2 sin ψ; (ϕ; ψ) W.
Подставив такие x, y, z в уравнение цилиндра, получаем
W = {|ϕ| < ψ < π/2}.
Дифференцируем x, y, z ïî ψ è ϕ :
~rψ = 2 |
− sin ψ sin ϕ |
; |
~rϕ = 2 |
|
cos ψ cos ϕ |
. |
|
sin ψ cos ϕ |
|
|
|
− cos ψ sin ϕ |
|
|
− cos ψ |
|
0 |
Теперь нетрудно вычислить detG(ϕ, ψ) = |
|
0 |
4 cos2 |
ψ |
= 16 cos2 ψ. Èòàê, |
|||||
|
|
|
|
|
|
4 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
π/2 |
ψ |
π/2 |
|
|
|
|
π/2 |
|
||
S(Σ) = Z |
dψ Z |
4 cos ψdϕ = Z |
8ψ cos ψdψ = 8(ψ sin ψ + cos ψ) 0 |
= 4π − 8. |
||||||
0 |
− |
ψ |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
61
