Scribd
Cargar
689 vistas4 páginas

Manu1 U2 A2

Este documento analiza la estabilidad numérica de diferentes algoritmos. Explica que usar la función f(x)=√x+1-√x para valores grandes de x no es estable debido a que su condicionamiento aumenta con x. También analiza un algoritmo para calcular g(x)=(e^x -1)/x, el cual es inestable para valores pequeños de x, y propone cambiar la variable para mejorar la estabilidad.

Cargado por

MaríaLuzPérez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
689 vistas4 páginas

Manu1 U2 A2

Este documento analiza la estabilidad numérica de diferentes algoritmos. Explica que usar la función f(x)=√x+1-√x para valores grandes de x no es estable debido a que su condicionamiento aumenta con x. También analiza un algoritmo para calcular g(x)=(e^x -1)/x, el cual es inestable para valores pequeños de x, y propone cambiar la variable para mejorar la estabilidad.

Cargado por

MaríaLuzPérez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

ANALISIS NUMERICO 1

UNIDAD 2
SISTEMA NUMERICO Y
ANALISIS DE ERROR
ACTIVIDAD 2
ESTABILIDAD DE SISTEMAS

ACTIVIDAD # 2.- ESTABILIDAD DE


SISTEMAS.

1.- Considera la funcin

f (x)= x +1 x . Al calcularla para valores

normales no parece haber ningn problema, pero considera usarla para


4
valores grandes de x , digamos x 10 .
a Obtn la expresin para el condicionamiento de
de

f ( x )= x1 x
u=x +1

u=x

du=1

du=1

f ' ( x )=

1
1

2 x+ 1 2 x

| |

f ' (x ) x
k=
=
f (x )

1
1
x
x

2 x+1 2 x
2 x +1 2 x
k=
=
=
x +1 x
x+1 x
xx
1
=
x +1 x 2 x+1 x

f ( x) dependiendo

b Trabajando con una precisin de 4 cifras significativas, explica en tus


propias palabras por qu la expresin

f s ( x )=

1
x+ 1+ x

es ms

estable.

.Tip: Considera comparar el error absoluto y relativo de valores grandes.


Intenta determinar en qu paso f (x) aumenta su condicionamiento.
'

f s ( x )=

1
2

( x1+1 1x )

Condicionamiento=

x
2

||

1
1
x x x+ 1

( x+1
)
2 ( x+ 1 x )
x
=
=
1
x+1+ x

||

x ( x x+1 ) ( x+ 1+ x )
=
2 ( x+1 x )

1
x +1+ x

x ( xx1 ) =
2 x +1

| | | |

x
1
x
=
2 x+ 1 2 x +1

Podemos notar que en este nuevo clculo, el


condicionamiento de la funcin, nos d el mismo resultado;
en el caso de que se asignan valores grandes, podemos ver
cmo se comporta el error absoluto para dichos valores.
La funcin

g (x )=

e x 1
= x i /(i+ 1)!
x
i=0

se puede calcular, ingenuamente, mediante el siguiente algoritmo


Si x=0

g=1

si no

g=(exp(x )1)/x

fin

a Dnde consideras que este algoritmo se vuelve inestable y por qu?

En g=(exp(x )1)/x

porque en cada caso involucrado la acumulacin

del error ir creciendo.


El algoritmo se vuelve inestable cuando se toman valores pequeos,
dado que al aumentar los valores en los datos a computar, el nmero
fraccionario que se obtiene es mucho menor y se acerca cada vez a 1 al
truncar el valor, provocando con la resta que se colapse a cero.

b Cmo consideras que el cambio de variable en el algoritmo siguiente


contribuye a su estabilidad?
y=exp ( x)

Si y=1

g=1

si no

g=( y1)/ log( y )

Tip: Compara un par de tablas con valores entre

10

fi n
16

y 10

con 10, 15,

y 20 cifras significativas para cada algoritmo para que puedas observar el


efecto numricamente.

En la siguiente tabla con cifras significativas, podemos


observar la manera en la que el valor , en tanto ms
pequeo, har que los algoritmos parezca que se van
alejando; por lo cual, es recomendable cambiar la
exponencial por un logaritmo con el fin de evitar la
inestabilidad o que los cambios no sean tan marcados.
Cifras
significativa
s
10

1er algoritmo

2do algoritmo

1.0000000050

1.000000005

15

1.00000005000000

1.000000004999999

20

1.00000000500000001666

1.0000000049999999991666

También podría gustarte