0% encontró este documento útil (0 votos)
98 vistas15 páginas

Actividad 1 Resolución de PVIs Con Métodos Numéricos - After 06

Descargar como pdf o txt
Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1/ 15

Actividad 1

Resolución de PVIs con métodos numéricos

Dra. Alicia Cordero


Neus Garrido
Dr. Juan R. Torregrosa

Máster en Ingenierı́a Matemática y Computación


Escuela Superior en Ingenierı́a y Tecnologı́a

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 1 / 15


Contenido

1 Caracterı́sticas del trabajo

2 Problemas
Problema 1
Problema 2

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 2 / 15


1

Caracterı́sticas del trabajo

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 3 / 15


Caracterı́sticas del trabajo

Actividad individual
Entrega en Word
Uso del editor de ecuaciones
No utilizar capturas de pantalla
Revisar Gráficas.pdf
Resultados con 6 decimales
10 páginas, Calibri 11, Interlineado 1.5
Criterios de evaluación
Presentación: 1.5 puntos
Problema 1: 4.5 puntos
Problema 2: 4 puntos

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 4 / 15


2

Problemas

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 5 / 15


Problemas

Se trata de dos ejercicios en los que deberás programar en Matlab las funciones
necesarias para la resolución de los problemas de valor inicial, por medio de los métodos
numéricos descritos en teorı́a.

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 6 / 15


Contenidos

1 Caracterı́sticas del trabajo

2 Problemas
Problema 1
Problema 2

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 7 / 15


Problema 1

PVI a través de ED de segundo orden (1 de 3)

Un sistema resonante de muelles sobre el que se


ejerce una fuerza externa periódica se modela 𝑘 = 25
mediante la ecuación 4sin(5𝑡)
x00 (t) = 4 sin(5t) − 25x(t), x(0) = x0 (0) = 0. 𝑚=1

a) (1 punto) Transforma el PVI en un sistema de ecuaciones diferenciales de primer orden.


Escribe una función PVI1.m que implemente el sistema de ecuaciones diferenciales de primer
orden y copia el código en este apartado.
Solución
  0 function dX=PVI1(t,X)
x1 = · · · x1 = · · · ...

x2 = · · · x02 = · · · end

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 8 / 15


Problema 1

PVI a través de ED de segundo orden (2 de 3)


b) (1.25 puntos) Usa el método de Heun de orden 2 para resolver el PVI en el intervalo [0, 2]
con 40 subintervalos. Representa la solución x(t) para t ∈ [0, 2]. Indica en una tabla los
valores de x(t) para t ∈ {0, 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2}.
Solución
Heun 40 subintervalos
1

>> [t,X]=Heun...
0.5
x(t)

0
t 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00
-0.5
x(t)
-1
0 0.5 1 1.5 2
t

c) (1.25 puntos) Usa el método de Runge-Kutta de orden 4 para resolver el PVI en el intervalo
[0, 2] con 40 subintervalos. Representa la solución x(t) para t ∈ [0, 2]. Indica en una tabla
los valores de x(t) para t ∈ {0, 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2}.
Solución
RK 40 subintervalos
1

0.5
>> [t,X]=RK4...
x(t)

0
t 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00
-0.5
x(t)
-1
0 0.5 1 1.5 2
t

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 9 / 15


Problemas

PVI a través de ED de segundo orden (3 de 3)


d) (1 punto) Proporciona una estimación numérica del orden de los métodos utilizados en los
apartados anteriores. Copia el código que has utilizado para realizar dicha estimación
numérica.

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 10 / 15


Contenidos

1 Caracterı́sticas del trabajo

2 Problemas
Problema 1
Problema 2

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 11 / 15


Problema 2

Sistema de tres PVIs a partir de EDs de primer orden (1 de 3)


El modelo de Lorenz-Haken para un láser de dos niveles viene dado por el sistema de ecuaciones
diferenciales 
 F 0 (t) = 3 (P (t) − F (t))
P 0 (t) = −P (t) + d(t)F (t)
 d0 (t) = 17 − d(t) − F (t)P (t)

donde F (t) es el campo, P (t) la polarización y d(t) la diferencia de polarizaciones.


a) Escribe una función PVI2.m que implemente el sistema de ecuaciones diferenciales de primer
orden y copia el código en este apartado.
Solución
function dFPd=PVI2(t,FPd)
...
end

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 12 / 15


Problemas

Sistema de tres PVIs a partir de EDs de primer orden (2 de 3)


b) Resuelve el sistema en el intervalo temporal t ∈ [0, 2] por el método de Runge-Kutta de
orden 4, tomando como paso h = 0.2 y como condiciones iniciales
F (0) = P (0) = 4.5, d(0) = 1. Representa las variables F , P y d en el intervalo t ∈ [0, 2] e
indica en una tabla los valores de F , P y d para t ∈ {0, 0.6, 1, 1.8, 2}.
Solución
2
F(t) >> [t,FPd]=RK4(...);
P(t)
1
d(t)
t 0.0 0.6 1.0 1.8 2.0
0
F (t)
-1 P (t)
-2
d(t)
0 0.5 1 1.5 2
t

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 13 / 15


Problemas

Sistema de tres PVIs a partir de EDs de primer orden (3 de 3)


c) Resuelve el sistema en el intervalo temporal t ∈ [0, 200] por el método de Runge-Kutta de
orden 4, tomando como paso h = 0.2 y como condiciones iniciales
F (0) = P (0) = 4.5, d(0) = 1. Representa las variables F , P y d en el intervalo
t ∈ [190, 200] e indica en una tabla los valores de F , P y d para t ∈ {190, 195, 200}.
Solución
2
F(t) >> [t,FPd]=RK4(...);
P(t)
1
d(t)
t 190 195 200
0
F (t)
-1 P (t)
-2
d(t)
190 195 200
t

A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 14 / 15


A Cordero + N Garrido + JR Torregrosa PER1071 (Otoño 2019) 15 / 15

También podría gustarte