RACES DE ECUACIONES

2010

RACES DE ECUACIONES

CONTENIDO

Definicin
Mtodos para la aproximacin de soluciones
1. Mtodo grafico
2. Cerrado o acotado:
a) Biseccin
b) Falsa Posicin
3. Abierto:
c) Secante
d) Newton-Raphson
e) Punto Fijo


RACES DE ECUACIONES

DEFINICIN

El objeto del clculo de las races de una ecuacin es determinar
los valores de x para los que se cumple:

f(x) = 0
Su importancia radica en que si podemos determinar las races de
una ecuacin tambin podemos determinar mximos y mnimos,
valores propios de matrices, resolver sistemas de ecuaciones
lineales y diferenciales, etc...

RACES DE ECUACIONES

Para resolver ecuaciones no lineales existen varios mtodos
numricos que los podemos clasificar as:

Cerrado o acotado:
(requiere de dos valores de
x que encierren la raz)

Biseccin

Falsa posicin
Abierto: (requiere de
uno o dos valores de x,
pero no necesariamente
encierran la raz)
Punto fijo
Newton-Raphson
Secante
Mtodo grafico
RACES DE ECUACIONES





La mayora de los mtodos utilizados para el clculo de las races de
una ecuacin son iterativos y se basan en modelos de aproximaciones
sucesivas. Estos mtodos trabajan del siguiente modo: a partir de una
primera aproximacin al valor de la raz, determinamos una
aproximacin mejor aplicando una determinada regla de clculo y as
sucesivamente hasta que se determine el valor de la raz con el grado
de aproximacin deseado.


MTODO GRAFICO

Consiste en graficar una funcin y determinar visualmente donde corta
el eje x. En y= f(x), establece el valor de x para el cual f(x)=0.

f(x)
a b
x
1. Si en un intervalo {a,b} cerrado se cumple que : no
existen races reales en el intervalo, pues y=f(x) no toca el eje x, por el
contrario pueden encontrarse una o ms races imaginarias.
f(a).f(b)>0
MTODO GRAFICO

2. Si en un intervalo {a,b} cerrado se cumple que :
Entonces existen dos races reales
f(a).f(b)>0
f(x)
a b
x
MTODO GRAFICO


3. Si en un intervalo {a,b} cerrado se cumple que :
da la certeza de encontrar una sola raz real en el intervalo.
f(a).f(b)<0
x
f(x)
a b

MTODO GRAFICO


4. Si en un intervalo {a,b} cerrado se cumple que :
hay ms de dos races.
f(a).f(b)<0
f(x)
a b
x

MTODO GRAFICO

5. Tambin puede existir una funcin , para la que existe una raz
real doble en x=0 , que no es apreciable por el mtodo grfico, pues
la ecuacin es tangente al eje x .
f(x)
a b
x

MTODO DE BISECCIN


Este mtodo, tambin conocido como mtodo de particin del
intervalo, parte de una ecuacin algebraica o trascendental f(x) y
un intervalo [xi, xs], tal que f(xi) y f(xs) tienen signos contrarios, es
decir, tal que existe por lo menos una raz en ese intervalo.
MTODO DE BISECCIN

El mtodo consiste en considerar un intervalo (xi,xs) en el que
se garantice que la funcin tiene raz.
El segmento se bisecta, tomando el punto de biseccin xr como
aproximacin de la raz.


se identifica luego en cual de los dos intervalos esta la raz.
Si: f(xi).f(xr)<0 la raiz esta en el intervalo {xi,xr} y xs=xr
Si: f(xr).f(xs)<0 la raiz esta en el intervalo {xr,xs} y xi=xr

el proceso se repite n veces, hasta que el punto de biseccin xr
coincide prcticamente con el valor exacto de la raz.

( )
2
s i
r
x x
x
+
=
MTODO DE BISECCIN
f(x)
f(xi)
f(xr)
f(xs)
xi xr xs
xi=xr
x
MTODO DE BISECCIN
El mtodo de biseccion se puede frenar:
con el numero mximo de iteraciones
cuando se alcanza el % de error

0
2
x A
1
2
x A
2
2
x A
3
2
x A
n
x
2
A
x A
xi xs
i s
x x x = A
n
i s
s
x x
2

