You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
3. EJERCICIOS 11<br />
Como conclusión el algoritmo es malo. Pero observemos que no lo es el problema en sí mismo.<br />
Como alternativas podemos calcular En por integración numérica o bien hacer el siguiente truco.<br />
Observemos que<br />
y como<br />
1<br />
En ≤<br />
0<br />
En−1 =<br />
1 − En<br />
n<br />
x n dx = 1<br />
→ 0<br />
n + 1<br />
podemos empezar de E20 ∼ 0 e ir hacia atrás usando En−1 = 1−En<br />
n . Este algoritmo es estable<br />
(el error en cada paso se multiplica por algo menor que uno).<br />
Como conclusión, los ejemplos analizados en esta sección muestran la diferencia entre el caso<br />
en el cual la amplificación de los errores de redondeo se debe a que el problema está “mal<br />
condicionado” o “mal planteado” y el caso en el que dicha amplificación se debe al uso de un<br />
“algoritmo inestable”. Es fundamental distinguir entre ambos casos y, por otra parte, encontrar<br />
las causas de la propagación indebida de errores con el objeto de mejorar los algoritmos.<br />
3. Ejercicios<br />
1. Utilizando el método de redondeo:<br />
a) Hallar el número de máquina más próximo a 125,6 y a= 126 si trabaja con<br />
Base 10 y mantisa de 2 dígitos.<br />
Base 2 y mantisa de 8 dígitos.<br />
b) Verificar para x = 125,6, la conocida cota para el error relativo<br />
| x − fl(x)<br />
| ≤ ɛ<br />
x<br />
si ɛ = 1/2β 1−d donde β es la base y d la longitud de la mantisa.<br />
c) ¿Cuál es, en cada caso, el valor que da la máquina como resultado de las operaciones<br />
126 + 125,6 y 126 − 125,6? ¿Cuál es el error relativo de estos resultados?<br />
2. Utilizando el método de truncamiento:<br />
a) Rehacer el Ejercicio 1, con el ɛ correspondiente, es decir: ɛ = β −d+1 , donde β y d<br />
son como antes.<br />
b) Demostrar que, en este caso, ɛ es el menor número de máquina tal que 1 + ɛ = 1.<br />
¿Cuánto da β + ɛ?<br />
3. Mostrar que fl(x) tiene (para ambos métodos) una escritura de la forma<br />
fl(x) = x(1 + δx)<br />
donde |δx| ≤ ɛ. (Usar la cota para el error relativo).<br />
4. Pérdida de dígitos significativos:<br />
a) Si x, y ≥ 0 demostrar que<br />
<br />
x<br />
+ y − fl(fl(x) + fl(y))<br />
<br />
<br />
≤ 2ɛ + ɛ<br />
x + y<br />
2 .