Raı́ces de ecuaciones no lineales

Curso: Métodos Numéricos en Ingenierı́a

Profesor: Dr. José A. Otero Hernández
Correo: [email protected]
web: http://metodosnumericoscem.weebly.com
Universidad: ITESM CEM
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

TÓPICOS

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Tópicos

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué son los métodos abiertos?

Los métodos abiertos para la determinación de raı́z
requieren únicamente de un solo valor de inicio (xi ) o que
comiencen con dos valores; pero que no encierren la raı́z.
Los métodos abiertos emplean una fórmula para calcular
la raı́z.
Los métodos abiertos a veces divergen o se alejan de la
raı́z verdadera a medida que avanza el cálculo.
Cuando los métodos abiertos convergen lo hacen mucho
más rápido que los métodos cerrados.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué son los métodos abiertos?

Los métodos abiertos para la determinación de raı́z
requieren únicamente de un solo valor de inicio (xi ) o que
comiencen con dos valores; pero que no encierren la raı́z.
Los métodos abiertos emplean una fórmula para calcular
la raı́z.
Los métodos abiertos a veces divergen o se alejan de la
raı́z verdadera a medida que avanza el cálculo.
Cuando los métodos abiertos convergen lo hacen mucho
más rápido que los métodos cerrados.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué son los métodos abiertos?

Los métodos abiertos para la determinación de raı́z
requieren únicamente de un solo valor de inicio (xi ) o que
comiencen con dos valores; pero que no encierren la raı́z.
Los métodos abiertos emplean una fórmula para calcular
la raı́z.
Los métodos abiertos a veces divergen o se alejan de la
raı́z verdadera a medida que avanza el cálculo.
Cuando los métodos abiertos convergen lo hacen mucho
más rápido que los métodos cerrados.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué son los métodos abiertos?

Los métodos abiertos para la determinación de raı́z
requieren únicamente de un solo valor de inicio (xi ) o que
comiencen con dos valores; pero que no encierren la raı́z.
Los métodos abiertos emplean una fórmula para calcular
la raı́z.
Los métodos abiertos a veces divergen o se alejan de la
raı́z verdadera a medida que avanza el cálculo.
Cuando los métodos abiertos convergen lo hacen mucho
más rápido que los métodos cerrados.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué son los métodos abiertos?

Los métodos abiertos para la determinación de raı́z
requieren únicamente de un solo valor de inicio (xi ) o que
comiencen con dos valores; pero que no encierren la raı́z.
Los métodos abiertos emplean una fórmula para calcular
la raı́z.
Los métodos abiertos a veces divergen o se alejan de la
raı́z verdadera a medida que avanza el cálculo.
Cuando los métodos abiertos convergen lo hacen mucho
más rápido que los métodos cerrados.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué métodos abiertos estudiaremos?

Métodos de iteración simple de punto fijo,
Método de Newton-Raphson,
Método de la secante,
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué métodos abiertos estudiaremos?

Métodos de iteración simple de punto fijo,
Método de Newton-Raphson,
Método de la secante,
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué métodos abiertos estudiaremos?

Métodos de iteración simple de punto fijo,
Método de Newton-Raphson,
Método de la secante,
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Métodos abiertos

¿Qué métodos abiertos estudiaremos?

Métodos de iteración simple de punto fijo,
Método de Newton-Raphson,
Método de la secante,
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Tópicos

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Métodos de iteración simple de punto fijo

Toda ecuación f (x) = 0 puede transformarse de tal modo que
la variable x esté del lado izquierdo de la ecuación.

f (x) = 0 → x = g(x)

De esta forma, dado un valor inicial para la raı́z xi , la ecuación

x = g(x) se utiliza para obtener una nueva aproximación de la
raı́z xi+1 , esto es:
xi+1 = g(xi )
Se llega a una fórmula iterativa para obtener la raı́z. El error
aproximado se puede obtener:
 
 xi+1 − xi 
εa =   100%
xi+1 
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Ejemplos
x2 +3
x2 − 2x + 3 = 0 → x = 2 , por tanto:
x2i +3
xi+1 = 2

sin(x) = 0 → x = sin(x) + x, por tanto:

xi+1 = sin(xi ) + xi

e−x − x = 0 → x = e−x , por tanto:

xi+1 = e−xi
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Ejemplos
x2 +3
x2 − 2x + 3 = 0 → x = 2 , por tanto:
x2i +3
xi+1 = 2

sin(x) = 0 → x = sin(x) + x, por tanto:

xi+1 = sin(xi ) + xi

e−x − x = 0 → x = e−x , por tanto:

xi+1 = e−xi
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Ejemplos
x2 +3
x2 − 2x + 3 = 0 → x = 2 , por tanto:
x2i +3
xi+1 = 2

