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

Основы высшей математики для инженеров 2009

.pdf
Скачиваний:
226
Добавлен:
12.03.2016
Размер:
10.04 Mб
Скачать

5.6. Тройной интеграл

261

П р и м е р. Вычислить DDDzdV , где V — пирамида, ограниченная плос

V

костями x 0, y 0, z 0, x y z 1.

11 x 1 x y

 

 

DDDzdV Ddx Ddy

Dzdz

 

 

 

 

 

 

V

0

0

0

 

 

 

 

 

1

1 x

(1 x y)2

1

(1 x)3

 

1

 

 

 

Ddx D

 

dy D

dx

.

 

 

 

 

2

6

24

Рис. 5.14

0

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

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

Пусть задана точка М в системе координат

OXYZ. Спроектируем ее на плоскость OXY. Тогда

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

определив полярные координаты r и ( точки N в

 

точки М. Числа r,

(,

плоскости OXY и аппликатуM62z

 

z называются

цилиндрическими

координатами

точки М. Переход от декартовых координат к ци

линдрическим

осуществляется

по формулам

(рис. 5.15)

 

 

Рис. 5.15

 

 

 

 

x r cos (,

 

 

y r sin (,

(*)

 

z z.

 

 

Пусть теперь проекция N области интегрирования V на плос кость OXY такова, что интеграл по ней удобнее брать в полярных ко ординатах, тогда делают замену переменных (*) и

 

* r2(()

h(r cos (, r sin ()

DDD f (x, y, z)dV Dd( Drdr

D f (r cos (, r sin (, z)dz,

V

) r1(()

g(r cos (, r sin ()

поскольку здесь dV rd(drdz.

262 Глава 5. Интегральное исчисление функций многих переменных

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

Положение точки в пространстве можно задать также числами (, >, r, где ( — полярный угол проекции точки M на плоскость OXY, > — угол, который составляет радиус вектор точки М с плоско стью OXY (отклонение радиуса вектора от плоскости OXY), r — ради ус вектор точки М. Здесь

0 ( 2',

' > ' , 2 2

0 + .

 

Декартовы координаты точки М связаны с ее

Рис. 5.16

сферическими координатами, очевидно, соотноше

ниями

x r cos > cos (,

y r cos > sin (, z r sin >,

 

рис. 5.16.

 

Можно показать, что dxdydz r2 cos > d(d>dr и

DDD f (x, y, z)dxdydz DDD f [r cos > cos (, r cos > sin (, r sin >]r2 cos > d(d>dr.

V

V

M62

 

 

П р и м е р. Переходя к сферическим координатам, вычислить

DDDx2 y2 z2 dxdydz,

V

где V — шар радиуса R.

Решение. Для шара пределы изменения сферических коорди

нат будут следующими:

 

 

0 ( 2',

' > ' ,

0 r R.

22

Поэтому

 

 

 

'

 

 

 

 

2'

 

 

 

R

 

 

2

 

 

 

DDD

x2 y2 z2 dxdydz Dd( Dd>Drr2 cos > dr 'R4.

V

0

'

0

 

 

 

 

 

 

 

 

 

 

 

 

 

2

5.7. Приложения тройного интеграла

263

5.7. Приложения тройного интеграла

Вычисление объема тела. Вычисление массы тела. Статические момен ты и центр тяжести тела.

5.7.1. Вычисление объема тела

Если подынтегральная функция в тройном интеграле тождествен но равна 1, то есть f (x, y, z) J 1, то

DDD f (x, y, z)dV DDDdV V,

VV

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

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

П р и м е р. Вычислить объем тела, ограниченного поверхностями

 

M62

 

z 1 x2

y2,

z

0.

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2'

 

 

1

 

 

1 r

 

 

 

 

 

 

 

 

 

 

 

 

 

V DDDdxdydz Dd(Drdr

 

Ddz

 

 

 

 

 

0

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

1

 

2

 

 

r

4

 

 

 

1

'

 

 

 

 

 

 

 

 

 

 

 

2'Dr(1 r2)dr 2'

r

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

0

2

 

4

 

0

2

 

Рис. 5.17

 

 

 

 

 

 

 

 

 

5.7.2. Вычисление массы тела

Пусть в пространстве OXYZ задано материальное тело. Рассмотрим его малую часть -V , содержащую точку P(x, y, z). Пусть масса этой части равна -m. Средней плотностью тела называется отношение его

массы к объему, т. е. -m . Если существует предел этого отношения

-V

при стягивании тела в точку P(x, y, z), то он называется плотностью тела в точке Р. Плотность является функцией x, y, z и обозначается +(x, y, z).

264 Глава 5. Интегральное исчисление функций многих переменных

Разобьем теперь наше тело на части -Vi, в каждой из них выберем точку Pi(xi, yi, zi) и положим плотность кусочка -Vi постоянной и рав ной +(xi, yi, zi). Тогда

-mi = +(xi, yi, zi)-Vi,

nn

m F-mi = F+(xi, yi, zi)-Vi.

1 1

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

m DDD+(x, y, z)dz.

V

П р и м е р. Вычислить массу прямого кругового цилиндра высо ты Н и радиуса R, если плотность в каждой его точке пропорциональ на расстоянию r от этой точки до оси цилиндра с коэффициентом

пропорциональности +0.

M62

 

 

Решение. Выберем

 

 

 

 

систему координат так, как показано на рис. 5.18.

 

 

 

m +0DDDrdV.

 

 

 

 

V

 

 

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

 

натах:

 

 

 

 

 

2'

R

H

2'R3H +0.

 

 

m +0 Dd(Drdr Dzdz

Рис. 5.18

 

0

0

0

3

 

 

5.7.3. Статические моменты и центр тяжести тела

Статическим моментом точки относительно плоскости называет ся произведение ее массы m на расстояние до плоскости. Пусть зада но тело с плотностью L(x, y, z). Тогда статические моменты Sxy, Sxz, S yz его относительно плоскостей OXY , OXZ, OYZ вычисляются по следующим формулам:

5.7. Приложения тройного интеграла

265

Sxy DDDz+(x, y, z)dV,

V

Sxz DDD y+(x, y, z)dV,

V

Syz DDD x+(x, y, z)dV.

V

Координаты центра масс x, y, z материального тела определяются формулами

 

 

S yz

 

DDD x+(x, y, z)dV

 

 

 

 

V

,

x

 

 

 

mDDD+(x, y, z)dV

 

 

 

 

 

 

V

 

 

 

 

 

 

 

DDD y+(x, y, z)dV

 

 

 

 

MV62

 

y

 

Sxz

 

V

,

 

 

 

 

m

 

DDD+(x, y, z)dV

 

 

 

 

 

Sxy

 

DDDz+(x, y, z)dV

 

 

 

 

 

 

V

.

z

 

 

 

 

 

mDDD+(x, y, z)dV

V

Если тело однородное, т. е. +(x, y, z) J +0, то

DDD xdV

 

 

V

 

,

x

 

 

 

 

 

V

DDD ydV

 

 

V

 

,

y

 

 

 

 

 

V

DDDzdV

 

 

 

V

 

,

z

 

 

 

 

 

 

 

V

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

266Глава 5. Интегральное исчисление функций многих переменных

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

Задачи, приводящие к криволинейному интегралу. Определение криво линейного интеграла. Теорема существования. Вычисление криволи нейного интеграла. Свойства криволинейного интеграла.

5.8.1. Задачи, приводящие к криволинейному интегралу

Пусть точка P(x, y) движется вдоль некоторой кривой L от точки М к точке N. К точке Р приложена сила F(P), которая меняется по вели чине и направлению при перемещении точки Р, т. е. является функ цией точки Р. Вычислим работу силы F по перемещению точки Р от М к N. Для этого разобьем дугу MN на n частей точками M0 M,

M1, …, Mn N и обозначим через -S i вектор

Mi 1Mi.

Величину

силы

 

F

 

в

точке

~ ~

 

 

 

Fi. Тогда ска

Pi(xi, yi) [Mi 1, Mi] обозначим

M62

(Fi, -S i)

 

 

 

лярное

произведение

можно

рас

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

боты силы F вдоль дуги M

T M

:

 

 

 

 

i 1

 

i

 

 

Рис. 5.19

Ai = (Fi, -S i).

 

 

 

 

Пусть F(P) F(x, y) X(x, y)i Y(x, y) j, где X(x, y) и Y(x, y) — проекции вектора F(P) на оси OX и OY. Вектор -S i можно в координатной форме представить так:

 

-S i -xi i -yi j.

Тогда

 

 

 

~

~

~

~

Ai = X(xi, yi)

-xi Y(xi, yi)-yi,

n

n

~ ~

(1)

 

 

~ ~

A FAi = F[X(xi, yi)-xi Y(xi, yi)-yi].

1

1

 

 

Если существует предел суммы, стоящей справа в (1), при -xi 0 и -yi 0, то этот предел выражает работу силы F по перемещению точки Р вдоль кривой L от точки М к точке N.

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

267

5.8.2. Определение криволинейного интеграла. Теорема существования

Если существует предел суммы (1) при n , -xi, -yi 0, кото рый не зависит ни от способа разбиения дуги на части, ни от выбора

~ ~

точек Pi(xi, yi), то он называется криволинейным интегралом второго рода или криволинейным интегралом по координатам и обозначается:

 

 

 

 

 

(N )

 

 

 

 

D X(x, y)dx Y(x, y)dy D X(x, y)dx Y(x, y)dy.

 

L

 

 

 

(M )

 

 

 

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

 

 

 

 

 

 

 

 

n

~

~

~

~

 

 

 

lim

 

]

 

X(x, y)dx Y(x, y)dy.

[X(x

, y )-x

Y(x

, y )-y

D

n

F

i

i i

i

i i

 

-xi 0

i 1

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

-yi 0

Если бы кривая располагаласьM62в пространстве, а функция F была соответственно функцией координат x, y, z точки Р, то мы пришли бы к интегралу от трех функций:

D X(x, y, z)dx Y(x, y, z)dy Z(x, y, z)dz.

L

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

Наряду с этим при решении таких задач, как вычисление длины плоской или пространственной кривой линии, вычислении массы ма териальной кривой и многих других, вводится понятие и определение криволинейного интеграла первого рода, при этом подынтегральная функция умножается на дифференциал ds длины дуги кривой L, по длине которой проводится интегрирование:

I D f (M)ds,

M(x, y, z) L.

(L)

 

Пусть кривая L задана в параметрической форме уравнениями x ((t), y >(t), где t [), *].

268 Глава 5. Интегральное исчисление функций многих переменных

Теорема существования. Пусть функции ((t) и >(t) непрерывно диффе2 ренцируемы на [), *], а функции X [((t), >(t)], Y [((t), >(t)] непрерывны на [), *] как функции t, тогда для функций X(x, y) и Y(x, y) существует пре2 дел интегральной суммы (1), который не зависит от способа разбиения

~ ~

дуги и от способа выбора точек xi, yi.

Криволинейный интеграл впервые был рассмотрен Клеро (1743).

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

Пусть точки ) t0, t1, t2 , …, tn * таковы, что

x0 ((t0) (()),

y0 >(t0) >()),

x1 ((t1),

y1 >(t1),

C C C C C C C C C C C

C C C C C C C C C C C

xn ((tn) ((*);

yn >(tn) >(*).

FX(M62x~i, ~yi)-xi FY(~xi, ~yi)-yi.

Интегральную сумму по дуге L

запишем в виде

n

n

1

1

Имеем: -xi ((ti) ((ti 1) ( (Pi)-ti, где можем выбирать произвольно, то будем

~

образом, что xi ((Pi). Тогда

~

ti 1 Pi ti. Т. к. точки xi мы считать их выбранными таким

nn

~

~

(2)

FX(xi, yi)-xi FX [((Pi), >(Pi)] ((Pi)-ti.

1

1

 

Справа стоит интегральная сумма для непрерывной функции X [((t), >(t)]( (t) на отрезке [), *].

Переходя к пределу в обеих частях равенства (2) при -xi 0, по

лучим

*

D X(x, y)dx D X [((t), >(t)]( (t)dt.

L)

Аналогично

*

DY(x, y)dy DY [((t), >(t)]> (t)dt.

L)

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

269

Окончательно получаем формулу для вычисления криволинейного интеграла по дуге L:

 

*

D X(x, y)dx Y(x, y)dy D{X [((t), >(t)]((t) Y(((t), >(t)]>(t)}dt.

L

)

Если кривая задана в пространстве уравнениями x ((t), y >(t),

z K(t) при t [), *],

то

D X(x, y, z)dx Y(x, y, z)dy Z(x, y, z)dz

L

*

D{X[((t), >(t), K(t)]((t) Y [((t), >(t), K(t)]>(t) Z[((t), >(t), K(t)]K(t)}dt.

)

5.8.4. Свойства криволинейного интеграла

(N ) (M )

1. D D.

(M ) (N ) M62

Значение криволинейного интеграла зависит от направления, в кото ром проходится кривая L. При изменении направления знак интегра ла меняется на противоположный.

2.Если кривая L состоит из нескольких кривых L1, L2 , …, Lk ,

иинтеграл по каждой кривой Li существует, то существует и интеграл по кривой L, и он равен

D D D D.

L L1 L2

Lk

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

(N )

J D2 xydx y2 dy z2 dz,

(M )

где MN — один виток винтовой линии x cost, y sin t, z 2t от точки M(1, 0, 0) до точки N(1, 0, 4').

Решение. Здесь t меняется от 0 до 2'.

2'

J D[2 cost sin t( sin t) sin2 t cost 4t2 2]dt

 

0

 

 

 

 

 

 

 

 

 

 

2'

 

 

 

 

2'

 

 

 

sin3 t

8

 

 

 

 

64

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

3

 

 

3

 

D( sin

t cost 8t

 

)dt

 

 

 

t

 

 

 

 

 

 

' .

 

 

 

 

 

 

0

 

 

 

3

3

 

 

 

 

0

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

270 Глава 5. Интегральное исчисление функций многих переменных

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

J D(x2 y2)dx 2 xydy

 

MN

 

 

 

 

 

 

 

 

 

вдоль параболы y x2

от точки (–1, 1) до точки (1, 1).

Решение. Здесь x x, y x2 и x можно рассматривать как параметр.

1

1

3

5

 

 

1

4

 

 

J D(x2 x4 2 x x2 2 x)dx D(x2

4x4)dx

x

4

x

 

 

2

.

 

 

 

1

1

3

 

5

 

1 15

 

 

 

 

5.9. Формула Грина – Остроградского. Независимость криволинейного интеграла от пути интегрирования

Формула Грина – Остроградского. Независимость криволинейного ин теграла от пути интегрирования

5.9.1. ФормулаM62Остроградского – Грина

Рассмотрим в плоскости XOY область N, ограниченную некоторой кривой L. Пусть в этой области заданы функции P(x, y) и Q(x, y), не прерывные вместе со своими частными производ

ными HP и HQ . Тогда имеет место следующая

 

Hx

 

Hy

 

 

 

формула:

 

 

 

 

 

 

HQ

HP

 

Рис. 5.20

DD

Hx

 

Hy

dN DP(x, y)dx Q(x, y)dy,

 

 

 

 

L

 

N

 

 

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

Эта формула была получена знаменитым русским математиком М. В. Остроградским (1801–1861) и английским физиком и математи ком Д. Грином (1793–1841). Она является частным случаем более об щей формулы, полученной Остроградским, с которой мы познако мимся позже.

Доказательство. Предположим, что область N имеет вид, изобра женный на рис. 5.21.