Ecuaciones no lineales Métodos
Cerrados 1
Competencias
Comprender la interpretación gráfica de una raíz.
Entender la diferencia entre los métodos cerrados y los métodos abiertos para la localización de
las raíces.
Utilizar los métodos cerrados de la bisección y falsa posición para encontrar las raíces de una
ecuación no lineal.
Aplicar los métodos abiertos para encontrar las raíces de una ecuación no lineal.
Identificar cuando una función es contractiva.
2
Raíz de una función
Muchos problemas de ciencia e ingeniería son expresados como:
Dada una función continua 𝑓(𝑥),
encuentre el valor de 𝑟 tal que 𝑓(𝑟) = 0
Estos problemas se denominan problemas de búsqueda de raíz.
3
Raíces de una ecuación
Sea f (x) una función de una variable real. Cualquier número r para el cual f (r) = 0 se llama
cero de la función.
Raíces (ceros)
r1 r2 r3
4
Raíces de ecuaciones: métodos
Raíces de
Ecuaciones
Métodos Cerrados Gráficos Métodos abiertos
Bisección Punto fijo
Secante
Falsa Posición Newton Raphson
5
Métodos gráficos
• Un método simple para obtener una aproximación a la raíz
de la ecuación f (x) = 0 consiste en graficar la función y
observar dónde cruza el eje x. Este punto representa el
valor de x para el cual f (x) = 0, ofrece una aproximación
inicial de la raíz.
6
Métodos gráficos
No hay raíces
Una sola raíz
Dos raíces
Tres raíces
7
Raíces de ecuaciones no lineales
Raíces de
Ecuaciones
Métodos Cerrados Gráficos Métodos abiertos
Bisección Punto fijo
Newton Raphson
Falsa Posición Secante
8
Métodos cerrados
• Los métodos cerrados o de intervalos, son aquellos que necesitan dos
valores iniciales para la raíz y dichos valores deben encerrar o estar a
ambos lados de la raíz.
• La idea es reducir sistemáticamente el tamaño del intervalo y así
converger a la respuesta correcta.
9
Método de la bisección
• El método de bisección es el más simple de todos los métodos para
encontrar una raíz de una ecuación no lineal. Se comienza con un
intervalo que contiene una raíz y se divide en una mitad izquierda y
una mitad derecha y se decide cuál de las dos mitades contiene la raíz
con el intervalo seleccionado se continua con la división. Se hace esto
hasta que el intervalo que contiene la raíz se reduce lo suficiente para
cumplir con el nivel de precisión requerido.
10
Teorema del valor intermedio
• Sea f(x) definida en un intervalo [a,b].
f(a)
• Teorema del valor intermedio:
Si la función es continua y f(a) y f(b)
tienen diferentes signos entonces la a b
función tiene al menos un cero en el
intervalo [a,b]. f(b)
11
Método de la bisección: algoritmo
• Paso 1: Encuentre un intervalo [a,b] donde f(a)f(b) < 0, o en otras palabras donde
f(x) cambie de signos entre f(a) y f(b)
f(a)
x b
a
f(b)
12
Método de la bisección: algoritmo
• Paso 2, estime la raíz de la ecuación f(x)=0, se determina
𝑎+𝑏
𝑥= f(a)
2
x b
a
f(x)
f(b)
13
Método de la bisección: algoritmo
• Paso 3, Verifique
➢ Si f(a)f(x) < 0, entonces la raíz se encuentra dentro del
subintervalo izquierdo. Entonces b=x.
➢ Si f(a)f(x) > 0, entonces la raíz se encuentra dentro del f(a)
subintervalo derecho. Entonces a=x.
x b
➢ Si f(a)f(x) = 0, la raíz es igual a x; termina el cálculo.
a
f(b)
14
Método de la bisección: criterio de finalización
• Paso 4, Encuentre el error relativo
𝑥𝑛𝑢𝑒𝑣𝑜 −𝑥𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝜀𝑎 = *100
𝑥𝑛𝑢𝑒𝑣𝑜
Vaya al Paso 2 usando
nuevas conjeturas
Si
superiores e inferiores.
Es a s ?
No Detener el algoritmo
15
Método de la bisección: criterio de finalización
• Una ventaja del método de la bisección es que después de cada iteración el intervalo
se reduce a la mitad, después de 2 iteraciones, el intervalo original se ha reducido 2𝑛
veces. Por lo tanto, si el tamaño del intervalo original es 𝑎 y el criterio de
convergencia aplicado es la diferencia del valor absoluto de dos 𝑥 consecutivos es 𝜀,
entonces se requieren n iteraciones, donde n se calcula con la igualdad
(𝑏 − 𝑎)
𝑛
≤𝜀
2
despejando
ln 𝑏 − 𝑎 − ln(ε)
𝑛=
ln(2)
16
Método de la bisección: algoritmo
𝑎+𝑏 f(x)f(a) > 0 -> NO
𝑥=
2 f(x)f(a) < 0 -> b→ x
f(a)
b
x
xb
a f(x)
f(x)
f(b)
f(b)
17
Método de la bisección: Ejemplo
Encuentre la raíz de la ecuación 8-4.5(x-sinx) con una tolerancia menor ℇ=0.03
Solución:
➢ Se asume como límites a=2 y b= 3
➢ Se verifica que haya cambio de signo
➢ f(2) = 8-4.5*(2-sen(2))= 3.091
➢ f(3) =8-4.5*(3-sen(3))= -4.86
Por lo tanto, al menos una raíz se encuentra entre a y b
Método de la bisección: Ejemplo
El gráfico demuestra que es correcto el cambio de signo
0
2 2.2 2.4 2.6 2.8 3 3.2
-2
-4
-6
-8
-10
19
Ejemplo
Iteración 1
𝑎+𝑏 2+3
La estimación de la raíz es 𝑥= = = 2.5
2 2
𝑓 𝑥 = 𝑓 2.5 = 8 − 4.5 ∗ 2.5 − 𝑠𝑒𝑛𝑜 2.5 = −0.5568
𝑓 𝑎 𝑓 𝑥 = 𝑓 2 𝑓 2.5 = 3.091 −0.5568 < 0 SI
𝑓 𝑏 𝑓 𝑥 = 𝑓 3 𝑓 2.5 = −4.86 −0.5568 < 0 NO
Por lo tanto, la raíz está entre a y x , es decir, entre 2 y 2.5. Por lo tanto, los límites inferior y
superior nuevos son
𝑎 = 2, 𝑏 = 2.5
En este punto, el error aproximado relativo absoluto no se puede calcular ya que no tenemos
una aproximación previa.
20
Método de la bisección: Ejemplo
f(a)
3
a x b
0
2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3
-1 f(X)
-2
f(b)
-3
-4
-5
-6
21
Ejemplo
Iteración 2
𝑎 + 𝑏 2 + 2.5
La estimación de la raíz es 𝑥= = = 2.25
2 2
𝑓 𝑥 = 𝑓 2.25 = 8 − 4.5 ∗ 2.25 − 𝑠𝑒𝑛𝑜 2.25 = 1.3763
𝑓 𝑎 𝑓 𝑥 = 𝑓 2 𝑓 2.25 = 3.091 1.37 < 0 𝑁𝑂
𝑓 𝑏 𝑓 𝑥 = 𝑓 2.5 𝑓 2.25 = −0.5568 1,3763 < 0 SI
Por lo tanto, la raíz está entre b y x , es decir, entre 2.5 y 2.25. Por lo tanto, los límites inferior
y superior nuevos son
𝑎 = 2.25, 𝑏 = 2.5
22
Método de la bisección: Ejemplo
23
Ejemplo
El error aproximado relativo absoluto al final de la iteración 2 es
𝑥 𝑎𝑐𝑡𝑢𝑎𝑙 − 𝑥 𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
∈𝑎 = 𝑎𝑐𝑡𝑢𝑎𝑙
× 100
𝑥
2.25 − 2.5
= × 100
2.25
= 11.11%
24
Ejemplo
Iteracion 3
𝑎 + 𝑏 2.25 + 2.5
La estimación de la raíz es 𝑥= = = 2.375
2 2
𝑓 𝑥 = 𝑓 2.375 = 8 − 4.5 ∗ 2.375 − 𝑠𝑒𝑛𝑜 2.375 = 0.4340
𝑓 𝑎 𝑓 𝑥 = 𝑓 2.25 𝑓 2.375 = 1.376 0.4340 < 0 NO
𝑓 𝑏 𝑓 𝑥 = 𝑓 2.5 𝑓 2.375 = −0.5568 0.4340 < 0
Por lo tanto, la raíz está entre b y x , es decir, entre 2.5 y 2.375. Por lo tanto, los límites inferior y
superior nuevos son
𝑎 = 2.375, 𝑏 = 2.5
25
Método de la bisección: Ejemplo
2.5
1.5 f(a)
0.5 f(x)
a x b
0
2.25 2.3 2.35 2.4 2.45 2.5
-0.5
f(b)
-1
-1.5
-2
26
Método de la bisección: Ejemplo
➢ Se repite el algoritmo hasta que el error relativo sea menor o igual a ℇ
a b x f(a) f(b) f(x) f(a)*f(x) f(b)*f(x) error
2,00000 3,00000 2,50000 3,09 -4,86 -0,56 -1,72 2,71
2,00000 2,50000 2,25000 3,09 -0,56 1,38 4,26 -0,77 11,11
2,25000 2,50000 2,37500 1,38 -0,56 0,43 0,60 -0,24 5,26
2,37500 2,50000 2,43750 0,43 -0,56 -0,06 -0,02 0,03 2,56
2,37500 2,43750 2,40625 0,43 -0,06 0,19 0,08 -0,01 1,30
2,40625 2,43750 2,42188 0,19 -0,06 0,07 0,01 0,00 0,65
2,42188 2,43750 2,42969 0,07 -0,06 0,01 0,00 0,00 0,32
2,42188 2,42969 2,42578 0,07 0,01 0,04 0,00 0,00 0,16
La raíz es 2.42578 ± 0.00391
27
Método de la falsa posición
• Si una raíz real está limitada por a
y b de f (x) = 0, entonces se puede
aproximar la solución haciendo
una interpolación lineal entre los f(b)
puntos [a, f (a)] y [b, f (b)]. La
intersección de esta línea con el
eje de las x representa una mejor
aproximación de la raíz
xr
𝑓(𝑎) 𝑓(𝑏) a
=
𝑥𝑟 − 𝑎 𝑥𝑟 − 𝑏 b
f(a)
28
Método de la falsa posición
1- Encuentre un intervalo [a,b] donde f(x) = 0
2- Estime el valor de la raíz con la siguiente formula
𝑓 𝑏 (𝑎 − 𝑏) 𝑓(𝑎) 𝑓(𝑏)
𝑥=𝑏− =
𝑓(𝑎) − 𝑓(𝑏) 𝑥𝑟 − 𝑎 𝑥𝑟 − 𝑏
3- Evalúe f(x),
➢ Si f(a)f(x) < 0, entonces la raíz se encuentra dentro del subintervalo izquierdo. Entonces b=x.
➢ Si f(a)f(x) > 0, entonces la raíz se encuentra dentro del subintervalo derecho. Entonces a=x.
➢ Si f(a)f(x) = 0, la raíz es igual a x; termina el cálculo.
29