MTODOS NUMRICOS Solucin de una ecuacin no lineal

MÉTODOS NUMÉRICOS Solución de una ecuación no lineal “Cálculo de raíces” Profesor: José Andrés Vázquez Flores Otoño 2020

Búsqueda de raíces Consiste en obtener una raíz x de una ecuación de la forma f(x) = 0 para una función f dada. Al número x se le llama también cero de f.

Algoritmo de bisección o método de búsqueda binaria Supongamos que f es una función continua definida en el intervalo [a, b] con f(a) y f(b) de signos diferentes. De acuerdo con el teorema de valor intermedio, existe un número p en (a, b) tal que f(p) = 0. Este procedimiento se aplica en el caso en que f(a) y f(p) tengan signos diferentes y exista más de una raíz en el intervalo (a, b), por razones de simplicidad suponemos que la raíz de este intervalo es única. El método requiere dividir varias veces a la mitad los subintervalos de [a, b] y en cada paso, localizar la mitad que contenga a p.

Algoritmo de bisección o método de búsqueda binaria Para empezar, supongamos que a 1=a y b 1=b, y sea p 1 el punto medio de [a, b] es decir; p 1 = (1/2)(a 1 + b 1) Si f(p 1) = 0, entonces p = p 1; de no ser así, entonces f(p 1) tiene el mismo signo que f(a 1) o f(b 1). Si f(a 1) y f(p 1) tienen el mismo signo, entonces p ϵ (p 1, b 1) y tomamos a 2=p 1 y b 2=b 1. Si f(a 1) y f(p 1) tienen signos opuestos, entonces p ϵ (a 1, p 1) y tomamos a 2=a 1 y b 2=p 1. Después volvemos a aplicar este proceso al intervalo [a 2, b 2]. Esto nos da el siguiente algoritmo.

Algoritmo de bisección o método de búsqueda binaria Para encontrar una solución de f(x)=0 dada la función continua f en el intervalo [a, b], donde f(a) y f(b) tienen signos opuestos: ENTRADA extremos a, b; tolerancia TOL; máximo numero de iteraciones No. SALIDA solución aproximada p o mensaje de fracaso. Paso 1 Tomar i = 1 Paso 2 Mientras que i ≤ No seguir Pasos 3 -6 Paso 3 Tomar p = a + (b - a) / 2 (Calcular pi) Paso 4 Si f(p) = 0 ó (b - a) / 2 < TOL entonces SALIDA (p); (Procedimiento completado satisfactoriamente) PARAR Paso 5 Tomar i = i + 1 Paso 6 Si f(a)f(p) > 0 entonces tomar a = p (Calcular ai, bi) si no tomar b = p Paso 7 SALIDA (“El método fracaso después de No iteraciones, No =“, No) (Procedimiento completado sin éxito) PARAR

Algoritmo de bisección o método de búsqueda binaria

Algoritmo de bisección o método de búsqueda binaria f(x) = x 3 + 4 x 2 – 10 n 1 2 3 4 5 6 7 8 9 10 11 12 13 an 1. 0 1. 25 1. 3125 1. 34375 1. 359375 1. 36328125 1. 364257813 1. 364746094 1. 364990235 [ä, b] = [1, 2] bn 2. 0 1. 5 1. 375 1. 3671875 1. 365234375 │p n-1 – pn │ / │pn │ <10 -4 pn 1. 5 1. 25 1. 375 1. 3125 1. 34375 1. 359375 1. 3671875 1. 36328125 1. 365234375 1. 364257813 1. 364746094 1. 364990235 1. 365112305 f(pn) 2. 375 -1. 79687 0. 16211 -0. 84839 -0. 35098 -0. 09641 0. 03236 -0. 03215 0. 000072 -0. 01605 -0. 00799 -0. 00396 -0. 00194

MÉTODO DE NEWTON RAPHSON f(x) x

MÉTODO DE NEWTON RAPHSON Consiste en elegir un punto inicial cualquiera x 1 como aproximación de la raíz.

MÉTODO DE NEWTON RAPHSON f(x) f(x 1) x 1 x

MÉTODO DE NEWTON RAPHSON Consiste en elegir un punto inicial cualquiera x 1 como aproximación de la raíz y obtener el valor de la función por ese punto. Trazar una recta tangente a la función por ese punto.

MÉTODO DE NEWTON RAPHSON f(x) f(x 1) x 1 x 2 x

