Operaciones básicas con Matlab

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

3.1 Operaciones aritméticas

Matlab en cuanto a operaciones básicas se refiere, actúa de manera equivalente a

cualquier calculadora, contando con una alta capacidad de procesamiento de
información. Por lo tanto, la velocidad de cómputo depende de los procesadores del
equipo donde está instalado el programa de Matlab.
Las operaciones aritméticas indicadas se evalúan por orden prioritario, es decir,
primero las potencias, luego las multiplicaciones, le sigue las divisiones y finalmente
las sumas y restas.

Suma + a+b
Diferencia - a-b
Multiplicación * a*b
División / a/b
Elevado a una potencia ^ a^n

Para el caso de operaciones algebraicas que tienen igual prioridad, Matlab evalúa
las operaciones de izquierda a derecha, de manera equivalente a como se realizaría
en papel. Por ejemplo, una suma:

>> 5+9 < enter >

El resultado es desplegado y se guarda en la variable ans (answer).

ans =
14

Para las demás operaciones, se tiene:

Resta Multiplicación División Potencia

>> 45-28 >> 8*4 >> 89/5 >> 12^2
ans = ans = ans = ans =
17 32 17.8000 144

Para operaciones combinadas. Obtener el valor de la siguiente expresión:

1
310 + 8 + 10000
0.7258 − 2(1 − √5)
En Matlab se teclea

>> (3^10+1/8+10000)/(0.7258-2*(1-sqrt(5)))
ans =
2.1592e+04

Como ocurre con las calculadores científicas, la notación 2.1592e+04 significa

21592.

En Matlab existen una vasta serie de instrucciones relacionadas con las

operaciones matemáticas, que van desde las más básico hasta lo más complejo, tal
como se resume en el siguiente cuadro:

Operación Instrucción Descripción

a+b plus(ab) Suma

a-b minus(a,b) Resta

-a uminus(a) Menos unitario

+a uplus(a) Más unitario

a.*b times(a,b) Multiplicación por elemento

a*b mtimes(a,b) Multiplicación matricial

a./b rdivide(a,b) División derecha por elemento

a.\b ldivide(a,b) División izquierda por elemento

a/b mrdivide(a,b) División matricial derecha

 a\b mldivide(a,b) División matricial izquierda
a.^b power(a,b) Potencia por elemento

a^b mpower(a,b) Potencia matricial

a<b lt(a,b) Menor que

a>b gt(a,b) Mayor que

a <= b le(a,b) Menor o igual a

a >= b ge(a,b) Mayor o igual a

a ~= b ne(a,b) Distinto de

a == b eq(a,b) Igual a

a&b and(a,b) Y lógico

a|b or(a,b) O lógico

~a not(a) No lógico

a:d:b colon(a,d,b) Operador dos puntos

a:b colon(a,b)

a’ ctranspose(a) Traspuesta conjugada compleja

a.’ transpose(a) Matriz transpuesta

display(a) Visualización pantalla

[a b] horzcat(a,b,...) Concatenación horizontal

[a; b] vertcat(a,b,...) Concatenación vertical

a(s1,s2,...sn) subsref(a,s) Referencia por subíndices

a(s1,...,sn)= subsasgn(a,s,b) Asignación por subíndices

b

(a) subsindex(a) Conversión al ser índice

3.2 Prefijos del Sistema Internacional

Los prefijos pertenecientes al SI los fija la Oficina Internacional de Pesas y

Medidas (Bureau International des Poids et Mesures), de acuerdo con el cuadro
siguiente:

10n Prefijo Símbolo Escala larga

10 24 Yotta Y Cuatrillón
1021 Zetta Z Mil trillones
 1018 Exa E Trillón
1015 Peta P Mil billones
1012 Tera T Billón
109 Giga G Mil millones
106 Mega M Millón
103 Kilo K Mil / Millar
102 Hecto H Cien / Centena
101 Deca Da Diez / Decena
1 ninguno Uno / Unidad
10−1 deci d Décimo
10−2 centi c Centésimo
10−3 mili m Milésimo
10−6 micro µ Millonésimo
10−9 nano n Milmillonésimo
10−12 pico p Billonésimo
10−15 femto f Milbillonésimo
10−18 atto a Trillonésimo
10−21 zepto z Miltrillonésimo
10−24 yocto y Cuatrillonésimo

Es importante conocer y usar correctamente los prefijos del sistema internacional,

porque en ciencias e ingeniería se emplean de forma regular, y si el uso es erróneo
puede conllevar a obtener datos inconsistentes. Para el caso de Matlab la notación
en base 10 se simboliza por la letra e. Es importante que se familiarice con este
símbolo para no incurrir en errores como confundir esta letra con la función
exponencial.

3.3 Unidades básicas del Sistema Internacional

Es importante recordar las unidades básicas del SI para el desarrollo de los

ejercicios planteados en cada capítulo, las cuales se citan en el siguiente cuadro.

Magnitud Nombre Símbolo

Longitud metro m
Masa kilogramo kg
Tiempo segundo s
Intensidad de corriente eléctrica ampere A
Temperatura termodinámica kelvin K
Cantidad de sustancia mol mol
Intensidad luminosa candela cd
Unidades del SI derivadas, estas se expresan a partir de unidades básicas y
suplementarias, como se observa en el siguiente cuadro:

Magnitud Nombre Símbolo

Superficie metro cuadrado m2
Volumen metro cúbico m3
Velocidad metro por segundo m/s
Aceleración metro por segundo cuadrado m/s2
Número de ondas metro a la potencia menos uno m-1
Densidad kilogramo por metro cúbico kg/m3
Velocidad angular radián por segundo rad/s
Aceleración angular radián por segundo cuadrado rad/s2

Ejemplos

2.5x105 en Matlab es 2.5e5 o 2.5*10^5

4.68 x10-8 en Matlab es 4.68e-8 o 4.68*10^-8

7Mm en Matlab se digita 7e6 m

9.85 µs en Matlab se digita 9.85e-6 s

Para operaciones como 4.7x106*1.25x102 en Matlab se digita:

>> 4.7e6*1.25e2
ans =
587500000

3.4 Práctica 1

Resolver los ejercicios en el command Window y los problemas en un script:

1. 6-17+9(25+8)9-5/6+6-82

2. 7+15(5+19)+ (9-4)(6+3)/8

