Métodos Numéricos, Sesión 1: Teoría del Error.
Contenidos.
a) Introducción.
b) Error de truncamiento.
c) Error de propagación.
d) Error absoluto.
e) Error relativo.
1. Errores
1.1. Introducción
La aritmética realizada con una calculadora o computadora es diferente a la Aritmética que
se imparte en los cursos de álgebra y cálculo. En el mundo matemático, la aritmética en R, es de
esperarse:
2+5 = 7
√ 3
3
2 = 2
sin embargo, con la aritmética computacional, esperamos resultados exactos para
2·2 = 4
23 = 8
√ 2
pero no obtendremos exactamente 3 = 3, en el mundo computacional, cada número represen-
table sólo tiene un número jo y nito de dígitos. Esto signi ca que, por ejemplo, sólo los números
racionales, e incluso no todos ellos, se pueden representar de forma exacta. Entonces, en muchos
caso, esta aritmética mecánica es satisfactoria y pasa sin importancia o preocupación, pero algunas
veces surgen problemas debido a su discrepancia.
1.2. Errores de truncamiento, de propagación
El error que se produce cuando se utiliza una calculadora o computadora para realizar cálculos
con números reales recibe el nombre de error de redondeo
Métodos Numéricos 1
1.3. Punto otante
La forma de trabajar el números de máquina, se denomina punto otante decimal normalizada:
±0, d1 d2 d3 . . . dk × 10n = y
1 ≤ d1 ≤ 9
0 ≤ di ≤ 9
La forma de punto otante y , denominada f l(y), se obtiene terminando la mantisa de y en k cifras
decimales. Tenemos dos formar de realizarlo:
1.3.1. Trucamiento o de corte
Consiste en cortar los dígitos dk+1 dk+2 . . .:
f l(y) = 0, d1 d2 d3 . . . dk × 10n
1.3.2. Redondeo
Suma 5 × 10n−(k+1) a y , luego trunca el resultado para obtener:
f l(y) = 0, a1 a2 a3 . . . ak × 10n
Por redondeo, si dk+1 ≥ 5, sumamos 1 a dk , si dk+1 < 5 truncamos todo excepto los primeros k
dígitos.
Ejemplo 1 La expansión decimal de π es in nita de la forma π = 3,14159265 . . ., escrito en forma
normalizada:
π = 0,314159265 . . . × 101
a) El formato de punto otante de π usando el recorte de cinco digitos es
f l(π) = 0,31415 × 101 = 3,1415
b) El sexto digito de la expansión decimal de π es un 9 , por lo que el formato de punto otante
de π con redondeo de cinco digitos es
f l(π) = (0,31415 + 0,00001) × 101 = 3,1416
Métodos Numéricos 2
1.4. Error Absoluto, Error Relativo
La siguiente de nición describe tres métodos para medir errores de aproximación. En general, el
error relativo es una mejor medición de precisión que el error absoluto porque considera el tamaño
del número que se va a aproximar.
De nición 1 Suponga que p∗ es una
∗
aproximación a p. El error real es p − p∗ , el error absoluto
es |p − p∗ | , y el error relativo es |p−p |
|p|
, siempre y cuando p 6= 0
1.1 Determine los errores real, absoluto y relativo al aproximar p con p∗ cuando:
a) p = 0,3000 × 101 y p∗ = 0,3100 × 101
b) p = 0,3000 × 10−3 y p∗ = 0,3100 × 10−3
c) p = 0,3000 × 104 y p∗ = 0,3100 × 104
De nición 2 Se dice que el número p∗ se aproxima a p para t dígitos signi cativos (o cifras) si t
es el entero no negativo más grande para el que
|p − p∗ |
≤ 5 × 10−t
|p|
A menudo, el término digitos signi cativos se usa para describir vagamente el número de digitos
decimales que parecen ser exactos. La de nición es más precisa y provee un concepto continuo.
1.5. Aritmética de dígitos nitos
Además de la representación inexacta de números, la aritmética que se efectúa en una compu-
tadora no es exacta. La aritmética implica manipular dígitos binarios mediante diferentes opera-
ciones de cambio, o lógicas. Puesto que la mecánica real de estas operaciones no es pertinente para
esta presentación, debemos idear una aproximación propia para aritmética computacional. A pe-
sar de que nuestra aritmética no proporcionará el panorama exacto, es su ciente para explicar los
problemas que se presentan. Piense que las representaciones de punto otante f (x) y f l(y) están
dadas para los números reales x y y y que los símbolos ⊕, , ⊗, y ⊕ representan operaciones de
máquina de suma, resta, multiplicación y división, respectivamente. Supondremos una aritmética
de dígitos nitos provista por
x ⊕ y = f l(f l(x) + f l(y)), x ⊗ y = f l(f l(x) × f l(y))
x y = f l(f l(x) − f l(y)), x ⊕ y = f l(f l(x) ÷ f l(y))
Esta aritmética corresponde a realizar aritmética exacta en las representaciones de punto otante
de x y y y, después, convertir el resultado exacto a su representación de punto otante de dígitos
nitos.
Métodos Numéricos 3
2. Ejercicios
2.1 Suponga que x = 75 y y = 13 . Utilice el corte de cinco digitos para calcular x + y, x − y, x × y
y x÷y
2.2 La fórmula cuadrática establece que las raíces de ax2 + bx + c = 0, cuando a 6= 0, son
√ √
−b + b2 − 4ac −b − b2 − 4ac
x1 = y x2 =
2a 2a
Considere esta fórmula aplicada a la ecuación x2 +62,10x+1 = 0, cuyas raices son aproximadamente
x1 = −0,01610723 y x2 = −62,08390
Utilice aritmética de cuatro dígitos para determinar las raíces.
2.3 Considere f (x) = x3 − 6,1x2 + 3,2x + 1,5 en x = 4,71 con aritmética de tres dígitos. El
resultado exacto de la evaluación es
Exacto: f (4,71) = 104,487111 − 135,32301 + 15,072 + 1,5 = −14,263899
a) Aplicar truncamiento.
b) Aplicar redondeo.
c) Medir los errores relativos obtenidos.
d) Reescribir de forma anidada;
f (x) = x3 − 6,1x2 + 3,2x + 1,5 = ((x − 6,1)x + 3,2)x + 1,5
y realizar nuevamente los cálculos anteriores. Analizar sus resultados.
3. Actividades para el estudiante.
3.1 Utilice redondeo aritmético a tres cifras para el cálculo de:
3
−10π + 6e −
62
Calcule el error absoluto y el error relativo.
Métodos Numéricos 4
3.2 Utilice redondeo aritmético a tres cifras para los siguientes cálculos. Calcule además error
absoluto y relativo con el valor exacto a por lo menos cinco cifras.
a. 133 + 0,921
b. 133 − 0,499
c. (121 − 0,327) − 119
d. (121 − 119) − 0,327
13
− 67
e. 2e−5,4
14
f. −10π 3
+ 6e − 62
g. 29 · 97
π− 22
h. 1 7
17
3.3 Los tres primeros términos no nulos de la serie de McLaurin para la función arco tangente
son:
1 1
x − x3 + x5
3 5
Calcule el error absoluto y el error relativo en la siguiente aproximación de π usando el polinomio
en lugar de la función arco tangente:
1 1
4 arctan + arctan
2 3
3.4 Considere:
1 3 3
− +
3 11 20
Calcular:
a) En forma exacta.
b) Truncamiento a tres cifras.
c) Con redondeo a tres cifras.
d) Calcule los errores relativos.
3.5 Use redondeo a 4 cifras para determinar las aproximaciones de la ecuación cuadrática:
1 2 123 1
x − x+ =0
3 4 6
Calcule los errores relativos.
Métodos Numéricos 5
3.6 Suponga que dos puntos (x0 , y0 ) y (x1 , y1 ) están en una línea recta, con y1 6= y0 . Se tienen
dos fórmulas para determinar la abscisa al origen de la recta:
x0 y1 − x1 y0
x =
y1 − y0
(x1 − x0 )y0
x = x0 −
y1 − y0
a) Pruebe que ambas son correctas algebraícamente.
b) Utilice los datos (x0 , y0 ) = (1,31, 3,24) y (1,93, 4,76) y el redondeo artimético a tres cifras
para calcular al abcisa al origen de ambas formas. ¾Cual es el mejor y por qué?
Métodos Numéricos 6
4. Respuestas
1.1 a) Para p = 0,3000 × 101 y p∗ = 0,3100 × 101 , el error real es −0,1, el error absoluto es
0,1y el error relativo es 0,3333 × 10−1
b) Para p = 0,3000 × 10−3 y p∗ = 0,3100 × 10−3 , el error real es −0,1 × 10−4 , el error absoluto
es 0,1 × 10−4 y el error relativo es 0,3333 × 10−1
c) Para p = 0,3000 × 104 y p∗ = 0,3100 × 104, el error real es −0,1 × 103 , el error ab soluto
es 0,1 × 103 y, de nuevo, el error relativo es 0,3333 × 10−1
2.1 Observe que
5
x= = 0.714285
7
1
y= = 0.3
3
implica que los valores de corte de cinco digitos de x y y son
f l(x) = 0,71428 × 100 y f l(y) = 0,33333 × 100
Por lo tanto,
x ⊕ y = f l(f l(x) + f l(y)) = f l 0,71428 × 100 + 0,33333 × 100
= f l 1,04761 × 100 = 0,10476 × 101
El valor verdadero es x + y = 57 + 31 = 22
21
, por lo que tenemos
22
Error absoluto = − 0,10476 × 10 = 0,190 × 10−4
1
21
y
0,190 × 10−4
Error relativo = = 0,182 × 10−4
22/21
2.2 Usaremos otra vez la aritmética de redondeo de cuatro digitos en los cálculos para determinar
la raíz. En esta ecuación, b2 es mucho más grande que 4ac, por lo que el numerador en el cálculo
para x1 implica la resta de números casi iguales. Ya que
√ p
b2 − 4ac = (62,10)2 − (4,000)(1,000)(1,000)
p √
= 3856. − 4,000 = 3852. = 62,06
Métodos Numéricos 7
tenemos
−62,10 + 62,06 −0,04000
f l (x1 ) = = = −0,02000
2,000 2,000
una aproximación de ciente a x1 = −0,01611, con un error relativo grande
| − 0,01611 + 0,02000|
≈ 2,4 × 10−1
| − 0,01611|
√
Por otro lado, el cálculo para x2 implica la suma de los números casi iguales −b y − b2 − 4ac.
Esto no presenta problemas debido a que
−62,10 − 62,06 −124,2
f l (x2 ) = = = −62,10
2,000 2,000
tiene un error relativo pequeño
| − 62,08 + 62,10|
≈ 3,2 × 10−4
| − 62,08|
2.3 a) Tres dígitos (corte): f (4,71) = ((104. − 134.) + 15,0) + 1,5 = −13,5.
b) Tres dígitos (redondeo): f (4,71) = ((105. − 135.) + 15,1) + 1,5 = −13,4.
c) Corte: −14,263899+13,5 ≈ 0,05, y redondeo: −14,263899+13,4 ≈ 0,06
−14,263899 −14,263899
d) Ahora, la aritmética de corte de tres dígitos produce
f (4,71) = ((4,71 − 6,1)4,71 + 3,2)4,71 + 1,5 = ((−1,39)(4,71) + 3,2)4,71 + 1,5
= (−6,54 + 3,2)4,71 + 1,5 = (−3,34)4,71 + 1,5 = −15,7 + 1,5 = −14,2
De manera similar, ahora obtenemos una respuesta de redondeo de tres digitos de −14,3.
Los nuevos errores relativos son Tres dígitos (corte): −14,263899+14,2 ≈ 0,0045, Tres dígitos
−14,263899
(redondeo): −14,263899 ≈ 0,0025
−14,263899+14,3
El anidado ha reducido el error relativo para la aproximación de corte a menos de 10 % del
valor obtenido al inicio. Para la aproximación de redondeo, la mejora ha sido todavía más
drástica; el error, en este caso, se ha reducido más de 95 por ciento.
Los polinomios siempre deberian expresarse en forma anidada antes de realizar una evaluación
porque esta forma minimiza el número de cálculos aritméticos. La disminución del error en
la ilustración se debe a la reducción de los cálculos de cuatro multiplicaciones y tres sumas
a dos multiplicaciones y tres sumas. Una forma de disminuir el error de redondeo es reducir
el número de cálculos.
Métodos Numéricos 8
3.2
Aproximacin Error absoluto. Error relativo.
a. 134 0,079 5,90 × 10−4
b. 133 0,499 3,77 × 10−3
c. 2,00 0,327 0,195
d. 1,67 0,003 1,79 × 10−3
e. 1,80 0,154 0,0786
f. −15,1 0,0546 3,60 × 10−3
g. 0,286 2,86 × 10−4 10−3
h. 0,00 0,0215 1,00
3.6 La primera fórmula da −0,00658, y la segunda −0,0100. El valor real a tres dígitos es −0,0116
Métodos Numéricos 9
