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

Вища математика / В.М.-iнтеграли

.pdf
Скачиваний:
53
Добавлен:
12.02.2016
Размер:
1.29 Mб
Скачать

площині XOY. Якщо область D обмежена двома кривими, полярні рівняння яких ρ=ρ1(ϕ) і ρ=ρ2(ϕ) (ρ1(ϕ)≤ρ2(ϕ)) і променями ϕ=ϕ1, ϕ=ϕ2, то

 

 

 

ϕ2

ρ2 (ϕ )

∫∫ f (x, y)dxdy =∫∫ f ×cosϕ, ρ ×sin ϕ )ρd ρdϕ = dϕ

f (ρ cosϕ , ρ sin ϕ )ρd ρ .

D

σ

 

ϕ1

ρ1 (ϕ )

 

Якщо область містить початок координат, то

 

 

ρ (ϕ )

 

 

∫∫ f (x, y)dxdy = dϕ f (ρ cosϕ, ρ sin ϕ)ρd ρ , де ρ=ρ(ϕ) -полярне рівняння

 

D

0

0

 

кривої, що обмежує область D.

 

 

Полярні координати зручно використовувати, якщо область є круг або

його частина.

Обчислити ∫∫ ydxdy , якщо область D обмежена верхньою

 

Приклад 5.

дугою кола x2+y2=ax

D

 

і віссю OX.

 

Розв’ язання. Введемо полярні координати x=ρcosϕ, y=ρsinϕ, тоді рівняння

кола прийме вид: ρ=acosϕ. Кут ( змінюється від 0 до π

) (Рис. 12.6).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

При кожному фіксованому значенні ϕ ρ

 

 

 

 

 

 

 

 

 

 

ρ=acosϕ

 

змінюється від 0 до ρ=acosϕ. Тоді за

 

 

 

 

 

 

 

 

 

 

 

 

 

формулою ( 12.1. 6) маємо:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X, ρ

 

 

 

 

 

 

 

 

 

 

π

 

 

 

a cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

2a

 

 

 

 

∫∫ ydxdy =dϕ

 

ρ 2 sin ϕd ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

0

 

 

 

 

0

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 12.6.

 

 

 

 

 

 

 

 

 

π

 

 

 

 

a cosϕ

 

 

 

 

 

ρ

 

 

a cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

3

 

 

a

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= sin ϕdϕ

ρ 2 d ρ = sin ϕ

 

 

 

 

dϕ = =

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

0

 

 

 

0

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приклад 6. Обчислити I = ∫∫ 1-

 

 

-

y

dxdy , де область обмежена еліпсом.

a

2

2

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Розв’ язання. Введемо так звані узагальнені полярні координати,

 

 

покладаючи x=aρcosϕ, y=bρsinϕ. Якобіан перетворення

 

 

 

 

 

 

 

 

 

 

I =

 

a cosϕ

 

- aρ sin ϕ

 

= abρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b sin ϕ

 

bρ cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кут

 

ϕ міняється від

0 до 2π. Рівняння еліпса в узагальнених полярних

координатах ρ=1, тому ρ змінюється від 0 до 1;

 

 

 

 

 

 

 

 

 

 

 

 

 

1-

x2

-

y2

 

= 1- ρ 2 cos2 ϕ - ρ 2 sin2 ϕ = 1- ρ 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

(1- ρ 2 )2

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I = ab dϕ

1- ρ

 

ρd ρ

= ab ×

-

 

 

 

 

×

 

 

 

 

 

 

 

 

 

=

 

π ab .

 

 

 

 

 

 

 

 

 

 

 

2

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приклад 7.

У теорії ймовірностей і математичній статистиці використовується

інтеграл Пуассона I = ex2 dx .

0

- 180 -

Первісна для підінтегральної функції f (x) = ex2 не виражається в елементарних функціях, однак даний невласний інтеграл може бути обчислений за допомогою подвійного інтеграла.

Оскільки визначений інтеграл не залежить від позначення змінної інтегрування, запишемо:

∞ ∞ ∞ ∞

I 2 = ex2 dx × ey2 dy = ∫ ∫ e−( x2 + y2 )dxdy .

0 0 0 0

Областю інтегрування є перший квадрант системи координат (Рис.12.7)

ϕ=π/2

Y

