Articulo de referencia

Método de Ridders

En análisis numérico , el método de Ridders es un algoritmo de búsqueda de raíces basado en el método de la falsa posición y el uso de una función exponencial para aproximar suc...

En análisis numérico , el método de Ridders es un algoritmo de búsqueda de raíces basado en el método de la falsa posición y el uso de una función exponencial para aproximar sucesivamente una raíz de una función continua.F(incógnita){\displaystyle f(x)}El método se debe a C.  Ridders. [ 1 ] [ 2 ]

El método de Ridders es más simple que el método de Muller o el método de Brent , pero con un rendimiento similar. [ 3 ] La fórmula siguiente converge cuadráticamente cuando la función se comporta bien, lo que implica que el número de dígitos significativos adicionales encontrados en cada paso se duplica aproximadamente; pero la función debe evaluarse dos veces para cada paso, por lo que el orden general de convergencia del método con respecto a las evaluaciones de la función en lugar de con respecto al número de iteraciones es2{\displaystyle {\sqrt {2}}}Si la función no se comporta correctamente, la raíz permanece entre corchetes y la longitud del intervalo de corchetes se reduce al menos a la mitad en cada iteración, por lo que se garantiza la convergencia.

Método

Dados dos valores de la variable independiente,incógnita0{\displaystyle x_{0}}yincógnita2{\displaystyle x_{2}}, que se encuentran en dos lados diferentes de la raíz que se busca, de modo queF(incógnita0)F(incógnita2)<0{\displaystyle f(x_{0})f(x_{2})<0}El método comienza evaluando la función en el punto medio. incógnita1=(incógnita0+incógnita2)/2{\displaystyle x_{1}=(x_{0}+x_{2})/2}Entonces se encuentra la función exponencial única.miaincógnita{\displaystyle e^{ax}}de tal manera que funcioneh(incógnita)=F(incógnita)miaincógnita{\displaystyle h(x)=f(x)e^{ax}}Satisfaceh(incógnita1)=(h(incógnita0)+h(incógnita2))/2{\displaystyle h(x_{1})=(h(x_{0})+h(x_{2}))/2}. Específicamente, el parámetroa{\displaystyle a}está determinado por

mia(incógnita1incógnita0)=F(incógnita1)firmar[F(incógnita0)]F(incógnita1)2F(incógnita0)F(incógnita2)F(incógnita2).{\displaystyle e^{a(x_{1}-x_{0})}={\frac {f(x_{1})-\operatorname {sign} [f(x_{0})]{\sqrt {f(x_{1})^{2}-f(x_{0})f(x_{2})}}}{f(x_{2})}}.}

Luego se aplica el método de posición falsa a los puntos.(incógnita0,h(incógnita0)){\displaystyle (x_{0},h(x_{0}))}y(incógnita2,h(incógnita2)){\displaystyle (x_{2},h(x_{2}))}, lo que conduce a un nuevo valorincógnita3{\displaystyle x_{3}}entreincógnita0{\displaystyle x_{0}}yincógnita2{\displaystyle x_{2}},

incógnita3=incógnita1+(incógnita1incógnita0)firmar[F(incógnita0)]F(incógnita1)F(incógnita1)2F(incógnita0)F(incógnita2),{\displaystyle x_{3}=x_{1}+(x_{1}-x_{0}){\frac {\operatorname {sign} [f(x_{0})]f(x_{1})}{\sqrt {f(x_{1})^{2}-f(x_{0})f(x_{2})}}},}

que se utilizará como uno de los dos valores de delimitación en el siguiente paso de la iteración. El otro valor de delimitación se toma comoincógnita1{\displaystyle x_{1}}siF(incógnita1)F(incógnita3)<0{\displaystyle f(x_{1})f(x_{3})<0}(lo cual será cierto en el caso de buen comportamiento), o de lo contrario, cualquiera de las siguientes opciones:incógnita0{\displaystyle x_{0}}yincógnita2{\displaystyle x_{2}}tiene un valor de función de signo opuesto aF(incógnita3).{\displaystyle f(x_{3}).}El procedimiento iterativo puede finalizarse cuando se alcanza la precisión deseada.

El método de Ridders es un algoritmo numérico para encontrar raíces que se utiliza para hallar un cero de una función continua. Es un método iterativo que combina el método de la falsa posición con una función exponencial para converger hacia una raíz. Generalmente se considera que este método es más sencillo y su rendimiento es similar al de otros métodos populares como los de Muller o Brent.

Referencias

  1. Ridders, C. (1979). "Un nuevo algoritmo para calcular una sola raíz de una función continua real". IEEE Transactions on Circuits and Systems . 26 (11): 979– 980. doi : 10.1109/TCS.1979.1084580 .
  2. Kiusalaas, Jaan (2010). Métodos numéricos en ingeniería con Python (2.ª ed.). Cambridge University Press. págs. 146–150 . ISBN   978-0-521-19132-6.
  3. Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). «Sección 9.2.1. Método de Ridders» . Numerical Recipes : The Art of Scientific Computing (3.ª ed.). Nueva York: Cambridge University Press. ISBN  978-0-521-88068-8.