Tema 2- Parte 1: Aproximaciones numéricas

Métodos computacionales/
Métodos Numéricos

Bibliografía:
Métodos Numéricos para Ingenieros.- Chapra y Canale. Ed. McGraw Hill
Interamericana. 2007.
Análisis Numérico. Burden y Faires. 7ma. Edición-.- Cencage Learning
Editores. 2009.-

INTRODUCCIÓN (0)
• La aritmética realizada por una computadora es diferente
de la que se utiliza en el Algebra y en el Calculo.
• Matemática tradicional: números con una cantidad
infinita de cifras.
• Por ej: √ 3 : Como el nro. positivo único que cuando se
multiplica por si mismo produce el entero 3.
• Computación digital: la representación de todo nro.
tiene un nro. finito, fijo de cifras.
• Como no se representa con un nro. finito de cifras -> se
proporciona una evaluación aproximada dentro de la
maquina, una cuyo cuadrado no es exactamente 3, pero
será suficientemente cercana para que sea aceptable en la
mayoría de los casos.
INTRODUCCIÓN (1)
• Los números utilizados en los distintos algoritmos son números reales
o complejos, y son concebidos como fracciones decimales infinitas.
• Con fines computacionales: Deben aproximarse mediante otra forma
de números, conocidos como FRACCIONES TERMINALES
FINITAS, que tienen un número finito de cifras decimales.
• Se introducen entonces situaciones :
• 1) Al reemplazar un número real por un número racional; o sea por
una fracción terminal finita.
• 2) + diferencia entre el sistema real y su modelo matemático.
• 3) + imposibilidad de procesamiento infinito ->
• Estas causas entre otras:
• Producen diferencias entre los resultados verdaderos,
obtenidos del sistema real y aquellos derivados del cálculo,
mediante la aplicación de algún método numérico o
computacional sobre un modelo matemático determinado.

INTRODUCCIÓN (2)

• La eficiencia en el cálculo de la solución numérica

— Depende de:
** la facilidad de implementación del algoritmo
** de las características especiales y limitaciones de los
instrumentos de cálculo (los computadores).
** Solo un subconjunto relativamente pequeño, del sistema de
los nros. reales se usa para representar todos los nros reales
** Este contiene solo números racionales, positivos y negativos.
• En gral, al emplear estos instrumentos de cálculo se introducen
errores.
• ERRORES: diferencias inevitables en toda aplicación
numérica y que, se propagan a lo largo del resto del
proceso de cálculo.

• Objetivo: Estudio de los errores y su influencia a lo

largo del procesamiento.
Números exactos y aproximados
• Clasificación al efecto del calculo numérico o
aproximado:
• CONSTANTES
— Absolutas
— Relativas
• VARIABLES
• NUMEROS EXACTOS:
— No sufren modificación POR CAUSAS OPERATIVAS al ser
utilizados como parámetros o variables de algún modelo
matemático.
— Ej. Enteros cortos ( 2 bytes o 16 bits) ; enteros largos y reales
de simple precisión (4 bytes o 32 BITS) ; reales doble
precisión ( 64 BITS)
— Rango de variación de enteros cortos: ( -2 15; 2 15 - 1 )
— Rango de variación de enteros largos: ( -2 31; 2 31 - 1 )

Números aproximados (I)

• Causas operativas:
• Producidas por el hardware en los casos en que le
resulta imposible soportar al número en su totalidad.
• P.ej.: un nro que necesita mas de 64 bits para ser
almacenado: 1/3; 157 99 ; e, Pi; etc.
• Estos números poseen infinitas cifras decimales o al
menos son muy grandes.
• La máquina es quien intrínsecamente produce el error.
• Causa de errores en computadora: Diferencia que
inevitablemente existe entre un numero a representar
y su real representación en la computadora.
Números aproximados (II)
• La mayor parte de los números que se utilizan en
una computadora no son exactos.
• La representación de los mismos no es continua
• Considerar la expresión:
• Ej: L = 2 Pi R (*)
• (*) tres números con características diferentes
• 2: nro exacto
• PI: computacionalmente solo puede hacerse uso de
un cdad. limitada de dígitos.
• R: proviene de una medición que depende de la
exactitud del instrumento.

