ez

nd
ná
er
Análisis numérico

H
Soluciones de ecuaciones de una variable

l
be
Dra. Anabel Hernández Ramı́rez
na
A
25 de enero de 2022
.
ra
D
 ez
El problema de la búsqueda de raı́ces

nd
ná
er
H
Consiste en obtener una raı́z, o solución, de una ecuación de la

l
forma f (x ) = 0 para una función dada. Una raı́z de esta ecuación

be
es llamada también un cero de la función f .
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

2 / 40
 ez
nd
ná
er
Método de bisección

H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

3 / 40
 ez
Preliminares

nd
ná
Teorema 1

er
Si f denota una función continua sobre un intervalo cerrado [a, b]

H
para el cual y si k es cualquier número entre f (a) y f (b), entonces
existe por lo menos un número x0 entre a y b tal que f (x0 ) = k.

l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

4 / 40
 ez
Método de bisección

nd
ná
• Se basa en el teorema del valor intermedio.

er
• Supongamos que f es una función continua definida en el

H
intervalo [a, b] con f (a) y f (b) de signos diferentes. De
acuerdo con el Teorema 1, existe un número p en (a, b) tal

l
que f (p) = 0.

be
• Supongamos (por razones de simplicidad) que la raı́z de este
intervalo es única.
na
A
• El método requiere dividir varias veces a la mitad (o bisectar)
los subintervalos de [a, b] y, en cada paso, localizar la mitad
.
ra

que contenga a p.
D

Dra. Anabel Hernández Ramı́rez

5 / 40
 ez
Pseudocódigo

nd
ná
Para empezar, supongamos que a1 = a y b1 = b, y sea p1 el punto
medio de [a, b]; es decir

er
b1 − a1 a1 + b1

H
p1 = a1 + = .
2 2

l
• Si f (p1 ) = 0, entonces p = p1 y listo.

be
• Si f (p1 ) 6= 0, entonces f (p1 ) tiene el mismo signo que f (a1 ) o
f (b1 ). na
• Si f (p1 ) y f (a1 ) tienen el mismo signo, entonces p ∈ (p1 , b1 ) y
A
entonces a2 = p1 y b2 = b1 .
• Si f (p1 ) y f (a1 ) tienen signos opuestos, entonces p ∈ (a1 , p1 )
.
ra

y entonces a2 = a1 y b2 = p1 .
D

Después volvemos a aplicar el proceso al intervalo [a2 , b2 ].

Dra. Anabel Hernández Ramı́rez

6 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

7 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

8 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

9 / 40
 ez
Ejemplo 1

nd
ná
er
Demuestre que f (x ) = x 2 − 2 tiene una raı́z en [0,2] y utilice el

H
método de bisección para determinar una aproximación a la raı́z

l
con una tolerancia de |f (x )| < 10−1 = 0.1 y |f (x )| < 10−2 =

be
0.01. Verifica que los resultados son cercanos a la raı́z, evaluando la
na
aproximación a la raı́z en la función f .
¿Qué pasarı́a si se cambia el intervalo por [2,4]?
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

10 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

11 / 40
 ez
Ejemplo 2

nd
ná
er
H
Demuestre que f (x ) = x 3 + 4x 2 − 10 tiene una raı́z en [1,2] y utilice

l
el método de bisección para determinar una aproximación a la raı́z

be
con una tolerancia de |f (x )| < 10−4 .
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

12 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

13 / 40
 ez
Ejemplo 3

nd
ná
er
H
√
Demuestre que f (x ) = x − cos(x ) tiene una raı́z en [0,1] y utilice

l
el método de bisección para determinar una aproximación a la raı́z

be
con una tolerancia de |f (x )| < 10−3 .
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

14 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

15 / 40
 ez
Ejercicios

nd
ná
er
Use el método de bisección para encontrar una solución aproximada
con una tolerancia de |f (x )| < 10−5 para las siguientes ecuaciones

H
NO LINEALES.

l
1 x − 2−x = 0 para 0 ≤ x ≤ 1.

be
2 e x − x 2 + 3x + 2 = 0 para −1 ≤ x ≤ 0.
3 na
2x cos(2x ) − (x + 1)2 = 0 para −3 ≤ x ≤ −2 y −1 ≤ x ≤ 0.
A
4 x cos x − 2x 2 + 3x − 1 = 0 para 0.2 ≤ x ≤ 0.3 y
1.2 ≤ x ≤ 1.3.
.
ra
D

Dra. Anabel Hernández Ramı́rez

16 / 40
 ez
Respuestas ejercicios

