114   Capítulo 4   Manipulación de matrices MATLAB

TALLER FINAL - PARTE A:

Realizar en los grupos de laboratorio los ejemplos 4.2 y 4.3 usando el comando meshgrid, para el
ejemplo 4.2 realizar la tabla como se muestra en los ejercicios de práctica 4.2

Cell arrays:
Vectores o matrices de celdas (cell arrays), que son vectores o matrices cuyos elementos pueden ser
cualquier otro tipo de dato. En un arreglo ordinario todos sus elementos son números o cadenas de
caracteres. Sin embargo, en un arreglo de celdas, el primer elemento puede ser un número; el segundo
una matriz; el tercero una cadena de caracteres; el cuarto una estructura, etc.

4.2 PROBLEMAS CON DOS VARIABLES

Todos los cálculos realizados hasta el momento han usado sólo una variable. Desde luego,
la mayoría de los fenómenos físicos puede variar con muchos factores diferentes. En esta
sección se considera cómo realizar los mismos cálculos cuando las variables se representan
mediante vectores.
Considere los siguientes enunciados MATLAB:

Dado que x y y son escalares, es un cálculo sencillo: x  y = 15, o

Ahora vea lo que ocurre si x es una matriz y y todavía es un escalar:

regresa cinco valores de x. Dado que y todavía es un escalar con sólo un valor (5),

04 Moore.indd 114 16/5/07 13:46:37

 Sección 4.2   Problemas con dos variables   115

regresa

Todo esto todavía es revisión. Pero, ¿qué ocurre si ahora y es un vector? Entonces

regresa un enunciado de error

Este enunciado de error (las dimensiones internas de matriz deben concordar) le re- Idea clave: cuando
cuerda que el asterisco es el operador para multiplicación matricial, que no es lo que se quie- se formulan problemas
re. Se quiere el operador punto-asterisco (.*), que realizará una multiplicación elemento por con dos variables, las
elemento. Sin embargo, los dos vectores, x y y, necesitarán tener la misma longitud para este dimensiones de la matriz
propósito. En consecuencia, deben concordar.

crea un nuevo vector y con cinco elementos igualmente espaciados:

No obstante, aunque esta solución funciona, el resultado probablemente no es lo que en reali-

dad quiere. Puede pensar en los resultados como en la diagonal en una matriz (tabla 4.1).
¿Y si quiere conocer el resultado para el elemento 3 del vector x y el elemento 5 del
vector y? Obviamente, este enfoque no da todas las posibles respuestas. Se quiere una matriz
bidimensional de respuestas que corresponda a todas las combinaciones de x y y. Con la fi-
nalidad de que su respuesta, A, sea una matriz bidimensional, los vectores de entrada deben
ser matrices bidimensionales. MATLAB tiene una función interna llamada meshgrid, que le
ayudará a lograr esto, e incluso x y y no tienen que ser del mismo tamaño.
Primero cambie y de nuevo a un vector de tres elementos:

04 Moore.indd 115 16/5/07 13:46:46

 116   Capítulo 4   Manipulación de matrices MATLAB

Luego se usará meshgrid para crear una nueva versión bidimensional de x y y que se llamarán
new_x y new_y:

Idea clave: use la El comando meshgrid toma los dos vectores de entrada y crea dos matrices bidimen-
función meshgrid para sionales. Cada una de las matrices resultantes tiene el mismo número de filas y columnas. El
mapear dos variables número de columnas se determina por el número de elementos en el vector x, y el número de
unidimensionales en filas se determina mediante el número de elementos en el vector y. Esta operación se llama
variables bidimensionales mapeo de vectores en un arreglo bidimensional:
de igual tamaño.

Note que todas las filas en new_x son las mismas y todas las columnas en new_y son las
mismas. Ahora es posible multiplicar new_x por new_y y obtener la retícula bidimensional de
resultados que realmente se quiere:

Ejercicio de práctica 4.2

1. El área de un rectángulo es longitud por ancho (área = longitud 3 ancho).
Encuentre las áreas de los rectángulos con longitudes de 1, 3 y 5 cm y con
anchos de 2, 4, 6 y 8 cm. (Debe tener 12 respuestas.)
2. El volumen de un cilindro es volumen = pr2h. Encuentre el volumen de los
contenedores cilíndricos con radios desde 0 hasta 12 m y alturas desde 10
hasta 20 m. Aumente la dimensión del radio por 3 metros y la altura por
2 metros conforme abarca los dos rangos.

