2022-2 Analisis Numerico EXAMEN UNAM

Examen de Conocimientos Generales
Análisis Numérico
Semestre 2022-2
Julio 2022
Instrucciones: Resuelva todos los ejercicios explicando con detalle su procedimiento.

Calificación: Suma de puntos dividida por 3. Requiere al menos 18 puntos para aprobar.
Duración: 4 horas.
1 Representación de los números de punto flotante. Los números reales se pueden

expresar en potencias de 10. Por ejemplo:
∞
34 1 0
X 3
= 11.3333 . . . = 1 × 10 + 1 × 10 + .
3 k=1
10k
En general, un número real x se representa en base decimal como:

∞
n
!
X
k
X bn
x = ±(an · · · a1 a0 .b1 b2 b3 . . .)10 = ± ak 10 + ,
k=1 k=1
10k
donde ak , bk ∈ {0, 1, . . . , 9}. La primera suma es la parte entera, mientras que la

segunda suma es la parte fraccionaria. Los dı́gitos son los coeficientes ak y bk que
acompañan a las potencias de 10. En particular, se le pide verificar esta representa-
ción para los números enteros.
a) (2 puntos) Sea M ∈ Z. Pruebe que existen d1 , . . . , dp ∈ {0, . . . , 9} tales que:
|M | = d1 10p−1 + · · · + dp−1 101 + dp 100 .
Los números reales se redondean con una cantidad finita de dı́gitos para represen-
tarlos en una computadora. La parte fraccionaria se trunca a una suma finita. Por
ejemplo, si usamos solo cuatro dı́gitos después del punto decimal, entonces x = 34/3
se aproxima por
3 3 3 3
b = 11.3333 = 1 × 101 + 1 × 100 +
x 1
+ 2 + 3 + 4.
10 10 10 10
No es recomendable fijar la posición del punto decimal entre 100 y 10−1 para
redondear números reales con distintas potencias de 10, ya que el error relativo de
redondeo puede cambiar drásticamente. Por ejemplo, si x = 5.9901 e y = 0.0059901
1
se truncan usando tres dı́gitos después del punto decimal, entonces x se aproxima
por xb = 5.99 con un error relativo |x − x b|/|x| ≈ 1.67 × 10−5 , mientras que y se
aproxima por yb = 0.005 con un error relativo |y − yb|/|y| ≈ 0.16.
La posición del punto se cambia junto a la potencia más grande de 10 para evitar
cambios abruptos en el error relativo de redondeo. Por ejemplo: y = 0.0059901
se escribe como y = 5.9901 × 10−3 , se trunca a tres dı́gitos significativos como
y = 5.99 × 10−3 , y su error relativo de redondeo es |y − yb|/|y| ≈ 1.67 × 10−5 . Se
puede probar que un número representado por una cantidad finita de dı́gitos se puede
escribir en notación cientı́fica.
b) (2 puntos) Un número en base decimal con p dı́gitos significativos es de la forma:

x = M · 10e−p+1 ,
donde M, e ∈ Z tales que |M | < 10p . Pruebe que existen d1 , . . . , dp ∈ {0, . . . , 9}
tales que:
d2 dp
x = ± d1 + + · · · + p−1 · 10e .
10 10
Esto se escribe en notación cientı́fica como:
x = ±(d1 .d2 · · · dp )10 × 10e .
2 Matrices Definidas Positivas. Algunas matrices simétricas dependen de parámetros.

Estas matrices pueden ser definidas positivas solamente para algunos valores de los
parámetros.
Sean a, b > 0. Considere la matriz
b −1 −1 · · · −1
 
 −1 a 0 ··· 0 

. .. .. 
A =  −1 0 .. . . ,
 
 . .. . .
 ..

. . a 0 
−1 0 ··· 0 a n×n
El problema es hallar condiciones suficientes sobre los valores de a y b de modo que
A sea una matriz definida positiva.
Una manera para identificar fácilmente los valores deseados de a y b es intercambiar
los renglones y columnas de la matriz A de modo que A se transforme en la matriz:
0 ··· 0 −1
 
a
. .. .. 
 0 a .. . . 

A
b= .
 .. . . . . .. .
.. .

 0 
 0 ··· 0 a −1 
−1 · · · · · · −1 b n×n
2
Si se obtienen valores de a y b para los que A
b es definida positiva, entonces se puede
probar que la matriz A es definida positiva con esos mismos valores.
a) (1 punto) Halle una matriz de permutación P de tamaño n × n tal que
A = P T AP.
b
b) (2 puntos) Tome n = 5, a = n − 1 y b = 1. Pruebe o refute que A