Números aproximados (II)

L= 2x 3,27 X 3 =19,62

L= 2x 3,27 X 3,1 =20,274

L= 2x 3,27 X 3,14 =20,5356 ͌

20,54
L= 2x 3,27 X 3,142 =20,548 68

L= 2x 3,27 X 3,1416 =20,546 064

L= 2x 3,27 X 3, 14159 =20,545 998 6

El mejor resultado es el obtenido cuando las aproximaciones de los elementos que intervienen
en una expresión, son homogéneas.
 Errores en los
Cálculos Científicos

Clasificación de los errores

• Definición: “Se denomina con el término genérico de
error, a la diferencia que existe entre el valor
verdadero de una magnitud determinada y otro valor
aproximado de ella.”
— E=Xv–X c

• Causas de ERROR:
• Aproximación Matemática forzada a la realidad física:
— ABSTRACCIÓN
— INHERENCIA
• Imposibilidad de realizar el cálculo en forma exacta
— TRUNCAMIENTO ( Modificación de la solución respecto de su
formulación)
— REDONDEO ( Falta de exactitud en las operaciones aritméticas
elementales.)
Clasificación de los errores
• La resolución de un problema pasa por etapas
• Todas ellas aportan al error
• ABSTRACCIÓN
— Diferencia entre el fenómeno real o verdadero (
sistema) y su descripción analítica simplificada (
modelo)
— Son introducidos por una única vez, al principio del
procesamiento.

Clasificación de los errores

— INHERENCIA:
— Datos introducidos desde el comienzo del procesamiento
del modelo.
— Las diferencias existentes entre los datos de entrada o
parámetros, respecto de sus verdaderos valores
— Generalmente son desconocidos, puede determinarse un
máximo valor estimativo denominado COTA DE ERROR.
— Son inevitables desde todo punto de vista y se introducen
por única vez al principio del procesamiento del modelo
matemático.
Clasificación de los errores
— TRUNCAMIENTO:
— Se produce a lo largo del procesamiento.
— Se produce por la diferencia causada por despreciar
en el desarrollo de una serie infinita, los términos de
orden n+1 en adelante.
— Se disminuye incrementando el nro. de iteraciones.
— La magnitud del error depende del tamaño que se
establezca para el incremento (h o ∆x ).
— Se expresan los errores en función de los incrementos
: E ≈ ( h n ).
— El error máximo cometido (o cota del error) es del
orden de h n , -> E nunca será > que el valor
absoluto de h n .
Clasificación de los errores
— REDONDEO:
— Proviene del modo en que los números son tratados.
— Manual: Cálculos realizados con números racionales,
expresados en notación decimal.
— Por computadora: Se valen de la notación científica.
— Se soporta un número determinado de cifras
significativas
— Producido por la limitación de los números a una
cierta cantidad de dígitos significativos.
— Sólo puede ser minimizado mediante el uso de una
mayor cantidad de dígitos decimales en los cálculos.

Clasificación de los errores

• Los resultados de muchas operaciones aritméticas tienen
más cifras de las que se puede almacenar y hay que
aproximarlos eliminando las cifras menos significativas. A
este proceso se llama "redondeo“.
• Los errores de redondeo son inevitables, pero
controlables
— En muchas ocasiones son poco significativos y no tienen ninguna
importancia
— Sin embargo, en algunos problemas pueden llegar a destruir por
completo el significado de un resultado. Conviene detectar estos
casos y tomar las medidas adecuadas
— Unos errores de redondeo catastróficos pueden ser consecuencia
de un problema difícil, de un mal algoritmo, o de ambas cosas a
la vez.
 Errores de redondeo y Aritmetica de
computadoras (0)

• El uso de dígitos binarios tiende a encubrir las

dificultades de cálculo que ocurren cuando se usa un
conjunto finito de números de máquina para
representar a todos los números reales.
• Para explicar los problemas que pueden surgir, se
considera que los números de máquina se
representan en la forma de punto flotante
decimal normalizada.
• En una computadora no se pueden poner infinitos
dígitos. Se trabaja solo con números de desarrollo
finito y de una longitud dada.

