METDO DE LA BISECCIN
Jhon Emanuel Morocho Guaman
[email protected]@est.ups.edu.ec
@est.ups.edu.ec
@est.ups.edu.ec
UNIVERSIDAD POLITECNICA SALESIANA
Resumen
Abstract
En el presente informe se encuentra detallado
el mtodo de la biseccin que es un algoritmo de
bsqueda de races que trabaja dividiendo
el intervalo a la mitad y seleccionando el
subintervalo que tiene la raz, que es uno de los
mtodos ms sencillos y de fcil intuicin para
resolver ecuaciones en una variable. Se basa en el
teorema del valor intermedio (TVI).
Para
ello
presentamos
un
anlisis
y
comprobaciones de las races extradas mediante el
mtodo de la biseccin, donde demostraremos que
los anlisis obtenidos matemticamente tambin
podemos obtener mediante un programa
(MATLAB) de una manera mucho ms rpida la
raz, para eso presentaremos con su respectivos
algoritmos y cdigos
The present report the bisection method which is a
root finding algorithm works by dividing the interval
in half and selecting the subinterval having the root,
which is one of the most simple and easy methods
intuition is detailed solving equations in one variable.
It is based on the intermediate value theorem (TVI).
Therefore we present an analysis of the roots and
checks drawn by the bisection method, which
demonstrate that mathematically analyzes obtained
can also be obtained through a program (MATLAB)
in a much more effective and fast way to present it
with their respective algorithms and codes
Palabras Claves:
Keywords:
Biseccin, Algoritmo, Variables.
Bisection Algorithm, Variables.
Objetivos
Desarrollar y analizar de una funcin sus
races aplicando el mtodo de la biseccin
El usuario debe ser capaz de desarrollar un
programa en matlab el cual encuentre de
manera correcta las races de diferentes
polinomios utilizando los mtodos de
biseccin
Comparar los resultados obtenidos de las
dos formas realizadas analticamente y
programado en matlab.
El objetivo de esta prctica es la
implementacin del mtodo de biseccin,
y la comparacin de ambos mtodos para
la resolucin de algunas ecuaciones no
lineales
figura1
Esta es una consecuencia del teorema del valor
intermedio para funciones continuas, que
establece que si f es continua en [a,b] y si k es un
nmero entre f(a) y f(b) , entonces existe por lo
menos
un c (a,b)
tal
que f(c)=k.
(para el caso en que f(a)f(b)<0 se escoge k=0,
luego f(c)=0, c (a,b)).
1. Introduccin
El mtodo de biseccin, conocido tambin como
de corte binario, de participacin en dos intervalos
iguales o mtodo Bolzano, es un mtodo de
bsqueda incremental en el que el intervalo se
divide siempre en dos. Si la funcin cambia su
signo sobre un intervalo, se evala el valor de la
funcin en el punto medio. La posicin de la raz
se determina situndola en el punto medio del
subintervalo dentro del cual ocurre un cambio. El
proceso
se
repite
hasta
obtener
una
mejor aproximacin.
2. Marco terico
2.1.
Mtodo de la Biseccin
El mtodo de la biseccin o corte binario es un
mtodo de bsqueda incremental que divide el
intervalo siempre en 2. Si la funcin cambia de
signo sobre un intervalo, se evala el valor de la
funcin en el punto medio. La posicin de la raz
se determina situndola en el punto medio del
subintervalo donde exista cambio de signo.
El proceso se
repite
hasta
mejorar
la
aproximacin.
El mtodo de biseccin consiste en dividir el
intervalo en 2 subintervalos de igual magnitud,
reteniendo el subintervalo en donde f cambia de
signo, para conservar al menos una raz o cero, y
repetir el proceso varias veces.
2.2.
En que se basa el mtodo de la
biseccin
El Mtodo de Biseccin se basa en la bsqueda
incremental donde el intervalo se divide siempre
en dos.
Si la funcin previamente planteada cambia de
signo sobre un intervalo, se evala el valor de la
funcin en el punto medio de este intervalo.
La posicin de la raz se determina situndola en
el punto medio del subintervalo dentro del cual
ocurre un cambio de signo. El proceso se repite
hasta obtener una mejor aproximacin.
2.3.
El mtodo consiste en lo siguiente:
Debe existir seguridad sobre la
continuidad de la funcin f(x) en el
intervalo [a, b]
A continuacin se verifica
que
Se calcula el punto medio m del
intervalo [a,b] y se evala f(m) si ese
valor es igual a cero, ya hemos
encontrado la raz buscada
En caso de que no lo sea, verificamos
si f(m) tiene signo opuesto con f(a) o
con f(b)
Se redefine el intervalo [a, b] como
Si
y
).
pasa a ser el punto medio (
a ser el punto medio (
).
Paso 4
Si f (Xa) f (Xb)=0 Error = | Xpm Xpm 1 |
< Tolerancia, donde Xpm es el punto medio de la
iteracin actual y Xpm 1 es el punto medio de
la iteracin anterior.
Al cumplirse la condicin del Paso 4, la raz o
solucin es el ltimo punto medio que se obtuvo.
Para el error relativo porcentual se tiene la
siguiente frmula:
sucesivamente encerrando la solucin
en un intervalo cada vez ms
pequeo, hasta alcanzar la precisin
deseada
Algoritmo
la
Si
, entonces la
solucin o raz est fuera del intervalo
entre Xa y el punto medio, y
pasa
determinado en cul de estos
intervalos ocurre un cambio de signo
Con este nuevo intervalo se contina
3.
Paso 1
Elegir los valores iniciales Xa y Xb, de tal forma
de que la funcin cambie de signo:
DESARROLLO DEL EJERCICIO
Donde tenemos la funcin
en el intervalo de
1) Obtenemos la grfica y la tabla de
datos
Paso 2
La primera aproximacin a la raz se determina
con la frmula del punto medio de esta forma:
Paso 3
Realizar las siguientes evaluaciones
determinar el intervalo de la raz:
entonces
solucin o raz est entre
[a, m] [m, b] segn se haya
x
-1
-0,5
0
0,5
1
para
F(x)
-4
-2,31
-2
-1,562
2
Figura2
1) Desarrollo del mtodo de la biseccin
Debemos considerar algunas condiciones en la
formulacin
4.
DESARROLLO DE CDIGOS EN
MATLAB
Procedemos a codificar en matlab
1)
2)
3)
4)
clear all
clc
% Metodo de biseccion %
fprintf ('/// Metodo de biseccion
/// \n ')
syms x;
5) f=input ('Ingrese la funcion f=')
6) ezplot (f)
7) grid on;
8)
9) xl=input('ingrese extremo
izquierdo del intervalo xl=')
10) xu=input(' ingrese extremo
derecho del intervalo xu=')
11) es=input('ingrese el error
deseado ea= ')
12)
13) x=xl;
14) fxl=eval(f)
15) x=xu;
16) fxu=eval(f)
17) if fxl*fxu>0;
18)
disp ('no hay cambio de signo
en los extremos')
19)
break
20) end
21)
22) ea=1000;
23)
24) i=0;
25) fprintf('interaciones raiz
error \n')
26)
27) while ea>es;
xr=((xl+xu)/2);
28)
x=xr;
29)
fxr=eval (f);
30)
31)
if fxr*fxl<0;
32)
xu=xr;
33)
elseif fxr*fxl>0;
34)
xl=xr;
else
35)
disp('la raiz es x= ')
Tambin
1
2
3
4
5
6
0,5
0,75
0,75
0,75
0,782
0,798
-0,418
0,596
0,049
-0,191
-0,070
-0,007
1
1
0,875
0,813
0,813
0,813
0,75
0,875
0,813
0,782
0,798
0,806
14,28
7,62
3,96
2
0.99
(-1,562)*( -0,418)=0,65350
(-0,418)*( 0,596)=-0,24920
(-0,418)*( 0,049)=-0,0202
(-0,418)*( -0,191)=0,080
(-0,191)*( -0,070)=0,0013
(-0,070)*( -0,007)=0,0005
36)
37)
38)
39)
40)
41)
42)
43)
xr
break
end
if i>1
ea=abs(((xrxrant)/xr)*100);
end
xrant=xr;
i=i+1;
fprintf('%6.0f %6.5f
%6.5f\n',i,xr,ea)
end
44)
45)
46) if ea<=es
47)
disp('la raiz es x=')
48)
xr
49)
disp('error es ea=')
50)
i
51) end
52) %FINAL%
Figura 3
5.3 Y a continuacin nos da
correspondiente.
Interacciones
raz
error
1
0.00000
1000.00000
2
0.50000
1000.00000
3
0.75000
33.33333
4
0.87500
14.28571
5
0.81250
7.69231
6
0.78125
4.00000
7
0.79688
1.96078
8
0.80469
0.97087
9
0.80859
0.48309
10
0.80664
0.24213
11
0.80762
0.12092
12
0.80713
0.06050
13
0.80688
0.03026
14
0.80701
0.01513
15
0.80695
0.00756
16
0.80692
0.00378
17
0.80690
0.00189
18
0.80689
0.00095
La raz es x= xr = 0.8069
Procedimiento
Luego de ingresar el algoritmo del mtodo de la
biseccin
correspondiente
en
matlab,
procedemos a mandar a correr y el programa nos
pedir que ingresemos la funcin y su intervalo
para calcular la raz, a continuacin tenemos
dicho procedimiento.
/// Metodo de biseccion ///
Ingrese la funcion f= x^4+3*x^3-2
f = x^4 + 3*x^3 - 2
ingrese extremo izquierdo del intervalo xl= 1
xl = -1
ingrese extremo derecho del intervalo xu= 1
xu = 1
ingrese el error deseado ea= 0.001
Luego de ingresar procedemos a dar un enter y
este nos dar primero la grfica de la funcin.
la
raz
5. Conclusiones
El mtodo de la Biseccin converge lentamente,
lo que genera la propagacin de error por la
cantidad de operaciones e iteraciones necesaria
para que el mtodo converja.
6. Bibliografa
[1] Libro de metodo numericos CHAPRA 5ta
edicion.
[2] W. Rafferty, Ground antennas in NASAs
deep space telecommunications, Proc. IEEE
vol. 82, pp. 636-640, 1994.
[3] http://www.monografias.com/trabajos43/
metodo-biseccion/metodobiseccion.shtml#ixzz3Z1WplbDx
[4] http://portales.puj.edu.co/objetosdeaprendizaj
e/Online/OA10/capitulo5/5.htm
Figura b
7. Anexos
Figura c
Figura a
Figura d
Figura e