ϕ=0

0

 

X

 

 

 

 

 

 

 

 

Рис.12.7

 

 

 

 

 

 

 

 

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

0 £ ϕ £ π , 0 < ρ < ¥, x2 + y2 = ρ 2 .

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

π

 

 

π

 

A

π

 

 

lim

a

 

 

 

 

 

 

 

2

 

 

 

 

 

1

I 2

= dϕ e− ρ 2 ρd ρ = ϕ

 

 

2

lim

e− ρ 2 ρd ρ =

-

e− ρ 2 d (-ρ 2 )=

 

 

 

 

 

0

0

 

 

0

A→+∞

0

2

 

2 A→+∞

0

π lim

(eA2 e0 ) = π ,

 

де lim eA2 =0 ,

 

 

 

 

 

 

 

 

 

 

 

4 A→+∞

4

 

 

 

A→∞

 

 

 

 

 

Отже, інтеграл Пуассона дорівнює:

 

 

 

π

 

ex2 dx =

.

 

0

2

 

 

 

 

12.2. Застосування подвійних інтегралів

Обчислення об'ємів тіл.

Циліндричне тіло, обмежене знизу областю D площини ХОУ, зверху –

поверхнею z = f (x, y) ,

збоку – циліндричною поверхнею

з

твірними,

паралельними осі ОZ і

з напрямною – межею області

D,

має об’єм

V = ∫∫ f ( x, y)dxdy.

 

 

 

D

Приклад. Подвійним інтегруванням знайти об'єм тіла, обмеженого

циліндрами y =

 

, y = 2

 

 

x

x і площинами z = 0, x + z = 6 (рис. 12.8).

Розвязання. Тіло

обмежене зверху площиною x + z = 6, знизу

площиною z = 0 і двома циліндрами y = x та y = 2 x , що проектують його на площину xOy в область D, обмежену прямою x = 6 і параболами y = x та y = 2 x (рис. 12.9).

- 181 -

z

0

y = x

y = 2 x

y

D

y

y = 2 x

D

y = x

0

6

x

x Рис. 12.8

Рис. 12.9

При цьому змінна x змінюється від 0 до 6; при будь-якому значенні х із зазначеного проміжку x < y < 2 x . Об'єм знаходимо за формулою

V = ∫∫ f (x, y)dxdy. У даному випадку z = f (x, y) = 6 − x . Тоді

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

2

 

 

x

 

6

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

V = ∫∫(6 − x)dxdy = (6 − x)dx

dy = (6 − x)y

 

2

 

 

 

dx = (6 − x)(2

 

 

 

)dx =

 

 

 

x

 

 

 

x

x

 

x

 

 

 

 

 

 

 

D

0

 

 

 

 

x

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

3 2

 

 

 

5 2

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

72

 

 

 

48

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (6

 

x x x )dx = 6

 

 

 

 

 

 

 

 

 

 

= 24 6 −

 

6 =

 

 

.

 

 

 

3 2

5 2

5

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приклад. Подвійним інтегруванням знайти об'єм тіла, обмеженого

циліндрами

x2 + y 2 = x

та

 

 

x2 + y 2 = 2x ,

 

 

параболоїдом

z = x2 + y 2 і

площинами x + y = 0, x – y = 0, z = 0 .

Розв’язання. Дане циліндричне тіло обмежене зверху поверхнею

z= x2 + y 2 . Об'єм його

v= ∫∫ f (x, y)dxdy = ∫∫zdxdy = ∫∫(x2 + y 2 )dxdy .

D

D

D

Побудуємо область D, що

зображує

проекцію тіла на площину xOy

(рис.12.10). Область D обмежена колами

 

y

y=x

 

 

ρ = 2 cos ϕ

 

D

 

0 0.5 1

2

x

ρ = cos ϕ

y=-x

Рис. 12.10

x 2 + y 2 = x

- 182 -

 

 

 

1

2

+ y

2

=

 

1 2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

або x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 + y 2 = 2x або (x −1)2 + y 2 = 1 і прямими y = - x, y = x.

 

У полярних координатах рівняння кола має вигляд: ρ = cos ϕ і ρ = 2 cos ϕ ;

рівняння прямих: ϕ = − π і ϕ = π . Знаходимо:

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

4

 

 

4

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

