Deducción: convergencia y orden del punto fijo

Desarrollar ϕ\phi por Taylor alrededor del punto fijo produce la ecuación del error del método y demuestra a la vez el criterio ϕ(α)<1|\phi'(\alpha)|<1 y el teorema del orden.

Taylor de φ alrededor del punto fijo

  1. Sea ek=xkαe_k=x_k-\alpha. Desarrollamos ϕ(xk)\phi(x_k) por Taylor en torno a α\alpha y usamos que ϕ(α)=α\phi(\alpha)=\alpha:

    xk+1=ϕ(xk)=α+ϕ(α)ek+ϕ(α)2ek2+ϕ(α)6ek3+x_{k+1}=\phi(x_k)=\alpha+\phi'(\alpha)\,e_k+\frac{\phi''(\alpha)}{2}\,e_k^2+\frac{\phi'''(\alpha)}{6}\,e_k^3+\cdots
  2. Restando α\alpha queda la ecuación del error de la iteración de punto fijo:

    ek+1=ϕ(α)ek+ϕ(α)2ek2+ϕ(α)6ek3+e_{k+1}=\phi'(\alpha)\,e_k+\frac{\phi''(\alpha)}{2}\,e_k^2+\frac{\phi'''(\alpha)}{6}\,e_k^3+\cdots
  3. Si ϕ(α)0\phi'(\alpha)\ne 0, el término dominante es lineal: ek+1ϕ(α)eke_{k+1}\approx\phi'(\alpha)e_k. Los errores se contraen si ϕ(α)<1|\phi'(\alpha)|<1 (convergencia lineal con factor ϕ(α)|\phi'(\alpha)|) y crecen si ϕ(α)>1|\phi'(\alpha)|>1: ese es el criterio de convergencia local.

  4. Si además ϕ(α)=ϕ(α)==ϕ(p1)(α)=0\phi'(\alpha)=\phi''(\alpha)=\dots=\phi^{(p-1)}(\alpha)=0 y ϕ(p)(α)0\phi^{(p)}(\alpha)\ne 0, todos los términos anteriores desaparecen y el primero superviviente fija el orden: es el teorema del orden de un método de punto fijo.

    ek+1=ϕ(p)(α)p!ekp+O(ekp+1)e_{k+1}=\frac{\phi^{(p)}(\alpha)}{p!}\,e_k^p+\mathcal{O}\bigl(e_k^{p+1}\bigr)
  5. Aplicación inmediata: para Newton, ϕ=xff\phi=x-\frac{f}{f'} y ϕ=ff(f)2\phi'=\frac{f f''}{(f')^2}, que se anula en una raíz simple (f(α)=0f(\alpha)=0). Como en general ϕ(α)0\phi''(\alpha)\ne 0, Newton tiene orden 2, en acuerdo con la demostración directa.