Corrección Primer Parcial, Teoría y Práctica, Análisis Numérico I, Semestre I08

Universidad Mayor de San Simón Hans Müller Santa Cruz
Facultad de Ciencias y Tecnologı́a Departamento de Matemáticas

Primer Parcial de Análisis Numérico I 30 de abril de 2008
Nombre y Apellido . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Carnet de Identidad . . . . . . . . . . . . . . . . . . . Firma . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Indicaciones: En hojas en blanco, escriba con letra clara el desarrollo de las preguntas que está respondiendo, indicando claramente
a que pregunta corrresponde.
Parte I (30 puntos)

1. ¿Qué es una operación elemental?
Es una operación matemática que puede ser realizada por un dispositivo de cálculo predeterminado,
en un tiempo finito.
2. Defina problema.
Un problema P(x) es darse un número finito de datos x ∈ Rn y obtener un único resultado P(x).
3. Indique qué es error del método o aproximación.
Cuando se remplaza un problema P(x) por un problema P̃(x), resoluble por algoritmo, se llama error
del método a la diferencia
P̃(x) − P(x).
4. Defina estabilidad en el sentido retrógrado.
Sean P(x) un problema, ŷ la solución numérica por un algoritmo dado, x̂ tal que ŷ = P(x̂), la estabilidad
retrógrada del algoritmo está dada por los valores c > 0, más pequeños que satisfacen
|x̂ − x|
≤ c · eps.
|x|
5. Defina condición de problema.

Sean P(x) un problema dado, la condición de P(x) es el real κ más pequeño tal que
|x̄ − x| |P x̄ − P(x)|
≤ eps ⇒ ≤ κ · eps.
|x| |P(x)|
6. Defina diferencia dividida.

Dada una sucesión (x0 , y1 ), (x2 , y2 ), · · · de puntos, con los xi diferentes entre si. Se define las diferencias
divididas de manera recursiva:
Orden cero: y[xk ] = yk .
y[xi1 , · · · , xin , xin+1 ] − y[xi0 , · · · , xin ]
Orden n + 1: y[xi0 , xi1 , · · · , xin , xin+1 ] =
xin+1 − xi0
7. Enuncie la fórmula de Newton para polinomios de interpolación.
Sea (x0 , y0 ), (x1 , y1 ), . . . , (xn , yn ), con los xi diferentes entre si; entonces el polinomio de interpolación
por estos puntos, está dado por
n
X k−1
Y
p(x) = y[x0 , . . . , xk ] (x − xi ).
k=0 i=0
8. ¿Cuando se interpola una función f , utilizando un polinomio de grado menor o igual a n, la diferencia
dividida de orden n es igual a qué?
f (n) (ξ)
y[x0 , x1 , . . . , xn ] = ,
n!
donde ξ ∈ (mı́n{xi }, máx{xi }).
9. Enuncie el teorema del error de un polinomio de interpolación.
Sea f : [a, b] → R lo suficientemente lisa, {x0 , x1 , . . . , xn } ⊂ [a, b] una colección de puntos diferentes
entre si, p(x) el polinomio de interpolación que satisface p(xi ) = f (xi ), i = 0, . . . , n. Entonces, para
todo x ∈ [a, b], existe ξ ∈ (mı́n{x, xi }, máx{x, xi }) tal que
f (n+1) (ξ)
f (x) − p(x) = (x − x0 )(x − x1 ) · · · (x − xn ).
(n + 1)!
10. Defina los polinomios de Chebichef.
El n-simo polinomio de Chebichef, Pn (x), es el polinomio de grado n, que sobre el intervalo [−1, 1]
satisface:
Pn (x) = cos(n arc cos(x)).
Parte II (70 puntos)

