Capítulo 6

Números complejos en Matlab

Msc. Ing. Jairo E. Márquez D.

6.1 Números complejos

Los números complejos ℂ son una extensión de los números reales ℝ y forman
el mínimo cuerpo algebraicamente cerrado que los contiene, es decir, ℝ ⊂ ℂ. Los
números complejos incluyen todas las raíces de los polinomios, a diferencia de los
números reales.

Todo número complejo se puede representar como la suma de un número real Re

y un número imaginario Im (𝑍 = 𝑥 + 𝑦𝑖), o en forma polar, donde “x” se asocia a un
número real, mientras “yi” a un número imaginario, tal como se ilustra en la figura.

Los números complejos son la herramienta

matemática empleada ampliamente en áreas del
álgebra, análisis, variable compleja, ecuaciones
diferenciales y cálculo integral entre otros. A nivel
de la ingeniería se puede encontrar en la
aerodinámica, hidrodinámica, electromagnetismo,
electrónica, telecomunicaciones, física cuántica,
etc.

Matlab dispone de la función complex, que crea

un número complejo a partir de dos argumentos,
la parte real e imaginaria. Así:
complex(7,3)
ans =
1.00 + 3.0000i

complex(10,-9)
ans =
10.0000 - 9.0000i

Se puede trabajar con complejos en su forma polar y exponencial respectivamente,

𝒁 = 𝝆(𝒄𝒐𝒔𝜽 + 𝒊𝒔𝒆𝒏𝜽)
𝒁𝒆 = 𝝆𝒆𝒊𝜽
Donde 𝜌 = |𝑍|

En términos generales en Matlab se procede a definir la magnitud y el ángulo, para

luego, declarar la ecuación Z, así.

x = a + bi
magnitud = abs(x)
angulo = angle(x)

z = magnitud*cos(angulo) + i*magnitud*sen(angulo)
z= magnitud*exp(i*angulo)

6.2 Funciones de números complejos

Instrucciones Descripción
real(Z) Parte real de un dato complejo.
imag(Z) Parte imaginaria de un dato complejo.
complex(x,y) Complejo construido a partir de las partes
real (x) e imaginaria (yi).
abs(Z) Módulo de un dato complejo, valor absoluto
de un dato real.
conj(Z) Complejo conjugado de un dato complejo
angle(Z) Argumento o ángulo de un dato complejo.
Ángulo de fase.
sign(Z) Función signo si z es complejo y no-nulo,
z/|z|.
isreal(Z) Determina si un dato es real o complejo.
unwrap(x) Ángulo de fase separado.
 cplxpair(x) Clasifica los números en pares complejos
conjugados.
[theta,rho]=cart2pol(x,y) Pasa de coordenadas cartesianas a polares.
Para el caso tridimensional, las
componentes x e y, corresponden a los
[theta,rho,z]=cart2pol(x,y,z)
elementos vectoriales y z indica la
coordenada cilíndrica.
[X,Y]= pol2cart(theta,rho) Pasa de coordenadas polares a cartesianas.
Transforma coordenadas cilíndricas a
cartesianas.
[X,Y,Z]=
pol2cart(theta,rho,Z)

Ejemplos

1. Sea Z=3+4i, hallar rho, theta, Z, Ze:

Z=3+4i
rho=abs(Z), theta=angle(Z);
Z=rho*(cos(theta)+i*sin(theta))
Ze=rho*exp(i*theta);

Z =
3.0000 + 4.0000i

rho =
5

Ze =
3.0000 + 4.0000i

2. Sea Z= 5-8i determinar la tabla de funciones de números complejos:

Z= 5-8i;
real(Z)

ans =
5

imag(Z)
ans =
-8

abs(Z)
ans =
9.4340
conj(Z) % o simplemente Z’
ans =
5.0000 + 8.0000i

angle(Z)
ans =
-1.0122

sign(Z)
ans =
0.5300 - 0.8480i

isreal(Z)
ans =
0

Ángulo de z en radianes.

theta = atan2(imag(Z),real(Z))
theta =
-1.0122

Para obtener el ángulo en grados:

Z= 5-8i;
theta = atan2d(imag(z),real(z))
theta =
36.8699

3. Sea z1=4-6i y z2=2+9i, determinar z =z1 + z2, z = z1 - z2, z = z1*z2 y z = z1/z2