sin(x) = 0 → x = sin(x) + x, por tanto:

xi+1 = sin(xi ) + xi

e−x − x = 0 → x = e−x , por tanto:

xi+1 = e−xi
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Convergencia y divergencia
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Convergencia y divergencia
El métodos de iteración simple de punto fijo establece que:

xi+1 = g(xi )

Considerando que la solución verdadera xr

xr = g(xr )

Restando
xr − xi+1 = g(xr ) − g(xi )
El teorema del valor medio, aplicado a nuestro problema,
establece que si una función g(x) y su derivada son continuas
en el intervalo xi 6 x 6 xr , existe al menos un valor de x = ξ
0
dentro del intervalo para el que: g (ξ) = g(xxrr)−g(x
−xi
i)
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Convergencia y divergencia
Finalmente:
0
xr − xi+1 = g (ξ)(xr − xi )
0
Ev,i+1 = g (ξ)Ev,i

Convergencia y divergencia
0
Si |g | < 1, entonces los errores disminuyen con cada
iteración.
0
Si |g | > 1, entonces los errores aumentan con cada
iteración.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
function p u n t o f i j o v 1 ( g , x0 , EE)
% p u n t o f i j o V 1 : Nombre de l a f u n c i o n
% V a l o r e s de e n t r a d a
% g : f u n c i o n matematica de e n t r a d a ( x = g ( x ) )
% x0 : V a l o r de i n i c i a l
% EE : E r r o r Estimado
% V a l o r e s de s a l i d a
% Raiz
% E r r o r Aproximado
% IM : I t e r a c i o n Maxima
IM =1;
x ( IM ) =x0 ;
x ( IM +1)=g ( x ( IM ) ) ;
EA( IM +1)=abs ( ( x ( IM +1)−x ( IM ) ) / x ( IM +1) ) ∗100;
while EA( IM +1)>EE
IM=IM +1;
x ( IM +1)=g ( x ( IM ) ) ;
EA( IM +1)=abs ( ( x ( IM +1)−x ( IM ) ) / x ( IM +1) ) ∗100;
end
S a l i d a 1 = [ ’ I t e r a c i o n Maxima= ’ , num2str ( IM ) ] ;
S a l i d a 2 = [ x ( 2 : s i z e ( x , 2 ) ) ’ EA ( 2 : s i z e ( x , 2 ) ) ’ ] ;
disp ( ’ ’ )
disp ( S a l i d a 1 )
disp ( ’ ’ )
disp ( ’ Raiz E r r o r Apro ’ )
disp ( S a l i d a 2 )
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> p u n t o f i j o v 1 (@( x ) exp(−x ) , 0 , 0 . 0 1 )

I t e r a c i o n Maxima=19

Raiz E r r o r Apro
1.0000 100.0000
0.3679 171.8282
0.6922 46.8536
0.5005 38.3091
0.6062 17.4468
0.5454 11.1566
0.5796 5.9034
0.5601 3.4809
0.5711 1.9308
0.5649 1.1089
0.5684 0.6244
0.5664 0.3556
0.5676 0.2012
0.5669 0.1143
0.5673 0.0648
0.5671 0.0367
0.5672 0.0208
0.5671 0.0118
0.5672 0.0067
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> p u n t o f i j o v 1 (@( x ) exp(−x ) , 0 , 0 . 0 0 1 )

I t e r a c i o n Maxima=23

Raiz E r r o r Apro
1.0000 100.0000
0.3679 171.8282
0.6922 46.8536
0.5005 38.3091
0.6062 17.4468
0.5454 11.1566
0.5796 5.9034
0.5601 3.4809
0.5711 1.9308
0.5649 1.1089
0.5684 0.6244
0.5664 0.3556
0.5676 0.2012
0.5669 0.1143
0.5673 0.0648
0.5671 0.0367
0.5672 0.0208
0.5671 0.0118
0.5672 0.0067
0.5671 0.0038
0.5671 0.0022
0.5671 0.0012
0.5671 0.0007
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Determinación del coeficiente de arrastre (c)

Calcular los ceros de la función:
gm  c

f(c) = 1 − e− m t − v
c
Considerando m = 68.1 kg, v = 40 m/s y t = 10 s.
667.38
1 − e−0.146843 c − 40 = 0


f (c) =
c
667.38
1 − e−0.146843 c = c


40
667.38
1 − e−0.146843 c


g (c) =
40
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> p u n t o f i j o v 1 (@( x ) 6 6 7 . 3 8 / 4 0 ∗ ( 1 − exp(−0.146843∗x ) ) , 1 4 , 0 . 0 0 1 )

I t e r a c i o n Maxima=8

