Laboratorio 1
Problemas de contorno unidimensionales, métodos de disparo
Métodos Numéricos Aplicados II
Eva García Villalba
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 1 / 14
Introducción
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 2 / 14
Índice
1 Objetivos
2 Conceptos previos: Disparo para PF no lineales
Condiciones Dirichlet
3 Laboratorio I
Criterios de calificación
Descripción del laboratorio
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 3 / 14
Objetivos
Objetivos
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 4 / 14
Objetivos
Objetivos
→ Aprender a describir el problema de disparo con esquema de Secante y esquema
de Newton.
→ Implementar código para los métodos de disparo no lineal.
→ Aplicar los métodos numéricos para resolver PFs con condiciones no Dirichlet.
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 5 / 14
Conceptos previos: Disparo para PF no lineales
Conceptos previos:
Disparo para PF no lineales
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 6 / 14
Conceptos previos: Disparo para PF no lineales Condiciones Dirichlet
Disparo para PF no lineales
y ′′ = F (x, y, y ′ ), x ∈ [a, b], y(a) = α, y(b) = β
1. PF −→ PVIs
′′
y = F (x, y, y ′ ), y ′′ = F (x, y, y ′ ),
⇒ y(x, t0 ) ⇒ y(x, t1 )
y(a) = α, y ′ (a) = t0 . y(a) = α, y ′ (a) = t1 .
2. Definimos F (tk ) = y(b, tk ) − β, k = 1, 2, . . .
3. ¿|F (tk )| < tol?
NO Obtenemos un nuevo t = tk+1 por los métodos de la Secante o el de Newton:
F (tk )(tk − tk−1 ) F (tk )
S: tk+1 = tk − , N: tk+1 = tk − ′
F (tk ) − F (tk−1 ) F (tk )
y resolvemos el PVI ′′
y = F (x, y, y ′ ),
y(a) = α, y ′ (a) = tk+1 .
Volvemos a 3.
SÍ Solución: y(x, tk )
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 7 / 14
Laboratorio I
Laboratorio I
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 8 / 14
Laboratorio I Criterios de calificación
Criterios de calificación
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 9 / 14
Laboratorio I Descripción del laboratorio
Problema de frontera
Descripción
Con esta actividad vas a conseguir poner en práctica los conceptos relacionados con
los problemas de frontera unidimensionales estudiados en la asignatura.
Concretamente, se aplicará el método de disparo con Secante y con Newton a un
problema de frontera cuyas condiciones no son de tipo Dirichlet.
Consideramos el siguiente problema de frontera:
y ′′ + y ′ y − 2 = 2xy, x ∈ [1, 2],
1 17
y ′ (1) − y(1) = , y ′ (2) + y(2) = .
2 2
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 10 / 14
Laboratorio I Descripción del laboratorio
Problema de frontera
Apartado (a) (1.5 pts)
Describe el método de disparo para este problema, utilizando el esquema de
secante para determinar los valores del parámetro t.
−→ Los copia y pega de la teoría se penalizarán.
−→ Explicar el método con esquema de Secante para este problema en concreto.
−→ Mínimo: El PVI, la transformación en sistema, la ecuación no lineal a resolver, la expresión
iterativa del método de Secante y porqué lo usamos y el criterio de parada.
Apartado (b) (2 pts)
Aproxima la solución del problema mediante el método del apartado (a),
tomando 10 subintervalos, una tolerancia de 10−8 y un valor inicial del parámetro
t = 0. Indica el número de iteraciones y el último valor de t.
−→ Si se piden 10 subintervalos debemos esperar una salida de la aproximación y con 11
nodos. Si no es así, se penalizará.
−→ t1 lo definimos como se hace normalmente.
−→ Mínimo: Tabla de nodos y aproximación de y e y ′ , número de iteraciones y t.
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 11 / 14
Laboratorio I Descripción del laboratorio
Problema de frontera
Apartado (c) (1.5 pts)
Describe el método de disparo para este problema, utilizando el esquema de
Newton para determinar los valores del parámetro t.
−→ Los copia y pega de la teoría se penalizarán.
−→ Explicar el método con esquema de Newton para este problema en concreto.
−→ Mínimo: Los PVIs a resolver, la transformación en sistema, la ecuación no lineal a resolver,
la expresión iterativa del método de Newton y porqué lo usamos y el criterio de parada.
Apartado (d) (2 pts)
Aproxima la solución del problema mediante el método del apartado (c),
tomando 10 subintervalos, una tolerancia de 10−8 y un valor inicial del parámetro
t = 0. Indica el número de iteraciones y el último valor de t.
−→ Si se piden 10 subintervalos debemos esperar una salida de la aproximación y con 11
nodos. Si no es así, se penalizará.
−→ Mínimo: Tabla de nodos y aproximación de y e y ′ , número de iteraciones y t.
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 12 / 14
Laboratorio I Descripción del laboratorio
Problema de frontera
Apartado (e) (1.5 pts)
Compara las aproximaciones obtenidas en los apartados (b) y (d), tanto gráfica
como numéricamente.
−→ Gráfica (0.75 pts). Mínimo: Título de la gráfica, leyenda y labels.
−→ Tabla (0.75 pts). Mínimo: nodos, aproximación por Secante y aproximación por Newton.
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 13 / 14
Laboratorio I Descripción del laboratorio
Métodos Numéricos Aplicados II Laboratorio 1 Curso 2023-2024 14 / 14