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

Лекции по вышмату за 1 курс

.pdf
Скачиваний:
110
Добавлен:
15.03.2015
Размер:
1.53 Mб
Скачать

Если в качестве векторной функции взять некоторую переменную си-

ëó ~

~

F

= F (x; y; z), то интеграл можно интерпретировать как работу, ко-

торую совершает эта сила по перемещению материальной точки из точки A (x ( ) ; y ( ) ; z ( )) в точку B (x ( ) ; y ( ) ; z ( )) по заданной кривой l:

87Правило вычисления криволинейных интегралов второго рода и их свойства

1.При заданных условиях на кривую и векторную функцию криволинейный интеграл второго рода существует (без доказательства).

запишем скалярное произведение в

 

 

 

 

 

 

n

 

e

e

 

e

 

 

 

 

 

 

 

 

=

 

P

 

 

 

 

 

2. Рассмотрим интегральную сумму n

 

 

~a

 

x

ti

; y ti

; z

ti

~ri

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~a x ti ; y

 

 

 

 

координатной форме:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ti ; z ti ~ri =

 

 

 

 

 

 

 

 

 

ãäå

 

i

 

 

i

 

 

i

ei

 

e

 

 

i

;ey ti

; z ti

yi

+R x ti

; y ti

 

; z ti

i

e

e

e

 

 

 

 

 

 

e

 

 

 

 

e

 

e

 

 

e

 

e

e

 

= P x t ; y t ; z t

 

x +Q x t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z ;

xi = x (ti) x (ti 1) ; yi = y (ti) y (ti 1) ; zi = z (ti) z (ti 1) :

 

 

 

 

Тогда интегральную сумму можно разбить на три слагаемых n =

 

 

 

x + y + z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

n, ãäå

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = P x t ; y t ; z t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

Xi

e

e

e

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xi

 

 

 

 

 

 

; y ti

; z ti

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ny = Q x ti

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

e

e

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è nz = R x ti

 

; y ti

; z ti

zi:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xi

 

e

e

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что, если каждая из сумм nx; ny; nz имеет предел в выше- указанном смысле, то существует предел общей интегральной суммы и, если обозначить

Ix = lim nx;

n!0

Iy = lim ny;

n!0

Iz = lim nz ;

n!0

òî I = Ix + Iy + Iz:

Числа Ix; Iy; Iz будем называть криволинейными интегралами второго рода от функций, соответственно, P (x; y; z) ; Q (x; y; z) ; R (x; y; z) по кривой l:

271

Замечание 42 Очевидно, что

Таким образом, криволинейный интеграл второго рода общего вида можно разбить на три слагаемых

Zl

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

P (x; y; z) dx+Zl

Q (x; y; z) dy+Zl

R (x; y; z) dz:

Теперь для вычисления интеграла общего вида достаточно научиться

 

вычислять каждое из этих слагаемых.

 

 

 

 

 

 

ется пределом

 

R

 

 

 

 

 

 

 

 

 

 

Рассмотрим интеграл

l P (x; y; z) dx. Как уже было сказано, он явля-

 

 

 

 

интегральной суммы

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

x ;

 

 

 

 

 

x = P x t ; y t ; z t

 

 

 

 

 

 

n

Xi

e e e

 

i

 

 

 

 

 

 

 

i

 

 

i

i

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

ãäå

xi = x (ti) x (ti 1)

, а точка ~

выбирается произвольным обра-

 

 

 

 

 

ti

 

зом на промежутке [ti 1; ti]. Так как функция x = x (t) удовлетворя-

 

ет условиям теоремы Лагранжа на этом промежутке, то существует

 

точка ti 2 [ti 1; ti], в которой x (ti) x (ti 1)

 

= x0 (ti) (ti ti 1) =

 

x (ti) ti:

Если положить ~

 

 

 

 

ti = ti, то интегральная сумма примет вид

n

P x ti

; y ti ; z ti

xi =

n

nx =

P (x (ti) ; y (ti) ; z (ti)) x0 (ti) ti;

X

e e e

Xi

i=1

 

 

 

=1

а эта последняя сумма является интегральной суммой для определенного интеграла

Z

P (x (t) ; y (t) ; z (t)) x0 (t) dt:

Таким образом,

ZZ

P (x; y; z) dx =

P (x (t) ; y (t) ; z (t)) x0 (t) dt:

(86)