Raiz E r r o r Apro
1 4 . 54 9 0 3.7738
1 4 . 71 4 5 1.1241
1 4 . 76 1 7 0.3203
1 4 . 77 5 0 0.0900
1 4 . 77 8 8 0.0252
1 4 . 77 9 8 0.0071
1 4 . 78 0 1 0.0020
1 4 . 78 0 2 0.0006
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> p u n t o f i j o v 1 (@( x ) exp(−x ) 6 6 7 . 3 8 / x ∗(1 − exp(−0.146843∗x ) )−40+x , 1 4 , 0 . 0 0 1 )

I t e r a c i o n Maxima=180

Raiz E r r o r Apro
1 5 . 56 8 7 10.0760
1 4 . 07 7 8 10.5907
1 5 . 48 5 7 9.0919
1 4 . 14 7 4 9.4596
1 5 . 41 2 3 8.2065
.................
.................
.................

1 4 . 78 0 1 0.0012
1 4 . 78 0 3 0.0011
1 4 . 78 0 1 0.0011
1 4 . 78 0 3 0.0010
1 4 . 78 0 1 0.0010
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Tópicos

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Fórmula
De la figura se puede apreciar que:

0 f (xi ) − 0
f (xi ) =
xi − xi+1

Despejando xi+1 se llega a:

f (xi )
xi+1 = xi − →Fórmula de Newton-Raphson
f 0 (xi )
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Análisis de errores
La expansión de la serie te Taylor se puede escribir como:
00
0 f (ξ)
f (xi+1 ) = f (xi ) + f (xi )(xi+1 − xi ) + (xi+1 − xi )2 + · · ·
2!
donde ξ se encuentra dentro del intervalo [xi xi+1 ]. Truncando
la serie de Taylor después del término primera derivada, se
obtiene:
0
f (xi+1 ) = f (xi ) + f (xi )(xi+1 − xi )
Como en la intersección con el eje x, f (xi+1 ) = 0, entonces de
la expresión anterior se llega a la fórmula de Newton-Raphson.
Por otro lado, evaluando la serie de Taylor xi+1 = xr (valor
verdadero) y como f (xr ) = 0
00
0 f (ξ)
0 = f (xi ) + f (xi )(xr − xi ) + (xr − xi )2
2!
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Análisis de errores
Restando las expresiones rojas:
00
0 f (ξ)
0 = f (xi )(xr − xi+1 ) + (xr − xi )2
2!
00
f (ξ) 2
0
0 = f (xi )Ev,i+1 +
Ev,i
2!
Finalmente, tanto xi como ξ se deberán aproximar a la raı́z xr :
00
f (xr ) 2
Ev,i+1 = − E
2!f 0 (xr ) v,i

El error es proporcional al cuadrado del error anterior, por lo

cual tenemos una convergencia cuadrática.
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Desventajas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
function newtonraphsonv1 ( f , x0 , EE)
% newtonraphsonv1 : Nombre de l a f u n c i o n
% V a l o r e s de e n t r a d a
% f : f u n c i o n matematica de e n t r a d a ( f ( x ) )
% x0 : V a l o r de i n i c i a l
% EE : E r r o r Estimado
% V a l o r e s de s a l i d a
% Raiz
% E r r o r Aproximado
% IM : I t e r a c i o n Maxima
syms x ;
Df=@( xx ) subs ( d i f f ( f , x ) , x , xx ) ; %Derivada de f : Funcion matematica
IM =1;
r ( IM ) =x0 ;
r ( IM +1)= r ( IM )−f ( r ( IM ) ) / Df ( r ( IM ) ) ;
EA( IM +1)=abs ( ( r ( IM +1)−r ( IM ) ) / r ( IM +1) ) ∗100;
while EA( IM +1)>EE
IM=IM +1;
r ( IM +1)= r ( IM )−f ( r ( IM ) ) / Df ( r ( IM ) ) ;
EA( IM +1)=abs ( ( r ( IM +1)−r ( IM ) ) / r ( IM +1) ) ∗100;
end
S a l i d a 1 = [ ’ I t e r a c i o n Maxima= ’ , num2str ( IM ) ] ;
S a l i d a 2 = [ r ( 2 : s i z e ( r , 2 ) ) ’ EA ( 2 : s i z e ( r , 2 ) ) ’ ] ;
disp ( ’ ’ )
disp ( S a l i d a 1 )
disp ( ’ ’ )
disp ( ’ Raiz E r r o r Apro ’ )
disp ( S a l i d a 2 )
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> newtonraphsonv1 (@( x ) 6 6 7 . 3 8 / x ∗(1 − exp(−0.146843∗x ) ) −40 , 14 , 0.001)

I t e r a c i o n Maxima=3

