Raíces de ecuaciones · Definición
Iteración de punto fijo
Antes de empezar necesitas conocer
Muchas ecuaciones no lineales pueden reescribirse como un punto fijo de una función. Este enfoque convierte el problema de encontrar raíces en una sucesión que debe converger.
Definición
Método iterativo que aproxima una solución de la ecuación \( x = g(x) \) generando la sucesión \( x_{n+1} = g(x_n) \) a partir de un valor inicial \( x_0 \).
Un punto fijo \( p \) de \( g \) cumple \( g(p) = p \). Si \( g \) es continua y la sucesión \( x_n \) converge, su límite es un punto fijo. El teorema del punto fijo establece que, si \( g \) envía \( [a,b] \) en sí mismo y \( |g'(x)| \le k < 1 \) en ese intervalo, existe un único punto fijo y la iteración converge a él para cualquier \( x_0 \in [a,b] \). La convergencia es lineal con constante \( |g'(p)| \) cuando \( g'(p) \neq 0 \).
Para resolver \( x^2 - 2 = 0 \) puede usarse \( g(x) = \frac{x + 2}{x + 1} \) porque \( g(x) = x \) equivale a \( x^2 = 2 \). Con \( x_0 = 1 \) se obtiene \( x_1 = 1.5 \), \( x_2 \approx 1.4 \) y \( x_3 \approx 1.4167 \), que se acerca a \( \sqrt{2} \approx 1.4142 \). Cada paso se calcula evaluando \( g \) en el valor anterior.
Un error típico es elegir una \( g \) con \( |g'(p)| \ge 1 \) y esperar convergencia. Por ejemplo, \( g(x) = x^2 - 2 \) tiene puntos fijos en \( x = 2 \) y \( x = -1 \), pero \( |g'(2)| = 4 \), así que la iteración desde valores cercanos a 2 se aleja de la raíz. Antes de iterar conviene verificar la condición sobre la derivada.
Para pensar
Si \( |g'(p)| > 1 \) y \( x_0 \) está muy cerca de \( p \), ¿qué ocurre con la iteración?
El error se multiplica por \( |g'(p)| > 1 \) en cada paso, así que la sucesión se aleja de \( p \), salvo que \( x_0 = p \).
Ejercicios 10
- ¿Cuál es la fórmula general de la iteración de punto fijo?
- ¿Qué condición se pide sobre \( g \) para que se cumpla \( g(p) = p \) en un punto fijo \( p \)?
- Si \( |g'(p)| < 1 \), ¿qué tipo de convergencia se espera cerca de \( p \)?
- ¿Qué establece el teorema del punto fijo sobre la existencia de un punto fijo?
- ¿Cuál es el valor de \( x_1 \) si \( g(x) = \cos(x) \) y \( x_0 = 0 \)?
- Si \( g(x) = \frac{x^2 + 2}{3} \) y \( x_0 = 1 \), ¿cuál es \( x_2 \)?
- ¿Cuál es la constante de convergencia lineal cuando \( g'(p) = -0.5 \)?
- Para \( g(x) = \sqrt{x + 6} \), ¿cuál es el punto fijo positivo?
- Si \( g(x) = \frac{x}{2} + 1 \), ¿cuál es el límite de \( x_n \) partiendo de cualquier \( x_0 \)?
- Se quiere resolver \( x^3 - x - 1 = 0 \) con \( g(x) = \sqrt[3]{x + 1} \). ¿Converge la iteración desde \( x_0 = 1 \)? Justifica con \( |g'(p)| \) en la raíz real \( p \approx 1.3247 \).
Después de esto
Referencias
- Richard L. Burden, J. Douglas Faires, Annette M. Burden, Numerical Analysis, 10th edition, Cengage Learning, 2016.
- Lloyd N. Trefethen, David Bau III, Numerical Linear Algebra, 1st edition, SIAM, 1997.