2 cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 cos ϕ

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

ρ4

 

 

v = dϕ

 

ρ3dρ = (внаслідок симетрії)

= 2

 

dϕ =

 

 

π

 

cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

4

cos

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

4

(16cos4

ϕ − cos4 ϕ)dϕ =

15

4

 

1 + cos 2ϕ

2

 

 

 

 

 

2

 

2

 

 

 

 

2

 

 

dϕ =

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

15

4

 

 

 

 

 

 

 

 

 

)dϕ =

15

4

 

 

+ 2cos 2ϕ +

1 + cos 4ϕ

 

8

(1 + 2cos 2ϕ + cos2

 

8

1

 

2

dϕ =

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

15

 

3

ϕ + sin 2ϕ +

sin 4ϕ 4

=

15 3π

 

 

 

 

 

 

 

 

 

 

8

 

 

2

 

 

8

 

 

 

 

 

8

 

+ 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приклад. Знайти об'єм тіла, вирізаного з кулі радіусом R, прямим

круговим циліндром радіусом R/2, твірна якого проходить через центр кулі.

 

 

 

 

z

 

 

 

 

 

 

 

 

 

Розвязання. Помістимо

початок

 

 

 

 

 

 

 

 

 

 

координат у центр кулі, вісь OZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

направимо уздовж твірної циліндра, а

 

 

 

 

 

 

 

 

 

 

 

 

вісь

OX

 

 

 

уздовж

діаметра основи

 

 

 

 

 

 

 

 

 

 

 

y

циліндра (рис. 12.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рівняння

сфери

має

вигляд:

 

 

 

 

 

 

 

 

 

 

 

 

x 2 + y 2 + z2 = R2 ,

тоді

в

першому

 

 

0

 

 

R/2

 

 

R

октанті

 

z =

 

R2 x2 y 2

і внаслідок

 

 

 

 

 

 

x симетрії об’єм тіла, вирізаного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

циліндром з кулі, буде рівним:

 

 

Рис. 12.11.

 

 

 

 

 

 

 

 

 

V = 4∫∫

R2 x2

y2 dxdy.

Проекція тіла на площину хОу збігається з кругом x2 + y2 Rx . Щоб

обчислити

отриманий

інтеграл, зручно

 

перейти до

 

полярних

координат.

 

 

 

 

 

 

 

 

 

 

 

 

- 183 -

 

 

 

 

 

 

 

 

 

 

 

 

 

Рівняння кола x2 + y2 = Rx в полярних координатах має вид ρ = R cos ϕ . Кут

ϕ змінюється від 0 до π/2 (враховуємо

симетрію),

ρ змінюється в межах

0 ≤ ρ ≤ R cos ϕ . Переходячи до полярних координат, одержимо:

 

 

 

 

 

 

 

 

 

 

π

 

R cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V = 4dϕ

R2 − ρ2

ρdρ ;

 

 

 

 

 

 

 

π

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π / 2

 

 

 

V

1

2

2

 

2

 

2

3 / 2

 

R cos ϕ

 

 

R3

 

3

 

 

 

 

 

 

 

 

 

 

= −

 

 

(R

 

− ρ

 

)

 

0

 

dϕ = −

 

(sin

 

ϕ − 1)dϕ =

4

2

3

 

 

 

3

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R3 π / 2

(1 − cos2 ϕ ) d (cosϕ ) +

π R3

 

R3

 

cos3 ϕ

=

 

 

 

=

 

cosϕ −

 

 

3

 

6

3

 

 

0

 

 

 

π

 

 

 

 

 

3

 

 

V =

4

R3

2

 

 

 

 

 

 

 

Отже,

 

 

 

.

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Обчислення маси неоднорідної пластини.

Пластина, що займає область D у площині ХОУ має масу

m = ∫∫ ρ (x, y)dxdy.

π / 2

+

π R3

=

R3

π

2

 

0

 

 

 

 

 

6

3

3

 

 

 

2

 

 

і має щільність ρ ( x, y ) ,

D

Обчислення статичних моментів і моментів інерції пластини.

Статичні моменти пластини відносно осей ОХ і ОУ відповідно

рівні

Мx = ∫∫ y ρ(x, y)dxdy ,

Мy = ∫∫ x ρ (x, y)dxdy.

D

D

Момент інерції пластини відносно початку координат визначається за формулою I0 = ∫∫(x2 + y2 )ρ(x, y)dxdy.

D

Моменти інерції пластини відносно осей ОХ і ОУ будуть відповідно

Ix = ∫∫ y2 ρ (x, y)dxdy ,

I y = ∫∫ x2 ρ (x, y)dxdy.

D

D

Знаходження центра ваги пластини.

Координати центра ваги xc та yc неоднорідної пластини дорівнюють,

відповідно, відношенням статичних моментів відносно осей ОУ й ОХ до маси пластини:

x =

M y

=

∫∫ xρ (x, y)dxdy

; y =

M

 

=

∫∫ yρ(x, y)dxdy

;

D

x

D

 

 

M

∫∫ ρ (x, y)dxdy

 

 

∫∫ ρ (x, y)dxdy

c

 

c

M

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

D

 

Якщо пластина однорідна:

- 184 -

 

∫∫xdxdy

 

∫∫ ydxdy

 

xc =

D

; yc =

D

.

S

S

 

 

 

Приклад. Знайти момент інерції квадрата зі стороною а, поверхнева щільність якого пропорційна відстані до однієї із сторін квадрата, відносно вершини, що належить даній стороні.

Розвязання. Нехай квадрат розташований у площині xOy; одна з його вершин належить початку координат, а дві інші збігаються з осями координат. Відзначимо, що момент інерції не залежить від вибору системи координат. Шуканий момент інерції дорівнює моменту інерції квадрата щодо початку координат з поверхневою щільністю kx, де k – коефіцієнт пропорційності (беремо поверхневу щільність, пропорційну відстані до осі Oy). Тоді

 

 

 

 

2

 

 

 

2

 

 

 

a

 

a

 

2

 

2

 

 

 

 

a

 

2

 

 

 

y

3

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I0 = ∫∫(x

 

 

+ y

 

)kxdxdy = k xdx(x

 

+ y

 

)dy = k x x

 

 

y

+

 

 

 

 

 

 

 

dx =

 

 

 

 

 

 

3

 

 

 

 

D

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

2

 

 

a

3

 

3

 

a

3

 

 

 

x

4

 

 

x

2

a

3

 

 

 

5ka

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= k x x

 

a +

 

 

 

dx = k x

 

a + x

 

 

dx = k a

 

 

+

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

.

 

 

3

 

3

 

 

 

 

6

 

 

 

 

 

12

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

4

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

12.3. Потрійні інтеграли і їх обчислення в декартовій

системі координат

Нехай у правильній, замкненій обмеженій області V задана обмежена

функція f (x, y, z) . Потрійним інтегралом від функції

f (x, y, z) по області V

називається:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

lim

f (ξi , ηi i ) Vi = ∫∫∫ f (x, y, z)dV = ∫∫∫ f (x, y, z)dxdydz .

 

λ→0 i =1

 

 

 

 

 

V

V

 

 

 

 

 

 

 

 

 

 

 

 

Обчислення потрійного інтеграла в

Z

 

 

 

 

 

 

 

z =z2 (x,y)

декартовій системі координат

 

 

 

 

 

 

 

 

 

 

зводиться до обчислення подвійного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

інтеграла по проекції D об’єму V на

 

 

 

 

 

 

 

 

 

 

будь-яку координатну площину (у

 

 

 

 

 

 

 

 

 

 

даному випадку ХОУ) і

 

 

 

 

 

 

 

 

 

z = z1(x,y)

внутрішнього інтеграла по третій

 

 

 

 

 

 

 

 

змінній (змінна z) (рис.12.12).

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Y

Внутрішній інтеграл береться від

 

 

 

 

 

 

 

 

 

нижньої межі z = z1,у) області V до

a

 

 

 

 

 

 

 

 

 

її верхньої межі z = z2,у)

 

 

 

 

 

 

 

 