z1=4-6i;
z2=2+9i;
z=z1+z2
z =
6.0000 + 3.0000i

z=z1-z2
z =
2.0000 -15.0000i

z=z1*z2
z =
62.0000 +24.0000i

z=z1/z2
z =
-0.5412 - 0.5647i
6.3 Gráfica de números complejos

Para graficar un número complejo se emplean ejes de coordenadas

rectangulares, partiendo del hecho que los números ±a±bi, se representan por un
punto cuyas coordenadas son (a,b). Tome en cuenta, que se debe marcar con algún
símbolo los puntos de las coordenadas para que estas sean visibles en la gráfica.
Por ejemplo, cada punto está marcado por un círculo negro ('ko')

En términos generales la sintaxis es:

plot(a1,b1,'ko',a2,b2,'ko',...,an,bn,'ko')

Tome en cuenta que la instrucción dada, aplica para una gráfica con números
reales.

Ejemplos

1. Graficar los números (4 + 2i), (-3 + 5i), (-5-1i), (-2+i) y (1+4i).

plot(4,2,'ko',-3,5,'ko',-5,-1,'ko',-2,1,'ko',1,4,'ko')

Figura. Distribución de parejas ordenadas en un plano cartesiano complejo

Se puede definir los límites de la gráfica para observar mejor los puntos que
representan las coordenadas, empleando la estructura [xmin xmax ymin ymax]

De igual manera se colocan las etiquetas para explicar mejor la distribución de los
puntos en el plano imaginario y real.
plot (4,2,'ko',-3,5,'ko',-5,-1,'ko',-2,1,'ko',1,4,'ko') % k indica el color negro
xlabel ('x (eje real)'),ylabel ('y (eje imaginario)')
title('Gráfica compleja')
grid on
axis([-6 5 -2 6])

Figura. Distribución de parejas ordenadas en un plano cartesiano complejo, al cual se le ha colocado

las etiquetas o nombres a los planos x e y, al igual que una grilla para facilitar la ubicación de las
parejas. Nótese que se han ajustado los ejes por medio de la instrucción axis, con el fin de mejorar
la observación de las parejas ordenadas en el plano cartesiano.

2. Gráfica de funciones en un plano complejo. Para este caso en particular se

procede a definir el intervalo y emplear las funciones real e imag para
representar los datos. Por ejemplo. Sean las funciones

−𝑥2 −𝑥2
𝑧1 = 𝑥 𝑒 𝑧2 = 2𝑥 𝑒

Graficar en el plano complejo en el intervalo [-3,3] con particiones de 0.2.

Personalice las curvas y etiquete la gráfica con su respectiva grilla.

x=-3:0.2:3;
z1 = x.^exp(-x.^2);
z2 = 2*x.^exp(-x.^2);
real_z1 = real(z1);
imag_z1 = imag(z1);
real_z2 = real(z2);
imag_z2 = imag(z2);
plot(real_z1,imag_z1,'r*',real_z2,imag_z2,'bo')
xlabel('R'),ylabel('I'), grid
title('plano complejo')
 Figura. Gráfica de funciones en un plano cartesiano complejo.

Practica 1.

Dada la tabla, organice los datos en un plano cartesiano real.

Para la segunda tabla crear parejas ordenadas entre valor propio y %de la varianza;
valor propio acumulativo y % acumulativo.
Para la tercera tabla crear parejas ordenadas entre las variables tamaño máx (mm)
y Volumen agregado con cada columna.

Para la cuarta tabla crear parejas ordenadas entre las variables retiros anuales de
agua y agricultura, retiros anuales de agua e industria, retiros anuales de agua y
doméstico.
Control de ejes coordenados

El definir los ejes en una gráfica es importante cuando se desea resaltar detalles de
la misma, o simplemente suprimirlos para observar mejor la curva de la función. La
instrucción estándar es axis, que ajusta las escalas de los ejes x e y según el
intervalo definido.

Las expresiones de la instrucción axis son las siguientes:

- axis: Muestra como salida un vector que define los valores:

axis([xmin, xmax, ymin, ymax])

- axis(axis): Mantiene los ejes de la gráfica en sus valores actuales, donde se

espera la adición de otras gráficas por medio del comando hold on.

- axis('ij'): Emplea el concepto de ejes de pantalla, en la que se especifica el

origen de la gráfica en la esquina superior izquierda, y el eje j describe la
dirección vertical descendente.