3. (4.225+5.14578)2-(51.3)(645)

4. 12.2658/45863 - 3.98752/2.3587+ 46-54+67 - 4.358.5452

5. (6+20)-5(8-7/6) + (7/9-3/5)6+12(53-24)

6. (4.575+2.508)3/54 + 3.256 (8.581.256)5

7. 4.7M*3.881T/7.5n

8. 4.5e-5*5.471e8+3.12e3-2.4578e2

9. 2.254 x10-9 0.89775 x10-8(4x1020)6

10. (0.45682658+2.2568)6/(3.65877+5.325684)2

11. (5.32e3*8.256e-2+1.214e-4*(-2.2455))5

48 45−19 25
12. 65 + 6+15 - + 13+17
38

13. 6 Mm+0.025 Gm+2.5 Km

14. 27.5 mm+8.397 µm

15. (5.7854 ns+0.085 µs)/42.135fs

16. 1024/1018+106 - 10−18 *10−15

17. La cantidad de bienes y servicios BS disponibles en una economía nacional

puede ser calculado mediante la siguiente ecuación:

𝐵𝑆 = 𝐶 + 𝐼 + 𝑋

Donde C es el consumo total (privado y público); I, inversión pública y privada;

X, valor de todas las exportaciones y M, valor de todas las importaciones.

Si en un país el consumo total es de 40.5217 T€, la inversión total es de 31.259

G€, el valor de exportaciones 0.8945 T€ y el valor de importaciones 987.578 M€.

a. Determine la cantidad de bienes y servicios.

b. Diferencia entre exportaciones e importaciones.

18. Un cuerpo describe un movimiento uniforme acelerado, el cual está

representado por la ecuación cinemática,

𝑉𝑓 = 𝑉𝑖 + 𝑎𝑡

Donde 𝑉𝑓 y 𝑉𝑖 son la velocidad final e inicial respectivamente en m/s, a, es la

aceleración en m/s2, y t el tiempo en segundos.
 a. Si la velocidad final e inicial es de 15 m/s y 8 m/s respectivamente, determine
los valores de la aceleración para tiempos de 2, 3, 4 y 5 segundos.
b. Cuál es el valor de la aceleración para una diferencia de velocidades de 60
Km/h en un tiempo de 0.12 minutos.

19. Determine el perímetro del terreno rectangular, según los datos de la figura,
cuando x equivale a 10.1253 m
5+x

2x-1

20. La pérdida de potencia P (Watts) en un circuito está representado por la

siguiente ecuación:
𝑉2
𝑃=
𝑅

Donde V es el voltaje (volts) y R la resistencia eléctrica (Ohmios=Ω)

a. Si R = 1.51 KΩ y P = 5.82 W, cuál es el valor del voltaje.

b. Si P = 0.5678 MW y V = 0.253 KV, cuál es el valor de la resistencia R.

21. Calcular la densidad crítica ρc (en kg/m3) estimada del universo mediante la
fórmula:
3𝐻𝑜2
𝜌𝑐 =
8𝜋𝐺

Donde Ho =67.8± 0.7 Km/sMpc y G la constante de gravitación universal.

22. La humedad específica (𝐻𝑒 ) es la relación entre la masa de vapor de agua y la

masa de aire húmedo, la cual se calcula por la fórmula:

0.622𝑃𝑣
𝐻𝑒 =
𝑃 − 0.378𝑃𝑣

Donde P es la presión atmosférica medida en hectopascales (Hpa); Pv es la

presión de vapor en Hpa; He se mide en un valor porcentual. Determine la
humedad específica cuando P = 2HPa y Pv toma valores de 0.9, 1.2 y 1.5 HPa.

23. La densidad de la Tierra D(h) (en g/cm3) a una altitud h metros bajo la superficie
terrestre, se puede aproximar con la fórmula: (Fuente. Swokowski, Cole, Álgebra
y trigonometría con geometría analítica. 12ª. Edición. Cengage Learning, pp.
308)

𝐷(ℎ) = 2.84 + 𝑎ℎ + 𝑏ℎ2 − 𝑐ℎ3

 Donde a = 1.4x10-3, b=2.49x10-6, c=2.19x10-9 y 0 ≤ ℎ ≤ 1000, determinar
D(150), D(400), D(800), D(1000).

24. La longitud l, de cada lado de cualquier pirámide cuadrangular se determina por

medio de la ecuación:
3𝑉
𝑙=√
ℎ

Donde V, representa el volumen (en pies3) de la pirámide y h la altura en pies.

Si el volumen de una pirámide es de 49.750.987.54477 pies 3, cuál es su altura

si la longitud de uno de sus lados es de 643.87891 pies.

25. En un movimiento uniforme y rectilíneo la distancia expresada en forma escalar

es:

𝑥(𝑡) = 𝑥𝑜 + 𝑣𝑜 𝑡

a. Determine la posición final de un objeto que ha recorrido inicialmente 0.6 Km

en 2 minutos con una velocidad inicial de 15 Km/h.
b. La diferencia de posición de un objeto es 120 m, determine la velocidad de
partida si el tiempo de medición es de 1.35 minutos.

26. La altura que alcanza un cuerpo que describe un movimiento parabólico está
descrita por la ecuación:
ℎ(𝑡) = −17𝑡 2 + 522𝑡 + 74

Determine la altura al cabo de 0, 1.247 y 3.712 segundos.

27. Para calcular la potencia eléctrica (Watt) mínima de un aerogenerador se emplea

la fórmula:
1
𝑃𝑒 = 𝐶𝑝 𝜌𝜋𝑅 2 𝑉 3
2

Donde 𝐶𝑝 = 0.351, 𝜌 es la densidad del aire (1.25 Kg/m3), R, radio del

aerogenerador (m) y V la velocidad de rotación (m/s).

a. Determine la potencia cuando R=1.4 y la velocidad es de 15.2 m/s.

b. Si la potencia de un aerogenerador es de 2386.2 W, determine la velocidad
de rotación cuando el radio es de 1.5 m.

28. Las fórmulas de la densidad 𝜌 y peso específico 𝛾 están representados por las
siguientes fórmulas:
 𝛾 𝑃
𝜌=𝑔 𝛾 = 𝑇𝑅

Donde g es la gravedad 9.81 m/s2; P, presión dada en kg/m2; T, temperatura en