(передбачається, що область є

 

 

 

y =y1(x)

 

 

 

 

 

 

D

 

b

 

 

 

 

 

 

 

 

y =y2(x)

правильною в напрямі осі Oz):

z2 (x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

Рис.12.12

 

 

 

 

∫∫∫ f (x, y, z)dxdydz =∫∫dxdy

f (x, y, z)d

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

D

z1 ( x, y)

 

 

 

- 185 -

 

 

 

Враховуючи правила обчислення подвійного інтеграла, останню формулу можна переписати таким чином:

 

b

y2 ( x)

z2 ( x, y)

∫∫∫ f (x, y, z)dxdydz = dx

dy

f (x, y, z)dz.

V

a

y1 ( x)

z1 ( x, y)

Якщо область V є неправильною, то її розбивають на скінченне число правильних областей і обчисляюють інтеграл, використовуючи властивість адитивності потрійного інтеграла.

Якщо областю інтегрування є прямокутний паралелепіпед, обмежений площинами х = а, х = b < b), у = с, у = d < d), z = m, z = n (m < n), то межі інтегрування будуть сталими, тобто

b d n

∫∫∫ f (x, y, z)dV = dxdy f (x, y, z)dz.

V

a c m

У випадку, якщо f (x, y, z) ≡ 1,