- axis('xy'): Emplea los ejes cartesianos estándar, donde el origen queda

definido en la esquina inferior izquierda, y el eje y vertical ascendente.

- axis('equal'): Permite crear un escalado equivalente en ambos ejes.

 - axis('square'): Muestra la ventana de la gráfica en forma cuadrada.

- axis('image'): Redimensiona la ventana, de tal manera que se tiene las

proporciones de la imagen manteniendo el escalado de los ejes.

- axis(’normal’): Elimina las restricciones establecidas por las instrucciones

’equal’ y ’square’.

- axis(’off’): Elimina etiquetas, números y ejes.

- axis('on'): Restablece etiquetas, números y ejes.

Ejemplo

Sea la función 𝑦 = 𝑥𝑠𝑒𝑛(2𝑥)cos(𝑥) graficar en el intervalo [-π,4π], con particiones

de 0.1. Emplear las diversas variantes de axis:

x=-pi:0.1:4*pi;
y=x.*sin(2*x).*cos(x);
plot(x,y,’k’), title('y=x.*sin(2*x).*cos(x)')
xlabel('x'),ylabel('y')

ans =
-4 14 -10 8

Figura 6.5
Figura 6.6

axis axis('ij')

Figura 6.7 Figura 6.8

 axis('xy') axis('equal')

Figura 6.9 Figura 6.10

axis('square') axis('image')

axis('off')

Dado el siguiente código, coloque etiquetas y aplique los diferentes tipos de control
de ejes.
t=-6*pi:0.1:6*pi;
x=2*cos(2*t)-pi/2;
y=2*sin(2*t)+pi/2;;
z=x+i*y;
A=1./z.^2;
a1=real(A); b1=imag(A); plot(a1,b1)

6.4 Transformación de coordenadas cartesianas a esféricas

El sistema de coordenadas esféricas se emplea para calcular la posición de un

punto mediante una distancia y dos ángulos. La sintaxis en Matlab es:

[azimuth,elevation,r] = cart2sph(X,Y,Z)
Esta instrucción transforma las coordenadas cartesianas almacenadas en los
elementos de las matrices X, Y y Z en coordenadas esféricas. El azimut y la
elevación son desplazamientos angulares expresados en radianes.

El azimut es el ángulo en sentido antihorario en el plano xy medido desde el eje x

positivo. Es decir, es el ángulo que forma un cuerpo celeste con el Norte, cuya
medida se realiza en sentido de rotación de las agujas de un reloj alrededor del
horizonte del observador. Por ejemplo, un cuerpo celeste que se encuentra al Norte,
su azimut es 0º, al Este es 90º, al Sur 180º y al Oeste 270º.

El ángulo de elevación, se traza desde el plano xy, y r es la distancia desde el

origen a un punto P, tal como de observa en la figura:

Tenga en cuenta que las Matrices X, Y, y Z deben ser del mismo tamaño (o
cualquiera de ellos puede ser escalar).

Ejemplos:

1. Sea z = 5 - 8i, determinar los ángulos de elevación y azimutal (en grados), al

igual que la distancia con respecto a la coordenada z. Esta transformación
consiste en pasar de coordenadas cartesianas a esféricas.

x=5; y=-8; z=x+y

[azimuth,elevation,r] = cart2sph(x,y,z)
azimuth =
-1.0122

elevation =
-0.3079

r =
9.8995

Para transformar el azimuth y elevación en grados, multiplicar por 180/π:

azimuth = -57.9948

elevation = -17.6414

2. Dado el número complejo A= 2 + 3i, graficar el radiovector.

A=2+3i;
compass(A); % genera un radiovector

3. Para graficar la pareja ordenada (2,3) que conforma el radiovector, se puede

expresar de la siguiente manera:

A=2+3i;
polar(angle(A), abs(A),'*k'); % k representa el color negro
 4. Dada la ecuación r(Ɵ), graficar en el intervalo [0,2π] con particiones de 0.1.

𝑟 = 2cos(2𝜃)

th=0:0.01:2*pi;
r=2*cos(2*th);
polar(th,r)
title('r=2cos(2th)')

Ejercicios resueltos

1. Sea la ecuación 𝑧 2 + 2𝑧 + 3. Determinar las raíces según la ecuación

cuadrática estándar.

a=1; b=2; c=3;

