Meto 2 PDF
Meto 2 PDF
Meto 2 PDF
luego
64 64 16 16
Asi que,
Como
r 0 v = f{ 3 x i) = t{ ~1) = 28
0
y
entonces
20 3) =fl ( 8
(r0v)e(r 0 s)=fl ( 8-"4 14) 28
=16
10 METOOOS NUMERICOS
Entonces
1.2 ERR ORES DE REDONDEO
Sabemos que todo numero real X;t 0 puede escribirse en la forma decimal normalizada
siguiente
c) x - - 123456789
Supongamos, para simplificar el analisis de los errores de redondeo, que nuestro conjunto de .. (.1234567
d) Y=.000021 3475
ii. Redondeando el numero x: En este caso
(213475)
si 0 ~ a t + 1 < 5
Entonces
Ejemplo 1.1 Supongamos t = 5 Y usemos las reglas de redondeo y cortado para encontrar
el representante de punto flotante decimal en cada uno de los siguientes casas :
Entonces
(27182) x 10 1, cortando
fl(e) ::
1
\ (27183) x 10 , redondeando ( ya que as = 8 > 5 )
Capitulo 1. ERRORES DE REDONDEO Y ESTABILIDAD 11
Entonces
( .314 15) x 1 0 1,
cortando
fl(rr) =
(.3 1416) x 10\ redondeando
c) X = - 123456789 (racional)
= - (.1 23456789) x 10 9 (forma decimal normalizada)
Entonces
_(.12345) x 10 9 , cortando
fl( X) =
-(.12346) x 10 9 , redondeando
d) Y=.0000213475 (racional)
= (.213475) x 10-4 (forma decimal normalizada)
Entonces
(.21347) x 10 -4 , cortando
fl(y) =
(.21348) x 10 -4 , redondeando
Entonces
(.66666) x 10° , cortando
fl(z) =
1(.66667) x 10 ° , redondeando •
Hay varias formas acostumbradas para medir errores de aproximaci6n ; algunas de elias se
dan en la siguiente definici6n.
12 METODOS NUMERICOS
b) Si 5 s a t+1 s 9 , entonces
Definici6n 1.1 Sea x· una aproximaci6n de un numero real x. EI error de x con
respecto a x es E= x - x· ; el error absoluto de x· con respecto a x es E = x - x·1 y el I asique
error relativo de x· con ..specto a x, "0, es Er " IXI ~~. I. Tambi"'n se define el error
I (.a ,a2 ···a,a,." ,,)
1
Er ='---- - - - , / -
E =I x - fl( x) I y Er =
Ix-fl(x)
Ix I ,
I 0
X;o;
•
Va vimos que el error de redondeo puede depender del tamal'\o del numero, pues los
numeros de punto flotante no estan distribuidos de manera uniforme en la recta real; desde
este punto de vista el error relativo es una mejor medida del error de redondeo que el error
absoluto.
.1
Estimemos la menor cota superior para el error relativo cuando un numero real . x ;0; 0 es
aproximado por su representante de punto fiotante, fl( x), en una aritmetica decimal de t
digitos.
Sea
Observe, en el trabajo
yentonces
Er =~----~----------------~----~
I(.a,a2 .. ·atat+, ···) x lO n I
I (.aI+ 1a t +2 " ') x lO n- t I
I(.a,a2 .. ·atat+ ",,) x 10 n
I (.at+ at+2·..) I
1 x lO - t
I (.a , a2.. ·atal+ , "') I
Capitulo 1. ERRORES DE REDONDEO Y ESTABILIDAD 13
b) Si 5 ~ a l +1 ~ 9 , entonces
l1li X. EI error de x' con
I
• x es E = x - x'i y el
asique
n
I (.a,a 2.· a l a l +,··.) x 10
n
- [(.a,a 2.. a l ) x lO +
n l
x lO - ] to I
I (.a t+ 1a t+ 2 · ·) - 1.0 I -I
• = x 10
.5
del numero, pues los
en la recta real; desde
~I ( . )l x 10
.a 1a 2·.. a ta +1···
l
-t
, yaque .al + ,al+ 2 ·· · ~ ·5
aritmetica decimal de t-
i
posici6n t + 1
Ix- fl(x) I t
Er = < 5 x 10
Ix I
Ejemplo 1.2 Encuentre el error absoluto y el error relativo de x' con respecto a x, en cada
uno de los siguientes casos:
14 METODOS NUMERICOS
2 2 2 1
E = I (.5) x 10 - (.5 1) x 10 1= 1- (·01) x 10 1=(·1) x 10 =1.0
.02 = 2%
(.5) x 10-3 .5 50
Er = 6
== 1
.02 = 2% •
(.5) x 10 (.5) x 10 50
.00005 < I y- y" 1=
Este ejemplo nos muestra que el error relativo es invariante aJ cambio de escala y se usa
y nuevamente, y" aproXlml
Teniendo en cuenta la menor cota superior para el error relativo usando redondeo , se define
e/ concepto de cifras significativas.
Deflnici6n 1.2 Se dice que e/ numero x· aproxima con sus primeros t-digitos 0 cifras
significativas al numero x * 0 , si t es el mayor entero no negativo para e/ cual
x - x" 1
Ix I
Los t-dfgitos significativos, a que se refiere esta deflnici6n, son 105 primeros t-dfgitos en la
Deflnici6n 1.3 Se dice que el numero x· aproxima con sus primeras k-cifras decimaJes
exactas a/ numero x, si k es el mayor entero no negativo tal que
C a pitulo 1. ERRORES DE REDONDEO Y ESTABILIDAD 15
Las k cifras decimales exactas, a que se refiere esta definicion, son las primeras k cifras
contadas a partir del punto decimal en x" , cuando x" se escribe en forma decimal. V
Los dos conceptos anteriores pueden aparecer definidos de manera distinta en otros textos.
Aqu i se usaran las definiciones dadas.
I
.00005 < x - x"1= .0001 03 < .0005 = 5 x 10- 4 < 5 x 10- 3 < 5 x 10- 2 < 5 x 10 1
asi que k =3 es el mayor entero no negativo ta l que 1 .003451- .003348 1:<; 5 x 1O - (k+1) .
Luego .0 03 348 aproxima a .0034 51 con sus tres primeras cifras decimales exactas , que son
en este caso 0, 0 y 3.
y nuevamente , y" aproxima a y con su s primeras tres cifras decimales exactas, que son , par
supuesto, 0, 0 y 3,
.000103 1
,005 < Er = .029,< .05 = 5 x 10- 2 < 5 x 10 -
.003451
1 .003451- .003348 1 t
!.-.-- - -- - - - - ' - < 5 x 10
1.0034511 '
y por tanto x' aproxima a x con sus primeros 2-digitos significativos que son 3 y 3 (Por que?)
Con cuantas cifras significativas aproxima y' a y? •
1
Ejemplo 1.4 Con cuantas cifras significativas aproxima .333 a '3?
Como
I i-·333
J
~ -.333
1 ~3_ _1= 11- .999 1= .001
Ii I 3
16 M~TOOOS NUM~RICOS
Y .0005 < .001 < .005 = 5 x 10- 3 < 5 x 10- 2 < 5 x 10- 1 , entonces t = 3 es el mayor entero no
por tanto I error rela
negativo tal que
--.333
1 I
\ 3 5 10- t
I~ I < x
1
Por 10 tanto .333 aproxima a 3" con 3 cifras significativas. Observe que .333 es el numero
con respecto a
numeros "caSI Ig
en aritmetica de punto flotante decimal con redondeo a tres digitos que representa a 3"' • operacl nes
Ejemplo 1.5 D6nde debe estar x· para que aproxime a 1000 con 4 cifras significativas?
1000 - x·
1000
1<5x 10- 4
1000 - x·
1000
I~ 5 x 10
5
La desigualdad i) tiene como soluci6n 999.5 < x· < 1000.5 Y la desigualdad ii) tiene como
obtene 5
Se sa be que fl(x) Y fl(y) aproximan a x e y, respectivamente, con todas sus seis · cifras
significativas (Verifiquelo).
Ahora ,
x - y=-1.53x10 - 6 = -.153 x 10 - 5
= fl(-2.0 x 10 - 6
) = fl( - .2 x 10-5 )= -2 x 10-
Capitulo 1. ERRORES DE RED ONDEO Y ESTABILIDAD 17
5
1- .2 x 10- - (- .153 x 10-
5
) 1_ .047 _
.307 ... < .5 = 5 x 10- '
1-.153 x 10- 5 1 - .153
Luego x e y aproxima al valor exacto x - y con unicamente una cifra significativa (1), asf
que hubo perdida de 5 cifras significativas (fl(x) , fl(Y) ten /a n cada uno 6 cifras signiticativas
co n respecto a x e y, respectivamente); 10 anterior sugie re que debe evitarse la resta de
numeros "casi iguales". Como ejercicio, revise en 'el mismo ejemplo, que pasa con las
operaciones (fJ , ® y 8
x2 - 4002x + 80 = 0
Si hacemos los calcu los para Xl y x 2 , usando aritmetica decimal con redondeo a 4 digitos,
obtenemos
x 2
= - -- - -2- - - - .1500
2 2 2
La deficiencia en la estimaci6n de xjl se debe a que los numeros 400.2 y ~(400.2)2 - 320
son nu meros muy cercanos entre sf (en una aritmetica finita con redondeo a 4 dfgitos). En
este caso se consigue una aproximacion mas exacta para X2, aumentando la precision de la
aritmetica 0 "racionalizando el numerador"
= 80l I 2
400.2 + ,,(400.2)2 - 320
1~=
x1
, donde c es el termino constante en la ecuaci6n
los primeros n + 1 I,••m-.,.
2
x + bx + C = 0 , obtenemos
• 80 80
x 2 = -- = - - = .2000
x~ 400.1
que coincide con el valor exacto de x 2 , en este caso . Supongamos que qUI"'"
C6mo resolve ria la ecuaci6n x 2 + 400.2x + 80 = 0 , usando aritmetica decimal con redondeo a
cuatro digitos, si quiere intentar evitar la perdida de cifras significativas en el calculo de las
raices? *
\ Ejercicio 1.1 Elabore un programa de computador que resuelva la ecuaci6n cuadratica
general ax 2 + bx + C = 0 (aun en el caso de ra ices complejas), usando aritmetica finita y que
intente evitar la perdida de cifras significativas en el calculo de las raices . *
Ejemplo 1.7 Recordemos que para todo x E R
xn x 2 x3 xn
eX = I
00
EI polinomio
x2 x3
xn
Pn(x) = 1 + x + - + --+.,.+
2! 3! n!
Se sabe que
con
n+1
Rn{x;O) = eE, _ ( x) para algun S entre 0 y x
n + 11
Capitulo 1. ERROR ES DE REDONDEO Y ESTABILIDAD 19
o tambiem
Observe que Rn(X;O) no es otra cosa que el residuo en la serie de Taylor cuando se to man
los primeros n + 1 terminos. A Rn(X;O) se Ie lIamara error de truncamiento 0 de formula al
5
Cual es la aproximaci6n que se obtiene para e- 5 y e , si se tra baja en una aritmetica (de
punta flotante) decimal con redondeo a 4 digitos?
5
Las aproximaciones correspondientes a e- 5 y e aparecen en la TABLA 1.3.
De acuerdo con los resultados de la TABLA 1.3, en una aritmetica decimal con redondeo a 4
n (-5/
dfgitos, e-s "" .9993 x 10- 2 (Ia suma Ik~O --
kl
se estabilizo en n = 22) Y e 5 "" 148.4 (Ia
n k
suma I
k=O
~
kl
se estabilizo en n = 14 ) .
Una forma mas adecuada de calcu lar e- 5 es calcular para la aritmetica de punto
e5
flotante decimal con redondeo a cuatro digitos
20 METODOS NUMERICOS
14
.7001 x 10- 1
.2771 x 10 - 1
148.4
15
-.2333 x 10- 1
.4380x10- 2 148.4
16
.7294 x 10-2
.1167x10 - 1
17
-.2145 x 10 - 2
.9525 x 10- 2
18
.5959 x 10- 3
.1012x10 - 1
19
-.1568 x 10- 3
.9963 x 10 - 2
20
.3920 x 10-4
.1000 X 10- 1
21
-.9333 x 10- 5
.9991 x 10-2
22
.2121 x 10- 5
.9993 x 10-2
23
-.4611 x 10-6
.9993 x 10- 2
TABLA 1.3
Los ejemplos 1.6 y 1.7 anteriores, muestran como un algoritmo mal concebido puede
conducir a una respuesta defectuosa de un problema perfectamente bien planteado. La
deficiencia fue corregida cambiando el algoritmo.
,
In = Jxnex-'dx, n = 1,2,3,..
o
Usando aritmetica (de punto flotante) decimal con redondeo a 6 digitos y la formula de
recurrencia In = 1- nln_, ' obtenemos
I, "" .367879 = I; , 12 "" .264242 = I;, 13 "" 207274 = I;, 14 "" .170904 = I~ ,
15"" .145480=1; , 16 "" .127120 = 1 ~, 17'" .110160=1;, Ie'" .118720 = 1 ~,
19 "" - .0684800 = I~
g
Es claro que el valor -.0684800 ('" 19) es incorrecto, pues x e X- ' es continua y positiva
sobre el intervalo (0,1) . Que caus6 este resultado? Observe que unica mente hay error de
1
redondeo en el calculo de I" donde - fue redondeado a 6 digitos sign ificativos. Como la
e
f6rmula de recurrencia obtenida en la integraci6n por partes es exacta para la aritmetica real ,
entonces no hay error de f6rmula y as! el error en 19 es debido en su totalidad al error de
redondeo en I, . EI error inicial fue E "" 4.412 x 10 - 7 .
entonces 12 - I; = - 2 E.
Ahora ,
22 METODOS NUMERICOS
De donde
19 -I~ :0:362880(4.412
7
X 10- ),., .160102656
19 = -.0684800+.160102656 = .091622656
Observe que el error absoluto, debido a los calculos, crece a medida que n aumenta, y es
mucho mas grande que el valor real (en valor absoluto) que se esta aproximando (se puede
ver que si E es el error inicial, entonces el error despues de n pasos es
1
In = 1- nln-1' n = 2,3, ... con 11 =
e
es inestable.
1-1
In- 1 = __ n , n = N, ... ,3,2
entonces en cada paso del calculo el error en In es dividido por n . Asi que, si comenzamos
con un valor para algun In con n » 1, Y trabajamos hacia atras, cualquier error inicial 0
errores de redondeo que ocurran estaran decreciendo en cada paso. Este es un ejemplo de
algorltmo astable.
Por ejemplo, si aproximamos 120 por 0 y usamos el valor 0 como un valor inicial , entonces
cometemos un error inicial E tal que 0 ~ E ~ ~ ; este error es multiplicado por _1_ al calcular
21 20
1,9 , as; que el error en el calculo de 1,9 , que es _1_ ,es tal que
20
1 1 1
o ~ -E ~
20 2021
1 1 1 1 1 1 1 -8 -8
16 17 20 16 17 20 21
10 que garantiza una precision de por 10 menos 6 cifras decimales exactas para los valores
calculados de 1'5 ,... ,1 9 ,
a medida que n aumenta , y es Haciendo los calculos para 120 ,... ,19 , obte nemos
toe estti aproximando (se puede
despu~s de n pasos es 120 ::::: .0000000000 , 1,9::::: .05000000000 , 1'8 '" .0500000000,
1
117 ::::: .0527777778 , 1'6 : : : .05571895425, 1,5 '" .05901 756536,
que 0 < In ~ - - ) . En
n+1
1'4 '" .06273216231, 1'3 '" .06694770269 , 1'2 '" .07177325364 ,
Supongamos que Xk = (~)Xk-l - (~)Xk-2 = (~r para 2:0; k < n, y veamos que
Si comlJ)ll:Ift'II.
xn - -(~)x
6 n-l -(~)x
6 n-2 _(~)n
- 3 .
Luego
x
n- (5)
- -6x n-l - -
6 x
n-2 -
- -
3
(1) (1) n para todo n = 0,1, ...
Analogamente, se puede verificar que la sucesi6n definida en ii) es igual a la sucesi6n {Pn} n
Si usamos aritmetica decimal con redondeo a 7 digitos para calcular los primeros terminos de
las sucesi6nes {Pn} n' {xn t Y {y n} n' se obtienen los resultados que se muestran en la
TABLA 1.4 siguiente.
n
Pn
. x~ y~
0 1.000000 1.000000 1.000000
.3333333 .3333333 .3333333
2 .1111111 .1111111 .1111111
3 -1
4
5 -2
6 -2
7 -3
8 -3
9 -4
10 -4
11 -5
12 -5
13 -5
14 -6
15 -6
16
17 -5
18
19 -5
20 -5
TABLA 1.4