> c
) 2 (
) ( ) (
In
In x In
n
s
c A
=
Para estimar el numero mximo
de iteraciones tenemos
Donde: x = longitud del intervalo
n= numero de iteraciones
error


n
i s
s
x x
2

> c
MTODO DE BISECCIN

Cuando se alcanza el %E

1. Alcanza la tolerancia exigida.



2. Alcanza el error relativo porcentual verdadero.



% 100
actual
r
anterior
r
actual
r
a
x
x x
E

=
% 100
verdadero
r
actual
r
verdadero
r
t
x
x x
E

=
s a
E c s
t a
E E <

MTODO DE BISECCIN

Ventajas:

Siempre converge.

til como aproximacin inicial de otros mtodos.

Desventajas:

No tiene en cuenta la magnitud de los valores de la funcin
en las aproximaciones calculadas xn, solo tiene en cuenta el
signo de f(x), lo que hace que una aproximacin intermedia,
mejor que la respuesta final, pase desapercibida.

Convergencia lenta.
iter xi xr xs f(xi) f(xr) f(xs) f(xi).f(xr) f(xs).f(xr) error
1 0,5 0,75 1 0,10653066 -0,27763345 -0,63212056 -0,02957647 0,17549781
2 0,5 0,625 0,75 0,10653066 -0,08973857 -0,27763345 -0,00955991 0,02491443 20
3 0,5 0,5625 0,625 0,10653066 0,00728282 -0,08973857 0,00077584 -0,00065355 11,1111111
4 0,5625 0,59375 0,625 0,00728282 -0,04149755 -0,08973857 -0,00030222 0,00372393 5,26315789
5 0,5625 0,578125 0,59375 0,00728282 -0,01717584 -0,04149755 -0,00012509 0,00071276 2,7027027
6 0,5625 0,5703125 0,578125 0,00728282 -0,00496376 -0,01717584 -3,615E-05 8,5257E-05 1,36986301
7 0,5625 0,56640625 0,5703125 0,00728282 0,0011552 -0,00496376 8,4131E-06 -5,7341E-06 0,68965517
8 0,56640625 0,56835938 0,5703125 0,0011552 -0,00190536 -0,00496376 -2,2011E-06 9,4577E-06 0,34364261
9 0,56640625 0,56738281 0,56835938 0,0011552 -0,00037535 -0,00190536 -4,336E-07 7,1518E-07 0,17211704
10 0,56640625 0,56689453 0,56738281 0,0011552 0,00038986 -0,00037535 4,5037E-07 -1,4633E-07 0,08613264
11 0,56689453 0,56713867 0,56738281 0,00038986 7,2379E-06 -0,00037535 2,8218E-09 -2,7167E-09 0,04304778

MTODO DE BISECCIN

Ejemplo:
Calcule la raz de:
FALSA POSICIN

Este mtodo, como en el mtodo de la biseccin, parte de dos puntos que
rodean a la raz f(x) = 0, es decir, dos puntos xi y xs tales que f(xi)f(xs) < 0.
La siguiente aproximacin, xr, se calcula como la interseccin con el eje X
de la recta que une ambos puntos empleando la ecuacin




La asignacin del nuevo intervalo de bsqueda se realiza como en el
mtodo de la biseccin: entre ambos intervalos, [xi,xr] y [xr,xs], se toma
aquel que cumpla f(xi)f(xr) < 0 ; f(xr)f(xs) < 0.
( )
) ( ) (
) (
s i
s i s
s r
x f x f
x x x f
x x

=
FALSA POSICIN

xi xr xs
f(x)
Raz falsa
Raz verdadera
FALSA POSICIN

iter xi xr xs f(xi) f(xr) f(xs) f(xi).f(xr) f(xs).f(xr) error
1 0,5 0,57211161 1 0,10653066 -0,00777908 -0,63212056 -0,00082871 0,00491732
2 0,5 0,56720422 0,57211161 0,10653066 -9,5491E-05 -0,00777908 -1,0173E-05 7,4283E-07 0,86518885
3 0,5 0,56714404 0,56720422 0,10653066 -1,172E-06 -9,5491E-05 -1,2485E-07 1,1192E-10 0,01061207
4 0,5 0,5671433 0,56714404 0,10653066 -1,4385E-08 -1,172E-06 -1,5324E-09 1,6859E-14 0,00013025
Ejemplo:
Calcule la raz de:
MTODO DE PUNTO FIJO

Usando el concepto de replantear la forma original del
problema:
Si

) ( ) ( x g x x f + =
x
Tal que
0 ) ( = x f
x
Tal que
0 ) ( = + x g x
) (x g x =
) (
1 i i
x g x =
+
MTODO DE PUNTO FIJO


Se pueden presentar cuatro situaciones al momento de buscar la raz.
1. Que y
) ( ) (
'
1
'
2
x f x f s
0 ) (
'
2
> x f
f2(x)
f1(x)
f(x)
x
RAIZ
xi
solucin monotnicamente convergente (mayor
acercamiento a la raz)
MTODO DE PUNTO FIJO


2. Que y
) ( ) (
'
1
'
2
x f x f s
0 ) (
'
2
< x f
f2(x)
f1(x)
f(x)
x
RAIZ
xi
solucin oscilatoriamente convergente (mayor acercamiento de manera
oscilatoria a la raz)
MTODO DE PUNTO FIJO

3. Que y
) ( ) (
'
1
'
2
x f x f >
0 ) (
'
2
> x f
f2(x)
f1(x)
f(x)
x
RAIZ
xi
solucin monotnicamente divergente (mayor alejamiento de
la raz)
MTODO DE PUNTO FIJO

4. Que y
) ( ) (
'
1
'
2
x f x f >
0 ) (
'
2
< x f
f2(x)
f1(x)
f(x)
x
RAIZ
xi
solucin oscilatoriamente divergente (mayor alejamiento de manera
oscilatoria de la raz)
MTODO DE PUNTO FIJO

iter xi x(i+1) f(xi) f(x(i+1) f(xi).f(x(i+1) error
1 0,5 0,60653066 0,10653066 -0,06129145 -0,00652942
2 0,60653066 0,54523921 -0,06129145 0,03446388 -0,00211234 17,5639365
3 0,54523921 0,57970309 0,03446388 -0,01963847 -0,00067682 11,2412032
4 0,57970309 0,56006463 -0,01963847 0,01110752 -0,00021813 5,94509212
5 0,56006463 0,57117215 0,01110752 -0,0063092 -7,008E-05 3,50646443
6 0,57117215 0,56486295 -0,0063092 0,0035751 -2,2556E-05 1,94468884
7 0,56486295 0,56843805 0,0035751 -0,00202859 -7,2524E-06 1,11694386
8 0,56843805 0,56640945 -0,00202859 0,00115018 -2,3333E-06 0,62893408
9 0,56640945 0,56755963 0,00115018 -0,00065242 -7,504E-07 0,35814989
10 0,56755963 0,56690721 -0,00065242 0,00036998 -2,4138E-07 0,20265386
11 0,56690721 0,5672772 0,00036998 -0,00020984 -7,7639E-08 0,11508432
12 0,5672772 0,56706735 -0,00020984 0,00011901 -2,4973E-08 0,06522085
13 0,56706735 0,56718636 0,00011901 -6,7496E-05 -8,0326E-09 0,03700515
14 0,56718636 0,56711886 -6,7496E-05 3,8279E-05 -2,5837E-09 0,02098221
Ejemplo:
Calcule la raz de:

MTODO DE NEWTON RAPHSON

Es uno de los mtodos mas usados en la ingeniera, por llegar al
resultado del problema de forma mas rpida.
Se basa en trazar rectas tangentes que toman la forma de la
funcin por medio de su primera derivada.



Se usa la proyeccin de la recta tangente para encontrar el valor
aproximado de la raz.


) (
) (
'
1
x f
x f
x x
i
i i
=
+
1
'
0 ) (
) (
+

=
i i
i
x x
x f
x f
MTODO DE NEWTON RAPHSON


f(x)
f(xi)
f(xi+1)
Xi+1
xi
RAIZ

MTODO DE NEWTON RAPHSON


Desventajas:

Lenta convergencia debida a la naturaleza de una funcin en
particular.

Cuando un punto de inflexin, f(x) = 0, ocurre en la vecindad de
una raz.
No existe un criterio general de convergencia.
Tener un valor suficientemente cercano a la raz.
Apoyarse de herramientas grficas.
Conocimiento del problema fsico.

Evaluacin de la derivada.

MTODO DE NEWTON RAPHSON


iter xi f(xi) f``(xi) x(i+1) f(xi+1) error
1 0,5 0,10653066 -1,60653066 0,566311 0,00130451
2 0,566311 0,00130451 -1,56761551 0,56714317 1,9648E-07 11,709291
3 0,56714317 1,9648E-07 -1,56714336 0,56714329 4,4409E-15 0,14672871
Ejemplo:
Calcule la raz de:
MTODO DE SECANTE


Este mtodo, a diferencia del de biseccin y regla falsa, casi nunca
falla ya que solo requiere de 2 puntos al principio, y despus el
mismo mtodo se va retroalimentando.

Lo que hace bsicamente es ir tirando rectas secantes a la curva
de la ecuacin que se tiene originalmente, y va checando la
interseccin de esas rectas con el eje de las X para ver si es la raz
que se busca.



MTODO DE SECANTE


El mtodo se define por:




Como se puede ver, este mtodo necesitar dos aproximaciones iniciales
de la raz para poder inducir una pendiente inicial.

) ( ) (
) )( (
1
1
1

=
i i
i i i
i i
x f x f
x x x f
x x
Despejando
MTODO DE SECANTE
Xi-1
xi
Xi+1
f(xi-1)
f(xi)
A
E
B
D
C
x
AD
AC
DE
BC
=
MTODO DE SECANTE
iter x(i-1) xi x(i+1) f(xi-1) f(xi) f(xi+1) f(xi-1).f(xi) f(xi+1).f(xi) error
1 0,5 1 0,57211161 0,10653066 -0,63212056 -0,00777908 -0,06734022 0,00491732
2 0,5 0,57211161 0,56720422 0,10653066 -0,00777908 -9,5491E-05 -0,00082871 7,4283E-07 74,7910687
3 0,5 0,56720422 0,57211161 0,10653066 -9,5491E-05 -0,00777908 -1,0173E-05 7,4283E-07 0,86518885
Ejemplo:
Calcule la raz de:
BIBLIOGRAFA

1. http://www.uv.es/diaz/mn/node17.html
2. http://www.scribd.com/doc/15638680/Metodo-de-NewtonRaphson
3. http://www.virtualum.edu.co/antiguo/metnum/raices/metgraf.htm
4. http://aprendeenlinea.udea.edu.co/lms/moodle/course/view.php?id=
229
5. http://www.mitecnologico.com/Main/MetodoDeLaSecante

6. SANTAFE, Elkin R. Elementos bsicos de modelamiento
matemtico. Clases-Universidad Industrial de Santander Ao-
2009.