MÉTODO DE NEWTON RAPHSON Consiste en elegir un punto inicial cualquiera x 1 como aproximación de la raíz. Obtener el valor de la función por ese punto y trazar una recta tangente a la función por ese punto. El punto de intersección de esta recta con el eje de las abscisas (xr, 0), constituye una segunda aproximación de la raíz.

MÉTODO DE NEWTON RAPHSON f(x) f(x 1) f(x 2) x 1 x 2 x

MÉTODO DE NEWTON RAPHSON Consiste en elegir un punto inicial cualquiera x 1 como aproximación de la raíz. Obtener el valor de la función por ese punto y trazar una recta tangente a la función por ese punto. El punto de intersección de esta recta con el eje de las abscisas (xr, 0), constituye una segunda aproximación de la raíz. El proceso se repite n veces hasta que el punto de intersección xn coincide prácticamente con el valor exacto de la raíz.

MÉTODO DE NEWTON RAPHSON f(x) f(x 1) f(x 2) x 1 x 2 x

MÉTODO DE NEWTON RAPHSON Aunque el método trabaja bien, no existe garantía de convergencia.

MÉTODO DE NEWTON RAPHSON Supóngase que la función f es continuamente diferenciable dos veces en el intervalo [a, b]; o sea f ε C 2[a, b]. Sea ẋ ε [a, b] una aprox. A p ɟ f´ (ẋ)≠ 0 y | ẋ-p| es “pequeño”. Consideremos el Pol. De Taylor de primer grado para f(x) alrededor de ẋ. donde �� (x) esta entre x y ẋ. Dado que f (p)=0 , con x=p tenemos:

MÉTODO DE NEWTON RAPHSON Derivando el método de Newton suponiendo que, como |p- ẋ| es tan pequeño el termino que contiene (p- ẋ)2 es mucho menor entonces podemos tomar: Despejando p de esta ecuación obtenemos: Esto prepara el método de Newton-Raphson, el cual comienza con una aproximación inicial p 0 y genera la sucesión {pn} definida como para n ≥ 1

ALGORITMO DE NEWTON RAPHSON Para encontrar una solución de f(x)=0 dada una aproximación inicial p 0. ENTRADA aproximación inicial p 0 ; tolerancia TOL; máximo numero de iteraciones No. SALIDA solución aproximada p o mensaje de fracaso. Paso 1 Tomar i = 1 Paso 2 Mientras que i ≤ No seguir Pasos 3 -6 Paso 3 Tomar p = p 0 – f(p 0) / f´(p 0) (Calcular pi) Paso 4 Si │p – p 0 │ < TOL entonces SALIDA (p); (Procedimiento completado satisfactoriamente) PARAR Paso 5 Tomar i = i + 1 Paso 6 Tomar p 0 = p (Redefinición de p 0, q 0, p 1, q 1) Paso 7 SALIDA (“El método fracaso después de No iteraciones, No =“, No) (Procedimiento completado sin éxito) PARAR

Método de la Secante 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. Por ejemplo: f(x)=x 23 xcos(2 x), entonces f´(x)= 2 x 3 xcos(2 x)+x 23 x(cos(2 x))ln 3 -2 x 23 xsen(2 x) La cual es extremadamente tediosa de evaluar

Método de la Secante El método se basa en obtener la ecuación de la recta que pasa por los puntos (xn− 1, f(xn− 1)) y (xn, f(xn)). A dicha recta se le llama secante por cortar la gráfica de la función. En la imagen anterior se toman los puntos iniciales x 0 y x 1, se construye una línea por los puntos (x 0, f(x 0)) y (x 1, f(x 1)). En forma punto-pendiente.

Método de la Secante Esta línea tiene la siguiente ecuación Posteriormente se escoge como siguiente elemento de la relación de recurrencia, xn+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 xn y xn-1).

Derivación del Método de la Secante Método de Newton Aproximación de la derivada Sustitución de la ecuación 2 en la ecuación 1 nos da el Método de la Secante (1) (2)

