Tema2 Ejercicios
Tema2 Ejercicios
Tema2 Ejercicios
—Ejercicios—
Ejercicio 1. Resuelve los siguientes PVI por medio de la transformada de Laplace. Representa gráficamente
la solución.
( (
y 0 + y = f (t) 1, 0 ≤ t < 1
a) donde f (t) =
y(0) = 0 0, caso contrario
( (
y 0 + y = f (t) 1, 1 ≤ t < 2
b) donde f (t) =
y(0) = 0 0, caso contrario
( (
y 00 + 4y = f (t) sin(t), 3π2 ≤t
c) 0
donde f (t) =
y(0) = 0, y (0) = 1 0, caso contrario
(
00 0
y + 4y + 3y = 1 − u(t − 1) + u(t − 2) − u(t − 4)
d)
y(0) = 0, y 0 (0) = 1
28 de diciembre de 2020
Ejercicio 5. Considera la ecuación
y 0 = t2 − y 2
en el intervalo [0, 5]. La utilizaremos para testear algunos métodos.
a) En primer lugar, resuelve esta ecuación utilizando métodos explícitos que ya tengas implementados
como Runge-Kutta IV y Crank-Nicolson, utilizando para ello las condiciones iniciales y(0) = −0.675 y
y(0) = −0.676. ode45 también te puede servir para este experimento. ¿Qué observas? ¿Qué deduces
de ello?
En adelante, se fija para la ecuación anterior la condición inicial y(0) = −0.675, queriéndose determinar el
valor de la solución en t = 5. Para ello, se presentan las siguientes tablas de Butcher.
0
1 1 1 1 1 1
3 3 4 4 0 2 2
3 2 1
I. 2
3 − 13 1 II. 4 4 4 III. 1 − 12 1
2
1 1 −1 1 1 1 1 1
2 2 2 2
1 3 3 1
8 8 8 8
b) Para cada tabla, escribe las ecuaciones correspondientes como métodos de Runge-Kutta. Clasifica
estos métodos como explícitos, diagonalmente implícitos o completamente implícitos.
c) Ejecuta 3 iteraciones completas del primer método, 2 del segundo y 1 del tercero, con paso de
integración h = 0.1. Calcula, cuando sea posible el error local y el global. Puedes apoyarte en
soluciones numéricas de métodos ya implementados con un paso relativamente más pequeño h =
0.001.
d) Implementa los métodos y determina el orden numérico de estos (4, 2 y 2, respectivamente), ana-
lizando el error global en t = 5. Puedes hacerlo frente a una solución considerada «exacta» o con
errores aproximados. Comprueba cuánto se reduce el error cuando se divide el paso por 2 (es decir,
se multiplica por 2 el número de iteraciones, N = 2, 4, 8, 16, . . .).
Ejercicio 6. Utilizando desarrollos de Taylor, calcula el error local del método de Crank-Nicolson. ¿Cuál es,
por tanto, el orden del método?
Ejercicio 7 ([1], 25.16). El oscilador armónico suele representarse como una masa puntual en el extremo
inferior de un muelle fijado por su extremo superior. La ecuación
d2 x dx
m 2
+c + kx = 0
dt dt
obtenida gracias a la ley de Hook, modela el desplazamiento vertical x (m) respecto a la posición de
equilibrio (x = 0) en función del tiempo t (s), siendo m (kg) la masa, c el coeficiente de amortiguamiento
(Ns/m), k la elasticidad (N/m) del muelle.
Utiliza un método numérico para simular posición y velocidad de este sistema cuando x(0) = 1 m, x0 (0) = 0
m/s, m = 20 kg, k = 20 N/m, en los siguientes casos:
a) Infra-amortiguado: c = 5 Ns/m (oscilamiento).
b) Amortiguado crítico: c = 40 Ns/m
c) Sobre-amortiguado: c = 200 Ns/m (sin oscilamiento).
Nota: el tipo amortiguamiento está relacionado con la forma de la solución, siendo el caso crítico c2 −4mk =
0.
Ejercicio 8 ([1], 25.26). Tres «puentistas» están atados entre sí. Cada cuerda-puentista se modela como un
oscilador con respecto a su punto de anclaje. Tomando en cuenta los equilibrios entre las fuerzas ejercidas
28 de diciembre de 2020
(ley de Hook) por los sistemas adyacentes, se define el siguiente sistema de ecuaciones diferenciales
d2 x1
m1 = m1 g − k1 (x1 − 0) + k2 (x2 − x1 )
dt2
d2 x2
m2 2 = m2 g − k2 (x2 − x1 ) + k3 (x3 − x2 )
dt
d2 x3
m3 2 = m1 g − k3 (x3 − x2 )
dt
donde mi es la masa del i-ésimo saltador; ki la constante de la cuerda, xi el desplazamiento (respecto a la
posición de equilibrio) y g = 9.81 m/s2 es la aceleración de la gravedad.
Simula el sistema suponiendo que las posiciones y velocidades iniciales son nulas para los tres saltadores
con lo siguientes datos: m1 = 60 kg, m2 = 70 kg, m3 = 80 kg, k1 = 50 N/m, k2 = 100 N/m, k3 = 50
N/m.
Referencias
[1] Steven C. Chapra y Raymond P. Canale. Métodos numéricos para ingenieros (7ª edición). McGraw-Hill, Education, 2015.
28 de diciembre de 2020