2Eva_2022PAOII_T2 EDO – población de protestantes en una sociedad

2da Evaluación 2022-2023 PAO II. 24/Enero/2023

Tema 1. (35 puntos) protestantismoEn el libro titulado “Looking at History Through Mathematics”, Rashevsky propone un modelo que se puede relacionar con el “protestantismo” en el siglo XVI como una reacción y denuncia de abusos impuestos sobre la sociedad de la época.

En un modelo de Rashevsky modificado con la ecuación logística de Verhulst, la población x(t) de individuos en la sociedad para cada año t, con tasas de natalidad b=0.02 y mortalidad d=0.015, cambia según la ecuación:

\frac{\delta}{\delta t}x(t) = b x(t) - d (x(t))^2 x(0)=1

La cantidad de individuos “protestantes” y(t) en la población se incrementa según la ecuación diferencial compuesta de dos términos.

\frac{\delta}{\delta t}y(t) = b y(t) - d (y(t))^2 +r b (x(t)-y(t)) y(0)=0.01

El primer término supone que todas familias de padre y madre “protestantes” tienen hijos que también se identifican como tales.

El segundo término supone que una porción r = 0.1 de jóvenes descendientes de los “conformistas” al meditar sobre la situación actual, los hechos y los argumentos de protesta se convierten a “protestantes”.

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

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

c.       Usando el algoritmo, aproxime la solución entre t=0 a t=200 años, adjunte sus resultados en la evaluación.

d.       Realice una observación sobre el crecimiento de población y(t) a lo largo del tiempo.

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

Referencia: Burden 5.2 Ejercicio 17 p276, Rashevsky, MIT 1968. pp102-110, Protestantismo https://es.wikipedia.org/wiki/Protestantismo. 3Eva_IIT2014_T2 Crecimiento demográfico. http://blog.espol.edu.ec/analisisnumerico/3eva_iit2014_t2-crecimiento-demografico/

La Reforma protestante y Lutero. Academia Play. 27 agosto 2019

 

2Eva_2022PAOII_T1 Altura de cohete en 30 segundos

2da Evaluación 2022-2023 PAO II. 24/Enero/2023

Tema 1. (30 puntos) La velocidad hacia arriba de un cohete se calcula con la fórmula:

v = u \ln\Big(\frac{m_0}{m_0-qt}\Big) - gt

Donde:https://www.debate.com.mx/Las-increibles-imagenes-del-lanzamiento-del-cohete-mas-potente-del-mundo-l201802060004.html
v   = velocidad hacia arriba,
u   = 1800 m/s, velocidad a que se expele el combustible en relación con el cohete,
m0 = 160 000 kg, masa inicial del cohete en el tiempo t = 0,
q    = 2 500 kg/s,  tasa de consumo de combustible y
g    = 9.8 m/s2, aceleración de la gravedad

Para determinar la altura alcanzada por el cohete en un vuelo de 30 segundos desarrolle la parte analítica con los siguientes métodos y compare los resultados.

a. Utilice la regla de Simpson, en el planteamiento incluya la cantidad de tramos o segmentos a usar

b. Use el método de cuadratura de Gauss para la misma cantidad de segmentos que el literal anterior

c. Compare y comente los resultados, sobre los errores entre los métodos.

Rúbrica: Planteamiento de tramos (5 puntos), integral con Simpson (10 puntos), cuadratura de Gauss (10 puntos), literal c (5 puntos).

Referencia: Chapra ejercicio 24.46 p701. NASA y SpaceX realizan con éxito el despegue del primer vuelo de EE. UU. hacia la Estación Espacial Internacional en nueve años. EFE 30 mayo 2020 https://youtu.be/npcgpQUKAbg

 

 

2Eva_2022PAOI_T3 EDP parabólica barra enfriada en centro

2da Evaluación 2022-2023 PAO I. 30/Agosto/2022

Tema 3. (40 puntos) Use el método de diferencias progresivas para aproximar la solución de la siguiente ecuación diferencial parcial parabólica:

\frac{\partial U}{\partial t} - \frac{1}{9} \frac{\partial ^2 U}{\partial x^2} = 0 0 \leq x \leq 2, t>0

