Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul...

14
Curs de Matematici Speciale - Cadastru + IC German… a 2016 Curs 5: Integrale triple Fie func‚ tia f : D R 3 ! R; continu… a‚ si not… am integrala tripl… a pe D astfel I = ZZZ D f (x; y; z) dxdydz: 1. Dac… a D este paralelipipedul dreptunghic [a; b] [c; d; ] [e; f ] atunci ZZZ D f (x; y; z) dxdydz = Z f e " Z b a Z d c f (x; y; z) dy ! dx # dz: Exemplul 1: S… a se calculeze ZZZ cos (x + y + z) dxdydz; n cazul n care este cubul h 0; 2 i h 0; 2 i h 0; 2 i : Solu‚ tie: Avem 1

Transcript of Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul...

Page 1: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Curs de Matematici Speciale - Cadastru + IC Germana 2016

Curs 5: Integrale triple

Fie functia f : D ⊂ R3 → R, continua si notam integrala tripla pe D astfel

I =

∫∫∫D

f (x, y, z) dxdydz.

1. Daca D este paralelipipedul dreptunghic

[a,b]× [c,d, ]× [e, f ]

atunci ∫∫∫D

f (x, y, z) dxdydz =

∫ f

e

[∫ b

a

(∫ d

c

f (x, y, z) dy

)dx

]dz.

Exemplul 1: Sa se calculeze∫∫∫

Ω

cos (x+ y + z) dxdydz, în cazul în care

Ω este cubul [0,π

2

]×[0,π

2

]×[0,π

2

].

Solutie: Avem

1

Page 2: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

I =

∫ π

2

0

∫ π

2

0

∫ π

2

0

cos (x+ y + z) dy

dx

dz =

=

∫ π

2

0

∫ π

2

0

[sin(x+ z +

π

2

)− sin (x+ z)

]dx

dz =

=√

2

∫ π

2

0

∫ π

2

0

cos(x+ z +

π

4

)dx

dz =

=√

2

∫ π

2

0

[sin

(z +

4

)− sin

(z +

π

4

)]dz =

= 2

∫ π

2

0

cos(z +

π

2

)dz = −2

∫ π

2

0

sin z =

= 2 cos z |π20 = 2 (0− 1) = −2.

Facem observatia ca daca functia este cu variabile separate, adica

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

atunci ∫∫∫D

f (x, y, z) dxdydz =

∫ b

a

u (x) dx ·∫ d

c

v (y) dy ·∫ f

e

w (z) dz.

Exemplul 2: Calculati integrala tripla

I =

∫∫∫D

zdxdydz,

undeD : [1, 2]× [1, 3]× [0, 4] .

Solutie:

2

Page 3: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Observam ca domeniul este un paralelipiped dreptunghic de dimensiuni 1, 2si 4, iar functia este cu variabile separate. Atunci avem:

I =

∫∫∫V

zdxdydz =

∫ 2

1

dx ·∫ 3

1

dy ·∫ 4

0

zdz = 1 · 2 · z2

2|40= 16.

2. Daca

D = (x, y, z) | x ∈ [a, b] , u (x) ≤ y ≤ v (x) , α (x, y) ≤ z ≤ β (x, y)

atunci∫∫∫D

f (x, y, z) dxdydz =

∫ b

a

[∫ v(x)

u(x)

(∫ β(x,y)

α(x,y)

f (x, y, z) dz

)dy

]dx.

3. DacaD = (x, y, z) | z ∈ [a, b] , (x, y) ∈ Ez

atunci ∫∫∫D

f (x, y, z) dxdydz =

∫ b

a

∫∫Ez

f (x, y, z) dxdy

dz.

4. Daca

D = (x, y, z) | (x, y) ∈ E, α (x, y) ≤ z ≤ β (x, y)

3

Page 4: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

atunci ∫∫∫D

f (x, y, z) dxdydz =

∫∫E

(∫ β(x,y)

α(x,y)

f (x, y, z) dz

)dxdy.

Exemplul 3: Calculati I =

∫∫∫V

ydxdydz unde V este limitat de planele:

x = 0, y = 0, z = 0, 2x+ y + 2z = 2.Solutie:Ecuatia planului prin taieturi este:

2x+ y + 2z = 2 |: 2⇔ x

1+y

2+z

1= 1

si avem figura

Din ecuatia planului 2x+ y + 2z = 2 avem z = 1− x− y

2de unde rezulta

I =

∫∫∫V

ydxdydz =

∫∫D

