4 Determinantes de matrices cuadradas

Definición 4.1 (Determinante).

Por inducción en el tamaño de A𝔐n(𝕂)\displaystyle A\in\mathfrak{M}_{n}(\mathbb{K}).

  •  

    Si a=1\displaystyle a=1: det((a11))=a11\displaystyle\det((a_{11}))=a_{11}.

  •  

    Supongamos que sabemos calcular determinantes de matrices de orden n1\displaystyle n-1 y sea A=(aij𝔐n(𝕂))\displaystyle A=(a_{ij}\in\mathfrak{M}_{n}(\mathbb{K})). Llamamos (ij)esimo\displaystyle(ij)-esimo adjunto de la matriz A\displaystyle A al numero

    αij=(1)i+1det(Aij)\displaystyle\alpha_{ij}=(-1)^{i+1}\det(A_{ij})

    donde Aij\displaystyle A_{ij} es la matriz de orden n1\displaystyle n-1 obtenida eliminando la fila i\displaystyle i y la columna j\displaystyle j de A\displaystyle A. Definimos el determinante de A\displaystyle A como

    det(A)=a11α11+a21α21++an1αn1\displaystyle\det(A)=a_{11}\alpha_{11}+a_{21}\alpha_{21}+\cdots+a_{n1}\alpha_{% n1}

    A esta formula se le llama desarrollo del determinante por adjuntos de la primera columna.

Lema 4.1.

El determinante de la matriz identidad es igual a 1\displaystyle 1.

Demostración.

Por induccion en el orden n\displaystyle n de la matriz identidad.

  •  

    Si n=1\displaystyle n=1, det(I1)=det(1)=1\displaystyle\det(I_{1})=\det(1)=1.

  •  

    Supongamos que det(In1)=1\displaystyle\det(I_{n-1})=1 y veamos que det(In)=1\displaystyle\det(I_{n})=1.

    Por una parte, el adjunto (11)\displaystyle(11)-esimo de In\displaystyle I_{n} es α11=(1)1+1det(In1)\displaystyle\alpha_{11}=(-1)^{1+1}\det(I_{n-1}), asi que por la hipotesis de induccion, α11=1\displaystyle\alpha_{11}=1.

    Ademas, las entradas (21),(31),,(n1)\displaystyle(21),(31),\ldots,(n1) de la matriz identidad son todas iguales a cero, asi que no hace falta calcular α21,,αn1\displaystyle\alpha_{21},\ldots,\alpha_{n1}. Sustituyendo en la formula del desarrollo del determinante por adjuntos de la primera columna

    det(In)=1α11+0α21++0αn1=1.\displaystyle\det(I_{n})=1\alpha_{11}+0\alpha_{21}+\cdots+0\alpha_{n1}=1.

Proposición 4.1 (Propiedades de los determinantes).

  1. 1.

    Si A,A,A′′\displaystyle A,A^{\prime},A^{\prime\prime} son tres matrices identicas de orden n\displaystyle n salvo en que la fila i\displaystyle i de A\displaystyle A es la suma de la fila i\displaystyle i de A\displaystyle A^{\prime} y la fila i\displaystyle i de A′′\displaystyle A^{\prime\prime} entonces

    det(A)=det(A)+det(A′′)\displaystyle\det(A)=\det(A^{\prime})+\det(A^{\prime\prime})
  2. 2.

    Si dos filas de A\displaystyle A son iguales, entonces det(A)=0\displaystyle\det(A)=0.

  3. 3.

    Si se intercambian dos filas de A\displaystyle A, entonces el determinante cambia de signo.

  4. 4.

    Si multiplicamos una fila de A\displaystyle A por un escalar λK\displaystyle\lambda\in K entonces el determinante de la matriz obtenida es igual a λdet(A)\displaystyle\lambda\det(A).

  5. 5.

    El determinante de A\displaystyle A no varia si a una fila de A\displaystyle A le sumamos otra fila multiplicada por un escalar.

  6. 6.

    A\displaystyle A es invertible si y solo si det(A)0\displaystyle\det(A)\neq 0.

  7. 7.

    det(AB)=det(A)det(B)\displaystyle\det(A\cdot B)=\det(A)\cdot\det(B).

  8. 8.

    det(A)=det(At)\displaystyle\det(A)=\det(A^{t}). En particular, las propiedades (1)(5)\displaystyle(1)-(5) anteriores se cumplen tambien si se cambia la palabra “fila” por la palabra “columna”.

  9. 9.

    Se puede desarrollar el determinante por cualquier columna o cualquier fila de A\displaystyle A.

Demostración.

