0% encontró este documento útil (0 votos)
161 vistas38 páginas

Errores y Raíces en Métodos Numéricos

Este documento trata sobre los errores y las raíces de ecuaciones. Explica que los métodos numéricos proveen aproximaciones a las soluciones exactas y que todo método numérico introduce un error. Describe dos tipos de errores: errores de redondeo debido a la representación finita de los números en un computador, y errores de truncamiento introducidos por los métodos de aproximación numérica. También presenta fórmulas para calcular diferentes tipos de errores y un ejemplo numérico del problema del paracaidista.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
161 vistas38 páginas

Errores y Raíces en Métodos Numéricos

Este documento trata sobre los errores y las raíces de ecuaciones. Explica que los métodos numéricos proveen aproximaciones a las soluciones exactas y que todo método numérico introduce un error. Describe dos tipos de errores: errores de redondeo debido a la representación finita de los números en un computador, y errores de truncamiento introducidos por los métodos de aproximación numérica. También presenta fórmulas para calcular diferentes tipos de errores y un ejemplo numérico del problema del paracaidista.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

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

También podría gustarte