Análisis Numérico I

PRACTICA 1a Ejercicio 3. En 1862 el físico Foucault, utilizando un

espejo giratorio, calculó en 298 000 km/s la velocidad
Errores. de la luz. Aceptando como exacta la velocidad de
299 776 km/s, estimar el error absoluto y el error
¿Cuánto es 2 + 2?
relativo cometido por Foucault. Por otra parte, la
Ingeniero: 3.999999. determinación de la constante universal h realizada
Físico: 4.00004 ± 0.00006. por Planck en 1913 dio el valor de 6.41×10−27 erg seg.
Matemático: La solución existe y Adoptando el valor de 6.623×10−27 erg seg, estimar el
es única... falta acotarla. error absoluto y relativo cometido por Planck. ¿Cual
medida es más precisa?

Definiciones básicas. Sea x e un valor aproximado

Número de decimales correctos, dígitos signi-
de una cantidad cuyo valor exacto es x.
ficativos y redondeo. Cuando se indica el número
El error absoluto de x
e se define como de dígitos de un valor numérico no se cuentan los posi-
bles ceros localizados al comienzo del número, puesto
∆=x−x
e. que los mismos solo ayudan a indicar la localización
del punto decimal. Por el contrario, si se indica el
El error relativo de x
e, como número de decimales, entonces debe contarse aque-
llos ceros iniciales que estén a la derecha del punto
∆ x−x decimal. Por ejemplo, el número 0.00147 está dado
(x 6= 0).
e
=
x x con tres dígitos pero tiene cinco decimales, mientras
que el número 13.24 está dado con cuatro dígitos y
∆
Y el error porcentual como 100 . solo dos decimales.
x
Se dice que el número x e aproxima a x con t
decimales correctos si t es el entero no negativo
+ En general no conocemos el valor verdadero y por
más grande para el cual
lo tanto, tampoco el error real. La situación típica
es que tenemos una estimación del error dada por 1
una cota positiva al tamaño máximo del error. Es |x − x
e| < × 10−t .
2
importante remarcar la diferencia entre el error, el
cual puede ser positivo o negativo, y una cota de error Se dice que el número x e aproxima a x con t
positiva  para la magnitud del error absoluto o δ dígitos o cifras significativas si t es el entero no
para la magnitud del error relativo: negativo más grande para el cual
|x − xe| 1
|∆| < × 101−t = 5 × 10−t .
|∆| ≤ , ≤ δ. |x| 2
|x|

En la práctica, sólo tales cotas pueden ser determina- + En virtud de estas definiciones, el número de
das. decimales correctos da una idea de la magnitud del
error absoluto, mientras que el número de dígitos
+ En la práctica suele utilizarse la notación x = xe± significativos da una idea de la magnitud del error
para indicar que |x − xe| ≤ . relativo.

+ Usualmente, el error relativo no puede ser cal- Dado un número x existen dos maneras de redon-
culado directamente puesto que x es una cantidad dear el número a, digamos, t decimales:
desconocida. Sin embargo, en tanto el error relativo
En el truncamiento o cortado simplemente se des-
sea pequeño, se puede estimar su magnitud a partir
cartan todos los decimales ubicados a la derecha
de considerar |∆|/|x| ' |∆|/|ex|.
del t−ésimo decimal para obtener x e. Procediendo
de esta manera el error cometido tiene, sistemá-
Ejercicio 1. Estimar los errores absolutos, relativos ticamente, el mismo signo que el número y su
y porcentuales que se comenten al tomar como valores magnitud está acotada por |x − x e| ≤ 10−t .
de π: En el redondeo simétrico, si el (t + 1)−ésimo de-
cimal es menor que 5 dejamos el t-ésimo decimal
a) 22/7, b) 333/106, c) 355/113. sin cambiar y eliminamos los siguientes, pero si
el (t + 1)−ésimo decimal es mayor o igual que
Ejercicio 2. Calcular el error absoluto que se comete 5 sumamos 1 al t-ésimo decimal. Los errores ob-
cuando se adopta como valor del área de un rectángulo tenidos resultarán positivos o negativos casi por
de lados a y b, el área del cuadrado que tiene como igual, y una cota para el error cometido está dada
lado el promedio de a y b. por |x − xe| ≤ 21 10−t .

Práctica 1a 1
 Análisis Numérico I