04 Moore.indd 116 16/5/07 13:46:52

Solución ejercicio de práctica 4.2
1. Área

1 long = [1 3 5];
2 ancho = [2 4 6 8];
3 [n_long, n_ancho] = meshgrid(long, ancho)

Command Window
n_long =
1.00 3.00 5.00
1.00 3.00 5.00
1.00 3.00 5.00
1.00 3.00 5.00
n_ancho =
2.00 2.00 2.00
4.00 4.00 4.00
6.00 6.00 6.00
8.00 8.00 8.00

4 area = n_long.*n_ancho

Command Window
area =
2.00 6.00 10.00
4.00 12.00 20.00
6.00 18.00 30.00
8.00 24.00 40.00

5 table1 = num2cell(area)

Command Window
table1 =
[2.00] [ 6.00] [10.00]
[4.00] [12.00] [20.00]
[6.00] [18.00] [30.00]
[8.00] [24.00] [40.00]

6 tableA = {'ancho\long' '1' '3' '5';...

'2' table1{1,1:end};...
'4' table1{2,1:end};...
'6' table1{3,1:end};...
'8' table1{4,1:end}}

Command Window
tableA =
'ancho\long' '1' '3' '5'
'2' [2.00] [ 6.00] [10.00]
'4' [4.00] [12.00] [20.00]
'6' [6.00] [18.00] [30.00]
'8' [8.00] [24.00] [40.00]

2. Volumen

1 radio = 0:3:12;
 2 altura = 10:2:20;
3 [n_radio, n_altura] = meshgrid(radio, altura)

Command Window
n_radio =
0 3.00 6.00 9.00 12.00
0 3.00 6.00 9.00 12.00
0 3.00 6.00 9.00 12.00
0 3.00 6.00 9.00 12.00
0 3.00 6.00 9.00 12.00
0 3.00 6.00 9.00 12.00
n_altura =
10.00 10.00 10.00 10.00 10.00
12.00 12.00 12.00 12.00 12.00
14.00 14.00 14.00 14.00 14.00
16.00 16.00 16.00 16.00 16.00
18.00 18.00 18.00 18.00 18.00
20.00 20.00 20.00 20.00 20.00

4 volumen = pi*n_radio.^2.*n_altura

Command Window
volumen =
0 282.74 1130.97 2544.69 4523.89
0 339.29 1357.17 3053.63 5428.67
0 395.84 1583.36 3562.57 6333.45
0 452.39 1809.56 4071.50 7238.23
0 508.94 2035.75 4580.44 8143.01
0 565.49 2261.95 5089.38 9047.79

5 table1 = num2cell(volumen)

Command Window
table1 =
[0] [282.74] [1130.97] [2544.69] [4523.89]
[0] [339.29] [1357.17] [3053.63] [5428.67]
[0] [395.84] [1583.36] [3562.57] [6333.45]
[0] [452.39] [1809.56] [4071.50] [7238.23]
[0] [508.94] [2035.75] [4580.44] [8143.01]
[0] [565.49] [2261.95] [5089.38] [9047.79]

6 tableV = {'altura\radio' '0' '3' '6' '9' '12';...

'10' table1{1, 1:end};...
'12' table1{2, 1:end};...
'14' table1{3, 1:end};...
'16' table1{4, 1:end};...
'18' table1{5, 1:end};...
'20' table1{6, 1:end}}

Command Window
tableV =
'altura\radio' '0' '3' '6' '9' '12'
'10' [0] [282.74] [1130.97] [2544.69] [4523.89]
'12' [0] [339.29] [1357.17] [3053.63] [5428.67]
'14' [0] [395.84] [1583.36] [3562.57] [6333.45]
'16' [0] [452.39] [1809.56] [4071.50] [7238.23]
'18' [0] [508.94] [2035.75] [4580.44] [8143.01]
'20' [0] [565.49] [2261.95] [5089.38] [9047.79]
 Sección 4.2   Problemas con dos variables   117

EJEMPLO 4.2
Distancia al horizonte
Probablemente ha experimentado estar de pie en lo alto de una colina o montaña y sentido que
puede ver hasta el infinito. ¿Realmente cuán lejos puede ver? Depende de la altura de la mon-
taña y del radio de la Tierra, como se muestra en la figura 4.2. La distancia hasta el horizonte
es muy diferente en la Luna que en la Tierra, porque el radio es diferente para cada una.
Con el teorema de Pitágoras se ve que

A partir de esta última expresión, encuentre la distancia hasta el horizonte en la Tierra y

