Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MatAnal2

.pdf
Скачиваний:
91
Добавлен:
20.05.2015
Размер:
1.88 Mб
Скачать
si, то в случае

23. Криволинейные интегралы

225

 

 

Физический смысл криволинейного интеграла первого рода.

Воспользуемся представлением (23.1) и интеграл справа запишем как пре-

дел интегральных сумм. Тогда получим

Z

N

X

= d(Π)→0 i=1 f (x1 i) , . . . , xn i)) si,

f (x) ds

lim

где ξi [si−1, si], si

= si − si−1, Π – разбиение [0, S], т. е. 0 = s0 <

s1 < · · · < sN = S. Поскольку разбиение отрезка [0, S] порождает разбиение кривой на кривые si−1,si (i = 1, . . . , N ) длины неотрицательной функции f выражение f (x1 i) , . . . , xn i)) si можно интерпретировать как приближенное значение массы кривой si−1,si , а

саму интегральную сумму – как приближенное значение массы кривой

с плотностью распределения массы f . Поэтому сам криволинейный

интеграл можно интерпретировать как массу кривой с плотностью распределения массы f .

Пример. Вычислить массу эллипса

x2

+

y2

= 1 (a > 0, b > 0), если

2

2

 

 

 

 

 

 

 

 

 

a

 

b

плотность распределения массы равна f (x, y) = |y|.

Параметрическое представление эллипса можно записать в виде

x = a cos t,

 

y = b sin t

(0 ≤ t ≤ 2π).

Масса кривой выражается криволинейным интегралом

Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R f (x, y) ds.

ds = |r0(t)| dt = q

 

 

 

 

dt = p

 

 

 

 

dt,

[x0(t)]2 + [y0(t)]2

a2 sin2 t + b2 cos2 t

f (x(t), y(t)) = |y(t)| = |b sin t|.

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z f (x, y) ds = b Z0

 

 

 

p

 

 

 

 

 

 

 

 

 

 

2π | sin t| a2 sin2 t + b2 cos2 t dt =

 

π

 

p

 

 

 

 

 

 

 

 

 

 

= b Z02π

| sin t|

 

a2 + (b2 − a2) cos2 t dt =

= 2b Z

 

sin tp

 

dt =

 

a2 + (b2 − a2) cos2 t

0

226

 

 

 

Четвертый семестр

 

 

 

 

 

 

1

 

 

 

 

 

= 2b Z−1 p

 

dz = 4b Z0

1 p

 

dz.

a2 + (b2 − a2) z2

a2 + (b2 − a2) z2

Упражнение. Вычислите самостоятельно последний интеграл. Для этого рекомендуется рассмотреть случаи a < b, a > b и a = b.

23.3Криволинейные интегралы второго рода

Областью в Rn называется открытое связное множество. Замкнутой об-

ластью называется объединение области и ее границы. Пусть область

Ω Rn. Предположим, что каждой точке x Ω поставлен в соответствие вектор F (x) Rn. Это можно записать с помощью n функций в

следующем виде:

F (x) = ¡ϕ1 ¡x1, . . . , xn¢ , . . . , ϕn ¡x1, . . . , xn¢¢ ≡ (ϕ1(x), . . . , ϕn(x)) .

В этом случае говорят, что на Ω задано векторное поле F . Если функции

ϕi ¡x1, . . . , xn¢ (i = 1, . . . , n) непрерывные, то векторное поле F называется непрерывным. Если все ϕi (i = 1, . . . , n) непрерывно дифференцируемы, то поле F называется непрерывно дифференцируемым в области Ω.

Пусть в области Ω Rn задано непрерывное векторное поле F = (ϕ1, . . . , ϕn), а r = r(t) (α ≤ t ≤ β) определяет кусочно гладкую кривую , лежащую в области Ω. Криволинейным интегралом II рода от векторного поля F вдоль кривой называется интеграл

ZZ β

(F, dr) ≡ ((ϕ1 (x1(t), . . . , xn(t)) , . . . , ϕn (x1(t), . . . , xn(t))) , r0(t)) dt ≡