Las propiedades 1), 2) y 4) se demuestran por inducción. Veamos la demostración del resto de propiedades:

  1. 3.

    Veamos que esta propiedad se deduce de 1) y 2). Expresamos A\displaystyle A por filas:

    A=(AiAj)\displaystyle A=\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)

    y construimos la matriz auxiliar

    A=(Ai+AiAj+Ai)\displaystyle A^{\prime}=\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}+A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}+A_{i}\\ \hline\cr\vdots\\ \end{array}\right)

    Entonces, por 1) y 2),

    0=1)det(A)=2)det(AiAi)=1)0+det(AiAj)+det(AjAi)+det(AjAj)=1)0\displaystyle 0=_{1)}\det(A^{\prime})=_{2)}\underbrace{\det\left(\begin{array}% []{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \end{array}\right)}_{=_{1)}0}+\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)+\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \end{array}\right)+\underbrace{\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)}_{=_{1)}0}

    de donde

    det(A)=det(AiAj)=det(AjAi)\displaystyle\det(A)=\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)=-\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \end{array}\right)
  2. 5.

    Si expresamos A\displaystyle A por filas

    A=(AiAj),\displaystyle A=\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right),

    tomamos λ𝕂\displaystyle\lambda\in\mathbb{K} y vamos a calcular el determinante de la matriz que en la fila i\displaystyle i tiene la fila i\displaystyle i de A\displaystyle A mas λ\displaystyle\lambda por la fila j\displaystyle j:

    det(Ai+λAjAj)=1),4)det(AiAj)+λdet(AjAj)=2)det(AiAj)\displaystyle\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}+\lambda A_{j}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)=_{1),4)}\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)+\lambda\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)=_{2)}\det\left(\begin{array}[]{c}\vdots\\ \hline\cr A_{i}\\ \hline\cr\vdots\\ \hline\cr A_{j}\\ \hline\cr\vdots\\ \end{array}\right)
  3. 6.

    Utilizando las propiedades anteriores y que det(In)=1\displaystyle\det(I_{n})=1, vamos a calcular cuánto valen los determinantes de las matrices elementales Pij,Pij(t)(ij)\displaystyle P_{ij},P_{ij}(t)(i\neq j) y Qi(s)\displaystyle Q_{i}(s):

    •  

      Como Pij\displaystyle P_{ij}, ij\displaystyle i\neq j, multiplicada por delante de la matriz identidad, intercambia las filas i\displaystyle i y j\displaystyle j de la matriz identidad,

      1=3)det(PijIn)=det(Pij).\displaystyle-1=_{3)}\det(P_{ij}I_{n})=\det(P_{ij}).
    •  

      Como Pij(λ)\displaystyle P_{ij}(\lambda), multiplicada por delante de la matriz identidad, suma a la fila i\displaystyle i la fila j\displaystyle j de la matriz identidad multiplicada por λ\displaystyle\lambda, y eso no afecta al valor del determinante,

      1=4)det(Pij(λ)In)=det(Pij(λ)).\displaystyle 1=_{4)}\det(P_{ij}(\lambda)I_{n})=\det(P_{ij}(\lambda)).
    •  

      Como Qi(s)\displaystyle Q_{i}(s), multiplicada por delante de la matriz identidad, multiplica la fila i\displaystyle i de la matriz identidad por s\displaystyle s, el determinante de la matriz resultante es sdet(In)\displaystyle s\det(I_{n}), asi que

      s=5)det(Qi(s)In)=det(Qi(s)).\displaystyle s=_{5)}\det(Q_{i}(s)I_{n})=\det(Q_{i}(s)).

    Las propiedades 3), 4) y 5) se pueden escribir usando matrices elementales. Si A\displaystyle A es una matriz cualquiera, ij\displaystyle i\neq j,

    •  

      det(PijA)=3)det(A)=det(Pij)det(A)\displaystyle\det(P_{ij}A)=_{3)}-\det(A)=\det(P_{ij})\det(A)

    •  

      det(PijA)=4)det(A)=det(Pij(λ))det(A)\displaystyle\det(P_{ij}A)=_{4)}\det(A)=\det(P_{ij}(\lambda))\det(A)

    •  

      det(Qi(s)A)=5)sdet(A)=det(Qi(s))det(A)\displaystyle\det(Q_{i}(s)A)=_{5)}s\det(A)=\det(Q_{i}(s))\det(A)

    y llegamos a la concluson de que si E\displaystyle E denota cualquiera de las matrices elementales y A\displaystyle A es una matriz en general

    det(EA)=det(E)det(A)\displaystyle\det(EA)=\det(E)\det(A)

    Vamos a usar esta propiedad para probar que A\displaystyle A es invertible si y solo si el determinante de A\displaystyle A es no nulo.

    )\displaystyle\Rightarrow) Sabemos que A\displaystyle A es invertible si y solo si A=E1E2Ek\displaystyle A=E_{1}\cdot E_{2}\cdots E_{k} para ciertas matrices elementales E1,,Ek\displaystyle E_{1},\ldots,E_{k}. Por tanto,

    det(A)=det(E1E2Ek)=det(E1)det(E2Ek)=det(E1)0det(E2)0det(Ek)00\displaystyle\det(A)=\det(E_{1}\cdot E_{2}\cdots E_{k})=\det(E_{1})\det(E_{2}% \cdots E_{k})=\underbrace{\det(E_{1})}_{\neq 0}\underbrace{\det(E_{2})}_{\neq 0% }\cdots\underbrace{\det(E_{k})}_{\neq 0}\neq 0

    )\displaystyle\Leftarrow) Veamos que si A\displaystyle A no es invertible, entonces el determinante de A\displaystyle A es cero: si A\displaystyle A no es invertible, cuando apliquemos el método de Gauss-Jordan, que consiste en ir multiplicando A\displaystyle A por delante por matrices elementales hasta encontrar una matriz escalonada reducida, nos encontraremos con una matriz escalonada reducida R\displaystyle R que tiene toda una fila de ceros. Por la propiedad 4), si una matriz tiene toda una fila de ceros, su determinante es cero. Así

    0=4)det(R)=det(E1EkA)=det(E1)0det(Ek)0det(A)\displaystyle 0=_{4)}\det(R)=\det(E_{1}\cdots E_{k}A)=\underbrace{\det(E_{1})}% _{\neq 0}\cdots\underbrace{\det(E_{k})}_{\neq 0}\det(A)

    de donde se obtiene que det(A)=0\displaystyle\det(A)=0.

  4. 7.

    Hasta ahora hemos probado que det(EA)=det(E)det(A)\displaystyle\det(EA)=\det(E)\det(A), para E\displaystyle E una matriz elemental. Veamos que en general det(AB)=det(A)det(B)\displaystyle\det(AB)=\det(A)\det(B), para A\displaystyle A y B\displaystyle B dos matrices cualquiera. Distinguimos dos casos:

    1. a)

      Si det(A)=0\displaystyle\det(A)=0 o det(B)=0\displaystyle\det(B)=0, entonces, por 6) A\displaystyle A o B\displaystyle B no son invertibles, con lo que AB\displaystyle AB no es invertible y, de nuevo, por 6)\displaystyle 6), 0=det(AB)=det(A)det(B)\displaystyle 0=\det(AB)=\det(A)\det(B)

    2. b)

      Si det(A)0\displaystyle\det(A)\neq 0 y det(B)0\displaystyle\det(B)\neq 0, entonces, por 6)\displaystyle 6), tanto A\displaystyle A como B\displaystyle B son invertibles (producto cada una de ellas por matrices elementales). Supongamos que A=E1Ek\displaystyle A=E_{1}\cdots E_{k}, B=E1Es\displaystyle B=E^{\prime}_{1}\cdots E^{\prime}_{s}. Así, por la conclusión anterior,

      det(AB)\displaystyle\det(AB) =det(E1EkE1Es)\displaystyle=\det(E_{1}\cdots E_{k}E^{\prime}_{1}\cdots E^{\prime}_{s})
      =det(E1)det(Ek)det(E1)det(Es)\displaystyle=\det(E_{1})\cdots\det(E_{k})\det(E^{\prime}_{1})\cdots\det(E^{% \prime}_{s})
      =det(A)det(B).\displaystyle=\det(A)\det(B).
  5. 8.

    Distinguimos dos casos:

    1. a)

      Si det(A)=0\displaystyle\det(A)=0, entonces, por 6),\displaystyle 6), A\displaystyle A no es invertible, con lo que At\displaystyle A^{t} tampoco es invertible y, por 6)\displaystyle 6), det(At)=0\displaystyle\det(A^{t})=0.

    2. b)

      Si det(A)0\displaystyle\det(A)\neq 0, entonces, por 6)\displaystyle 6), A\displaystyle A es invertible (A\displaystyle A es producto de matrices elementales A=E1Ek\displaystyle A=E_{1}\cdots E_{k}). Para cada matriz elemental, es fácil comprobar que det(E)=det(Et)\displaystyle\det(E)=\det(E^{t}), así que

      det(At)\displaystyle\det(A^{t}) =det((E1Ek)t)=det(EktE1t)\displaystyle=\det((E_{1}\cdots E_{k})^{t})=\det(E^{t}_{k}\cdots E^{t}_{1})
      =det(Ekt)det(E1t)=det(Ek)det(E1)\displaystyle=\det(E^{t}_{k})\cdots\det(E^{t}_{1})=\det(E_{k})\cdots\det(E_{1})
      =det(E1)det(Ek)=det(E1Ek)\displaystyle=\det(E_{1})\cdots\det(E_{k})=\det(E_{1}\cdots E_{k})
      =det(A)\displaystyle=\det(A)
  6. 9.

    Se sigue de 3) y 8).