Integrais Múltiplos em Coordenadas Não Cartesianas
Integrais Múltiplos em Coordenadas Não Cartesianas
52 / 100
3.3 Integração em coordenadas não cartesianas
3.3.1 Coordenadas polares e integração dupla
Coordenadas polares
Coordenadas cilı́ndricas
Coordenadas esféricas
53 / 100
Problema
Estudar integração dupla em regiões da forma
y
y
x2 + y 2 ≤ 1
x2 + y 2 = 4
D1
D2
x
x2 + y 2 = 1
x2 + y 2 = z 2
x y
x
Coordenadas polares
55 / 100
Coordenadas cartesianas vs Coordenadas polares
Coordenadas polares
Coordenadas cartesianas
y P = (x, y) P = (r, θ)
O x
I origem do referencial O, um
I origem do referencial O e dois eixo e um ângulo
eixos
I r é a distância a O
I x distância na horizontal a O
I y distância na vertical a O I θ ângulo entre o eixo polar e
a horizontal
56 / 100
Coordenadas cartesianas vs Coordenadas polares
y P = (x, y)
Coordenadas
r Coordenadas polares
cartesianas
P = (r, θ)
P = (x, y)
θ
x
I Da trigonometria do retângulo vem
x = r cos θ r ∈ [0 + ∞[
y = r sen θ θ∈R
I Logo
p
x2 + y 2 = (r cos θ)2 + (r sen θ)2 = r2 =⇒ r = x2 + y 2
e para x 6= 0
y
= tan θ
x
57 / 100
I Para passar de coordenadas polares a cartesianas
x = r cos θ
r ∈ [0, +∞[
θ ∈ [0, 2π[
y = r sen θ
tan θ = y ,
x 6= 0
x
58 / 100
Observação
59 / 100
Exemplo
π
I As coordenadas cartesianas de (r, θ) = 7, são
3
√
π π 7 7 3
(x, y) = (r cos θ, r sen θ) = 7 cos , 7 sen = , .
3 3 2 2
√ π
I As coordenadas polares de (x, y) = (3, 3) são (r, θ) = 3 2, , pois
4
√ √ √
2 2
r = 3 + 3 = 18 = 3 2
θ = arctan 3 = arctan 1 = π ,
para θ ∈ [0, 2π[.
3 4
π
I As coordenadas polares de (x, y) = (0, 2) são (r, θ) = 2, .
2
60 / 100
Exemplo
I A circunferência de equação x2 + y 2 = 4 é descrita em coorde-
nadas polares pelas equações
r = 2, 0 ≤ θ < 2π.
x2 + y 2 = 1
61 / 100
Retângulo polar
Estas regiões são casos particulares de “retângulos polares”
R = {(r, θ) : a ≤ r ≤ b, α ≤ θ ≤ β}.
62 / 100
Integração dupla em coordenadas polares
R = {(r, θ) : a ≤ r ≤ b, α ≤ θ ≤ β}.
63 / 100
I Subdividimos [a, b] em n subintervalos e [α, β] em k subintervalos:
a = r0 < r1 < · · · < rn−1 < rn = b e α = θ0 < θ1 < · · · < θk−1 < θk = β.
b θj+1
I O “centro” de Rij é, em coordenadas polares,
1 1
ri∗ = (ri + ri+1 ), θj∗ = (θi + θi+1 ). θj
2 2
I A área do retângulo polar Rij é
1 2 1 a
∆Aij = ri+1 ∆θj − ri2 ∆θj = ri∗ ∆ri ∆θj , ri+1
2 2 ri
onde ∆θj = θj+1 − θj e ∆rj = rj+1 − rj . β
α
f (x∗i , yj∗ )∆Aij = f (ri∗ cos θj∗ , ri∗ sen θj∗ )ri∗ ∆ri ∆θj .
65 / 100
I Quando n, k −→ ∞ o valor da soma de Riemann de g é o integral
definido de g em R e
ZZ Z βZ b
g(r, θ) dr dθ = g(r, θ) dr dθ.
R α a
66 / 100
Observação
I A integração dupla em coordenadas polares goza das mesmas propriedades
que a integração dupla em coordenadas retangulares.
67 / 100
Exercı́cio
Calcular ZZ
x2 + y 2 dx dy
D
onde
D = {(x, y) ∈ R2 : 1 ≤ x2 + y 2 ≤ 4, x ≥ 0, y ≥ 0}.
y
Descrição de D em coordenadas
cartesianas:
x2 + y 2 = 4 p p
D 0 ≤ x ≤ 1, 1 − x2 ≤ y ≤ 4 − x2
p
x2 + y 2 = 1 1 ≤ x ≤ 2, 0 ≤ y ≤ 4 − x2
x
Resolução.
π
ZZ Z
2
Z 2
2 2
(r cos θ)2 + (r sen θ)2 r dr dθ
x + y dx dy =
D 0 1
68 / 100
Resolução (cont.).
π
ZZ Z
2
Z 2
x2 + y 2 dx dy = (r cos θ)2 + (r sen θ)2 r dr dθ
D 0 1
π
Z
2
Z 2
r2 cos2 θ + sen2 θ r dr dθ
=
0 1
π π r=2
2
r4
Z Z Z
2 2
3
= r dr dθ = dθ
0 1 0 4 r=1
Z π Z π
2 1 2 15
= 4− dθ = dθ
0 4 0 4
π
15 r= 2
15π
= θ = .
4 θ=0 8
x = r cos θ, y = r sen θ
e temos sempre x2 + y 2 = r2 .
69 / 100
Exercı́cio
Determinar o volume do sólido limitado pelo plano z = 0 e o pelo
parabolóide z = 1 − x2 − y 2 .
Resolução.
Primeiro observe-se que a projeção do parabolóide com o plano
z = 0 é a cı́rculo D de equação x2 + y 2 ≤ 1.
O volume do sólido é dado por
ZZ ZZ
2 2
1 − (x2 + y 2 ) dx dy
V = 1 − x − y dA =
D D
70 / 100
Resolução (cont.).
Usando coordenadas polares para o cálculo do integral, vem
ZZ
1 − (x2 + y 2 ) dx dy
V =
D
Z 2π Z 1
1 − r2 r dr dθ
=
0 0
Z 2π Z 1
r − r3 dr dθ
=
0 0
2π
r2 r4 r=1
Z
= − dθ
0 2 4 r=0
Z 2π θ=2π
1 1 π
= dθ = θ = .
0 4 4 θ=0 2
Se tivéssemos usado coordenadas retangulares em vez de polares,
terı́amos obtido o integral
ZZ Z 1 Z √1−x2
2 2
1 − x2 − y 2 dy dx
1 − x − y dA = √
D −1 − 1−x2
que não é fácil de calcular. 71 / 100
Exercı́cioZ Z
Calcular (3x + y) dA onde
D
D = {(x, y) ∈ R2 : 1 ≤ x2 + y ≤ 4, y ≥ 0}.
Resolução.
ZZ Z πZ 2
(3x + y) dx dy = 3r cos θ + r sen θ r dr dθ
D 0 1
Z πZ 2
3r2 cos θ + r2 sen θ dr dθ
=
0 1
π r=2
r3
Z
3
= r cos θ + sen θ dθ
0 3 r=1
Z π
7
= 7 cos θ + sen θ dθ
0 3
θ=π
7 14
= 7 sen θ + (− cos θ) = .
3 θ=0 3
72 / 100
3.3.2 Coordenadas cilı́ndricas e integração tripla
Coordenadas cilı́ndricas
73 / 100
Coordenadas cartesianas vs Coordenadas cilı́ndricas
z
P = (x, y, z)
Coordenadas
Coordenadas cilı́ndricas
cartesianas
P = (r, θ, z)
P = (x, y, z)
y
r
θ
x
(x, y, 0)
I Para passar de coordenadas cilı́ndricas a cartesianas
x = r cos θ r ∈ [0, +∞[
y = r sen θ θ ∈ [0, 2π[
z=z z ∈ R.
74 / 100
Observação
75 / 100
Exemplo
π
I As coordenadas cartesianas de (r, θ, z) = 7, , 5 são
√ ! 3
7 7 3
(x, y, z) = , ,5 .
2 2
I As coordenadas
√ πcilı́ndricas
de (x, y, z) = (3, 3, 1) são
(r, θ, z) = 3 2, , 1 .
4
0 ≤ r ≤ 2, 0 ≤ θ < 2π, 0 ≤ z ≤ 3.
76 / 100
Exemplo
O cone (metade) definido por
n p o
E = (x, y, z) ∈ R3 : x2 + y 2 ≤ z ≤ 2
p
(limitado inferiormente pela superfı́cie cónica z = x2 + y 2 e superiormente
pelo plano z = 2) é descrito em coordenadas cilı́ndricas da forma
0 ≤ r ≤ 2, 0 ≤ θ < 2π, r ≤ z ≤ 2.
p
Note-se que z = x2 + y 2 =⇒ z 2 = x2 + y 2 .
Para z = 2 temos a circuferência x2 + y 2 = 4 e a projeção do cone no plano xy
é o cı́rculo x2 + y 2 ≤ 4.
p √
Note-se também que podemos escrever z ≥ x2 + y 2 = r2 = r.
77 / 100
Integração tripla em coordenadas cilı́ndricas
Suponhamos que se pretende calcular o integral da função f contı́nua em U ⊂ R3
(dado em coordenadas retangulares)
ZZZ
f (x, y, z) dV,
U
onde U é uma região da forma
z z
x2 + y 2 = a 2
x2 + y 2 = z 2
y y
x x
Em coordenadas cilı́ndricas
U1 = {(r, θ, z) : 0 ≤ r ≤ a, 0 ≤ θ < 2π, 0 ≤ z ≤ h}
U2 = {(r, θ, z) : 0 ≤ r ≤ h, 0 ≤ θ < 2π, r ≤ z ≤ h}
78 / 100
I Consideremos o caso geral de f estar definida em U , uma
região do tipo I de R3
z U
z = γ2 (x, y)
z = γ1 (x, y)
x D
79 / 100
I Da teoria de integração sobre uma região do tipo I de R3
sabe-se que
ZZZ Z Z "Z #
γ2 (x,y)
f (x, y, z) dV = f (x, y, z) dz dA
U D γ1 (x,y)
80 / 100
[Mudança para coordenadas cilı́ndricas num integral triplo]
Se f é contı́nua numa região U ⊂ R3 que pode ser descrita em
coordenadas cilı́ndricas por
α≤ θ ≤β
µ1 (θ) ≤ r ≤ µ2 (θ)
então
ZZZ
f (x, y, z) dV =
U
Z β Z µ2 (θ) Z γ2 (r cos θ,r sen θ)
= f (r cos θ, r sen θ, z) r dz dr dθ.
α µ1 (θ) γ1 (r cos θ,r sen θ)
81 / 100
Exemplo
Calcular
ZZZ
(x2 + y 2 ) dV z
U x2 + y 2 = 1
sendo 3
n o
U = (x, y, z) ∈ R3 : x2 + y 2 ≤ 1, −2 ≤ z ≤ 3 .
y
x
-2
82 / 100
Resolução.
Fazendo a mudança para coordenadas cilı́ndricas,
x = r cos θ, y = r sen θ, z = z,
temos
ZZZ Z 2π Z 1Z 3
2 2
(r cos θ)2 + (r sen θ)2 r dz dr dθ
(x + y ) dV =
U 0 0 −2
Z 2π Z 1Z 3
= r2 r dz dr dθ
0 0 −2
Z 2π Z 1Z 3
= r3 dz dr dθ
0 0 −2
Z 2π Z 1 z=3
3
= r z dr dθ
0 0 z=−2
Z 2π Z 1 Z 2π r=1
3 5 4
= 5r dr dθ = r dθ
0 0 0 4 r=0
Z 2π θ=2π
5 5 5
= dθ = θ = π.
0 4 4 θ=0 2
83 / 100
3.3.3 Coordenadas esféricas e integração tripla
Coordenadas esféricas
84 / 100
Coordenadas cartesianas vs Coordenadas esféricas
P = (x, y, z) P = (ρ, θ, φ)
y
−−→
r ρ = kOP k
θ
x (x, y, 0)
85 / 100
I Para passar de coordenadas esféricas a cartesianas
x = ρ sen φ cos θ ρ ∈ [0, +∞[
y = ρ sen φ sen θ , θ ∈ [0, 2π[
z = ρ cos φ φ ∈ [0, π]
86 / 100
Observação
87 / 100
Exemplo
1. As coordenadas esféricas do ponto (x, y, z) = (4, 0, 0) são
π
(ρ, θ, φ) = 4, 0, , obtidas da forma
2
p √
ρ = x2 + y 2 + z 2 = 42 + 0 + 0 = 4
y 0
tan θ = = = 0 =⇒ θ = 0
x 4
z 0 π
cos φ = = = 0 =⇒ φ =
ρ 4 2
π π
2. As coordenadas cartesianas do ponto (ρ, θ, φ) = 4, , são
√ √ √ 6 4
(x, y, z) = ( 6, 2, 2 2), obtidas fazendo
√ √
√
π π 2 3
x = ρ sen φ cos θ = 4 sen cos = 4 · · = 6
4 6 √2 2
π π 2 1
y = ρ sen φ sen θ = 4 sen sen = 4 · ·
4√ 6 2 2
√
z = ρ cos φ = 4 cos π = 4 2 = 2 2
4 2
3. As coordenadas
√ π esféricas do ponto (x, y, z) = (0, 1, 1) são
π
(ρ, θ, φ) = 2, , .
2 4 88 / 100
Exemplo
1. Uma superfı́cie esférica de equação x2 + y 2 + z 2 = a2 , com
a > 0, em coordenadas esféricas descreve-se simplesmente por
ρ = a, θ ∈ [0, 2π[ , φ ∈ [0, π],
e a esfera de inequação x2 + y 2 + z 2 ≤ a2 por
ρ ≤ a, θ ∈ [0, 2π[ , φ ∈ [0, π].
z
y
y
x
x
E = {(ρ, θ, φ) : a ≤ ρ ≤ b, α ≤ θ ≤ β, c ≤ φ ≤ d},
onde a ≥ 0, β − α ≤ 2π e d − c ≤ π.
1
Uma cunha esférica é a um “paralelipı́pedo” em coordenadas esféricas
90 / 100
Seja P = {(ρ, θ, φ) : a ≤ ρ ≤ b, α ≤ θ ≤ β, c ≤ φ ≤ d} e f : P ⊂ R3 −→ R.
I A uma subdivisão do intervalo [a, b] em n subintervalos, de [α, β] em m
subintervalos e de [c, d] em l subintervalos, associamos uma subdivisão da
cunha esférica E em n × m × l cunhas esféricas
Eijk = [ρi , ρi+1 ] × [θj , θj+1 ] × [φk , φk+1 ].
I Uma aproximação para o volume de Eijk
é dada por
ρ2i sen φ
∆Vijk = e ek ∆ρi ∆θj ∆φk
onde (e
ρi , e
θj , e
ρk ) é um ponto em Eijk .
I Para este ponto considerem-se as suas coordenadas retangulares
(x∗i , yj∗ , zk∗ ) = (e
ρi sen φ
ej cos e
θk , e
ρi sen φ
ej sen e
θk , e
ρi cos φ
ek ).
91 / 100
I Consideremos a soma
n−1
X m−1
XXl−1
ρi sen φ
f (e ej cos e
θk , e
ρi sen φ
ej sen e
θk , e
ρi cos φ ρ2i sen φ
ek ) e ek ∆ρi ∆θj ∆φk
i=0 j=0 k=0
a ≤ ρ ≤ b, α ≤ θ ≤ β, c≤φ≤d
93 / 100
Observação
94 / 100
Exemplo
Calcular ZZZ
2 +y 2 +z 2 )3/2
e(x dV
U
onde U é a esfera unitária,
U = (x, y, z) ∈ R3 : x2 + y 2 + z 2 ≤ 1 .
Resolução.
Vamos usar coordenadas esféricas:
U = {(ρ, θ, φ) : 0 ≤ ρ ≤ 1, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π} .
95 / 100
Resolução (cont.).
ZZZ
3/2
Z π Z 2π Z 1
(x2 +y 2 +z 2 ) 2 )3/2
e dV = e(ρ ρ2 sen φ dρ dθ dφ
U 0 0 0
Z π Z 2π Z 1
3
= ρ2 eρ sen φ dρ dθ dφ
0 0 0
π 2π 3 ρ=1
eρ
Z Z
= sen φ dθ dφ
0 0 3 ρ=0
Z π Z 2π
e−1
= 3 sen φ dθ dφ
0 0
Z π θ=2π
e−1
= 3 sen φ θ dφ
Z0 π θ=0
2π(e−1)
= 3 sen φ dφ
0
φ=π
2π(e−1) 4π(e−1)
= 3 − cos φ = 3
φ=0
96 / 100
Resolução (cont.).
Observe-se que teria sido muito difı́cil avaliar o integral deste
exemplo sem coordenadas esféricas. Em coordenadas retangulares
o integral teria ficado
Z Z √
1 1−x2 Z √ 1−x2 −y 2
2 +y 2 +z 2 )3/2
√ √ e(x dz dy dx.
−1 − 1−x2 − 1−x2 −y 2
97 / 100
Exemplo
Determinar o volume de um cone de gelado limitado superiormente
pela parte superior da superfı́cie esférica 2 2 2
p x + y + z = 2 e inferi-
ormente pela superfı́cie cónica z = x2 + y 2 .
98 / 100
Resolução.
Interseção das superfı́cies:
2
x +p y2 + z2 = 2
=⇒
z = x2 + y 2
2
x + y2 + z2 = 2
=⇒ =⇒
z 2 = x2 + y 2
2
x + y2 = 1
=⇒
z=1