α

Zβ

1 (x1(t), . . . , xn(t)) x01(t) + · · · + ϕn (x1(t), . . . , xn(t)) x0n(t)] dt ≡

α

Z

 

ϕ1 dx1 + · · · + ϕn dxn.

В частности, в R2 обозначим F = (P (x, y), Q(x, y)), : r = r(t) = (x(t), y(t)) (α ≤ t ≤ β). Тогда получим

ZZ

(F, dr) ≡ P (x, y) dx + Q(x, y) dy =

 

 

23. Криволинейные интегралы

227

 

 

Zβ

=[P (x(t), y(t))x0(t) + Q(x(t), y(t))y0(t)] dt.

α

В этой записи первые два интеграла – суть обозначения криволинейного интеграла II рода, а интеграл справа – определение криволинейного интеграла II рода.

Аналогично, в R3 имеем F = (P (x, y, z), Q(x, y, z), R(x, y, z)), : r = r(t) = (x(t), y(t), z(t)) (α ≤ t ≤ β),

ZZ

(F, dr) =

P (x, y, z) dx + Q(x, y, z) dy + R(x, y, z) dz =

 

 

Zβ

=[P (x(t), y(t), z(t))x0(t) + Q(x(t), y(t), z(t))y0(t)+

α

+ R(x(t), y(t), z(t))z0(t)] dt.

Свойства криволинейных интегралов II рода.

Мы будем рассматривать эти свойства в R2 (изменения на случай произвольного n ≥ 2 очевидны).

1. Криволинейный интеграл II рода не зависит от способа парамет-

ризации кривой.

Доказательство. Если : r = r(t) (α ≤ t ≤ β) и : ρ = ρ(τ ) (a ≤ τ ≤ b), то t = t(τ ), t(a) = α, t(b) = β и t – кусочно непрерывно дифференцируемая функция переменной τ . Тогда

ZZ β

(F, dr) = [P (x(t), y(t))x0(t) + Q(x(t), y(t))y0(t)] dt =

 

α

Zb

=[P (x(t(τ )), y(t(τ )))x0(t(τ )) + Q(x(t(τ )), y(t(τ )))y0(t(τ ))] t0(τ ) dτ =

a

Z b Z

= [P (¯x(τ ), y¯(τ )) x¯0(τ ) + Q (¯x(τ ), y¯(τ )) y¯0(τ )] dτ = (F, dρ),

a

где r(t) = (x(t), y(t)) α ≤ t ≤ β), ρ(τ ) = (¯x(τ ), y¯(τ )) (a ≤ τ ≤ b).

Замечание. Напомним, что r = r(t) и ρ = ρ(τ ) – два параметрических задания кривой , если они, кроме того, что t = t(τ ), сохраняют ориентацию кривой. Другими словами, доказательство свойства 1 имеет силу

228 Четвертый семестр

лишь в том случае, когда представления r = r(t) и ρ = ρ(τ ) определяют одну и ту же кривую и имеют одну и ту же ориентацию.

2. При изменении ориентации кривой на противоположную криво-

линейный интеграл II рода меняет знак на противоположный, т. е.

Z Z

(F, dr) = − (F, dr).

Доказательство. Пусть : r = r(t) (α ≤ t ≤ β), : ρ = r(α+β −t) (α ≤ t ≤ β). Тогда ρ0(t) = −r0(α + β − t) и

Z

Z β

(F, dr) =

(F (ρ(t)), ρ0(t)) dt =

α

Zβ

=− (F (r(α + β − t)), r0(α + β − t)) dt =

Z αβ Z

= − (F (r(τ )), r0(τ )) dτ = − (F, dr).

α

3. Криволинейный интеграл II рода аддитивен относительно кри-

вой .

Доказательство этого свойства такое же, как и для интегралов I рода,

и мы его опускаем.

 

 

 

 

Пример 1. Вычислить криволинейный интеграл

y dx − x dy, где

 

2

2

начинается в точке (1, 0) и

