Добавил:
Я уверяю Вас, мне можно доверить огнестрельное оружие Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

А. Р. Лакерник. Высшая математика.Краткий курс

.pdf
Скачиваний:
5
Добавлен:
10.10.2023
Размер:
3.65 Mб
Скачать

VII. Кратные и криволинейные интегралы. Теория поля

Интеграл Эйлера–Пуассона

Теперь можно доказать приведенную в разд. 13.4 формулу для вы - числения интеграла Эйлера–Пуассона:

+∞

x2

p

 

 

 

 

ò e dx =

2 .

(19.13)

0

 

 

 

 

Сначала преобразуем исходный интеграл следующим образо м:

+∞

 

 

 

æ a

 

 

 

ö

1

 

a

 

 

 

a

2

ò

ex2dx = lim

ò

ex2dx = lim ç

ò

ex2dx

×

ò

ey2dy ÷

=

a→∞

 

a→∞ ç

 

 

÷

 

0

 

0

 

è

0

 

 

0

ø

 

 

æ

 

 

 

 

1

 

 

 

 

 

 

a

a

 

ö2

 

 

 

 

 

= ç lim→∞ òex2dx × òey2dy ÷ .

 

 

 

 

ç a

0

0

 

÷

 

 

 

 

 

è

 

 

ø

 

 

 

 

Произведение интегралов в круглых скобках равно

y

 

 

 

 

òò e

x2

e

y2

dxdy = òò e

(x2+ y2)

dxdy,

 

 

 

 

 

 

 

 

 

 

 

 

 

(D)

 

 

 

(D)

 

 

 

(D)

 

 

где областью (D) является заштрихованный

 

 

 

квадрат на рис. 113.

 

 

 

 

 

 

 

 

0

a

 

 

Действительно,прирасстановкепределовпо

2

x

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

 

a

 

Ðèñ. 113

 

 

ралов:

 

 

 

 

 

 

 

 

 

 

 

a

a

 

 

 

a

a

 

 

òò e(x2+ y2)dxdy = òdx òex2ey2dy = òex2dx òey2dy.

 

(D)

 

 

0

0

 

 

 

0

0

 

a

В последнем интеграле внутренняя часть òey2dy не зависит от х ,

0

следовательно, ее можно вынести за знак внешнего интеграла. Тогда òò e(x2+ y2)dxdy = aòey2dy × aòex2dx , что мы и утверждали.

(D)

0

0

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

 

 

 

æ

 

 

1

 

+∞

 

 

 

 

ö2

 

ò

ex2dx

=

ç

lim

òò

e(x2+ y2)dxdy ÷

(19.14)

 

ça→+∞

÷ .

0

 

 

è

 

(D)

ø

 

19. Кратные интегралы

Учитывая положительность подынтегральной функции и обо зна- чая изображенные на рис. 113 четверти кругов радиусов а и a 2 через (D1) è (D2) соответственно, имеем

òò e(x2+ y2)dxdy £ òò e(x2+ y2)dxdy £ òò e(x2+ y2)dxdy.

(19.15)

(D1)

(D)

(D2)

 

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

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2

 

2)

2

a

 

2

 

 

p

æ

 

1 ö

a

 

2

 

 

 

p

 

r 2

 

a

 

òò

e

+ y

dxdy = òdjòe

r

rdr =

-

òe

r

d(-r

2

) = -

e

 

=

 

 

 

 

 

 

ç

÷

 

 

 

 

 

 

0

(D )

 

 

 

 

0

0

 

 

 

 

2

è

 

2 ø

0

 

 

 

 

 

4

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= p

(1- ea2 );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

òò e(x2+ y2)dxdy = p4 (1-e2a2 ).

(D2)

Для нахождения предела среднего интеграла неравенства (19.15) применим теорему 2.12:

 

 

 

 

 

 

 

 

теоремаа

 

 

 

 

 

 

 

 

 

 

 

p

2

p

 

p

 

 

p

2.

 

 

 

2+ 2

 

 

p

 

 

2

2.12

a→+∞ òò

 

 

.

a→+∞ 4

(1- e a ) =

4

;

a→+∞ 4

(1- e 2a

) =

4

Þ

e (x

 

y

)dxdy =

4

lim

 

lim

 

lim

 

 

 

 

 

 

 

 

 

 

 

 

 

(D)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

2

 

p

 

 

 

p

 

 

Отсюда по формуле (19.14) получаем: ò e

dx =

=

 

 

 

 

 

x

4

 

 

.

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

19.7. Замена переменных в тройном интеграле

Рассмотрим òòò f (x,y,z)dxdydz , где область (T) ограничена гладкой

(T )

или кусочно-гладкой поверхностью, а функция u = f (x,y,z) непрерывна в этой области. Пусть x = x(u,v,w), y = y(u,v,w), z = z(u,v,w) , где эти функции осуществляют взаимно однозначное и в обе стор оны непрерывно дифференцируемое соответствие между точками о бластиT() пространства xyz и точками некоторой области (T ¢) пространства uvw.

398

399

VII. Кратные и криволинейные интегралы. Теория поля

При этом границе одной из этих областей соответствует гра ница другой, и наоборот. Аналогично результатам разд. 19.5, получим, чт о

òòò f (x,y,z)dxdydz =

(T )

= òòò f (x(u,v,w),y(u,v,w),z(u,v,w))

 

I(u,v,w)

 

dudvdw,

(19.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(T )

 

 

 

 

 

 

 

 

 

 

 

 

 

где якобиан

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

u

 

u

 

 

u

 

 

 

I(u,v,w) =

x

 

y

 

z

 

 

(19.17)

v

 

v

 

v

 

 

 

x

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

w

 

 

 

w

 

 

 

 

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

V = I(u,v,w)Du Dv Dw ).

Рассмотрим два частных случая формул (19.16),(19.17).

Тройной интеграл в цилиндрических координатах

Цилиндрические координаты представляют собой соединени е полярных координат на плоскости 0xy с обычной декартовой аппликатой z (рис. 114).

На рис. 114 точка M имеет координаты (r,j,z) , где r и ϕ – полярные координаты точки Р – проекции M на плоскость 0xy, z – обычная аппликата точки M. Формулы перехода к цилиндрическим координа-

там имеют вид: x = r cosϕ,

y = r sinϕ,

z = z . По формуле (19.17) нахо-

дим якобиан этих функций:

 

 

 

 

 

 

 

cosj

sin j

0

 

 

cosj

sinj

 

 

 

 

 

 

 

 

I(r,j,z) =

-r sin j

r cos j

0

 

= 1×

=

 

0

0

1

 

 

-r sinj

r cosj

 

 

 

 

 

 

 

 

= r cos2 ϕ + r sin2 ϕ = r.

 

(19.18)

19. Кратные интегралы

z

M (r, j, z)

z

0

y j r

P

x

Ðèñ. 114

К цилиндрическим координатам стоит переходить, если урав нение границы области или, может быть, подынтегральная функция содержит выражение x2 + y2 = r2.

Пример. Вычислить òòòzdxdydz , где область (T) ограничена поверхно-

 

 

(T )

 

стями y = x, y =

x

, x2 + y2 = x, z = 0,

z = x2 + y2 .

3

 

 

 

Ð å ø å í è å

Первые три поверхности – цилиндрические, параллельные ос иz0; они проектируются на плоскость 0xy в границу области, изображенной на рис. 115. Рссмотрим уравнения этих границ:

x

2

+ y

2

æ

1

ö2

2

=

1

;

 

 

= x Û ç x -

2

÷

+ y

 

4

 

 

 

 

è

ø

 

 

 

 

y= x: k =tgj =1Þj = p

;

 

1

1

1

4

 

y= x : k

2

=tgj =

1 Þj = p.

3

 

2

3

2

6

 

 

 

 

Перейдем к цилиндрическим координатам. Уравнение окружности примет вид r 2 = r cosj Ю r = cosj, уравнение верхней «крышки» преобразуется к z = r2 = r , и интеграл примет вид

y

0

 

1

 

1 x

 

 

2

 

 

Ðèñ. 115

 

400

401

VII. Кратные и криволинейные интегралы. Теория поля

 

π

 

 

 

π

 

 

 

π

 

4

cosϕ

r

1

4

cosϕ

 

1

4

òòòzdxdydz = òdj

ò

rdr ò zdz =

ò dj

ò

r3dr =

òcos4 j dj =

2

8

(T )

π

0

0

 

π

0

 

 

π

 

6

 

 

 

6

 

 

 

6

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

1

 

4

 

 

 

 

 

 

2

 

 

 

1

4 æ

 

 

 

1+cos4j ö

 

=

 

 

ò

(1+cos2j)

 

dj=

 

 

 

ç1+2cos2j

+

 

 

÷dj=

32

 

32

2

 

 

 

 

 

 

 

 

 

 

òè

 

 

 

ø

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

p

 

 

 

p

 

 

 

 

p ö

 

2p+16 – 9 3

 

 

=

 

1 ç 3

j

4

+sin 2j

4

+

1

sin4j

4 ÷

=

.

 

32

ç

2

p

p

8

p ÷

 

512

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

6

 

 

 

6

 

 

 

 

6 ø

 

 

 

 

 

 

Тройной интеграл в сферических координатах

В сферических координатах положение точки M в пространстве определяется тремя числами, ρ, ϕ, θ (рис. 116): r – это расстояние от точки М до начала координат (r ³ 0 ); j – это тот же угол, что в цилиндрических координатах, т.е. угол от положительного направ ления оси 0x до вектора OP ( ϕ (−π,π]); q – это угол между положительным направлением оси 0z и вектором OM ( θ [0,π]). Угол j показывает, в какую сторону вектор OM отклоняется от оси 0z, а угол q показывает, как велико такое отклонение.

Формулы перехода к сферическим координатам имеют вид:

x = OP cosj = rsin qcosj; y = OPsinj = rsinqsinj; z = rcosq.

z

M (r,j, q)

q q

 

r

z

 

 

0

y

 

 

 

j

 

y

 

 

 

 

P

x

Ðèñ. 116

 

 

 

19. Кратные интегралы

По формуле (19.17) определяем якобиан этих функций:

 

¶x

¶y

¶z

 

 

 

 

 

 

 

 

 

 

 

 

¶r

¶r

¶r

 

sin qcosj

sin qsinj

cosq

 

 

 

 

 

¶x

¶y

¶z

 

 

 

I(r,j,q) =

=

-rsin qsinj

rsin qcosj

0

.

¶j

¶j

¶j

 

 

rcosqcosj

rcosqsinj

-rsinq

 

 

¶x

¶y

¶z

 

 

 

 

 

 

 

 

 

¶q

¶q

¶q

 

 

 

 

 

Раскладывая этот определитель по элементам второй строк и, находим:

I(r,j,q) = rsinqsinj

 

sinqsinj cosq

 

+rsinqcosj

 

sinqcosj cosq

 

=

 

 

 

 

 

 

rcosqsinj -rsinq

 

 

 

rcosqcosj -rsinq

 

 

= rsin qsinj(-rsin2 qsin j - rcos2 qsin j) + rsinqcosj(-rsin2 qcosj -

-rcos2 qcosj) = -r2 sin qsin2 j ×1- r2 sinqcos2 j ×1 = -r2 sin q ×1.

Значит,

 

I(r, j, q)

 

= r2 sin q.

(19.19)

 

 

К сферическим координатам стоит переходить, когда уравне ние границы области или подынтегральная функция содержит вы ражение x2 + y2 + z2 , которое в этих координатах преобразуется следующим образом:

x2 + y2 + z2 = r2 sin2 qcos2 j + r2 sin2 qsin2 j + r2 cos2 q =

= ρ2 sin2 θ + ρ2 cos2 θ = ρ2.

Пример. Найти объем тела, ограниченного

конусом z2 = x2 + y2 (z ³ 0) и сферой x2 + y2 + z2 = 1 z (z ³ 0) (ðèñ. 117).

Ð å ø å í è å

 

 

 

Сечение конуса плоскостью 0yz (х = 0) име-

 

 

åò âèä z2 = y2 Ы z = ±y , значит, максимальное

0

y

p

 

 

 

 

 

 

значение угла q равно 4

(так как угловой коэф-

 

 

фициент прямых z = ±y

равен k = ±tg p

= ±1 ) è

 

 

объем тела

4

 

x

 

 

 

Ðèñ. 117

 

402

403

VII. Кратные и криволинейные интегралы. Теория поля

 

 

 

π

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

4

 

π

 

1

 

 

1

2p æ

 

 

ö

 

 

 

 

 

 

4

r3

 

2

 

 

 

 

 

 

 

 

 

 

V =

òòò

dxdydz =

ò

dq

ò

dj

ò

r2 sin qdr = 2p(-cos q)

 

 

=

 

ç1

-

 

÷.

 

3

3

2

 

 

 

 

 

 

 

ç

 

÷

 

(T )

 

0

 

−π

 

0

 

0

 

0

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

20. КРИВОЛИНЕЙНЫЕ ИНТЕГРАЛЫ

20.1. Криволинейный интеграл первого рода

Определение 20.1. Пусть в пространстве задана кривая (AB) и пусть на ней задана некоторая функция f (M ) = f (x,y,z), где M (x,y,z) (AB ). Разобьем кривую (AB ) на части точками Ai (xi ,yi ,zi ). На каждой части

 

~ ~ ~

 

– äëè-

(Ai Ai+1) возьмем произвольную точку Mi ( xi , yi , zi ). Пусть si

íà äóãè (Ai

Ai+1). Составим интегральную сумму σ = lim

n1

si. Îáî-

å f (Mi )

 

λ→0

i=0

 

 

 

 

значим λ максимальную поi длину дуги (Ai

Ai+1): λ = max si . Еслисуще-

 

i

ствует предел интегральных сумм при λ → 0, который не зависит от разбиения(AB )начасти(AiAi+1)иотвыбораточек Mi (Ai Ai+1),тоэтотпредел называется криволинейным интегралом первого рода от функции f (M ) =

= f (x,y,z) по кривой (AB ) и обозначается

ò

f (M )ds =

ò

f (x,y,z)ds.

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

 

(AB)

(AB)

 

 

 

 

 

 

 

 

 

n1

 

 

ò f (M )ds =

ò f (x,y,z)ds =

λ→lim0

å f (Mi )

si ,

(20.1)

(AB)

(AB)

 

i=0

 

 

еслиэтотпределсуществуетинезависитотвыбораточекAi èMi (рис.118). Определения криволинейного и обычного интегралов анало гичны, поэтому совпадают и одинаково доказываются их свойства. П риведем

некоторые из этих свойств:

A = A0

A1

B = An

 

 

A2

Ai

Ai = 1

Mi

Ðèñ. 118

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

1)

ò

αf (M )ds = α ò f (M )ds , если интеграл справа существует;

 

(AB)

 

(AB)

 

2)

ò

[ f (M ) + g(M )]ds = ò

f (M )ds + ò g(M )ds, если интегралы

 

(AB)

 

 

(AB)

(AB)

справа существуют;

 

 

 

3)

ò

f (M )ds =

ò

f (M )ds + ò f (M )ds , åñëè (AB) = (AC) (CB)

 

(AB)

 

(AC)

 

(CB)

и все эти три интеграла существуют;

4)

ò

f (M )ds =

ò

f (M )ds при условии существования этих ин-

 

(AB)

 

(BA)

 

 

тегралов (это свойство следует из определения 20.1).

Вычисление криволинейного интеграла первого рода

Пусть кривая (AB) задана параметрическими уравнениями x = x (t), y = y (t), z = z (t), t [α,β], α < β, ãäå x = x(t), y = y(t), z = z(t) непрерывно

дифференцируемы на отрезке

α,β

]

, а функция

f (x,y,z)

непрерывна

 

[

 

 

вдоль (AB ).

 

 

 

 

 

Можно доказать, что при этих условиях интеграл ò

f (M )ds ñó-

(AB)

ществует. Не вдаваясь в это доказательство, получим формулу для вычисления интеграла. Пусть точки Ai è Ai+1 соответствуют значениям параметра ti è ti+1, тогда запишем si по аналогии со случаем двух переменных х и y (см. формулу (10.14)):

si = tiò+1 [x(t)]2 + [y(t)]2 + [z(t)]2dt.

ti

Применяякпоследнемуинтегралутеорему10.1осреднем,полу чаем

 

 

 

s i =

[x

( ti )]

+ [y (ti

)]

+ [z ( ti

)]

 

t i ,

 

 

 

 

 

 

 

 

~

2

~

2

~

 

2

 

 

 

ãäå

~

],

t i

=t i+1 –t i

. Тогда

 

 

 

 

 

 

 

ti [t i , t i+1

 

 

 

 

 

 

 

 

ò

 

 

0

n–1

 

 

~

 

~

 

 

~

 

 

 

 

 

 

λ→lim åf (M i

 

2

 

 

2

+[z (ti )]

2

 

 

 

f (x,y,z )ds=

) [x (ti )]

+[y (ti )]

 

t i .

( AB )

 

 

i=0

 

 

 

 

 

 

 

 

 

 

404

405

VII. Кратные и криволинейные интегралы. Теория поля

Так как интеграл ò f (M )ds существует, то для его вычисления

(AB)

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

соответствует значению параметра ~ , тогда эта сумма является интег- ti

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

b

ò f ( x,y,z) ds= ò f ( x(t),y(t),z(t)) [ x¢(t)]2 +[ y¢(t)] 2 +[ z¢(t)]2 dt. (20.2)

( AB) a

(в этой формуле интеграл справа существует в силу непреры вности подынтегральной функции).

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

Определение 20.2. Пусть в пространстве задана кривая (AB ) и пусть на ней задана некоторая функция P (M ) = P (x,y,z), где M (x,y,z) (AB ). Разобьем кривую (AB ) на части точками Ai (xi ,yi ,zi ). На каждой части

~ ~ ~

=

(Ai Ai+1) возьмем произвольную точку Mi ( xi , yi , zi ). Пусть xi

= xi+1 – xi , ò.å. xi – это проекция дуги (Ai Ai+1) на ось 0x (длина этой про-

 

n1

 

екции со знаком). Составим интегральную сумму σ = lim

åP(Mi )

xi .

λ→0

i=0

 

Обозначим через λ максимальную по i длину дуги (Ai Ai+1): λ = max

si .

 

i

 

Если существует предел наших интегральных сумм при λ → 0 , который

не зависит от разбиения (AB) на части (Ai Ai+1) и от выбора точек

Mi (Ai Ai+1), то этот предел называется криволинейным интегралом вто-

рого рода от функции P (M ) =P (x,y,z) по кривой (AB ) по переменной х и обозначается ò P(x,y,z)dx .

(AB)

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

 

 

n1

 

n1

 

ò

P(x,y,z)dx = λ→lim0

åP(Mi )

xi = λ→lim0

åP (xi ,yi ,zi ) xi ,

(20.3)

(AB)

 

i=0

 

i=0

 

если этот предел существует и не зависит от выбора точек Ai è Mi (ñì. ðèñ. 118).

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

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

Аналогично определяются

 

 

n1

 

 

 

n

 

ò

Q(x,y,z)dy = lim å Q(Mi ) y i

,

ò

R(x,y,z)dz = lim åR(Mi )

z i

 

λ→0 i=0

{

 

λ→0 i=0

{

(AB)

 

y i+1 –yi

 

(AB)

 

zi+1 –z i

è ò

Pdx +Qdy + Rdz как сумма таких интегралов:

 

(AB)

 

 

 

 

 

 

 

 

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

 

 

 

( AB)

 

 

 

 

 

 

 

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

 

 

 

( AB)

( AB)

 

 

( AB)

 

Как и в предыдущем разделе свойства криволинейных интегр алов второго рода аналогичны свойствам обычных интегралов. А и менно:

1)

ò αP(x,y,z)dx = α ò P(x,y,z)dx, если интеграл справа суще-

ствует;

(AB)

(AB)

 

 

 

 

 

 

2) ò [P1(x,y,z) + P2(x,y,z)]dx = ò P1(x,y,z)dx + ò P2(x,y,z)dx,

 

(AB)

 

(AB)

(AB)

если интегралы справа существуют;

 

3)

ò P(x,y,z)dx =

ò P(x,y,z)dx + ò

P(x,y,z)dx, åñëè (AB ) =

 

(AB)

(AC)

(CB)

 

(AC ) (CB ) (рис. 119) и все эти три интеграла существуют;

4) из определения 20.2 также следует, что в отличие от интеграла

первого рода ò P(x,y,z)dx = − ò

P(x,y,z )dz (при составлении интег-

 

(AB)

(BA)

 

B

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

C

 

 

 

 

Ai è Ai+1 меняются местами, и xi

– xi+1 =

 

= –(xi+1 – xi)).

 

 

A

 

 

Ðèñ. 119

406

407

VII. Кратные и криволинейные интегралы. Теория поля

С учетом этого свойство 3 будет верным независимо от взаим ного расположения точек A, B и C .

Естественно, все перечисленные свойства справедливы и дл я интегралов по переменным y и z, и для суммы таких интегралов.

D

 

Если контур (L) замкнут, то криволи-

 

нейный интеграл вдоль (L) в заданном на-

 

B

 

правлении определяется так: пусть A (L),

 

 

 

 

B (L) – произвольные точки, тогда для

 

 

интеграла по любой переменной и для

 

 

суммы этих интегралов по определению

A

ò

=

ò

+

ò

(ðèñ. 120).

 

C

(L)

(ACB) (BDA)

 

Ðèñ. 120

 

При другом выборе точек значение ин-

 

теграла не изменится:

ò

+ ò

= ò

+ ò

+ ò

+ ò ;

(ACB)

(BDA)

(AC)

(CB )

(BD) (DA)

ò + ò = ò

+ ò

+ ò

 

+ ò

Þ

ò + ò = ò + ò .

(CBD) (DAC) (CB)

(BD)

(DA) (AC)

(ACB) (BDA) (CBD) (DAC)

Если (L) – замкнутый контур, то под интегралом ò Pdx +Qdy + Rdz

 

 

 

 

 

 

(L)

при отсутствии указания на направление обхода контура по нимают ин-

 

теграл, взятый в положительном направлении

 

(так, чтобы при движении вдоль контура бли-

 

жайшая часть области, ограниченной этим кон-

 

туром, оставалась бы слева, что в простых слу-

 

чаях равносильно обходу контура против часо-

Ðèñ. 121

вой стрелки – рис. 121).

В случае замкнутого контура иногда также пишут

òPdx+Qdy+ Rdz.

( L )

Вычисление криволинейного интеграла второго рода

Теорема 20.1. Пусть кривая (AB) задана параметрическими уравнениями x = x(t), y = y(t), z = z(t) , t О[a,b] и при изменении парамет-

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

ðà t îò α до b кривая описывается в направлении от A к B; при t = α получаем точку A, при t = b – точку B. Пусть также функции P(x,y,z), Q(x,y,z), R(x,y,z) непрерывны вдоль кривой (AB), а функции х(t), у(t), z(t) – непрерывно дифференцируемы на отрезке [a,b] . Тогда криволинейный интеграл существует и

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

(AB)

b

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

a

(20.4)

b

 

= ò[ P( x(t),y(t),z(t)) x¢(t)+Q( x(t),y(t),z(t)) y¢(t)+ R( x(t),y(t),z(t)) z¢(t)] dt

a

 

(т.е. для вычисления надо и в подынтегральных функциях, и под знаками дифференциалов заменить x, y, z на их выражения через параметр t).

¡ Докажем одну часть этого равенства (остальные части дока зываются аналогично):

 

b

 

òP( x,y,z) dx= òP( x(t),y(t),z(t)) x¢(t)dt.

(20.5)

( AB)

a

 

В соответствии с определением 20.2 левая часть этого равенст ва

 

 

 

n1

 

 

 

 

 

n1

 

 

 

ò

P(x,y,z)dx = limσ = lim

å

P(M

 

)

x

= lim

å

P (x ,y ,z

 

) x .

λ→0

λ→0

 

i

 

i

λ→0

i i

i

i

(AB)

 

 

i=0

 

 

 

 

 

i=0

 

 

 

Докажем, что этот предел равен правой части равенства, т.е. при достаточно малых λ интегральные суммы σ будут отличаться от интегралавправойчасти формулы (20.5) меньше,чемналюбоенапередзаданноечисло.

Зададим произвольное число ε > 0 . Пусть в интегральной сумме σ

точки Ai è Mi соответствуют значениям параметра ti è

~

. Тогда

ti

n–1

 

 

 

 

 

~

~

~

))Dxi,

 

 

