Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
MA3.pdf
Скачиваний:
13
Добавлен:
09.02.2015
Размер:
979 Кб
Скачать
☆

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

f(x, y) = x

Теорема 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]