– дуга окружности x

 

+ y = 1, которая

 

 

R

заканчивается в точке (0, 1).

Параметрическое представление кривой имеет вид : x = cos t, y =

sin t ¡0 ≤ t ≤

π

. Поэтому

 

 

 

2

¢

π/2

π/2

π

Z y dx − x dy = Z0

[sin t(− sin t) − cos t · cos t] dt = − Z0

dt = −

 

.

2

Пример 2. Пусть в предыдущем примере – отрезок, начинающийся

в точке (1, 0) и заканчивающийся в точке (0, 1). Тогда параметрическое

представление имеет вид : x = 1 − t, y = t (0 ≤ t ≤ 1) и поэтому

Z Z 1

y dx − x dy = [t(−1) − (1 − t) · 1] dt = −1.

 

0

23. Криволинейные интегралы

 

 

 

 

 

 

 

 

229

 

 

 

 

 

 

Пример 3. Вычислить

R

y dx + x dy, где :

x = cos t,

y = sin t

(0 ≤ t ≤ π/2). Имеем

 

 

 

 

 

 

 

 

 

 

Z y dx + x dy = Z0

π/2

[sin t(− sin t) + cos t · cos t] dt =

 

 

= Z0

π/2 ¡cos2 t − sin2 t¢ dt = Z0

π/2 cos 2t dt = 0.

 

 

 

4. Вычислить

R

y dx + x dy, где : x = 1

t, y = t (0

t

1).

Имеем

 

 

 

 

 

 

 

 

 

 

Z Z 1 Z 1

y dx + x dy = [t(−1) + (1 − t) · 1] dt = (1 − 2t) dt = 0.

0 0

Замечание. Как в первых двух примерах, так и в двух последних мы вычисляли интегралы от одних и тех же векторных полей вдоль разных кривых, начала и концы которых совпадают. В первой паре примеров интегралы оказались разными, а во второй паре – одинаковыми. Ниже мы увидим, что совпадение интегралов во второй паре примеров не является случайностью.

23.4Формула Грина

Напомним, что областью в R2 называется открытое связное множество.

Замкнутой областью называется объединение области и ее границы. Для дальнейшего нам понадобится следующая теорема, которую мы принимаем без доказательства.

Теорема Жордана. Каждая простая, непрерывная, замкнутая кривая в R2 разделяет всю плоскость на две области – ограниченную и

неограниченную, причем сама кривая является границей для обеих областей.

Простую непрерывную замкнутую кривую в R2 называют простым контуром. Говорят, что простой контур γ ограничивает область Ω R2, если Ω – ограниченная область, существование которой гарантируется

теоремой Жордана.

230

Четвертый семестр

 

 

Область Ω R2 называется односвязной, если для любого простого контура γ Ω область, ограниченная контуром γ, содержится в Ω.

Введем понятие ориентации простого контура . Будем говорить, что простой контур ориентирован положительно, если при обходе контура ограничиваемая им область остается слева. Противоположно ори-

ентированный простой контур называют отрицательно ориентированным и обозначают через . Иначе говоря, кусочно гладкий простой контур ориентирован положительно, если в каждой точке простого контура ,

в которой существует касательный вектор, он образует с вектором внутренней нормали правую пару.

Замечане. Данное определение ориентации простого контура носит интуитивный характер. Можно дать и строгое формальное определение ориентации, но мы не будем на этом останавливаться.

Следующая важная теорема устанавливает связь между интегралом по замкнутому контуру и двойным интегралом по области, ограниченной этим контуром.

Теорема (формула Грина). Пусть функции P (x, y) и Q(x, y) непрерывно дифференцируемы в односвязной области Ω R2, а простой кусочно гладкий контур Ω ограничивает область G Ω. Тогда спра-

ведливо равенство

 

 

 

 

 

 

Z∂G P (x, y) dx + Q(x, y) dy = Z ZG ·

∂Q

∂P

¸ dxdy,

(23.2)

 

 

 

∂x

∂y

где ∂G – положительно ориентированная граница области G.