то потрійний інтеграл чисельно дорівнює

об'єму області інтегрування.

Приклад. Обчислити об'єм тіла, обмеженого поверхнями z = x2 + y 2 , z = 2(x2 + y 2 ), y = x, y = x2 .

Розвязання. Тіло обмежене площиною y = x, циліндром y = x2 і

параболоїдами обертання z = x2 + y 2 , z = 2(x2 + y 2 ). Нехай D - проекція тіла на площину xOy; у даному випадку проектуючими поверхнями є площина

y

 

 

y = x

1

 

y = x2

 

D

 

0 1

x

Рис. 12.13.

 

y = x і циліндр y = x2 (рис. 12.13). Тоді об'єм

- 186 -

 

1

x

2(x 2 + y 2 )

1

x

 

 

2(x 2 + y

2 )

1

x

V = ∫∫∫dv = dx dy

dz

= dx

z

 

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

 

 

x 2 + y 2

 

 

 

v

0

x 2

x 2 + y 2

0

x 2

 

 

 

 

0

x 2

1

 

 

3

 

 

x

 

 

 

 

2

y +

y

 

 

 

x

 

3

 

 

x 2

 

 

 

 

 

0

 

 

 

 

 

 

1

 

4x

3

 

x

6

 

3

 

dx =

 

x4

 

dx =

.

 

 

 

 

 

0

 

3

 

3

 

35

 

 

 

 

 

 

 

 

 

 

Приклад. Обчислити масу тіла, обмеженого циліндром x2 = 2y і площинами z = 0, 2y + z = 2, якщо в кожній його точці об'ємна щільність чисельно дорівнює аплікаті цієї точки.

Розвязання. Циліндричне тіло (рис. 12.14) обмежено зверху площиною z = 2 - 2y, що перетинається з площиною z = 0 по прямій y = 1. Маса тіла, що займає область V, обчислюється через потрійний інтеграл: m = ∫∫∫ ρ (x, y, z)dxdydz , де ρ (x, y, z,) - об'ємна щільність.

V

В нашій задачі ρ (x, y, z) = z , тому

 

1

 

2 y

 

2−2 y

1

 

2 y

 

1

 

 

m = ∫∫∫zdxdydz = dy

 

dx

zdz = 2(1 − y)2 dy

 

dx = 4(1 − y)2

 

 

2 ydy

V

0

 

 

0

0

 

 

0

 

 

2 y

2 y

 

 

z

0

1

y

x

Рис.12.14

12.4. Потрійні інтеграли і їх обчислення в циліндричній і сферичній системах координат

Заміна змінних у потрійному інтегралі: нехай функція f (x, y, z) неперервна в області V і формули.

 

x = x (u,v, w); y = y (u,v, w), z = z (u,v, w)

встановлюють

взаємно однозначну відповідність між точками M (x, y, z)

області V і точками M '(u,v, w) деякої області V ', тоді

∫∫∫ f (x, y, z)dxdydz = ∫∫∫ f (x (u,v, w), y (u,v, w), z (u,v, w))

 

J

 

dudvdw ,

 

 

 

 

V

V '

- 187 -

де J - абсолютне значення якобіана.

У циліндричній системі координат положення точки визначається полярними координатами ϕ, ρ та аплікатою z (рис. 12.15), а формули, що зв'язують прямокутні і циліндричні координати мають вигляд: x = ρ cos ϕ ; y = ρ sin ϕ ; z = z. Модуль Якобіана дорівнює J = ρ , тоді