n
Y (x − xj )
1. (30 puntos) Sean x0 < x1 < · · · < xn y li (x) = .
j=0
(xi − xj )
j6=i
Verificar que la función
n
X
λn (x) = |li (x)|
i=0
tiene un solo máximo local sobre cada [xi−1 , xi ].
Denotemos por Ij = [xj−1 , xj ] el intervalo de estudio. Puesto que cada uno de los li (x) tiene sus n
ceros diferentes y éstos son los nudos de la subdivisión, li (x) no cambia de signo en cada uno de los
subintervalos [xk−1 , xk ], en particular sobre Ij . En consecuencia, sobre Ij , se tiene, ya sea |li (x)| = li (x),
o ya sea |li (x)| = −li (x), con lo que se justifica la validez de la indicación.
Para i ∈ {0, 1, . . . , n} determinemos el valor de i en el subintervalo Ij . Separando los i, en aquellos
que son menores o iguales a j − 1 y aquellos que son mayores o iguales a j, se obtiene:
(−1)j−1−i si i ≤ j − 1,

i = (1)
(−1)i−j si i ≥ j.
Por consiguiente, λn (x) es igual sobre el intervalo Ij a
j−1
X n
X
ln (x) = (−1)j−1−i li (x) + (−1)i−j li (x). (2)
i=0 i=j
La identidad (2), muestra claramente que sobre Ij , λn (x) es un polinomio de grado menor o igual a n.
Por consiguiente, prolongamos λn (x) dado por (2) a todo R, denotando la prolongación por T (x).
Remplazando en los xi , el polinomio T (x) de grado menor o igual a n, da los siguientes valores:
(−1)j−1−i si i ≤ j − 1,

T (xi ) = (3)
(−1)i−j si i ≥ j.
Ahora mostremos que T (x) admite un único máximo local so-

bre Ij . Para tal efecto, suficiente mostrar que sobre (xj−1,xj )
existe un único ξ, tal que T 0 (ξ) = 0, además T 0 (x) > 0 sobre
(xj−1 , ξ) y T 0 (x) < 0 sobre (ξ, xj ), (condiciones suficientes pa- 1
ra la existencia y unicidad de máximo local, cf. Cálculo I).
Remarcamos que T (xj−1 ) = T (xj ) = 1, aplicando el Teorema
de Rolle, deducimos la existencia de c ∈ (xj−1 , xj ), con lo que
T (x) tiene al menos un punto crı́tico sobre Ij , mostraremos
que es el único y que no hay otro más en este intervalo.
xj−1 ξ xj
Aplicando la relación (3), el polinomio T (x) cambia de signo en los intervalos [xi−1 , xi ] para 1 ≤ i ≤ i−1
2
y j + 1 ≤ i ≤ n. Como T (x) es una función continua, existen ζi ∈ (xi−1 , xi ) tales que T (ζi ) = 0, para
1 ≤ i ≤ i − 1 y j + 1 ≤ i ≤ n. Por lo tanto existen al menos j − 1 ceros de T (x) menores a xj−1 y
n − j ceros de T (x) mayores a xj , lo que hacen n − 1 ceros; con lo que T (x) es un polinomio de grado
mayor o igual a n − 1. Sabemos que T (x) es un polinomio de grado menor o igual a n. Tenemos dos
posibilidades, que T (x) sea de grado n − 1 o T (x) sea de grado n.
Supongamos que el grado de T (x) es n−1. Para los intervalos [ζi−1 , ζi ], para 2 ≤ i ≤ j−1 y j+1 ≤ i ≤ n,
se tiene T (ζi−1 ) = T (ζi ) = 0, podemos aplicar el Teorema de Roller, con lo que existe al menos un
ci ∈ (ζi−1 , ζi ) con T (ci ). Por lo tanto existen al menos n − 3 ceros de T 0 (x) fuera de Ij , incluyendo
c ∈ Ij , se tiene al menos n − 2 ceros de T 0 (x), como el grado de T 0 (x) es n − 2, se tiene necesariamente
que c ∈ Ij con T 0 (c) es el único número con esta propiedad.
Si el grado de T (x) es n, faltarı́a encontran un cero. Ubiquemos el cero que falta. Utilizando (3), se
obtiene
−1 si n par,

T (x0 )T (xn ) = (−1) (−1)
j−1 n−j
= (−1) n−1
= (4)
1 si n impar.
Ahora bien, como T (x) es un polinomio de grado n, se tiene
1 si n par,