s=åP(x(ti

), y(ti

), z(ti

 

 

i=0

ãäå xi = xi+1 xi = x(ti+1) x(ti ).

408

409

VII. Кратные и криволинейные интегралы. Теория поля

 

 

ti +1

 

 

 

 

 

 

Используя равенство

ò

 

ti+1

 

 

 

 

 

 

 

x (t)dt = x(t)

 

ti = x(ti+1) x(ti

) = xi , ïåðå-

 

 

 

 

 

ti

 

 

 

 

 

 

пишем предыдущую формулу в виде

 

 

 

 

 

n–1

~ ~ ~

ti +1

n–1 ti +1

~ ~

 

~

s= åP (x (ti ), y (ti ), z (ti )) ò x¢( t) dt= å

 

òP (x(ti ), y (ti ), z

(ti ))x¢(t) dt.

i=0

 

t

i=0

 

t

 

 

 

 

 

i

 

 

i

 

 

 

Преобразуем интеграл в правой части формулы (20.5) следующим образом:

β

n1

ti+1

òP(x(t),y(t),z(t))x¢(t)dt = å

ò

P(x(t ),y(t ),z(t ))x ¢(t )dt

α

i=0

t

 

 

 

i

 

и рассмотрим разность двух полученных выражений:

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

σ − òP(x(t),y(t),z(t))x (t)dt =

 

 

 

 

 

α

 

 

n–1 ti +1

~

~

 

~

n–1 ti +1

=å

 

 

 

 

òP(x(ti ),

y(ti ),

z(ti ))x¢( t) dt –å òP( x( t), y( t), z( t) ) x¢( t) dt =

i=0

t

i

 

 

 

 

i=0 t

 

 

 

 

 

 

i

 

 

n–1 ti +1

~

 

 

~ ~

 

=å ò[P(x(ti ),

y(ti ), z(ti )) – P( x( t), y( t), z( t) )] x¢( t) dt.

i=0 ti

Отсюда из свойств абсолютных величин и свойств интеграло в получаем

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s - òP(x(t),y(t),z(t))x¢(t)dt

£

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n–1

ti +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£å

ò

 

~

~

~

 

 

 

 

x¢(t)

 

dt

.

 

(20.6)

 

 

 

 

 

 

P(x(ti

), y(ti

), z(ti )) – P(x(t), y(t), z(t))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=0

t

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[

 

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как по условию функция x(t) непрерывна на отрезке α,β

 

, òî îíà

ограничена на этом отрезке:

 

 

K , где K – некоторая постоянная.

 

 

 

x (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[

]

 

Сложная функция P(x(t), y(t), z(t)) непрерывна на отрезке α,β

 

,

тогдапотеоремеКантора3.9онаравномернонепрерывнанаэтомотрезке. Исходя из этого можно выбрать l столь малым, что на каждом отрезке разбиения [ti , ti+1] модуль разности значений функции в двух про-

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

извольных точках не будет превосходить любое наперед зад анное по- e

ложительное число, за которое в нашем случае возьмем K b - a . Тогда по формуле (20.6) получим

 

β

 

 

n1

 

e

 

 

 

 

 

e

n1

 

 

 

 

 

 

 

 

 

s - òP(x(t), y(t),z(t))x¢(t)dt

 

£

å

 

K

 

ti+1 - ti

 

=

å

 

ti+1 - ti

 

.

 

 

 

 

 

 

 

 

 

K

b - a

b - a

 

 

α

 

 

i=0

 

 

 

 

 

i=0

 

 

 

 

Так как при движении вдоль кривой (АВ) параметр t либо растет, либо уменьшается, то все слагаемые последней суммы имеют один и тот же знак:

n1

 

 

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

å

 

ti+1 - ti

 

=

 

å(ti+1 - ti )

=

 

 

(t1 – t0)+(t

2 – t1)+...+ (tn – tn –1)

 

 

 

 

 

 

 

 

=

tn – t0

=

 

 

i=0

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

b - a

 

Þ

 

s - òP(x(t),y(t),z(t))x¢(t)dt

£

 

 

b - a

 

= e.x

 

 

 

 

 

 

 

 

 

 

 

 

 

b - a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1. Вычислить ò xdx + ydy + (x + y -1)dz, где L – отрезок прямой

от точки

 

 

 

 

 

 

 

 

 

 

 

 

 

(L)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А (1, 1, 1) до точки В (2, 3, 4).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения прямой имеют вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x -1

=

y -1

=

z -1

Þ

x -1

=

y -1

=

z -1

= t Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 -1 3-1 4 -1

1

2

 

 

 

3

 

 

 

 

 

 

 

 

 

x = t +1, y = 2t +1, z = 3t +1 – параметрические уравнения этой прямой.

Ïðè

 

 

t = 0 получаем x = 1, y = 1, z = 1 , т.е. точку А; при t = 1 получаем

x = 2, y = 3, z = 4 , т.е. точку В. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

1

ò xdx + ydy + (x + y -1)dz = ò(t +1)dt + (2t +1)d(2t +1) +(t +1 +2t +1 -1)d(3t +1) =

(L)

0

= ò1

(t + 1)dt + (2t +1)2dt +(3t +1)3dt = ò1

(t +1 +4t +2 +9t +3)dt =ò1(14t +6)dt =

0

 

 

 

 

 

0

 

0

 

=14

t2

 

 

1

+6t

 

1

= 7 + 6 =13.

 

 

 

 

 

 

 

 

2

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

410

411

VII. Кратные и криволинейные интегралы. Теория поля

y

1

y = x 2

y = x

0

1 x

 

Ðèñ. 122

 

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

ò

ydx + xdy , ãäå

 

 

(L)

 

 

 

 

 

 

L – путь от точки (0, 0) до точки (1, 1) вдоль:

à) y = x ; á) y = x2 (ðèñ. 122).

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

à) y = x: берем за параметр x [0,1], тогда

ò

ydx + xdy = ò1 xdx + xdx =2ò1 xdx = 2

x2

 

 

1

= 1;

 

 

(L)

0

0

2

 

0

 

 

 

á) y = x2 : опять берем за параметр x [0,1] , тогда

 

1

1

1

x3

 

 

1

 

 

ò

ydx + xdy = òx2dx + xdx2 = ò

(x2 + 2x2 )dx =3òx2dx = 3

 

 

=1.

3

(L)

0

0

0

 

 

0

 

 

В обоих случаях получили один и тот же ответ. Возникает воп рос: не является ли этот факт правилом независимости криволиней ного интеграла (второго рода) от формы пути интегрирования? Приведен ный ниже пример показывает, что, вообще говоря, такого правила быть не может.

Пример 3. Для тех же путей, что в примере 2, вычислить ò 2ydx + xdy.

à) y = x : ò 2ydx

(L)

á) y = x2 :

1

ò 2ydx + xdy = ò

(L)

0

Ð å ø å í è å

 

 

 

 

 

 

(L)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

x2

 

 

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ xdy = ò2xdx + xdx =3òxdx = 3

 

 

=

;

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

0

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

x3

 

 

 

4.

 

 

 

 

 

 

 

2x2dx + xdx2 = ò

(2x2 + 2x2 )dx =4 òx2dx = 4

 

 

 

=

 

 

 

0

 

 

0

 

3

 

 

0

 

3

 

 

 

 

 

 

Однакоприопределенныхусловияхуказаннаянезависимос тькриволинейногоинтегралаотформыпутиинтегрированиявсежеи меетместо.

20.3. Условия независимости криволинейного интеграла второго рода от формы пути интегрирования

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

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

целиком принадлежащей множеству; открытое связное множе ство на-

зывается областью.

Далее будем рассматривать криволинейный интеграл

