Práctica 1, Análisis Numérico I
Práctica 1, Análisis Numérico I
Práctica 1, Análisis Numérico I
M1 Supóngase que, el dispositivo de cálculo, con el que se cuenta, puede efectuar la división con resto; es
decir, para a, b enteros no negativos, con b 6= 0, el dispositivo calcula p, r satisfaciendo:
a = pb + r y 0 ≤ r < b.
bn = an ;
bi = ai + x0 bi+1 , i = n − 1, . . . , 1, 0.
x0 = 1,
xn+1 = xn (2 − xn v).
Mostrar que para p > 0, grande, y q ≥ 0, muy pequeño, este algoritmo es numéricamente inestable.
Utilizando la relación λ1 λ2 = q, encontrar un algoritmo que es numéricamente estable en el sentido de
backward analysis.
M12 Asumiendo la Tierra esférica, dos puntos sobre esta superficie pueden ser expresados en coordenadas
cartesianas
pi = (xi , yi , zi ) = (r cos αi cos βi , r sin αi cos βi , r sin βi ), i = 1, 2,
donde r es el radio de la Tierra y αi , βi son las longitudes y latitudes de los puntos pi , respectivamente.
Si
pt p2
cos δ = 1 2 = cos(α1 − α2 ) cos β1 cos β2 + sin β1 sin β2 ,
r
entonces rδ es la distancia entre estos dos puntos.
(a) Mostrar que utilizando la función arc cos para determinar δ de la expresión de más arriba no es
numéricamente estable.
(b) Derive una expresión numéricamente estable para δ.
M13 Mostrar como evaluar las siguientes expresiones de una manera numéricamente estable:
1
1+2x − 1−x para |x| 1,
q q 1+x
x + x1 = x − x1 para x 1.
2
M14 Supóngase que la librerı́a de programas de una computadora contiene una subrutina que evalua arcsin y
con una precisión en punto flotante de t decimales. Considérese la relación
x
arctan x = arcsin √ .
1 + x2
Determine para que valores de x la fórmula permite calcular establemente arctan.
M15 Para un x y un entero k, el valor de cos kx puede ser calculado recursivamente utilizando para m =
1, 2, . . . , k − 1 la fórmula
Por consiguiente, el cálculo de cos kx puede ser efectuado con una sola evaluación de cos x. Planteando
c1 = cos x y c0 = 1, deduzca el algoritmo para evaluar cos kx. Analize la estabilidad del algoritmo para
|x| =
6 0 pequeños.
M16 Para un x dado y un entero k grande, cos kx y sin kx pueden ser evaluados recursivamente, utilizando
las fórmulas:
cos mx = cos x cos(m − 1)x − sin x sin(m − 1)x;
sin mx = sin x cos(m − 1)x + cos x sin(m − 1)x, m = 1, . . . , k.
Planteando cm = cos mx, sm = sin mx, estudiar la estabilidad del algoritmo propuesto por las relaciones
recursivas, se debe notar que para el cálculo solo se requiere una evaluación de sin x y otra de cos x.
M17 Para |x| = 6 0 pequeño, el algoritmo de M16 puede ser mejorado en su estabilidad, definiendo las
diferencias
dcm+1 = cm+1 − cm ,
dsm+1 = sm+1 − sm .
y para m = 1, 2, . . . , k: