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

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

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

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

251

Рис. 5.8

Рис. 5.9

Каждой точке P(xo, yo) из области D соответствует точка P (uo, vo) из области D . Рассмотрим в области D линии u const, v const. В области D этим линиям будут соответствовать, вообще говоря, неко торые кривые. Будем считать, что D разбита на прямоугольники ли ниями u const и v const. (При этом прямоугольники, прилегающие к границе области, мы не рассматриваем.) Тогда область D разобьется на некоторые криволинейные четырехугольники. Рассмотрим в плос

также через S и S . M62

кости D прямоугольную площадку -S , ограниченную прямыми

u const, u -u const, v const, v -v const, и соответствующую ей нелинейную площадку -S в D. Площади этих площадок обозначим

- - Очевидно, что

-S -u -v.

Пусть в области D задана непрерывная функция z f (x, y). Тогда в D также задана непрерывная функция

z F(u, v) f [((u, v), >(u, v)].

Если мы рассматриваем интегральные суммы для f (x, y) в D, то, оче видно,

nn

F f (xi, yi)-Si FF(ui, vi)-Si.

(2)

i 1

i 1

 

Если бы удалось выразить -S через -S , то мы выразили бы инте гральную сумму для f через интегральную сумму для F.

Для этого вычислим -S, т. е. площадь криволинейного четырех угольника в области D. Координаты его вершин:

P1(x1, y1):

x1 ((u, v),

y1 >(u, v);

P2(x2, y2):

x2 ((u, v -v),

y2 >(u, v -v);

P3(x3, y3):

x3 ((u -u, v -v),

(3)

y3 >(u -u, v -v);

P4(x4, y4):

x4 ((u -u, v),

y4 >(u -u, v).

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

При вычислении площади -S будем считать линии P1P2 , P2 P3, P3P4 и P4P1 попарно параллельными, а приращение функции будем заменять ее дифференциалом, т. е. пренебрегать бесконечно малыми более вы сокого порядка, чем -u и -v. Тогда формулы (3) примут вид

x1 ((u, v),

 

 

y1 >(u, v),

 

 

x2

((u, v)

H(

-v,

 

 

y2

>(u, v)

H>

-v,

 

 

 

 

 

 

 

 

 

 

 

 

Hv

 

 

 

 

Hv

 

 

x3

((u, v)

H(

-u

H(

-v,

y3

>(u, v)

H>

-u

H>

-v,

 

 

 

 

 

 

 

 

Hu

Hv

 

 

Hu

Hv

x4

((u, v)

H(

-u,

 

 

y4

>(u, v)

H>

-u.

 

 

 

 

 

 

 

 

 

 

 

Hu

 

 

 

 

Hu

 

 

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

 

 

 

 

 

 

 

 

 

 

M62

 

 

 

 

 

 

 

изведения [P2 P1, P2 P3]:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-S = |(x3 x )(y3 y1) (x3 x1)(y3 y2)|

 

 

H(

 

 

H>

 

H>

 

 

 

 

 

H(

 

H(

 

 

 

H>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-u

 

 

-u

 

 

 

-v

 

 

 

-u

 

 

-v

 

 

 

-u

 

Hu

 

 

Hv

 

Hv

 

 

 

 

 

 

 

Hu

 

 

 

 

 

 

Hu

 

 

 

 

Hu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H(

 

H(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H(

 

H>

-u-v

H(

 

H>

-v-u

 

 

 

Hu Hv

 

 

-u-v.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hu Hv

 

 

Hv

 

 

Hu

 

 

 

 

 

 

 

H>

 

H>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hu

 

Hv

 

 

 

 

 

Определитель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

 

u

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

v

 

 

 

 

 

 

 

 

 

 

 

называется якобианом функций ((u, v), >(u, v). Итак, -S = | J| -S . При -S 0 это равенство становится точным. С учетом последнего соот ношения равенство (2) приобретает вид

nn

F f (xi, yi)-Si = FF(ui, vi)| J | -Si.

i 1 i 1

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

253

После перехода к пределу при n и diam -Si 0 получим:

 

DD f (x, y)dxdy DDF (u, v)| J | dudv.

(4)

(D)

(D )

 

 

Формула (4) позволяет свести вычисление интеграла по области D к

вычислению интеграла по области D , что часто упрощает задачу.

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

 

Наиболее часто в приложениях встречается замена, позволяющая

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

ласти интегрирования D пересекается прямой, проходящей через на

чало координат, не более чем в двух точках, и ограничена кривыми

+ E1(B), + E2(B).

 

 

 

 

Пусть в области D

 

задана функция

f (x, y)

 

и нужно вычислить

 

 

 

 

DD f (x, y)dxdy.

 

 

D

 

 

 

 

Координаты x, y любой точки области D связа

 

ны с ее полярными

координатами + и B сле

 

 

M62

 

дующим образом:

 

 

Рис. 5.10

 

x + cosB,

 

y + sin B.

(6)

 

Формулы (5) определяют замену переменных с якобианом J:

 

 

 

Hx

 

Hx

 

 

 

 

 

 

 

 

 

 

 

| J |

 

HB

 

H+

 

 

 

+ sin B

cosB

 

 

 

+(sin2 B cos2 B)

 

+.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hy

 

 

+ cosB

sin B

 

 

 

 

 

Hy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HB

 

H+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в полярных координатах получим

*E2(B)

DD f (x, y)dxdy DdB

D[f (+ cosB, + sin B)]+d+.

(D)

)

E1(B)

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

DD4 x2 y2 dxdy,

(D)

если D — круг радиуса 2 с центром в начале координат.

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

Сделаем замену x + cosB,

y + sin B. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2'

2

 

 

 

 

 

2

 

3

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4

)

2

 

 

 

16'

 

 

 

2

2

 

 

2

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

DD 4

x

y dxdy DdBD

4

+ +d+ 2'

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

3

 

(D)

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.5. Приложения двойного интеграла

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

Мы уже отметили два приложения двойного интеграла. 1. Объем тела

V DD f (x, y)dxdy.

(N)

2. Площадь области

 

 

 

 

 

 

 

 

N DDdN DD+ d+dB .

 

(N)

 

 

 

(N)

 

 

M62

5.5.1. Масса плоской фигуры

Пусть в плоскости XOY задана материальная площадка N с поверх ностной плотностью L(x, y). Под плотностью L(x, y) здесь понимается предел, к которому стремится отношение массы -m элемента площад ки, окружающей точку (x, y) к площади этого элемента площадки, ко

гда diam -N 0:

L(x, y) lim -m .

diam -N 0 -N

Разобьем площадку N на n частей -N1, -N2 , …, -Nn; в каждой из площадок -Ni выберем точку Pi(xi, yi) и будем считать плотность на

всей площадке -Ni постоянной и равной L(xi, yi). Тогда

 

-mi = L(xi, yi)-Ni,

(1)

где -mi — масса площадки -Ni. Пользуясь (1), запишем

 

nn

m F-mi = FL(xi, yi)-Ni.

(2)

i 1

i 1

 

5.5. Приложения двойного интеграла

255

За массу m пластинки примем предел, к которому стремится инте

гральная сумма (2) при стремлении diam -Ni

0 и при n :

 

 

n

 

 

 

m

lim

FL(xi, yi)-Ni DDL(x, y)dN.

 

n

i 1

(

N

)

 

diam -Ni 0

 

 

 

 

 

 

 

 

5.5.2. Статические моменты, центр тяжести плоской фигуры

Под статическим моментом Sx материальной точки P(x, y), распо ложенной в плоскости XOY и имеющей массу m, относительно оси OX понимают произведение массы точки на ее ординату:

Sx m y.

Аналогично

S y m x.

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

Пусть вновь задана материальнаяM62площадка N с поверхностной плотностью L(x, y). Для определения ее статических моментов посту пим следующим образом: разобьем ее на части -Ni; в каждой из этих частей выберем точку с координатами Pi(xi, yi) и вычислим в ней плот ность L(xi, yi). За массу кусочка -Ni приближенно примем

-mi = L(xi, yi)-Ni.

Заменим каждую площадку -Ni материальной точкой (xi, yi) и мас сой -mi, равной массе площадки. Статические моменты этой точки относительно осей дадут приближенные значения статических мо

ментов -Sxi и -S yi площадки -Ni:

-Sxi = yi -mi = yi L(xi, yi)-Ni; -S yi = xi -mi = xi L(xi, yi)-Ni.

Так как статический момент всей площадки равен сумме статических моментов малых площадок -Ni, то мы получим следующие прибли

женные равенства:

n

Sx = Fyi L(xi, yi)-Ni,

1

n

S y = Fxi L(xi, yi)-Ni.

1

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

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

Sx

lim

 

n

 

diam -Ni

S y

lim

 

n

 

diam -Ni

n

Fyi L(xi, yi)-Ni DD yL(x, y)dN,

 

i 1

(

N

0

 

)

 

 

 

 

n

 

 

 

Fxi L(xi, yi)-Ni DD xL(x, y)dN.

 

i 1

 

N

0

 

( )

 

 

 

Как известно из механики, координаты x и y центра тяжести мате риальной пластинки определяются по следующим формулам:

 

 

S y

,

 

 

 

S

x

 

x

y

 

 

 

 

mm

 

 

 

DD xL(x, y)dN

 

 

 

 

 

 

DD yL(x, y)dN

 

 

 

N

 

 

 

 

 

;

 

 

 

 

 

 

N

 

.

x

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M620

 

 

 

 

 

DDL(x, y)dN

 

 

 

 

 

 

DDL(x, y)dN

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

Если пластинка однородная, т. е. L(x, y) J L

, то

 

 

 

 

 

 

L0DD xdN

 

 

 

DD xdN

 

 

 

 

 

 

 

 

 

N

 

N

 

 

;

 

 

 

 

 

x

 

 

 

 

 

 

 

L0DDdN

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DD ydN

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

.

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

П р и м е р. Найти координаты центра тяжести однородной пло щадки, ограниченной параболой y 4 x2 и осью OX.

Решение. Так как фигура симметрична относительно оси OY, то x 0. Так как пластинка однородная, то по формуле (3)

 

 

 

 

 

2 4 x2

1

2

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DD ydxdy

 

Ddx

D ydy

D(4 x2)2 dx

D(4 8x2 x4)dx

 

256

 

 

 

 

 

 

 

 

2

 

2

 

 

8

 

 

 

 

 

 

 

 

15

 

 

 

 

 

 

N

 

2

0

 

2

 

 

2

 

 

 

 

 

 

 

.

y

 

 

 

 

 

 

 

 

4 x2

 

2

 

 

 

 

 

 

2

32

 

 

 

 

 

DDdxdy

2

 

 

2

 

 

 

3

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

N

 

Ddx

Ddy

 

 

D(4

x )dx

 

 

4x

1

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

2

0

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.5. Приложения двойного интеграла

 

 

257

5.5.3. Площадь поверхности

 

 

 

 

Пусть в области N задана непрерывная функция f (x, y). Нужно оп

ределить площадь поверхности, определяемой уравнением z f (x, y);

f (x, y) предполагается непрерывной вместе с частными производными

1 го порядка.

 

 

 

 

 

 

 

 

 

 

Лемма. Пусть N — площадь проекции некоторой плоской фигуры с

площадью O на некоторую плоскость,

тогда

 

 

 

 

 

 

N O cos (,

 

 

 

 

 

 

где ( — угол между плоскостью фигуры и плоскостью проекции.

 

Доказательство. Этот факт известен для треугольников и, следова

тельно, он справедлив для многоугольников. Для доказательства лем

мы достаточно вписать в нашу фигуру много

 

 

 

 

угольник и перейти к пределу при неограни

 

 

 

 

ченном увеличении числа его сторон.

 

 

 

 

 

 

Для вычисления

площади

поверхности

 

 

 

 

S разобьем N на n площадок. Через -Si

 

 

 

 

обозначим ту часть поверхности, которая

 

 

 

 

проектируется в -Ni. Выберем

в -Ni точку

 

 

 

 

 

 

координатой

 

 

 

 

Pi(xi, yi) и через точку M62Mi с

 

 

 

 

zi f (xi, yi) проведем к поверхности S каса

 

 

 

 

тельную плоскость O. Рассмотрим ту ее

 

 

 

 

часть, которая проектируется в -Ni. Обозна

 

 

 

 

чим ее -Oi. По лемме

 

 

 

 

 

 

 

Рис. 5.11

 

-O

-Ni ,

 

 

 

 

 

 

 

 

 

i

cos Li

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Li — угол между площадкой -Ni и -Oi, который равен углу между

нормалью к поверхности и осью OZ. Согласно 4.9.2,

 

cos Li

 

 

1

 

 

,

 

 

 

 

 

 

 

 

 

 

1 f 2(x , y ) f 2

(x , y )

 

 

 

 

x

 

i

i

y

i

i

 

 

 

поэтому

 

 

 

 

 

 

 

 

 

 

-O 1 f 2(x , y ) f 2

(x , y ) -N

i

.

(4)

i

x

i

i

y

i

i

 

 

За площадь S поверхности примем предел

 

 

 

 

 

 

n

 

b

 

 

S lim

F-Oi

lim

F 1 fx 2(xi, yi) fy2(xi, yi) -Ni.

n

 

n

 

 

 

diam -Ni 0 i 1

diam -Ni 0 i 1

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

Это выражение есть интегральная сумма для функции 1 fx 2 fy2 ,

которая интегрируема, так как fx и fy по предположению непрерыв ны. Следовательно,

S DD1 fx 2(x, y) fy2(x, y) dN.

N

П р и м е р. Определить площадь поверхности той части параболоида вращения z x2 y2 , которая вырезается цилиндром x2 y2 4. Здесь f (x, y) x2 y2 , N — окружность с центром в точке (0, 0) радиуса 2:

 

 

 

 

 

fx 2 x,

 

fy 2 y,

 

S DD 1 4x2 4y2 dxdy.

Вычисление

интеграла

 

 

N

полярных

координатах:

x + cosB,

 

проведем в

 

y + sin B.

 

 

 

M62

 

2'

2

 

 

 

 

 

 

'

 

 

 

 

 

 

 

S DdB D 1 4+2 +d+

 

 

 

 

 

 

(17

17 1).

Рис. 5.12

 

 

0

0

 

 

 

 

 

 

 

 

 

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

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

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

Пусть в пространстве задана область V, объем которой равен V. Пусть в каждой точке этой области определена функция u f (P)

f (x, y, z). Выполним следующие действия.

1.Разобьем область V на n частей -V1, -V2 , …, -Vn таких, что

n

F-Vi V.

i1

2.В каждом малом элементе -Vi выберем точку Pi(xi, yi, zi).

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

259

3. Составим сумму

nn

F f (Pi)-Vi F f (xi, yi, zi)-Vi,

11

которая называется интегральной суммой. 4. Рассмотрим

 

n

lim

F f (xi, yi, zi)-Vi.

n

1

diam -Vi 0

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

DDD f (x, y, z)dV.

V

Теорема существования. Для всякой функции u f (x, y, z), непрерыв2 ной в замкнутой ограниченной области V, существует тройной интеграл, т. е. существует предел интегральныхM62сумм при неограниченном увеличе2 нии числа малых элементов -Vi и при стягивании каждого из них в точ2 ку, который не зависит ни от способа разбиения объема V на части -Vi, ни от способа выбора точек Pi(xi, yi, zi).

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

1.DDDkf (x, y, z)dV kDDD f (x, y, z)dV.

V V

2. DDD[f (x, y, z) ((x, y, z)]dV DDD f (x, y, z)dV DDD((x, y, z)dV.

V V V

3. Если в области интегрирования f (x, y, z) 0, то

DDD f (x, y, z)dV 0.

V

4.Если в области интегрирования f (x, y, z) ((x, y, z), то

DDD f (x, y, z)dV DDD((x, y, z)dV.

VV

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

5. Теорема о среднем значении. Если функция f (x, y, z) непрерывна в замкнутой ограниченной области V, то в этой области существует та2 кая точка P0(x0, y0, z0), что

DDD f (x, y, z)dV f (x0, y0, z0) V,

V

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

6. Свойство аддитивности.

DDD fdV DDD fdV DDD fdV DDD fdV,

V

V1

V2

Vk

где V1 V2 Vk V .

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

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

сверху и снизу, соответственно, поверхностями z h(x, y), z g(x, y).

Пусть это тело проектируется на плоскости XOY в область N, ограни

ченную кривыми

y (1(xM62), y (2(x) [(2(x) (1(x)] и

прямыми x a,

x b. Тогда, как

можно

показать, в

случае

непрерывной функции

f (x, y, z)

 

 

h(x, y)

 

 

 

 

 

 

 

 

DDD f (x, y, z)dV DDdN

D f (x, y, z)dz.

 

 

V

N

g(x, y)

 

 

 

 

Для того, чтобы вычислить тройной инте

 

грал, нужно сначала вычислить интеграл

 

 

 

h(x, y)

 

 

 

 

 

D f (x, y, z)dz,

 

 

 

 

g(x, y)

 

 

 

считая x и y постоянными. При этом получим

 

функцию от x и y. Ее нужно проинтегрировать

 

по области N. Пользуясь формулой для вычис

 

ления двойного интеграла, получим

Рис. 5.13

 

 

b

(2(x) h(x, y)

 

 

DDD f (x, y, z)dz Ddx Ddy

D f (x, y, z)dx.

V

a (1(x) g (x, y)