qwertyuiopasdfghjklzxcvbnmqwe
rtyuiopasdfghjklzxcvbnmqwertyu
iopasdfghjklzxcvbnmqwertyuiopa
sdfghjklzxcvbnmqwertyuiopasdfg
INSTITUTO POLITECNICO NACIONAL
hjklzxcvbnmqwertyuiopasdfghjkl
ESCUELA SUPERIOR DE INGENIERIA MECANICA Y ELECTRICA
zxcvbnmqwertyuiopasdfghjklzxcv
Alumno: Jorge Chávez Torres
bnmqwertyuiopasdfghjklzxcvbnm
Materia: Análisis Numérico
qwertyuiopasdfghjklzxcvbnmqwe
rtyuiopasdfghjklzxcvbnmqwertyu
Profesor: María Alicia Ramírez Cruz
iopasdfghjklzxcvbnmqwertyuiopa
Practica N°2 : Bisección
sdfghjklzxcvbnmqwertyuiopasdfg
hjklzxcvbnmqwertyuiopasdfghjkl
zxcvbnmqwertyuiopasdfghjklzxcv
bnmqwertyuiopasdfghjklzxcvbnm
qwertyuiopasdfghjklzxcvbnmqwe
rtyuiopasdfghjklzxcvbnmrtyuiopa
Introducción
El método de la bisectriz es un método numérico sencillo, muy versátil para encontrar una
raíz real en un intervalo en el que existe una raíz.
Su singular ventaja consiste en que funciona incluso con funciones no analíticas; sin
embargo, sólo se debe utilizar el método después de un análisis gráfico.
Supongáse que hay una raíz de f(x) = 0 en un intervalo entre;
x = a y x = c,
denotado por [a, c] o, de forma equivalente, por c ³ x ³ a.
El método de la bisectriz se basa en el hecho de que, cuando un intervalo [a, c] tiene una
raíz, el signo de y(x) en los extremos es distinto, o sea, f(a) f(c) < 0.
El primer paso de este método consiste en bisectar el intervalo [a, c] en dos mitades: [a, b] y
[b, c], donde;
b=(a+b)/2
Si se verifican los signos de f(a)f(b) y f(b)f(c), se sabe en que mitad del intervalo se
encuentra la raíz.
De hecho, si 0 ³ f(a)f(b), el intervalo [a, b], que incluye x = a y x = b, contiene a la raíz; de
lo contrario, la raíz esta en el otro intervalo, [b, c].
A continuación, se bisecta de nuevo el intervalo que contiene a la raíz.
Al repetir este procedimiento, el tamaño del intervalo que contiene a la raíz se hará cada
vez más pequeño.
En cada paso se toma el punto medio del intervalo como la aproximación más actualizada a
la raíz.
La iteración se detiene cuando el tamaño de la mitad del intervalo es menor que una
tolerancia dada.
El tamaño del intervalo después de n pasos de iteración es
(c0- a0)/2n
donde a0 y c0 son valores iniciales, de modo que el numerador es el tamaño de intervalo
inicial.
La ecuación anterior representa el máximo error que existe cuando se aproxima la raíz con
el n-ésimo punto medio. Por tanto, si la tolerancia del error es t , el número de pasos de
iteración necesarios es el entero n más pequeño que satisface
t ³( c0- a0)/2n
De forma equivalente, n ³ log[(c0- a0)/ t ]/log(2) donde t es la tolerancia.
Desarrollo
Algoritmo:
1. Elija valores Iniciales para “a” y “b” de forma tal que lea función cambie de signo sobre
el intervalo. Esto se puede verificar asegurándose de que:
f(a)*f(b) < 0
2. La primera aproximación a la raíz se determina con la fórmula:
xn = (a + b) / 2
3. Realizar las siguientes evaluaciones para determinar en que subintervalo se encuentra la
raíz:
f(a)*f(xn ) < 0 Entonces b = xn
f(a)*f(xn) > 0 Entonces a = xn
f(a)*f(xn) = 0 Entonces xn Es la Raíz
4. Calcule la nueva aproximación
xn+1 = (a + b) / 2
5. Evaluar la aproximación relativa
| (xn+1 - xn ) / xn+1 | < Tolerancia
No. (Falso) Repetir el paso 3, 4 y 5
Sí. (Verdadero) Entonces xn+1 Es la Raíz
Codigo:
#include<stdio.h>
#include<conio.h>
#include<math.h>
#include<stdlib.h>
double funcion(float X); //declaracion de la funcion.
int main(void){
int SI=1;
do{
system("cls" );
float Xa=0, Xb=0,error=0,Xn=0,Xr=0,promedio=0,p=0;
int cont=0;
int I;
printf("PRACTICA #2" );
printf(" METODO DE BISECCION" );
printf("\nintroduce Xa: " );
scanf("%f",&Xa);
printf("\nintroduce Xb: " );
scanf("%f",&Xb);
printf("Numero de iteraciones: ");
scanf("%d",&I);
printf("\nXA XB Xr f(Xa) f(Xb) f(Xr) error\n" );
Xr=(Xb+Xa)/2;
printf("\n%.3f %.3f %.3f %.3f %.3f %.3f ",Xa,Xb,Xr,funcion(Xa),funcion(Xb),funcion(Xr));
do{
promedio=funcion(Xa)*funcion(Xr);
if(promedio==0)
printf("\nLa raiz es %f",Xb);
else
if(promedio<0)
Xb= Xr;
else
Xa= Xr;
Xn=(Xb+Xa)/2;
error=fabs((Xn-Xr)/Xn)*100;
printf("\n%.3f %.3f %.3f %.3f %.3f %.3f %.3f
%",Xa,Xb,Xn,funcion(Xa),funcion(Xb),funcion(Xr),error);
Xr=Xn;
cont++;
}while(cont<I);
printf("\n");
printf("\nDESEAS REALIZAR LA OPERACION OTRA VEZ? SI[1] NO[2] " );
scanf("%d",&SI);
}while(SI<=1);
system("PAUSE" );
return 0;
}
double funcion(float X){ //funcion que desp es llamada dentro del main.
return exp(X)+ log(X);
}
Conclusion:
El método lo veo muy sencillo en la forma de aplicar pero me cuesta ubicar las partes a operar
después de la 3° iteración ya que me revuelven tantos números y eso lo resuelvo con el programa
ya que con dar el número de iteraciones este las despliega en una tabla y es más fácil ubicar un
valor o valores específicos.
Bibliografia:
Título Métodos numéricos para ingenieros
Autores Steven C. Chapra, Raymond P. Canale
Edición 5
Editor McGraw-Hill, 2007
ISBN 9701061144, 9789701061145
N.º de páginas 977 páginas
