Teoria Mayo 2022 Final Solucion Es
Teoria Mayo 2022 Final Solucion Es
Teoria Mayo 2022 Final Solucion Es
Métodos Numéricos
Examen final.
a) (0.5 puntos) Dados a = 1.3524 y b = 1.3528, calcule el punto medio de [a, b] utilizando
aritmética de 4 dı́gitos (a.k.a. aritmética de coma fija con 3 decimales) obtenidos mediante
redondeo. Comente el resultado.
b) (0.5 puntos) Resuelva la ecuación de segundo grado
x2 + 1016 x + 1 = 0
De la función f (x) = 12 e2x − 9x. se sabe que la ecuación f (x) = 0 tiene exactamente dos soluciones
reales c1 y c2 tales que c1 ∈ [0, 1/2] y c2 ∈ [3/2, 2].
a) (1 punto) Demuestre que es posible aproximar la solución c1 mediante el método de punto fijo
eligiendo como función de iteración g(x) = e2x /18 y como semilla cualquier punto del intervalo
[0, 1/2]. ¿Cuántas iteraciones son necesarias para asegurar que dicha raı́z c1 se aproxima con
un error absoluto menor que 10−4 tomando como semilla x0 = 0?
b) (1 punto) Aproxime la solución c2 mediante una iteración del metodo de la secante tomando
como semillas x0 = 3/2 y x1 la segunda iterada del metodo de bisección.
Solución:
a) Tenemos que probar que hay convergencia global del método del punto fijo en el intervalo
[0, 1/2] con función de iteración g(x).
Primero, es necesario comprobar que, en [0, 1/2], el punto fijo de g(x) coincide con la raı́z de
f (x).
1
g(x) = x ⇔ e2x = 18x ⇔ e2x = 9x ⇔ f (x) = 0
2
Notar que g(x) es indefinidamente derivable en todo R y con derivadas continuas, en particular,
g ′ (x) = e2x /9 > 0.
Tenemos que comprobar también la condición de mapeo g([0, 1/2]) ⊂ [0, 1/2] y la condición
sobre la derivada |g ′ (x)| ≤ L < 1 en [0, 1/2]. Siendo g monótona es suficiente verificar la
condición de mapeo en los extremos del intervalo:
Y se comprueba que g(0), g(1/2) ∈ [0, 1/2]. Para la condicón sobre la derivada tenemos que
y L ≈ 0.8210 < 1.
Por lo tanto, por el teorema de convergencia global el punto fijo, se puede garantizar que g
tiene un único punto fijo en [0, 1/2] y que el método iterativo de punto fijo converge a dicho
punto fijo para cualquier aproximación inicial x0 ∈ [0, 1/2].
Ln
|xn − c1 | ≤ |x1 − x0 |
1−L
donde x0 = 0, x1 = g(x0 ) = 1/18, L = e2 /9 ≈ 0.8210. Entonces se tiene que cumplir
Ln
|x1 − x0 | < 10−4
1−L
de donde n > 40.76. Por lo tanto son necesarias 41 iteraciones.
b) Vamos a calcular la semilla x1 mediante dos iteraciones del método de bisección. Ponemos
[a0 , b0 ] = [3/2, 2] = [1.5, 2] y llamaremos x̃n las aproximaciones obtenidas por el método de
bisección.
Se inicializa a0 = 1.5 y b0 = 2,
Para n = 0, x̃1 = a0 +b2
0
= 1.5+2
2 = 1.75 y f (x̃1 ) = 0.8077 > 0
Como f (a0 )f (x1 ) < 0, se toma a1 = a0 = 1.5 y b1 = x1 = 1.75
Para n = 1, x̃2 = a1 +b2
1
= 1.5+1.75
2 = 1.6250
La primara iteración del método de la secante es
x1 − x0
x2 = x1 − f (x1 )
f (x1 ) − f (x0 )
2
3. Ejercicio Tema 3 (2 puntos)
Solución:
donde P A = LU , con
1 0 0 1 0 0
P = 0 0 1 , L= 0 1 0
0 1 0 1/7 16/21 1
7 5 1
y U = 0 3 6
0 0 −26/7
Para obtener un intervalo del error relativo cometido se tiene en cuenta la siguiente desigualdad:
3
Tomando normas (||r||∞ = 0.09, ||b||∞ = 3) resulta que:
es estrictamente diagonal dominante por filas, y, por lo tanto, los métodos de Jacobi y Gauss-
Seidel son convergentes cuando se aplican a la resolución de cualquier sistema de la forma
Ax = b.
a) (0.5 puntos) Primero pensamos en usar cuatro nodos equiespaciados. Escriba explı́citamente los
nodos y obtenga una cota razonable del máximo error de interpolación si usamos el polinomio
de interpolación de Lagrange.
b) Cambiamos de opinión y definimos los tres nodos xi = (i − 2)π/4 para i = 1, 2, 3.
i) (0.25 puntos) Calcule la base de Lagrange relativa a los nodos x1 , x2 y x3 .
ii) (0.25 puntos) Use esta base para obtener una expresión para p2 (x), el polinomio de inter-
polación de Lagrange de f (x) en estos tres nodos.
c) (1 punto) Por último, abandonamos la interpolación. Encuentre a ∈ R tal que la función
y = a · (2 − cosh x) ajuste los datos (xi , cos xi )i=1,2,3 del apartado (b) en el sentido de los
mı́nimos cuadrados.
Nota: La función cosh x suele estar disponible de forma directa en la mayorı́a de las calculadoras. En
cualquier caso, se puede calcular como cosh x = 0.5(e−x +ex ). Use al menos 5 cifras en todos los cálculos
intermedios y tres en la respuesta final. Por ejemplo, cosh(π/4) = 1.3246.
Solución:
a) Cuatro nodos significa n = 3 huecos. Como la longitud del intervalo es π, tenemos que h = π/3.
Por tanto,
π π π π
x1 = − , x2 = − , x3 = , x4 = .
2 6 6 2
Como la función es el coseno, trivialmente M4 = 1, ası́ que una cota del error es
M4 4 1 π 4
|E3 (x)| ≤ h = ≈ 5.01 × 10−2 .
24 24 3
x(x − π/4) 8
ℓ1 (x) = = 2 x(x − π/4).
−π/4(−π/4 − π/4) π
4
√
ii) En esta base, y usando que y1 = y3 = 2/2 e y2 = 1, el polinomio de interpolación de
Lagrange es
√ √
4 2 16 4 2
p2 (x) = 2 x(x − π/4) − 2 (x + π/4)(x − π/4) + 2 (x + π/4)x.
π π π
c) Tenemos que resolver en el sentido de los mı́nimos cuadrados el sistema incompatible
√
(2 − cosh(−π/4))a = 2/2
(2 − cosh(0))a = 1
√
(2 − cosh(π/4))a = 2/2
(Nótese que, aunque las ecuaciones primera y tercera son la misma, ¡no podemos deshacernos
de ninguna de ellas!). Este es un sistema de la forma Xa = y, donde
√
2 − cosh(−π/4) 0.6754 2/2 0.7071
X= 1 ≈ 1 , y = √1 ≈ 1
2 − cosh(π/4) 0.6754 2/2 0.7071
La solución en el sentido de los mı́nimos cuadrados de este sistema es la única solución de las
ecuaciones normales X T Xa = X T y. En este caso:
X T X = 0.67542 + 1 + 0.67542 = (1.9123),
a) (0.75 puntos) Una forma de deducir fórmulas de derivación numérica es calcular un polinomio
interpolador y aproximar la derivada que nos interese mediante el la derivada del polinomio
interpolador.
Sea f (x) una función de clase C 2 en un entorno de un punto c ∈ R y consideremos un tamaño
de paso h > 0. Calcule p2 (x), el polinomio interpolador de f en los nodos c − h, c y c + h y
deduzca una fórmula para aproximar f ′′ (c) usando p′′2 (c).
b) (0.75 puntos) En el intervalo [−1, 1], obtenga la fórmula de Newton-Côtes cerrada de 4 nodos.
¿Cuál es su grado de precisión?
Solution.
a) Para calcular el polinomio de interpolación, usaremos la forma de Newton. Los nodos son
c − h, c, c + h y el árbol de diferencias divididas es
c − h f (c − h)
f (c)−f (c−h)
c f (c) h
f (c+h)−f (c) f (c+h)−2f (c)+f (c−h)
c + h f (c + h) h 2h2
Luego
f (c) − f (c − h) f (c + h) − 2f (c) + f (c − h)
p2 (x) = f (c − h) + (x − c + h) + (x − c + h)(x − c)
h 2h2
f (c+h)−2f (c)+f (c−h)
Como el único factor de x2 es 2h2
, tenemos que
f (c + h) − 2f (c) + f (c − h)
f ′′ (c) ≈ p′′2 (c) =
h2
5
b) Como la fórmula es cerrada y los nodos son equiespaciados, tenemos que h = 2/3, y los nodos
son t1 = −1, t2 = −1/3, t3 = 1/3 y t4 = 1. Como la fórmula de Newton-Côtes es de tipo
interpolatorio, si grado de precisión será al menos 3, ası́ que para calcular los pesos imponemos
que sea exacta para f (t) = tn , con n = 0, 1, 2, 3.
Z 1
1dt = 2 = ω1 + ω2 + ω3 + ω4
−1
Z 1
1 1
tdt = 0 = −ω1 − ω2 + ω3 + ω4 .
−1 3 3
Z 1
2 1 1
t2 dt = = ω1 + ω2 + ω3 + ω4 .
−1 3 9 9
Z 1
1 1
t3 dt = 0 = ω1 − ω2 + ω3 + ω4 .
−1 27 27
Tenemos un sistema de 4 ecuaciones con 4 incógnitas. Podemos aplicar, por ejemplo, el método
de Gauss para resolverlo. Otra posibilidad es fijarse en que las ecuaciones segunda y cuarta
son
1 1
−ω1 − ω2 + ω3 + ω4 = 0
3 3
1 1
−ω1 − ω2 + ω3 + ω4 = 0
27 27
Restando ambas ecuaciones tenemos que
8 8
− ω2 + ω3 = 0
27 27
y por tanto ω2 = ω3 . Multiplicando la segunda por 9 y restando, se tiene que 8ω1 − 8ω4 = 0,
luego ω1 = ω4 . (Esto pasa en general y se suele decir que para nodos simétricos los pesos
coinciden). Aprovechando esto en las ecuaciones primera y tercera obtenemos
ω1 + ω2 = 1
1 1
ω1 + ω2 =
9 3
de donde ω1 = 41 , ω2 = 3
4 y la fórmula es
1 3 3 1
Q(f, −1, 1) = f (−1) + f (−1/3) + f (1/3) + f (1)
4 4 4 4
Veamos si es exacta para f (t) = t4 . Por un lado
Z 1
2
t4 dt = .
−1 5
Por el otro
1 3 1 3 1 1 14
Q(t4 , −1, 1) =
+ + + = .
4 4 81 4 81 4 27
Luego el grado de precisión es exáctamente 3.
6
a) (0.75 puntos) Aproxime y(1) usando el método de Euler con dos pasos.
b) (0.75 puntos) Usando el método de Heun con dos pasos se ha obtenido la aproximación y(t2 ) ≈
y2H = 5.7267. Utilice este resultado para estimar el error absoluto cometido al aproximar y(1)
con el método de Euler con dos pasos. ¿Cuántas iteraciones del método de Euler utilizarı́a
para conseguir que el error absoluto sea menor que 10−3 ?
Solución.
a) En el problema de valor inicial del enunciado la ecuación diferencial ordinaria es de orden uno,
y podemos escribirlo en la forma
′
y (t) = f (t, y(t)) en (0, 1]
y(0) = y0
donde f (t, y) = −0.6 y 1/3 y y0 = 8. Vamos a utilizar N = 2 pasos para aproximar la solución
en T = 1, de forma que t0 = 0, t1 = t0 + h y t2 = t0 + 2h = 1; para que esto último se cumpla,
necesitamos tomar el paso h = 1−0 2 = 0.5 y entonces t0 = 0, t1 = 0.5 y t2 = 1.
Si aplicamos el método de Euler, obtenemos
y(t0 ) = y0E = y0 = 8;
y(t1 ) ≈ y1E = y0E + f (t0 , y0E ) = 8 − 0.6 81/3 = 8 − 0.6 2 = 6.8;
y(t2 ) ≈ y2E = y1E + f (t1 , y1E ) = 6.8 − 0.6 6.81/3 ≈ 5.6633.
b) Como el método de Heun es de orden mayor que el de Euler, podemos usar y(t2 ) ≈ y2H para
estimar el error cometido en la aproximaciı́on y(t2 ) ≈ y2E . De esta forma,
Como el método de Euler es de orden 1, esperamos que el error absoluto se comporte lineal-
mente con respecto al paso h; i.e. el error absoluto cometido al aproximar nos sugiere tomar la
constante C = 0.1268. En consecuencia, para garantizar que g(h) = 0.1268 h < 10−3 usarı́amos
h = 1/N < 7.8864 10−3 , es decir, N > 126.801. Si optamos por el menor valor que lo satisface,
necesitamos N = 127 pasos.