Rl Q (x; y; z) dy = R Q (x (t) ; y (t) ; z (t)) y0 (t) dt è Rl R (x; y; z) dz = R R (x (t) ; y (t) ; z (t)) z0 (t) dt:

Формула для интеграла общего вида будет выглядеть следующим образом

Z

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

l

=(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

(87)

Замечание 43 Если кривая лежит в плоскости ОХУ, то она может быть задана уравнением y = f (x) ; x 2 [a; b]. Тогда в качестве

параметра можно взять переменную и параметрические уравнения

272

кривой будут x = t; y = f (t) ; t 2 [a; b]. В этом случае формула (86) для вычисления интеграла будет иметь вид

Z

Z b

 

P (x; y) dx =

P (x; f (x)) dx:

(88)

la

Формула для вычисления интеграла общего вида будет

Z

Z b

P (x; y) dx + Q (x; y) dy = (P (x; f (x)) + Q (x; f (x)) f0 (x)) dx

l

a

(89)

3. Криволинейный интеграл второго рода обладает свойством линейно-

ñòè, ò.å.

Z

~

C1 ~a (x; y; z) + C2 b (x; y; z) d~r:

l

4. Криволинейный интеграл второго рода аддитивен относительно кривой интегрирования, т.е.

Z Z Z

~a (x; y; z)d~r = ~a (x; y; z)d~r + ~a (x; y; z)d~r;

l l1 l2

если кривая l разбита на две непересекающиеся части l1 è l2:

5. Криволинейный интеграл второго рода зависит от направления на кривой. При изменении направления величина интеграла меняет знак,

ò.å.

Z

~a (x; y; z) d~r = Z

~a (x; y; z) d~r:

Это непосредственно следует из определения интегральной суммы для

такого интеграла.

 

 

 

 

 

 

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

Rl

xdy

 

ydx, где l - дуга параболы y = x2 между

точками (0; 0) и (2; 4):

 

 

 

 

Решение Для вычисления воспользуемся формулой (89). Для этого подставим туда y = x2 и заметим, что при движении по кривой переменная

меняется от

= 0 äî

= 2.

 

 

 

 

 

 

2

 

2

2

x2dx = 3:

Zl xdy ydx = Z0

x x2 0 x2 dx = Z0

2x2 x2 dx = Z0

 

 

 

 

 

 

8

 

ñòè (x 1)2

+ (y 1)2 = 1,R

лежащая в полуплоскости x + y 1. Направ-

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

l (x + y) dx + (x y) dy, где l - часть окружно-

ление обхода - против часовой стрелки.

Решение

Запишем уравнения окружности в параметрической форме:

x 1 = cos t

y 1 = sin t

273

èëè x = 1 + cos t; y = 1 + sin t.

 

Интегрирование происходит по той части окружности, которая получа-

 

 

 

 

 

 

 

ется при изменении t от 2

 

äî :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вычисления интеграла воспользуемся формулой (87). Тогда

 

 

 

 

 

 

 

 

 

 

 

Zl

(x + y) dx + (x y) dy = Z /2

((2 + cos t + sin t) ( sin t) + (cos t sin t) cos t) dt =

 

 

 

 

 

= Z /2

2 sin t 2 sin t cos t + cos2 t sin2 t dt = Z /2

( 2 sin t sin 2t + cos 2t) dt =:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2 cos t + 2 cos 2t + 2 sin 2t

 

 

= 2 + 1 = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx + 2yzdy

 

 

x2dz, где l - кривая x =

 

 

 

 

 

 

 

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

l

y2

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t; y = t ; z = t ; t 2 [0; 1]

. Обход кривой в направлении возрастания пара-

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

метра.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение Воспользуемся формулой (87). Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

Zl

y2 z2 dx + 2yzdy x2dz = Z0

 

t4 t6 + 2t5 2t t2 3t2 dt = Z0

2t4

+ 3t6

dt =

:

 

 

 

 

35

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

 

y2

 

z2

 

dx + z2

 

 

x2

 

dy +

2

2 y

2

dz, ãäå l

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

 

2

 

 

ñ

 

 

 

 

 

 

 

- кривая, которая образуется при пересечении сферы

x

+ y + z = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

координатными плоскостями, при условии, что x 0; y 0; z

0. Íà-

 

 

 

 

 

 

 

правление обхода таково, что движение в плоскости ОХУ происходит от

 

 

 

 

 

 

 

 

точки (1; 0; 0) до точки (0; 1; 0):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение Кривая, по которой берется интеграл, является кусочно-гладкой,

 

 

 

 

 

 

поэтому разобьем ее на три гладкие части, каждая из которых лежит в од-

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

ОХУ через l1, часть лежащую в плоскости OYZ - через l2 и часть, лежащую

 

 

 

 

 

 

 

в плоскости OXZ - через l3:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда l1 определяется уравнениями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z = 0

 

 

+ z2

= 1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

+ y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что равносильно уравнению x2 + y2

= 1 при условиях x 0 и y 0 в

 

 

 

 

 

 

 

плоскости ОХУ. Очевидно, что l1 - четверть окружности радиуса 1 с цен-

 

 

 

 

 

 

 

тром в начале координат. Параметрические уравнения этой кривой будут

 

 

 

 

 

 

 

x = cos t; y = sin t; 0 t 2 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь воспользуемся формулой (87):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zl1

y2

z2

dx +

z2 x2

dy +

 

x2

y2

dz = Z0

/2

 

sin2 t ( sin t) + cos2 t

cos t

 

dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/2

sin3 t + cos3 t

dt = Z0

/2 sin2 t d (cos t) Z0

/2 cos2 t d (sin t) = Z0

/2

1 cos2 t

d (cos t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos3 t

 

 

 

 

 

sin3 t

/

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

/2

 

1 sin2 t

d (sin t) =

 

cos t

 

sin t +

0

2

 

 

 

1 +

1

 

4

 

 

 

Z0

 

 

 

 

3

 

3

 

= 1 + 3

3 =

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

274

~a ~ =

Также вычисляются интегралы по кривым l2 è l3, причем значения этих интегралов такие же как и вычисленный интеграл по кривой l1. Поэтому

Z

y2 z2 dx + z2 x2 dy + x2 y2 dz = 4:

l

88Связь между криволинейными интегралами первого и второго рода

Теорема 86 Пусть l - кусочно-гладкая кривая и ~a = ~a(x; y; z) = (P (x; y; z) ; Q (x; y; z) ; R (x; y; z))

- векторная функция, заданная в некоторой области, содержащей кривую l и непрерывная на этой кривой. Обозначим через ~ единичный вектор ка-

сательной к кривой l, причем направление ~ соответствует направлению движения по кривой. Тогда имеет место равенство:

 

 

Zl

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

(~a ~) ds;

 

 

 

(90)

 

 

 

 

 

 

 

где ~a ~ - скалярное произведение векторов ~a и ~:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство Пусть кривая l задана уравнениями x = x (t) ; y =

 

 

 

 

 

 

 

y (t) ; z = z (t) ; t 2 [t0; T ]. Тогда вектор касательной будет иметь вид (x0 (t) ; y0 (t) ; z0 (t))

 

 

 

 

и орт этого вектора равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ = 0

 

 

 

 

x0 (t)

 

 

;

 

 

 

y0 (t)

 

 

 

 

 

;

 

 

 

 

z0 (t)

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(y))

2

2

+ (z0 (y))

2

q(x0 (y))

2

+ (y0 (y))

2

+ (z0

(y))

2

q(x0

(y))

2

 

2

(y))

