3Eva_2024PAOI_T3 Salto Bungee y método de Jacobi

3ra Evaluación 2024-2025 PAO I. 17/Septiembre/2024

Tema 3. (35 puntos) Para el salto del Bungee y la tabla del ejercicio del tema 2, realice el sistema de ecuaciones para el polinomio de interpolación usando los tiempos en el vector ts = [0, 0.75,1.375, 2.55]

datos usando h=1/8
ti vi
0 0.0000
0.25 2.4479
0.5 4.8849
0.75 7.3001
1 9.6832
1.375 13.1763
1.75 16.5451
2.125 19.7641
2.4 22.0193
2.55 23.2075
  1. Plantee el sistema de ecuaciones usando los cuatro puntos datos en el vector ts
  2. Presente el sistema de ecuaciones en su forma matricial y muestre la matriz aumentada y pivoteada
  3. Desarrolle el ejercicio usando el método de Jacobi, para tres iteraciones. Justifique el vector inicial X0
  4. Comente sobre la convergencia del ejercicio y adjunte los archivos para algoritmo.py y resultados.txt

Rúbrica: Planteamiento (5 puntos), Matriz aumentada y pivoteada (5 puntos), iteraciones(15 puntos), error por iteración (5 puntos), literal d (5 puntos)

3Eva_2024PAOI_T2 Salto Bungee interpolar velocidad primer tramo

3ra Evaluación 2024-2025 PAO I. 17/Septiembre/2024

Tema 2. (30 puntos) Bungee Jumping 02
Para el salto del Bungee del ejercicio del tema anterior se toman lecturas con un sensor de velocidad sujetado a la persona.

De la tabla de datos obtenida, se observa que los tamaños de paso en tiempo no siempre son equidistantes.
Se requiere encontrar un polinomio interpolación de al menos grado 3.

datos usando h=1/8
ti vi
0 0.0000
0.25 2.4479
0.5 4.8849
0.75 7.3001
1 9.6832
1.375 13.1763
1.75 16.5451
2.125 19.7641
2.4 22.0193
2.55 23.2075

a. Describa el planteamiento del ejercicio, selección de puntos, expresiones usando el método de Lagrange.

b. Resuelva el sistema usando los algoritmos correspondientes.

c. Presente el polinomio obtenido, simplificado y grafique verificando que P(x) pase por los puntos de muestra.

d. Use el resultado P(x) para estimar el error con los valores de vi que no fueron usados para la interpolación.

e. Adjunte los archivos para algoritmo.py, resultados.txt y grafica.png

Rúbrica: literal a (10 puntos), literal b (5 puntos), literal c (5 puntos), literal d (5 puntos), literal e (5 puntos)

Referencia: Chapra. capítulo 28. Ejercicio 28.41 p852.

3Eva_2024PAOI_T1 EDO Salto Bungee extiende y estira toda la cuerda

3ra Evaluación 2024-2025 PAO I. 17/Septiembre/2024

Tema 1. (35 puntos) Salto Bungee 01
La ecuación diferencial para la velocidad de alguien que practica el salto del bungee es diferente según si el saltador ha caído una distancia en la que la cuerda está extendida por completo y comienza a estirarse o encogerse.

Si la distancia recorrida es menor que la longitud de la cuerda, el saltador sólo está sujeto a las fuerzas gravitacional y de arrastre de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 y \leq L

Salto Bungee 02Suponga que las condiciones iniciales son:

y(0) =0

\frac{dy(0)}{dt} = 0

Una vez que la cuerda comienza a estirarse, también deben incluirse las fuerzas del resorte y del amortiguamiento de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 -\frac{k}{m}(y-L) - \frac{\gamma}{m}( \frac{dy}{dt}\Big) y \gt L
dy/dt m/s velocidad (v)
t s tiempo
g 9.8 m/s2 gravedad
cd 0.25 kg/m coeficiente de arrastre
m 68.1 Kg masa
L 30 m Longitud de la cuerda
k 40 N/m constante de resorte de la cuerda
γ 8 N s/m coeficiente de amortiguamiento de la cuerda
signo(v) función que devuelve –1, 0 y 1, para v negativa, cero y positiva, respectivamente