b es una matriz
definida positiva.
c) (2 puntos) Sea n = 7. Halle condiciones sobre los valores de a y b de modo que
Ab sea una matriz definida positiva. Sugerencia: Factorize la matriz A.
b
d ) (1 punto) Pruebe que si Ab es definida positiva, entonces A es definida positiva.
3 Curvas Parametrizadas por Polinomios Cúbicos. El diseño asistido por computadora

requiere de curvas suaves que puedan modificarse de manera fácil y rápida. Una
familia de curvas que pueden emplearse está parametrizada por polinomios cúbicos
de Hermite.
Dados los puntos (x0 , y0 ) y (x1 , y1 ), se desea obtener una curva parametrizada:
φ(t) = (x(t), y(t)) : [0, 1] → R2 ,
donde x e y son polinomios cúbicos que conectan ambos puntos:
x(0) = x0 , x(1) = x1 ,
y(0) = y0 , y(1) = y1 .
La recta tangente a la curva φ en el punto (x0 , y0 ) puede especificarse mediante un

punto guı́a (x0 + α0 , y0 + β0 ), donde
α0 = x0 (0), β0 = y 0 (0).
Análogamente, la recta tangente a la curva φ en el punto (x1 , y1 ) puede especificarse

mediante un punto guı́a (x1 − α1 , y1 − β1 ), donde
α1 = x0 (1), β1 = y 0 (1).
a) (2 puntos) Dados αi , βi , i = 0, 1. Halle expresiones para los polinomios cúbicos

de Hermite x e y en la base canónica {1, t, t2 , t3 }.
b) Sea φ la curva parametrizada por polı́nomios cúbicos de Hermite que une
(x0 , y0 ) = (0, 0) con (x1 , y1 ) = (0, 1) mediante puntos guı́a:
(x0 + α0 , y0 + β0 ) = (2, 2), (x0 + α1 , y0 + β1 ) = (2, −1).
3
Figura 1: Curva parametrizada φ
i. (1 punto) Pruebe que la curva φ atraviesa el eje horizontal.

ii. (1 punto) Dibuje la curva φ.
4 Integración Múltiple por Cuadratura Gaussiana. Se desea aproximar el volumen deba-

jo de una función gaussiana de variables x e y sobre el cuadrado unitario [0, 1] × [0, 1]
mediante una regla de cuadratura:
Z 1Z 1 2 4
x + y2

1 X
exp − dxdy ≈ wi p(xi , yi ), (1)
2π 0 0 2 i=1
donde p(x, y) es un polinomio bicuadrático.

a) (3 puntos) Halle pesos wbi y nodos xi tales que para cualquier polonimio q de
hasta grado tres se cumpla:
Z 1
q(x)dx = w
b1 q(x1 ) + w
b2 q(x2 ) + w
b3 q(x3 ).
−1
Desarrolle su procedimiento.
b) (1 punto) Encuentre el polinomio cuadrático q que interpola la función:
2
exp − x2
f (x) = √ ,
2π
en x1 = −1, x2 = 0 y x3 = 1.
c) (2 puntos) Sea p(x, y) = q(x)q(y). Use la aproximación:
2
x + y2

1
exp − ≈ p(x, y), 0 ≤ x, y ≤ 1,
2π 2
para obtener los pesos wi y los nodos (xi , yi ) de la regla de cuadratura (1).
Sugerencia: Use un cambio de variable adecuado.
4
5 Teorema de Punto Fijo. Sea φ : [a, b] → [a, b] una función en C 1 ([a, b]) tal que:
máx |φ0 (x)| < 1. (2)
a≤x≤b
Entonces, dado x0 ∈ [a, b], la sucesión:
xk+1 = φ(xk ), k ∈ N,
converge a un único punto α ∈ [a, b] tal que φ(α) = α.

Use las hipótesis del Teorema de Punto Fijo para probar lo siguiente:
a) (1 punto) La sucesión {xk } converge.

b) (1 punto) Si 0 < φ0 (α) < 1 y x0 > α, entonces xk > α para k ∈ N.
c) (1 punto) Si −1 < φ0 (α) < 0, entonces xk − α y xk+1 − α tienen signos opuestos
para k ∈ N.
La cota sobre la pendiente es una condición suficiente para garantizar que la

iteración de punto fijo {xk } converga. Si la cota (2) no se cumple en todos los puntos
del intervalo, la sucesión {xk } puede ser divergente.
d ) (1 punto) Pruebe que si |φ0 (α)| > 1, entonces {xk } no puede converger.
6 Ajuste de una sección cónica por mı́nimos cuadrados. Considere la familia de fun-
ciones bicuadráticas ϕα : R2 → R dadas por
ϕα (x, y) = α1 x2 + α2 y 2 + α3 xy + α4 x + α5 y + 1,
donde
α5 ∈ R5 .

