M. y Det
M. y Det
M. y Det
MATRICES Y DETERMINANTES
6.1. Introduccion
Las matrices y los determinantes son herramientas del algebra que facilitan el ordenamiento de datos, as como su manejo. Los conceptos de matriz y todos los relacionados fueron desarrollados basicamente en el siglo XIX por matematicos como los ingleses J.J. Sylvester y Arthur Cayley y el irlandes William Hamilton. Las matrices se encuentran en aquellos ambitos en los que se trabaja con datos regularmente ordenados y aparecen en situaciones propias de las Ciencias Sociales , Economicas y Biologicas.
Filas de la matriz
a11 a21 . . .
a12 a22 . . .
a13 a23 . . .
... ... .. .
a1n a2n . . .
amn
am 1 am 2 am 3 ...
Columnas de la matriz A
Abreviadamente se puede expresar A = (aij). Cada elemento de la matriz lleva dos subndices. El primero de ellos i, indica la fila en la que se encuentra el elemento, y el segundo, j, la columna. As el elemento a23 esta en la fila 2 y columna 3. Las matrices siempre se representaran con letras mayusculas. Ejemplos: Son ejemplos de matrices los siguientes:
A tiene 2 filas y 2 columnas, diremos que su tamano es 2 x 2.Que elemento es a21?. B tiene 2 filas y 3 columnas, diremos que su tamano es 2 x 3.Que elemento es b23?. C tiene 4 filas y 3 columnas, diremos que su tamano es 4 x 3.Que elemento es c42?. En general, si una matriz A tiene m filas y n columnas, diremos que su tamano o d imension es m x n (se lee m por n), siempre en primer lugar el n de filas y en segundo lugar el de columnas.
82
es una matriz nula de tamano 2 x 5. 2. Se llama matriz fila a la que solo tiene una fila, es decir su dimension es 1x n. Por ejemplo, es una matriz fila de tamano 1 x 4. 3. Se llama matriz columna a la que solo consta de una columna, es decir su dimension sera m x 1, como por ejemplo:
x 1.
Una matriz es cuadrada cuando tiene el mismo numero de filas que de columnas, es decir su dimension es n x n. La matriz ( 23 14 ) del primer ejemplo anterior es cuadrada de tamano 2 x 2 o simplemente de orden 2. Otro ejemplo de matriz cuadrada es:
de orden 3. Dentro de las matrices cuadradas llamaremos diagonal principal a la formada por los elementos a11,a22,a33,...,ann, siendo la matriz:
En la matriz D del ejemplo anterior, su diagonal principal estara formada por 1, 5, 0. Se llama traza de la matriz a la suma de los elementos de la diagonal. Es decir, Traza (A)=a11 + a22 + a33 + ... + ann , y en el caso de D, Traza (D)= 1+5+0 = 6. La diagonal secundaria es la formada por los elementos a1n,a2,n1,a3,n 2,...,an1. En la matriz D estara formada por 3, 5, -3. Una clase especial de matrices cuadradas son las matrices triangulares. Una matriz es triangular superior si todos los elementos por debajo de la diagonal principal son nulos y triangular inferior si son nulos todos los elementos situados por encima de dicha diagonal. Son ejemplos de estas matrices:
Triangular inferior
Si una matriz es a la vez triangular superior e inferior, solo tiene elementos en la diagonal principal. Una matriz de este tipo se denomina matriz diagonal. Un ejemplo de matriz diagonal sera: 1 G 0 0 0 45 0 0 3 0 0 0 0 0 0 0 Por ultimo, si una matriz diagonal tiene en su diagonal principal solo unos, se denomina matriz unidad o identidad. Se suelen representar por In , donde n es el orden o tamano de la matriz. Algunas matrices identidad son:
ud
Estas matrices se denominan matrices de informacion, y simplemente recogen los datos numericos del problema en cuestion. Otras matrices son las llamadas matrices de relacion, que indican si ciertos elementos estan o no relacionados entre s. En general, la existencia de relacion se expresa con un 1 en la matriz y la ausencia de dicha relacion de expresa con un 0. Estas matrices se utilizan cuando queremos trasladar la informacion dada por un grafo y expresarla numericamente. En Matematicas, un grafo es una coleccion cualquiera de puntos conectados por lineas. Existen muchos tipos de grafos. Entre ellos, podemos destacar: * Grafo simple: Es el grafo que no contiene ciclos, es decir, lineas que unan un punto consigo mismo, ni lineas paralelas, es decir, lineas que conectan el mismo par de puntos. * Grafo dirigido: Es el grafo que indica un sentido de recorrido de cada linea, mediante una flecha. Estos tipos de grafo pueden verse en la figura:
Figura 6.1: Grafo, Grafo simple y Grafo dirigido. Relacionadas con los grafos se pueden definir algunas matrices. Entre todas ellas, nosotros nos fijaremos en la llamada matriz de adyacencia, que es aquella formada por ceros y unos exclusivamente, de tal forma que: * un 1 en el lugar (i,j) expresa la posibilidad de ir desde el punto de la fila i hasta el punto de la columna j mediante una linea que los una directamente. * un 0 en el lugar (i,j) expresa la imposibilidad de ir del primer punto al segundo mediante una linea que los una directamente. La matriz de adyacencia del grafo dirigido de la figura anterior sera: A A 0 B 0 C D0 Ejercicio 1) Escribe las correspondientes matrices de adyacencia de los grafos: B C D 1 0 0 1 0 0 10 1 0 0 00
A A 0 B D0
B 0C 1
1 1 11 0 0 00 0 0 1 1 2) Dibuja los grafos dirigidos que correspondan a las matrices de adyacencia: A B C A B 101 C 000
Si las matrices tienen diferente tamano, no se pueden sumar o restar entre s. Propiedades de la suma (y diferencia) de matrices: a) Conmutativa: A + B = B + A b) Asociativa: A + (B + C) = (A + B) + C c) Elemento neutro: La matriz nula del tamano correspondiente. d) Elemento opuesto de A: La matriz -A, que resulta de cambiar de signo a los elementos de A. Ejemplo: Si porque:
0 2 + 3 9
3x2
0 4 3
3x2
2 9
3x2
= 0
0 0
Ejercicios: 1. Las exportaciones, en millones de euros, de 3 pases A, B, C a otros tres X, Y, Z, en los anos 2000 y 2001 vienen dadas por las matrices:
Calcula y expresa en forma de matriz el total de exportaciones para el conjunto de los dos anos. Cuantos millones ha exportado el pas B al Z en total? Calcula el incremento de las exportaciones del ano 2000 al 2001 con los datos del ejemplo anterior. 2. Calcula x, y, z en la suma:
3.
Propiedades: a) Distributiva respecto de la suma de matrices: k (A + B) = k A + k B b) Distributiva respecto de la suma de numeros: (k + d) A= k A + d A c) Asociativa: k (d A)=(k d) A d) Elemento neutro, el numero 1: 1 A=A Ejercicios: 1. Si , halla una matriz X que verifique la ecuacion: 2X4A=B 2. Determina las matrices X y Y sabiendo que:
At 0 7 2 1
Evidentemente, si A es una matriz de tamano m x n, su traspuesta At tendra tamano n x m, pues el numero de columnas pasa a ser el de filas y viceversa. Si la matriz A es cuadrada, su traspuesta tendra el mismo tamano. Propiedades: a) (At)t = A, es decir, la traspuesta de la traspuesta es la matriz inicial. b) (A + B)t = At + Bt c) (k A)t = k At En base a esta nueva operacion, podemos definir otras dos clases de matrices, que son: Matriz simetrica, que es aquella para la que se cumple que At = A, por ejemplo la matriz:
es simetrica (compruebalo). En una matriz simetrica, los elementos son simetricos respecto a la diagonal principal. Ejercicio: Puede ser simetrica una matriz que no sea cuadrada?Por que?. Matriz antisimetrica, es aquella para la que se cumple que At = A. Por ejemplo:
es antisimetrica ( comprueba ). En una matriz antisimetrica, los elementos de la diagonal principal son siempre nulos (por que?), y los restantes son opuestos respecto a dicha diagonal. Ejercicios: 1. 2. Dadas las matrices Obtener las matrices X e Y que verifiquen los sistemas: calcula 3At Bt.
primero comprobamos que se puede realizar el producto A B, pues el n de columnas de A es 4 y el n de filas de B tambien es 4, y el resultado, segun lo dicho sera una matriz de tamano 2 x 3, tiene 2 filas y 3 columnas:
Solo nos falta completar los elementos de la matriz producto. Para ello, seguimos la regla anterior: El elemento de la fila 1 y columna 1 de A B proviene de multiplicar elemento a elemento la fila 1 de A por la columna 1 de B y sumar, es decir: (3) 0 + 2 1 + 1 2 + 4 3 = 0 + 2 + 2 + 12 = 16 El elemento de la fila 1 y columna 2 de A B proviene de multiplicar elemento a elemento la fila 1 de A y la columna 2 de B y sumar: (3) (4) + 2 (2) + 1 0 + 4 2 = 12 4 + 0 + 8 = 16 El elemento de la fila 1 y columna 3 de A B proviene de multiplicar elemento a elemento la fila 1 de A y la columna 3 de B y sumar: (3) 1 + 2 1 + 1 2 + 4 1 = 3 + 2 + 2 + 4 = 5
3. 4.
Lo mismo si
calcula: A + B,3A 4B,A B,A D,B C,C D,At C,Dt At,Bt A,Dt D,D Dt Propiedades del producto de matrices a) Asociativa: A (B C) = (A B) C b) Distributiva respecto de la suma: A (B + C) = A B + A C B + C) A = B A + C A c) Elemento neutro, la matriz identidad correpondiente, si A es m x n: A In = A Im A = A d) En general el producto de matrices no es conmutativo
Pueden verse ejemplos en los ejercicios anteriores. Esta es una propiedad muy importante. e) El producto de dos matrices no nulas A y B puede dar lugar a una matriz nula:
CAP ITULO 6. MATRICES Y DETERMINANTES 91 Se dice que el conjunto de las matrices con la operacion producto tiene divisores de cero, es decir, hay matrices no nulas cuyo producto es nulo. Ejercicios: 1. Si A y B son dos matrices cuadradas del mismo orden, son ciertas las propiedades siguientes, que son ciertas para las operaciones con numeros reales?: a) (A + B)2 = A2 + B2 + 2 A B b) (A B)2 = A2 + B2 2 A B c) (A + B) (A B) = A2 B2 2. 3. Determina los valores de a y b de la matriz Que matrices conmutan con la matriz ?. para que A2 = A.
CAP ITULO 6. MATRICES Y DETERMINANTES 92 2) No todas las matrices cuadradas no nulas tienen matriz inversa (sea lo que sea, por analoga con los numeros). Definamos, en primer lugar, el termino de matriz inversa: Dada una matriz cuadrada de orden n , A, se dice que A es invertible (o que posee inversa o que es no singular o que es regular ), si existe otra matriz del mismo orden, denominada matriz inversa de A y representada por A1 y tal que: A A1 = In y A1 A = In Si A no tiene inversa, se dice que es singular o no invertible. Si una matriz tiene inversa, dicha matriz inversa es unica (solo hay una). Para calcular dicha matriz inversa, podemos utilizar dos vas:
Es decir, hemos de resolver un sistema de 4 ecuaciones con 4 incognitas, aunque en realidad son 2 sistemas de dos ingonitas cada uno (uno con x y z y otro con y y t). Resolviendo el sistema se obtiene que
Se puede comprobar que tambien se cumple que A1 A = I2, luego A es invertible, tiene inversa. Si el sistema no tiene solucion, la matriz no tiene inversa. Por ejemplo, en el caso en que , del mismo modo :
Y por ejemplo de 2x+2z=0 se obtiene x = -z, si se sustituye en la primera ecuacion es -z+z=1, es decir 0 = 1 (imposible). El sistema no tiene solucion. Por tanto A no es invertible, es singular. Este metodo directo solo se suele utilizar para matrices cuadradas de tamano 2, puesto que para las de tamano 3 obtenemos un sistemas de 9 ecuaciones con 9 incognitas! que realmente es difcil de resolver.
ii) Se hace la matriztriangular superior (es decir, hacemos ceros por debajo de la diagonal principal) usando transformaciones elementales en filas. La mejor forma de realizar esto es hacer cero los elementos por debajo de la diagonal en la primera columna usando la fila 1. Luego, hacer cero los elementos por debajo de la diagonal en la segunda columna usando la fila 2, y as sucesivamente. En nuestro caso, basta sumar la fila 2 con la fila 1, y se obtiene:
iii) Una vez hecha la matriz triangular superior, se hace la matriz triangular inferior, haciendo ceros a los elementos por encima de la diagonal. El proceso es parecido al anterior: Hacer cero los elementos por encima de la diagonal en la ultima columna usando la ultima fila. Luego, hacer cero los elementos por encima de la diagonal en la penultima columna usando la penumtima fila, y as sucesivamente. En nuestro caso:
iv) Ya tenemos una matriz diagonal. Lo unico que falta es dividir a cada fila entre el numero adecuado para obtener unos en la diagonal principal, es decir, para obtener la matriz identidad en la parte izquierda:
v) Una vez se tiene la matriz identidad en la parte de la izquierda, la parte derecha es la matriz inversa, es decir, llegamos a:
matriz que habamos obtenido antes por el metodo directo. Si al realizar el metodo de Gauss-Jordan en algun momento alguna fila es de ceros, la matriz no tiene inversa. Cuanto mayor sea el orden de la matriz, mejor es este metodo frente al directo. Veamos otro ejemplo: Calcular la inversa de la matriz Siguiendo los pasos anteriores: por el metodo de Gauss-Jordan.
Como aparece una fila de ceros, la matriz A no tiene inversa. Ejercicios: 1. Calcular por el metodo de Gauss-Jordan la inversa de las matrices:
2.
CAP ITULO 6. MATRICES Y DETERMINANTES 95 rapida la inversa de una matriz diagonal cualquiera?.
a) b)
, Rg(A)=1 ,solo una fila distinta de cero. , Rg(B)=2 hay 2 filas no nulas. c)
Los ejemplos anteriores ponen de manifiesto que el rango de cualquier matriz siempre es menor o igual que el numero de filas de la matriz. De hecho se verifica que el rango de cualquier matriz siempre es menor o igual que su numero de filas y de columnas, pues el proceso para hacer el metodo de Gauss se puede hacer indistintamente mediante operaciones elementales en filas o en columnas. Esto permite, antes de calcular el rango de una matriz, saber entre que valores va a estar ese rango. Por ejemplo, en el caso c) del ejemplo, como la matriz es 3x3 , el rango solo puede ser 0, 1, 2 o 3 , no hay otras posibilidades.
CAP ITULO 6. MATRICES Y DETERMINANTES 96 En el caso del apartado d), como la matriz es 2 x 3, el rango solo puede ser 0,1 o 2. (De hecho, podemos reducir esto algo mas , pues una matriz solo tiene rango cero si es la matriz nula).Resumiendo: Propiedad: Si A es una matriz de tamano m x n no nula se cumple que: 1 Rg(A) min{m,n} Ejemplo: Calcular en funcion de k el rango de la matriz:
Aplicando Gauss,
Ahora es evidente que si k-6=0, la ultima fila es nula. Por tanto, si k=6, la ultima fila es nula y el rango de A es 1, Rg(A)=1, mientras que si k-6 es distinto de cero, es decir si k es distinto de 6, hay 2 filas no nulas y el rango de A es 2, Rg(A)=2. Resumiendo: Si k = 6, entonces Rg(A)=2 Si k=6, entonces Rg(A)=1 La siguiente propiedad permite relacionar el concepto de rango con el de matriz inversa visto anteriormente: Propiedad: Una matriz cuadrada A tiene inversa Rg(A) es maximo. Ejercicios: 1. Calcula el rango de A segun los valores de k: inversa?. 2. Calcula el rango de las matrices: .Para que valores de k tiene A
6.8. Determinantes
Introduciremos a continuacion el concepto de determinante asociado a una matriz cuadrada. Este concepto permite simplificar operaciones matriciales tales como el calculo del rango o de la matriz inversa. Definicion:
CAP ITULO 6. MATRICES Y DETERMINANTES 97 Si es una matriz 2 x 2 se define el determinante de la matriz A, y se expresa como det(A) o bien | A|, como el numero:
b) Para definir determinantes de matrices de orden mayor que 2 es necesario introducir previamente algunos conceptos. Dada una matriz cuadrada A de orden n, definimos el menor complementario de un elemento de A,aij , como el determinante de la matriz que se obtiene al suprimir la fila i y la columna j en la que se encuentra dicho elemento aij. Se representa por Mij. Ejemplo: En la matriz elementos de la primera fila son: Menor complementario de -2: , los menores complementarios de cada uno de los
Menor complementario de 4: Menor complementario de 5: Y as sucesivamente. Ejercicio: Obtener los restantes menores complementarios de los elementos de la matriz A. Estrechamente ligado al concepto de menor complementario se encuentra el de adjunto de una matriz. Dada una matriz cuadrada A de orden n, definimos el adjunto de un elemento aij de A como el numero: Aij = (1)i+j Mij es decir, no es mas que el menor complementario correspondiente acompanado de un signo mas o menos dependiendo de la fila y la columna en la que se encuentre el elemento en cuestion. Por ejemplo, para la matriz anterior, los adjuntos de los elementos de la primera fila son: Adjunto de -2:A11 = (1)1+1 M11 = 1 14 = 14 (coincide con el menor complementario) Adjunto de 4:A12 = (1)1+2 M12 = (1)21 = 21 (menor complementario con signo cambiado) Adjunto de 5:A13 = (1)1+3 M13 = 1 21 = 21 (coincide con el menor complementario). Ejercicio: Obtener los restantes adjuntos de los elementos de la matriz A.
CAP ITULO 6. MATRICES Y DETERMINANTES 98 En general puede saberse si el signo del menor complementario y del adjunto coinciden o no utilizando una sencilla regla grafica, por ejemplo, para matrices 3 x 3 y 4 x 4 basta fijarse en las matrices:
donde el + significa que el adjunto coincide con el menor complementario y el - indica que tienen signo contrario. Una vez vistos estos conceptos se puede definir ya: Definicion: Dada una matriz cuadrada A de tamano n se define su determinante como la suma del producto de los elementos de una linea cualquiera de la matriz (fila o columna) elegida, por sus correpondientes adjuntos. Se puede demostrar, aunque dicha demostracion excede los contenidos del curso, que el valor del determinante no depende de la fila o columna elegida para calcularlo. Ejemplo: Para la matriz ,aplicando la definicion, si elegimos la fila tercera queda:
Ejercicio: Calcula, desarrollando por la fila que tu elijas los determinantes de las matrices: 1 1 1 1 1 1 1 1 0 1 0 1 1 1 0 3 4 1 3 2 5 1 3 5 7 1 0 8 7 0 0 3 20 3 4
21 23
CAP ITULO 6. MATRICES Y DETERMINANTES 99 En el caso de las matrices cuadradas de orden 3, esta regla facilita el calculo de dichos determinantes. Si la matriz es , entonces el determinante de A se calcula mediante la resta
de dos expresiones obtenidas del siguiente modo: Llamaremos sumandos positivos a los obtenidos al multiplicar: - Los elementos de la diagonal principal,a11 a22 a33. - Los elementos de la linea paralela superior a la diagonal principal por el elemento aislado de la esquina inferior izquierda:a12 a23 a31. - Los elementos de la linea paralela inferior a la diagonal principal por el elemento aislado de la esquina superior derecha:a21 a32 a13. Graficamente:
Llamaremos sumandos negativos a los obtenidos al multiplicar: - Los elementos de la diagonal secundaria,a13 a22 a31. - Los elementos de la linea paralela superior a la diagonal secundaria por el elemento aislado de la esquina inferior derecha: a12 a21 a33. - Los elementos de la linea paralela inferior a la diagonal secundaria por el elemento aislado de la esquina superior izquierda: a32 a23 a11. Graficamente:
Y entonces det (A)= Sumandos positivos - Sumandos negativos. Por ejemplo, en el caso de la matriz anterior:
CAP ITULO 6. MATRICES Y DETERMINANTES 100 , se tiene que aplicando la regla de Sarrus: det(A)=(-2) 7 2+4 3 (-3)+6 5 0-(3 7 5+0 (-2) (-3)+6 4 2)=-28-36-105-48=-217. Ejercicio: Comprobar, mediante la regla de Sarrus, los determinantes de orden 3 obtenidos en el ejercicio anterior.
4.
Si multiplicamos todos los elementos de una linea de un determinante por un numero, el determinante queda multiplicado por ese numero. Por ejemplo:
pero
5.
Si a una linea de una matriz se le suma otra linea multiplicada por un numero, el determinante no cambia. Esta propiedad permite utilizar un metodo mas sencillo para calcular determinantes de orden mayor que 3.
6.
7.
CAP ITULO 6. MATRICES Y DETERMINANTES 101 Una estrategia a tener en cuenta en este caso de determinantes de orden 4 o superior, o incluso de orden 3 si la matriz es compleja, es el metodo de hacer ceros, puesto que el valor del determinante no vara al realizar a la matriz ciertas transformaciones elementales en filas,como indica la propiedad 5 anterior, si bien hemos de ser cuidadosos al aplicar dicha propiedad. As pues la mejor forma de calcular un determinante es hacer ceros en una fila o columna y desarrollar por dicha fila o columna, porque entonces solo tendremos que calcular un adjunto. Por ejemplo, si calculamos:
Como hemos dicho, hemos de tener especial cuidado al aplicar esta regla con determinantes, puesto que no podemos hacer las mismas operaciones que con las matrices, lo que puede confundir. Por ejemplo, si queremos calcular el determinante: mediante
Si hiciesemos ceros en la primera columna, y desarrollasemos nos debera dar lo mismo. Ahora bien,podemos hacer cero el 4 de la primera columna mediante: .
no obtenemos lo mismo, porque hemos multiplicado la fila sustituida por un numero y eso altera el valor del determinante. Luego la fila a sustituir conviene no multiplicarla, como en el primer ejemplo, puesto que si no nos damos cuenta, podemos variar el valor del determinante.
donde Adj(A) denota la matriz adjunta de A, es decir, aquella que se obtiene de sustituir cada elemento de A por su adjunto. Ejemplo: Calcular, si es posible, la inversa de la matriz .
3 1 1 1 1
32 21 1 1 3
Por tanto,
(Adj(A
Y entonces, se obtiene:
A =2
2 21 12
2 1
Como es nulo, el rango de la matriz NO es 2. Menores de orden 1 hay 4, por ejemplo |1| = 1, que es no nulo,luego el rango de la matriz es Rg(A)=1 (el tamano de dicho menor complementario). b) Solo hay un menor de orden 2, que es:
Como no es nulo, el rango de la matriz es complementario). c) Solo hay un menor de orden 3, que es:
Como es nulo, podemos asegurar que el rango NO es 3. Menores de orden 2 hay 9. Calculando alguno:
resulta que es no nulo, luego el rango es Rg(C)=2 (el tamano de dicho menor complementario). d) El menor mas grande que podemos formar es de orden 2. Hay 3 de ellos:
Son todos nulos, luego el rango NO es 2. Menores de orden 1 hay 6, y por ejemplo nulo, luego el rango es Rg(D)=1. Ejercicio Calcula,utilizando los determinantes, el rango de las matrices:
= 0, es no