Errores de redondeo y Aritmética de

computadoras (I)

En consecuencia los números de máquina

decimales con k dígitos serán de la forma:

0 ≤ di ≤ 9 Para cada i = 2,. . ., k.

M1 ≤ n ≤ M2 ; el exp. n (orden del nro.) estará
limitado a cierto rango.

Los números k, M1 y M2 dependen de la maquina.

 Errores de redondeo y Aritmetica de
computadoras (II)
• Cualquier número real positivo y puede ser normalizado a

Si y está dentro del rango numérico de la máquina, la

forma de punto flotante de y, que se representará
por fl (y), se obtiene terminando la mantisa de y en
k cifras decimales.

Existen dos formas de llevar a cabo tal terminación.

Errores de redondeo y Aritmetica de

computadoras (III)

Método 1) Un método es simplemente cortar

los dígitos dk+1, dk+2. . .. y se obtiene:

Este método es bastante preciso y se llama truncar el

número.
 Errores de redondeo y Aritmetica de
computadoras (IV)
• Método 2) Agregar 5x10 a y y después
n-(k+1)

truncar para que resulte un número de la

forma:

• Este último método comúnmente se designa por

redondeo del número.
• Si dk+1 ≥ 5, se agrega 1 a dk para obtener fl (y);
esto es, redondeamos hacia arriba.
• Si d k+1 < 5, simplemente se trunca luego de los
primeros k dígitos; se redondea así hacia abajo.

Errores de redondeo y Aritmetica de

computadoras (V)
• EJEMPLO 1: El número pi (π) tiene un desarrollo decimal
infinito de la forma π = 3.14159265. . .
• Escrito en forma decimal normalizada, se tiene:

π = 0.314159265...........x101

Método 1) La forma de punto flotante de cinco

dígitos de π utilizando truncamiento es:
fl ( )  0.31415 *101  3.1415
Método 2) Dado que el sexto digito de la expansión
decimal de π es 9, la forma de π con redondeo a cinco
dígitos es:
fl ( )  (0.31415  0.00001) *101  3.1416 Dk+1 ≥ 5
 Errores de redondeo y Aritmetica de
computadoras (VI)
• El error que resulta al reemplazar un número, por su forma
de punto flotante, se llama error de redondeo (sin que
importe si se usa el método de redondeo o de
truncamiento).
• Como medir los errores de aproximación: Definiciones:
Si p* es una aproximación a p,
El error absoluto (Ea) es: |p - p*| . Desventaja?

El error relativo (Er) es: p  p * siempre que p ≠ 0.

p
Este error permite normalizar el error respecto al
valor verdadero.
p  p * * 100 % = E p
Error porcentual =
p

Errores de redondeo y Aritmetica de

computadoras (VII)
• Considere los errores absoluto y relativo al representar p
por p* en el siguiente ejemplo.
p p*
0.3000 x 101 0.3100 x 101

0.3000x 10-3 0.3100 x 10-3

0.3000x 10 4 0.3100 x 104

Analice los errores obtenidos!!

 Errores de redondeo y Aritmetica de
computadoras (VII)
• Considere los errores absoluto y relativo al representar p
por p* en el siguiente ejemplo.
p p* Absoluto Relativo
0.3000 x 101 0.3100 x 101

0.3000x 10-3 0.3100 x 10-3

0.3000x 10 4 0.3100 x 104

Errores de redondeo y Aritmetica de

computadoras (VII)
• Considere los errores absoluto y relativo al representar p
por p* en el siguiente ejemplo.
p p* Absoluto Relativo
0.3000 x 101 0.3100 x 101 0.1 0.3333 x 10-1

0.3000x 10-3 0.3100 x 10-3 0.1 x 10-4 0.3333 x 10-1

0.3000x 10 4 0.3100 x 104 0.1 x 103 0.3333 x 10-1.

El error relativo es una medida de mayor significación.

El error absoluto puede ser puesto en 2do. termino
 Errores de redondeo y Aritmetica de
computadoras (VIII)

Representación de los números en la máquina:

la de punto flotante fl(y) de un número y tiene el error relativo

y  fl ( y )
y
Si se emplean k cifras decimales y el truncamiento para la
representación en la máquina de:

y  0.d1d 2 ......d k d k 1..... *10n

Entonces:

Errores de redondeo y Aritmetica de

computadoras (VIII)
Representación de los números en la máquina:
Se observa que fl(y) como punto flotante para el número y tiene
el error relativo:
y  fl ( y )
y
Si se emplean k cifras decimales y el truncamiento para la
representación en la máquina de:

y  0.d1d 2 ......d k d k 1..... *10n

Entonces:

y  fl ( y ) 0.d1d 2 ...d k d k 1... *10 n  0.d1d 2 ...d k *10 n


y 0.d1d 2 ...d k *10 n
 Errores de redondeo y Aritmetica de
computadoras (IX)

y  fl ( y ) 1 Cota de error relativo

 *10  k  10  k 1 por truncamiento
y 0 .1

De manera similar, una Cota para el error relativo

cuando se usa aritmética con redondeo a k dígitos
es
0.5 x 10-k+1.

Errores de redondeo y Aritmetica de

computadoras (IX)

Dado que d1 ≠ O, el valor mínimo del denominador

es 0.1. El numerador está acotado superiormente
por 1. Entonces ,

y  fl ( y ) 1 Cota de error relativo

 *10  k  10  k 1 por truncamiento
y 0 .1

De manera similar, una Cota para el error relativo

cuando se usa aritmética con redondeo a k dígitos
es 0.5 x 10-k+1.
 Errores de redondeo y Aritmetica de
computadoras (X)
• Importante: las cotas para el error relativo, cuando se usa
aritmética de k dígitos, son independientes del número que
se representa. Esto se debe a la forma en que los números
de máquina están distribuidos a lo largo de la recta real.
• Debido a la forma exponencial de la característica, la misma
cantidad de números de máquina decimales se emplean para
representar cada uno de los intervalos [0.1, 1], [1, 10] y [10,
100].
• Dentro de los límites de la máquina, la cdad de números
decimales de máquina en [10n , 10n+1] es constante para
todos los enteros n.

Errores de redondeo y Aritmética de computadoras

(XI)
• Además de la representación imprecisa de números, la
aritmética realizada en una computadora no es exacta.
• Las operaciones aritméticas generalmente implican
manipular dígitos binarios mediante diversos corrimientos u
operaciones lógicas.
• Dado que la mecánica real de estas operaciones no tiene que
ver con esa representación, es necesario utilizar una
aproximación apropiada a la aritmética de computadora.
• La pérdida de precisión debida al error de redondeo puede
con frecuencia ser evitada por una cuidadosa serie de
operaciones o por una reformulación del problema,
• Se plantea a continuación:
 Errores de redondeo y Aritmetica de
computadoras (XII)
• Considere que la representación de punto flotante fl(x) y
fl(y) esté dada para los números reales x e y, y que los
símbolos +,-,x,/, representan las operaciones de adición,
sustracción, multiplicación y división de máquina o en la
computadora, respectivamente.
• Supondremos una aritmética con un nro. finito de cifras
dadas por
x+ y = fl (fl(x) + fl( y)), x * y = fl (fl(x) * fl( y)),
x - y = fl (fl(x) - fl( y)), x / y = fl(fl(x) / fl( y))

Esto corresponde a realizar aritmética exacta con las

representaciones de punto flotante de x e y ;
para luego convertir el resultado exacto a su
representación de pto. flotante con un nro. finito
de cifras.

Errores de redondeo y Aritmetica de

computadoras (XII)
• Considere que la representación de punto flotante fl(x) y
fl(y) esté dada para los números reales x y y, y que los
símbolos +,-,x,/, representan las operaciones de adición,
sustracción, multiplicación y división de máquina o en la
computadora, respectivamente.
• Supondremos una aritmética de dígitos finitos dada por
x+ y = fl (fl(x) + fl( y)), x * y = fl (fl(x) * fl( y)),
x - y = fl (fl(x) - fl( y)), x / y = fl(fl(x) / fl( y))

Esta aritmética corresponde a realizar aritmética

exacta en las representaciones de punto flotante
de x e y ; y luego convertir el resultado exacto a
su representación de pto. flotante para dígitos
finitos.
 Errores de redondeo y Aritmetica de