α = α1 · · ·
Se desea ajustar una sección cónica:
ϕα (x, y) = 0,
a los puntos (xi , yi ), i = 1, . . . , 8 dados en la Tabla 1 minimizando la suma de

cuadrados:
X8
J(α) = ϕα (xi , yi )2 .
i=1
5
i xi yi
1 √ 0 1
2 2 − 3 1/2
3 √2 0
4 2 + 3 1/2
5 √4 1
6 2 + 3 3/2
7 √2 2
8 2 − 3 3/2
Tabla 1: Puntos (xi , yi ), i = 1, . . . , 8.
a) (2 puntos) Encuentre una matriz A ∈ R8×5 y un vector b ∈ R8 tales que
J(α) = kAα − bk22 .
b) (2 puntos) Supónga que:

A = QR,
donde Q una matriz de tamaño 8×5 con columnas ortonormales y R una matriz
triangular superior de tamaño 5 × 5 tal que todos sus elementos en la diagonal
principal son distintos de cero. Pruebe que el mı́nimo de J está dado por la
solución del sistema de ecuaciones:
Rα = Qt b.
c) (1 punto) Obtenga las ecuaciones normales correspondientes y compruebe que

T
αmı́n = 1/4 1 0 −1 −2
es la solución de esas ecuaciones.

d ) (1 punto) Dibuje la sección cónica ϕα (x, y) = 0 que se obtiene con α = αmı́n .

2022-2 Analisis Numerico EXAMEN UNAM

Cargado por

Copyright:

Formatos disponibles

2022-2 Analisis Numerico EXAMEN UNAM

Cargado por

Información del documento

Descripción original:

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

2022-2 Analisis Numerico EXAMEN UNAM

Cargado por

Copyright:

Formatos disponibles

Examen de Conocimientos Generales

Instrucciones: Resuelva todos los ejercicios explicando con detalle su procedimiento.

1 Representación de los números de punto flotante. Los números reales se pueden

En general, un número real x se representa en base decimal como:

donde ak , bk ∈ {0, 1, . . . , 9}. La primera suma es la parte entera, mientras que la

a) (2 puntos) Sea M ∈ Z. Pruebe que existen d1 , . . . , dp ∈ {0, . . . , 9} tales que:

|M | = d1 10p−1 + · · · + dp−1 101 + dp 100 .

b) (2 puntos) Un número en base decimal con p dı́gitos significativos es de la forma:

2 Matrices Definidas Positivas. Algunas matrices simétricas dependen de parámetros.

a) (1 punto) Halle una matriz de permutación P de tamaño n × n tal que

b) (2 puntos) Tome n = 5, a = n − 1 y b = 1. Pruebe o refute que A

3 Curvas Parametrizadas por Polinomios Cúbicos. El diseño asistido por computadora

φ(t) = (x(t), y(t)) : [0, 1] → R2 ,

donde x e y son polinomios cúbicos que conectan ambos puntos:

La recta tangente a la curva φ en el punto (x0 , y0 ) puede especificarse mediante un

Análogamente, la recta tangente a la curva φ en el punto (x1 , y1 ) puede especificarse

a) (2 puntos) Dados αi , βi , i = 0, 1. Halle expresiones para los polinomios cúbicos

(x0 + α0 , y0 + β0 ) = (2, 2), (x0 + α1 , y0 + β1 ) = (2, −1).

i. (1 punto) Pruebe que la curva φ atraviesa el eje horizontal.

4 Integración Múltiple por Cuadratura Gaussiana. Se desea aproximar el volumen deba-

donde p(x, y) es un polinomio bicuadrático.

Entonces, dado x0 ∈ [a, b], la sucesión:

converge a un único punto α ∈ [a, b] tal que φ(α) = α.

a) (1 punto) La sucesión {xk } converge.

La cota sobre la pendiente es una condición suficiente para garantizar que la

a los puntos (xi , yi ), i = 1, . . . , 8 dados en la Tabla 1 minimizando la suma de

a) (2 puntos) Encuentre una matriz A ∈ R8×5 y un vector b ∈ R8 tales que

J(α) = kAα − bk22 .

b) (2 puntos) Supónga que:

c) (1 punto) Obtenga las ecuaciones normales correspondientes y compruebe que

es la solución de esas ecuaciones.

También podría gustarte