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

Письменный Д.Т. Конспект лекций по высшей математике - Полный курс

.pdf
Скачиваний:
26441
Добавлен:
10.03.2016
Размер:
20.05 Mб
Скачать

Переходя к полярным координатам, находим:

V = ~JJ(7 -

r 2 )r · drdr.p -

JJ(1 + r2 )r · drdr.p =

 

 

 

D

 

 

 

D

 

 

 

 

 

1

211"

1

 

 

211"

1

 

 

 

= 3

Jdr.p J(7r -

r 3 ) dr -

Jdr.p J(r + r 3 ) dr =

 

 

о

о

 

 

о

о

 

 

 

= ! (~ - !) .r.p,211" -

(! + !) .r.pl21r =

13 • 271" -

~. 27Г = ~71".

3 2 4

о

2

4

о

12

4

3

 

При.мер 53.4.

Найти массу, статические мо­

 

 

 

менты Sx и Sy и координаты центра тяжести фигу­

 

у

 

 

1

 

ры, лежащей в первой четверти, ограниченной эл-

 

 

 

 

 

2

=

1 и координатными осями (см.

 

 

 

липсом 4 + у2

 

 

2 х

рис. 224). Поверхностная плотность в каждой точке

 

о

 

 

 

фигуры пропорциональна произведению координат

 

Рис. 224

 

точки.

Q Решение: По формуле (53.6) находим массу пластинки. По условию,

'У= "((х;у) = k · ху, где k - коэффициент пропорциональности.

т = Jj kxy dx dy = k j2 х·dx

 

J уdy = 2k j2 хdx · у2lp;4

=

 

 

 

р;

 

 

 

 

 

 

 

 

 

 

 

 

D

О

 

 

О

 

 

 

 

 

О

 

 

О

 

 

 

 

 

= - · -

2

х(

4 -

х

2

) dx

= -

2

- -

4

 

2

= -

 

 

 

 

 

 

 

 

 

1

 

k

1 /

 

 

 

 

 

k (

 

х )

 

k

 

2

4

 

 

 

 

 

 

 

8

 

 

4

 

о

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим статические моменты пластинки:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г:-7

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

v·-4

у

dy = ... =

15k,

 

 

Sx = JJу · kxy dx dy = k Jхdx

!

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

4

 

 

 

 

D

О

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

2

 

 

v·-4

уdy = ... =

 

 

 

 

 

Sy = jj х · kxy dx dy = k j

х2 dx

 

!

 

185 k.

 

 

 

D

О

 

 

 

 

о

 

 

 

 

 

 

 

 

 

Находим координаты центра тяжести пластинки, используя формулы

х

 

-

Si. и у

-

Si.. х

-

16 у

-

8

 

с

-

т

с -

т .

с -

15,

с -

15.

 

390

§ 54. ТРОЙНОЙ ИНТЕГРАЛ

54.1. Основные понятия

Обобщением определенного интеграла на случай функции трех пе­

ременных является так называемый «тройной интеграл».

Теория тройного интеграла аналогична теории двойного интегра­

ла. Поэтому изложим ее в несколько сокращенном виде.

Пусть в замкнутой области V пространства Oxyz задана 11епре­

рывная функция и= f(x; у; z). Разбив область V сеткой поверхностей

на п частей Vi (i = 1, п) и выбрав в каждой из них произвольную точ­ n

ку Mi(xi; Yii zi), составим интегральную сумму 2:: f (xi; Yii Zi)дVi для

i=l

функции f (х; у; z) по области V (здесь дVi -объем элементарной обла-

сти Vi).

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

область» Vi стягивается в точку (т. е. диаметр области di стремится к нулю, т.е. di-+ О), то его называют троi1'Н:ым интегралом от функции

и= f(x; у; z) по области V и обозначают

 

jjj f(x;y;z) ·dxdydz (или

jjj f(x;y;z)dv).

v

v

Таким образом, по определению, имеем:

 

n

 

jjj f(x;y;z) ·dxdydz = JhIIJo Lf(xi;Yi;zi)дVi =

 

v

(maxd;-+0) i=l

(54.1)

 

