Método de la secante

En análisis numérico el método de la secante es un método para encontrar los ceros de una función de forma iterativa.

Es una variación del método de Newton-Raphson donde en vez de calcular la derivada de la función en el punto de estudio, teniendo en mente la definición de derivada, se aproxima la pendiente a la recta que une la función evaluada en el punto de estudio y en el punto de la iteración anterior. Este método es de especial interés cuando el coste computacional de derivar la función de estudio y evaluarla es demasiado elevado, por lo que el método de Newton no resulta atractivo.

En otras palabras, el método de la secante es un algoritmo de la raíz de investigación que utiliza una serie de raíces de las líneas secantes para aproximar mejor la raíz de una función f. El método de la secante se puede considerar como una aproximación en diferencias finitas del método de Newton-Raphson. Sin embargo, este método fue desarrollado independientemente de este último.

El método se define por la relación de recurrencia:

Como se puede ver, este método necesitará dos aproximaciones iniciales de la raíz para poder inducir una pendiente inicial.

El método se basa en obtener la ecuación de la recta que pasa por los puntos (x_n−1, f(x_n−1)) y (x_n, f(x_n)). A dicha recta se le llama secante por cortar la gráfica de la función. En la imagen de arriba a la derecha se toman los puntos iniciales x₀ y x₁, se construye una línea por los puntos (x₀, f(x₀)) y (x₁, f(x₁)). En forma punto-pendiente, esta línea tiene la ecuación mostrada anteriormente. Posteriormente se escoge como siguiente elemento de la relación de recurrencia, x_n+1, la intersección de la recta secante con el eje de abscisas obteniendo la fórmula, y un nuevo valor. Seguimos este proceso, hasta llegar a un nivel suficientemente alto de precisión (una diferencia lo suficientemente pequeñas entre x_n y x_n-1).

El orden de convergencia de este método, en un punto cercano a la solución, es ${displaystyle varphi }$ donde

es el número áureo, por lo que se trata de una convergencia superlineal inferior a la del método de Newton-Raphson. En caso de que la aproximación inicial sea demasiado lejana o la raíz no sea simple, este método no asegura la convergencia y tiene un comportamiento similar al de Newton-Raphson.

El método de bisección necesita de muchas iteraciones comparado con el método de la secante, ya que el proceso que este sigue es mucho más preciso que el de bisección, el cual solo divide por mitades sucesivamente hasta dar con un valor aproximado al real y por consecuente conlleva un número significativamente mayor de iteraciones.

El método de la regla falsa utiliza la misma fórmula que el método de la secante. Sin embargo, no se aplica la fórmula en x_n−1 y x_n, como el método de la secante, pero en x_n y en la última iteración x_k tal que f(x_k) y f(x_n) tiene un signo diferente. Esto significa que el método de regla falsa siempre converge.

La fórmula de recurrencia del método de la secante se puede derivar de la fórmula para el método de Newton-Raphson:

utilizando la aproximación de diferencias finitas:

Si comparamos el método de Newton-Raphson con el método de la secante, vemos que el método de Newton-Raphson converge más rápido (para 2 en contra α ≈ 1,6). Sin embargo, el método de Newton-Raphson requiere la evaluación de ambos f y su derivada en cada paso, mientras que el método de la secante sólo requiere la evaluación de f. Por lo tanto, el método de la secante puede muy bien ser más rápido en la práctica.

Utilice el método de la secante para encontrar una raíz real de la ecuación polinomial: F(x)=x³+2x²+10x-20=0.

Utilizando la ecuación:

Obtenemos:

Y mediante x₀=0 y x₁=1 se calcula x₂

Los valores posteriores son los siguientes:

Ahí tenemos el resultado, cuando ${displaystyle |x_{n+1}-x_{n}|}$ ≤ ${displaystyle epsilon =10^{-3}}$

Comprobando el resultado graficando la función utilizando software obtenemos:

Si bien no se converge a la raíz tan rápido como resolviéndolo utilizando el método Newton-Raphson, la velocidad de convergencia no es tan lenta como resolviéndolo por el método de punto fijo; entonces se tiene para este ejemplo una velocidad de convergencia intermedia.

Programa escrito en Java correspondiente al ejemplo f(x) = exp( x ) - ( 2 * x^3 )

Programa escrito en Fortran 90 correspondiente al ejemplo f(x) = x³ + 2x² + 10x - 20

Para compilar en GNU/Linux con compilador de GNU, se escribe en una terminal:

Programa escrito en Matlab para ejecutar el método de la secante.

Alternativa al código anterior.

Antonio Nieves, Federico C. Domínguez. "Métodos numéricos aplicados a la ingeniería" (Tercera Edición). Editorial Patria.