Доказательство. Докажем сначала равенство (23.2) для случая, ко-

гда область G элементарна, т. е.

G = {(x, y) : a ≤ x ≤ b, ϕ(x) ≤ y ≤ ψ(x)} =

= {(x, y) : c ≤ y ≤ d, α(y) ≤ x ≤ β(y)} ,

где функции ϕ, ψ, α и β – кусочно непрерывно дифференцируемы.

23. Криволинейные интегралы

231

 

 

 

y

 

d

H

F

 

 

E

 

K

 

 

 

G

 

 

D

 

A

 

c

B

C

 

 

 

x

0

a

b

В этом случае, применяя формулу сведения двойного интеграла к по-

вторному, получим

 

 

 

 

 

∂P

 

b

ψ(x) ∂P

 

Z ZG

 

dxdy = − Za

 

dx Zϕ(x)

 

(x, y) dy =

∂y

 

∂y

Z b Z b Z b

= − [P (x, ψ(x)) − P (x, ϕ(x))] dx = P (x, ϕ(x)) dx− P (x, ψ(x)) dx =

a a a

Z Z Z

= P (x, y) dx + P (x, y) dx + P (x, y) dx+

ABCD DE EF HK

Z Z

+ P (x, y) dx = P (x, y) dx.

KA ∂G

Аналогично имеем

 

Z ZG

∂Q

dxdy = Z∂G Q(x, y) dy.

 

 

∂x

Складывая эти два равенства, получаем (23.2).

Предположим теперь, что область G ограничена кусочно гладкой кривой ∂G и такова, что ее можно разбить кусочно гладкой кривой на две

элементарные области G1

и G2. Введем обозначения i = ∂G ∩ ∂Gi (i =

1, 2). Ясно, что ∂G1 = 1

и ∂G2 = 2 −. Применяя к каждой из

областей G1 и G2 уже доказанную формулу Грина, получаем

 

 

∂Q

∂P

Z

ZG1

·

 

 

¸ dxdy =

∂y

∂x

232

 

 

 

 

 

 

 

 

 

 

 

 

Четвертый семестр

 

 

 

 

 

 

 

 

 

= Z∂G1 P dx + Q dy = Z 1 P dx + Q dy + Z P dx + Q dy

и

 

 

 

 

 

 

 

 

 

 

 

 

¸ dxdy =

 

 

 

 

 

∂Q

 

∂P

 

 

 

Z

ZG2

·

 

 

 

 

 

 

 

 

∂y

 

 

∂x

 

 

= Z∂G2

P dx + Q dy = Z 2

P dx + Q dy + Z P dx + Q dy.

Поэтому

Z

ZG ·

 

 

 

 

 

 

 

¸ dxdy =

 

 

∂Q

 

∂P

 

 

 

 

 

 

 

 

∂y

 

∂x

= Z 1

P dx + Q dy + Z P dx + Q dy + Z P dx + Q dy = Z∂G P dx + Q dy.

Далее, применяя метод математической индукции, можно доказать формулу Грина для случая, когда область G с помощью n − 1 кусочно гладких кривых 1, . . . , n−1 разбивается на n элементарных областей

G1, . . . , Gn. В частности, любой многоугольник можно разбить на конеч-

ное число треугольников, т. е. элементарных областей. Таким образом, формула Грина справедлива и для любого многоугольника.

В общем случае произвольной области G, ограниченной кусочно гладким контуром , формула Грина получается, если контур аппроксимировать замкнутой ломаной или, что то же самое, область G аппроксими-

ровать многоугольником.

Формула Грина может быть обобщена и на случай многосвязных областей. Точнее, пусть область G ограничена внешним контуром и n

контурами γ1, . . . , γn. Тогда формально имеет место равенство (23.2), где

R R Pn R

под символом ∂G понимается + i=1 γi .

Положим в формуле Грина Q(x, y) = x, P (x, y) = −y. Тогда получим

Z Z Z

dxdy = 1 x dy − y dx.

G 2 ∂G