(∫ 1−x−y2

0

ydz

)dxdy =

=

∫∫D

(yz |1−x−

y2

0

)dxdy =

∫∫D

y(

1− x− y

2

)dxdy

4

Page 5: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

unde D este proiectia piramidei (V ) în planul xOy si avem figura

si scriem:

I =

∫∫D

ydxdy −∫∫D

xydxdy − 1

2

∫∫D

y2dxdy =

=

∫ 1

0

(∫ 2−2x

0

ydy

)dx−

∫ 1

0

(∫ 2−2x

0

xydy

)dx− 1

2

∫ 1

0

(∫ 2−2x

0

y2dy

)dx =

= 2

∫ 1

0

(1− x)2dx− 2

∫ 1

0

x (1− x)2dx− 1

2· 8

3

∫ 1

0

(1− x)3dx =

=2

3− 1

6− 1

3=

1

6.

5. Schimbarea de variabilaDaca x = x (u, v, w)

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

; (u, v, w) ∈ ∆; x, y, z ∈ C1∆;

D (x, y, z)

D (u, v, w)6= 0

atunci ∫∫∫D

f (x, y, z) dxdydz =

=

∫∫∫∆

f (x (u, v, w) , y (u, v, w) , z (u, v, w)) ·∣∣∣∣D (x, y, z)

D (u, v, w)

∣∣∣∣ dudvdw.Cazuri particulare de schimbari de variabilaa) Coordonatele sferice x = ρ cos θ sinϕ

y = ρ sin θ sinϕz = ρ cosϕ

ρ ≥ 0θ ∈ [0, 2π]ϕ ∈ [0, π]

∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = ρ2 sinϕ.

5

Page 6: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Exemplul 4: Calculati integrala tripla I =

∫∫∫V

zdxdydz, unde

V : x2 + y2 + (z − 1)2 ≤ 4, z ≥ 0.

Solutie:

Domeniul este o semisfera (z ≥ 0) centrata în (0, 0, 1) si de raza 2, deci vomfolosi coordonatele sferice: x = ρ cos θ sinϕ

y = ρ sin θ sinϕz = 1 + ρ cosϕ

ρ ∈ [0, 2]θ ∈ [0, 2π]

ϕ ∈z≥0

[0,π

2

] ∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = ρ2 sinϕ

de unde rezulta

∫∫∫V

zdxdydz =

∫∫∫∆

(1 + ρ cosϕ) ρ2 sinϕdρdθdϕ =

=

∫∫ ∫[0,2]×[0,2π]×[0, π2 ]

ρ2 sinϕdρdθdϕ+

+

∫∫ ∫[0,2]×[0,2π]×[0, π2 ]

ρ3 sinϕ cosϕdρdθdϕ.

6

Page 7: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Deoarece functia din ambele integrale este cu variabile separate, iar domeniuleste un paralelipiped dreptunghic, avem:

I =

∫ 2

0

ρ2dρ ·∫ 2π

0

dθ ·∫ π

2

0

sinϕdϕ+

+

∫ 2

0

ρ3dρ ·∫ 2π

0

dθ ·∫ π

2

0

sinϕ cosϕdϕ

=16π

3+ 4π =

28π

3.

b) Coordonatele sferice generalizate x = aρ cos θ sinϕy = bρ sin θ sinϕz = cρ cosϕ

ρ ≥ 0θ ∈ [0, 2π]ϕ ∈ [0, π]

∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = abcρ2 sinϕ.

Exemplul 5: Calculati integrala tripla∫∫∫

V

zdxdydz, unde

V :x2

4+y2

9+ z2 ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0.

Solutie:

Domeniul este o optime de elipsoid (primul octant) de semiaxe 2, 3, 1, decivom folosi coordonatele sferice generalizate

x = 2 cos θ sinϕy = 3 sin θ sinϕz = 1 · ρ cosϕ

ρ ∈ [0, 1]

θx≥0∈y≥0

[0,π

2

]ϕ ∈z≥0

[0,π

2

]∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = 6ρ2 sinϕ

si avem

7

Page 8: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

∫∫∫∆

ρ cosϕ · 6ρ2 sinϕdρdθdϕ =

∫∫ ∫[0,1]×[0, π2 ]×[0, π2 ]

6ρ3 sinϕ cosϕdρdθdϕ =

= 6 ·∫ 1

0

ρ3dρ ·∫ π

2

0

dθ ·∫ π

2

0

sinϕ cosϕdϕ =

=3π

8.

c) Coordontele cilindrice x = aρ cos θy = bρ sin θz = z