nd
ná
er
El método de bisección da como resultado:

H
1 p̃ = 0.6412353515625.
p̃ = −0.68084716796875.

l
2

be
3 p̃ = −2.191314697265625 para −3 ≤ x ≤ −2 y
p̃ = −0.79815673828125 para −1 ≤ x ≤ 0.
na
4 p̃ = 0.29755859375000004 para 0.2 ≤ x ≤ 0.3 y
A
p̃ = 1.2566406250000002 para 1.2 ≤ x ≤ 1.3.
.
ra
D

Dra. Anabel Hernández Ramı́rez

17 / 40
 ez
Desventajas y ventajas del método

nd
ná
er
1 El método de bisección converge lentamente (es decir, N

H
puede volverse muy grande antes de que |p − pn | sea bastante
pequeño) e inadvertidamente se puede desechar una buena

l
be
aproximación intermedia.
2 El procedimiento siempre converge en una solución y, por esta
na
razón, a menudo sirve para inciar otros métodos más
A
eficientes.
.
ra
D

Dra. Anabel Hernández Ramı́rez

18 / 40
 ez
Cota del error de aproximación

nd
ná
er
Teorema 2

H
Supongamos que f ∈ C[a, b] y f (a) · f (b). El método de bisección

l
genera una sucesión {pn }∞
n=1 que aproxima un cero p de f tal que

be
b−a
na
|pn − p| ≤
2n
, cuando n ≥ 1.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

19 / 40
 ez
nd
ná
er
Método de Newton

H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

20 / 40
 ez
Preliminares: Polinomio de Taylor de

nd
grado n-ésimo de f en a.

ná
er
H
Definición 1

l
El polinomio de Taylor de n−ésimo grado de f en a, se define como

be
n
f (i) (a)
na
X

n=i
i!
(x − a)i
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

21 / 40
 ez
Preliminares. Ejemplo.

nd
Considera f (x ) = ex , entonces

ná
x2 x2 x3
T1 (x ) = 1 + x , T2 (x ) = 1 + x + , T3 (x ) = 1 + x + + .

er
2! 2! 3!

H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

22 / 40
 ez
Preliminares. Residuo.

nd
ná
er
Teorema 3

H
Si f (x ) = T (x ) + R(x ) donde T (x ) es el polinomio de Taylor de

l
grado n de f en a y

be
lı́m R(x ) = 0
n→∞
na
para |x − a| < R entonces f es igual a la suma de sus series de
A
Taylor en el intervalo |x − a| < R.
.
ra
D

Dra. Anabel Hernández Ramı́rez

23 / 40
 ez
Preliminares. Desigualdad de Taylor.

nd
ná
er
Teorema 4

H
Si |f (n) (x )| ≤ M para |x − a| ≤ d entonces el residuo Rn (x ) de la

l
serie de Taylor cumple con la desigualdad

be
M
|Rn (x )| ≤ na
(n + 1)!
|x − a|n+1 para |x − a| ≤ d.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

24 / 40
 ez
Método de Newton

nd
ná
er
El método de Newton (o método de Newton-Raphson) es una técni-

H
ca de iteración funcional pn = g(pn−1 ) para la cual

l
be
f (pn−1 )
g(pn−1 ) = pn−1 − , para n ≥ 1.
na f 0 (pn−1 )
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

25 / 40
 ez
Ejemplos

nd
ná
er
Use el método de Newton para encontrar una solución aproximada

H
con una tolerancia de |f (x )| < tol y una aproximación inicial para
una solución de f (x ) = 0 denotada por x0 .

l
be
1 x 2 − 2 = 0 para tol = 1 × 10−6 y x0 = 1.5.
2 na
x 3 − x 2 − 1 = 0 para tol = 1 × 10−10 y x0 = 1.
3 f (x ) = x 1/3 para tol = 1 × 10−2 = 0.01 y x0 = 0.1.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

26 / 40
 ez
Ejercicios

nd
ná
er
Use el método de Newton para encontrar una solución aproximada

H
con una tolerancia de |f (x )| < tol y una aproximación inicial para
una solución de f (x ) = 0 denotada por x0 .

l
be
1 x 2 − 2 = 0 para tol = 1 × 10−6 y x0 = 1.5.
2 na
x 3 − x 2 − 1 = 0 para tol = 1 × 10−10 y x0 = 1.
3 f (x ) = x 1/3 para tol = 1 × 10−2 = 0.01 y x0 = 0.1.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

27 / 40
 ez
nd
ná
er
Análisis del error para métodos

H
iterativos

l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

28 / 40
 ez
