Errores
Raı́ces de Ecuaciones
Errores y Raı́ces de Ecuaciones
Omar Andrés Cuervo Fernández
Universidad de la Sabana
Agosto 2020
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Contenido
1 Errores
Tipos de errores
Errores representativos
Cifras significativas decimales
2 Raı́ces de Ecuaciones
Métodos de solución
Métodos Cerrados
Métodos Abiertos
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores
Observaciones
1 Los métodos numéricos proveen una aproximación de la
solución exacta, llamada solución analı́tica (si existe).
2 Toda aproximación posee un error propio al método.
3 Ası́ una solución no exista, se puede estudiar un tipo de error.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Tipos de Errores
Existen dos tipos de errores principales
1 Errores de redondeo: Son aquellos dados por cálculos
realizados con números que no son almacenados por la
calculadora.
e = 2,718281828459045235360 · · ·
2 Errores de truncamiento: Dados por los métodos númericos
como métodos de aproximación a soluciones matemáticas
exactas. Ej: Problema del paracaidista.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Tipos de Errores
Existen dos tipos de errores principales
1 Errores de redondeo: Son aquellos dados por cálculos
realizados con números que no son almacenados por la
calculadora.
e = 2,718281828459045235360 · · ·
2 Errores de truncamiento: Dados por los métodos númericos
como métodos de aproximación a soluciones matemáticas
exactas. Ej: Problema del paracaidista.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores de redondeo
Ejemplos
Insertamos los números 13 y 17 en la lı́nea de comandos de Matlab y
en la calculadora. Comparamos:
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores de redondeo
Números decimales finitos en el computador
Observaciones
1 El comando format long en Matlab, nos permite observar
15 cifras decimales del número. A la última cifra significativa
se le hace el redondeo.
2 No todo número real se puede representar en el computador.
Los decimales infinitos, se deben redondear (a 15 cifras
decimales). Al número redondeado se le llama notación
flotante del número fl(x).
Ejemplo Determinar en Matlab fl(π).
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores de Redondeo
Error de redondeo
El error de redondeo que se genera al utilizar el número fl(x) en
lugar de x es
|x − fl(x)| M
εredondeo = ≤ , (1)
|x| 2
donde M es el épsilon de la máquina. Escribir el comando eps en
la lı́nea de comandos de Matlab.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores de Redondeo
Observación Importante
Se debe tener cuidado con la propagación del error de rendon-
deo. En un número muy grande de cálculos, el error se puede
volver significativo.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Errores
Definimos los errores dados por las aproximaciones:
1 Si se conoce el valor verdadero
1 Error verdadero: Ev = valor verdadero - valor aproximado.
2 Error relativo porcentual verdadero:
|valor verdadero - valor aproximado|
εv = × 100 %
valor verdadero
2 Si se desconoce el valor verdadero
Error relativo porcentual aproximado:
|aproximación actual - aproximación anterior|
εa = × 100 %
aproximación actual
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Ejemplo del paracaidista
Ejercicio
Crear una rutina (o hacer en calculadora) para determinar una tabla
del error relativo porcentual verdadero y el error relativo porcen-
tual aproximado del ejemplo del paracaidista, con los dos pasos de
tiempo.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Ejemplo del paracaidista
Resultados
t 0 2 4 6 8 10
S. Analı́tica 0 15.42 24.77 30.45 33.89 35.98
S. Numérica 0 17.15 26.79 32.22 35.27 36.99
E. R. P. Verdadero 0 11.19 8.14 5.81 4.07 2.80
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Cifras Significativas
Teorema (Scarborough-1966)
Si el error relativo porcentual (verdadero o aproximado) es menor
que el valor de tolerancia εs = (0,5×102−n ) %, entonces el resultado
es correcto en al menos n cifras significativas.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Tipos de errores
Errores
Errores representativos
Raı́ces de Ecuaciones
Cifras significativas decimales
Ejercicios
1 Con la serie de Taylor de f (x) = ln(x) alrededor de x = 1,
determinar el número de términos de la serie para aproximar el
valor de ln(2) con 5 cifras significativas.
2 Desarrollar un programa donde se agregue términos de la serie
de Maclaurin de la función f (x) = e x , calculando el valor de
error relativo verdadero hasta que sea menor al 0,01 % de e 4
(leer ejercicio 4.2, página 103 del libro). Repita el ejercicio
para un error relativo aproximado menor a 0,01 %. ¿Cuál
converge más rápido?
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Contenido
1 Errores
Tipos de errores
Errores representativos
Cifras significativas decimales
2 Raı́ces de Ecuaciones
Métodos de solución
Métodos Cerrados
Métodos Abiertos
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Raı́ces de Ecuaciones
Algunas Observaciones
Polinomios de grado mayor a 2, es complicado determinar la
raı́z de la ecuación.
En algunas funciones trascendentales (no algebráicas), no
existe método analı́tico para determinar la solución de la
ecuación. Ej: f (x) = e −x + x
Las variables dentro de las ecuaciones que no se pueden
despejar se dicen que están en forma implı́cita. Ej: Despejar c
de la ecuación
mg c
v (t) = 1 − e− m t .
c
Las raı́ces de una función pueden ser reales o complejas. Ej:
x 2 + 1 = 0. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Raı́ces de Ecuaciones
Algunas Observaciones
Polinomios de grado mayor a 2, es complicado determinar la
raı́z de la ecuación.
En algunas funciones trascendentales (no algebráicas), no
existe método analı́tico para determinar la solución de la
ecuación. Ej: f (x) = e −x + x
Las variables dentro de las ecuaciones que no se pueden
despejar se dicen que están en forma implı́cita. Ej: Despejar c
de la ecuación
mg c
v (t) = 1 − e− m t .
c
Las raı́ces de una función pueden ser reales o complejas. Ej:
x 2 + 1 = 0. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Raı́ces de Ecuaciones
Algunas Observaciones
Polinomios de grado mayor a 2, es complicado determinar la
raı́z de la ecuación.
En algunas funciones trascendentales (no algebráicas), no
existe método analı́tico para determinar la solución de la
ecuación. Ej: f (x) = e −x + x
Las variables dentro de las ecuaciones que no se pueden
despejar se dicen que están en forma implı́cita. Ej: Despejar c
de la ecuación
mg c
v (t) = 1 − e− m t .
c
Las raı́ces de una función pueden ser reales o complejas. Ej:
x 2 + 1 = 0. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Raı́ces de Ecuaciones
Algunas Observaciones
Polinomios de grado mayor a 2, es complicado determinar la
raı́z de la ecuación.
En algunas funciones trascendentales (no algebráicas), no
existe método analı́tico para determinar la solución de la
ecuación. Ej: f (x) = e −x + x
Las variables dentro de las ecuaciones que no se pueden
despejar se dicen que están en forma implı́cita. Ej: Despejar c
de la ecuación
mg c
v (t) = 1 − e− m t .
c
Las raı́ces de una función pueden ser reales o complejas. Ej:
x 2 + 1 = 0. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos de Solución
Tipos de Métodos
Son de dos tipos:
1 Métodos Cerrados: Utilizan un intervalo donde se encuentre
la raı́z para su ejecución.
1 Método de la Bisección.
2 Método de la Falsa Posición.
2 Métodos Abiertos: Requieren uno o varios puntos iniciales
cercanos a la raı́z.
1 Método de Newton-Rhapson.
2 Método de la Secante.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos de Solución
Tipos de Métodos
Son de dos tipos:
1 Métodos Cerrados: Utilizan un intervalo donde se encuentre
la raı́z para su ejecución.
1 Método de la Bisección.
2 Método de la Falsa Posición.
2 Métodos Abiertos: Requieren uno o varios puntos iniciales
cercanos a la raı́z.
1 Método de Newton-Rhapson.
2 Método de la Secante.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Bisección
Conociendo un intervalo [a, b] donde pertenece la raı́z (se verifica
haciendo f (a)f (b) < 0) aplicamos el siguiente algoritmo de 2 pasos:
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Bisección
1 Una aproximación de la raı́z es el punto medio del intervalo
a+b
xm = .
2
2 Determinar en qué subintervalo está la raı́z:
1 Si f (a)f (xm ) = 0, entonces la raı́z es xm .
2 Si f (a)f (xm ) < 0, entonces la raı́z está en el subintervalo
[a, xm ].
3 Si f (a)f (xm ) > 0, entonces la raı́z está en el subintervalo
[xm , b].
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Bisección
1 Con el comando ezplot podemos obtener un gráfico rápido
de una función en un intervalo [a, b]. Ej: ezplot(’x 2 − 1’,
[-1,1]).
2 Este método siempre es convergente, es decir que en cada
iteración se acerca al valor de la raı́z.
3 Generalmente se trabaja con situaciones donde no se conoce
el valor real de la raı́z, por lo que debemos utilizar el error
aproximado para verificar que el método es convergente.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Bisección
1 Con el comando ezplot podemos obtener un gráfico rápido
de una función en un intervalo [a, b]. Ej: ezplot(’x 2 − 1’,
[-1,1]).
2 Este método siempre es convergente, es decir que en cada
iteración se acerca al valor de la raı́z.
3 Generalmente se trabaja con situaciones donde no se conoce
el valor real de la raı́z, por lo que debemos utilizar el error
aproximado para verificar que el método es convergente.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Bisección
1 Con el comando ezplot podemos obtener un gráfico rápido
de una función en un intervalo [a, b]. Ej: ezplot(’x 2 − 1’,
[-1,1]).
2 Este método siempre es convergente, es decir que en cada
iteración se acerca al valor de la raı́z.
3 Generalmente se trabaja con situaciones donde no se conoce
el valor real de la raı́z, por lo que debemos utilizar el error
aproximado para verificar que el método es convergente.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de la Bisección
Ejemplo
Aproximar la raı́z de la función f (x) = e −x −x por medio del método
de la Bisección. Grafique la función para determinar un intervalo
inicial.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de la Bisección
Ejercicio
Implementar un algoritmo en Matlab donde determine la raı́z de
la función f (x) = e −x − x por medio del método de la Bisección.
Grafique la función, determine el intervalo inicial y pare el cálculo
cuando el error sea menor a 0,05 %.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Cerrados
Método de la Falsa Posición
También conocido como el método de la interpolación lineal, apro-
vecha la magnitud de los valores f (a) y f (b). Tiene como fórmulas
bf (a) − af (b)
xaprox = ,
f (a) − f (b)
ó
f (b)(a − b)
xaprox = b − .
f (a) − f (b)
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de la Falsa Posición
Ejemplo
Aproximar la raı́z de la función f (x) = e −x −x por medio del método
de la Bisección. Grafique la función para determinar un intervalo
inicial.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de la Falsa Posición
Ejemplo
Implementar un algoritmo en Matlab donde determine la raı́z de la
función f (x) = e −x − x por medio del método de la falsa posición.
Grafique la función, determine el intervalo inicial y pare el cálculo
cuando el error sea menor a 0,05 %. Compare los resultados con los
dados por el método de la Bisección.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Abiertos
Método de Newton-Rhapson
Método donde se toma un valor inicial xi aproximado a la raı́z y se
traza una tangente a la curva en el punto (xi , f (xi )).
El punto donde la tangente toca el eje x (nombrado xi+1 ) es una
aproximación de la raı́z. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Abiertos
Método de Newton-Rhapson
Observaciones
El valor de xi+1 se deduce de la gráfica
f (xi )
xi+1 = xi − .
f 0 (xi )
El método de Newton-Raphson (en general los métodos
abiertos) no siempre son convergentes. Esto depende de la
elección del punto inicial xi .
En caso de ser convergentes los métodos abiertos, lo hacen
más rápido que los métodos cerrados. Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de Newton-Rhapson
Ejemplo
1 Con una calculadora, determinar el valor de la raı́z de
f (x) = e −x − x, tomando como punto inicial x0 = 1. Hacer la
iteración 3 veces y mirar como va variando el error relativo.
Compare con el resultado de 3 iteraciones del método de la
bisección y falsa posición, tomando como intervalo inicial
[0, 1].
2 Implementar un algoritmo en Matlab donde determine la raı́z
de la función f (x) = e −x − x por medio del método de
Newton-Raphson. Pare el cálculo cuando el error sea menor a
0,05 %. Compare los resultados con los dados por el método
de la bisección y falsa posición.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de Newton-Rhapson
Ejemplo
1 .
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Abiertos
Método de la Secante
En el método de Newton Raphson, cambiamos el término f 0 (xi ) por
su término aproximado
f (xi ) − f (xi−1 )
f 0 (xi ) ≈
xi − xi−1
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Métodos Abiertos
Método de la Secante
Reemplazando en la formula tenemos
f (xi )(xi − xi−1 )
xi+1 = xi − .
f (xi ) − f (xi−1 )
Observemos que se necesitan dos puntos anteriores (xi y xi−1 ) para
poder calcular el punto actual xi+1 . Se necesitan dos puntos iniciales
para ejecutar el algoritmo.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos
Métodos de solución
Errores
Métodos Cerrados
Raı́ces de Ecuaciones
Métodos Abiertos
Método de la Secante
Ejemplo
Implementar un algoritmo en Matlab donde determine la raı́z de la
función f (x) = e −x − x por medio del método de Secante. Pare el
cálculo cuando el error sea menor a 0,05 %. Compare los resultados
con los dados por el método de la bisección, falsa posición y Newton-
Raphson.
Usabana
Omar Andrés Cuervo Fernández Métodos Numéricos