En conocimiento que la primera ecuación es válida solo hasta tc=2.55, L=30mts, v= 23.20752.
Encuentre el tiempo td cuando se alcanza la longitud MÁXIMA de la cuerda extendida y estirada por completo, es decir y>L, con velocidad = 0. (solo 2da ecuación)

a. Realice el planteamiento del ejercicio usando Runge-Kutta de 2do Orden

b. Desarrolle tres iteraciones para y(t) con tamaño de paso h=0.5

c. Usando el algoritmo, aproxime la solución para y en el intervalo entre [0,tc], adjunte sus resultados.txt

d. Indique el valor de td, muestre cómo mejorar la precisión y realice sus observaciones sobre los resultados.

Observación: dy/dt = v, función signo(v) en Numpy:np.sign(v)

Rúbrica: literal a (5 puntos), literal b (15 puntos), literal c resultados.txt y grafica.png (10 puntos), literal d (5 puntos),

Referencia: [1] Chapra. capítulo 28. Ejercicio 28.41 p852.

[2] Extreme Bungy Jumping with Cliff Jump Shenanigans! Play On in New Zealand! 4K! – devinsupertramp. 23 mar 2015.

 

2Eva_2024PAOI_T3 EDP Parabólica

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 3. (30 puntos)

Para la siguiente Ecuación Diferencial Parcial con b = 2, resuelva usando las condiciones mostradas

\frac{\partial ^2 u}{\partial x^2} = b \frac{\partial u}{\partial t}
0 < x < 1  0 < t < 0.5
Condiciones iniciales:  u(x,0)=0
Condiciones de frontera: u(0,t)=1
u(1,t)= 2

Utilice diferencias finitas centradas y hacia adelante para las variables independientes x,t

a. Plantee las ecuaciones para usar un método numérico en un nodo i,j

b. Realice la gráfica de malla,

c. Desarrolle y obtenga el modelo discreto para u(xi,tj)

d. Realice al menos tres iteraciones en el eje tiempo.

e. Estime el error de u(xi,tj) y adjunte los archivos del algoritmo y resultados.

Rúbrica: Aproximación de las derivadas parciales (5 puntos), construcción de la malla (5), desarrollo de iteraciones (15), literal e (5 puntos)

Referencia: EDP Parabólicas. Chapra & Canale. 5ta Ed. Ejercicio 30.15. P.904

2Eva_2024PAOI_T2 Salto Bungee longitud total de cuerda

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 2. (40 puntos) Bungee Jumping 02
Para el salto del Bungee del ejercicio del tema anterior se toman lecturas con un sensor de velocidad sujetado a la persona.

2.1 De la tabla de datos obtenida, se observa que los tamaños de paso en tiempo no siempre son equidistantes.
Se requiere encontrar la distancia recorrida en el intervalo de [0,2.55] usando fórmulas de integración compuestas.

ti vi
0 0.0000
0.25 2.4479
0.5 4.8849
0.75 7.3001
1 9.6832
1.375 13.1763
1.75 16.5451
2.125 19.7641
2.4 22.0193
2.55 23.2075

2.2 Usando los datos de la tabla para el intervalo [2.55, 5.175] donde la velocidad de la caída de la persona al primer salto ha llegado a casi cero, o antes del primer rebote, se ha obtenido un polinomio de interpolación:

v = -3.979t2 + 21.557t – 5.3997

Obtenga la distancia recorrida en el segundo intervalo usando el método de Cuadratura de Gauss.

a. Realice el planteamiento de las ecuaciones para cada sección del ejercicio.

b. Describa el criterio usado para determinar el número de tramos usado en cada caso.

c. Desarrolle las expresiones completas del ejercicio para cada sección.

d. Encuentre la distancia total (profundidad) alcanzada por la persona al dar el salto.

Rúbrica: literal a 2.1 (5 puntos), a 2.2 (5 puntos) literal b (5 puntos), literal c 2.1 (10 puntos), c 2.2 (10 puntos), literal d (5 puntos)

Referencia:[1] Chapra. capítulo 28. Ejercicio 28.41 p852.

[2] Extreme Bungy Jumping with Cliff Jump Shenanigans! Play On in New Zealand! 4K! – devinsupertramp. 23 mar 2015.