grados Kelvin; R, constante de los gases, 29.3 m/oK

Determine la densidad del aíre a 30 oC, donde la presión es de 10336 Kg/m2.

29. El rendimiento térmico viene dado por la ecuación:

𝑇1 −𝑇2
𝑅= ∗ 100%, 𝑇1 > 𝑇2
𝑇1

a. Hallar el rendimiento ideal de una máquina térmica que funciona entre 230.45
oC y 70.51 oC. Recuerde que debe pasar grados centígrados a kelvin.

b. Cuál debe ser la temperatura del foco caliente para que el rendimiento sea
del 50.5 %. Tome el valor de 𝑇2 del ejercicio anterior.

30. El principio de Arquímedes afirma que un cuerpo total o parcialmente sumergido

en un fluido en reposo, recibe un empuje de abajo hacia arriba igual al peso del
volumen del fluido que desaloja.

Si el empuje E que experimenta un cuerpo es de 12.5 KN, que se encuentra en

agua de mar (Densidad 𝜌 = 1027 Kg/m3), ¿Cuál es el volumen desplazado?

𝐸 = 𝜌𝑔𝑉

Donde g es la aceleración de la gravedad (9.81 m/s2).

31. Un divisor de voltaje es una configuración de circuito eléctrico que reparte el

voltaje de una fuente (Vin) entre una o más resistencias (R1, R2) conectadas en
serie (una a continuación de otra), tal como se muestra en la figura 3.2.

𝑅2
𝑉𝑜𝑢𝑡 = 𝑉
𝑅1 +𝑅2 𝑖𝑛

a. Determinar el voltaje de salida cuando las resistencias

son de 4 KΩ y 5.5 KΩ (R2), y voltaje de entrada es de
70.5231 Voltios.
Figura 3.2

b. Determinar la resistencia R1, si el voltaje de salida es 1/3 del voltaje de

entrada y la resistencia R2 es 3.58 KΩ.

32. Para levantar una plataforma de 10.12 toneladas se utiliza un gato hidráulico. Si
en el pistón actúa una presión de 14.52 Kg/cm 2, cuál es el área del pistón.
 𝑃𝑒𝑠𝑜
𝑃𝑝𝑖𝑠𝑡ó𝑛 =
Á𝑟𝑒𝑎

Si la masa de un pistón es de 850 gramos, distribuida sobre una superficie de 5

cm por 6 cm, ¿qué presión ejerce en unidades de pascales, atmósferas y bares?

33. El área de superficie corporal S de una persona (en pies 2) se puede aproximar
con:
𝑆 = 0.1091𝑤 0.425 ℎ0.725

Donde h es la estatura medida en pulgadas y w el peso en libras.

Hallar S para una persona que mide 6.4 pies de alto y pesa 165.78 libras.

34. La máxima elevación de las mareas, es el doble en valor absoluto, de la mínima

elevación. De modo que, la diferencia entre altura máxima de la bajamar y la
pleamar es,
𝑀 𝑅 3
∆ℎ𝑚á𝑥 = ( ) 𝑅
2𝑀𝑇 𝑟

Donde M es la masa del cuerpo perturbador, que en este caso es el Sol

(1.98·1030 kg) y la Luna (7.35·1022 kg); MT = 5.98·1024 kg, es la masa de la Tierra;
El radio de la Tierra R=6.37·106 m.

La distancia media entre el centro de la Tierra y el centro de la Luna es 384.4x106

m.
La distancia media entre el centro de la Tierra y el centro del Sol es 149.6x109
m.

Determine los valores máximos para:

a. Las mareas producidas por la Luna.

b. Las mareas producidas por el Sol.

35. Calcular el valor máximo de las fuerzas de marea debida a la Luna F L y al Sol
FS, según las siguientes ecuaciones:

2𝐺𝑀𝐿 𝑅 2𝐺𝑀𝑆 𝑅
𝐹𝐿 = 𝐹𝑆 =
𝑟2 𝑟 𝑟2 𝑟

Donde G es la constante de gravitación universal que equivale a: 6.67x10-

11 Nm2/kg2

Halle el cociente entre estas dos fuerzas FL/FS, que indica que las fuerzas de
marea son muy pequeñas comparadas con la fuerza de atracción terrestre sobre
 un objeto de masa m situado en su superficie, aunque no se desestiman sus
efectos tanto en el clima como en la vida misma.

36. Dado el sistema de cuerpos m1 y m2 ligados por una

cuerda tal como se muestra en la figura 3.3, cuya
tensión es transmitida a lo largo de todos los puntos
y es igual en los dos extremos. El sistema se mueve
por la acción del peso del cuerpo 2, por lo cual la
aceleración total del sistema y la tensión son
respectivamente:

Figura 3.3
𝑚2 𝑔 𝑚1 𝑚2
𝑎=𝑚 ; 𝑇=𝑚 𝑔
1 +𝑚2 1 +𝑚2

a. Si las masas m1=2.3 Kg y m2 = 6.89 Kg, determine la aceleración y tensión

del sistema (g=9.81 m/s2).
b. Si la aceleración de los cuerpos es de 4 m/s2, cuál es la tensión de la cuerda,
si m2 = 7 Kg.
c. Si la tensión de la cuerda es de 25 N, con qué aceleración se mueve el
sistema, en la que se sabe que m1=2.85 Kg

37. Tomando como referencia el punto anterior, la aceleración y tensión cuando el

sistema de cuerpos está sujeto a fricción 𝜇 , se tiene:
𝑚2 −𝜇𝑚1 𝑚 𝑚
𝑎= 𝑔; 𝑇 = 𝑚 1+𝑚2 (1 + 𝜇)𝑔
𝑚1 +𝑚2 1 2

a. Calcular la aceleración y tensión cuando el coeficiente de rozamiento es

0.456. Tome los valores de las masas del punto anterior y compare los
resultados.
b. Cuál es el coeficiente de rozamiento del sistema de masas, cuya tensión es
de 37.8512 N y las masas son de 4.2 Kg y 8.3 Kg.
c. Cuál es la aceleración según los datos del punto b.

38. La Fuerza de Casimir por unidad de área Fc/A para placas ideales y
perfectamente conductoras con vacío entre ambas es:

𝐹𝑐 ℏ𝑐𝜋 2
=−
𝐴 240𝑎4