Definición. Orden de convergencia.

nd
ná
er
Suponga que {pn }∞
n=0 es una sucesión que converge a p, con

H
pn = p para todo n. Si existen constantes positivas λ y α tales que

l
|pn+1 − p|

be
lı́m = λ,
n→∞ |pn − p|α

entonces {pn }∞
na
n=0 converge a p de orden α, con un error asintótico
A
constante λ.
.
ra
D

Dra. Anabel Hernández Ramı́rez

29 / 40
 ez
nd
ná
• Una técnica iterativa de la forma pn = g(pn−1 ) se dice de

er
orden α si la sucesión {pn }∞
n=0 converge a la solución
p = g(p) de orden α.

H
• En general, una sucesión con orden de convergencia de orden

l
mayor converge más rápidamente que una de orden más bajo.

be
Dos tipos de orden son importantes:
Si α = 1 (y λ < 1), la sucesión es linealmente convergente.
1
2
na
Si α = 2, la sucesión es cuadráticamente convergente.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

30 / 40
 ez
Preliminares. Raı́ces múltiples.

nd
ná
er
H
El método de Newton, generalmente, dará problemas si f 0 (p) = 0

l
cuando f (p) = 0. Para examinar estas dificultades en más detalle,

be
se darán a continuación algunas definiciones.
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

31 / 40
 ez
Definición. Cero de multiplicidad m.

nd
ná
er
H
Una solución p de f (x ) = 0 es un cero de multiplicidad m de f ,

l
si para x 6= p, podemos escribir f (x ) = (x − p)m q(x ), donde

be
lı́m q(x ) 6= 0.
x →p
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

32 / 40
 ez
Teoremas

nd
ná
er
1 La función f ∈ C 1 [a, b] tiene un cero simple en p ∈ (a, b), si y

H
sólo si, f (p) = 0 pero f 0 (p) 6= 0.

l
2 La función f ∈ C m [a, b] tiene un cero de multiplicidad m en

be
p ∈ (a, b), si y sólo si,
na
0 = f (p) = f 0 (p) = f 00 (p) = · · · = f (m−1) (p), pero
f (m) (p) 6= 0.
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

33 / 40
 ez
Ejemplo 1

nd
ná
er
H
Sea f (x ) = e x − x − 1. (a) Muestre que la función f tiene un cero

l
de multiplicidad 2 en x = 0. (b) Muestre que el método de Newton

be
con p0 = 1 converge a 0 pero no cuadráticamente.
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

34 / 40
 ez
nd
ná
er
H
l
be
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

35 / 40
 ez
Tabla de Ejemplo 1

nd
n pn
0 1.0

ná
1 0.58198
2 0.31906

er
3 0.16800

H
4 0.08635
5 0.04380

l
be
6 0.02206
7 0.01107
8 0.005545 na
2.7750 × 10−3
A
9
10 1.3881 × 10−3
.

6.9411 × 10−4
ra

11
12 3.4703 × 10−4
D

13 1.7416 × 10−4
14 8.8041 × 10−5
Dra. Anabel Hernández Ramı́rez

15 4.2610 × 10−5
36 / 40
 ez
Método de Newton modificado

nd
ná
er
El método de Newton modificado (o método de Newton-Raphson

H
modificado) es una técnica de iteración funcional pn = g(pn−1 ) para

l
la cual

be
f (pn−1 )f 0 (pn−1 )
g(pn−1 ) = pn−1 − na , para n ≥ 1.
(f 0 (pn−1 ))2 − f (pn−1 )f 00 (pn−1 )
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

37 / 40
 ez
Ejemplo 2

nd
ná
er
H
En el Ejemplo 1, se muestra que f (x ) = e x − x − 1 tiene un ce-

l
ro de multiplicidad 2 en x = 0 y que el método de Newton con

be
p0 = 1 converge a 0 pero no cuadráticamente. Muestre que Newton
modificado mejora la tasa de convergencia.
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

38 / 40
 ez
Tabla de Ejemplo 2

nd
ná
er
n pn

H
0 1.0
1 −2.3421061 × 10−1

l
−8.4582788 × 10−3

be
2
3 −1.1889524 × 10−5
4 na
−6.8638230 × 10−6
5 −2.8085217 × 10−7
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

39 / 40
 ez
Bibliografı́a

nd
ná
er
H
Burden, Richard L. y Faires, J. Douglas. Análisis numérico. CEN-

l
be
GAGE Learning. Novena edición. 2011.
na
A
.
ra
D

Dra. Anabel Hernández Ramı́rez

40 / 40