Raiz E r r o r Apro
1 4 . 75 6 6 5.1274
1 4 . 78 0 2 0.1593
1 4 . 78 0 2 0.0001
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> newtonraphsonv1 (@( x ) x ˆ 1 0 − 1 , 0 . 5 , 0 . 0 0 1 )

I t e r a c i o n Maxima=43

Raiz E r r o r Apro
5 1 . 65 0 0 99.0319
4 6 . 48 5 0 11.1111
4 1 . 83 6 5 11.1111
3 7 . 65 2 9 11.1111
3 3 . 88 7 6 11.1111
.................
.................
.................
1.0237 5.8305
1.0023 2.1299
1.0000 0.2292
1.0000 0.0024
1.0000 0.0000
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Tópicos

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Presentación del método

Fórmula
Un problema potencial en la implementación del método de
Newton-Raphson es la evaluación de la derivada. La derivada
se podrı́a aproximar a la secante, esto serı́a:

0 f (xi−1 ) − f (xi )
f (xi ) ≈
xi−1 − xi

Esta aproximación se conoce como diferencia finita hacia

atrás. Sustituyendo esta expresión en la fórmula de
Newton-Raphson:

(xi−1 − xi )f (xi )
xi+1 = xi −
f (xi−1 ) − f (xi )

La expresión anterior es la fórmula del método de la secante.

Aquı́ se necesitan dos valores de inicio: xi y xi−1
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
function secantev1 ( f , x0 , x1 , EE)
% secantev1 : Nombre de l a f u n c i o n
% V a l o r e s de e n t r a d a
% f : f u n c i o n matematica de e n t r a d a ( f ( x ) )
% x0 : V a l o r de i n i c i a l ( x i −1)
% x1 : V a l o r de i n i c i a l ( xi )
% EE : E r r o r Estimado
% V a l o r e s de s a l i d a
% Raiz
% E r r o r Aproximado
% IM : I t e r a c i o n Maxima
IM =1;
x ( IM ) =x0 ;
x ( IM +1)=x1 ;
x ( IM +2)=x ( IM +1)−f ( x ( IM +1) ) ∗( x ( IM )−x ( IM +1) ) / ( f ( x ( IM ) )−f ( x ( IM +1) ) ) ;
EA( IM +2)=abs ( ( x ( IM +2)−x ( IM +1) ) / x ( IM +2) ) ∗100;
while EA( IM +2)>EE
IM=IM +1;
x ( IM +2)=x ( IM +1)−f ( x ( IM +1) ) ∗( x ( IM )−x ( IM +1) ) / ( f ( x ( IM ) )−f ( x ( IM +1) ) ) ;
EA( IM +2)=abs ( ( x ( IM +2)−x ( IM +1) ) / x ( IM +2) ) ∗100;
end
S a l i d a 1 = [ ’ I t e r a c i o n Maxima= ’ , num2str ( IM ) ] ;
S a l i d a 2 = [ x ( 3 : s i z e ( x , 2 ) ) ’ EA ( 3 : s i z e ( x , 2 ) ) ’ ] ;
disp ( ’ ’ )
disp ( S a l i d a 1 )
disp ( ’ ’ )
disp ( ’ Raiz E r r o r Apro ’ )
disp ( S a l i d a 2 )
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Programa MATLAB

Programa MATLAB
>> secantev1 (@( x ) 6 6 7 . 3 8 / x ∗(1 − exp(−0.146843∗x ) ) − 4 0 , 1 5 , 1 6 , 0 . 0 0 1 )

I t e r a c i o n Maxima=4

Raiz E r r o r Apro
1 4 . 76 9 6 8.3306
1 4 . 78 0 7 0.0752
1 4 . 78 0 2 0.0035
1 4 . 78 0 2 0.0000
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Tópicos

1 INTRODUCCIÓN
Métodos abiertos
2 ITERACIÓN SIMPLE DE PUNTO FIJO
Presentación del método
Programa MATLAB
3 NEWTON-RAPHSON
Presentación del método
Programa MATLAB
4 SECANTE
Presentación del método
Programa MATLAB
5 Problemas
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Problema 1
Determine la raı́z real más grande
f (x) = 0.95x3 − 5.9x2 + 10.9x − 6:
a-) En forma gráfica.
b-) Con el uso del método de Newton-Raphson (tres
iteraciones, xi = 3.5).
c-) Con el uso de la secante (tres iteraciones,
xi−1 = 2.5 y xi = 3.5).
INTRODUCCIÓN ITERACIÓN SIMPLE DE PUNTO FIJO NEWTON-RAPHSON SECANTE Problemas

Problema 2
Utilice el método de Newton-Raphson para encontrar la raı́z de
f (x) = e−0.5x (4 − x) − 2. Utilice los valores iniciales de a) 2, b)
6 y c) 8. Explique el resultado.