Practicas de MatLab
Practicas de MatLab
Practicas de MatLab
⎡6 9 5 1⎤ ⎡4 8⎤
⎢8 7 2 ⎥
3⎥ ⎢3 7 ⎥⎥
001. Dadas las matrices C = ⎢ y D=⎢ , realizar en MATLAB las siguientes
⎢1 3 4 4⎥ ⎢2 3⎥
⎢ ⎥ ⎢ ⎥
⎣5 2 8 2⎦ ⎣5 1⎦
operaciones:
a. Crear la matriz E1 con las dos columnas centrales de C usando el operador dos puntos.
b. Crear la matriz E2 con las filas 1 y 2, y las columnas 2 y 3 de C usando el operador dos puntos.
c. Crear la matriz E3 colocando E1 y D lado por lado.
d. Encuentre el producto de C2,4 y D1,2
⎡7 2 ⎤ ⎡2 3⎤
003. Dada las matrices A = ⎢ ⎥ y B=⎢ ⎥ , calcular lo siguiente:
⎣2 1 ⎦ ⎣− 4 − 5⎦
a. A+B
b. A-B
c. 2*B
d. A/4
e. A.*B
f. B.*A
g. A*B
h. B*A
i. A.^2
j. A^2
k. A.^B
l. A./B
[
005. Dados los vectores x = [3 2 6 8] y y = 2 − 3 2 0 , ]
a. Obtenga la suma de los elementos de x y y
b. Obtenga un vector z cuyos componentes sean los elementos del vector x elevados a la potencia
especificada por cada elemento correspondiente en el vector y.
c. Dividir cada elemento de y para cada elemento correspondiente de x.
d. Obtenga un vector z cuyos componentes sean los elementos del vector x multiplicados por cada
elemento correspondiente del vector y.
e. Ejecutar la operación x y − z
T
007. Construir una matriz A de 2x3 cuyas filas son los 6 primeros impares consecutivos.
a. Anular el elemento (2,3)
b. Obtener la matriz B = A’
c. Construir una matriz C, formada por la matriz B y la matriz identidad de orden 3 adosada a su derecha
d. Construir una matriz D extrayendo las columnas impares de la matriz C
e. Construir una matriz E formada por la intersección de las dos primeras filas de C y sus columnas
tercera y quinta
f. Construir una matriz F formada por la intersección de las dos primeras filas y las tres últimas
columnas de la matriz C
g. Construir una matriz diagonal G tal que los elementos de su diagonal principal son los mismos que los
de la diagonal principal de D
h. Calcular el orden de la matriz C
a. A − [ x ' y' ]
b. [ x ; y' ]
c. [ x ; y]
d. A - 3
⎡ 1 0⎤ ⎡ 1 2 3⎤ ⎡ − 1 1⎤
011. Ingresar las matrices A = ⎢ ⎥ , B=⎢ ⎥ , C=⎢ ⎥y
⎣ − 1 2⎦ ⎣ − 1 0 4⎦ ⎣ 5 5⎦
⎡ − 1 1 3⎤
D=⎢ ⎥ , y evaluar cada una de las siguientes expresiones. Explicar cualquier mensaje de
⎣ 2 2 2⎦
error:
a. 2*A+3*C
b. A-4*D
c. B^2
d. B.^2
e. A*B
f. B*A
g. C*D
h. C.*D
i. A*B+D
⎡1 2 3 3⎤ ⎡1 2 3 3⎤
⎢2 2 0 ⎥
0⎥ ⎢− 1 2 0 0⎥⎥
014. Crear las matrices X = ⎢ , Y=⎢ . Usar operaciones entre
⎢5 6 3 1⎥ ⎢5 6 3 1⎥
⎢ ⎥ ⎢ ⎥
⎣0 2 0 5⎦ ⎣0 2 0 5⎦
arreglos para responder cada uno de los literales:
017. En una sola orden de MATLAB crear una matriz 3 x 5 cuyo único elemento sea el 7.
018. Con una sola orden de MATLAB crear una matriz aleatoria 4 x 4 de números reales entre -5 y 5.
020. Sea x = (0:pi/2:2*pi). Con una sola orden de MATLAB crear una matriz cuya primera fila es x, su
segunda fila es el seno de cada elemento de x y cuya tercera fila el coseno de cada elemento de x.
021. Definir un vector a formado por los cuatro primeros números impares y otro b formado por los
cuatro primeros números pares de varias formas distintas. Emplearlos para construir la matriz
⎡2 4 6 8 ⎤
⎢ 6 12 18 24⎥
A=⎢ ⎥
⎢10 20 30 40⎥
⎢ ⎥
⎣14 28 42 56 ⎦
022. Construya una matriz n x n, C = (cij)
a. Con cij = i·j;
b. Con cij = cos(i·j);
024. En una sola instrucción, cambiar todos los valores de la diagonal de una matriz cuadrada a cero.
025. En una sola instrucción, sustituir todos los valores de la diagonal de una matriz cuadrada por los
elementos de un vector dado.
026. Ordenar los elementos de una matriz del menor al mayor manteniendo su forma (indicación:
emplear la orden reshape)
027. De tres formas distintas (cada una en una sola instrucción), averiguar el número de elementos de una
matriz, de forma que al final tengamos un número.
001. Trabaje en los resultados de las siguientes operaciones, antes de verificar con MATLAB,
dados los vectores a = [-1 0 3]; y b = [0 3 1];
a. ˜a
b. a&b
c. a|b
d. xor(a, b)
e. a>0&b>0
f. a>0|b>0
g. ˜a>0
h. a + (˜ b)
i. a>˜b
j. ˜a>b
k. ˜ (a > b)
002. Determine los valores de verdad de las siguientes operaciones antes de verificar com MATLAB,
a. 1 & -1
b. 13 & ˜(-6)
c. 0 < -2|0
d. ˜[1 0 2] * 3
e. 0 <= 0.2 <= 0.4
f. 5>4>3
g. 2>3&1
003. Dados a = [1 0 2] y b = [0 2 2] determine los valores de las siguientes expresiones. Verifique sus
respuestas con MATLAB.
a. a ˜= b
b. a<b
c. a<b<a
d. a<b<b
e. a | (˜a)
f. b & (˜b)
g. a(˜(˜b))
h. a = b == a (determine el valor final de a)
[-4 0 5 -3 0 3 7 -1 6]
006. Crear una matriz A de 6×6 cuyos elementos son números aleatorios comprendidos -4 y 4. En una
sola instrucción, poner a cero todos los elementos negativos de una matriz.
007. Crear una matriz A de 6×6 cuyos elementos son números aleatorios comprendidos -3 y 3. En una
sola instrucción, poner a cero todos los elementos de la matriz que estén entre -1 y 1. (La conjunción
lógica es &).
008. Una compañía ofrece 7 niveles de salario (dólares): 12000, 15000, 18000, 24000, 35000, 50000,
70000. El número de empleados en cada nivel son, respectivamente: 3000, 2500, 1500, 1000, 400, 100 y
25. Escriba sentencias en la ventana de comandos para calcular lo siguiente:
a. El nivel de salario promedio. (Use el comando mean) (Rta: 32000)
b. El número de empleados sobre y bajo este salario medio. Use vectores lógicos para encontrar qué
niveles de salario están sobre y bajo el nivel promedio. Multiplique estos vectores lógicos elemento por
elemento por el vector empleados, y sume el resultado. (Rta: 525 sobre, 8000 bajo)
GRÁFICAS
Gráficos 2D.
sin(50t ) sin(50t )
b. x ( t ) = cos(t + π) , y( t ) = sin (t + π) , en el intervalo -10π ≤ t ≤ 10π
t t
x (x − 2)
c. f ( x ) = , en el intervalo -2≤x≤2
(x + 1)(x − 2)
⎧x 2 si x < 0
d. f ( x ) = ⎨ , en el intervalo -2≤x≤2
⎩− 1 si x ≥ 0
⎧− x si x < 0
⎪
e. f ( x ) = ⎨1 si 0 < x < 2 , en el intervalo -4≤x≤4
⎪x 2 si x > 2
⎩
⎧ 1− x si x < -1
⎪⎪
f. f ( x ) = ⎨1 − x si - 1 < x < 1 , en el intervalo -2≤x≤2
2
⎪
⎪⎩ x - 1 si x > 1
002. Dibujar las curvas paramétricas siguientes, usando los comandos plot y comet; dibujar además los
vectores velocidad utilizando el comando quiver.
003. Dibujar las gráficas de las siguientes funciones en coordenadas polares (Usar el comando
polar(theta,rho) ):
a. r = 7 − 7 sin (θ ) ; θ∈[-π,π]
b. r = sin (6θ ) ; θ∈[-π,π]
c. r = 5 cos(2θ ) ; θ∈[-π,π]
d. r = sin (6θ ) ; θ∈[-π,π]
004. Dados los problemas del ejercicio 4, utilizar el comando pol2cart(theta,rho) para cambiar las
coordenadas polares a coordenadas cartesianas. Usar después el comando plot para obtener las gráficas
en las nuevas coordenadas.
Gráficos 3D
001. Utilizando los comandos plot3 y comet3, dibujar las curvas que son imagen de las trayectorias
dadas. Represéntese algunos vectores velocidad utilizando el comando quiver3.
b. ( )
c(t ) = 4e t ,6t 4 , cos( t ) , donde 0 ≤ t ≤ 2π
c. ( )
c(t ) = cos 2 ( t ),3t − t 2 , t , donde 0 ≤ t ≤ 2π
d. ( ( ) )
c(t ) = sen e t , t ,4 − t 3 , donde 0 ≤ t ≤ 2
e. ( )
c(t ) = 4e t ,6t 4 , cos( t ) , donde 0 ≤ t ≤ 2π
⎧v = 1m/s
⎛ ⎛ vt ⎞ ⎛ vt ⎞ ⎞ ⎪
f. c(t ) = ⎜⎜ vt − r sin ⎜ ⎟, R − r cos⎜ ⎟ ⎟⎟, donde ⎨R = 1m
⎝ ⎝R⎠ ⎝ R ⎠⎠ ⎪r = {0,0.5,1,2}
⎩
002. Graficar las superficies y las curvas de nivel para cada una de las siguientes funciones (Use los
comandos mesh y contour),
⎧ xy
⎪ 2 (x, y) ≠ (0,0 )
(
c. f ( x , y) = ⎨ x + y 2 ) 1/ 2
para -1≤x≤1 y -1≤y≤1
⎪0 (x, y) = (0,0 )
⎩
d. z = x 2 + y 2 , - 2 ≤ x ≤ 2, - 2 ≤ y ≤ 2
e. z = − xy , - 2 ≤ x ≤ 2, - 2 ≤ y ≤ 2
f. z = x + y + 1 , - 2 ≤ x ≤ 2, - 2 ≤ y ≤ 2
2 2
x
g. z = , - 2 ≤ x ≤ 2, - 2 ≤ y ≤ 2
y
entre -1 y 1.
b. Graficar la curva: 3y + y − x = 5
3 3
c. Graficar la curva: x − y = x + y
2 2
( 2
)
2 2
(lemniscata)
d. Graficar la curva de nivel de la función f ( x , y) = y ln (x ) + x ln (y ) que contenga el punto (1,1).
x y2 z2
a. = +
4 4 9
2 2
y z x2
b. + = 1+
9 4 16
c. 4 x − 3y + 2z = 0
2 2 2
x 2 y2 z2
d. + + =1
9 2 9
005. Utilizando el comando sphere , dibujar una esfera de radio 2, y de centro (1,-3,5)
a. R ( t ) = t , t∈[-1,1]
b. R ( t ) = t , t∈[-1,1]
2
007. Utilizando los comandos sphere y cylinder dibujar la intersección entre una esfera de radio 2 con
centro en el origen; y tres cilindros de radio 1, con ejes de simetría a lo largo de los ejes x , y , z, y de
longitud 2.
x+y
008. Dada la función V(x, y ) = ,
x 2 + y2
a. Utilizando los comandos grad y quiver, esbozar el campo gradiente − gradV para el intervalo
-1≤x≤1 , -1≤y≤1.
PROGRAMACIÓN
001. Escribir un archivo m que permita determinar si un número entero dado es par o impar. (Utilizar el
operador mod)
002. Utilizando el operador relacional > (mayor que), escribir un archivo.m que permita definir si un
número a es mayor que un número b. El programa debe admitir ingresar los números a y b, e imprimir el
resultado a es mayor que b, o a es menor que b, o a es igual a b.
004. Crear un archivo.m que calcule las raíces de la ecuación: a·x2 + b·x + c = 0. Teniendo en cuenta los
siguientes casos:
− b ± b 2 − 4·a·c
xi =
2a
La expresión d=b2-4ac se denomina discriminante.
Indicar con literales adecuados los datos a introducir, así como los resultados obtenidos.
005. Escribir un programa que permita imprimir sobre la ventana de comandos los números impares del
1 al 20.
a. Utilizar la sentencia for…end
b. Utilizar la sentencia while…end
006. Escribir un archivo m que de como resultado la suma de los 100 primeros números naturales.
a. Utilizar la sentencia for…end
b. Utilizar la sentencia while … end
007. Crear un archivo m que de como resultado la suma de los números pares comprendidos entre dos
números enteros a y b ingresados por el usuario.
Para empezar la secuencia, se deben asignar valores de partida a1 y a2. Escribir un archivo de función
para calcular la secuencia, asignando como parámetros de entrada a1 y a2. No obstante, la primera vez
que | an | > 1000, el cómputo de valores deberá finalizar.
a. Utilizar la sentencia de control for…end
b. Utilizar la sentencia de control while…end
009. Implementar un programa que permita evaluar el factorial de un número entero positivo. Por
ejemplo
Factorial(5) = 5*4*3*2*1 = 120
Factorial(3) = 3*2*1 = 6
010. Implementar un programa que permita determinar si un número entero positivo es primo.
011. Implementar un programa que de como resultado un vector formado por los números primos
menores que un número dado n (positivo y entero).
012. Vamos a suponer que nuestro problema es sumar los términos de una serie infinita, ∑a n , en el
computador. No obstante, es posible sumar únicamente un número finito de términos para aproximar la
−6
suma exacta de la serie. Especificamos el criterio a n < 10 como test de convergencia. Escriba
archivos de función para calcular las siguientes series:
a. ln (2) = 1 −
1 1 1 1
+ − + − ... = ∑
∞
(− 1) n +1
2 3 4 5 n =1 n
b.
π 1 1 1
= 1 − + − + ... = ∑
∞
(− 1) n +1
4 3 5 7 n =1 2n − 1
∞ ∞
1 1 1 1
c. e = ∑ a n = ∑ = + + + ...
n =0 n = 0 n! 0! 1! 2!
013. Escribir un programa que lea una cadena de n caracteres e imprima el resultado que se obtiene cada
vez que se realice una rotación de un carácter a la derecha sobre dicha cadena. El proceso finalizará
cuando se haya obtenido nuevamente la cadena de caracteres original. Por ejemplo:
HOLA AHOL LAHO OLAH HOLA
014. Escribir un programa que permita ingresar una palabra, y posteriormente imprimir la palabra con las
letras en orden inverso. Por ejemplo:
Se ingresa: INGENIERIA
Se imprime: AIREINEGNI
015. Escribir un programa que dibuje un triángulo de n filas, empleando el caracter asterisco. Crear el
programa utilizando:
a. Sentencias for…end
b. Sentencias while…end
B* BNB *BN
NB *B*B NB
BNB*BN BN
NB*B* B NB
B*BNB * BN
*BNBNB * B
BNBNBNB*
NBNBNBNB
La caída de presión del flujo de un fluido en un tubo puede expresarse mediante la forma adimensional
del factor de fricción de Darcy, definida por la ecuación:
⎛1 ⎞⎛ L ⎞
∆P* = (P1 + ρgh ent ) − (P2 + ρgh sal ) = f ⎜ ρV 2 ⎟⎜ ⎟
⎝2 ⎠⎝ D ⎠
El factor de fricción de Darcy f, depende del régimen de flujo (laminar o turbulento). Para flujo
turbulento, Colebrook, encontró en experimentos en tubos revestidos con un material áspero con una
altura física (de la rugosidad) ε, que el factor de fricción f medido, se relaciona con el número de
Reynolds (Re) y la razón rugosidad a diámetro del tubo (ε/D) mediante la fórmula:
1 ⎛ε/D 2.51 ⎞
= −2 log10 ⎜ + ⎟ si Re ≥ 2300 (Flujo turbulento)
f ⎝ 3.7 Re f ⎠
Analíticamente se puede demostrar que para flujo laminar en un tubo circular, el factor de fricción de
Darcy, se relaciona con el número de Reynolds mediante la fórmula:
64
f = si Re<2300 (flujo laminar)
Re
a. Escribir un archivo de función para calcular el factor de fricción de Darcy, dados como parámetros de
entrada: el número de Reynolds (Re), la rugosidad de la tubería (ε) y el diámetro del tubo (D). Para
despejar f de la ecuación de Colebrook, utilice la relación recursiva:
−2
⎛ ⎛ε/D 2.51 ⎞⎟ ⎞⎟
f i +1 = ⎜ − 2 log10 ⎜ +
⎜ ⎜ 3.7 Re f ⎟ ⎟
⎝ ⎝ i ⎠⎠
Para empezar el cálculo iterativo, especificar dentro del programa un valor de partida para f (por ejemplo
f i +1 − f i
f = 0.001). Para finalizar las iteraciones utilizar el criterio de convergencia error = < 10 −6 o
fi
un número máximo de 1000 iteraciones (el criterio que se alcance primero).
d. Bajo las condiciones de flujo especificadas en el literal c, calcúlese el caudal para el cual la caída de
presión estática sea cero.
018. Cálculo del caudal en una tubería para una pérdida dada de presión.
En muchos de los casos de flujos en tuberías, el gasto volumétrico por la tubería no está controlado sino
que crece hasta alcanzar un nivel que concuerde con la caída de presión disponible. Por ejemplo al abrir
por completo el grifo en un lavamanos, el caudal aumenta hasta que la caída de presión en la tubería
iguala la diferencia de presión entre el suministro de agua y la presión atmosférica. El procedimiento
para calcular el gasto Q es el siguiente:
∆P * (Pent + ρgh ent ) − (Psal + ρgh sal )
1. Se calcula el cambio de carga hidráulica ∆h f = =
ρg ρg
g(∆h f )D 2
2. Suponiendo que el flujo es laminar (Re<2300), se calcula V a partir de la ecuación V =
32νL
y el número de Reynolds por la ecuación Re = VD / ν .
πD 2 V
3. Si Re<2300, el flujo es laminar y Q = .
4
⎛ 2g (∆h f )D 3 ⎞
4. Si el flujo no es laminar (Re>2300), se calcula f Re = ⎜ ⎟⎟ y el número de Reynolds
⎜ ν 2
L
⎝ ⎠
b. Una tubería de acero (ε = 1E(-5) m) comercial con diámetro de 8 in = 0.2032 m, conduce agua (ρ =
1000 kg/m3, ν = 1E(-6) m2/s) desde la salida de un tanque de almacenamiento (a una altura de 100 m) a
una distancia de 2 km (a una altura de 22 m) desde el tanque de almacenamiento. Calcule el gasto
volumétrico (Q) cuando la caída de presión estática Pent – Psal es cero.
CÁLCULO SIMBÓLICO
Álgebra simbólica
⎧3x + 4 y + 5z = 2
⎪
a. ⎨2 x − 3y + 7 z = −1
⎪x − 6 y + z = 3
⎩
⎧f ( x, y) = y − 4 x 2 + 3 = 0
⎪
b. ⎨ x2
⎪g ( x, y) = + y2 −1 = 0
⎩ 4
⎧f ( x , y) = y − sin( x ) = 0
c. ⎨
⎩g ( x , y) = y − x + 1 = 0
2
⎧⎪f ( x, y) = y − x 2 + 1 = 0
d. ⎨
⎪⎩g( x , y) = y − x 2 + 1 = 0
e. f ( x, y) = x − exp(xy) = 0 , despejar y.
f. ax + bx + c = 0 , despejar x.
2
Límites.
sin (x )
a. Lim
x →0 x
1 + cos(x )
b. Lim
x →−π x+π
2 −x
c. Lim x e
x →∞
1
d. Lim−
x →1 x −1
⎛1⎞
e. Lim+ sin ⎜ ⎟
x →0 ⎝x⎠
f. Lim
(x + y )2 − (x − y )2
( x , y )→ (0 , 0 ) xy
cos(xy ) − 1
g. Lim
( x , y )→ (0 , 0 ) x 2 y2
sin( xy)
h. Lim
( x , y )→(0 , 0 ) xy
e xy
i. Lim
( x , y )→ (0 , 0 ) x + 1
cos( x ) − 1 − ( x 2 / 2)
j. Lim
( x , y )→ (0 , 0 ) x 4 + y4
xyz
k. Lim
( x , y , z )→ ( 0 , 0 , 0 ) x 2 + y 2 + z 2
Derivación
d
a. (acrtg(sinh (x )))
dx
b.
d
dx
(
tanh −1 (sin (x )) )
d ⎛ 1 ⎛b ⎞⎞
c. ⎜⎜ arctg⎜ tg (x )⎟ ⎟⎟
dx ⎝ ab ⎝a ⎠⎠
d ⎛ x ⎛ x ⎞ ⎞⎟
d. ⎜ − arcsin ⎜ ⎟⎟
dx ⎜⎝ a 2 − x 2 ⎝ a ⎠⎠
d ⎛ ⎛ x − a ⎞⎞
e. ⎜⎜ (x − a ) 2ax − x 2 + a 2 arcsin⎜ ⎟ ⎟⎟
dx ⎝ ⎝ a ⎠⎠
∂ ⎛ x 2 + y2 ⎞
a. ⎜ ⎟
∂y ⎜⎝ x 2 − y 2 ⎟⎠
∂ 2 ⎛⎜ xy ⎞
⎟
b.
∂x∂y ⎜ x 2 + y 2 ⎟
⎝ ⎠
∂ ⎛
2
⎛ x ⎞⎞
c. ⎜ x arctan ⎜⎜ ⎟⎟ ⎟⎟
∂x 2 ⎜⎝ ⎝ y ⎠⎠
d.
∂3
∂x ∂y
2
(
A cos x 2 + y 2( ))
∂ ⎛
e.
⎛A
⎜⎜ A 2 + B 2 exp⎜ − x 2 − y 2
∂y ⎝
( )⎞⎟ ⎞⎟⎟
⎝B ⎠⎠
mM
003. Demostrar que el potencial gravitatorio de Newton V = −G satisface la ecuación
x 2 + y2 + z2
de Laplace:
∂2V ∂2V ∂2V
+ + =0
∂x 2 ∂y 2 ∂z 2
004. Escribir un archivo de función llamado div para calcular la divergencia de un campo vectorial en
forma simbólica. Como parámetro de entrada se especificará el campo vectorial con nomenclatura
matricial: F = [Fx , Fy , Fz ] . Donde F : ℜ 3 a ℜ 3 .
La divergencia de un campo vectorial en coordenadas cartesianas es:
Obsérvese que (∇ • F ) : ℜ 3 a ℜ
005. Escribir un archivo de función llamado rot para calcular el rotacional de un campo vectorial en
forma simbólica. Como parámetro de entrada se especificará el campo vectorial con nomenclatura
matricial: F = [Fx , Fy , Fz ] . Donde F : ℜ 3 a ℜ 3 .
El rotacional de un campo vectorial en coordenadas cartesianas es:
Obsérvese que (∇ × F ) : ℜ 3 a ℜ 3
006. Escribir un archivo de función llamado gradiente para calcular el gradiente de un campo escalar en
forma simbólica. Como parámetro de entrada se especificará el campo escalar f en coordenadas
cartesianas. Donde f : ℜ a ℜ .
3
⎡ ∂f ∂f ∂f ⎤
grad(f ) = ∇f = ⎢ , , ⎥
⎣ ∂y ∂z ∂z ⎦
∂f ˆ ∂f ˆ ∂f ˆ
≡ i+ j+ k
∂x ∂y ∂z
007. Utilizando los comandos div y rot de los problemas 4 y 5, calcular la divergencia y el rotacional de
los siguientes campos vectoriales:
a. F( x, y, z) = xi + yj + zk
b. F( x , y, z ) = yzi + xzj + xyk
yzi − xzj + xyk
c. F( x , y, z ) =
x 2 + y2 + z2
d. F( x, y, z) = r , siendo r = xˆi + yj + zkˆ ≡ [x, y, z ]
r
e. F( x , y, z ) = 3 , siendo r = xˆi + yj + zkˆ ≡ [x , y, z ] , r = (r • r )
1/ 2
r
008. Utilizando el comando gradiente del problema 6, calcular,
⎛1⎞
a. ∇⎜ , siendo r = x + y + z
2 2 2
n ⎟
⎝ ⎠
r
b. (
grad x 2 + y 2 + z 2 )
⎛ ⎞
mM
c. FG = −∇⎜ − G ⎟
⎜ 2 ⎟
x +y +z ⎠
2 2
⎝
⎛ 1 qQ ⎞
d. FE = −∇⎜ ⎟
⎜ 4πε 0 2 ⎟
x +y +z ⎠
2 2
⎝
(
e. ∇ x y + y z
2 2 2 2
)
009. Sea f ( x, y, z) = x 2 y 2 + y 2 z 2 . Verificar que ∇ × ∇f = 0
a. F • (∇ × G )
b. F • (∇f )
c. F × (∇f )
d. ∇ 2 f = ∇ • ∇f
Sumatorias
n
⎛1 1 ⎞
a. ∑ ⎜⎝ k − 1 + k ⎟⎠
k =1
∞
1
b. ∑n
n =1
2
n
c. ∑k
k =1
2
∞
d. ∑r
k =0
k
n
xk
e. ∑
k = 0 k!
; la función factorial no trabaja con entradas simbólicas, pero se puede usar sym(‘k!’) o la
∞
función Γ(x ) = ∫e
− t x −1
t dt llamada gamma en MATLAB, que satisface Γ(k + 1) = k!
0
∞
1
f. ∑ (z − k )
k = −∞
2
Integración
∫ x sin (x )dx
2
b.
∫e
−x2
e. dx
−∞
1 a
x
f. ∫∫
0 0 a 2 − y2
dydx
3 x 2 +1
g. ∫ ∫ xydydx
0 − x 2 +1
1 e2 x
h. ∫ ∫ x ln(y)dydx
0 ex
π/2
y cos(xy )dxdy
arcsin( y ) / y
i. ∫ ∫
0 0
π / 2 arcsin( y ) / y
j. ∫ ∫ y cos(xy )dxdy
0 0
π/2 3 sin( x )
k. ∫ ∫ x (1 + y)dydx .
0 sin( x )
2 3⎛⎜ 4 − x 2 ⎞⎟ / 2 ) ⎛ 5 ⎞
l. ∫∫
0
⎝ ⎠
−3⎛⎜ 4 − x 2 ⎞⎟ / 2
⎝ ⎠
⎜
⎝ 2+x
+ y 3 ⎟dydx .
⎠
1 x y
m. ∫ ∫ ∫ (y + xz )dzdydx
0 0 0
1 z y
∫ ∫ ∫ xy
2
n. z 3 dxdydz
0 0 0
1 y x/ 3
x
o. ∫∫ ∫
0 0 0 x + z2
dzdxdy
2
2 z 2
∫ ∫ ∫ yz dxdydz
2
p.
1 1 1/ y
∞ y
∫ ∫ xe
− y3
q. dxdy
0 0
002. El potencial gravitacional a una distancia R de la región comprendida entre dos esferas concéntricas
de radios ρ1 y ρ 2 , con ρ 2 > ρ1 , viene dado por la expresión:
ρ2 π 2 π
ρ 2 sin φdθdφdρ
V(0,0, R ) = ∫∫ ∫ ρ 2 sin 2 φ + (ρ cos φ − R )
2
ρ1 0 0
Evaluar la integral triple y encontrar una fórmula simplificada para el cálculo del potencial gravitacional.
Series de Taylor
001. Encontrar el polinomio de Taylor de grado n (Usar el comando taylor), alrededor del punto x0, para
las siguientes funciones:
a. f ( x ) = e , n = 6, x 0 = 0
x
b. f ( x ) = sin( x ), n = 4 y 6, x 0 = 0
c. f ( x ) = sin( x ), n = 6, x 0 = 2
d. f ( x ) = tan( x ), n = 6, x 0 = 0
e. f ( x ) = ln( x ), n = 4, x 0 = 1
f. f ( x ) = erf ( x ), n = 8, x 0 = 0
002. Escribir un archivo m de función llamado taylor1 para encontrar la serie de Taylor de primer orden,
alrededor del punto x0, de una función f : D ⊂ ℜ a ℜ . Recuérdese que la serie de Taylor de primer
3
f (x ) = f (x 0 ) + h∇f (x 0 ) + ...
Donde,
⎡ x 1 − x 1, 0 ⎤
⎢ ⎥
⎢ x 2 − x 2,0 ⎥
⎢x − x ⎥ ⎡ ∂f ⎤
h T = (x − x 0 ) T = ⎢ 3 3, 0
⎥ , ∇f (x 0 ) = ⎢ (x 0 ), ∂f (x 0 ),..., ∂f (x 0 )⎥
⎢ . ⎥ ⎣ ∂x 1 ∂x 2 ∂x n ⎦
⎢ ⎥
⎢ . ⎥
⎢⎣ x n − x n ,0 ⎥⎦
003. Escribir un archivo de función llamado taylor2 para encontrar la serie de Taylor de segundo orden,
alrededor del punto x0, de una función f : D ⊂ ℜ a ℜ . Recuérdese que la serie de Taylor de
3
1 T
f (x ) = f (x 0 ) + h∇f (x 0 ) + h Hf (x 0 )h + ...
2!
Donde
⎡ x 1 − x 1, 0 ⎤
⎢ ⎥
⎢ x 2 − x 2,0 ⎥
⎢x − x ⎥ ⎡ ∂f ⎤
h T = (x − x 0 ) T = ⎢ 3 3, 0
⎥ , ∇f (x 0 ) = ⎢ (x 0 ), ∂f (x 0 ),..., ∂f (x 0 )⎥
⎢ . ⎥ ⎣ ∂x 1 ∂x 2 ∂x n ⎦
⎢ ⎥
⎢ . ⎥
⎢⎣ x n − x n ,0 ⎥⎦
⎡ ∂ 2f ∂ 2f ∂ 2f ⎤
⎢ (x 0 ) (x 0 ) ..... (x 0 )⎥
⎢ ∂x21
2
∂x 1∂x 2 ∂x 1∂x n ⎥
⎢ ∂ f ∂ f ∂ 2f
(x 0 )⎥⎥
2
⎢ ∂x ∂x (x 0 ) 2
(x 0 ) .....
∂x 2 ∂x n
∂x 2
Hf (x 0 ) = ⎢ 2 1 ⎥
⎢ . . . . ⎥
⎢ . . . . ⎥
⎢ 2 ⎥
⎢ ∂ f (x ) ∂ f ∂ f
2 2
⎢⎣ ∂x n ∂x 1 0 ∂x n ∂x 2 0
( x ) ..... 2
(x 0 ) ⎥
⎥⎦
∂x n
004. Encontrar las series de taylor de primero y segundo orden alrededor del punto x0 , dadas las
funciones:
1
c. f ( x , y) = , x 0 = (0,0)
x + y2 + 1
2
d. f ( x , y) = e
− x 2 − y2
cos(xy ), x 0 = (0,0 )
( x −1)
cos(y ), x 0 = (1,0)
2
e. f ( x , y) = e
005. La velocidad de una reacción química está dada por la siguiente expresión
⎛ E ⎞ 2
r (C A , C B , T ) = k 0 exp⎜ − ⎟C A C B
⎝ RT ⎠
a. Obtener una aproximación lineal de ésta función (Serie de Taylor de primer orden), alrededor de los
valores de referencia CA0, CB0 y T0.
b. Obtener una aproximación cuadrática de ésta función (Serie de Taylor de segundo orden), alrededor
de los valores de referencia CA0, CB0 y T0.
006. La ecuación de Antoine para calcular la presión de vapor de una sustancia pura , p0, como una
función de la temperatura (T), está dada por:
⎛ B ⎞
p 0 (T ) = A exp⎜ − ⎟
⎝ T +C⎠
donde A, B, C son constantes.
a. Obtener una aproximación lineal de ésta función (Serie de Taylor de primer orden), alrededor de los
valores de referencia T0.
b. Obtener una aproximación cuadrática de ésta función (Serie de Taylor de segundo orden), alrededor
de los valores de referencia T0.
007. La densidad de un gas ideal (ρ), en función de la presión (p) y la temperatura (T), está dada por la
expresión
Mp
ρ(T, p ) =
RT
Donde M es el peso molecular del gas y R es la constante universal de los gases.
a. Obtener una aproximación lineal de ésta función (Serie de Taylor de primer orden), alrededor de los
valores de referencia p0 y T0.
b. Obtener una aproximación cuadrática de ésta función (Serie de Taylor de segundo orden), alrededor
de los valores de referencia p0 y T0.
Ecuaciones diferenciales
001. Resolver simbólicamente las siguientes ecuaciones diferenciales (Usar el comando dsolve).
dy
a. Ecuación logística: = Ay − By 2
dx
dy ⎛1 ⎞ 4
b. cos (x ) + 3y = 1, y⎜ π ⎟ =
2
dx ⎝4 ⎠ 3
dy
c. + ky = e 2 kx
dx
dy ⎛1 ⎞
d. + 4 y cot (2 x ) = 6 cos(2x ), y⎜ π ⎟ = 2
dx ⎝4 ⎠
dy
e. 2 y + y 2 sin (x ) = sin (x ), y(0) = 2
dx
e − x sinh (x )
3
dy
f. +x y= 2
dx 3y 2
dy
g. + y = y 2 , y(0) = -1
dx
h.
d2y
dx 2
+ 2k
dy
dx
( ) dy
+ k 2 + ω 2 y = 0, y(0) = 1, (0) = −k
dx
2
2 d y dy dy
i. x 2
− 4x + 6 y = 0, y(1) = 1, (1) = 0
dx dx dx
d2y dy 1
j. 2
+9 = cos(x ) + cos(3x )
dx dx 3
3
d y dy
k. + =0
dx 3 dx
d3y d2y dy dy d2y
l. − 9 + 27 − 27 y = 0, y(0) = 3, (0 ) = 10, =5
dx 3 dx 2 dx dx dx 2
d4y d2y dy d2y d3y
m. − 26 2 + 251y = 0, y(0) = 1, (0) = −1, 2 (0) = 5, 3 (0) = −10
dx 4 dx dx dx dx
⎧ dy1
⎪⎪ dt = 3y 2
n. ⎨
⎪ dy 2 = 12 y
⎪⎩ dt 1
⎧ dy1
⎪ dt = 2 y1 + 8y 2 − 4 y 3
⎪
⎪ dy
o. ⎨ 2 = −4 y1 − 10 y 2 + 2 y 3
⎪ dt
⎪ dy 3
⎪ dt = −4 y1 − 4 y 2 − 4 y 3
⎩
⎧ dy1
⎪⎪ dt = −2 y 2 + 4 t
⎨
dy
p. ⎪ 2 = 2 y1 − 2 t
⎪⎩ dt
1
y1 (0) = 4, y 2 (0) =
2
⎧ dy1
⎪⎪ dt = − y1 + 5y 2
⎨
q. ⎪ dy 2
= − y 1 + 3y 2
⎪⎩ dt
y1 (0) = 7, y 2 (0) = 2
Problemas de aplicación
Las distribuciones de velocidad y presión para el flujo reptante alrededor de una esfera son:
⎡ 3 ⎛ R ⎞ 1 ⎛ R ⎞3 ⎤
v r = v ∞ ⎢1 − ⎜ ⎟ + ⎜ ⎟ ⎥ cos θ
⎢⎣ 2 ⎝ r ⎠ 2 ⎝ r ⎠ ⎥⎦
⎡ 3⎛R⎞ 1⎛R⎞
3
⎤
v θ = v ∞ ⎢− 1 + ⎜ ⎟+ ⎜ ⎟ ⎥ sin θ
⎢⎣ 4⎝ r ⎠ 4⎝ r ⎠ ⎥⎦
vφ = 0
2
3 µv ∞ ⎛R⎞
p = p 0 − ρg(r cos θ) − ⎜ ⎟ cos θ
2 R ⎝r⎠
En la última ecuación, la cantidad p0 es la presión en el plano z=0, lejos de la esfera. El término -ρgz es
la presión hidrostática resultante del peso del fluido, y el término que contiene v∞ es la contribución del
movimiento del fluido.
⎡ ∂v ⎤
τ rr = −µ ⎢2 r ⎥
⎣ ∂r ⎦
⎡ ∂ ⎛ v ⎞ 1 ∂v r ⎤
τ rθ = −µ ⎢r ⎜ θ ⎟ + ⎥
⎣ ∂r ⎝ r ⎠ r ∂θ ⎦
b. Determinar la fuerza ejercida por el flujo de fluido sobre la esfera. Debido a la simetría alrededor del
eje z, la fuerza resultante estará en la dirección del eje z. Empezar con la integración de la fuerza normal:
En cada punto sobre la superficie de la esfera, el fluido ejerce una fuerza por unidad de área
− (p + τ rr ) r = R sobre el sólido, actuando en dirección normal a la superficie. La componente z de la
fuerza es − (p + τ rr ) r = R (cos θ ) . Multiplicamos la presión por elemento diferencial de superficie
R 2 sin θdθdφ para obtener la fuerza sobre el elemento de superficie. Calcular la fuerza normal
resultante en la dirección z:
∫ (− (p + τ ) )
2π π
Fn = ∫ rr r =R
cos θ R 2 sin θdθdφ
0 0
c. Integración de la fuerza tangencial: En cada punto de la superficie del sólido hay también un esfuerzo
cortante actuando tangencialmente. La fuerza tangencial resultante en la dirección z es:
∫ (τ )
2π π
Ft = ∫ rθ r = R sin θ R 2 sin θdθdφ
0 0
MÉTODOS NUMÉRICOS
∫e
sin( x )
a. dx
0
1
b. ∫
0
x 3 + 1dx
+∞
∫e
−x2
c. dx
−∞
Ecuaciones no lineales
(
b. x + y
2 2 2
) = 2xy y y = x 3
⎛ a ⎞
⎜ P + 2 ⎟(v − b ) = RT
⎝ v ⎠