z1=(-b+sqrt(b^2-4*a*c))/(2*a)
z1 =
-1.0000 + 1.4142i

z2=(-b-sqrt(b^2-4*a*c))/(2*a)
z2 =
-1.0000 - 1.4142i

a*z1^2+b*z1+c % prueba de la raíz x1

ans =
-4.4409e-16

a*z2^2+b*z2+c % prueba de la raíz x2

ans =
-4.4409e-16

25
√2 √2
2. Desarrollar (− + 𝑖)
2 2

(-sqrt(2)/2+sqrt(2)/2i)^25
ans =
-0.7071 - 0.7071i

3. Sea z1 = 2+3i, z2 = 4+2i, z3 = 3+i, z4 = iπ. Calcular:

a. A=z1-z2
b. B=z1*z2
c. C=z1/z2
d. D=z3^2
e. E=sqrt(z1)
f. F= 𝑒 𝑧4
g. G=z1’ (conjugada)

Solución

z1 = 2+3*i;
z2 =4+2*i;
z3 =3+i;
z4 =pi*i;
A=z1-z2
B=z1*z2
C=z1/z2
D=z3^2
E=sqrt(z1)
F=exp(z4)
G=z1'

A =
-2.0000 + 1.0000i
B =
2.0000 +16.0000i

C =
0.7000 + 0.4000i

D =
8.0000 + 6.0000i

E =
1.6741 + 0.8960i

F =
-1.0000 + 0.0000i

G =
a. - 3.0000i

4. Potencias de números complejos: Dado z=i.

a. Calcular v=Zn para n=1,...,40

b. Graficar dichos valores en el plano complejo.
c. Repetir los apartados anteriores con z=1.1i y z=1.1268i
d. ¿Qué ocurre con w cuando n crece indefinidamente?

%(a)

for n=1:40
v(n)=i^n;
end
% b) Grafica

plot(v)

% (c) Repetir los resultados cuando z=1.1i y z=1.1268i-1

for n=1:40
v(n)=(1.1*i)^n;
end
plot(v)
for n=1:40
v(n)=(1.1268*i-1)^n;
end
plot(v)

Cuando n tiende a infinito: se deja al lector que realice las pruebas.

6.8 Función Compass

La función compass toma las entradas de los datos y los pasa en formato
cartesiano, pero emite los gráficos polares. Esta función representa la longitud de
cada flecha corresponde a la magnitud de un elemento de los datos y su dirección
de orientación indica el ángulo de los datos complejos.

La función crea radiovectores que salen del origen de los ejes en un sistema de
coordenadas polares. Para ilustrar esta función, se crea un conjunto de flechas que
aumentan de tamaño en sentido contrario a las agujas del reloj.

t = 1:10;
r = t.*exp(i*t*36*(pi/180));
compass(r)
Graficar la función 𝑧 = 𝑐𝑜𝑠𝜃 + 𝑖𝑠𝑒𝑛𝜃, evaluada en el intervalo [-π, π] con particiones
de π/5.

th=-pi:pi/5:pi;
x=cos(th); y=sin(th);
z=x+y*i; compass(z)

Practica 2

1. Si z1=2-3i y z2=-7+8i. Determinar las siguientes operaciones:

a. z1z2
b. z1+z2
c. z1/z2
d. z1-z2
e. z1^2
 f. (z1^z2+z2/z1)z2’-z1’
g. √𝑧1𝑧2
h. exp(i*pi)
i. z1^3
j. (z2’)^10
k. 2z1+3z2-5z1’z2
l. (z1’/z2)^3
m. 𝑒 𝑧1
n. tanh(z1’*z2’)
ñ. z1/z2’

2. Graficar los puntos (-7+2i), (3+5i), (-2-4i), (1+3i), (3i), (6-2i) y (4). Ajuste los ejes
para mejorar la presentación de la gráfica. Colocar etiquetas, título y grilla.

3. Graficar los complejos dados en forma cartesiana, polar y trazar sus respectivos
radiovectores, todos en una misma ventana:

a. 10i
b. 2
c. 3+5i
d. 7-7i
e. 15+6i
f. -9-2i
g. -10+7i
h. sen30-icos30
i. –tan45 + isec25
j. 3cot60+2icsc40
k. 5cos50-5isen80
l. 𝑒 2𝜋𝑖
m. 2𝑒 3𝜋𝑖
5𝜋
n. 5𝑒 7 𝑖
o. √5 − 2√2𝑖