ò Pdx + Qdy + Rdz, в котором функции P(x,y,z), Q(x, y,z), R(x,y,z)

(AB)

непрерывны в некоторой области (Т) трехмерного пространства. Теорема 20.2. Для того чтобы ò Pdx +Qdy + Rdz при любых кри-

(AB)

вых (АВ) (Т) не зависел от формы пути интегрирования (АВ), а зависел только от положения начальной и конечной точек А и В, необходи-

мо и достаточно, чтобы

ò

Pdx +Qdy + Rdz = 0 для любого замкнутого

контура (L) (T).

 

(AB)

 

 

 

 

 

 

¡ Необходимость. Пусть наш ин-

D

теграл не зависит от формы пути ин-

B

тегрирования. Тогда для произволь-

 

ного замкнутого контура (L) (T),

 

изображенного на рис. 123, имеем

C

 

 

 

A

 

ò Pdx +Qdy + Rdz =

ò

Pdx +Qdy + Rdz +

Ðèñ. 123

 

(L)

(ACB)

 

 

+ ò Pdx + Qdy + Rdz =

ò

Pdx + Qdy + Rdz

ò Pdx + Qdy + Rdz = 0,

(BDA)

 

(ACB)

(ADB)

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

выполнении условий теоремы

ò Pdx + Qdy + Rdz =

ò Pdx + Qdy + Rdz.

(ACB)

(ADB)

Для этого докажем, что разность левой и правой частей этог о равенства есть 0:

ò Pdx +Qdy + Rdz ò Pdx +Qdy + Rdz = ò Pdx + Qdy + Rdz +

(ACB) (ADB) (ACB)

+ ò

Pdx + Qdy + Rdz = ò Pdx +Qdy + Rdz = 0,

(BDA)

(ACBDA)

как интеграл по замкнутому контуру. x

412

413

(AB)
ò

VII. Кратные и криволинейные интегралы. Теория поля

Теорема 20.3. Для того чтобы Pdx +Qdy + Rdz не зависел от

формы пути интегрирования для любой кривой (АВ )М (Т), необходимо и достаточно, чтобы выражение Pdx +Qdy + Rdz являлось в области (Т) дифференциалом отнекоторойфункциитрехпеременных u(x,y,z), т.е. чтобы существовала такая функция u(x,y,z) , что du(x,y,z) =

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

¡Необходимость. Пусть криволинейный интеграл

ò Pdx +Qdy + Rdz

(AB)

не зависит от формы пути интегрирования. Тогда введем фун кцию

(x,y,z )

u(x,y,z )= ò Pdx +Qdy + Rdz ,

(x0,y0,z0)

ãäå (x0,y0,z0) М (T ) – фиксированная точка и криволинейный интеграл берется по любому пути в (Т), соединяющему (x0,y0,z0) и (х, у, z). Тогда, считая, что путь в точку (x + Dx,y,z) проходит через точку (х, у, z), имеем

 

 

(xx,y,z)

Dxu = u(x + Dx,y,z) - u(x,y,z) =

ò

Pdx +Qdy + Rdz -

 

 

(x0,y0,z0)

(x,y,z)

(xx,y,z)

 

- ò Pdx +Qdy +Rdz =

 

ò

Pdx +Qdy + Rdz.

(x0,y0,z0)

 

(x,y,z )

 

Последний интеграл берется по любому пути. Возьмем за так ой путь отрезок, соединяющий точки (х, у, z) и (x + Dx,y,z) ; его уравнения: x = t, t О[x,x + Dx]; y = y, z = z, – постоянные, следовательно, dy = dz = 0. Отсюда

( x+ Dx,y,z)

x+ Dx

теорема 10.1

о среднем

 

 

ò

Pdx+Qdy+ Rdz=

ò

P( t,y,z) dt

= P( c,y,z) ( x+ Dx – x) ,

 

 

 

 

 

 

14243

 

 

 

( x,y,z)

 

x

 

 

 

 

 

Dx

 

где с – точка между х и x + Dx, следовательно

 

 

 

 

 

 

 

 

 

 

 

 

 

P(x,y,z ) -

 

 

 

 

 

 

 

 

 

 

 

P (x,y,z) –

 

 

Dxu

 

 

 

 

 

Dxu

 

 

непре-

 

 

 

 

¶u(x,y,z ) =

 

 

непрерывнавна

 

 

= P(c,y,z )Þ

lim

= lim P (c,y ,z )

=

P (x

,y ,z ).

 

 

 

 

Dx

 

¶x

 

x0

Dx

x0

 

 

 

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

Аналогично

¶u(x,y,z) =Q(x,y,z) è

¶u(x,y,z) = R(x,y,z). Следова-

тельно,

 

¶y

 

 

¶z

 

 

 

 

 

du =

¶u dx +

¶u dy +

¶u dz = Pdx +Qdy + Rdz.

 

 

¶x

¶y

¶z

 

Достаточность. Пусть Pdx +Qdy + Rdz = du и кривая (АВ) задана параметрическими уравнениями x = x(t),y = y(t),z = z(t),t О[a,b] . Òîã-

¶u

 

¶u

 

¶u

 

äà P = ¶x

, Q =

¶y

, R =

¶z

è

 

 

 

 

b

 

òPdx+Qdy+ Rdz = ò[ P( x(t),y(t),z(t)) x¢(t)+Q( x(t),y(t),z(t)) y¢(t)+

( AB)

 

 

 

a

 

 

 

β é

+R(x(t),y(t),z(t))z¢(t)]dt = ò ê

 

 