lı́m signo(T (x)) · signo(T (−x)) = (5)
x→+∞ −1 si n impar.
Comparando (4) y (5), se deduce que existe un cero ζ de T (x) que es menor a x0 o mayor a xn .
Aplicando el mismo razonamiento que el caso en que T (x) era de grado n − 1, obtenemos que T 0 (x)
tiene n − 2 ceros fuera de Ij , agregando c de Ij , tenemos exactamente los n − 1 ceros de T 0 (x). Con lo
que está demostrado que c es el único punto crı́tico de T (x) en el intervalo Ij .
cj−1 ζj cj
ζj−1 xj−1 c xj
Solo faltarı́a mostrar que T (x) presenta un máximo local en c y para este efecto, tomamos en cuenta que
T 0 (x) tiene raices simples y consecuencia T 0 (x) no cambia de signo entre dos ceros sucesivos. Además
en la aplicación del Teorema de Rolle, tanto en la determinación de ceros de T (x), como de T 0 (x), se
puede remarcar que entre el primer ci menor a c o el primer ci mayor a c, existe un ζi , cero de T (x).
Por lo tanto si, ζi ∈ (ci , c), ζi < xj−1 , T (x) es creciente y T 0 (x) > 0 o si ζi ∈ (c, ci ), xj < ζi , T (x) es
decreciente y T 0 (x) < 0.
2. (20 puntos) Hallar el polinomio de Hermite p(x) que satisface:
p(0) = 1, p(1) = 0, p0 (1) = 1, p(2) = 0, p0 (2) = −1, p(3) = 1.
Aplicamos el esquema de diferencias divididas
3
0 1
&
% −1
&
1 0 % 2
& &
% 1 % − 23
& &
1 0 % −1 %
3
4
& & &
% 0 % 0 % 0
& &
2 0 % −1 %
3
4
& &
% −1 %
3
2
&
2 0 % 2
&
% 1
3 1
La fórmula de Newton, proveniente del esquema, es:

3 3
p(x) = 1 − x + 2x(x − 1) − x(x − 1)2 + x(x − 1)2 (x − 2).
2 4
3. (20 puntos) Las raices del polinomio x2 − 2px − q = 0 son:

p p
λ1 = p + p2 + q, λ2 = p − p2 + q.
Mostrar que para p > 0, grande, y q ≥ 0, muy pequeño, la segunda fórmula es numéricamente inestable.
Utilizando la relación λ1 λ2 = q, encontrar un algoritmo que sea numéricamente estable
En el cálculo de la segunda fórmula observamos que
p
p ≈ p2 + q,
por lo que la substracción p

p− p2 + q
es una operación mal condicionada. Proponemos el algoritmo
p
λ1 = p + p2 + q,
q
λ2 =
λ1
Todas las operaciones involucradas en estas fórmulas son operaciones bien condicionadas.

Corrección Primer Parcial, Teoría y Práctica, Análisis Numérico I, Semestre I08

Cargado por

Copyright:

Formatos disponibles

Corrección Primer Parcial, Teoría y Práctica, Análisis Numérico I, Semestre I08

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Corrección Primer Parcial, Teoría y Práctica, Análisis Numérico I, Semestre I08

Cargado por

Copyright:

Formatos disponibles

Universidad Mayor de San Simón Hans Müller Santa Cruz

Facultad de Ciencias y Tecnologı́a Departamento de Matemáticas

Parte I (30 puntos)

5. Defina condición de problema.

6. Defina diferencia dividida.

Parte II (70 puntos)

Ahora mostremos que T (x) admite un único máximo local so-

2. (20 puntos) Hallar el polinomio de Hermite p(x) que satisface:

p(0) = 1, p(1) = 0, p0 (1) = 1, p(2) = 0, p0 (2) = −1, p(3) = 1.

Aplicamos el esquema de diferencias divididas

La fórmula de Newton, proveniente del esquema, es:

3. (20 puntos) Las raices del polinomio x2 − 2px − q = 0 son:

por lo que la substracción p

También podría gustarte