2Eva_2024PAOI_T1 EDO Salto Bungee tiempo extiende cuerda

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 1. (30 puntos) Salto Bungee 01
La ecuación diferencial para la velocidad de alguien que practica el salto del bungee es diferente según si el saltador ha caído una distancia en la que la cuerda está extendida por completo y comienza a estirarse o encogerse.

Si la distancia recorrida es menor que la longitud de la cuerda, el saltador sólo está sujeto a las fuerzas gravitacional y de arrastre de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 y \leq L

Salto Bungee 02Suponga que las condiciones iniciales son:

y(0) =0

\frac{dy(0)}{dt} = 0

Una vez que la cuerda comienza a estirarse, también deben incluirse las fuerzas del resorte y del amortiguamiento de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 -\frac{k}{m}(y-L) - \frac{\gamma}{m}( \frac{dy}{dt}\Big) y \gt L
dy/dt m/s velocidad (v)
t s tiempo
g 9.8 m/s2 gravedad
cd 0.25 kg/m coeficiente de arrastre
m 68.1 Kg masa
L 30 m Longitud de la cuerda
k 40 N/m constante de resorte de la cuerda
γ 8 N s/m coeficiente de amortiguamiento de la cuerda
signo(v) función que devuelve –1, 0 y 1, para v negativa, cero y positiva, respectivamente

Encuentre el tiempo tc y la velocidad de la persona cuando se alcanza la longitud de la cuerda extendida y sin estirar (30 m), es decir y<L, aún se entra cayendo signo(v)=1. (solo primera ecuación)

a. Realice el planteamiento del ejercicio usando Runge-Kutta de 2do Orden

b. Desarrolle tres iteraciones para y(t) con tamaño de paso h=0.5

c. Usando el algoritmo, aproxime la solución para y en el intervalo entre [0,tc], adjunte sus resultados.txt

d. Indique el valor de tc, muestre cómo mejorar la precisión y realice sus observaciones sobre los resultados.

Observación: dy/dt = v, función signo(v) en Numpy:np.sign(v)

Rúbrica: literal a (5 puntos), literal b (15 puntos), literal c resultados.txt y grafica.png (5 puntos), literal d (5 puntos),

Referencia: [1] Chapra. capítulo 28. Ejercicio 28.41 p852.

[2] Extreme Bungy Jumping with Cliff Jump Shenanigans! Play On in New Zealand! 4K! – devinsupertramp. 23 mar 2015.

 

1Eva_2024PAOI_T3 Tasas de natalidad de reemplazo en Ecuador

1ra Evaluación 2024-2025 PAO I. 2/Julio/2024

Tema 3 (30 puntos) El  acelerado crecimiento de la población de mayor edad tendrá un fuerte impacto en los gastos futuros del Instituto de Seguridad Social (IESS), una entidad a la que actualmente (2024) no le alcanzan sus ingresos para pagar las pensiones a sus jubilados y otros rubros. [1] trabajadores vs Pensionistas

Un motivo a considerar, podría ser la caída de la natalidad es un fenómeno generalizado en los países desarrollados que tienen tasa de hijos por mujer en edad fértil por debajo de la tasa de reemplazo poblacional de 2.1

Ejemplos de estos casos son Corea del Sur, Japón y China [2,3], países donde la cantidad de personas que trabajan y aportan al seguro social tiende a ser cada vez menor respecto a los pensionistas (jubilados).

Para el caso de Ecuador a fin de realizar un análisis preliminar, se requiere disponer de un modelo matemático que permita estimar cuando ser alcanzaría esta tasa de natalidad.

Promedio Hijos por mujer en Ecuador. INEC 29 Agosto 2019 [4]
tasa 7.32 6.39 5.58 4.89 4.31 3.85 3.50 3.26 3.03 2.79 2.54 2.35
año 1965 1970 1975 1980 1985 1990 1995 2000 2005 2010 2015 2020
k 0 1 2 3 4 5 6 7 8 9 10 11

Realice un modelo de interpolación polinómica usando los datos de los años 1965, 1980, 1995 y 2010.

a. Describa el planteamiento del ejercicio, justificando el grado del polinomio seleccionado

b. Realice el del sistema de ecuaciones en su forma matricial y muestre la matriz aumentada