Donde ℏ es la constante reducida de Planck; c es la velocidad de la luz; a es

la distancia entre dos placas. La fuerza es negativa, indicando el carácter
atractivo de la misma: disminuyendo la distancia entre placas, la energía es
reducida.
 Determine la Fuerza de Casimir por unidad de área cuando la distancia entre las
placas es de 1 μm, 10 μm, 1 nm, 10 nm.

39. La densidad atmosférica terrestre D(h) (en kg/m 3) a una altitud de h metros se
puede aproximar con la fórmula: (Fuente. Swokowski, Cole, Álgebra y
trigonometría con geometría analítica. 12ª. Edición. Cengage Learning, pp. 308)

𝐷(ℎ) = 1.2 − 𝑎ℎ + 𝑏ℎ2 − 𝑐ℎ3

Donde a = 1.096x10-4, b=3.42x10-9, c=3.6x10-14 y 0 ≤ ℎ ≤ 30000.

Determinar D(10000), D(19000), D(23000) y D(30000).

40. El número de carga de rodamiento hidrodinámico está dado por la ecuación:

𝜋𝜀√𝜋 2 (1 − 𝜀 2 ) + 16𝜀 2
𝑁𝐻 =
(1 − 𝜀 2 )2

Donde 𝜀 representa el coeficiente de excentricidad. Calcular NH cuando 𝜀 = 0.81

y 𝜀 = 0.75.

41. Área o sección resistente efectiva. La sección resistente para los tornillos
métricos se puede determinar por:
𝜋
𝐴𝑠 = (𝑑 − 0.938194𝑃)2
4

Donde d es el diámetro nominal de la rosca del tornillo y P el paso de la rosca.

Complete el cuadro según los datos suministrados:

Referencias DIN Paso P (mm) Diámetro As (mm2)

ISO 13-1 nominal (mm)

M1 0.25 0.693
M1,2 0.25 0.893
M1,6 0.35 1.170
M2 0.4 1.509
M2,3 0.45 1.948
M3 0.5 2.387
M4 0.7 3.242
M5 0.8 4.134
Fuente: Rosca Métrica ISO DIN 13 (M).

42. El periodo de oscilación de una estructura, está definido como el tiempo

necesario para que un edificio realice un ciclo completo en vibración libre:
 𝑇 = 𝐶𝑍 𝛼

Donde Z es la altura total del edificio, C y α son valores empíricos que dependen
del tipo de estructura. Así, C = 0,047, α = 0,9 para edificios de pórticos de
concreto, C = 0,072, α = 0,8 para edificios aporticados de acero y C = 0,073, α
= 0,75 para edificios de acero con pórticos arriostrados.

Determine y compare el periodo de oscilación para edificios cuya altura es de 20

metros, conforme a los datos suministrados de C y α.

43. El área total de un cilindro circular recto es: 𝐴 = 2𝜋𝑟(ℎ + 𝑟)

Donde r es el radio de la circunferencia de la base del cilindro y h la altura.

Si el área de un depósito cilíndrico es de 100,81 m 2, determine el diámetro
óptimo para una altura de 6.1251 m.

44. Calcule la velocidad de propagación de una onda electromagnética en el vacío,

empleando la ecuación:
1
𝑐=
√𝜇𝑜 𝜀𝑜

Donde μo y εo representan las constante de permeabilidad magnética y

permitividad eléctrica respectivamente, que equivalen a 1.3x10 -6 N/A2 y 8.9x10-
12 C2/Nm2

Usar las constantes dadas y calcular la resistencia eléctrica del vacío empleando
la siguiente ecuación:
𝜇𝑜
𝑅=√
𝜀𝑜

45. Efecto invernadero. El promedio de aumento de temperatura de la superficie del

planeta Tierra debido al efecto invernadero se calcula con la ecuación:

21.1 3 127.01 2 1293

𝑇(𝑡) = 𝑡 − 𝑡 + 𝑡
5. 106 1. 106 5. 104

Donde 0 ≤ 𝑡 ≤ 60, y 𝑡 = 0 corresponde a 1980. Calcular la temperatura en

grados celcius para los años 1995, 2000, 2011, 2016, 2018, 2020 y 2035.

46. El cambio de energía cinética ∆𝐸 de un sistema se expresa mediante la

ecuación:
1
∆𝐸 = 𝑚[𝑣22 − 𝑣12 ]
2

Donde m es la masa en Kg; v2 es la velocidad final, y v1, la velocidad inicial.

 a. Determine el cambio de energía cinética para un cuerpo de masa (m) 0.2 g
que inicialmente se desplazaba a 50 km/h y luego la triplicó.
b. Si el cambio de energía cinética es de 13.75 J para un cuerpo de masa de
100 g con una velocidad final de 30 m/s. ¿Cuál era su velocidad inicial?
c. El cambio de energía cinética de un cuerpo es de 1293.6 J, donde sus
velocidades son de 60 y 38 m/s. ¿Cuál es la masa del cuerpo?

47. Un cable de un puente colgante tiene sus soportes en el mismo nivel, separados
a una distancia de L pies. Los soportes están a “a” pies por encima del punto
mínimo del cable, tal como se observa en la figura 3.5.

Figura 3.5. Cable colgante de un puente que emplea dos soportes como punto de apoyo.

Si el peso del cable es despreciable pero el puente tiene un peso uniforme de w

libras por pies. Se deducen las siguientes fórmulas:

a. La tensión en libras para un cable en el punto más bajo es:

𝑊𝐿2
𝑇𝑐 =
8𝑎

b. La tensión en libras de los soportes es:

𝑊𝐿 2
𝑇𝑠 = √𝐿 + 16𝑎2
8𝑎

Determinar las tensiones del cable y soporte, cuando los soportes equidistan a
350.8912 pies, donde la distancia de estos está por encima del punto mínimo
del cable a 48.5812 pies. El peso del puente es de 5.12x10 6 libras
aproximadamente.

3.5 Operadores relacionales

Los operadores relacionales se emplean para realizar comparaciones entre dos

valores numéricos. Por consiguiente, si la comparación es correcta, la expresión
se considera como verdadera, dando como resultado 1, en caso contrario, falsa y
resultado 0.
>> true >> false
ans = ans =
1 0
Cuadro de operadores relacionales