Ejercicio 4. Mostrar que 0.001234 ± 0.000004 tiene Problema directo e inverso de errores. La de-
cinco decimales correctos y tres dígitos significati- terminación del error de un valor obtenido como re-
vos, mientras que 0.001234 ± 0.000008 tiene cuatro sultado de operaciones sobre cantidades cuyos errores
decimales correctos y dos dígitos significativos. se conocen, constituye el problema directo del cálculo
de errores. Por otra parte, la estimación de los erro-
Ejercicio 5. Redondear π a cuatro decimales em-
res de los valores que intervienen en una cantidad
pleando el cortado y el redondeo simétrico.
de manera de que ésta sea calculada con un error
prefijado, constituye el problema inverso del cálculo
Propagación de errores. Nos interesa analizar de errores. La solución de estos problemas se apoyan
como se propaga el error al efectuar operaciones arit- en las fórmulas de propagación de errores.
méticas con aproximaciones.
Ejercicio 11. Determinar el valor y el error come-
Ejercicio 6. Mostrar que en la adición y sustracción,
tido en el cálculo del volumen de un cubo de lado
una cota para el error absoluto del resultado está dada
l = 7.3 cm medido con un error menor que medio milí-
por la suma de las cotas de los errores absolutos de
metro. ¿Cuántos decimales correctos tiene el volumen
los operandos
calculado? ¿cuántas cifras significativas?
|∆(x ± y)| ≤ |∆x| + |∆y|,
Ejercicio 12. Dados π y e redondeados a seis de-
mientras que en la multiplicación y división, una cota cimales, ¿cuantos dígitos significativos y decimales
del error relativo está dada, aproximadamente, por correctos tendrá el cociente π/e?
la suma de las cotas de los errores relativos de los
operandos, Ejercicio 13. Supongamos que se han determinado
|∆(x · y)| |∆x| |∆y| los siguientes valores para dos lados y un ángulo de
. + , un triángulo ABC:
|x · y| |x| |y|
|∆(x/y)| |∆x| |∆y| b = 465.24 m, c = 404.63 m, α = 48◦ 250 1000 .
. + .
|x/y| |x| |y|
Los lados se han determinado con errores absolutos
Ejercicio 7. Mostrar que el error relativo en la di- inferiores a 0.05 m y el ángulo con un error absoluto
ferencia de dos números aproximadamente iguales inferior a 1000 . El área del triángulo se calcula según
puede ser muy grande, aún cuando los dos números la fórmula A = 1 b c sin α. Estimar el valor de dicha
2
tengan errores pequeños (fenómeno de cancelación). área y el error cometido.
−4
Considerar como ejemplo x = 0.5764 ± 0.5 × 10 e
y = 0.5763 ± 0.5 × 10−4 . Ejercicio 14. Se desea calcular el volumen de un
paralelepípedo rectángulo, cuyas aristas miden apro-
Consideremos ahora el caso general de una relación ximadamente x = 35 cm, y = 40 cm, z = 45 cm, con
funcional z = z(x1 , x2 , . . . , xn ) para ciertas varia- un error menor que 50 cm3 . ¿Con qué aproximación
bles x1 , x2 , . . . , xn . Una cota para el error absoluto deben medirse las aristas?
∆z = z(x1 , x2 , . . . , xn ) − z(e x1 , x
e2 , . . . , x
en ) al estimar
z con el valor aproximado x ei de xi con error ∆xi , Ejercicio 15. Se desea calcular con una aproxima-
para i = 1, . . . , n, está dada por la siguiente fórmula ción del 0.1 % la superficie de un círculo cuyo radio
fundamental del cálculo de errores: mide aproximadamente 25 cm. ¿Con qué aproxima-
n  ción debe medirse el radio y cuantos decimales de π
 ∂z 

será necesario considerar?
X
|∆z| . |∆xi |.
∂x
 
i x1 ,e
(e x2 ,...,e
xn )
i=1

+ La cantidad |∂z/∂xi | puede ser interpretada co- Error estadístico. En la práctica, la cota para
mo una medida de la sensibilidad de z frente a una |∆z| obtenida por la fórmula fundamental de pro-
perturbación ∆xi en el argumento xi . pagación de errores suele dar una sobreestimación
muy pesimista del error puesto que ésta contempla el
Ejercicio 8. Mostrar que si se consideran las rela- peor caso posible donde los errores ∆xi contribuyen
ciones funcionales z = x + y, z = x − y, z = x · y con el mismo signo y sus magnitudes son iguales sus
y z = x/y la fórmula fundamental conduce a las cotas. Por esta razón, la cota para |∆z| así obtenida
fórmulas de propagación de errores ya enunciadas. se denomina error máximo. Una estimación alternati-
va de la magnitud del error resulta de la aplicación
Ejercicio 9. Mostrar que el error relativo de la raíz
de métodos estadísticos, y tiene su mayor aplicación
cuadrada de un número es inferior a la mitad de una
en los errores de observación que se comenten al de-
cota superior del error relativo del radicando.
terminar el valor númerico de una magnitud física.
Ejercicio 10. Mostrar que una cota para error rela- Los errores de observación se clasifican en errores
tivo en el volumen de una esfera de radio R es tres sistemáticos y errores aleatorios. Un error sistemá-
veces el error relativo con que ha sido medido el radio tico es un error intrínsico al método de observación
R. en sí, por lo que estará siempre presente y tendrá

Práctica 1a 2
 Análisis Numérico I

una magnitud constante en la aplicación de dicho

método. Tales errores no pueden ser modelados en
forma probabilística, pero si pueden ser analizados
con la fórmula fundamental de propagación de erro-
res. Por otra parte, los errores aleatorios se deben a
un gran número de causas irregulares que el obser-
vador no puede controlar por completo, de manera
que fluctuarán de forma irregular de una observación
a otra. Tal clase de errores pueden ser tratados con
métodos estadísticos. Especificamente, asumiendo que
los errores aleatorios ∆xi pueden considerarse como
variables aleatorias independientes con valor medio
cero y desviaciones estándar 1 , 2 , . . . , n , entonces el
error estadístico para z = z(x1 , x2 , . . . , xn ) está dado
por la fórmula:

n 
!1/2
X ∂z 2 2
' i .
i=1
∂xi

Ejercicio 16. Considere la suma z = x1 + x2 + · · · +

xn .
a) Mostrar que si cada xi tiene una cota de error
absoluto ε entonces el error máximo para z es nε.
b) Mostrar que si los errores en cada xi pueden mo-
delarse estadísticamente como variables aleatorias
independientes con desviación√estándar ε entonces
el error estadístico para z es n ε.

Práctica 1a 3