c. Resuelva el sistema usando los algoritmos correspondientes.

d. Presente el polinomio obtenido y grafique verificando que P(x) pase por los puntos de muestra.

e. Use el resultado P(x) para estimar la tasa de hijos por mujer para el año 2020 y calcule el error.

f. Estime el año en que se alcanza la tasa mínima de reemplazo de 2.1

Adjunte los archivos de resultados, algoritmos y gráficas realizados para el ejercicio

Rúbrica: literal a (5 puntos), literal b (5 puntos), literal c (5 puntos), literal d (5 puntos), literal e (5 puntos) ), literal f (5 puntos).

Referencia:  [1] IESS: Envejecimiento de la población disparará gastos; Peña plantea subir aportes. Primicias.ec: 3 Marzo 2024. https://www.primicias.ec/noticias/economia/envejecimiento-poblacion-iess-aportes/
[2] «Emergencia nacional» en Corea del Sur: por qué las mujeres surcoreanas no están teniendo hijos. BBC News Mundo. 30 marzo 2024.

[3] Por qué China amplió a 3 el número de hijos que pueden tener las parejas. BBC News Mundo. 4 junio 2021.
https://www.youtube.com/watch?v=8kErwjPKwjY .
[4] 5 datos sobre población del Ecuador. INEC Ecuador. 29 Agosto 2019. Min 0:45 https://www.youtube.com/watch?v=wjTZNfykmZU [5] ¿Por qué los PAÍSES RICOS se enfrentan al COLAPSO DEMOGRÁFICO? VisualEconomik. 11 Octubre 2022.

import numpy as np
tasa = [7.32, 6.39, 5.58, 4.89, 4.31, 3.85, 3.50, 3.26, 3.03, 2.79, 2.54, 2.35]
anio = [1965, 1970, 1975, 1980, 1985, 1990, 1995, 2000, 2005, 2010, 2015, 2020]
k    = [   0,     1,   2,    3,    4,   5,     6,    7,    8,    9,   10,   11]
cual = [0,3,6,9]
tasamin = 2.1
fi = np.take(tasa,cual) # datos seleccionados

1Eva_2024PAOI_T2 Temperatura en nodos de placa cuadrada

1ra Evaluación 2024-2025 PAO I. 2/Julio/2024

Tema 2 (40 puntos) La distribución de temperatura en estado estable en una placa cuadrada caliente está modelada por la ecuación de Laplace [1], cuya solución en su forma iterativa cuando el factor
(Δy)2/(Δx) = 1 se interpreta como:

“La temperatura en los nodos de la malla de una placa se puede calcular con el promedio de las temperaturas de los 4 nodos vecinos de la izquierda, derecha, arriba y abajo” [2].

Considere placa cuadrada de 4.5 cm de lado tiene la temperatura en los nodos de los bordes como se indica en la figura. Placa Cuadrada Calentada Nodos01

a) Plantee el sistema de ecuaciones para encontrar los valores en los nodos a, b, c, d. Use la solución descrita para la ecuación de Laplace.

b) Presente la matriz aumentada y Muestre los pasos detallados para el pivoteo parcial por filas.

c) Desarrolle las expresiones para resolver mediante el método de Gauss-Seidel. Considere para el vector inicial Xo, valores intermedios entre las temperaturas de los bordes de la placa.

d) Realice al menos 3 iteraciones, indicando el error por iteración.

e) Analice la convergencia del método, número de condición y resultados obtenidos.

Adjunte los archivos del algoritmo y resultados de computadora utilizados.

Rúbrica: Literal a (5 puntos), literal b (5 puntos), literal c (5 puntos), literal d (15 puntos). literal e (5 puntos) Adjuntos (5 puntos)

Referencia: [1] Ejercicio 12.39 p339 Steven C. Chapra. Numerical Methods 7th Edition.
[2] Ecuaciones Elípticas. Método iterativo. http://blog.espol.edu.ec/analisisnumerico/edp-elipticas-metodo-iterativo/

1Eva_2024PAOI_T1 Vaciado de reservorio semiesférico

1ra Evaluación 2024-2025 PAO I. 2/Julio/2024