∫∫∫ f (x, y, z)dxdydz = ∫∫∫ f (ρcosϕ,ρsin ϕ, zdρdϕdz.

V

 

 

V '

 

z

 

 

z

 

 

 

М

θ М

 

 

 

 

 

0

ρ

z

0 ρ

 

ϕ

y

ϕ

y

 

 

 

x

 

 

x

 

 

Рис. 12.15

Рис. 12.16

 

Якщо областю інтегрування є круговий циліндр із віссю Oz, то потрійний інтеграл за цією областю в циліндричній системі координат матиме сталі межі по всіх змінних, тобто

R H

∫∫∫ f (x, y, z)dV = dϕρdρ f (ρcos ϕ,ρsin ϕ, z)dz .

V

0 0

0

Сферичні координати ρ,ϕ,θ зв'язані з прямокутними координатами співвідношеннями:

x = ρ sin θ cos ϕ, y = ρ sin θ sin ϕ, z = ρ cos θ , J = ρ2 sin θ .

Перехід у потрійному інтегралі до сферичних координат здійснюється за формулою:

∫∫∫ f (x, y, z)dxdydz = ∫∫∫ f (ρsin θcos ϕ,ρsin θsin ϕ,ρcos θ)ρ2 sin θdρdϕdθ .

V V

Очевидно, якщо областю інтегрування є куля з центром на початку координат і радіусом R, то потрійний інтеграл за цією областю в сферичній системі координат матиме сталі межі інтегрування по всіх змінних, тобто

πR

∫∫∫f (x, y, z)dV = sin θdθ dϕρ2 f (ρsin θcos ϕ,ρsin θsin ϕ,ρcos θ)dρ .

V

0

0 0

Нижче на конкретних прикладах проілюстровані правила для розставлення меж інтегрування в циліндричній і сферичній системах координат і показані їх геометричні та фізичні застосування.

- 188 -

Приклад. Обчислити

∫∫∫xyzdxdydz , де V – частина області, яка обмежена

 

V

сферою х2 + у2 + z2 = 4

і параболоїдом х2 + у2 = 3z, розташована в першому

октанті (рис.12.17).

 

2

 

3

3 x2+y2=3

Рис.12.17

Розвязання. Перший спосіб. Обчислення інтеграла в декартовій системі координат. Перед тим, як проектувати об'єм V на площину ХОY, знайдемо лінію перетину сфери і параболоїда. Для цього розв’яжемо систему рівнянь

 

2

+ y

2

+ z

2

 

 

 

 

 

x

 

 

 

= 4

z 2

+ 3z − 4 = 0 z = 1 x2 + y2

= 3 ,

 

 

2

+ y

2

= 3z

 

 

x

 

 

 

 

 

 

 

 

 

 

 

тобто поверхні перетинаються по колу радіуса R = 3 , що лежить у площині z = 1. Поверхня проектується на площину хОу у чверть круга даного радіусу, що знаходиться в першій чверті.

Одержимо:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3− x 2

 

 

 

 

 

 

4− x 2 y 2

 

 

 

 

 

 

3− x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

4− x2 y

2

 

 

∫∫xyzdxdydz = xdx

 

 

 

 

 

 

ydy

 

 

 

 

 

zdz =

xdx

yz 2

 

 

dy =

 

 

 

 

 

 

 

 

 

 

 

 

2

 

x2 + y 2

 

 

D

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

x2 + y 2

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ y2 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3− x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

xdx

y 4

x2 y2

 

 

 

 

 

 

 

 

 

 

dy

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

 

x2 y2

y4

 

 

x4 y2

 

 

 

 

x2 y4

 

 

 

y6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 3− x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

x 2 y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

0

 

 

 

 

 

 

2

 

 

 

 

4

 

 

 

18

 

 

 

18

 

 

 

 

54

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

 

 

3

1 3

 

− 2 x 3 +

 

x 5

+

 

x 7

 

 

 

 

=

2 7

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

d x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

3 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 4

 

 

 

 

 

 

 

 

 

 

 

 

5 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Другий спосіб. Обчислення інтеграла в циліндричній системі координат:

- 189 -

Соседние файлы в папке Вища математика