α ë

+

¶u (x(t),y(t),z(t))y¢(t) +

¶u

 

¶y

¶z

¶u (x(t),y(t),z(t))x¢(t) + ¶x

(x(t),y(t),z(t))z¢(t)ùdt.

ú

û

Согласно правилу дифференцирования сложной функции (см. разд. 12.4) последнее выражение можно представить в виде

β

ò é d u(x(t),y(t

ê

α ëdt

Èòàê, åñëè ¶u

¶x

 

ù

 

 

½t = βu (x (t ), y (t ), z (t ))½t = α =

 

 

 

),z(t))ú dt = u(x(t),y(t),z(t))

 

 

 

û

 

 

 

 

 

 

 

 

= u (B ) u (A ).

 

dx +

¶u dy +

¶u dz = du , òî

 

 

¶y

¶z

 

ò

Pdx +Qdy + Rdz = u(B) -u(A) = u

 

B

(20.7)

 

 

A .

(AB)

 

 

 

 

 

 

Правая часть формулы (20.7) не зависит от пути из А в В, что доказывает независимость интеграла от формы пути интегриров анияx.

Заодно мы получили формулу для вычислений интеграла от по л- ного дифференциала и для примера 2 в разд. 20.2:

414

415

VII. Кратные и криволинейные интегралы. Теория поля

(1,1)

(1,1)

 

(1,1)

 

 

 

 

ò ydx + xdy =

ò d(xy) = (xy)

 

= 1 для любого пути.

(0,0)

(0,0)

 

(0,0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 20.4. Пусть функции P(x,y), Q(x,y),

¶P(x,y)

,

¶Q(x,y)

 

¶x

 

 

 

 

 

¶y

непрерывны в некоторой области (D) на плоскости 0xy. Для того чтобы интеграл ò Pdx +Qdy не зависел от формы пути интегрирования для

(AB)

любой кривой (АВ)М (D), необходимо, а в предположении односвязности области и достаточно, чтобы при (x,y)ОD

¶P(x,y)

=

¶Q(x,y).

(20.8)

 

¶y

¶x

 

Эта теорема сразу следует из теорем 20.3 и 14.2, которая утвержда - ет, что (20.8) есть необходимое, а в предположении односвязнос ти области и достаточное условие того, чтобы выражение P (x,y)dx + Q (x,y)dy являлось дифференциалом некоторой функции двух перемен ных. Напомним здесь, что область (D ) на плоскости называется односвязной, если вместе с каждым замкнутым самонепересекающимся кон туром она содержит и область, ограниченную этим контуром, т.е. односв язная область – это область без «дырок».

Для формулировки аналогичной теоремы для случая трех переменных приведем следующее определение.

Определение 20.3. Область в пространстве называется односвязной, если вместе с каждым самонепересекающимся замкнутым кон туром она содержит и некоторую поверхность, натянутую на этот конту р (т.е. поверхность, границей которой является наш контур).

Теорема 20.5. Пусть функции P(x,y,z), Q(x, y,z), R(x,y,z) и их

частные производные ¶P , ¶P , ¶Q , ¶Q , ¶R , ¶R непрерывны в неко-

¶y ¶z ¶x ¶z ¶x ¶y

торой области (Т ) трехмерного пространства. Для того чтобы

ò Pdx +Qdy + Rdz не зависел от формы пути интегрирования для

(AB)

любой кривой (АВ ) М (Т ), необходимо, а в предположении односвязности области и достаточно, чтобы при (x,y,z)О(T )

21. Теория поля

¶P =

¶Q ;

¶Q =

¶R ;

¶R =

¶P .

(20.9)

¶y

¶x

¶z

¶y

¶x

¶z

 

¡ Необходимость. Пусть криволинейный интеграл

ò Pdx +Qdy + Rdz

(AB)

не зависит от формы пути интегрирования. Тогда по теореме 20.3 выражение Pdx +Qdy + Rdz является дифференциалом отнекоторойфункции трех переменных u(x,y,z): Pdx +Qdy + Rdz = du Ю

P =

¶u , Q =

¶u , R =

¶u Þ

¶P =

2u

;

¶Q =

2u

.

 

 

 

¶x

¶y

¶z

¶y

¶y¶x

¶x

¶x¶y

Íî 2u = 2u по теореме (12.9) о смешанных производных. Отсюда

¶y¶x ¶x¶y

¶P = ¶Q . Аналогично проверяются другие равенства:

¶y ¶x

¶Q =

2u

;

¶R =

2u

Þ

¶Q =

¶R ;

 

¶y¶z

¶z

¶z¶y

¶y

 

¶z

¶y

¶R =

 

2u

;

¶P =

2u

Þ

¶R =

¶P .

 

 

¶z¶x

¶x

¶x¶z

¶z

 

¶x

¶z

Достаточность условий (20.9) будет проверена позднее (см. разд. 21.5). x

21. ТЕОРИЯ ПОЛЯ

21.1. Скалярное поле. Производная по направлению. Градиент

Определение 21.1. Пусть каждой точке М некоторой области (Т) трехмерного пространства соответствует скалярная (числ овая) величи- на U =U (M ). Тогда говорят, что в области (Т) задано скалярное поле.

Если положение точки М определять ее координатами в некоторой прямоугольной декартовой системе координат M (x,y,z) , то U =U (M ) =U (x,y,z). Будем полагать, что эта функция имеет непрерывные частные производные.

416

417