Representación
Operador Nombre significado
Matemática
< menor que a<b a menor que b
> mayor que a>b a mayor que b
== igual a a==b a igual a b
~= no igual a a~=b a no es igual a b
menor o igual a menor o igual
<= a<=b
a ab
mayor o igual a mayor o igual
>= a>=b
a ab

Ejemplo del uso de los operadores:

>> 2>3
ans =
0

El resultado (0) indica que la desigualdad es falsa porque 2 no es mayor a 3.

>> 2<3
ans =
1

El resultado (1) indica que la desigualdad es verdadera porque 2 es menor a 3.

Para los demás ejemplos puede deducir los estados de falso o verdadero según el
caso:

>> 2>=3 >> 2>=2 >> 2<=3 >> 2==2

ans = ans = ans = ans =
0 1 1 1
Devuelve Falso Devuelve Devuelve Devuelve
verdadero verdadero verdadero

Cuando es diferente

>> 2~=2
ans =
0

Devuelve Falso

Nota 5. El carácter ~ se obtiene al oprimir las teclas Alt + 126.

 Nota 6. No se debe confundir el operador asignación =, con el operador
relacional igual a ==

Ejemplo de operación de asignación,

>> x=5;
>> whos x

Name Size Bytes Class Attributes

x 1x1 8 double

Se asigna una variable x de valor cinco, cuyo atributo es real de tamaño 1x1,
compuesto por 8 bytes.

3.6 Práctica 2

Desarrolle las siguientes operaciones relacionales, observe los resultados.

1. 5<9
2. -5>-10
3. -5>8
4. 8>=9
5. -3<-2
6. 9==8
7. 10~=12
8. 9<=12
9. 7>9
10. 9<15
11. 8~=8
12. 2>=1
13. 3<=4
14. -7<-4
15. 9>=-5
16. -11<=-13
17. 6~=-9
18. -9==9
19. -5~=-5
20. (6-5) > (6+2)
21. 63 > 53
22. (18+25) < (36+28)
23. (5*9) <= (1589/12)
24. (15-6.32)4 >= (25.4587/10)
25. 53 >= 62
26. (-4)3 ~= -(4)3
27. (25+15)2 <= (1589*12)
28. 7 > 5 – 2x ≥ 3 con x = 1
29. √𝑥 + 2 > √2𝑥 + 3, con x = 2

3𝑥+5 4𝑥 2 −1
30. ≥ , con x = 3
𝑥 2 +7 5𝑥+1

𝑥+2
31. −10 < < 10 con x = 1, x = 4
3

32. (2𝑥 − 5)(4𝑥 + 7) ≤ (8𝑥 + 2)(𝑥 − 5) con x = -1

6𝑥+7
33. 5 ≥ > 2 con x = 2
2

𝑦+1 𝑦 2𝑦−1
34. −3 > 1+ con y = 2
4 6

35. En vigas continuas o apoyadas de ancho b y canto d, la distancia libre entre

los elementos rigidizadores laterales, 1, debe cumplir la doble condición:
𝑏2
1 < 60𝑏 ; 1 < 250 𝑑

En voladizos se debe cumplir

𝑏2
1 < 25𝑏 ; 1 < 100 𝑑

Calcular las condiciones cuando b=0.5 m, d=1.2 m

36. La desigualdad de Schwarz, se expresa como:

𝑥1 𝑦1 + 𝑥2 𝑦2 ≤ √𝑥12 + 𝑥22 √𝑦12 + 𝑦22

Verifique la desigualdad para x1 = 5, x2 = -2 y y1 = 3, y2 = 4

 37. Las cimentaciones se clasifican en dos clases, tal como se muestra a
continuación
Profundas superficiales

Determinar:

a. Tipo de cimentación cuya profundidad D es de 10 m y ancho B de 50 cm.

b. Tipo de cimentación cuya profundidad es de 50 cm y ancho de 2.5 m.
3.7 Operadores lógicos

Los operadores lógicos son:

▪ & (AND): Si dos expresiones x e y son verdaderas, entonces el resultado es

verdadero, por ende, su(s) salida(s) van a ser 1 y en caso contrario será cero,
es decir, falso.

▪ | (OR): Si alguna de las dos expresiones x e y es verdadera, entonces el

resultado es verdadero, por ende, su(s) salida(s) van a ser 1.

▪ ~ (NOT): Dada una expresión x, el operador cambia el estado de verdadero

a falso o viceversa.
Tablas de verdad de los operadores lógicos:

Operador lógico AND, &

x y x&y
1 1 1
1 0 0
0 1 0
0 0 0

Operador lógico OR, |

x y x|y
1 1 1
1 0 1
 0 1 1
0 0 0

Operador lógico NOT, ~

x ~x
1 0
0 1

Nota 7. Otras funciones equivalentes a los operadores lógicos que tiene

Matlab a disposición son:
▪ and(x,y) equivale a x&y
▪ or(x,y) equivale a x|y
▪ not(x) equivale a ~x

Ejemplos de operadores lógicos.

( 1 < 2 ) & ( 2 < 3) ( 1 < 2) & ( 2 < 1 )

ans = ans =
1 0
Verdadero Falso

Para la operación OR

>> ( 1 < 2 ) | ( 2 < 1 )

ans =
1
Verdadero

Para la operación NOT

>> ~ ( 2 < 1)
ans =
1
Verdadero

>> and(1==1,2==3)
ans =
0
Falso

>> 1==1 & 2==2

ans =
 1
Verdadero

>> 1==1 | 2==1

ans =
1
Verdadero

>> or(1==1, 2==1)

ans =
1
Verdadero

Para la operación OR exclusiva (XOR)

La puerta XOR es una puerta lógica digital que implementa el o exclusivo; es decir,
una salida verdadera (1) resulta si una, y solo una de las entradas a la puerta es
verdadera. Si ambas entradas son falsas (0) o ambas son verdaderas, resulta en
una salida falsa.

La XOR representa la función de la desigualdad, es decir, la salida es verdadera si

las entradas no son iguales, de otro modo el resultado es falso. Una manera de
recordar XOR es "uno o el otro, pero no ambos".

INPUT OUTPUT
A B A XOR B
0 0 0
0 1 1
1 0 1
1 1 0

Ejemplos

1. Revise con detenimiento los ejercicios resueltos:

xor(1==1,2==2) xor(1==1,2==2) xor(18<=17,1==2) xor(18<=17,2==2)