en la Luna, para montañas desde 0 hasta 8000 metros. (El monte Everest tiene 8850 metros de
alto.) El radio de la Tierra es 6378 km y el de la Luna es de 1737 km.
1. Establezca el problema.
Encontrar la distancia hasta el horizonte desde lo alto de una montaña en la Luna y en
la Tierra.
2. Describa las entradas y salidas.
Entrada
Radio de la Luna 1737 km
Radio de la Tierra 6378 km
Altura de las montañas 0 a 8000 metros
Salida
Distancia hasta el horizonte, en kilómetros
3. Desarrolle un ejemplo a mano.

Figura 4.2
Distancia hasta el
horizonte.

04 Moore.indd 117 16/5/07 13:46:56

 118   Capítulo 4   Manipulación de matrices MATLAB

Con el radio de la Tierra y una montaña de 8000 metros se obtiene

4. Desarrolle una solución MATLAB.

EJEMPLO 4.3
Caída libre
La ecuación general para la distancia que recorre un cuerpo en caída libre (sin tomar en cuenta
la fricción del aire) es

04 Moore.indd 118 16/5/07 13:47:02

 Sección 4.2   Problemas con dos variables   119

donde
d = distancia,
g = aceleración debida a la gravedad, y
t = tiempo.

Cuando un satélite orbita un planeta, está en caída libre. Muchas personas creen que,
cuando el transbordador espacial entra en órbita, deja detrás la gravedad; pero la gravedad es
lo que mantiene al transbordador en órbita. El transbordador (o cualquier satélite) en realidad
cae hacia la Tierra (figura 4.3). Si va lo suficientemente rápido de manera horizontal, perma-
nece en órbita; si va muy lentamente, golpea el suelo.
El valor de la constante g, la aceleración debida a la gravedad, depende de la masa del
planeta. En diferentes planetas, g tiene diferentes valores (tabla 4.2).
Encuentre qué tan lejos caería un objeto en tiempos desde 0 hasta 100 segundos en cada
uno de los planetas del sistema solar y en la Luna.
1. Establezca el problema.
Encontrar la distancia que recorre un objeto en caída libre en planetas con diferentes
gravedades.
2. Describa las entradas y salidas
Entrada Valor de g, la aceleración debida a la gravedad, en cada uno de los planetas
y la Luna
Tiempo = 0 a 100 s
Salida Distancias calculadas para cada planeta y la Luna
3. Desarrolle un ejemplo a mano.
d = ½ gt2, así que en Mercurio, a 100 segundos:
d = ½ 3 3.7 m/s2 3 1002 s2
d = 18,500 m

Figura 4.3
El transbordador espacial
cae constantemente hacia
la Tierra. (Cortesía de
NASA/Jet Propulsion
Laboratory.)

04 Moore.indd 119 16/5/07 13:47:05

 120   Capítulo 4   Manipulación de matrices MATLAB

Tabla 4.2 Aceleración debida a la gravedad en el sistema solar

4. Desarrolle una solución MATLAB.

04 Moore.indd 120 16/5/07 13:47:10

PARTE B:

1) MATLAB dispone de cuatro funciones básicas para crear gráficas 2-D. Estas
funciones se diferencian principalmente por el tipo de escala que utilizan en los ejes
de abscisas y de ordenadas. Consultar cuales son esas funciones y cómo se utilizan.

2) Existen además otras funciones orientadas a añadir títulos al gráfico, a cada uno de
los ejes, a dibujar una cuadrícula auxiliar, a introducir texto, etc. Consultar qué
funciones existen y cómo se utilizan.

3) Consultar cómo cambiar el grosor de una gráfica, el tamaño del texto, el color, si se
quiere una curva con linea continua o discontinua, cómo graficar varias curvas en una
sola, cómo poner leyendas, cómo poner o quitar la cuadricula, cómo generar una
matriz de gráficas.

4) Para construir gráficos en tres dimensiones primero se construye la superficie que

representa la magnitud de una función de dos variables; las dos variables se representan
en el plano horizontal y el valor de la función como una altura sobre dicho plano. Consultar las
funciones utilizadas para graficar superficies y cómo graficarlas usando la función meshgrid
para construir el plano XY.

5) Consultar las funciones que permiten definir las características de los gráficos 3-D como
se plantea en los puntos 2) y 3)

6) Crear un script en donde se ejempliquen las funciones consultadas, separar cada una de las
graficas 2-D y 3-D en secciones. Ustedes definen qué funciones desean graficar.