4. Sea 𝑥 = 2 + 2𝑖, 𝑦 = −7 − 3√3𝑖 calcular:

a. x/y
b. xyi
c. y3
d. x2/y90
e. y1/2
f. log(x)
g. 3x/y’
h. cos(xy)
3 𝑦′
i. √𝑥
 j. 𝑥𝑦 2

5𝜋 5𝜋
5. 𝑧 = 10 (𝑐𝑜𝑠 + 𝑖𝑠𝑒𝑛 4 ) Calcule la conjugada y determine los ángulos de
4
elevación y azimutal (en grados), al igual que la distancia con respecto a la
coordenada dada. Grafique la pareja ordenada con un círculo rojo en un plano
complejo con las debidas etiquetas y grilla.

3𝜋 3𝜋
6. 𝑧 = 9 (𝑐𝑜𝑠 − 𝑖𝑠𝑒𝑛 ) Calcule la magnitud, ángulo en radianes y grados.
7 7
Grafique en coordenadas polares.

𝜋 𝜋
7. 𝑧 = 2 (𝑐𝑜𝑠 + 𝑖𝑠𝑒𝑛 18) Calcule su conjugada y determine los ángulos de
18
elevación y azimutal (en grados), al igual que la distancia con respecto a la
coordenada dada.

8. 0.4(𝑐𝑜𝑠200𝑜 + 𝑖𝑠𝑒𝑛200𝑜 ) calcule su conjugada.

1+𝑐𝑜𝑠𝑥+𝑖𝑠𝑒𝑛𝑥 𝑛
9. Dada la ecuación: 𝑧 = ( )
1+𝑐𝑜𝑠𝑥−𝑖𝑠𝑒𝑛𝑥

Si x =30.5o, determine los valores de z cuando n = 0, 1, 2, 3, 4, 5, 6.

10. Sea la ecuación 3𝑧 2 − 5𝑧 + 9 = 0. Hallar los valores solución de z y pruebe que

cumplen con la condición de la ecuación.

11. Si a=5 y b=-7, evaluar:

𝑎 𝑏
𝑍 = (𝑎 + 𝑏𝑖) ( − 𝑖)
𝑎2 + 𝑏2 𝑎2 + 𝑏2

12. Sea la ecuación 2𝑧 2 + 3𝑧 + 7 = 0.

a. Hallar la magnitud
b. Hallar los ángulos en radianes y en grados.
c. z1*z2
d. (z’1/z2)^3

−3cosh −1(tan(𝑋)+𝑐𝑜𝑡ℎ(𝑥))+csch−1(1+senh −1 (𝑥))

13. Dada la función 𝑓 (𝑥 ) = calcular f(-π/7) y
cosh −1(2𝑥𝑠𝑒𝑐ℎ(𝑥))
f(π/5) como vector, luego determine las magnitudes y ángulos en radianes y
grados.

14. Evaluar las siguientes operaciones:

1−𝑖 5 𝑖 8 −𝑖 −8
𝑖 90 + 1 ( ) +1
𝑖 3−4𝑖

𝑖 𝑠𝑒𝑛(1+𝑖) (2𝑖 10 − 𝑖 15 )5 (2 + ln(𝑖))1/𝑖

 1−𝑖
(1 + 𝑖 )𝑖 𝑖 𝑙𝑛(1+𝑖) (1 + √3𝑖 )
3𝑖 5 +𝑖 −10
𝑖 70 + 𝑖 50 + 𝑖 25 + 𝑖
7𝑖

15. Dadas las impedancias Z1 = 3 - 5i y Z2 = 6 +3i dispuestas en paralelo en un

circuito. Determinar

a. Impedancia total:
𝑍1 𝑍2
𝑍=
𝑍1 + 𝑍2
b. Magnitud |𝑍| cuya unidad es el ohmio.
c. Ángulo de desfasamiento en radianes de la impedancia total 𝜃 = 𝑡𝑎𝑛−1 (𝑍)

16. El coeficiente de reflexión Гv de voltaje en una línea de transmisión está definido

por la fórmula:
𝑍𝐿 − 𝑍𝑜
Γ𝑣 =
𝑍𝐿 + 𝑍𝑜

Donde ZL Representa una impedancia de carga sobre la línea, y Z o la impedancia