ans = ans = ans = ans =

0 0 0 1

Falso Falso Falso Verdadero

 2. Dados los vectores x = [5 0 8 3 9], y= [2 8 7 0 1], realizar las operaciones
indicados en el cuadro a la izquierda:

x = [5 0 8 3 9]; y= [2 8 7 0 1];

n = x & y n = 1 0 1 0 1
m = ~(y | x) m = 0 0 0 0 0
p = xor(x,y) p = 0 1 0 1 0
x<y ans = 0 1 0 0 0
x>y ans = 1 0 1 1 1
x<5 ans = 0 1 0 1 0
q = x((x<y) & (x<3)) q = 0

3.8 Práctica 3

Desarrolle las siguientes operaciones lógicas y analice los resultados.

1. and(2==3,4==3) 2. and(5==5,3==3)
3. and(6>7,3>2) 4. ~3<~5
5. and(6>7,3<2) 6. and(5<8,8<10)
7. or(2==3,4==3) 8. or(5==5,3==3)
9. or(6>7,3>2) 10. or(5<8,8<10)
11. xor(2==3,4==3) 12. xor(5==5,3==3)
13. xor(6>7,3>2) 14. xor(5<8,8<10)
15. ~(6<5) 16. ~(9>7)
17. ~(8<=10) 18. ~(15>=11)
19. 8>5&6>9 20. 1<=2&2>6
21. 3>=8 | 1>=1 22. 20==8 | 6<=14
23. ~5 24. ~(-6)
24. (9>=5 | 3>=1)&(5>=1 | 4>=2)
25. (~1<4 & 5>=4) | (7>=10 & 4<=7)
26. ~(4<=7 & 5>=2)
27. (~4>=15 | 13>=5)&( ~5>=4 | 7>=9)
28. (11<~7 & ~5>=-5) | (-10<=19 & -8<=7) &(8>~5)
 29. (9>7 & 3>=1) | (17<=19 & 6<=5) &(2>4)
30. (-1<~2 & 10<-2) | (1>3 & ~6>7)
31. t=sin(2) < 0.8 & log(2) > 0.5
32. x=cos(pi)> 3.11 | sqrt(10) > 3
33. M=[true,false,true,true, true, false]
34. x=1:6; y=[1,0,2,5,9,4]; z = x > y
35. x=1:6; y=[1,0,2,5,9,4]; emplear la instrucción find(x>y)
36. Se tiene los vectores x = [-6,-4,0,3,7] y z = [-15,12,8,5,9], calcular las
siguientes operaciones:
a. z = (x < y) z = (x > y) z = (x ~= y)
b. z = (x = = y) z = (x > 3) z = ~y > x
c. z= x&y z= x|y z = xor(x,y)
d. z=~(y < x) z = (2x > -3y) z = -2x&y
e. z= 5(x|y) z = (x2 > y2)
Además de los operadores booleanos citados, existen varias funciones lógicas
integradas útiles, tales como:

Instrucción Descripción

any Verdadero si cualquier elemento de un vector es verdadero.

all Verdadero si todos los elementos de un vector son

verdaderos.

exist Verdadero si el argumento existe.

isempty Verdadero true para una matriz vacía.

isinf Verdadero para todos los elementos infinitos de una matriz.

isnan Devuelve una matriz A del mismo tamaño que contiene 1

lógico (true) donde los elementos de A son NaNs y 0 lógico
(false) donde no lo son.

find Encuentra índices de elementos no nulos de una matriz.

3.10 Formatos de salida

La función format se encarga de cambiar la precisión numérica de los datos de

salida, y se emplea para hacer cambios de visualización de salida de cálculos
numéricos según el grado de precisión deseado. Sus variantes son:
 Instrucción Descripción
format Por defecto. Lo mismo que short.

format short Formato de punto fijo con 4 dígitos decimales. El formato

de visualización numérica oscila entre 0.001 y 1000.

format long Formato de punto fijo con 15 dígitos decimales.

format short e Formato de punto flotante con 4 dígitos decimales

expresados en notación científica.

format long e Formato de punto flotante con 15 dígitos decimales

expresados en notación científica.
format short g Formato para los mejores 5 dígitos fijos o punto flotante.

format long g Formato para los mejores 15 dígitos fijos o punto flotante.

format hex Formato hexadecimal.

format bank Formato fijo para dólares y céntimos en la que se toman 2
dígitos decimales.

format rat Formato que transforma una operación a fraccionario. Se

emplea cuando se operan racionales o reales en general
y se desea que la respuesta sea fraccionaria.

format compact Elimina las líneas vacías para permitir mostrar mayor
información en la pantalla.

format loose Es lo opuesto a compact, añade líneas vacías.

Los símbolos +, - y espacio en blanco se emiten para

elementos positivos, negativos y de valor cero. No se
format + o -
tienen en cuenta las partes imaginarias.

Con esta información, el usuario puede controlar el número de decimales que desea
que aparezcan en pantalla. Para emplear un determinado número de decimales
Matlab establece por defecto cuatros dígitos.

Por ejemplo. Sea el valor de la constante pi:

Estándar

>> pi
ans =
3.1416

Forma extensa
>> format long; pi
ans =
3.141592653589793

Forma extensa representada con exponente:

>> format long e;pi

ans =
3.141592653589793e+000

Forma corta con potencia:

>> format short e;pi

ans =
3.1416e+000

>> format long g;pi

ans =
3.14159265358979

>> format short g;pi

ans =
3.1416

Forma corta, en la que se representa mediante una relación comando número.

>> format rat; pi

ans =
355/113

>> format rat;

>> 2/7+1/3-3/11
ans =
80/231

Para representar números racionales en diversos formatos:

>> format bank; pi

ans =
3.14

>> format compact; pi

ans =
3.14

>> format hex; pi

ans =
400921fb54442d18

>> format loose; pi

ans =
400921fb54442d18

Recapitulando

>> format long >> format short e >> format long e

>> x=1/13 >> x=1/13 >> x=1/13
x= x= x=
0.076923076923077 7.6923e-02 7.692307692307693e-02

>> format hex >> format bank >> format rat

>> x=1/13 >> x=1/13 >> x=1/13
x= x= x=
3fb3b13b13b13b14 0.08 1/13

Para establecer cifras decimales con la precisión deseada, se emplea el comando