computadoras (XIII)
• EJEMPLO: Dado que x = 1/3, y = 5/7, y que se utiliza
truncamiento de cinco cifras para los cálculos aritméticos
donde intervienen x e y. La Tabla da los valores de las
operaciones en computadora con

fl ( x)  0.33333 *100 fl ( y )  0.71428 *100

Operación Resultado Valor Real

x+y 0,10476*101 22/21

y-x 0,38095*100 8/21

x*y 0,23809*100 5/21

y/x 0,21428*101 15/7

Errores de redondeo y Aritmetica de

computadoras (XIV)

fl ( x)  0.33333 *100 fl ( y )  0.71428 *10 0

Opera Error Error
ción Resultado Valor Real Absoluto Relativo
x+y 0,10476*101 22/21

y-x 0,38095*100 8/21

x*y 0,23809*100 5/21

y/x 0,21428*101 15/7

 Errores de redondeo y Aritmetica de
computadoras (XIV)

fl ( x)  0.33333 *100 fl ( y )  0.71428 *10 0

Opera Error Error
ción Resultado Valor Real Absoluto Relativo
x+y 0,10476*101 22/21 0,190*10-4 0,182*10-4

y-x 0,38095*100 8/21 0,238*10-5 0,625*10-5

x*y 0,23809*100 5/21 0,524*10-5 0,220*10-4

y/x 0,21428*101 15/7 0,571*10-4 0,267*10-4

• Nótese que el máximo error relativo para las operaciones en

el ejemplo es 0.267 x 10-4, -> la aritmética produce
resultados satisfactorios a cinco dígitos. ( cifras significativas
1) ( cifras significativas 2 )

Errores de redondeo y Aritmetica de

computadoras (XV)

Si se tiene. u= 0.714251, v = 98765.9, y w = 0.11111 X 10-4

de tal forma que

fl (u )  0.71415 *10 0 Se muestran algunos

problemas que pueden
fl (v)  0.98765 *105
surgir con la aritmética
fl ( w)  0.11111*10  4 cuando se tiene una
cantidad finita de cifras.
 Errores de redondeo y Aritmética de
computadoras (III)

Error Error
Operac. Resultado Valor real absoluto relativo

y-u 0.30000 X 10-4 0.34714 X 10-4

(y-u)/w 0.27000 X l01 0.31243 X l01

(y-u)*v 0.29629 X l01 0.34285 X l01

u+v 0.98765 X l05 0.98766 X l05

Errores de redondeo y Aritmetica de

computadoras (III)
Error Error
Operac. Resultado Valor real absoluto relativo

y-u 0.30000 X 10-4 0.34714 X 10-4 0.471 X l0-5 0.136

(y-u)/w 0.27000 X l01 0.31243 X l01 0.424 0.136
(y-u)*v 0.29629 X l01 0.34285 X l01 0.465 0.136

u+v 0.98765 X l05 0.98766 X l05 0.161 X l01 0.163 X 10-4

• En particular, u+v nos dice que si se tiene que sumar

varios nros x1; x2; :::::xN conviene hacerlo de menor a
mayor (¿Por que?).
• La pérdida de precisión debida al error de redondeo puede
ser evitada con frecuencia por una cuidadosa serie de
operaciones o por una reformulación del problema.
 Cancelación de cifras significativas (I)
• Uno de los cálculos mas comunes que producen errores es
debido a la resta de números casi iguales.
• Suponga que dos números casi iguales x e y, con x > y, con
representaciones de k cifras
fl ( x )  0.d 1 d 2 ....... d p  p 1 p  2 ...... k x 10 n

y
fl ( y )  0.d1d 2 .......d p  p 1 p  2 ...... k x 10 n ,
La forma de punto flotante de x-y es:
fl ( fl ( x)  fl ( y ))  0. p 1 p  2 ...... k x 10 n -p ,
donde
0. p 1 p  2 ...... k  0. p 1 p  2 ...... k  0. p 1  p  2 ...... k

Cancelación de cifras significativas (II)

