MATRICES ELEMENTALES
Profesores Omar Dar o Saldarriaga Ort z Ivan Dar o G omez Hern an Giraldo
2009
Definici on Sea E una matriz de tama no n n, decimos que E es una matriz elemental si E se obtiene de la identidad al aplicar una operaci on elemental de la. Ejemplo 0 Las siguientes matrices son matrices elementales: A = 1 0 1 2 0 1 0 0 B = 0 1 0 C = 0 3 0 . 0 0 1 0 0 1 1 0 0 0 0, 1
Definici on Sea E una matriz de tama no n n, decimos que E es una matriz elemental si E se obtiene de la identidad al aplicar una operaci on elemental de la. Ejemplo 0 Las siguientes matrices son matrices elementales: A = 1 0 1 2 0 1 0 0 B = 0 1 0 C = 0 3 0 . 0 0 1 0 0 1 1 0 0 0 0, 1
Observaciones Como hay tres tipos diferentes de operaciones elementales, hay un n umero igual de tipos de matrices elementales entonces usaremos la siguiente notaci on. Eij denotar a la matriz elemental que se obtiene al intercambiar las las i y j de la matriz identidad. Eij (c) denotar a la matriz que se obtiene al sumar c veces la la i a la la j de la matriz identidad. Ei (c) la matriz que se obtiene al multiplicar por la constante c la la i de la matriz identidad. Ejemplo En el ejemplo anterior tenemos que A = E12 , B = E21 (2) y C = E2 (3).
Observaciones Como hay tres tipos diferentes de operaciones elementales, hay un n umero igual de tipos de matrices elementales entonces usaremos la siguiente notaci on. Eij denotar a la matriz elemental que se obtiene al intercambiar las las i y j de la matriz identidad. Eij (c) denotar a la matriz que se obtiene al sumar c veces la la i a la la j de la matriz identidad. Ei (c) la matriz que se obtiene al multiplicar por la constante c la la i de la matriz identidad. Ejemplo En el ejemplo anterior tenemos que A = E12 , B = E21 (2) y C = E2 (3).
Ejemplo Para ilustrar la de las operaciones elementales consideremos invertibilidad 0 1 3 la matriz A = 2 0 2 3 1 0 Teorema Toda matriz elemental es invertible y las inversas estan dadas por
1 Eij = Eij ,
Eij (c)1 = Eij (c)
Ei (c)1 = Ei
1 c
Ejemplo Para ilustrar la de las operaciones elementales consideremos invertibilidad 0 1 3 la matriz A = 2 0 2 3 1 0 Teorema Toda matriz elemental es invertible y las inversas estan dadas por
1 Eij = Eij ,
Eij (c)1 = Eij (c)
Ei (c)1 = Ei
1 c
Teorema Sea A una matriz de tama no m n y E una matriz elemental de tama no m m asociada a una operaci on elemental de la, el producto EA es la matriz que se obtiene la aplicar la operaci on elemental de la a la matriz A.
Ejemplo Para ilustrar la de las operaciones elementales consideremos invertibilidad 0 1 3 la matriz A = 2 0 2 3 1 0 Teorema Toda matriz elemental es invertible y las inversas estan dadas por
1 Eij = Eij ,
Eij (c)1 = Eij (c)
Ei (c)1 = Ei
1 c
Teorema Sea A una matriz de tama no m n y E una matriz elemental de tama no m m asociada a una operaci on elemental de la, el producto EA es la matriz que se obtiene la aplicar la operaci on elemental de la a la matriz A.
Corolario Toda matriz se puede expresar como el producto de un n umero nito de matrices elementales por una matriz en forma escalonada reducida. Mas concretamente, si A una matriz de tama no m n, existen matrices elementales E1 , . . . , Ek todas de tama no m m y una matriz escalonada reducida A tal que A = E1 Ek A . Ejemplo 1 Expresar la matriz A = 1 2 elementales por una matriz 1 0 1 2 como un producto de matrices 2 2 en forma escalonada reducida.
Corolario Toda matriz se puede expresar como el producto de un n umero nito de matrices elementales por una matriz en forma escalonada reducida. Mas concretamente, si A una matriz de tama no m n, existen matrices elementales E1 , . . . , Ek todas de tama no m m y una matriz escalonada reducida A tal que A = E1 Ek A . Ejemplo 1 Expresar la matriz A = 1 2 elementales por una matriz 1 0 1 2 como un producto de matrices 2 2 en forma escalonada reducida.
Inplementaci on al Matlab Consideremos las siguientes matrices elementales: >> E 1=[1 0 0;1 1 0;0 0 1]; E 2=[1 0 0;0 1 0;2 0 1]; E 3=[1 0 0;0 2 0;0 0 1]; E 4=[1 0 0;0 1 0;0 2 1] La matriz escalonada >> A =[1 1 0;00 1;0 0 0] >> A = E 1 E 2 E 3 E 4 A . Obtendremos la matriz de el ejemplo anterior. Lema Sea A una matriz de tama no n n en forma escalonada reducida, entonces A es invertible si y solo si A = I .
Inplementaci on al Matlab Consideremos las siguientes matrices elementales: >> E 1=[1 0 0;1 1 0;0 0 1]; E 2=[1 0 0;0 1 0;2 0 1]; E 3=[1 0 0;0 2 0;0 0 1]; E 4=[1 0 0;0 1 0;0 2 1] La matriz escalonada >> A =[1 1 0;00 1;0 0 0] >> A = E 1 E 2 E 3 E 4 A . Obtendremos la matriz de el ejemplo anterior. Lema Sea A una matriz de tama no n n en forma escalonada reducida, entonces A es invertible si y solo si A = I . Teorema Sea A una matriz de tama no n n, entonces A es invertible si y solo si A se puede escribir como un producto de matrices elementales.
Inplementaci on al Matlab Consideremos las siguientes matrices elementales: >> E 1=[1 0 0;1 1 0;0 0 1]; E 2=[1 0 0;0 1 0;2 0 1]; E 3=[1 0 0;0 2 0;0 0 1]; E 4=[1 0 0;0 1 0;0 2 1] La matriz escalonada >> A =[1 1 0;00 1;0 0 0] >> A = E 1 E 2 E 3 E 4 A . Obtendremos la matriz de el ejemplo anterior. Lema Sea A una matriz de tama no n n en forma escalonada reducida, entonces A es invertible si y solo si A = I . Teorema Sea A una matriz de tama no n n, entonces A es invertible si y solo si A se puede escribir como un producto de matrices elementales.
Ejemplo 0 Expresar la matriz A = 1 1 matrices elementales. 1 1 2 2 1 y su inversa como un producto de 2
Teorema (Algoritmo para calcular A1 ) Sea A una matriz invertible de tama no n n, al aplicar reducci on Gauss-Jordan a la matriz aunmentada A I obtenemos la matriz I A1 . M as a un, si B es una matriz de tama no n q , al aplicar reducci on Gauss-Jordan a la matriz aumentada A B obtenemos la matriz I A1 B .
Ejemplo 0 Expresar la matriz A = 1 1 matrices elementales. 1 1 2 2 1 y su inversa como un producto de 2
Teorema (Algoritmo para calcular A1 ) Sea A una matriz invertible de tama no n n, al aplicar reducci on Gauss-Jordan a la matriz aunmentada A I obtenemos la matriz I A1 . M as a un, si B es una matriz de tama no n q , al aplicar reducci on Gauss-Jordan a la matriz aumentada A B obtenemos la matriz I A1 B .
Inplementaci on al Matlab 0 1 2 1 1 calcular A1 . Para hallar la inversa de la matriz Sea A = 1 1 2 2 A en el Matlab se da por: >> AI = [0 1 2 1 0 0; 1 1 1 0 1 0; 1 2 2 0 0 1]; rref (AI ) Inplementaci on al Matlab 0 1 2 1 1 1 y B = 2 Sean A = 1 1 2 2 0 en Matlab : >> AB = [0 1 2 1 2; 1 1 1 2 3; 1 2
2 3, calcular A1 B . Inplementando 1 2 0 1]; rref (AB )
Inplementaci on al Matlab 0 1 2 1 1 calcular A1 . Para hallar la inversa de la matriz Sea A = 1 1 2 2 A en el Matlab se da por: >> AI = [0 1 2 1 0 0; 1 1 1 0 1 0; 1 2 2 0 0 1]; rref (AI ) Inplementaci on al Matlab 0 1 2 1 1 1 y B = 2 Sean A = 1 1 2 2 0 en Matlab : >> AB = [0 1 2 1 2; 1 1 1 2 3; 1 2
2 3, calcular A1 B . Inplementando 1 2 0 1]; rref (AB )
Observaciones Si en el teorema anterior las columnas de la matriz B son los vectores b1 , . . . , bq , entonces al aplicar reducci on Gauss-Jordan a la matriz aumentada A B = A b1 bq obtenemos la matriz la matriz A1 B = A1 b1 A1 bq , obteniendo soluciones simultaneas a los sistemas Ax = b1 , . . . , Ax = bq . Ejemplo Usar la observaci on anterior para resolver los sistemas x2 + 2x3 = 1 x1 + x2 + x3 = 2 x1 2x2 2x3 = 0 y x2 + 2x3 = 2 x1 + x2 + x3 = 3 . x1 2x2 2x3 = 1
Observaciones Si en el teorema anterior las columnas de la matriz B son los vectores b1 , . . . , bq , entonces al aplicar reducci on Gauss-Jordan a la matriz aumentada A B = A b1 bq obtenemos la matriz la matriz A1 B = A1 b1 A1 bq , obteniendo soluciones simultaneas a los sistemas Ax = b1 , . . . , Ax = bq . Ejemplo Usar la observaci on anterior para resolver los sistemas x2 + 2x3 = 1 x1 + x2 + x3 = 2 x1 2x2 2x3 = 0 y x2 + 2x3 = 2 x1 + x2 + x3 = 3 . x1 2x2 2x3 = 1
Definici on Sea A una matriz de tama no m n, 1. decimos que A tiene inversa a la izquierda si existe una matriz L de tama no n m tal que LA = In , 2. decimos que A tiene inversa a la derecha si existe una matriz R de tama no n m tal que AR = Im . Teorema Sea A una matriz de tama no m n, entonces las siguientes armaciones son equivalentes: 1. A tiene inversa a la derecha. 2. El sistema Ax = b tiene soluci on para cada b Rm . 3. rango(A) = m = # de las de A, 4. La funci on TA : Rn Rm denida por TA (x) = Ax es sobreyectiva.
Definici on Sea A una matriz de tama no m n, 1. decimos que A tiene inversa a la izquierda si existe una matriz L de tama no n m tal que LA = In , 2. decimos que A tiene inversa a la derecha si existe una matriz R de tama no n m tal que AR = Im . Teorema Sea A una matriz de tama no m n, entonces las siguientes armaciones son equivalentes: 1. A tiene inversa a la derecha. 2. El sistema Ax = b tiene soluci on para cada b Rm . 3. rango(A) = m = # de las de A, 4. La funci on TA : Rn Rm denida por TA (x) = Ax es sobreyectiva.
Ejemplo Determine si la matriz A = 1 2 2 tiene inversa a la derecha y en 1 0 4 caso armativo calcular una inversa a la derecha de A.
Teorema Sea A una matriz de tama no m n, entonces las siguientes armaciones son equivalentes: 1. A tiene inversa a la izquierda. 2. El sistema Ax = b tiene a lo sumo una soluci on para cada b Rm . 3. La funci on TA : Rn Rm denida por TA (x) = Ax es inyectiva. 4. El sistema Ax = n tiene soluci on u nica. 5. rango(A) = n = # de columnas de A.
Ejemplo Determine si la matriz A = 1 2 2 tiene inversa a la derecha y en 1 0 4 caso armativo calcular una inversa a la derecha de A.
Teorema Sea A una matriz de tama no m n, entonces las siguientes armaciones son equivalentes: 1. A tiene inversa a la izquierda. 2. El sistema Ax = b tiene a lo sumo una soluci on para cada b Rm . 3. La funci on TA : Rn Rm denida por TA (x) = Ax es inyectiva. 4. El sistema Ax = n tiene soluci on u nica. 5. rango(A) = n = # de columnas de A.
Observaciones El teorema provee un algorimo para calcular dicha matriz, el cual describimos a continuaci on 1. Se aplica reducci on Gauss-Jordan a la matriz 2. Se calculan las matrices elementales E1 , . . . , Ek asociadas a cada una de las operaciones elementales aplicadas en el paso 1. 3. Se calcula el producto Ek E1 . 4. La matriz L formada por las primeras n las de la matriz Ek E1 es la inversa a la izquierda. Inplementaci on al Matlab Determine si las matrices tienen inversa a la izquierda y encuentrela en caso armativo. 2 1 1 1 1 1 0 2 1 . a. A = 4 b. B = c. C = 0 2 1 1 2 1 2 1
Observaciones El teorema provee un algorimo para calcular dicha matriz, el cual describimos a continuaci on 1. Se aplica reducci on Gauss-Jordan a la matriz 2. Se calculan las matrices elementales E1 , . . . , Ek asociadas a cada una de las operaciones elementales aplicadas en el paso 1. 3. Se calcula el producto Ek E1 . 4. La matriz L formada por las primeras n las de la matriz Ek E1 es la inversa a la izquierda. Inplementaci on al Matlab Determine si las matrices tienen inversa a la izquierda y encuentrela en caso armativo. 2 1 1 1 1 1 0 2 1 . a. A = 4 b. B = c. C = 0 2 1 1 2 1 2 1
Inplementaci on al Matlab 1 1 1 . Calcular una inversa a la izquierda de C = 0 2 1 >> C t = [1, 0, 2; 1, 1, 1], e1 = [1; 0], e2 = [0; 1], r1 = Ae1, r2 = Ae2
Corollary (Caracterizaci on de una matriz invertible) Sea A una matriz cuadrada de tama no n n, entonces las siguientes armaciones son equivalentes: 1. A es invertible. 2. A tiene inversa a la izquierda. 3. El sistema Ax = b tiene a lo sumo una soluci on para cada b Rm . 4. La funci on TA : Rn Rm denida por TA (x) = Ax es inyectiva. 5. El sistema Ax = n tiene soluci on u nica. 6. rango(A) = n = # de columnas de A=# de las deA. 7. A tiene inversa a la derecha. 8. El sistema Ax = b tiene soluci on para cada b Rm . 9. rango(A) = m = # de las de A, 10. La funci on TA : Rn Rm denida por TA (x) = Ax es sobreyectiva. 11. La funci on TA : Rn Rn denida por TA (x) = Ax es biyectiva. 12. At es invertible. 13. A es un producto de matrices elementales.