= JJJ f(x;y;z)dv.

 

 

 

v

 

Здесь dv = dx dy dz -

элемент объема.

 

Теорема 54.1 (существования). Если функция и= f(x; у; z) непре­

рывна в ограниченной замкнутой области V, то предел интегральной

суммы (54.1) при п -+ оо и maxdi -+ О существует и не зависит

ни от способа разбиения области V на части, ни от выбора точек

Mi(xi; Yi; zi) в них.

Тройной интеграл обладает теми же свойствами, что и двойной

интеграл:

1. jjj c·f(x;y;z)dv=c· jjj f(x;y;z)dv, c-const.

v

v

391

2.JJJU1(x;y;z) ± f2(x;y;z))dv = v

=JJJ fi(x;y;z)dv± JJJ f2(x;y;z))dv.

v

 

v

 

 

3.

JJJ f(x;y;z)dv = JJJ f(x;y;z)dv + JJJ J(x;y;z)dv, если V

=

 

V

V1

Vz

 

= V1 UV2 ,

а пересечение Vi

и V2 состоит из границы, их разделяющей.

4.

J JJ

j(x; у; z) dv ~ О,

если в области V функция f(x; у; z) ~

О.

 

v

 

 

 

Если в области интегрирования f(x; у; z) ~ 1i0(x; у; z), то и

 

 

 

JJJ f(x;y;z)dv ~ J JJ 1i0(x;y;z)dv.

 

 

 

v

v

 

5.

jJJ dv = V, так как в случае f(x;y;z) = 1 любая интегральная

 

V

n

 

 

сумма имеет вид ЕдVi = V и численно равна объему тела.

 

i=l

6.Оценка тройного интеграла:

т· V::::; JJJ f(x;y;z)dv::::; М ·V, v

где т и М - соответственно наименьшее и наибольшее значения функ­

ции f(x;y;z) в области V.

7. Теорема о среднем значении: если функция f(x;y; z) непрерыв­ на в замкнутой области V, то в этой области существует такая точка

Мо(хо; Уо; zo), что

J jJ f(x; у; z) dv = f(xo; Уо;zo) · V, v

где V - объем тела.

54.2. Вычисление тройного интеграла в Аекартовых

координатах

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

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

Пусть областью интегрирования V является тело, ограниченное

снизу поверхностью z =

z1 ( х; у), сверху - поверхностью z = z2 ( х; у),

причем z1 (x; у) и z2(x; у)

(z1 (х; у) ::::; z2(x; у)) - непрерывные функции в

замкнутой области D, являющейся проекцией тела на плоскость Оху

392

(см. рис. 225). Будем считать область V - прави.11ьно1J. в направле­

нии оси Oz: любая прямая, параллельная оси Oz, пересекает границу

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

области V функции f (х; у; z)

имеет место формула

 

JJJ f(x;y;z)dv= JJ(

z2(x;y)

 

J

f(x;y;z)dz)ds,

(54.2)

V

D

z1(x;y)

 

 

сводящая вычисление тройного интеграла к вычислению двойного ин­

теграла от однократного (доказательство формулы (54.2) не приводим).

При этом сначала вычисляется внутренний интеграл по переменной z при постоянных х и у в пределах изменения z. Нижней границей инте­ грала является аппликата точки А - точки входа прямой, параллель­

ной оси Оz в область V, т. е. z = z1 ( х; у); верхней границей - аппликата

точки В - точки выхода прямой из области V, т. е. z = z2 (x; у). Ре­ зультат вычисления этого интеграла есть функция двух переменных:

х и у.

у

 

 

. . . .. ' . ..

. . .. . . '

 

.

.................

 

 

...................

 

:::::::::::::::::n:::::::::::::::::

 

 

.. .. . . . .. . . . .

' . .

О

а

 

ьх

Рис. 225

 

Рис. 226

 

Если область D ограничена линиями х =а, х = Ь (а< Ь), у = ip1 (х) и у= ip2(x), где ip1(x) и ip2(x) - непрерывные на отрезке [а, Ь] функции, причем ip1(х) ~ ip2(x) (см. рис. 226), то, переходя от двойного интеграла

