Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Taller12 EDO 3

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

METODOS NUMERICOS 3006907

SEMESTRE 02, 2020, TALLER 12


Sistemas y problemas de orden superior

Rutinas disponibles:
En la colección de rutinas acompañantes del curso está el método de Runge-Kutta de orden 4 para
sistemas: rks4.m
En MATLAB: Para sistemas de ecuaciones diferenciales ordinarias están todas las rutinas que empiezan
por ode, por ejemplo, ode45.m, las cuales vienen preparadas para resolver sistemas y también pueden
resolver problemas escalares.

1. Considere un péndulo simple de longitud L = 2 p y suponga que la gravedad está dada por g =
32.17p/s2 . Con h = 0.1s, compare el ángulo θ que se obtiene para los siguientes problemas de valor
inicial en t = 0, 1 y 2 s.
d2 θ g
a) + senθ = 0, θ(0) = π/6, θ0 (0) = 0.
dt2 L
d2 θ g
b) + θ = 0, θ(0) = π/6, θ0 (0) = 0.
dt2 L
Al segundo problema se le llama linealización.
2. Un resorte resonante con función de forzamiento periódica se modela por medio de la ecuación
diferencial ordinaria de segundo orden

x00 (t) + 25x(t) = 8sen(5t), x(0) = 0, x0 (0) = 0.

Use el método de Runge-Kutta de cuarto orden para resolver esta ecuación en el intervalo [0, 2]
usando M = 40 pasos o sea h = 0.05.
3. Un circuito eléctrico de tipo RLC se modela por medio de la ecuación diferencial

Q00 (t) + 20Q0 (t) + 125Q(t) = 9sen(5t)

con Q(0) = 0 y Q0 (0) = 0. Use el método de Runge-Kutta de cuarto orden para resolver esta ecua-
ción en el intervalo [0, 2] usando M = 40 pasos o sea h = 0.05.

Nota: I(t) = Q0 (t) es la corriente en el tiempo t.


4. Modelo predador-presa: Sean x(t) y y(t) las poblaciones de conejos y de zorros respectivamente
en un tiempo t. El modelo predador-presa consiste en el sistema de ecuaciones diferenciales ordinarias

x0 (t) = Ax(t) − Bx(t)y(t) (1)


0
y (t) = Cx(t)y(t) − Dy(t). (2)

donde los parámetros tienen el siguiente significado:


A: Tasa de crecimiento de los conejos.
B: Tasa de muerte de los conejos, que depende tanto de la población de conejos como de la población
de zorros.
C: Tasa de crecimiento de los zorros, la cual depende de ambas especies.
D: Tasa de muerte de los zorros.

Considere los siguientes conjuntos de parámetros:


a) A = 2, B = 0.02, C = 0.0002, D = 0.8.
b) A = 1, B = 1, C = 1, D = 1.

1
Resuelva el sistema por el método de Runge-Kutta de cuarto orden en el intervalo [0, 10] con M = 50
con las siguientes condiciones iniciales:
a) x(0) = 3000 conejos, y(0) = 120 zorros.
b) x(0) = 5000 conejos, y(0) = 100 zorros.

Presente sus soluciones por medio de diagramas de fase.

También podría gustarte