de la línea, ambas representadas en ohmios (Ω).
Si ZL = 50.1+j50.1 y Zo = 47.5 Ω. Determinar:

a. Coeficiente de reflexión.
b. Magnitud del coeficiente de reflexión.
c. Ángulo del coeficiente de reflexión.
d. Si Гv = 057+j0.687 y Zo = 50.5 Ω. Cuál es el valor de ZL.

6.9 Función rose

Con la función rose (Rosa) se puede crear histogramas angulares que se dibujan
en coordenadas polares. Usando la función rose (ángulo), determinará cuántos de
los ángulos (en radianes) caen dentro de un cuadrante dado. Por defecto hay 20
compartimientos uniformemente espaciados entre 0 y 2π.

El número de contenedores puede cambiarse usando la función rose (ángulo,

escalar), el escalar especifica el número de compartimientos que deben estar
espaciados entre 0 y 2π. También, se puede especificar los centros de los
compartimentos pasando un vector vc, a la función rosa, así: rose (ángulo, vc).

Ejemplos

1. Dado el vector x, graficar empleando la instrucción rose.

x=[0,2,4,6,8,10];
rose(x)

2. El siguiente código produce un gráfico de la rosa cuyos datos se distribuyen

normalmente en un ángulo cerca de 90º.

x=angle(exp(i*randn(1,500)))+pi/2;
rose(x)

3. Dado el vector x definido en el intervalo [0,2π], graficar empleando la

instrucción rose
x=0:2*pi;
rose(x)
 4. Dado el ángulo definido en el intervalo [0,5π], graficar para un vector x=10.

ang=0:5*pi;
x=10;
rose(ang,x)

5. Graficar la función y = sen(x+10), evaluada en el intervalo [0,2π], con

particiones de π/20.

x=0:pi/20:2*pi;
y=sin(x+10);
rose(x,y)
title('x=sin(x+10)')
6.10 Instrucciones para conversión de sistemas de coordenadas

La conversión de un sistema de coordenadas a otro, es un tema clave en cálculo y

física en general. Para ello, Matlab implementa funciones específicas, cuya sintaxis
se resumen en el siguiente cuadro:

Instrucciones Transforma

[theta,rho,Z]=cart2pol(X,Y,Z) Cartesianas a cilíndricas

[theta,rho]=cart2pol(X,Y) Cartesianas a polares

[theta,phi,R]=cart2sph(X,Y,Z) Cartesianas a esféricas

[X,Y,Z]=pol2cart(theta,rho,Z) Cilíndricas a cartesianas

[X,Y]=pol2cart(theta,rho) Polares a cartesianas

[x,y,z]=sph2cart(theta,phi,R) Esféricas a cartesianas

Las ecuaciones que relacionan las conversiones son:

Cilíndricas a rectangulares

𝑥 = 𝑟𝑐𝑜𝑠𝜃 𝑦 = 𝑟𝑠𝑒𝑛𝜃 𝑧=𝑧

Rectangulares a cilíndricas

𝑦
𝑟 = √𝑥 2 + 𝑦 2 𝜃 = 𝑎𝑟𝑐𝑡𝑎𝑛 ( )
𝑥
𝑧=𝑧

Esféricas a rectangulares

𝑥 = 𝜌𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 𝑦 = 𝜌𝑠𝑒𝑛𝜑𝑠𝑒𝑛𝜃 𝑧 = 𝜌𝑐𝑜𝑠𝜃

Rectangulares a esféricas
 𝑦
𝜌 = √𝑥 2 + 𝑦 2 + 𝑧 2 𝜃 = 𝑎𝑟𝑐𝑡𝑎𝑛 ( )
𝑥

Ejemplos

1. Transformar las coordenadas (-1,5,8) a cilíndricas

[theta,rho,Z]=cart2pol(-1,5,8)
theta =
1.7682

rho =
5.0990
Z =
8

2. Transformar las coordenadas cilíndricas (π/3,7,2) a cartesianas

[X,Y,Z]=pol2cart(pi/3,7,2)
X =
3.5000

Y =
6.0622

Z =
2

3. Transformar de coordenadas polares a cartesianas el punto (-3,π/5)

[X,Y]=pol2cart(-3,pi/5)
X =
-0.6220

Y =
-0.0887

4. Transformar de coordenadas esféricas a cartesianas el punto (π/3, π/5, 3)