по области D к повторному, получаем формулу

JJJ f(x;y;z)dxdydz=

Ь

<р2(х)

z2(x;y)

 

j dx

J dy

J

f(x;y;z)dz,

(54.3)

V

а

1 (х)

z1(x;y)

 

 

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

393

Заме'Чания.

1. Если область V более сложная, чем рассмотренная, то ее следует разбить на конечное число таких областей (правильных), к которым

можно применить формулу (54.3).

2. Порядок интегрирования в формуле

(54.3), при определенных условиях, может

быть иным.

Пример

54.1. Вычислить

JJJ + z) dxdydz, v

у

где V ограничена плоскостями х = О, у = О,

z = 1, х +у+ z = 2 (рис. 227).

Q Решение: Область V является правильной

Рис. 227

в направлении оси Oz (как, заметим, и в на-

правлении осей Ох и Оу). Ее проекция на плоскость Оху является пра­

вильной в направлении оси Оу (и оси Ох). Согласно формуле (54.3),

имеем:

Jjj(x+z)dxdydz= Jdx Jx dy

2J-y(x+z)dz=

v

о

о

1

 

 

1

1-х

z2 2-х-у

 

 

= Jdx

J dy · (xz + 2 ) /1

=

оо

1

1-х

 

(

2 - х2- у

)2

1

= Jdx

J (-

х2 - ху - х +

 

-

2) dy =

оо

= 1dx

(

2

у2

-

(2 - х - у)з

1

) ,1-х

=

 

ху - х у -

х2

6

- "2у

0

о

= ](х-х2-х2+хз_х(l;х)2 -~+ (2~х)з -~+~x)dx=

о

=а.х; _2 _х; +х44 -~(~2 _2.х; +х:)-~·х-(2;4х)4 )\:=

3

2

1

1

2

1

16

1

= 4- 3+ 4- 24 - 3 - 24 + 24 = 4.

394

54.З. Замена переменных в тройном интеграле. Вычисление тройного интеграла

в цилиндрических и сферических координатах

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

ных.

 

Пусть совершена подстановках =

ip(и;v;w), у= ф(и;v;w), z =

= х(и; v; w). Если эти функции имеют

в некоторой области V* про­

странства Ouvw непрерывные частные производные и отличный от ну-

ля определитель

дх

дх

дх

ди

дv

дw

I(u ·, v·, w) = !!.JL

l!JJ.

!!.JL

ди

дv

дw

дz

дz

дz

ди

дv

дw

то справедлива формула замени переменних в тройном интеграле:

jjj f(x;y;z) dxdydz = v

= jjj f(ip(и;v;w);ф(и;v;w);x(и;v;w)) · II(и;v;w)ldиdvdw. (54.4)

Здесь I(u;v;w) - определитель Якоби, или якобиан преобразования (примем без доказательства).

Для вычисления тройного интеграла

 

часто используют так называемые цилин­

z

 

дрические координаты.

 

Положение точки М(х; у; z) в простран­

стве Oxyz можно определить заданием трех

чисел r, <р, z, где r - длина радиуса-вектора

проекции точки М на плоскость Оху, <р - угол, образованный: этим радиусом-векто­

ром с осью Ох, z - аппликата точки М (см.

рис. 228).

Эти три числа (r, <р, z) называются ци­

линдри"tескими координатами точки М.

M(x;y;z)

M(r; ер; z) z

у

Рис. 228

Цилиндрические координаты точки связаны с ее декартовыми ко­

ординатами следующими соотношениями:

J х = r · cos <р, у = r · sin <р, z = z J