Con las condiciones iniciales de borde e iniciales:

U(0,t) = U(2,t) = 0, t>0 U(x,0) = \cos \Big( \frac{\pi}{2}(x-3)\Big) , 0 \leq x \leq 2

Aplique un método numérico para encontrar los valores de U(x,t) usando Δx = 1/3, Δt = 0.02 y muestre:

a. La grafica de malla
b. Ecuaciones de diferencias divididas  a usar
c. Encuentre las ecuaciones considerando las condiciones dadas en el problema.
d. Determine el valor de λ, agrupando las constantes durante el desarrollo, revise la convergencia del método.
e. Resuelva para tres pasos
f. Estime el error (solo plantear)
g. Usando el algoritmo, aproxime la solución para t=0.02 y t=0.1

Rúbrica: literal a (3 puntos), literal b (2 puntos), literal c (5 puntos), literal d (5 puntos), aplicación de condiciones iniciales (5 puntos), literal e (10 puntos), literal f (5 puntos). literal g, usando algoritmo (5 puntos)

Referencia: 2Eva_IT2017_T3 EDP parabólica http://blog.espol.edu.ec/analisisnumerico/2eva_it2017_t3-edp-parabolica/


2Eva_2022PAOI_T2 EDO de circuito RLC con interruptor intermedio

2da Evaluación 2022-2023 PAO I. 30/Agosto/2022

Tema 2. (30 puntos) El circuito de la figura 2a tiene el interruptor en posición cerrada por largo tiempo antes de t=0, con lo que la corriente en el inductor será de 2 Amperios, y(0)=2. Para t<0, el inductor opera como un conductor sin caída de voltaje, el capacitor está cargado a 10V y solo pasaría corriente por la resistencia de 5 Ohm.


En el tiempo t=0, el interruptor se abre de forma instantánea y el circuito cambia al modelo de la figura 2b.


La corriente del inductor y(t) para t≥0 está dada por la ecuación:

\frac{\delta}{\delta t}y(t) + 2 y(t) + 5 \int_{-\infty}^t y(\tau) \delta \tau = 10 \mu(t)

En t=0, luego de abrir el interruptor, los voltajes de la fuente y el capacitor son iguales. La corriente inicial sobre el resistor de 2 A genera un voltaje que se compensa con el voltaje del inductor pero en signo opuesto. Lo que implica que y’(0) = -4

V_{Inductor} = - V_{resistor} y'(0) = -4

Derive la expresión de corrientes y(t) para obtener una ecuación diferencial ordinaria.

a) Realice el planteamiento del problema usando el método de Runge-Kutta de 2do orden para 2da derivada

b) Desarrolle las expresiones para al menos tres iteraciones usando h=0.01

c) Estime el valor del error.

d) Muestre el resultado con el algoritmo para el intervalo t entre [0,5] segundos

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

Referencia: Lathi B.P. Green R. Linear Systems and Signals, 3rd Edition. ejemplo 4.13 p364

2Eva_2022PAOI_T1 Comparar integrales numéricos Simpson y Cuadratura de Gauss

2da Evaluación 2022-2023 PAO I. 30/Agosto/2022

Tema 1. (30 puntos) Determine el área bajo la curva dada por la expresión mostrada para el intervalo de x entre [0,3]:

A = \int_0^3 \frac{e^x \sin(x)}{1+x^2} \delta x

Desarrolle el ejercicio mostrando las expresiones completas para integración numérica usando:

a) Un método de Simpson aplicado al menos dos veces para el intervalo del integral. Determine el tamaño de paso propuesto y el número de puntos necesario para usar un solo método.

b) El método de Cuadratura de Gauss de dos puntos, usando dos tramos en el intervalo.

c) Estime el error de integración para los literales a y b. Compare los resultados obtenidos.

Rúbrica: Literal a. tamaño de paso (5 puntos) expresiones correctas y completas (10 puntos), literal b (10 puntos), literal c (5 puntos)

Referencia: Chapra 5Ed. ejercicio 22.14 p667