[x,y,z]=sph2cart(pi/3,pi/5,3)
x =
1.2135

y =
2.1019

z =
1.7634
Practica 3

1. Transformar de coordenadas esféricas a cartesianas el punto (2π/5, 3π/7, 10).

2. Transformar el punto de coordenadas (2,-5,9) a coordenadas cilíndricas,

esféricas y polares.

3. Transformar las coordenadas esféricas a cartesianas el punto (π/8,π/11,10) y

luego transformar a cilíndricas.

4. Representar las funciones dadas en gráfica independiente en el intervalo [1,5]

en forma de radiovectores.

a. 𝑟 = 𝑐𝑜𝑠(2𝑡)
b. 𝑟 = 2𝑡𝑠𝑒𝑛(2𝑡)
c. 𝑟 = −3𝑒 −2𝑡𝑖
d. 𝑟 = 3𝑡𝑖
e. 𝑟 = 𝑠𝑒𝑛(𝑡) + cos(𝑡) 𝑖
f. 𝑟 = 𝑡 2𝑖

5. Graficar en coordenadas polares y forma de radiovector la función 𝑧 = 2cos(𝜃 −

1
) + 𝑖𝑠𝑒𝑛(2𝜃), evaluada en el intervalo [-π, π] con particiones de π/10. Adicione
2
título a la gráfica.

6. La longitud l de la cuerda de un círculo de radio r se representa por la fórmula:

𝜃
𝑙(𝜃) = 2𝑟𝑠𝑒𝑛 ( )
2

Graficar l en polares y como radiovector, en el intervalo [0,2π] con particiones

de π/20, donde r es igual a 5.

7. La distancia d de un péndulo desde su posición de reposo en un tiempo t, se

representa por la ecuación.
2𝜋 2 0.35
𝑑 = −15.2𝑒 −0.5𝑡/30 𝑐𝑜𝑠 (√( ) − 3600 𝑡).
3

Graficar en coordenadas polares d en el intervalo [0,6].

8. Graficar la función dada como radiovector y polar en el intervalo [0,4π] con

particiones de π/10, luego pruebe con particiones de π/20, π/40 y π/100.

7𝜃
𝑟 = 3cos ( )
2
9. Graficar usando la función rose según la condición e intervalo dado:

a. x=10, 0 ≤ 𝜃 ≤ 5𝜋
b. −2𝜋 ≤ 𝜃 ≤ 5𝜋
c. −5𝜋/3 ≤ 𝜃 ≤ 3𝜋
d. x=5, 𝜋 ≤ 𝜃 ≤ 3𝜋
e. −𝜋/2 ≤ 𝜃 ≤ 2𝜋

10. Graficar en coordenadas polares las ecuaciones de la mariposa en el intervalo

[0,20π] con particiones de 0.01.

𝑟 = 𝑒 𝑐𝑜𝑠𝜃 − 2 ∗ cos(4𝜃)

𝜃
𝑟 = 𝑒 𝑐𝑜𝑠𝜃 − 2 ∗ cos(4𝜃) + 𝑐𝑜𝑠 5 ( )
12

11. Graficar en coordenadas polares en el intervalo [0,2π] con particiones de 0.01.

9
𝑟2 =
1 + 𝑠𝑒𝑛2 𝜃

12. Movimiento de cometas. La ecuación polar que se muestra a continuación,

describe el movimiento de cometas graficando sus trayectorias:

𝑑𝑝𝑒𝑟 (1 + 𝑒)
𝑟=
1 − 𝑒𝑐𝑜𝑠(𝜃)

Donde dper representa la distancia del perihelio medida en Unidades

astronómicas (UA); e es la excentricidad de la cónica.

a. Graficar en polar la órbita de saturno en el intervalo [0,2π] con particiones

de 0.01, donde e= 0.056 y dper = 9.006.

b. Graficar en polar la órbita del cometa Halley en el intervalo [0,2π] con

particiones de 0.01, donde e= 0.9673 y dper = 0.5871.
c. Graficar en polar la órbita del cometa Halley en el intervalo [-36, 36] con
particiones de 2. Luego con el intervalo [-24.24] con particiones de 1.6.

d. Determine qué tipo de órbita describe el cometa Enke: elíptica, parabólica

o hiperbólica. Emplear el intervalo [0,2π], con particiones de 0.01, donde e=
0.8499 y dper = 0.3317.