2

@q(x0

 

+ (y0 (y))

 

 

 

 

 

 

 

 

+ (y0 (y)) + (z0

 

A

Тогда в точках кривой l будет

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) q =

(x0 (t))2 + (y0 (t))2 + (z0 (t))2

P(M (t)) dx + Q (M (t)) dy + R (M (t)) dz

=q

(x0 (t))2 + (y0 (t))2 + (z0 (t))2

è

Zl

 

T

 

 

 

 

 

 

 

 

 

 

 

 

P (M (t)) x0((x0 (t))2 + (y0 (t))02 + (z0 (t))2

0

 

 

q(x0 (t))2

+ (y0 (t))2 + (z0 (t))2dt =

(~a ~) ds = Zt0

(t)

 

 

 

 

t) + Q (M (t)) y (t) + R (M (t)) z

 

 

 

 

 

 

T

 

q

 

 

 

 

 

 

 

 

 

 

= Zt0

(P (M (t)) x0 (t) + Q (M (t)) y0 (t) + R (M (t)) z0 (t)) dt = Zl

P (x; y; z) dx + Q (x; y; z) dy + R (x; y; z) d

 

В последних равенствах были использованы формулы для вычисления

 

 

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

 

 

 

 

 

 

 

 

Замечание 44 Формулу (90) можно записать в виде

 

 

 

 

 

 

 

Zl

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

(P (x; y; z) cos + Q (x; y; z) cos + R (x; y; z) cos ) ds;

где ; ; - углы между касательным вектором ~ и координатными осями.

 

 

275

Замечание 45 Если интегрирование ведется по кривой, лежащей в плоскости ОХУ, то формула (90) примет вид

Z Z

P (x; y) dx + Q (x; y) dy = (P (x; y) cos + Q (x; y) sin ) ds;

l

l

где - угол между касательным вектором и осью ОХ.

в интеграл первого рода. l - виток цилиндрической винтовой линииR

x =

Пример 250 Преобразовать криволинейный интеграл второго рода

l ydx + xdy + adz

R cos t; y = R sin t; z = at:

 

 

 

 

 

 

 

 

 

 

 

 

Решение Вектор касательной к данной кривой равен

 

 

 

(x0 (t) ; y0 (t) ; z0 (t)) = (

R sin t; R cos t; a) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда орт касательной

 

 

+ a2

; pR2 + a2 :

 

 

 

~ = pR2 + a2 ; pR2

 

 

 

 

 

R sin t

R cos t

 

 

a

 

 

 

По теореме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ds:

 

Zl

ydx + xdy + adz = Zl

 

 

 

 

+ xR cos t + a2

 

yR sin tp

 

 

 

 

 

R2 + a2

 

 

Если в последний интеграл подставить уравнения кривой, то оконча-

тельно получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zl

ydx + xdy + adz = p

 

Zl

 

 

 

 

R2 + a2

ds:

 

 

в интеграл второго рода, если и - углы, которыеRlобразует нормаль к

Пример 251 Преобразовать интеграл первого рода

(x cos + y cos ) ds

кривой с осями ОХ и ОУ.

Решение Если нормаль к кривой интегрирования составляет с координатными осями углы и , то будет выполнено одно из двух равенств: либо

+ = 2 , ëèáî = 2 . В том и другом случае будет верно соотношение cos = sin , т.е. координаты орта нормали будут (cos ; sin ) :

Тогда угол между одним из касательных векторов и осью ОХ будет

равен +

и, следовательно, координаты орта этого вектора

2

 

cos +

 

; sin

+

 

= ( sin ; cos ) :

 

 

 

 

2

2

Запишем подынтегральную функцию как скалярное произведение некоторого вектора на орт касательной:

x cos + y cos = x cos + y sin = ( y) ( sin ) + x cos :

Тогда по формуле (90) получим

Z Z

(x cos + y cos ) ds = ydx + xdy:

l l

Если взять касательный вектор, направленный в противоположную сто-

ðîíó, òî

Z Z

(x cos + y cos ) ds = ydx xdy:

l l

276

89 Теорема Грина и ее следствия

Определение 141 Область будем называть односвязной, если для

любого простого контура, лежащего в этой области, та часть плоско- ñòè 1, которая лежит внутри этого контура, целиком содержится в области , т.е. 1 :

Теорема 87 Пусть дана область односвязная и элементарная, как от-

носительно оси ОХ, так и относительно оси ОУ, и простой кусочногладкий контур - граница этой области. Тогда ее можно записать в

âèäå

= f(x; y) =a x b; g1 (x) y h1 (x)g

èëè â âèäå

= f(x; y) =c y d; g2 (y) x h2 (y)g ;

где функции g1 (x) ; h1 (x) ; g2 (y) ; h2 (y) непрерывные и кусочно-непрерывно

дифференцируемые.

 

 

и Q (x; y) непрерывны и

Предположим также, что функции P (x; y)

имеют непрерывные частные производные в области

и на ее границе.

Тогда имеет место формула

@x

@y

dxdy;

 

Z P (x; y) dx + Q (x; y) dy = ZZ

(91)

 

 

@Q

 

@P

 

где обход по контуру совершается так, чтобы область оставалась сле-

âà.

Формула (91) называется формулой Грина.

Доказательство Разобьем двойной интеграл в правой части формулы (91) на два слагаемых и каждое слагаемое превратим в повторный интеграл

 

@P

b

h1(x) @P

ZZ

 

dxdy = Za

 

dx Zg1(x)

 

dy =

@y

 

@y

b

 

b

P (x; g1 (x)) dx Za

b

Za

P (x; y)jy=g1(x) dx = Za

P (x; h1 (x)) d

 

y=h1

(x)

 

 

= Z

P (x; y) dx + Z

P (x; y) dx + Z

P (x; y) dx + Z

P (x; y) dx = Z P (x; y) dx

 

 

R

 

 

R

 

 

Добавленные интегралы P (x; y) dx è

P (x; y) dx равны нулю.

Аналогично доказывается, что

 

@Q dxdy =

 

Q (x; y) dy:

 

 

 

 

@x

 

 

 

 

 

RR

 

 

R

 

Таким образом, сформулированная теорема доказана.

Замечание 46 Если область не элементарна, но может быть разбита

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

(рис) Тогда

ZEF AE P (x; y) dx + Q (x; y) dy = ZZ

@x

@y dxdy;

 

 

@Q

@P

1

 

 

 

 

 

@x

@y dxdy;

ZABCA P (x; y) dx + Q (x; y) dy = ZZ

 

@Q

@P

2

 

 

 

 

 

277

ZEACDE P (x; y) dx + Q (x; y) dy = ZZ

@x

@y dxdy:

 

 

 

@Q

@P

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

RAE

REA è

RAC

 

CA равны нулю, получим требуемое.

 

 

 

 

Суммируя эти равенства и, учитывая, что суммы

 

 

+

 

+

R

Замечание 47 Формулу Грина можно применять и тогда, когда область не является односвязной, но может быть разбита не несколько односвязных областей. Например, пусть область , ограниченная кривыми и

разбита на две односвязные области 1 è 2. Тогда

 

 

Z 1[[ 1[ P (x; y) dx + Q (x; y) dy = ZZ

@x

 

@y dxdy

 

 

 

 

 

 

 

@Q

 

@P

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@y dxdy;

 

 

 

 

Z 2[[ 2[ P (x; y) dx + Q (x; y) dy = ZZ @x

 

 

è

 

 

 

@Q

@P

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

ãäå = 1 [ 2 è = 1 [ 2

 

 

 

 

 

 

RAB

+

RBA

= 0

è

CD +

DC = 0, получим, что

 

 

 

 

 

 

R

Складывая эти равенства и используя тот факт, что

 

 

 

R

 

 

 

 

 

 

 

 

 

 

ZZ

@x

@y

dxdy = Z [ P (x; y) dx + Q (x; y) dy = Z + Z ;

 

@Q

@P

 

где контуры и обходятся так, чтобы область оставалась слева.

Следствие 10 Положим P (x; y) = y; Q (x; y) = x. Тогда формула Грина

RRR

принимает вид xdy ydx =

2dxdy. Деля обе части этого равенство

лу для вычисления

 

RR

 

 

на 2 и, учитывая, что

dxdy равен площади области , получаем форму-

 

 

 

 

 

 

площади с помощью криволинейного интеграла второго

ðîäà:

 

S = 2 Z xdy ydx

(92)

 

 

 

1

 

 

Следствие 11 Если функции P (x; y) и Q (x; y) удовлетворяют условиям теоремы Грина и область односвязна, то для того, чтобы выполня-

R

лось равенство l P (x; y) dx + Q (x; y) dy = 0, где l произвольный простой,

кусочно-гладкий контур, лежащий в области, необходимо и достаточно, чтобы в каждой точке этой области выполнялось условие @P@y = @Q@x .

Доказательство Напишем формулу Грина для произвольного просто-

го, кусочно-гладкого контура, лежащего в области

 

l

P (x; y) dx + Q (x; y) dy =

@x @y dxdy;

Z

 

ZZD

 

@P

 

 

 

 

@Q

 

 

где D - область, ограниченная контуром l:

Тогда, достаточность следует из этой формулы очевидно.

278

Докажем необходимость, то есть докажем, что, если интеграл по про-

извольному простому, кусочно-гладкому контуру равен нулю, то в каждой точке области выполняется условие @P@y = @Q@x . Доказательство проведем

от противного. Предположим, что найдется точка в области , в которой

@P@y 6= @Q@x , например предположим , что в этой точке

@Q@x @P@y > 0. Тогда

в силу непрерывности производных

@Q

@P

 

 

@x è

@y , найдется такая окрестность

этой точки, в каждой точке которой @Q@x @P@y

> > 0 и, если взять простой

контур, лежащей в этой окрестности, то

 

 

 

l P (x; y) dx + Q (x; y) dy =

@x @y dxdy > SD > 0;

Z

ZZD

 

@P

 

 

 

@Q

 

 

что противоречит условию.

Пример 252 Вычислить криволинейный интеграл второго рода, применяя формулу Грина

Z

xy2dy x2ydx; l : x2 + y2 = a2:

l

Решение По формуле Грина

ZZZ

xy2dy x2ydx = y2 + x2 dxdy;

l

D

где D - область , ограниченная окружностью x2 + y2 = a2:

Для вычисления двойного интеграла перейдем к полярным координатам:

x = r cos '; y = r sin ':

Тогда

 

 

 

 

 

2

a

 

 

 

a4

 

ZZ

x2

+ y2

 

dxdy = ZZ

r3drd' = Z0

d' Z0

r3dr =

 

:

 

 

 

2

 

D

 

 

D0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Пример 253 Найти площадь, ограниченную кривой x =

3at

; y =

3at

; 0

 

1+t3

1+t3

t < +1:

Решение По формуле (91) имеем

S =

1

 

xdy ydx =

1

 

+1

 

 

 

3at

 

 

 

3at2

 

 

 

 

3at2

3at

=

9a2

0

+1

 

t2 + t5 dt

9a2 +1

2

 

2

Z

0

 

 

 

1 + t3 d

1 + t3

 

1 + t3 d

1 + t3

2

Z

 

(1 + t3)3 =

2

Z

0

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9a2

 

+1 d t3 + 1

 

 

 

3a2

1

 

 

 

+1

 

 

 

3a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

Z0

(t3 + 1)2

=

 

 

 

 

 

0

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

3

2

 

t3 + 1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 254 Выполнено ли условие следствия 2 для функций

P =

 

 

x

 

 

 

 

 

 

x2+y2

 

 

 

 

 

 

è Q = x2+y2 . Вычислите интеграл

Rl

2

 

2

 

 

по окружности x2 + y2 = 1

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

ydx

xdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и прокомментируйте результат.

x +y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

(1 +

279

Следовательно, условие

Решение Проверим, будет ли выполнено условие @P@y = @Q@x :

@P

=

@

 

 

y

 

=

 

x2 y2

;

@y

@y

x2 + y2

 

(x2 + y2)2

 

 

 

@Q

=

@

 

 

x

 

=

 

x2 y2

:

@x

@x

x2 + y2

 

 

(x2 + y2)2

 

 

 

 

@P@y = @Q@x выполняется во всех точках плоскости, кроме начала координат, где функции P (x; y) и Q (x; y) не определены.

R ydx xdy

Теперь вычислим интеграл l x2+y2 параметрическому заданию окружности

ydx xdy

2

dt = 2 :

Тогда Rl x2+y2

= R0

.Для его вычисления перейдем к x = cos t; y = sin t:

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

гладкий контур, ограничивающий

R

 

D

 

S

Пример 255 Вычислить интеграл

l (x cos + y cos ) ds, где l простой

 

область

 

, площадь которой равна

 

и и - углы между вектором внешней нормали к кривой l и осями координат. Направление обхода кривой - положительно.

Решение Если направление обхода кривой положительно, то орт касательного вектора будет иметь координаты ~ = ( sin ; cos ). Тогда, пе-

реходя от интеграла первого рода к интегралу второго рода и используя

R R

формулу (90) получим l (x cos + y cos ) ds = l xdy ydx = 2S:

90Независимость криволинейного интеграла от пути интегрирования

Теорема 88 Пусть в односвязной области заданы две непрерывные и непрерывно дифференцируемые функции P (x; y) и Q (x; y). Тогда следующие три утверждения равносильны

R

1.l P (x; y) dx + Q (x; y) dy = 0, где l - произвольный простой, кусочно-

гладкий контур, лежащий внутри области .

R

2. P (x; y) dx + Q (x; y) dy по любой простой, кусочно-гладкой кривой, лежащей в области и соединяющей точки А и В, не завит от этой кривой. (Зависит только от точек А и В.)

3. Существует функция U = U (x; y) такая, что P (x; y) = @U@x ; Q (x; y) =

@U

P (x; y) dx + Q (x; y) dy является полным

@y . (При этом выражение

дифференциалом этой функции.)

Доказательство Докажем сначала, что из 1) следует 2). Пусть две кусочно-гладкие кривые l1 è l2 соединяют точки и , лежат в области и не пересекаются. Тогда их объединение образует простой кусочно-гладкий контур l, также лежащий внутри , но ведущий сначала из в по первой кривой,

R

а затем назад из в по второй. По условию теоремы l P (x; y) dx + Q (x; y) dy =

280