PRÁ CTI CA COO RDENADAS C ILÍNDRI C AS Y ESFÉRI CAS C URS O 2016 - 2017
C ÁL CUL O I I
Prácticas Matla b
Práctica 3 (21/02/2017)
Objetivos
o Estudiar los sistemas de coordenadas cilíndricas y esféricas.
o Definir regiones de
de integrales ℝ3 en coordenadas cilíndricas y esféricas, para el planteamiento
triples.
o Calcular integrales triples en coordenadas cilíndricas y esféricas.
Coordenadas cilíndric as
Las coordenadas cilíndricas se obtienen utilizando coordenadas polares en uno de los planos
coordenados y manteniendo la tercera variable.
Paso de cilíndricas a cartesianas:
x = r cos
y = rsen
z = z
Paso de cartesianas a cilíndricas:
y
r x2 y2 arctg ; zz
x
;
( , signo signo y)
Las superficies de ecuación más sencilla en coordenadas cilíndricas son:
• r = a . Es el cilindro de eje 0Z y radio a ;
• = b . Es el semiplano que contiene al eje 0Z y forma ángulo b con el plano XZ
( x > 0 );
• z = c . Es un plano perpendicular al eje 0Z .
Llamaremos caja cilíndrica a la región del espacio acotada por las superficies anteriores.
PÁGINA MATLAB: COORDENADAS CILÍNDRICAS Y
Coordenadas esféric as
Paso de esféricas a cartesianas:
x = sen cos
y = sensen
z = cos
Paso de cartesianas a esféricas:
y
x2 y2 z2 arctg ;
; x
x2 y2
arctg
z
Las superficies de ecuación más sencilla en coordenadas esféricas son:
• = a . Es la esfera de centro el origen y radio a ;
• = b . Es el semiplano que contiene al eje 0Z y forma ángulo b con el plano XZ ,
( x > 0 );
• = c . Es un semicono de eje 0Z .
Llamaremos caja esférica a la región del espacio acotada por las superficies anteriores.
Herramienta para visualizar cajas cilíndric as y esféric as
Vete a la página
http://personales.unican.es/alvareze/CalculoWeb/CalculoII/integracion_multiple.html
y pincha en el enlace Coordenadas cilíndricas y esféricas
MATLAB: PRÁCTICA PÁGINA
Se abrirá una herramienta que mostrará, en la ventana de la derecha, los valores de las
coordenadas de un punto del espacio y, en el la ventana de la izquierda, sólidos determinados
por variación de las coordenadas en un cierto intervalo.
Ejercicios
Si H es el sólido limitado inferiormente por z =x2 y2 y superiormente por z = 3 , calcula:
El volumen de H
1
El valor de la integral I =
H (x y 2z)dV
El planteamiento se hará a mano y las integrales se resolverán con Matlab.
Resolución
Este es el ejercicio propuesto 14 del tema 1.
Se recomienda utilizar coordenadas cilíndricas ya que la proyección del sólido en z 0 es un
disco, y los límites de integración de la integral doble son constantes en este sistema de
coordenadas.
La representación del sólido se muestra en la figura
a) La ecuación del cono en coordenadas cilíndricas es
z = x2 y2 zr
b) La región cónica de la figura, expresada en coordenadas cilíndricas es:
H ( , r, z) / 0 2 , 0 r 3, r z 3
Trasladamos estos límites a la integral para calcular el volumen:
PÁGINA MATLAB: COORDENADAS CILÍNDRICAS Y
2 3 3
V rdzdrd
0 0 r
Resolvemos con Matlab:
syms t r z
V=int(int(int(r,z,r,3),r,0,3),t,0,2*pi)
%V=9pi
c) Planteamos la integral propuesta
2 3 3
I = (x y 2z)dV
H r cos r sen 2z rdzdrdt
0 0 r
Resolvemos con Matlab:
syms t r z
I=int(int(int((r*cos(t)+r*sin(t)-2*z)*r,z,r,3),r,0,3),t,0,2*pi)
%I=-81pi/2
Halla mediante una integral triple el volumen del sólido del primer
octante limitado por los planos coordenados, por el cono z2 = x2 y2 y por la esfera x2 y2
2 Calcula la densidad media de dicho sólido, sabiendo que su densidad de masa viene dada por la
z
x, y, z =
x2 y2
El planteamiento se hará a mano y las integrales se resolverán con Matlab.
Resolución
Este es el ejercicio propuesto 15 del tema 1.
Se recomienda utilizar coordenadas esféricas ya que el sólido es una caja esférica, y los límites
de integración son constantes en este sistema de coordenadas.
La representación del sólido se muestra en la figura
MATLAB: PRÁCTICA PÁGINA
Escribimos las ecuaciones del cono y de la esfera en coordenadas esféricas:
2 2 2 2 2 2 2
o Ecuación del cono: z = x y cos sen
4
o Ecuación de la esfera: 22
2
x y z =8
2 2 2
8
El sólido de la figura utilizando coordenadas esféricas es:
H , , / 0 , 0 , 02 2
2 4
La integral para obtener el volumen del sólido utilizando coordenadas esféricas es:
/2 2 2 /4
volumen H 2sen d d
0 0 0
d
Resolvemos con Matlab:
syms rho t phi
volumen=int(int(int(rho^2*sin(phi),phi,0,pi/4),rho,0,2*sqrt(2)),t,0,pi/2))
%volumen =(8*pi*(2^(1/2) - 1))/3
Nota: este volumen también podría calcularse utilizando coordenadas cilíndricas:
/2 2
8r 2
volumen H
0 0 r
r dz dr d
a) La expresión para calcular la densidad media es:
PÁGINA MATLAB: COORDENADAS CILÍNDRICAS Y
x, y, z dV
densidad media H
volumen H
La función densidad en coordenadas esféricas es:
z cos cos
x, y, z = x, y, z
x2 y2 sen sen
Y la
integral /2 2 2 /
4
2 cos d d
masa de H
0 0 0
d
Escribimos el código para calcular la densidad media:
masa=int(int(int(cos(phi)*rho^2,phi,0,pi/4),rho,0,2*sqrt(2)),t,0,pi/2)
densidadMedia=masa/volumen
%masa =(8*pi)/3
%densidadMedia =1/(2^(1/2) - 1)
Calcula mediante una integral triple el volumen del sólido H, limitado superiormente por el paraboloide z = 9
Encuentra la temperatura media de H, sabiendo que la temperatura en cada punto viene dada por la distancia del p
3
El planteamiento se hará a mano y las integrales se resolverán con Matlab.
Resolución
Este es el ejercicio propuesto 16 del tema 1.
MATLAB: PRÁCTICA PÁGINA
La representación del sólido se muestra en la figura
Se recomienda utilizar coordenadas cilíndricas ya que la proyección del sólido en z 0 es un
disco, y los límites de integración de la integral doble son constantes en este sistema de
coordenadas.
a) Escribimos las ecuaciones del paraboloide y el cilindro en coordenadas cilíndricas:
o Ecuación del
z = 9 x2 y2 z9r2
paraboloide:
o Ecuación del x2 y2 =1 r 1
cilindro:
La región de la figura, expresada en coordenadas cilíndricas es:
H ( , r, z) / 0 2 , 1 r 3, 0 z 9 r 2
La integral triple para calcular el volumen de H es:
2 3 9r
2
V rdzdrd
0 1 0
Resolvemos con Matlab:
syms t r z
V=int(int(int(r,z,0,9-r^2),r,1,3),t,0,2*pi)
%V=32pi
b) La función temperatura es
T (x, y, z) x2 y2 T (r, , z) r
Y la temperatura media
T x, y, z dV
temperatura media H
volumen H
Escribimos la integral de la temperatura en coordenadas cilíndricas
2 3 9r
2
I r 2dzdrdt
0 1 0
Escribimos el código para calcular la temperatura media:
I=int(int(int(r^2,z,0,9-r^2),r,1,3),t,0,2*pi)
Tm=I/V
%I =(296*pi)/5
%Tm =37/20
Podemos calcular también los puntos del sólido que tienen una temperatura igual a la
temperatura media:
PÁGINA MATLAB: COORDENADAS CILÍNDRICAS Y
T (r, , z) r 37
r 20
Es decir los puntos buscados pertenecen al cilindro de radio 37 / 20
El código en Matlab para dibujar el sólido, junto con los puntos que están a temperatura
media es:
%dibujo del paraboloide
t=linspace(0,2*pi,20);
r=linspace(1,3,20);
[T,R]=meshgrid(t,r);
X=R.*cos(T);Y=R.*sin(T);Z=9-R.^2;
surf(X,Y,Z); alpha 0.5;
%dibujo del cilindro interior
zv=[0,8];
[U V]=meshgrid(t,zv);
X1=cos(U);Y1=sin(U);Z1=V;
hold on
surf(X1,Y1,Z1)
%dibujo del cilindro cuyos puntos están a temperatura media
radio=37/20;
zv1=[0,9-radio^2];
[U1 V1]=meshgrid(t,zv1);
X2=radio*cos(U1);Y2=radio*sin(U1);Z2=V1;
surf(X2,Y2,Z2);
axis equal
hold off
Dibuja, utilizando Matlab la porción del cilindro x2 y2 = 2x entre z = 3 y z =3 y sobre la misma figura u
4
Calcula, utilizando el paquete simbólico de Matlab, el volumen del sólido H limitado entre las dos hojas de z 2 = 2 x
Calcula con Matlab, utilizando coordenadas cilíndricas, el volumen del sólido limitado por x2 y2 = x entre las d
Resolución
Este es el ejercicio propuesto 12 del tema 1.
a) La representación del sólido se muestra en la figura de la derecha
MATLAB: PRÁCTICA PÁGINA
Código para representar la figura de la izquierda con Matlab:
tv=linspace(0,2*pi,30);
zv=[-3 3];
[T,V]=meshgrid(tv,zv);
X=1+cos(T);Y=sin(T);Z=V;
surf(X,Y,Z)
axis equal
alpha 0.5
hold on
y=linspace(-1,1,20);
z=linspace(-2,2,30);
[Y1,Z1]=meshgrid(y,z);
surf(Z1.^2/2,Y1,Z1)
hold off
b) Expresamos el sólido H utilizando coordenadas cartesianas
H x, y, z / 0 x 2 , 2x x2
y 2x x2 , 2x z 2x
Teniendo en cuenta las simetrías respecto a y=0 y z=0, se sabe que el volumen de H es
4 veces el volumen del primer octante
2
2 x x2 2x 2 2 x x2
volumen H 4 dzdydx 2x dydx
0
0 0 0 0
Escribimos el código Matlab para calcular el volumen del sólido en coordenadas
cartesianas
syms x y z
V=4*int(int(int(1,z,0,sqrt(2*x)),y,0,sqrt(2*x-x^2)),x,0,2)
% V=128/15
La ecuación del disco
x2 y2 = 2x en polares es r 2 , por lo tanto el sólido
cos
en coordenadas cilíndricas es:
PÁGINA MATLAB: COORDENADAS CILÍNDRICAS Y
H ( , r, z)
, 0r2 z 2r cos
/
cos2 , 2 2r cos
Planteamos el volumen en este sistema de coordenadas, aplicando simetrías,
/2 2cos 2r cos / 2 2cos
volumen H 4 2r cos dr d
r dz dr d 4
0
0 0 0 0
El código para resolver la integral es:
syms r t z
V=4*int(int(int(r,z,0,sqrt(2*r*cos(t))),r,0,2*cos(t)),t,0,pi/2)
% V=128/15
c) La representación del sólido z 0 , se muestra en la figura de la derecha
para
La ecuación del disco x2 y2 = x en polares es r cos , por lo tanto el sólido en
coordenadas cilíndricas es:
H ( , r, z)
, 0 r cos , z 1 r 2
/
2 2 1 r2
Y el volumen, en este sistema de coordenadas, aplicando simetrías,
/2 cos 1r2 /2 2cos
volumen H 4 1 r2 dr d
r dz dr d 4
0
0 0 0 0
El código para resolver la integral es:
syms r t z
V=double(4*int(int(int(r,z,0,sqrt(1+r^2)),r,0, cos(t)),t,0,pi/2))
% V= 1.8338