Tema 1. (30 puntos) Un reservorio semiesférico de radio R = 3 m, está lleno agua hasta la altura h como se muestra en la figura. En la base tiene un tubo de salida con abertura de área a = 0.01 m2. El coeficiente de fricción hidráulico K = 0.85   g = 9.8 m/s2.
tanque semiesfera
Se requiere que el reservorio se vacíe en menos de tf = 15 min.
La ecuación diferencial para el vaciado de tanques es:

A(h)\frac{dh}{dt} = -Ka\sqrt{2gh}

Considerando la relación de la altura del líquido h con
respecto al radio R de la semiesfera: r^2 + \Big( R-h\Big)^2 = R^2
Y el Área del círculo en función de h: A(h) = \pi \Big(2Rh -h^2 \Big)

se encuentra la solución general en la expresión:

\frac{4}{3}R \Big(h_f^{3/2} - h_0^{3/2} \Big) - \frac{2}{5}\Big(h_f^{5/2} - h_0^{5/2} \Big) = \frac{Ka}{\pi}\sqrt{2g}\Big(t_0-t_f \Big)

Dado que se vacía el reservorio, hf = 0 y que el experimento inicia en t0=0, h0 =h y tf=t, la expresión se simplifica:


-\frac{4}{3}R h^{3/2} + \frac{2}{5} h^{5/2} = -\frac{Ka}{\pi} t \sqrt{2g}

a. Plantear el ejercicio para encontrar h para un t dado, muestre el intervalo de búsqueda y una gráfica.

b. Desarrolle usando el método de Newton-Raphson para tres iteraciones y tolerancia milimétrica.

c. Verifique el orden de convergencia y observe sus resultados usando el algoritmo.

Nota: la gravedad g tiene unidades de segundos, el tiempo de vaciado está en minutos.

Rúbrica: Planteamiento (5 puntos), iteraciones y error (15 puntos), análisis de la convergencia (5 puntos). observación de resultados, algoritmo y gráficas adjuntos (5 puntos).

Referencia:
[1] Ejercicio 25.21 p765 y 5.17. p143 Steven C. Chapra. Numerical Methods 7th Edition.
[2] Vaciado de un tanque semiesférico. Tiempo total de vaciado. Demostración y aplicación. Sebastian Rodriguez

3Eva_2023PAOII_T3 Volumen por solido de revolución de un peón

3ra Evaluación 2023-2024 PAO II. 15/Febrero/2024

Tema 3 (40 puntos) Los sólidos de revolución se generan al girar una región plana alrededor de un eje.un peon 3D

El volumen generado al girar la región de una función f(x) en el intervalo [a,b], se puede calcular como el volumen del disco de radio f(x) y anchura dx.

V = \int_{a}^{b} \pi (f(x))^2 dx

Calcule el volumen de revolución generado por la región sombreada y limitada los puntos de la tabla del tema anterior y la gráfica 2D mostrada. volumen de un Peon 2D de revolucion

Realice el ejercicio usando para los integrales el método de integración por Cuadratura de Gauss para al menos lo tres primeros intervalos.

xi=[ 0, 3, 5.  , 9.985 , 14.97 , 17.97, 40.04, 43.29, 51.6449, 60]
yi=[15,15,13.25,14.1552, 9.6768,  9.67,  4.64,  4.64,  8.9768, 0.]

Para el desarrollo de cada intervalo:
a. Realice el planteamiento de las formulas de volumen de sólido de revolución.
b. Desarrolle las expresiones completas con valores numéricos que permitan revisar sus operaciones.
c. Indique el resultado obtenido para cada integral.
d. Determine el volumen de revolución generado por la región sombreada presentada en la gráfica usando el algoritmo en Python.
e. Adjunte sus resultados.txt y algoritmos.py

Rúbrica: literal a (12 puntos), literal b (12 puntos), literal c (6 puntos), literal d (5 puntos), literal e (5 puntos)

Referencia: [1] Volúmenes de sólidos de revolución. Moisés Villena Muñoz. Capítulo 4 p78. https://www.dspace.espol.edu.ec/bitstream/123456789/4800/4/7417.pdf

[2]Curso Torno Madera. Práctica de realización de peón de ajedrez. Taller Escuela Pinocho. 21 octubre 2021