Algoritmo del Método de la Secante Para encontrar una solución de f(x)=0 dada la función continua f y unas aproximaciones iniciales p 0, p 1. ENTRADA aproximaciones iniciales p 0, p 1; tolerancia TOL; máximo numero de iteraciones No. SALIDA solución aproximada p o mensaje de fracaso. Paso 1 Tomar i = 2; q 0 = f(p 0); q 1 = f(p 1) Paso 2 Mientras que i ≤ No seguir Pasos 3 -6 Paso 3 Tomar p = p 1 – q 1(p 1 - p 0) / (q 1 - q 0) (Calcular pi) Paso 4 Si │p – p 1 │ < TOL entonces SALIDA (p); (Procedimiento completado satisfactoriamente) PARAR Paso 5 Tomar i = i + 1 Paso 6 Tomar p 0 = p 1 (Redefinición de p 0, q 0, p 1, q 1) q 0 = q 1 p 1 = p q 1 = f(p) Paso 7 SALIDA (“El método fracaso después de No iteraciones, No =“, No) (Procedimiento completado sin éxito) PARAR

Ejemplo

MÉTODO DE LA REGLA FALSA f(x) x

MÉTODO DE LA REGLA FALSA Consiste en considerar un intervalo (xi, xi+1) en el que se garantice que la función tiene raíz.

MÉTODO DE LA REGLA FALSA f(x) f(xi) xi f(xi+1) Xi+1 x

MÉTODO DE LA REGLA FALSA Consiste en considerar un intervalo (xi, xi+1) en el que se garantice que la función tiene raíz. Se traza una recta que une los puntos (xi, f(xi)), (xi+1, f(xi+1)).

MÉTODO DE LA REGLA FALSA f(x) f(xi) xi f(xi+1) Xi+1 x

MÉTODO DE LA REGLA FALSA Consiste en considerar un intervalo (xi, xi+1) en el que se garantice que la función tiene raíz. Se traza una recta que une los puntos (xi, f(xi)), (xi+1, f(xi+1)). Se obtiene el punto de intersección de esta recta con el eje de las abscisas: (xr, 0); se toma xr como aproximación de la raíz buscada.

MÉTODO DE LA REGLA FALSA f(x) f(xi) f(xr) f(xi+1) xi xr xi+1 x

MÉTODO DE LA REGLA FALSA Consiste en considerar un intervalo (xi, xi+1) en el que se garantice que la función tiene raíz. Se traza una recta que une los puntos (xi, f(xi)), (xi+1, f(xi+1)) y se obtiene el punto de intersección de esta recta con el eje de las abscisas: (xr, 0); se toma xr como aproximación de la raíz buscada. Se identifica luego en cuál de los dos intervalos está la raíz.

MÉTODO DE LA REGLA FALSA f(x) x i+1 = x r f(xi) f(xr) f(xi+1) xi xr xi+1 x

MÉTODO DE LA REGLA FALSA Consiste en considerar un intervalo (xi, xi+1) en el que se garantice que la función tiene raíz. Se traza una recta que une los puntos (xi, f(xi)), (xi+1, f(xi+1)) Se obtiene el punto de intersección de esta recta con el eje de las abscisas: (xr, 0); se toma xr como aproximación de la raíz buscada. Se identifica luego en cuál de los dos intervalos está la raíz. El proceso se repite n veces, hasta que el punto de intersección xr coincide prácticamente con el valor exacto de la raíz.

MÉTODO DE LA REGLA FALSA f(x) x r = x i+1 - f(xi) f(xr) f(xi+1) xi xr xi+1 f (x i+1)( x i+1 - x i ) f ( x i+1) - f ( x i ) x

ALGORITMO DE LA REGLA FALSA ENTRADA aproximaciones iniciales (p 0, p 1) = (a, b) ; tolerancia TOL; máximo numero de iteraciones No. SALIDA solución aproximada p o mensaje de fracaso. Paso 1 Tomar i = 2; q 0 = f(p 0); q 1 = f(p 1) Paso 2 Mientras que i ≤ No seguir Pasos 3 -6 Paso 3 Tomar p = p 1 – q 1(p 1 - p 0) / (q 1 - q 0) (Calcular pi) Paso 4 Si │p – p 1 │ < TOL entonces SALIDA (p); (Procedimiento completado satisfactoriamente) PARAR Paso 5 Tomar i = i + 1 , q = f(p) Paso 6 Si q * q 1 < 0 entonces Tomar p 0 = p (Redefinición de p 0, q 0, p 1, q 1) q 0 = q Si no p 1 = p q 1 = q Paso 7 SALIDA (“El método fracaso después de No iteraciones, No =“, No) (Procedimiento completado sin éxito) PARAR