C.4 El método secante
( \newcommand{\kernel}{\mathrm{null}\,}\)
Dejarf(x) ser una función continua. El método secante es una variante del método de Newton que evita el uso de la derivada def(x) — lo cual puede ser muy útil cuando se trata de la derivada no es fácil. Evita el uso de la derivada aproximandof′(x) porf(x+h)−f(x)h para algunos Esh. decir, aproxima la línea tangente af atx por una línea secante paraf que pasa porx. Para limitar el número de evaluaciones def(x) requerido, utiliza x=xn−1yx+h=xn. así es como funciona.
Supongamos que ya hemos encontradoxn. Entonces denotamos pory=F(x) la ecuación de la línea (secante) que pasa por(xn−1,f(xn−1))(xn,f(xn)) y y elegimosxn+1 ser el valor dex dondeF(x)=0.
La ecuación de la línea secante es
y=F(x)=f(xn−1)+f(xn)−f(xn−1)xn−xn−1(x−xn−1)
por lo quexn+1 está determinado por
\ begin {alinear*} & 0=F (x_ {n+1}) = f (x_ {n-1}) +\ frac {f (x_n) -f (x_ {n-1})} {x_n-x_ {n-1}} (x_ {n+1} -x_ {n-1})\\ &\ iff x_ {n+1} = x_ {n-1}} -\ frac {x_n-x_ {n-1}} {f (x_n) -f (x_ {n-1})} f (x_ {n-1})\ final {alinear*}
o, simplificando,
xn+1=xn−1f(xn)−xnf(xn−1)f(xn)−f(xn−1)
Por supuesto, para empezarn=1, necesitamos dos conjeturas iniciales,x0 yx1, para la raíz.
En este ejemplo calculamos, aproximadamente, la raíz cuadrada de dos aplicando el método secante a la ecuación
f(x)=x2−2=0
y compararemos los resultados del método secante con los resultados correspondientes del método de Newton. (Ver Ejemplo C.1.2.)
Ya quef′(x)=2x, (C.1.1) dice que, bajo el método de Newton, debemos aplicar iterativamente
xn+1=xn−f(xn)f′(xn)=xn−x2n−22xn=xn2+1xn
mientras que (C.4.1) dice que, bajo el método secante, debemos aplicar iterativamente (después de simplificar un poco el álgebra)
\ begin {alinear*} x_ {n+1} & =\ frac {x_ {n-1} f (x_n) - x_n f (x_ {n-1})} {f (x_n) -f (x_ {n-1})} =\ frac {x_ {n-1} [x_n^2-2] - x_n [x_ {n-1} ^2-2]} {x_n^2-x_ {n-1} ^2}\\ & =\ frac {x_ {n-1} x_n [x_n-x_ {n-1}] +2 [x_n-x_ {n-1}]} {x_n^2-x_ {n-1} ^2}\\ & =\ frac {x_ {n-1} x_n+2} {x_ {n-1} +x_n}\ final {alinear*}
Aquí están los resultados, comenzando con el método de Newtonx1=4 e iniciando el método secante conx0=4,x1=3. (Entonces le estamos dando un poco de ventaja al método secante).
\ begin {alignat*} {3} & & &\ text {método secante}\ qquad & &\ text {método de Newton}\\ & x_0\ quad & & & 4 & &\\ & x_1 & & & 3 & 4\ & x_2 & & 2 & & & 2.25\ & x_3 & & 1.6 & & 1.57\ & x_4 & & 1.444 & amp; 1.422\\ & x_5 & & & 1.4161 & & 1.414234\\ & x_6 & & & 1.414233 & & 1.414213562525\\ & x_7 & & 1.414213575 & & 1.414213562373095\ end {alignat*}
Para fines de comparación, la raíz cuadrada de2, a 15 decimales, es1.414213562373095. Así que el método secantex7 es preciso a 7 lugares decimales y el método de Newtonx7 es exacto a al menos 15 lugares decimales.
La ventaja que tiene el método secante sobre el método de Newton es que no utiliza la derivada def. Esto puede ser una ventaja sustancial, por ejemplo cuando la evaluación de la derivada es computacionalmente difícil o costosa. Por otro lado, el ejemplo anterior sugiere que el método secante no es tan rápido como el método de Newton. El siguiente subapartado muestra que efectivamente así es.