Método de Steffensen

El método de Steffensen (por Johan Frederik Steffensen) es un algoritmo para obtener los ceros de una función. El método de Steffensen se puede considerar como una combinación del método de punto fijo y del método de Aitken. Como el método de Aitken esencialmente acelera la convergencia de otro método, se puede definir este método como el método de punto fijo acelerado.

Ventajas[editar]

El método de Steffensen presenta una convergencia rápida y no requiere, como en el caso del método de Newton, la evaluación de derivada alguna. Presenta además, la ventaja adicional de que el proceso de iteración sólo necesita un punto inicial.

Otra ventaja del método de Steffensen es que -al igual que el de Newton- tiene convergencia cuadrática. Es decir, ambos métodos permiten encontrar las raíces de una función f "rápidamente" - en este caso rápidamente significa que en cada iteración, el número de dígitos correctos en la respuesta se duplica. Pero la fórmula para el método de Newton requiere la evaluación de la derivada de la función, el método de Steffensen no, por lo que este último puede ser programado para una función genérica, mientras que la función cumpla la restricción mencionada anteriormente.

El precio de la convergencia rápida es una doble evaluación de la función: tanto $f(x_{n})$ como $f(x_{n}+h)$ deben ser calculadas, lo que podría llevar un tiempo considerable dependiendo de la función f. Por comparación, el método de la secante sólo necesita una evaluación de la función por cada paso, así que con dos evaluaciones de la función del método de la secante se pueden hacer dos pasos, y esos dos pasos aumentan el número de dígitos correctos en un factor de 1,6. En un solo paso de tiempo el método de Steffensen (o de Newton) aumenta los dígitos correctos en un factor de 2, lo que es sólo un poco mejor.

Al igual que el método de Newton y otros métodos cuadráticamente convergentes, la debilidad fundamental en el método de Steffensen es la elección del valor inicial $x_{0}$ . Si el valor de $x_{0}$ no está "lo suficientemente cerca" de la solución, el método puede fallar y la secuencia de valores $x_{0},x_{1},x_{2},x_{3},\dots$ o bien puede oscilar entre dos valores, o bien divergir hacia infinito (ambas alternativas pueden suceder).

Se calcula el siguiente punto de iteración a partir de la expresión:

x_{n+1}=x_{n}-{[f(x_{n})]^{2} \over {f(x_{n}+f(x_{n}))-f(x_{n})}}

Algoritmo[editar]

Para una sucesión {x_n}, obtenida por el método del punto fijo x_n+1 = f(x_n), partimos de tres puntos:^{[cita requerida]}

y₀= f(x₀)

z₀= f(y₀)

donde x₀ es el punto inicial. Obteniendo así:

x₁ = x₀ – (y₀ –x₀)^1/2

z₀ – 2*y₀ – x₀

En forma general:

X_n+1 = x_n – (y_n – x_n)^1/2

z_n – 2* y_n – x_n

Donde si |x_n+1 – x_n| = error < Tol entonces se satisface la tolerancia.

Código Matlab[editar]

function [x,n_iter]=steffensen(x0,f,tol,nmax) 
x1=feval(f,x0);
x2=feval(f,x1);
x_nuevo=x0-(x1-x0)^2/(x2-2*x1+x0)
error= abs(x2-x_nuevo);
i=1;
	while(i<=nmax & error>=tol)
		x0=x_nuevo;
		x1=feval(f,x0);
        x2=feval(f,x1);
		x_nuevo=x0-(x1-x0)^2/(x2-2*x1+x0);
		error= abs(x0-x_nuevo);
		i=i+1;
	end
	n_iter=i-1;
	x=x0;
end