(r ~ О, 'РЕ [О; 2п], z Е IR).

395

Возьмем в качестве и; v, w цилиндрические координаты r, rp, z и вычислим якобиан преобразования:

 

дх

дх

дх

 

 

 

 

дr

&rp

дz

cos rp

-r sin ер

о

 

 

 

 

I(r·, m·, z) =

S!.1L

!!JL

S!.1L

sin rp

rcosrp

О =r;;::: О.

r

дr

д<р

дz

 

 

1

 

дz

дz

дz

о

о

 

дr д<р дz

 

Формула замены переменных (54.4) принимает вид

 

JJJ f(x;y;z)dxdydz = JJJ f(rcoscp;rsincp;z)rdrdcpdz.

(54.5)

v

 

Таким образом, вычисление тройного интеграла приводится к инте­

грированию по r, по rp и по z аналогично тому, как это делается в

декартовых координатах.

Заме-ч,ание. К цилиндрическим координатам бывает удобно перей­ ти в случае, если область интегрирования образована цилиндрической

поверхностью.

При.мер 54. 2. Вычислить JJJz · dx dy dz, где V - область, огра­

v

ниченная верхней частью конуса х2 + у2 = z2 и плоскостью z = 1.

Q Решение: На рис. 229 изображена область интегрирования V. Вы­

числим интеграл путем перехода к цилиндрическим координатам: х =

= r · coscp, у= r · sincp, z = z. Здесь dxdydz = r · drdcpdz. Уравнение

конуса примет вид r 2 cos2 cp+r2 sin2 ср = z 2 , т. е. z = r. Уравнение окруж­

ности х2 + у2 = 1 (границы области D) запишется так: r = 1. Новые

переменные изменяются в следующих пределах: r - от О до 1, ер - от

О до 2к, а z - от r до 1 (прямая, параллельная оси Oz, пересекающая

область D, входит в конус z = r и выходит из него на высоте z = 1). Таким образом, согласно формуле (54.5), получаем:

JJJz dx dy dz = JJJz · r · dr drp dz =

27Т

1

 

1

 

Jdcp Jr d1·

Jz dz =

 

V

V

 

 

О

О

 

r

 

27Т

1

2 1

2ir

1

1

2

 

= Jdcp Jr dr · z2 1r =

Jdcp Jr · ( 2 - r2 ) dr =

 

о

о

 

о

о

 

 

 

 

 

271"

r2

r4

1

1

27r

1 2ir

 

Jdcp · ( 4 - 8) lo= 8

Jdrp = 8 · 'Pjo

=

оо

Заметим, что, не переходя к цилиндрическим координатам, полу-

чим:

1

~

1

 

JJJz · dx dy dz = J dx

J

dy

J z · dz.

V

-1

-~

 

Jx2+y2

 

396

z

z=l

1 у

Рис. 229

z

М(х· у· z)

, M(pi <р; 8)

1

1

1

1

у

Рис. 230

Сфери'Ческими координата.ми точки М(х; у; z) пространства Oxyz

называется тройка чисел р, ip, (},где р - длина радиуса-вектора точки

М, ifJ - угол, образованный проекцией радиуса-вектора ОМ на плос­ кость Оху и осью Ох, (} - угол отклонения радиуса-вектора ОМ от

оси Oz (см. рис. 230).

Сферические координаты р, ip, (} связаны с декартовыми коорди­

натами х, у, z соотношениями:

1 х = pcosip ·sinO, у= psinip ·sinO, z = pcos(} 1

(р ~ о, о ~ ip ~ 27r, о~ (} ~ 1Г).

В некоторых случаях вычисление тройного интеграла удобно про­

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

 

 

COS ipSill (}

-

р sin ер sin (}

р cos ер cos (}

 

 

I(r; <р; z) = sin sin (}

р cos ер sin (}

р sin ер cos (}

 

 

 

 

 

cos (}

 

о

 

 

-psinO

 

 

.

.

(}

1 sin ерsin (}

рsin ерcos (}

1

+

 

 

= р sm ер sш

·

cos (}

. (}

 

 

 

 

 

 

 

 

 

-psm

 

 

 

 

 

 

 

 

 

 

.

 

 

1 cos ерsin (}

pcOSipCOS01-

 

 

 

 

 

+ р cos ер sш (} .

 

cos (}

-psin(}

-

= рsin ерsin О(-рsin ер sin2 (} - рsin epcos2 О)+

 

 

 

 

 

 

+ рcos ерsin О(-р cos ерsin2 (} -

р cos ерcos2 О) =

= - р2 sin2 ifJ sin (} · 1 -

р2 cos2 ер sin (}

· 1 = - р2 sin (} · 1 = - р2 sin (},

то

JJJ f(x;y;z)dxdydz = v

397

= JJJ f(pcosrpsinO; psinrpsinO; pcosO) · р2 sinO · dpdrpdO. (54.6) v•

Замечание. Переходить к сферическим координатам удобно, когда

область интегрирования V есть шар (уравнение его границы х2 + у2 + + z2 = R2 в сферических координатах имеет вид р = R) или его часть, а также если подынтегральная функция имеет вид f(x2 + у2 + z2 ).

Прuмер 54.3.

Вычислить

 

 

 

JJJ

dx · dy · dz l ,

 

 

 

v

1 + (х2 + у2 + z2)2 -

 

 

где V - шар х2 + у2 + z2 :::;

1.

 

 

Q Решение: Вычислим интеграл путем перехода к сферическим коор­

динатам: х = pcosrpsinO, у= psinrpsinO, z = pcosO. Тогда

 

 

dV = dx dy dz = р2 sin О dp drp dO.

 

 

Граница области V -

сфера и ее уравнение имеет вид р = 1, подынте­

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

+

(12 3 ; 2 ,

 

 

1

р )

т. е. ~l. Новые переменные изменяются в следующих пределах: р-

1

 

 

 

 

 

 

 

 

 

 

 

от О до 27Г, О -

 

 

 

 

 

 

 

 

 

от О до 1, rp -

от О до 7Г. Таким образом, согласно фор-

муле (54.6),

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

27Г

1

 

2

 

 

 

I = JJJ 1 + рз

· р2 sinOdpdrpd() = JsinOdO

Jdrp J1 : рз dp =

 

v

 

 

о

 

о

о

 

 

 

 

 

27Г

1

 

1

 

2,.-

1

 

 

 

= Jsin ОdO Jdrp ( Зln 11 + p3 I) j0 = Jsin ОdO

J

3ln 2 drp =

 

о

о

 

 

о

 

 

о

 

 

 

 

 

1

,.-

127Г

2

 

,.-

 

 

 

 

 

= Зln 2

Jsin () dO · rp 0

= ; ln 2 Jsin

ОdO =

 

 

 

 

 

о

 

 

о

 

 

 

 

 

 

 

 

 

 

27Г

1n2(-cos0)

!,.-

47Г

1n2.

 

 

 

 

= 3

0 =

3

 

54.4. Некоторые приложения тройного интеграла

Объем тела

Объем области V выражается формулой V = JJJdv или

v

V = JJJdx dy dz - в декартовых координатах, v

398

V

= JJJr dr di.p dz - в цилиндрических координатах,

 

v

V

= JJJр2 sin (} dp d'fJ d(} - в сферических координатах.

 

v

Масса тела

Масса тела т при заданной объемной плотности 'У вычисляется с

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

т = JJJ 1(x;y;z)dxdydz, v

где 1(х; у; z) - объемная плотность распределения массы в точке

M(x;y;z).

Статические моменты

Моменты Sxy, Bxz, Syz тела относительно координатных плоско­ стей Оху, Oxz, Oyz вычисляются по формулам

Sxy = JjJz · 1(х;у; z) dv,

Syz = JjJх · 1(х;у;z) dv,

v

v

Bxz= jjj y·1(x;y;z)dv. v

Центр тяжести тела

Координаты центра тяжести тела V находятся по формулам

Syz

Bxz

Sxy

Хе=-, Ус=-,

Zc = --.

т

т

т

Моменты инерции тела

 

 

Моменты инерции тела относительно координатных плоскостей

вычисляются по формулам

lyz = JJj х2 ·1(x;y;z) dv,

Ixy = JJJ z2 ·1(x;y;z)dv,

v

 

v

lxz= JJJ y2 ·1(x;y;z)dv, v

а моменты инерции относительно координатных осей:

lx = JjJ2 +z2)·1(х;у;z) dv,

ly = JJJ2 +z2)·1(х;у;z) dv,

v

v

lz = JJJ2 2) ·1(x;y;z)dv_. v

Пример 54.4. Найти объем тела, ограниченного поверхностями z = х2 + у2 и z = 1.

399