ρ ≥ 0θ ∈ [0, 2π]z ∈ [h1, h2]

D (x, y, z)

D (ρ, θ, z)= abρ.

Exemplul 6: Calculati I =

∫∫∫V

√x2 + y2zdxdxydz, unde V este marginit

de cilidrul x2 + y2 = 1 si de planele z = 0, z = 4.Solutie:Folosim coordonatele cilidrice si avem: x = ρ cos θ, ρ ∈ [0, 1]y = ρ sin θ, θ ∈ [0, 2π]z = z, z ∈ [0, 4]

,D (x, y, z)

D (ρ, θ, z)= ρ, ∆ = [0, 1]× [0, 2π]× [0, 4]

de unde rezulta

I =

∫∫∫V

√x2 + y2zdxdxydz =

∫∫∫∆

√ρ2z · ρdρdθdz =

=

1∫0

ρ2dρ ·2π∫0

dθ ·4∫

0

zdz =1

3· 2π · 8 =

16π

3.

6. Aplicatii geometrice si fizicea) Volumul unui domeniu masurabil, compact D ⊂ R3

V ol (D) =

∫∫∫D

dxdydz.

Exemplul 7:Calculati volumul corpului (V ) marginit de sfera x2+y2+z2 = 4 si interiorul

conului√x2 + y2 = z.

Solutie:

8

Page 9: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Din intersectia conului√x2 + y2 = z cu sfera x2 + y2 + z2 = 4 obtinem

z =√

2⇒ x2 + y2 = 2

adica un cerc centrat în (0, 0) si de raza r =√

2 si avem figura:

Din cele de mai sus observam ca ]α care este de fapt unghiul ϕ (din coordo-natele sferice, format cu Oz pozitiv) este de 45, deoarece cercul x2+y2 =

(√2)2

este la o înaltime de√

2 unitati. Folosind coordonatele sferice rezulta: x = ρ cos θ sinϕy = ρ sin θ sinϕz = ρ cosϕ

ρ ∈ [0, 2]θ ∈ [0, 2π]

ϕ ∈[0,π

4

] ∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = ρ2 sinϕ

V ol (V ) =

∫∫∫V

dxdydz =

∫∫∫∆

ρ2 sinϕdρdθdϕ =

=

∫∫ ∫[0,2]×[0,2π]×[0, π4 ]

ρ2 sinϕdρdθdϕ =

=

∫ 2

0

ρ2dρ ·∫ 2π

0

dθ ·∫ π

4

0

sinϕdϕ =

=8

3· 2π ·

(−√

2

2+ 1

)=

3

(2−√

2).

b) Masa (M) si centrul de greutate (xG, yG,zG) al unui corp material,

9

Page 10: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

care ocupa domeniul D din R3, de densitate ρ (x, y, z) cunoscuta

M =

∫∫∫D

ρ (x, y, z) dxdydz,

xG =1

M

∫∫∫D

xρ (x, y, z) dxdydz,

yG =1

M

∫∫∫D

yρ (x, y, z) dxdydz,

zG =1

M

∫∫∫D

zρ (x, y, z) dxdydz.

Exemplul 8: Calculati masa corpului marginit de portiunea de elipsoid x2 +y2

9+z2

4= 1

y ≥ 0

stiind ca densitatea lui variaza dupa legea ρ (x, y, z) = 3x+ 1.Solutie:

M =

∫∫∫V

(3x+ 1) dxdydz = 3

∫∫∫V

xdxdydz +

∫∫∫V

dxdydz.

10

Page 11: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Pentru calculul primei integrale vom folosi coordonatele sferice generalizate,iar cea de-a doua reprezinta formula volumului si avem: x = 1 · ρ cos θ sinϕ

y = 3 · ρ sin θ sinϕz = 2 · ρ cosϕ

ρ ∈ [0, 1]θ ∈ [0, π]ϕ ∈ [0, π]

∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = 6ρ2 sinϕ

M = 3

∫∫∫∆

ρ cos θ sinϕ · 6ρ2 sinϕdρdθdϕ+1

2· 4π · 1 · 3 · 2

3=

= 18

∫ 1

0

ρ3dρ ·∫ π

0

cos θdθ︸ ︷︷ ︸0

·∫ π

0

sin2 ϕdϕ+ 4π = 4π.

c) Momentele de inertie ale unui corp material, (al carui model matem-atic îl constituie domeniul compact si marginit D ⊂ R3 de densitate ρ (x, y, z)cunoscuta), fata de axele de coordonate, fata de origine si fata de planele decoordonate sunt:

I0x =

∫∫∫D

(y2 + z2

)ρ (x, y, z) dxdydz

I0y =

∫∫∫D

(x2 + z2

)ρ (x, y, z) dxdydz

I0z =

∫∫∫D

(x2 + y2

)ρ (x, y, z) dxdydz

I0 =

∫∫∫D

(x2 + y2 + z2

)ρ (x, y, z) dxdydz

Ix0y =

∫∫∫D

z2ρ (x, y, z) dxdydz

Iy0z =

∫∫∫D

x2ρ (x, y, z) dxdydz

Iz0x =

∫∫∫D

y2ρ (x, y, z) dxdydz.

Exemplul 9: Calculati momentul de inertie, în raport cu originea, pentruportiunea de sfera

(V ) :

x2 + y2 + z2 ≤ 9z ≥ 0

,

11

Page 12: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

densitatea de masa fiind ρ (x, y, z) = y.Solutie:

I0 =

∫∫∫V

(x2 + y2 + z2

)ρ (x, y, z) dxdydz =

=

∫∫∫V

(x2 + y2 + z2

)ydxdydz.

Pentru a calcula aceasta integrala vom folosi coordonatele sferice x = ρ cos θ sinϕy = ρ sin θ sinϕz = ρ cosϕ

ρ ∈ [0, 3]θ ∈ [0, 2π]

ϕ ∈[0,π

2

] ∣∣∣∣D (x, y, z)

D (ρ, θ, ϕ)

∣∣∣∣ = ρ2 sinϕ

si avem

I0 =

∫∫∫∆

ρ2 · ρ sin θ sinϕ · ρ2 sinϕdρdθdϕ =

=

∫∫ ∫[0,3]×[0,2π]×[0, π2 ]

ρ5 · sin θ · sin2 ϕdρdθdϕ =

=

∫ 3

0

ρ5dρ ·∫ 2π

0

sin θdθ︸ ︷︷ ︸0

·∫ π

2

0

sin2 ϕdϕ = 0.

12

Page 13: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

d) Potentialul newtonian (gravitational) al unui corp D într-un punctP (x0, y0, z0) dat din spatiu

U (x0, y0, z0) =

∫∫∫D

ρ (x, y, z) dxdydz√(x− x0)

2+ (y − y0)

2+ (z − z0)

2.

unde ρ (x, y, z) este densitatea într-un punct curent M (x, y, z) ale corpului D.e) Atractia exercitata de un corp material D asupra unui punct

material P (x0, y0, z0) de masa m

Fx = Km

∫∫∫D

(x− x0) ρ (x, y, z) dxdydz[(x− x0)

2+ (y − y0)

2+ (z − z0)

2] 32

,

Fy = Km

∫∫∫D

(y − y0) ρ (x, y, z) dxdydz[(x− x0)

2+ (y − y0)

2+ (z − z0)

2] 32

,

Fz = Km

∫∫∫D

(z − z0) ρ (x, y, z) dxdydz[(x− x0)

2+ (y − y0)

2+ (z − z0)

2] 32

,

unde K este constanta atractiei universale iar ρ (x, y, z) este densitatea într-unpunct curent M (x, y, z) ale corpului D.

13

Page 14: Curs de Matematici Speciale - Cadastru + IC German…a … · Dac…a Deste paralelipipedul dreptunghic [a;b] [c;d;] [e;f] atunci ZZZ D ... Observ…am c …a domeniul este un paralelipiped

Probleme propuse1. Sa se calculeze I =

∫∫∫V

ydxdydz dacã V este marginit de planele de

coordonate x = 0, y = 0, z = 0 si planul

6x+ 3y + 4z = 12.

2. Calculati ∫∫∫V

e−x2−y2−z2 · zdxdydz

undeV =

(x, y, z) ∈ R3 | x2 + y2 + z2 ≤ 9, z ≥ 0

.

3. Sa se calculeze ∫∫∫V

1√x2 + 9y2 +

(z − 1

2

)2 dxdydz,daca

V =

(x, y, z) ∈ R3 | x2 + 9y2 + z2 ≤ 3

4+ z

.

4. Trecând la coordonate sferice calculati urmatoarea integrala∫∫∫V

√x2 + y2 + z2dxdydz

unde V este domeniul limitat de sfera x2 + y2 + z2 = z.

Lector dr. mat. Hedrea Ioan CiprianDepartamentul de MatematicaUniversitatea Politehnica Timisoara

14