Biseccion

4.metodo de la biseccion

Document Transcript

223 RAÍCES REALES DE ECUACIONES NO-LINEALESSea f: R→R. Dada la ecuación f(x) = 0, se debe encontrar un valor real r tal que f(r) = 0.Entonces r es una raíz real de la ecuaciónSi no es posible obtener la raíz directamente, entonces se debe recurrir a los métodos numéricositerativos para calcular r en forma aproximada con alguna precisión controlada. Se han creadomuchos métodos numéricos para resolver este problema clásico, pero con el uso decomputadoras para el cálculo, conviene revisar solamente algunos de estos métodos que tengancaracterísticas significativamente diferentes.3.1 Método de la bisecciónSea f: R→R. Suponer que f es continua en [a, b], y que además f(a) y f(b) tienen signosdiferentes. Por continuidad, el intervalo (a, b) contendrá al menos una raíz real.El siguiente teorema establece la existencia de la raíz r:Teorema de Bolzano: Si una función f es continua en un intervalo [a, b] y f(a) tiene signodiferente que f(b), entonces existe por lo menos un punto r en (a, b) tal que f(r)=0.Si además f(x) no cambia de signo en el intervalo [a, b], entonces la solución es única.El método de la bisección es un método simple y convergente para calcular r. Consiste encalcular el punto medio c=(a+b)/2 del intervalo [a, b] y sustituirlo por el intervalo [c, b] ó [a, c]dependiendo de cual contiene a la raíz r. Este procedimiento se repite hasta que la distanciaentre a y b sea muy pequeña, entonces el último valor calculado c estará muy cerca de r.Interpretación gráfica del método de la bisecciónEn la figura se puede observar que luego de haber calculado c, para la siguiente iteración debesustituirse el intervalo [a, b] por [c, b] debido a que f(a) y f(c) tienen igual signo y por lo tanto laraíz estará en el intervalo [c, b]3.1.1 Convergencia del método de la bisecciónSean ai, bi, ci los valores de a, b, c en cada iteración i=1, 2, 3, . . . respectivamenteEl método de la bisección genera una sucesión de intervalos [a, b], [a1, b1], [a2, b2], …, [ai, bi]tales que a ≤ a1 ≤ a2 … ≤ ai constituyen una sucesión creciente y b ≥ b1 ≥ b2 …, ≥ bi unasucesión decreciente con ai < bi. Además por definición del método: ci, r ∈ [ai, bi] en cadaiteración i

23Sean di = bi – ai longitud del intervalo [ai, bi] en la iteración i=1, 2, 3, . . . d = b – a longitud del intervalo inicialRecorrido de las iteraciones Iteración Longitud del intervalo 1 d1 = d /2 2 2 d2 = d1/2 = d/2 3 3 d3 = d2/2 = d/2 4 4 d4 = d3/2 = d/2 ... ... i i di = d/2Entonces d → 0 ⇒ di → 0 ⇒ ai → bi ⇒ ci → r ⇒ ∃i>0 | ci − r |< ε para cualquier valor positivo ε 2i i→∞ i→∞ i→∞i→∞Suponer que se desea que el último valor calculado ci tenga precisión E = 0.001, entonces si elalgoritmo termina cuando bi – ai < E, se cumplirá que |ci – r| < E y ci será una aproximaciónpara r con un error menor que 0.0001Ejemplo. Calcule una raíz real de f(x) = x e - π = 0 en el intervalo [0, 2] con precisión 0.01 xLa función f es continua y además f(0)<0, f(2)>0, por lo tanto la ecuación f(x)=0 debe conteneralguna raíz real en el intervalo [0, 2]Cálculo manual para obtener la raíz con el método de la Bisección iteración a b c sign(f(a)) sign(f(c)) inicio 0 2 1 - - 1 1 2 1.5 - + 2 1 1.5 1.25 - + 3 1 1.25 1.125 - + 4 1 1.125 1.0625 - - 5 1.0625 1.125 1.0938 - + 6 1.0625 1.0938 1.0781 - + 7 1.0625 1.0781 1.0703 - - 8 1.0703 1.0781 1.0742En la última iteración se observa que el intervalo que contiene a la raíz se ha reducido a[1.0703, 1.0781], por lo tanto el último valor calculado de c = 1.0742 debe estar cerca de r conuna distancia no mayor a 0.01

243.1.2 Eficiencia del método de la bisecciónSuponer el caso más desfavorable, en el que r está muy cerca de uno de los extremos delintervalo [a, b]:Sean Ei = r − ci : error en la iteración i Ei+ 1= r − ci+ 1 : error en la iteración

...