23. Криволинейные интегралы

233

 

 

R R

Но поскольку G dxdy = m(G), то получили формулу, выражающую меру области G, ограниченной кусочно гладким контуром ∂G, через кри-

волинейный интеграл по этому контуру

Z

m(G) = 1 x dy − y dx.

2 ∂G

Пример. Вычислить площадь области, заключенной между кривой

 

3at

 

3at2

(0 ≤ t ≤ 1) и биссектрисой первой четверти.

 

 

x =

 

, y =

 

 

 

 

 

1+t3

 

1+t3

 

 

На данной кривой имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x dy − y dx = 9a2

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 + t3)2

 

 

 

 

 

 

 

 

Если же x = y, то x dy − y dx = 0. Таким образом, получаем

 

 

 

 

 

 

1

 

1

 

1

t2

1

 

1

 

3

 

 

 

 

 

Z0

 

 

 

m(G) =

 

Z∂G x dy − y dx =

 

9a2

 

dt =

 

 

9a2

 

 

=

 

a2.

 

2

2

(1 + t3)

2

6

4

23.5Потенциальные поля

Непрерывное векторное поле (P (x, y), Q(x, y)) называется потенциальным в области G R2, если существует непрерывно дифференцируемая функция U (x, y), заданная на G, такая, что

∂U

(x, y) = P (x, y),

∂U

(x, y) = Q(x, y).

∂x

∂y

 

 

Такая функция U называется потенциалом поля (P, Q). Другими словами, функция U называется потенциалом поля (P, Q), если

dU (x, y) = P (x, y) dx + Q(x, y) dy.

Следующая теорема содержит необходимое и достаточное условие потенциальности поля.

Теорема 1. Для того чтобы непрерывное поле (P (x, y), Q(x, y)) было потенциальным в области G, необходимо и достаточно, чтобы было

выполнено условие

Z

P (x, y) dx + Q(x, y) dy = 0

(23.3)

234 Четвертый семестр

для любой кусочно гладкой замкнутой кривой G.

Доказательство. Необходимость. Пусть поле (P, Q) потенциаль-

но, т. е. пусть существует такая функция U (x, y), что

∂U

(x, y) = P (x, y),

∂U

(x, y) = Q(x, y) ((x, y) G).

 

 

∂x

∂y

Далее, пусть : r = r(t) = (x(t), y(t)) (α ≤ t ≤ β) – произвольная кусочно гладкая, замкнутая кривая, лежащая в G. Тогда

Z Z β

P (x, y) dx + Q(x, y) dy = [P (x(t), y(t))x0(t) + Q(x(t), y(t))y0(t)] dt =

α

β

·

∂U

 

∂U

(x(t), y(t))y0(t)¸ dt =

= Zα

 

(x(t), y(t))x0(t) +

 

∂x

∂y

=Z β d U (x(t), y(t)) dt = U (x(β), y(β)) − U (x(α), y(α)) = 0.

αdt

Последнее равенство справедливо в силу условия r(α) = r(β), т. е. в силу замкнутости кривой .

Достаточность. Пусть выполнено условие (23.3). Покажем сначала,

R

что в этом случае криволинейный интеграл γ P (x, y) dx + Q(x, y) dy зависит лишь от начальной и от конечной точек кривой γ G и не зависит от самой кривой γ, соединяющей эти точки.

Итак, пусть γ1 : r = r1(t) (α ≤ t ≤ β) и γ2 : r = r2(τ ) (a ≤ τ ≤ b)

– две кусочно гладкие кривые, лежащие в G и такие, что r1(α) = r2(a), r1(β) = r2(b). Тогда кривая = γ1 2)является замкнутой, и поэтому, в силу условия (23.3),

Z Z Z

0 = P dx + Q dy = P dx + Q dy + P dx + Q dy =

γ1 (γ2)

= Zγ1

P dx + Q dy − Zγ2

P dx + Q dy.

Отсюда следует, что Rγ1

P dx + Q dy = Rγ2

P dx + Q dy.

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