• El nro. de pto. flotante utilizado para representar la diferencia
x- y, tiene a lo sumo: k- p cifras significativas
• La mayor parte de las computadoras a x-y le asignarán k
cifras ->> las últimas p se anularán o serán asignadas al azar.
• Así en los cálculos posteriores con x-y se tendrá el problema
de contar con k-p cifras significativas
• Si una representación con un número finito de cifras o un
calculo introduce error –> este aumenta al dividir entre un
nro. con magnitud pequeña
• Multiplicación: el error aumenta en forma equivalente al
multiplicar dicho error por un nro. con magnitud grande.
 Cancelación de cifras significativas (III)
Suponga que el nro. z tiene la aproximación con un nro. finito de
cifras z+δ, donde el error δ surge por la representación o por
un cálculo anterior.
z  fl ( z ) 
 fl    ( z   ) x10 n
  fl ( ) 
Así el error absoluto en esta aproximación,  x10
n

Suponga que dividimos entre ε = 10 –n, donde n > 0 .

Entonces:

Es el error absoluto original , δ multiplicado por el factor 10 n

Cancelación de cifras significativas (IV)

Sean p= 0.54617 y q = 0.54601. El valor exacto de la resta
R= p – q es = 0.00016.
1) Si la resta se realiza con aritmética a 4 cifras.
1) Redondear p y q. Hallar el error relativo de la resta.
2) Truncar, realizar la resta y obtener el error relativo.
Analizar si se pierde precisión en dichos cálculos.
 Cancelación de cifras significativas (IV)
Sean p= 0.54617 y q = 0.54601. El valor exacto de la resta
R= p – q es = 0.00016.
Valores redondeados
P = 0.5462 P = 0.54617 + 0.00005 = 0.54622 a 4
dígitos es 0.5462

Q = 0.5460 Q = 0.5460 + 0.00005 = 0.54605 a 4

dígitos es 0.5460
R = 0.5462- 0.5460 = 0.0002 (con valores
redondeados )

Er = 0.00016 – 0.0002 = -0.25

0.00016

Cancelación de cifras significativas (IV)

• b) Trabajar con aritmética con 4 dígitos, truncar P y Q,
restar y hallar el Error relativo de dicha recta
• Valores truncados
• P = 0.5461 Q= =.5460
• R = 0.5461- 0.5460 = 0.0001 (Con valores truncados)
•
• Er = 0.00016 – 0.0001 = 0.3750
• 0.00016
•
• El error relativo en valor absoluto es menor cuando se
redondean los valores que intervienen en la resta que
cuando se los trunca, por lo que se concluye que con
valores truncados se pierde precisión en los resultados.
Consideraciones aritméticas: Problemas Mal
Planteados (1)
• Diferencia de números parecidos o cancelación por
resta
—La resta entre números de una magnitud parecida
puede hacer perder varias cifras significativas La
solución es tratar de realizar las operaciones de otra
forma.
—Un caso común donde esto ocurre es en la
determinación de las raíces cuadráticas o parábola
usando la formula cuadrática:

 b  b 2  4ac
x1 y x 2 
2a

Consideraciones aritméticas: Problemas Mal

Planteados (2)

—En los casos donde b2 >= 4ac, la diferencia en el

numerador puede ser muy pequeña. En tales casos
la precisión doble reduce el problema. Además se
puede usar como alternativa la formula:

 2c
x1 y x 2 
 b  b 2  4ac

— Para estos casos se cambia la fórmula cuadrática

mediante la racionalización del numerador.
 Consideraciones aritméticas: Problemas Mal
Planteados (3)
• Orden de las operaciones de acumulación
— Si se suman primero los términos más pequeños se pierden
menos cifras significativas que si se empieza sumando los
términos de mayor valor
• Comparaciones
— Signo de números pequeños: Para saber si un número pequeño
es mayor que cero conviene establecer un valor límite (1e-12,
por ejemplo) pues los números próximos a cero pueden tener
un signo u otro según los errores de redondeo
— Comparación de números de punto flotante : Nunca se deben
comparar con el operador == directamente. Hay que ver si el
valor absoluto de su diferencia dividido por el valor absoluto del
número es menor que un determinado número pequeño (1e-12,
por ejemplo)