Práctica 2 (Métodos Numéricos)

De Cuba-Wiki
Revisión del 22:35 1 dic 2007 de 200.114.149.135 (discusión) (algunos ej interesantes)
(difs.) ← Revisión anterior | Revisión actual (difs.) | Revisión siguiente → (difs.)

Plantilla:Back

Ejercicio 9

Probar que si tiene todas sus submatrices principales no singulares entonces A tiene factorizacion LU sin pivoteo. Ademas esa factorizacion es unica.

Por induccion en n.
Casos base:
n=1,
n=2, como entonces puedo hacer gauss sin intercambio de filas en entonces se que existen unicos tal que

Paso inductivo:
Supongo que vale
donde son matrices de n-1 x n-1, columnas de n-1 elementos asi como son filas de n-1 elementos.
Quiero ver que vale
Es decir que mis incognitas son
Como los bloques son del mismo tama~no puedo multiplicar por bloques y queda
i) Despejo haciendo: .
ii) Despejo haciendo: .
iii) Despejo haciendo:.
i) y ii) tiene soluciones unicas porque y son inversibles.

Ejercicio 10

Supongamos que una matriz tiene factorizacion A = LU y que L y U son conocidas. Dar una algoritmo que calcule el elemento (i,j) de en aproximadamente flops. (Un flop es una operacion de punto flotante)

Si A = LU entonces y por lo tanto .
Calcular nos lleva O() porque se puede plantear directamente el sistema donde es el canonico con 1 en la posicion j y nos daria la columna j de . Se puede llevar a O() ya que se sabe de antemano que la columna j de tiene j ceros, es decir no es necesario calcularlos. Idem con .

Ejercicio 11

Sea y sea la matriz que se obtiene a partir de A por el metodo de eliminacion Gaussiana cuando las primeras k columnas ya han sido trianguladas. Usando propiedades de determinantes, probar que A es no singular si y solo si es no singular.

Nota: es la matriz identidad con los coeficientes que se usaron en la eliminacion gaussiana para poner un 0 en la posicion i,j.
Por ejemplo para una matriz de 4x4.
Asi,

Volviendo al ejercicio


Ejercicio 24

Sea tal que , siendo cualquier norma consistente.
a) Probar que .
b) Sea una matriz inversible y tal que . Probar que es inversible y vale


a) No se si el enunciado supone que (I+R) es inversible o hay que demostrarlo, por las dudas es asi:

entonces . Si I+R es inversible entonces Nu(I+R) = {0}
Sea .
Da que para que exista un x tal que y entonces lo cual ya sabemos que no pasa. Como cualquier vector se puede llevar a vector de norma 1 el unico x que cumple es 0 entonces es inversible.

Volviendo al ejercicio, por inspiracion divina se me ocurre que

Meto norma en los 2 lados


porque ya que es consistente. El unico numero que cumple esto es 1.

Entonces,

asi que no se puede pasar dividiendo.


b)
i) es inversible. Primero pruebo que es inversible por lo hecho en a). Para eso veo que :
inversible.
Pero entonces porque A es inversible, entonces es inversible.

ii)

en el ultimo paso use que (por enunciado) y asi aseguro que no divide por 0.