vpa, la sintaxis general es:

R = vpa (A)
R = vpa (A, d)

Calcular la relación de 1/3 con 5 dígitos significativos:

>> vpa (1/3, 5)

ans =
0.33333

Sin restricción alguna,

>> vpa(1/3)
ans =
0.33333333333333333333333333333333

Aunque existen variantes

Para este caso por defecto muestra 31 cifras:

>> vpa pi
ans =
3.1415926535897932384626433832795

Para establecer la precisión se tiene:

 >> vpa 'pi' 5 >> vpa 'pi' 7
ans = ans =
3.1416 3.141593

Otra forma de mostrar resultados exactos es con digits(n).

>> digits(25)
>> vpa(pi)
ans =
3.141592653589793238462643

>> digits(15); vpa (pi)

ans =
3.14159265358979

Para el caso de números enteros se tiene:

>> digits(24); vpa (5^35)

ans =
2910383045673370270040066.0

>> digits(30); vpa (5^35)

ans =
2910383045673370270040064.0

2 7 5
Sea la suma y resta de números racionales 3 + 11 − 9. Se puede evaluar empleando
la forma simbólica (sym) y luego situar la respuesta con una precisión de 18 dígitos.

>>sym(2/3+7/11+5/9)
ans =
184/99

>>digits(18)
vpa('2/3+7/11+5/9')
ans =
1.85858585858585859

Funciones de redondeo

Redondear un número significa reducir la cantidad de cifras manteniendo un

resultado más o menos exacto. Matlab posee varias funciones que permiten el
redondeo de un número. Como son:

Round(x). Redondea el valor de la variable x al entero más próximo.

>> round (-7.236587401215489301254789)

ans =
-7

>> round(8.25489)
ans =
8

Fix(x). Redondea el valor de x al entero más cercano por defecto.

>> fix (6.23510135479555130218738)

ans =
6

Cuando se declara una variable, esta se coloca en la función, así:

>> x=6.23510135479555130218738;
>> fix(x)
ans =
6
Floor(x). Redondea el valor de x al entero más cercano por defecto si es positivo, y
por exceso si es negativo.

>> floor (3.6589710124789654135)

ans =
3

>> floor (-3.6589710124789654135)

ans =
-4

Ceil(x). Redondea el valor de x al entero más próximo por exceso si es positivo, si

es negativo se aproxima por defecto.

>> ceil (9.3254987512568475555)

ans =
10

>> ceil (-9.3254987512568475555)

ans =
-9

En resumen, los comandos anteriores son:

Round(u) fix(u) floor(u) ceil(u)

>> round(u) >> fix(u) >> floor(u) >> ceil(u)
ans = ans = ans = ans =
3 2 2 3
Ejemplos

La distancia recorrida en línea recta por un vehículo es de 17.5879 metros, en un

tiempo de 12 segundos. Determine la rapidez empleando las diferentes funciones
de redondeo básicas, luego ajústelas a 5 dígitos.
𝑥
La fórmula de rapidez es: 𝑢 = al digitarse en Matlab, se tiene:
𝑡

>> x=17.5879;
>> t=12;
>> u=x/t;
>>vpa (u) % muestra el cálculo sin redondear
ans =

1.4656583333333335073

% luego use las demás funciones de redondeo, recuerde use vpa (u) para restablecer la respuesta.

>> round(u)
ans =

1
>> vpa(u)
ans =

1.4656583333333335073
>> fix(u)
ans =
1

>> vpa(u)
ans =
1.4656583333333335073

>> floor(u)
ans =
1

>> vpa(u)
ans =
1.4656583333333335073

>> ceil(u)
ans =
2

Para ajuste a 5 dígitos

 >> round(5);u >> fix(5);u >> floor(5);u >> ceil(5);u
u= u= u= u=
1.4657 1.4657 1.4657 1.4657

3.11 Práctica 4

1. Aproximar el número áureo empleando los diferentes formatos que dispone

Matlab y ordene en una tabla. Indague que es este número, en donde se
encuentra en la naturaleza y para que se emplea.

2. Responda a las siguientes preguntas:

a. ¿Es fix (4.71) igual que floor (4.71)?
b. ¿Es fix (5.4) igual que fix (-5.4)?
c. ¿Es fix (9.2) igual que floor (9.2)?
d. ¿Es fix (-7.2) igual que floor (-7.2)?
e. ¿Es fix (-2.15) igual que ceil (-2.15)?

En los siguientes ejercicios, aproximar a por lo menos cuatro cifras significativas.

3𝑥+5
3. Dada la función 𝑓(𝑥) = 𝑥 2 +3𝑥−5 . Calcular f(2.514455) y f(-1.1257), cuya
respuesta sea de 7 dígitos decimales.

4. La función polinomial h(t) representa la altura (en pies) de un segmento de

montaña rusa en un tiempo t:

ℎ(𝑡) = 0.00151𝑡 3 − 0.12587𝑡 2 + 3.698412𝑡 + 12.4587

a. Determine la altura para un tiempo t = 1.2554. Exprese la respuesta en forma

fraccionaria.
b. Determine la altura para un tiempo t = 5.012478. Exprese la respuesta con 9
dígitos decimales.
c. Determine la altura para un tiempo t = 12.895471. Exprese la respuesta en
notación científica.

5. La relación entre lecturas de escalas Celsius y Fahrenheit de temperatura se

representa mediante la siguiente ecuación:

5
𝐶 = (𝐹 − 32)
9
 a. Determine el equivalente en grados Celsius cuando la temperatura es de
58.26, 69.75, 85.691 y 125.01258 grados Fahrenheit. Exprese la
respuesta en formato compacto.
b. Determine el equivalente en grados Fahrenheit cuando la temperatura es
de 15, 30, 67, 89, 98 y 110 grados centígrados. Exprese la respuesta en
formato bancario.

6. El momento de inercia de un cilindro homogéneo alrededor de un eje

perpendicular a su eje de simetría que pasa por su centro de masas es:

1 2
ℎ2
𝐼 = 𝑀 (𝑟 + )
4 3

Donde M es la masa del cuerpo, r el radio y h la longitud. Exprese la respuesta

empleando un formato que aproxime por defecto.

a. Si la masa de un cuerpo cilíndrico es de 40.0851 Kg, determine su momento

de inercia si la longitud h es de 3.415 m y el diámetro es de 2.42 m.
b. Si el momento de inercia de un cuerpo cilíndrico es de 12.9103 Kgm 2,
determine la altura de este, si su radio es de 50.7 cm y la masa de 60 Kg.

7. La altura h (en pies) de la base de una nube se puede estimar usando la

ecuación:
ℎ = 227(𝑇𝑠 − 𝑇𝑝𝑟 )

Donde Ts es la temperatura del suelo y Tpr es la temperatura del punto de rocío.

Exprese la respuesta con máximo 3 decimales.

a. Si la temperatura es 72.521 oF y el punto de rocío es 55,854 oF, encuentre la

altura de la base de la nube.
b. Si el punto de rocío es 69.125 oF y la base de la nube está a 3500.213 pies,
estime la temperatura del suelo.

8. La presión a una distancia r del centro de la Tierra es:

𝑅
3𝐺𝑀2 𝑟 3𝐺𝑀2 2
𝑃𝑜 − 𝑃 = − ∫ 𝑑𝑟 = − (𝑅 − 𝑟 2 )
4𝜋𝑅 3 𝑅 3 8𝜋𝑅 6
𝑟

Donde p0=1.0131x105 Pa es la presión cuando r = R, es decir, la presión

atmosférica. En el centro de la Tierra r = 0, por lo que la presión P es igual a:

3𝐺𝑀2
𝑃 = 𝑃𝑜 +
8𝜋𝑅 4
 Obtenga el valor de P con R=6.37·106 m, M=5.98·1024 kg y G=6.67·10-11
Nm2/kg2. Aplique el o los métodos de aproximación que garantice el menor error
en el cálculo final.

9. La ecuación del periodo de un cuerpo que describe un Movimiento Armónico

Simple (MAS) a nivel gravitacional está descrita por la ecuación:
𝑅3
𝑇 = 2𝜋√
𝐺𝑀

Obtenga el valor de T con R=6.37·106 m, M=5.98·1024 kg y G=6.67·10-11

Nm2/kg2. Aplique el o los métodos de aproximación que garantice el menor error
en el cálculo final. Exprese la respuesta en horas, minutos y segundos.

10. El radio de un satélite geoestacionario se obtiene al suponer un periodo de órbita

de 24 horas. Calcule el radio y expréselo en m, Km, cm, mm. Ayuda. Tome la
ecuación del punto anterior y despeje R, teniendo en cuenta las unidades y
conversión de las mismas.

11. La circunferencia de una elipse puede ser aproximada por medio de la siguiente
ecuación:
𝐶 = 𝜋[3(𝑎 + 𝑏) − √(3𝑎 + 𝑏)(𝑎 + 3𝑏)]

Calcular la circunferencia C, cuando a = 12.125 cm y b = 15.784 cm. Exprese

la respuesta con 5 decimales de exactitud.

12. Un cuerpo esférico cayendo bajo su propio peso en un fluido tiene como
velocidad límite,
2𝑟 2 (𝜌 − 𝜌𝑓 )𝑔
𝑣𝐿 =
9𝜂

Donde 𝜌 es la densidad del cuerpo, 𝜌𝑓 es la densidad del fluido, 𝜂 es la

viscosidad del fluido y r es el radio del cuerpo.

Sea un glóbulo rojo, cuyo radio es de aproximadamente 2 µm y su densidad es

de 1.312x103 Kg/m3, la densidad del plasma sanguíneo es de 1.06 x103 Kg/m3.
Exprese la respuesta con 3 decimales de exactitud.

a. Cuál es la velocidad límite de los glóbulos rojos, si la viscosidad del

plasma es de 1.528x10-3 Ns/m2.
b. Que radio tendría un glóbulo rojo, si la velocidad límite es 0.91 veces la
velocidad del punto anterior.

13. El momento de inercia de la Tierra respecto del eje Z es;

 I= 0.345·MR2

El momento de inercia de una distribución esférica y homogénea de masa M es,

I = 0.40MR2

Siendo M=5.98·1024 kg la masa de la Tierra, y R=6.371·106 m el radio de la

misma. Calcular los momentos aplicando el o los métodos de aproximación que
garantice el menor error en el cálculo final.

14. La salida de potencia mecánica P de un músculo en contracción está dada por

la ecuación:
(1 − 𝑢)
𝑃 = 𝑘𝑢
𝑘+𝑢

Donde u = v/vmax, es la relación de velocidades, siendo v y vmax las velocidades

de acortamiento muscular estándar y máxima respectivamente. La constante k
varía entre 0.15 y 0.25. Determine P cuando la relación de velocidades es de
0.75 y k = 0.22; 0.67 y 0.18. Exprese la respuesta con 3 decimales de exactitud.

15. Debajo de la base de una nube, la temperatura del aire T (medida en oF) a una
altura h (en pies) se expresa de manera aproximada con la ecuación:

5.51
𝑇 = 𝑇𝑜 − ℎ ( )
1000

Donde 𝑇𝑜 es la temperatura al nivel del suelo, con dos 3 decimales de exactitud.

a. Determine la temperatura del aire a una altura de 1.25 millas si la

temperatura del suelo es de 72.2368 oF.
b. ¿A qué altura se alcanza la temperatura de congelación?

16. El costo en dólares f(p) de eliminar la contaminación de un río se aproxima a la

fórmula:
185000𝑝
𝑓(𝑝) = 100−𝑝 0 ≤ 𝑝 ≤ 100

Evalúe el costo de eliminar los porcentajes de contaminación de 35.47% y

72.0458%.

17. Los momentos inerciales de un cilindro circular de radio r y longitud L están

representados por las siguientes ecuaciones:
1 1
𝐼𝑥 = 2 𝑚𝑟 2 ; 𝐼𝑦 = 𝐼𝑧 = 12 𝑚(3𝑟 2 + 𝐿2 )
a. Calcular los momentos para un cilindro de masa m = 850.12 g, con un
diámetro de 4.2 cm y una longitud de 75.15 cm. Exprese la respuesta en
forma fraccionaria.

b. El momento inercial en la dirección z es de 16.687x10 -3 kgm2, determine el

momento 𝐼𝑥 cuando la masa de una varilla cilíndrica es de 0.8 kg y la longitud
es de 80 cm. Exprese la respuesta con 4 decimales de exactitud.