Articulo de referencia

Determinante

En matemáticas , el determinante es una función escalar de los elementos de una matriz cuadrada . El determinante de una matriz A se suele denotar como det( A ) , det A o | A | ...

En matemáticas , el determinante es una función escalar de los elementos de una matriz cuadrada . El determinante de una matriz A se suele denotar como det( A ) , det A o | A | . Su valor caracteriza algunas propiedades de la matriz y de la aplicación lineal representada, sobre una base dada , por la matriz. En particular, el determinante es distinto de cero si y solo si la matriz es invertible y la aplicación lineal correspondiente es un isomorfismo . Sin embargo, si el determinante es cero, la matriz se denomina singular, lo que significa que no tiene inversa.

El determinante está completamente determinado por las dos propiedades siguientes: el determinante de un producto de matrices es el producto de sus determinantes, y el determinante de una matriz triangular es el producto de sus elementos diagonales.

El determinante de una matriz de 2 × 2 es

|abdod|=adbdo,{\displaystyle {\begin{vmatrix}a&b\\c&d\end{vmatrix}}=ad-bc,}

y el determinante de una matriz de 3 × 3 es

|abdodmiFgramohi|=amii+bFgramo+dodhdomigramobdiaFh.{\displaystyle {\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}}=aei+bfg+cdh-ceg-bdi-afh.}

El determinante de una matriz n × n se puede definir de varias maneras equivalentes, siendo la más común la fórmula de Leibniz , que expresa el determinante como una suma denorte¡{\displaystyle n!}(el factorial de n ) productos con signo de entradas de matrices. Se puede calcular mediante la expansión de Laplace , que expresa el determinante como una combinación lineal de determinantes de submatrices, o mediante la eliminación gaussiana , que permite calcular una forma escalonada por filas con el mismo determinante, igual al producto de las entradas diagonales de la forma escalonada por filas.

Los determinantes también pueden definirse por algunas de sus propiedades. En concreto, el determinante es la única función definida en matrices de n × n que tiene las cuatro propiedades siguientes:

  1. El determinante de la matriz identidad es 1 .
  2. El intercambio de dos filas multiplica el determinante por −1 .
  3. Al multiplicar una fila por un número, el determinante se multiplica por ese número.
  4. Sumar un múltiplo de una fila a otra fila no cambia el determinante.

Las propiedades anteriores relativas a las filas (propiedades 2 a 4) pueden sustituirse por las declaraciones correspondientes relativas a las columnas.

El determinante es invariante bajo similitud de matrices . Esto implica que, dado un endomorfismo lineal de un espacio vectorial de dimensión finita , el determinante de la matriz que lo representa sobre una base no depende de la base elegida. Esto permite definir el determinante de un endomorfismo lineal, que no depende de la elección del sistema de coordenadas .

Los determinantes aparecen en todas las matemáticas. Por ejemplo, una matriz se usa a menudo para representar los coeficientes en un sistema de ecuaciones lineales , y los determinantes se pueden usar para resolver estas ecuaciones ( regla de Cramer ), aunque otros métodos de solución son computacionalmente mucho más eficientes. Los determinantes se usan para definir el polinomio característico de una matriz cuadrada, cuyas raíces son los autovalores . En geometría , el volumen n -dimensional con signo de un paralelepípedo n- dimensional se expresa mediante un determinante, y el determinante de un endomorfismo lineal determina cómo se transforman la orientación y el volumen n -dimensional bajo el endomorfismo. Esto se usa en cálculo con formas diferenciales exteriores y el determinante jacobiano , en particular para cambios de variables en integrales múltiples .

Matrices de dos por dos

El determinante de una matriz de 2 × 2(abdod){\displaystyle {\begin{pmatrix}a&b\\c&d\end{pmatrix}}}se denota ya sea por " det " o por barras verticales alrededor de la matriz, y se define como

det(abdod)=|abdod|=adbdo.{\displaystyle \det {\begin{pmatrix}a&b\\c&d\end{pmatrix}}={\begin{vmatrix}a&b\\c&d\end{vmatrix}}=ad-bc.}

Por ejemplo,

det(3714)=|3714|=(3(4))(71)=19.{\displaystyle \det {\begin{pmatrix}3&7\\1&-4\end{pmatrix}}={\begin{vmatrix}3&7\\1&{-4}\end{vmatrix}}=(3\cdot (-4))-(7\cdot 1)=-19.}

Primeras propiedades

El determinante tiene varias propiedades clave que pueden probarse mediante la evaluación directa de la definición para2×2{\displaystyle 2\times 2}matrices, y que siguen siendo válidas para determinantes de matrices más grandes. Son las siguientes: [ 1 ] primero, el determinante de la matriz identidad(1001){\displaystyle {\begin{pmatrix}1&0\\0&1\end{pmatrix}}}es 1. Segundo, el determinante es cero si dos filas son iguales:

|abab|=abba=0.{\displaystyle {\begin{vmatrix}a&b\\a&b\end{vmatrix}}=ab-ba=0.}

Esto se cumple de manera similar si las dos columnas son iguales. Además,

|ab+bdod+d|=a(d+d)(b+b)do=|abdod|+|abdod|.{\displaystyle {\begin{vmatrix}a&b+b'\\c&d+d'\end{vmatrix}}=a(d+d')-(b+b')c={\begin{vmatrix}a&b\\c&d\end{vmatrix}}+{\begin{vmatrix}a&b'\\c&d'\end{vmatrix}}.}

Finalmente, si alguna columna se multiplica por algún númeror{\displaystyle r}(es decir, todas las entradas de esa columna se multiplican por ese número), el determinante también se multiplica por ese número:

|rabrdod|=radbrdo=r(adbdo)=r|abdod|.{\displaystyle {\begin{vmatrix}r\cdot a&b\\r\cdot c&d\end{vmatrix}}=rad-brc=r(ad-bc)=r\cdot {\begin{vmatrix}a&b\\c&d\end{vmatrix}}.}

significado geométrico

El área del paralelogramo es el valor absoluto del determinante de la matriz formada por los vectores que representan los lados del paralelogramo.

Si las entradas de la matriz son números reales, la matriz A representa la transformación lineal que asigna los vectores base a las columnas de A. Las imágenes de los vectores base forman un paralelogramo que representa la imagen del cuadrado unitario bajo dicha transformación. El paralelogramo definido por las columnas de la matriz anterior es aquel cuyos vértices se encuentran en (0, 0) , ( a , c ) , ( a + b , c + d ) y ( b , d ) , como se muestra en el diagrama adjunto.

El valor absoluto de adbc es el área del paralelogramo y, por lo tanto, representa el factor de escala por el cual las áreas se transforman mediante A.

El valor absoluto del determinante junto con el signo se convierte en el área con signo del paralelogramo. El área con signo es igual que el área usual , excepto que es negativa cuando el ángulo entre el primer y el segundo vector que define el paralelogramo gira en sentido horario (que es opuesto a la dirección que se obtendría para la matriz identidad ).

Para demostrar que adbc es el área con signo, se puede considerar una matriz que contiene dos vectores u ≡ ( a , c ) y v ≡ ( b , d ) que representan los lados del paralelogramo. El área con signo se puede expresar como | u | | v | sin θ para el ángulo θ entre los vectores, que es simplemente la base por la altura, la longitud de un vector por la componente perpendicular del otro. Debido al seno , esto ya es el área con signo, pero se puede expresar de manera más conveniente usando el coseno del ángulo complementario a un vector perpendicular, por ejemplo u = (− c , a ) , de modo que | u | | v | cos θ se convierte en el área con signo en cuestión, que se puede determinar por el patrón del producto escalar para que sea igual a adbc según las siguientes ecuaciones:

Área señalizada=|||v|pecadoθ=|||v|porqueθ=(doa)(bd)=adbdo.{\displaystyle {\text{Área con signo}}=|{\boldsymbol {u}}|\,|{\boldsymbol {v}}|\,\sin \,\theta =\left|{\boldsymbol {u}}^{\perp }\right|\,\left|{\boldsymbol {v}}\right|\,\cos \,\theta '={\begin{pmatrix}-c\\a\end{pmatrix}}\cdot {\begin{pmatrix}b\\d\end{pmatrix}}=ad-bc.}

Así, el determinante proporciona el factor de escala del área y la orientación inducida por la transformación representada por A. Cuando el determinante es igual a uno, la transformación lineal definida por la matriz conserva el área y la orientación.

El volumen de este paralelepípedo es el valor absoluto del determinante de la matriz formada por las columnas construidas a partir de los vectores r1, r2 y r3.

Si una matriz real A de n × n se escribe en términos de sus vectores columnaA=[a1a2anorte]{\displaystyle A=\left[{\begin{array}{c|c|c|c}\mathbf {a} _{1}&\mathbf {a} _{2}&\cdots &\mathbf {a} _{n}\end{array}}\right]}, entonces

A(100)=a1,A(010)=a2,,A(001)=anorte.{\displaystyle A{\begin{pmatrix}1\\0\\\vdots \\0\end{pmatrix}}=\mathbf {a} _{1},\quad A{\begin{pmatrix}0\\1\\\vdots \\0\end{pmatrix}}=\mathbf {a} _{2},\quad \ldots ,\quad A{\begin{pmatrix}0\\0\\\vdots \\1\end{pmatrix}}=\mathbf {a} _{n}.}

Esto significa queA{\displaystyle A}mapea el cubo unitario n al paralelepípedo n -dimensional definido por los vectoresa1,a2,,anorte,{\displaystyle \mathbf {a} _{1},\mathbf {a} _{2},\ldots ,\mathbf {a} _{n},}la regiónPAG={do1a1++donorteanorte0doi1 i}{\displaystyle P=\left\{c_{1}\mathbf {a} _{1}+\cdots +c_{n}\mathbf {a} _{n}\mid 0\leq c_{i}\leq 1\ \forall i\right\}}({\textstyle \forall }representa "para todos" como símbolo lógico .

El determinante da el volumen n -dimensional con signo de este paralelepípedo,det(A)=±vol(PAG),{\displaystyle \det(A)=\pm {\text{vol}}(P),}y por lo tanto describe de manera más general el factor de escala de volumen n -dimensional de la transformación lineal producida por A. [ 2 ] (El signo muestra si la transformación conserva o invierte la orientación ). En particular, si el determinante es cero, entonces este paralelepípedo tiene volumen cero y no es completamente n -dimensional, lo que indica que la dimensión de la imagen de A es menor que n . Esto significa que A produce una transformación lineal que no es ni sobreyectiva ni biyectiva , y por lo tanto no es invertible.

Definición

Sea A una matriz cuadrada con n filas y n columnas, de modo que se puede escribir como

A=[a1,1a1,2a1,nortea2,1a2,2a2,norteanorte,1anorte,2anorte,norte].{\displaystyle A={\begin{bmatrix}a_{1,1}&a_{1,2}&\cdots &a_{1,n}\\a_{2,1}&a_{2,2}&\cdots &a_{2,n}\\\vdots &\vdots &\ddots &\vdots \\a_{n,1}&a_{n,2}&\cdots &a_{n,n}\end{bmatrix}}.}

Las entradasa1,1{\displaystyle a_{1,1}}etc. son, para muchos propósitos, números reales o complejos. Como se explica más adelante, el determinante también se define para matrices cuyas entradas están en un anillo conmutativo .

El determinante de A se denota por det( A ), o se puede denotar directamente en términos de las entradas de la matriz escribiendo barras en lugar de corchetes:

|a1,1a1,2a1,nortea2,1a2,2a2,norteanorte,1anorte,2anorte,norte|.{\displaystyle {\begin{vmatrix}a_{1,1}&a_{1,2}&\cdots &a_{1,n}\\a_{2,1}&a_{2,2}&\cdots &a_{2,n}\\\vdots &\vdots &\ddots &\vdots \\a_{n,1}&a_{n,2}&\cdots &a_{n,n}\end{vmatrix}}.}

Existen diversas formas equivalentes de definir el determinante de una matriz cuadrada A , es decir, una con el mismo número de filas y columnas: el determinante puede definirse mediante la fórmula de Leibniz , una fórmula explícita que involucra sumas de productos de ciertos elementos de la matriz. El determinante también puede caracterizarse como la función única que depende de los elementos de la matriz y que satisface ciertas propiedades. Este enfoque también puede utilizarse para calcular determinantes simplificando las matrices en cuestión.

Fórmula de Leibniz

matrices de 3 × 3

La fórmula de Leibniz para el determinante de una matriz de 3 × 3 es la siguiente:

|abdodmiFgramohi|=amii+bFgramo+dodhdomigramobdiaFh. {\displaystyle {\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}}=aei+bfg+cdh-ceg-bdi-afh.\ }

En esta expresión, cada término tiene un factor de cada fila, todos en columnas diferentes, ordenados de forma ascendente según la fila. Por ejemplo, bdi tiene b de la primera fila, segunda columna; d de la segunda fila, primera columna; e i de la tercera fila, tercera columna. Los signos se determinan por la cantidad de transposiciones de factores necesarias para ordenarlos de forma ascendente según sus columnas (dado que los términos están ordenados de izquierda a derecha de forma ascendente según la fila): positivo para un número par de transposiciones y negativo para un número impar. En el ejemplo de bdi , la transposición simple de bd a db da dbi, cuyos tres factores son de la primera, segunda y tercera columna respectivamente; este es un número impar de transposiciones, por lo que el término aparece con signo negativo.

El gobierno de Sarrus

La regla de Sarrus es una regla mnemotécnica para la forma expandida de este determinante: la suma de los productos de tres filas diagonales de elementos de la matriz, desde el noroeste hasta el sureste, menos la suma de los productos de tres filas diagonales de elementos, desde el suroeste hasta el noreste, cuando se escriben junto a él las copias de las dos primeras columnas de la matriz, como se muestra en la ilustración. Este método para calcular el determinante de una matriz de 3 × 3 no se aplica a dimensiones superiores.

matrices n × n

Generalizando lo anterior a dimensiones superiores, el determinante de unnorte×norte{\displaystyle n\times n}Una matriz es una expresión que involucra permutaciones y sus signaturas . Una permutación del conjunto{1,2,,norte}{\displaystyle \{1,2,\dots ,n\}}es una función biyectivaσ{\displaystyle \sigma }de este conjunto a sí mismo, con valoresσ(1),σ(2),,σ(norte){\displaystyle \sigma (1),\sigma (2),\ldots ,\sigma (n)}agotando todo el conjunto. El conjunto de todas esas permutaciones, llamado grupo simétrico , se denota comúnmenteSnorte{\displaystyle S_{n}}La firmasgn(σ){\displaystyle \operatorname {sgn}(\sigma )}de una permutaciónσ{\displaystyle \sigma }es+1,{\displaystyle +1,}si la permutación se puede obtener con un número par de transposiciones (intercambios de dos entradas); de lo contrario, es1.{\displaystyle -1.}

Dada una matriz

A=[a1,1a1,norteanorte,1anorte,norte],{\displaystyle A={\begin{bmatrix}a_{1,1}\ldots a_{1,n}\\\vdots \qquad \vdots \\a_{n,1}\ldots a_{n,n}\end{bmatrix}},}

La fórmula de Leibniz para su determinante es, utilizando la notación sigma para la suma,

det(A)=|a1,1a1,norteanorte,1anorte,norte|=σSnortesgn(σ)a1,σ(1)anorte,σ(norte).{\displaystyle \det(A)={\begin{vmatrix}a_{1,1}\ldots a_{1,n}\\\vdots \qquad \vdots \\a_{n,1}\ldots a_{n,n}\end{vmatrix}}=\sum _{\sigma \in S_{n}}\operatorname {sgn}(\sigma )a_{1,\sigma (1)}\cdots a_{n,\sigma (n)}.}

Utilizando la notación pi para el producto, esto se puede abreviar en

det(A)=σSnorte(sgn(σ)i=1norteai,σ(i)){\displaystyle \det(A)=\sum _{\sigma \in S_{n}}\left(\operatorname {sgn}(\sigma )\prod _{i=1}^{n}a_{i,\sigma (i)}\right)}.

El símbolo de Levi-Civitaεi1,,inorte{\displaystyle \varepsilon _{i_{1},\ldots ,i_{n}}}se define en las n - tuplas de enteros en{1,,norte}{\displaystyle \{1,\ldots ,n\}}como 0 si dos de los enteros son iguales, y en caso contrario como la signatura de la permutación definida por la n- tupla de enteros. Con el símbolo de Levi-Civita, la fórmula de Leibniz se convierte en

det(A)=i1,i2,,inorteεi1inortea1,i1anorte,inorte,{\displaystyle \det(A)=\sum _{i_{1},i_{2},\ldots ,i_{n}}\varepsilon _{i_{1}\cdots i_{n}}a_{1,i_{1}}\!\cdots a_{n,i_{n}},}

donde la suma se toma sobre todas las n -tuplas de enteros en{1,,norte}.{\displaystyle \{1,\ldots ,n\}.}[ 3 ] [ 4 ]

Propiedades

Caracterización del determinante

El determinante se puede caracterizar por las siguientes tres propiedades clave. Para enunciarlas, es conveniente considerar unnorte×norte{\displaystyle n\times n}matriz A como compuesta de sunorte{\displaystyle n}columnas, así denominadas como

A=(a1,,anorte),{\displaystyle A={\big (}a_{1},\dots ,a_{n}{\big )},}

donde el vector columnaai{\displaystyle a_{i}}(para cada i ) está compuesta por las entradas de la matriz en la i -ésima columna.

  1. det(I)=1{\displaystyle \det \left(I\right)=1}, dóndeI{\displaystyle I}es una matriz identidad .
  2. El determinante es multilineal : si la j -ésima columna de una matrizA{\displaystyle A}se escribe como una combinación linealaj=rv+w{\displaystyle a_{j}=r\cdot v+w}Si A tiene dos vectores columna v y w y un número r , entonces el determinante de A se puede expresar como una combinación lineal similar:
    |A|=|a1,,aj1,rv+w,aj+1,,anorte|=r|a1,,v,anorte|+|a1,,w,,anorte|{\displaystyle {\begin{aligned}|A|&={\big |}a_{1},\dots ,a_{j-1},r\cdot v+w,a_{j+1},\dots ,a_{n}|\\&=r\cdot |a_{1},\dots ,v,\dots a_{n}|+|a_{1},\dots ,w,\dots ,a_{n}|\end{aligned}}}
  3. El determinante es alternante : siempre que dos columnas de una matriz son idénticas, su determinante es 0:
    |a1,,v,,v,,anorte|=0.{\displaystyle |a_{1},\dots ,v,\dots ,v,\dots ,a_{n}|=0.}

Si el determinante se define utilizando la fórmula de Leibniz como se indicó anteriormente, estas tres propiedades se pueden demostrar mediante la inspección directa de dicha fórmula. Algunos autores también abordan el determinante directamente utilizando estas tres propiedades: se puede demostrar que existe exactamente una función que asigna a cualquiernorte×norte{\displaystyle n\times n}matriz A un número que satisface estas tres propiedades. [ 5 ] Esto también muestra que este enfoque más abstracto del determinante produce la misma definición que la que utiliza la fórmula de Leibniz.

Para comprobarlo, basta con expandir el determinante mediante multilinealidad en las columnas en una (enorme) combinación lineal de determinantes de matrices en la que cada columna es un vector base estándar . Estos determinantes son 0 (si las columnas son linealmente dependientes, según la propiedad 3) o ±1 (según las propiedades 1 y 3; el signo negativo aparece cuando las columnas se permutan según una permutación impar ), por lo que la combinación lineal da la expresión anterior en términos del símbolo de Levi-Civita. Si bien su apariencia es menos técnica, esta caracterización no puede reemplazar por completo la fórmula de Leibniz para definir el determinante, ya que sin ella no está clara la existencia de una función adecuada.

Consecuencias inmediatas

Estas reglas tienen varias consecuencias adicionales:

  • El determinante es una función homogénea , es decir,det(doA)=donortedet(A){\displaystyle \det(cA)=c^{n}\det(A)}(para unnorte×norte{\displaystyle n\times n}matrizA{\displaystyle A}).
  • Intercambiar cualquier par de columnas de una matriz multiplica su determinante por  −1. Esto se deduce de que el determinante es multilineal y alternante (propiedades 2 y 3 anteriores):|a1,,aj,ai,,anorte|=|a1,,ai,,aj,,anorte|.{\displaystyle |a_{1},\dots ,a_{j},\dots a_{i},\dots ,a_{n}|=-|a_{1},\dots ,a_{i},\dots ,a_{j},\dots ,a_{n}|.}Esta fórmula se puede aplicar de forma iterativa cuando se intercambian varias columnas. Por ejemplo|a3,a1,a2,a4,anorte|=|a1,a3,a2,a4,,anorte|=|a1,a2,a3,a4,,anorte|.{\displaystyle |a_{3},a_{1},a_{2},a_{4}\dots ,a_{n}|=-|a_{1},a_{3},a_{2},a_{4},\dots ,a_{n}|=|a_{1},a_{2},a_{3},a_{4},\dots ,a_{n}|.}En términos más generales, cualquier permutación de las columnas multiplica el determinante por el signo de la permutación.
  • Si alguna columna puede expresarse como una combinación lineal de las demás columnas (es decir, las columnas de la matriz forman un conjunto linealmente dependiente ), el determinante es 0. Como caso especial, esto incluye: si alguna columna tiene todos sus elementos iguales a cero, entonces el determinante de esa matriz es 0.
  • Sumar un múltiplo escalar de una columna a otra no altera el valor del determinante. Esto se debe a la multilinealidad y a la alternancia: por multilinealidad, el determinante cambia en un múltiplo del determinante de una matriz con dos columnas iguales, cuyo determinante es cero, ya que el determinante es alternante.
  • SiA{\displaystyle A}es una matriz triangular , es deciraij=0{\displaystyle a_{ij}=0}, cuando seai>j{\displaystyle i>j}o, alternativamente, cuandoi<j{\displaystyle i<j}, entonces su determinante es igual al producto de los elementos de la diagonal:det(A)=a11a22anortenorte=i=1norteaii.{\displaystyle \det(A)=a_{11}a_{22}\cdots a_{nn}=\prod _{i=1}^{n}a_{ii}.}De hecho, dicha matriz puede reducirse, añadiendo adecuadamente múltiplos de las columnas con menos entradas no nulas a aquellas con más entradas, a una matriz diagonal (sin cambiar el determinante). Para dicha matriz, el uso de la linealidad en cada columna reduce a la matriz identidad, en cuyo caso la fórmula enunciada se cumple por la primera propiedad que caracteriza a los determinantes. Alternativamente, esta fórmula también puede deducirse de la fórmula de Leibniz, ya que la única permutaciónσ{\displaystyle \sigma }La permutación identidad es aquella que aporta una contribución distinta de cero.

Ejemplo

Estas propiedades características y sus consecuencias enumeradas anteriormente son teóricamente significativas, pero también pueden utilizarse para calcular determinantes de matrices concretas. De hecho, la eliminación gaussiana puede aplicarse para transformar cualquier matriz en una matriz triangular superior, y los pasos de este algoritmo afectan al determinante de forma controlada. El siguiente ejemplo concreto ilustra el cálculo del determinante de la matrizA{\displaystyle A}utilizando ese método:

A=[212214331].{\displaystyle A={\begin{bmatrix}-2&-1&2\\2&1&4\\-3&3&-1\end{bmatrix}}.}

La combinación de estas igualdades da como resultado|A|=|mi|=(183(1))=54.{\displaystyle |A|=-|E|=-(18\cdot 3\cdot (-1))=54.}

Transponer

El determinante de la transpuesta deA{\displaystyle A}es igual al determinante de A :

det(AT)=det(A){\displaystyle \det \left(A^{\textsf {T}}\right)=\det(A)}.

Esto se puede comprobar examinando la fórmula de Leibniz. [ 6 ] Esto implica que en todas las propiedades mencionadas anteriormente, la palabra "columna" puede sustituirse por "fila". Por ejemplo, si consideramos una matriz n × n compuesta por n filas, el determinante es una función n -lineal.

Multiplicatividad y grupos matriciales

El determinante es una aplicación multiplicativa , es decir, para matrices cuadradas.A{\displaystyle A}yB{\displaystyle B}Si son de igual tamaño, el determinante de un producto matricial es igual al producto de sus determinantes:

det(AB)=det(A)det(B){\displaystyle \det(AB)=\det(A)\det(B)}

Este hecho clave puede probarse observando que, para una matriz fijaB{\displaystyle B}, ambos lados de la ecuación son alternantes y multilineales como una función que depende de las columnas deA{\displaystyle A}. Además, ambos toman el valordetB{\displaystyle \det B}cuandoA{\displaystyle A}es la matriz identidad. La caracterización única de los mapas multilineales alternados mencionada anteriormente demuestra, por lo tanto, esta afirmación. [ 7 ]

Una matrizA{\displaystyle A}Una matriz con entradas en un campo es invertible precisamente si su determinante es distinto de cero. Esto se deduce de la multiplicatividad del determinante y de la fórmula para la inversa que involucra la matriz adjunta mencionada a continuación. En este caso, el determinante de la matriz inversa viene dado por

det(A1)=1det(A)=[det(A)]1{\displaystyle \det \left(A^{-1}\right)={\frac {1}{\det(A)}}=[\det(A)]^{-1}}.

En particular, los productos e inversas de matrices con determinante distinto de cero (respectivamente, determinante uno) aún conservan esta propiedad. Por lo tanto, el conjunto de tales matrices (de tamaño fijo)norte{\displaystyle n}sobre un campoK{\displaystyle K}) forma un grupo conocido como el grupo lineal generalGLnorte(K){\displaystyle \operatorname {GL} _{n}(K)}(respectivamente, un subgrupo llamado grupo lineal especialSLnorte(K)GLnorte(K){\displaystyle \operatorname {SL} _{n}(K)\subset \operatorname {GL} _{n}(K)}. De manera más general, la palabra "especial" indica el subgrupo de otro grupo matricial de matrices de determinante uno. Ejemplos incluyen el grupo ortogonal especial (que si n es 2 o 3 consta de todas las matrices de rotación ), y el grupo unitario especial .

Debido a que el determinante respeta la multiplicación y los inversos, de hecho es un homomorfismo de grupo deGLnorte(K){\displaystyle \operatorname {GL} _{n}(K)}en el grupo multiplicativoK×{\displaystyle K^{\times }}de elementos distintos de cero deK{\displaystyle K}Este homomorfismo es sobreyectivo y su núcleo es SLnorte(K){\displaystyle \operatorname {SL} _{n}(K)}(las matrices con determinante uno). Por lo tanto, según el primer teorema de isomorfismo , esto demuestra queSLnorte(K){\displaystyle \operatorname {SL} _{n}(K)}es un subgrupo normal deGLnorte(K){\displaystyle \operatorname {GL} _{n}(K)}y que el grupo cocienteGLnorte(K)/SLnorte(K){\displaystyle \operatorname {GL} _{n}(K)/\operatorname {SL} _{n}(K)}es isomorfo aK×{\displaystyle K^{\times }}.

La fórmula de Cauchy-Binet es una generalización de la fórmula del producto para matrices rectangulares . Esta fórmula también puede reformularse como una fórmula multiplicativa para matrices compuestas cuyos elementos son los determinantes de todas las submatrices cuadráticas de una matriz dada. [ 8 ] [ 9 ]

expansión de Laplace

La expansión de Laplace expresa el determinante de una matriz.A{\displaystyle A}recursivamente en términos de determinantes de matrices más pequeñas, conocidas como sus menores . El menorMETROi,j{\displaystyle M_{i,j}}se define como el determinante de la(norte1)×(norte1){\displaystyle (n-1)\times (n-1)}matriz que resulta deA{\displaystyle A}eliminando eli{\displaystyle i}-ésima fila y laj{\displaystyle j}-ésima columna. La expresión(1)i+jMETROi,j{\displaystyle (-1)^{i+j}M_{i,j}}se conoce como cofactor . Por cadai{\displaystyle i}, uno tiene la igualdad

det(A)=j=1norte(1)i+jai,jMETROi,j,{\displaystyle \det(A)=\sum _{j=1}^{n}(-1)^{i+j}a_{i,j}M_{i,j},}

que se denomina expansión de Laplace a lo largo de la i- ésima fila . Por ejemplo, la expansión de Laplace a lo largo de la primera fila (i=1{\displaystyle i=1}) da como resultado la siguiente fórmula:

|abdodmiFgramohi|=a|miFhi|b|dFgramoi|+do|dmigramoh|{\displaystyle {\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}}=a{\begin{vmatrix}e&f\\h&i\end{vmatrix}}-b{\begin{vmatrix}d&f\\g&i\end{vmatrix}}+c{\begin{vmatrix}d&e\\g&h\end{vmatrix}}}

Desentrañando los determinantes de estos2×2{\displaystyle 2\times 2}Las matrices devuelven la fórmula de Leibniz mencionada anteriormente. De manera similar, la expansión de Laplace a lo largo de laj{\displaystyle j}La columna -ésima es la igualdad

det(A)=i=1norte(1)i+jai,jMETROi,j.{\displaystyle \det(A)=\sum _{i=1}^{n}(-1)^{i+j}a_{i,j}M_{i,j}.}

La expansión de Laplace se puede utilizar de forma iterativa para calcular determinantes, pero este método es ineficiente para matrices grandes. Sin embargo, es útil para calcular los determinantes de matrices altamente simétricas, como la matriz de Vandermonde.|1111incógnita1incógnita2incógnita3incógnitanorteincógnita12incógnita22incógnita32incógnitanorte2incógnita1norte1incógnita2norte1incógnita3norte1incógnitanortenorte1|=1i<jnorte(incógnitajincógnitai).{\displaystyle {\begin{vmatrix}1&1&1&\cdots &1\\x_{1}&x_{2}&x_{3}&\cdots &x_{n}\\x_{1}^{2}&x_{2}^{2}&x_{3}^{2}&\cdots &x_{n}^{2}\\\vdots &\vdots &\vdots &\ddots &\vdots \\x_{1}^{n-1}&x_{2}^{n-1}&x_{3}^{n-1}&\cdots &x_{n}^{n-1}\end{vmatrix}}=\prod _{1\leq i<j\leq n}\left(x_{j}-x_{i}\right).}La expansión de Laplace de n términos a lo largo de una fila o columna se puede generalizar para escribir un determinante de n x n como una suma de(nortek){\displaystyle {\tbinom {n}{k}}}términos , cada uno el producto del determinante de una submatriz k x k y el determinante de la submatriz complementaria ( n−k ) x ( n−k ).

matriz adjugada

La matriz adjugadaadj(A){\displaystyle \operatorname {adj} (A)}es la transpuesta de la matriz de los cofactores, es decir,

(adj(A))i,j=(1)i+jMETROji.{\displaystyle (\operatorname {adj} (A))_{i,j}=(-1)^{i+j}M_{ji}.}

Para cada matriz, se tiene [ 10 ]

(detA)I=AadjA=(adjA)A.{\displaystyle (\det A)I=A\operatorname {adj} A=(\operatorname {adj} A)\,A.}

Por lo tanto, la matriz adjugada puede utilizarse para expresar la inversa de una matriz no singular :

A1=1detAadjA.{\displaystyle A^{-1}={\frac {1}{\det A}}\operatorname {adj} A.}

Matrices de bloques

La fórmula para el determinante de un2×2{\displaystyle 2\times 2}La matriz anterior sigue siendo válida, bajo supuestos adicionales apropiados, para una matriz de bloques , es decir, una matriz compuesta por cuatro submatrices.A,B,do,D{\displaystyle A,B,C,D}de dimensiónmetro×metro{\displaystyle m\times m},metro×norte{\displaystyle m\times n},norte×metro{\displaystyle n\times m}ynorte×norte{\displaystyle n\times n}, respectivamente. La fórmula más sencilla de este tipo, que puede demostrarse utilizando la fórmula de Leibniz o una factorización que involucre el complemento de Schur , es

det(A0doD)=det(A)det(D)=det(AB0D).{\displaystyle \det {\begin{pmatrix}A&0\\C&D\end{pmatrix}}=\det(A)\det(D)=\det {\begin{pmatrix}A&B\\0&D\end{pmatrix}}.}

SiA{\displaystyle A}es invertible , entonces se deduce con los resultados de la sección sobre multiplicatividad que

det(ABdoD)=det(A)det(ABdoD)det(A1A1B0Inorte)=det(A1)=(detA)1=det(A)det(Imetro0doA1DdoA1B)=det(A)det(DdoA1B),{\displaystyle {\begin{aligned}\det {\begin{pmatrix}A&B\\C&D\end{pmatrix}}&=\det(A)\det {\begin{pmatrix}A&B\\C&D\end{pmatrix}}\underbrace {\det {\begin{pmatrix}A^{-1}&-A^{-1}B\\0&I_{n}\end{pmatrix}}} _{=\,\det(A^{-1})\,=\,(\det A)^{-1}}\\&=\det(A)\det {\begin{pmatrix}I_{m}&0\\CA^{-1}&D-CA^{-1}B\end{pmatrix}}\\&=\det(A)\det(D-CA^{-1}B),\end{aligned}}}

lo cual se simplifica adet(A)(DdoA1B){\displaystyle \det(A)(D-CA^{-1}B)}cuandoD{\displaystyle D}es un1×1{\displaystyle 1\times 1}matriz.

Se obtiene un resultado similar cuandoD{\displaystyle D}es invertible, es decir

det(ABdoD)=det(D)det(ABdoD)det(Imetro0D1doD1)=det(D1)=(detD)1=det(D)det(ABD1doBD10Inorte)=det(D)det(ABD1do).{\displaystyle {\begin{aligned}\det {\begin{pmatrix}A&B\\C&D\end{pmatrix}}&=\det(D)\det {\begin{pmatrix}A&B\\C&D\end{pmatrix}}\underbrace {\det {\begin{pmatrix}I_{m}&0\\-D^{-1}C&D^{-1}\end{pmatrix}}} _{=\,\det(D^{-1})\,=\,(\det D)^{-1}}\\&=\det(D)\det {\begin{pmatrix}A-BD^{-1}C&BD^{-1}\\0&I_{n}\end{pmatrix}}\\&=\det(D)\det(A-BD^{-1}C).\end{aligned}}}

Ambos resultados pueden combinarse para derivar el teorema del determinante de Sylvester , que también se enuncia a continuación.

Si los bloques son matrices cuadradas del mismo tamaño, se cumplen otras fórmulas. Por ejemplo, sido{\displaystyle C}yD{\displaystyle D}desplazamiento (es decir,doD=Ddo{\displaystyle CD=DC}), entonces [ 11 ]

det(ABdoD)=det(ADBdo).{\displaystyle \det {\begin{pmatrix}A&B\\C&D\end{pmatrix}}=\det(AD-BC).}

Esta fórmula se ha generalizado a matrices compuestas por más de2×2{\displaystyle 2\times 2}bloques, nuevamente bajo condiciones de conmutatividad apropiadas entre los bloques individuales. [ 12 ]

ParaA=D{\displaystyle A=D}yB=do{\displaystyle B=C}, la siguiente fórmula se cumple (incluso siA{\displaystyle A}yB{\displaystyle B}no se desplacen diariamente).

det(ABBA)=det(A+BBB+AA)=det(A+BB0AB)=det(A+B)det(AB).{\displaystyle \det {\begin{pmatrix}A&B\\B&A\end{pmatrix}}=\det {\begin{pmatrix}A+B&B\\B+A&A\end{pmatrix}}=\det {\begin{pmatrix}A+B&B\\0&A-B\end{pmatrix}}=\det(A+B)\det(A-B).}

Es posible calcular el determinante mediante matrices de bloques de forma rápida utilizando algoritmos rápidos de multiplicación de matrices en el tiempoO(norteω){\displaystyle O({n^{\omega }})}para 2.37ω<3{\displaystyle ~2.37\leq \omega <3}, por elLU{\displaystyle LU}descomposición. [ 13 ]

Teorema determinante de Sylvester

El teorema del determinante de Sylvester establece que para A , una matriz m × n , y B , una matriz n × m (de modo que A y B tienen dimensiones que permiten que se multipliquen en cualquier orden formando una matriz cuadrada):

det(Imetro+AB)=det(Inorte+BA),{\displaystyle \det \left(I_{\mathit {m}}+AB\right)=\det \left(I_{\mathit {n}}+BA\right),}

donde I m e I n son las matrices identidad de m × m y n × n , respectivamente.

De este resultado general se derivan varias consecuencias.

  1. Para el caso de un vector columna c y un vector fila r , cada uno con m componentes, la fórmula permite calcular rápidamente el determinante de una matriz que difiere de la matriz identidad en una matriz de rango 1:
    det(Imetro+dor)=1+rdo.{\displaystyle \det \left(I_{\mathit {m}}+cr\right)=1+rc.}
  2. De manera más general, [ 14 ] para cualquier matriz invertible m × m X ,
    det(incógnita+AB)=det(incógnita)det(Inorte+Bincógnita1A),{\displaystyle \det(X+AB)=\det(X)\det \left(I_{\mathit {n}}+BX^{-1}A\right),}
  3. Para un vector columna y fila como el anterior:
    det(incógnita+dor)=det(incógnita)det(1+rincógnita1do)=det(incógnita)+radj(incógnita)do.{\displaystyle \det(X+cr)=\det(X)\det \left(1+rX^{-1}c\right)=\det(X)+r\,\operatorname {adj} (X)\,c.}
  4. Para matrices cuadradasA{\displaystyle A}yB{\displaystyle B}del mismo tamaño, las matricesAB{\displaystyle AB}yBA{\displaystyle BA}tienen los mismos polinomios característicos (por lo tanto, los mismos valores propios).

Una generalización esdet(Z+AWB)=det(Z)det(W)det(W1+BZ1A){\displaystyle \det \left(Z+AWB\right)=\det \left(Z\right)\det \left(W\right)\det \left(W^{-1}+BZ^{-1}A\right)}(véase el lema del determinante de matrices ), donde Z es una matriz invertible de m × m y W es una matriz invertible de n × n .

Suma

El determinante de la sumaA+B{\displaystyle A+B}La relación entre dos matrices cuadradas del mismo tamaño no se puede expresar en general en términos de los determinantes de A y de B.

Sin embargo, para matrices semidefinidas positivasA{\displaystyle A},B{\displaystyle B}ydo{\displaystyle C}de igual tamaño, det(A+B+do)+det(do)det(A+do)+det(B+do),{\displaystyle \det(A+B+C)+\det(C)\geq \det(A+C)+\det(B+C){\text{,}}} con el corolario [ 15 ] [ 16 ]det(A+B)det(A)+det(B).{\displaystyle \det(A+B)\geq \det(A)+\det(B){\text{.}}}

El teorema de Brunn-Minkowski implica que la raíz n -ésima de un determinante es una función cóncava , cuando se restringe a funciones hermitianas definidas positivas.norte×norte{\displaystyle n\times n}matrices. [ 17 ] Por lo tanto, si A y B son hermíticas definidas positivasnorte×norte{\displaystyle n\times n}matrices, una tiene det(A+B)nortedet(A)norte+det(B)norte,{\displaystyle {\sqrt[{n}]{\det(A+B)}}\geq {\sqrt[{n}]{\det(A)}}+{\sqrt[{n}]{\det(B)}},}ya que la raíz enésima del determinante es una función homogénea .

Identidad de la suma para matrices de 2×2

Para el caso especial de2×2{\displaystyle 2\times 2}En matrices con entradas complejas, el determinante de la suma se puede escribir en términos de determinantes y trazas en la siguiente identidad:

det(A+B)=det(A)+det(B)+tr(A)tr(B)tr(AB).{\displaystyle \det(A+B)=\det(A)+\det(B)+{\text{tr}}(A){\text{tr}}(B)-{\text{tr}}(AB).}

Propiedades del determinante en relación con otras nociones

Autovalores y polinomio característico

El determinante está estrechamente relacionado con otros dos conceptos centrales del álgebra lineal: los valores propios y el polinomio característico de una matriz.A{\displaystyle A}frijolnorte×norte{\displaystyle n\times n}matriz con entradas complejas . Entonces, por el Teorema Fundamental del Álgebra,A{\displaystyle A}debe tener exactamente n valores propiosλ1,λ2,,λnorte{\displaystyle \lambda _{1},\lambda _{2},\ldots ,\lambda _{n}}. (Aquí se entiende que un valor propio con multiplicidad algebraica μ aparece μ veces en esta lista). Entonces, resulta que el determinante de A es igual al producto de estos valores propios,

det(A)=i=1norteλi=λ1λ2λnorte.{\displaystyle \det(A)=\prod _{i=1}^{n}\lambda _{i}=\lambda _{1}\lambda _{2}\cdots \lambda _{n}.}

El producto de todos los autovalores distintos de cero se denomina pseudodeterminante .

De esto se deduce inmediatamente que el determinante de una matrizA{\displaystyle A}es cero si y solo si0{\displaystyle 0}es un valor propio deA{\displaystyle A}. En otras palabras,A{\displaystyle A}es invertible si y solo si0{\displaystyle 0}no es un valor propio deA{\displaystyle A}.

El polinomio característico se define como [ 18 ]

χA(t)=det(tIA).{\displaystyle \chi _{A}(t)=\det(t\cdot I-A).}

Aquí,t{\displaystyle t}es la indeterminada del polinomio yI{\displaystyle I}es la matriz identidad del mismo tamaño queA{\displaystyle A}Mediante este polinomio, se pueden utilizar determinantes para hallar los valores propios de la matriz.A{\displaystyle A}: son precisamente las raíces de este polinomio, es decir, esos números complejos.λ{\displaystyle \lambda }de tal manera que

χA(λ)=0.{\displaystyle \chi _{A}(\lambda )=0.}

Una matriz hermitiana es definida positiva si todos sus valores propios son positivos. El criterio de Sylvester afirma que esto es equivalente a los determinantes de las submatrices.

Ak:=[a1,1a1,2a1,ka2,1a2,2a2,kak,1ak,2ak,k]{\displaystyle A_{k}:={\begin{bmatrix}a_{1,1}&a_{1,2}&\cdots &a_{1,k}\\a_{2,1}&a_{2,2}&\cdots &a_{2,k}\\\vdots &\vdots &\ddots &\vdots \\a_{k,1}&a_{k,2}&\cdots &a_{k,k}\end{bmatrix}}}

ser positivo, para todosk{\displaystyle k}entre1{\displaystyle 1}ynorte{\displaystyle n}. [ 19 ]

Rastro

La traza tr( A ) es por definición la suma de las entradas diagonales de A y también es igual a la suma de los valores propios. Por lo tanto, para matrices complejas A ,

det(exp(A))=exp(tr(A)){\displaystyle \det(\exp(A))=\exp(\operatorname {tr} (A))}

o, para matrices reales A ,

tr(A)=registro(det(exp(A))).{\displaystyle \operatorname {tr} (A)=\log(\det(\exp(A))).}

Aquí exp( A ) denota la exponencial matricial de A , porque cada valor propio λ de A corresponde al valor propio exp( λ ) de exp( A ). En particular, dado cualquier logaritmo de A , es decir, cualquier matriz L que satisfaga

exp(L)=A{\displaystyle \exp(L)=A}

El determinante de A viene dado por

det(A)=exp(tr(L)).{\displaystyle \det(A)=\exp(\operatorname {tr} (L)).}

Por ejemplo, para n = 2 , n = 3 y n = 4 , respectivamente,

det(A)=12((tr(A))2tr(A2)),det(A)=16((tr(A))33tr(A) tr(A2)+2tr(A3)),det(A)=124((tr(A))46tr(A2)(tr(A))2+3(tr(A2))2+8tr(A3) tr(A)6tr(A4)).{\displaystyle {\begin{aligned}\det(A)&={\frac {1}{2}}\left(\left(\operatorname {tr} (A)\right)^{2}-\operatorname {tr} \left(A^{2}\right)\right),\\\det(A)&={\frac {1}{6}}\left(\left(\operatorname {tr} (A)\right)^{3}-3\operatorname {tr} (A)~\operatorname {tr} \left(A^{2}\right)+2\operatorname {tr} \left(A^{3}\right)\right),\\\det(A)&={\frac {1}{24}}\left(\left(\operatorname {tr} (A)\right)^{4}-6\operatorname {tr} \left(A^{2}\right)\left(\operatorname {tr} (A)\right)^{2}+3\left(\operatorname {tr} \left(A^{2}\right)\right)^{2}+8\operatorname {tr} \left(A^{3}\right)~\operatorname {tr} (A)-6\operatorname {tr} \left(A^{4}\right)\right).\end{aligned}}}

cf. teorema de Cayley-Hamilton . Tales expresiones se pueden deducir de argumentos combinatorios, identidades de Newton o el algoritmo de Faddeev-LeVerrier . Es decir, para un n genérico , det A = (−1) n c 0 el término constante con signo del polinomio característico , determinado recursivamente a partir de

donorte=1;   donortemetro=1metrok=1metrodonortemetro+ktr(Ak)  (1metronorte) .{\displaystyle c_{n}=1;~~~c_{n-m}=-{\frac {1}{m}}\sum _{k=1}^{m}c_{n-m+k}\operatorname {tr} \left(A^{k}\right)~~(1\leq m\leq n)~.}

En el caso general, esto también se puede obtener de [ 20 ].

det(A)=k1,k2,,knorte0k1+2k2++norteknorte=nortel=1norte(1)kl+1lklkl¡tr(Al)kl,{\displaystyle \det(A)=\sum _{\begin{array}{c}k_{1},k_{2},\ldots ,k_{n}\geq 0\\k_{1}+2k_{2}+\cdots +nk_{n}=n\end{array}}\prod _{l=1}^{n}{\frac {(-1)^{k_{l}+1}}{l^{k_{l}}k_{l}!}}\operatorname {tr} \left(A^{l}\right)^{k_{l}},}

donde la suma se toma sobre el conjunto de todos los enteros k l ≥ 0 que satisfacen la ecuación

l=1nortelkl=norte.{\displaystyle \sum _{l=1}^{n}lk_{l}=n.}

La fórmula se puede expresar en términos del polinomio exponencial completo de Bell de n argumentos s l = −( l – 1)! tr( A l ) como

det(A)=(1)nortenorte¡Bnorte(s1,s2,,snorte).{\displaystyle \det(A)={\frac {(-1)^{n}}{n!}}B_{n}(s_{1},s_{2},\ldots ,s_{n}).}

Esta fórmula también se puede utilizar para hallar el determinante de una matriz A I J con índices multidimensionales I = ( i 1 , i 2 , ..., i r ) y J = ( j 1 , j 2 , ..., j r ) . El producto y la traza de dichas matrices se definen de forma natural como

(AB)JI=KAKIBJK,tr(A)=IAII.{\displaystyle (AB)_{J}^{I}=\sum _{K}A_{K}^{I}B_{J}^{K},\operatorname {tr} (A)=\sum _{I}A_{I}^{I}.}

Una importante identidad de dimensión arbitraria n se puede obtener a partir del desarrollo en serie de Mercator del logaritmo cuando el desarrollo converge. Si cada valor propio de A es menor que 1 en valor absoluto,

det(I+A)=k=01k¡(j=1(1)jjtr(Aj))k,{\displaystyle \det(I+A)=\sum _{k=0}^{\infty }{\frac {1}{k!}}\left(-\sum _{j=1}^{\infty }{\frac {(-1)^{j}}{j}}\operatorname {tr} \left(A^{j}\right)\right)^{k}\,,}

donde I es la matriz identidad. De forma más general, si

k=01k¡(j=1(1)jsjjtr(Aj))k,{\displaystyle \sum _{k=0}^{\infty }{\frac {1}{k!}}\left(-\sum _{j=1}^{\infty }{\frac {(-1)^{j}s^{j}}{j}}\operatorname {tr} \left(A^{j}\right)\right)^{k}\,,}

se expande como una serie de potencias formal en s entonces todos los coeficientes de s m para m > n son cero y el polinomio restante es det( I + sA ) .

Límites superior e inferior

Para una matriz definida positiva A , el operador traza proporciona los siguientes límites inferiores y superiores ajustados para el determinante logarítmico.

tr(IA1)registrodet(A)tr(AI){\displaystyle \operatorname {tr} \left(I-A^{-1}\right)\leq \log \det(A)\leq \operatorname {tr} (A-I)}

con igualdad si y solo si A = I. Esta relación se puede derivar mediante la fórmula de la divergencia de Kullback-Leibler entre dos distribuciones normales multivariadas .

También,

nortetr(A1)det(A)1norte1nortetr(A)1nortetr(A2).{\displaystyle {\frac {n}{\operatorname {tr} \left(A^{-1}\right)}}\leq \det(A)^{\frac {1}{n}}\leq {\frac {1}{n}}\operatorname {tr} (A)\leq {\sqrt {{\frac {1}{n}}\operatorname {tr} \left(A^{2}\right)}}.}

Estas desigualdades pueden demostrarse expresando las trazas y el determinante en términos de los autovalores. De este modo, representan el hecho bien conocido de que la media armónica es menor que la media geométrica , que es menor que la media aritmética , que a su vez es menor que la raíz cuadrática media .

Derivado

La fórmula de Leibniz muestra que el determinante de matrices cuadradas reales (o análogamente para matrices cuadradas complejas) es una función polinómica deRnorte×norte{\displaystyle \mathbf {R} ^{n\times n}}aR{\displaystyle \mathbf {R} }. En particular, es diferenciable en todas partes . Su derivada se puede expresar utilizando la fórmula de Jacobi : [ 21 ]

ddet(A)dα=tr(adj(A)dAdα).{\displaystyle {\frac {d\det(A)}{d\alpha }}=\operatorname {tr} \left(\operatorname {adj} (A){\frac {dA}{d\alpha }}\right).}

dóndeadj(A){\displaystyle \operatorname {adj} (A)}denota el adjugado deA{\displaystyle A}. En particular, siA{\displaystyle A}es invertible, tenemos

ddet(A)dα=det(A)tr(A1dAdα).{\displaystyle {\frac {d\det(A)}{d\alpha }}=\det(A)\operatorname {tr} \left(A^{-1}{\frac {dA}{d\alpha }}\right).}

Expresado en términos de las entradas deA{\displaystyle A}, estos son

det(A)Aij=adj(A)ji=det(A)(A1)ji.{\displaystyle {\frac {\partial \det(A)}{\partial A_{ij}}}=\operatorname {adj} (A)_{ji}=\det(A)\left(A^{-1}\right)_{ji}.}

Otra formulación equivalente es

det(A+ϵincógnita)det(A)=tr(adj(A)incógnita)ϵ+O(ϵ2)=det(A)tr(A1incógnita)ϵ+O(ϵ2){\displaystyle \det(A+\epsilon X)-\det(A)=\operatorname {tr} (\operatorname {adj} (A)X)\epsilon +O\left(\epsilon ^{2}\right)=\det(A)\operatorname {tr} \left(A^{-1}X\right)\epsilon +O\left(\epsilon ^{2}\right)},

utilizando la notación de la gran O. El caso especial dondeA=I{\displaystyle A=I}, la matriz identidad, produce

det(I+ϵincógnita)=1+tr(incógnita)ϵ+O(ϵ2).{\displaystyle \det(I+\epsilon X)=1+\operatorname {tr} (X)\epsilon +O\left(\epsilon ^{2}\right).}

Esta identidad se utiliza para describir álgebras de Lie asociadas a ciertos grupos de Lie de matrices . Por ejemplo, el grupo lineal especial.SLnorte{\displaystyle \operatorname {SL} _{n}}se define mediante la ecuacióndetA=1{\displaystyle \det A=1}La fórmula anterior muestra que su álgebra de Lie es el álgebra de Lie lineal especial .slnorte{\displaystyle {\mathfrak {sl}}_{n}}que consiste en aquellas matrices cuya traza es cero.

Escribir un3×3{\displaystyle 3\times 3}matriz comoA=[abdo]{\displaystyle A={\begin{bmatrix}a&b&c\end{bmatrix}}}dóndea,b,do{\displaystyle a,b,c}Si son vectores columna de longitud 3, entonces el gradiente sobre uno de los tres vectores se puede escribir como el producto vectorial de los otros dos:

adet(A)=b×dobdet(A)=do×adodet(A)=a×b.{\displaystyle {\begin{aligned}\nabla _{\mathbf {a} }\det(A)&=\mathbf {b} \times \mathbf {c} \\\nabla _{\mathbf {b} }\det(A)&=\mathbf {c} \times \mathbf {a} \\\nabla _{\mathbf {c} }\det(A)&=\mathbf {a} \times \mathbf {b} .\end{aligned}}}

Historia

Históricamente, los determinantes se utilizaron mucho antes que las matrices: un determinante se definió originalmente como una propiedad de un sistema de ecuaciones lineales . El determinante "determina" si el sistema tiene una solución única (lo cual ocurre precisamente si el determinante es distinto de cero). En Europa, Cardano expresó en 1545 las soluciones de sistemas lineales de dos ecuaciones mediante una entidad similar a un determinante. [ 22 ]

Los determinantes propiamente dichos surgieron por separado del trabajo de Seki Takakazu en 1683 en Japón y, paralelamente, del de Leibniz en 1693. [ 23 ] [ 24 ] [ 25 ] [ 26 ] Cramer (1750) enunció, sin demostración, la regla de Cramer. [ 27 ] Tanto Cramer como Bézout (1779) llegaron a los determinantes a partir de la cuestión de las curvas planas que pasan por un conjunto dado de puntos. [ 28 ]

Vandermonde (1771) fue el primero en reconocer los determinantes como funciones independientes. [ 24 ] Laplace (1772) dio el método general para expandir un determinante en términos de sus menores complementarios : Vandermonde ya había dado un caso especial. [ 29 ] Inmediatamente después, Lagrange (1773) trató determinantes de segundo y tercer orden y lo aplicó a cuestiones de la teoría de la eliminación ; demostró muchos casos especiales de identidades generales.

Gauss (1801) realizó el siguiente avance. Al igual que Lagrange, utilizó mucho los determinantes en la teoría de números . Introdujo la palabra "determinante" (Laplace había usado "resultante"), aunque no con el significado actual, sino más bien aplicada al discriminante de una forma cuadrática . [ 30 ] Gauss también llegó a la noción de determinantes recíprocos (inversos) y estuvo muy cerca del teorema de la multiplicación.

El siguiente contribuyente importante es Binet (1811, 1812), quien enunció formalmente el teorema relativo al producto de dos matrices de m columnas y n filas, que para el caso especial de m = n se reduce al teorema de la multiplicación. El mismo día (30 de noviembre de 1812) que Binet presentó su trabajo a la Academia, Cauchy también presentó uno sobre el tema. (Véase la fórmula de Cauchy-Binet ). En este, utilizó la palabra "determinante" en su sentido actual, [ 31 ] [ 32 ] resumió y simplificó lo que se sabía entonces sobre el tema, mejoró la notación y dio el teorema de la multiplicación con una demostración más satisfactoria que la de Binet. [ 24 ] [ 33 ] Con él comienza la teoría en su generalidad.

Jacobi (1841) utilizó el determinante funcional que Sylvester posteriormente denominó jacobiano . [ 34 ] En sus memorias publicadas en el Crelle's Journal de 1841, trata especialmente este tema, así como la clase de funciones alternantes que Sylvester denominó alternantes . Alrededor de la época de las últimas memorias de Jacobi, Sylvester (1839) y Cayley comenzaron su trabajo. Cayley (1841) introdujo la notación moderna para el determinante utilizando barras verticales. [ 35 ] [ 36 ]

El estudio de formas especiales de determinantes ha sido el resultado natural de la finalización de la teoría general. Los determinantes axisimétricos fueron estudiados por Lebesgue , Hesse y Sylvester; los determinantes persimétricos por Sylvester y Hankel ; los circulantes por Catalan , Spottiswoode , Glaisher y Scott; los determinantes sesgados y los pfaffianos , en conexión con la teoría de la transformación ortogonal , por Cayley; los continuos por Sylvester; los wronskianos (llamados así por Muir ) por Christoffel y Frobenius ; los determinantes compuestos por Sylvester, Reiss y Picquet; los jacobianos y hessianos por Sylvester; y los determinantes gauche simétricos por Trudi . De los libros de texto sobre el tema, el de Spottiswoode fue el primero. En Estados Unidos, Hanus (1886), Weld (1893) y Muir/Metzler (1933) publicaron tratados.

Aplicaciones

La regla de Cramer

Los determinantes pueden utilizarse para describir las soluciones de un sistema lineal de ecuaciones , escrito en forma matricial comoAincógnita=b{\displaystyle Ax=b}Esta ecuación tiene una solución única.incógnita{\displaystyle x}si y solo sidet(A){\displaystyle \det(A)}es distinto de cero. En este caso, la solución viene dada por la regla de Cramer :

incógnitai=det(Ai)det(A)i=1,2,3,,norte{\displaystyle x_{i}={\frac {\det(A_{i})}{\det(A)}}\qquad i=1,2,3,\ldots ,n}

dóndeAi{\displaystyle A_{i}}es la matriz formada al reemplazar lai{\displaystyle i}-ésima columna deA{\displaystyle A}por el vector columnab{\displaystyle b}. Esto se sigue inmediatamente mediante la expansión en columna del determinante, es decir

det(Ai)=det[a1banorte]{\displaystyle \det(A_{i})=\det {\begin{bmatrix}a_{1}&\ldots &b&\ldots &a_{n}\end{bmatrix}}}

=j=1norteincógnitajdet[a1ai1ajai+1anorte]=incógnitaidet(A){\displaystyle =\sum _{j=1}^{n}x_{j}\det {\begin{bmatrix}a_{1}&\ldots &a_{i-1}&a_{j}&a_{i+1}&\ldots &a_{n}\end{bmatrix}}=x_{i}\det(A)}

donde los vectoresaj{\displaystyle a_{j}}son las columnas de A. La regla también está implícita en la identidad.

Aadj(A)=adj(A)A=det(A)Inorte.{\displaystyle A\,\operatorname {adj} (A)=\operatorname {adj} (A)\,A=\det(A)\,I_{n}.}

La regla de Cramer se puede implementar enO(norte3){\displaystyle \operatorname {O} (n^{3})}tiempo, que es comparable a métodos más comunes de resolución de sistemas de ecuaciones lineales, como LU , QR o descomposición en valores singulares . [ 37 ]

Independencia lineal

Los determinantes pueden utilizarse para caracterizar vectores linealmente dependientes :detA{\displaystyle \det A}es cero si y solo si los vectores columna de la matrizA{\displaystyle A}son linealmente dependientes. [ 38 ] Por ejemplo, dados dos vectores linealmente independientesv1,v2R3{\displaystyle v_{1},v_{2}\in \mathbf {R} ^{3}}, un tercer vectorv3{\displaystyle v_{3}}se encuentra en el plano generado por los dos primeros vectores exactamente si el determinante de la3×3{\displaystyle 3\times 3}La matriz formada por los tres vectores es cero. La misma idea también se utiliza en la teoría de ecuaciones diferenciales : dadas las funcionesF1(incógnita),,Fnorte(incógnita){\displaystyle f_{1}(x),\dots ,f_{n}(x)}(se supone que esnorte1{\displaystyle n-1}veces diferenciable ), el Wronskiano se define como

W(F1,,Fnorte)(incógnita)=|F1(incógnita)F2(incógnita)Fnorte(incógnita)F1(incógnita)F2(incógnita)Fnorte(incógnita)F1(norte1)(incógnita)F2(norte1)(incógnita)Fnorte(norte1)(incógnita)|.{\displaystyle W(f_{1},\ldots ,f_{n})(x)={\begin{vmatrix}f_{1}(x)&f_{2}(x)&\cdots &f_{n}(x)\\f_{1}'(x)&f_{2}'(x)&\cdots &f_{n}'(x)\\\vdots &\vdots &\ddots &\vdots \\f_{1}^{(n-1)}(x)&f_{2}^{(n-1)}(x)&\cdots &f_{n}^{(n-1)}(x)\end{vmatrix}}.}

Es distinto de cero (para algunosincógnita{\displaystyle x}) en un intervalo especificado si y solo si las funciones dadas y todas sus derivadas hasta el ordennorte1{\displaystyle n-1}son linealmente independientes. Si se puede demostrar que el Wronskiano es cero en todo un intervalo, entonces, en el caso de funciones analíticas , esto implica que las funciones dadas son linealmente dependientes. Véase el Wronskiano y la independencia lineal . Otro uso del determinante es el resultante , que proporciona un criterio cuando dos polinomios tienen una raíz común . [ 39 ]

Producto cruzado

El cálculo de un producto vectorial es equivalente a encontrar el determinante formal de una matriz específica de 3x3. Específicamente, para vectoresa=a1i+a2j+a3k{\displaystyle \mathbf {a} =a_{1}\mathbf {i} +a_{2}\mathbf {j} +a_{3}\mathbf {k} }, yb=b1i+b2j+b3k{\displaystyle \mathbf {b} =b_{1}\mathbf {i} +b_{2}\mathbf {j} +b_{3}\mathbf {k} }:

a×b=det(ijka1a2a3b1b2b3){\displaystyle \mathbf {a\times b} =\det {\begin{pmatrix}\mathbf {i} &\mathbf {j} &\mathbf {k} \\a_{1}&a_{2}&a_{3}\\b_{1}&b_{2}&b_{3}\\\end{pmatrix}}}

Orientación de una base

El determinante puede entenderse como la asignación de un número a cada secuencia de n vectores en R n , utilizando la matriz cuadrada cuyas columnas son los vectores dados. El determinante será distinto de cero si y solo si la secuencia de vectores es una base para R n . En ese caso, el signo del determinante determina si la orientación de la base es consistente o opuesta a la orientación de la base estándar . En el caso de una base ortogonal, la magnitud del determinante es igual al producto de las longitudes de los vectores de la base. Por ejemplo, una matriz ortogonal con entradas en R n representa una base ortonormal en el espacio euclidiano , y por lo tanto tiene un determinante de ±1 (ya que todos los vectores tienen longitud 1). El determinante es +1 si y solo si la base tiene la misma orientación. Es −1 si y solo si la base tiene la orientación opuesta.

En términos más generales, si el determinante de A es positivo, A representa una transformación lineal que preserva la orientación (si A es una matriz ortogonal de 2 × 2 o 3 × 3 , se trata de una rotación ), mientras que si es negativo, A cambia la orientación de la base.

Volumen y determinante jacobiano

Como se indicó anteriormente, el valor absoluto del determinante de vectores reales es igual al volumen del paralelepípedo formado por dichos vectores. En consecuencia, siF:RnorteRnorte{\displaystyle f:\mathbf {R} ^{n}\to \mathbf {R} ^{n}}es el mapeo lineal dado por la multiplicación con una matrizA{\displaystyle A}, ySRnorte{\displaystyle S\subset \mathbf {R} ^{n}}es cualquier subconjunto medible , entonces el volumen deF(S){\displaystyle f(S)}es dado por|det(A)|{\displaystyle |\det(A)|}veces el volumen deS{\displaystyle S}. [ 40 ] De manera más general, si el mapa linealF:RnorteRmetro{\displaystyle f:\mathbf {R} ^{n}\to \mathbf {R} ^{m}}está representado por elmetro×norte{\displaystyle m\times n}matrizA{\displaystyle A}, entonces la relación entre elnorte{\displaystyle n}- volúmenes dimensionales deF(S){\displaystyle f(S)}yS{\displaystyle S}está dado por:

volumen(F(S))volumen(S)=det(ATA).{\displaystyle {\frac {\operatorname {volume} (f(S))}{\operatorname {volume} (S)}}={\sqrt {\det \left(A^{\textsf {T}}A\right)}}.}

Cuandometro<norte{\displaystyle m<n}Esto es cero.

Calculando el volumen del tetraedro delimitado por cuatro puntos, se pueden utilizar para identificar líneas oblicuas . El volumen de cualquier tetraedro, dados sus vérticesa,b,do,d{\displaystyle a,b,c,d},16|det(ab,bdo,dod)|{\displaystyle {\frac {1}{6}}\cdot |\det(a-b,b-c,c-d)|}o cualquier otra combinación de pares de vértices que formen un árbol de expansión sobre los vértices.

Un mapa no linealF:R2R2{\displaystyle f\colon \mathbf {R} ^{2}\to \mathbf {R} ^{2}}Se envía un pequeño cuadrado (izquierda, en rojo) a un paralelogramo deformado (derecha, en rojo). El jacobiano en un punto proporciona la mejor aproximación lineal del paralelogramo deformado cerca de ese punto (derecha, en blanco translúcido), y el determinante jacobiano da la razón entre el área del paralelogramo aproximado y la del cuadrado original.

Para una función diferenciable general , gran parte de lo anterior se extiende al considerar la matriz jacobiana de f . Para

F:RnorteRnorte,{\displaystyle f:\mathbf {R} ^{n}\rightarrow \mathbf {R} ^{n},}

La matriz jacobiana es la matriz n × n cuyas entradas vienen dadas por las derivadas parciales.

D(F)=(Fiincógnitaj)1i,jnorte.{\displaystyle D(f)=\left({\frac {\partial f_{i}}{\partial x_{j}}}\right)_{1\leq i,j\leq n}.}

Su determinante, el determinante jacobiano , aparece en la versión de mayor dimensión de la integración por sustitución : para funciones adecuadas f y un subconjunto abierto U de R n (el dominio de f ), la integral sobre f ( U ) de alguna otra función φ  : R nR m viene dada por

F(U)ϕ(v)dv=Uϕ(F())|det(DF)()|d.{\displaystyle \int _{f(U)}\phi (\mathbf {v} )\,d\mathbf {v} =\int _{U}\phi (f(\mathbf {u} ))\left|\det(\operatorname {D} f)(\mathbf {u} )\right|\,d\mathbf {u} .}

El jacobiano también aparece en el teorema de la función inversa .

Cuando se aplica al campo de la cartografía , el determinante puede utilizarse para medir la tasa de expansión de un mapa cerca de los polos. [ 41 ]

Áreas y colinealidad

El determinante proporciona una forma conveniente de calcular el área de un triángulo en el plano xy. El área de un triángulo cuyos vértices son(incógnita1,y1){\displaystyle (x_{1},y_{1})},(incógnita2,y2){\displaystyle (x_{2},y_{2})}y(incógnita3,y3){\displaystyle (x_{3},y_{3})}está dado por: [ 42 ]

±12det(incógnita1y11incógnita2y21incógnita3y31){\displaystyle \pm {\frac {1}{2}}\det {\begin{pmatrix}x_{1}&y_{1}&1\\x_{2}&y_{2}&1\\x_{3}&y_{3}&1\end{pmatrix}}}

donde está el letrero±{\displaystyle \pm }Se elige para obtener un área positiva. Una consecuencia inmediata de este resultado es que si el determinante es cero, no hay triángulo, lo que implica que los puntos son colineales . Por lo tanto, el determinante proporciona una herramienta útil para comprobar la colinealidad.

Dado que dos puntos identifican de forma única una línea en el plano xy, una de las filas del determinante puede ser reemplazada por cualquier punto arbitrario.(incógnita,y){\displaystyle (x,y)}expresar la ecuación general de una recta que pasa por puntos(incógnita1,y1){\displaystyle (x_{1},y_{1})}y(incógnita2,y2){\displaystyle (x_{2},y_{2})}como:

det(incógnitay1incógnita1y11incógnita2y21)=0.{\displaystyle \det {\begin{pmatrix}x&y&1\\x_{1}&y_{1}&1\\x_{2}&y_{2}&1\end{pmatrix}}=0.}

Aspectos algebraicos abstractos

Determinante de un endomorfismo

Las identidades anteriores relativas al determinante de productos e inversas de matrices implican que matrices semejantes tienen el mismo determinante: dos matrices A y B son semejantes si existe una matriz invertible X tal que A = X −1 BX . De hecho, al aplicar repetidamente las identidades anteriores se obtiene

det(A)=det(incógnita)1det(B)det(incógnita)=det(B)det(incógnita)1det(incógnita)=det(B).{\displaystyle \det(A)=\det(X)^{-1}\det(B)\det(X)=\det(B)\det(X)^{-1}\det(X)=\det(B).}

Por lo tanto, el determinante también se denomina invariante de similitud . El determinante de una transformación lineal

T:VV{\displaystyle T:V\to V}

Para algún espacio vectorial de dimensión finita V, se define como el determinante de la matriz que lo describe, con respecto a una elección arbitraria de base en V. Debido a la invariancia de similitud, este determinante es independiente de la elección de la base para V y, por lo tanto , solo depende del endomorfismo T.

Matrices cuadradas sobre anillos conmutativos

La definición anterior del determinante utilizando la regla de Leibniz funciona de manera más general cuando las entradas de la matriz son elementos de un anillo conmutativo.R{\displaystyle R}, tales como los números enterosZ{\displaystyle \mathbf {Z} }, a diferencia del campo de los números reales o complejos. Además, la caracterización del determinante como el único mapa multilineal alternante que satisfacedet(I)=1{\displaystyle \det(I)=1}sigue siendo válido, al igual que todas las propiedades que resultan de esa caracterización. [ 43 ]

Una matrizAEsteranorte×norte(R){\displaystyle A\in \operatorname {Mat} _{n\times n}(R)}es invertible (en el sentido de que existe una matriz inversa cuyas entradas están enR{\displaystyle R}) si y solo si su determinante es un elemento invertible enR{\displaystyle R}. [ 44 ] ParaR=Z{\displaystyle R=\mathbf {Z} }, esto significa que el determinante es +1 o −1. Dicha matriz se denomina unimodular .

Al ser el determinante multiplicativo, define un homomorfismo de grupo.

GLnorte(R)R×,{\displaystyle \operatorname {GL} _{n}(R)\rightarrow R^{\times },}

entre el grupo lineal general (el grupo de invertibles)norte×norte{\displaystyle n\times n}-matrices con entradas enR{\displaystyle R}) y el grupo multiplicativo de unidades enR{\displaystyle R}. Dado que respeta la multiplicación en ambos grupos, esta aplicación es un homomorfismo de grupos .

El determinante es una transformación natural.

Dado un homomorfismo de anillosF:RS{\displaystyle f:R\to S}Hay un mapaGLnorte(F):GLnorte(R)GLnorte(S){\displaystyle \operatorname {GL} _{n}(f):\operatorname {GL} _{n}(R)\to \operatorname {GL} _{n}(S)}dado al reemplazar todas las entradas enR{\displaystyle R}por sus imágenes bajoF{\displaystyle f}. El determinante respeta estos mapas, es decir, la identidad

F(det((ai,j)))=det((F(ai,j))){\displaystyle f(\det((a_{i,j})))=\det((f(a_{i,j})))}

se cumple. En otras palabras, el diagrama conmutativo mostrado conmuta.

Por ejemplo, el determinante del conjugado complejo de una matriz compleja (que también es el determinante de su transpuesta conjugada) es el conjugado complejo de su determinante, y para matrices enteras: la reducción módulometro{\displaystyle m}del determinante de dicha matriz es igual al determinante de la matriz reducida módulometro{\displaystyle m}(este último determinante se calcula utilizando aritmética modular ). En el lenguaje de la teoría de categorías , el determinante es una transformación natural entre los dos functores.GLnorte{\displaystyle \operatorname {GL} _{n}}y()×{\displaystyle (-)^{\times }}. [ 45 ] Añadiendo otra capa de abstracción, esto se captura diciendo que el determinante es un morfismo de grupos algebraicos , del grupo lineal general al grupo multiplicativo ,

det:GLnorteGRAMOmetro.{\displaystyle \det :\operatorname {GL} _{n}\to \mathbb {G} _{m}.}

Álgebra exterior

El determinante de una transformación linealT:VV{\displaystyle T:V\to V}de unnorte{\displaystyle n}espacio vectorial de dimensiónV{\displaystyle V}o, más generalmente, un módulo libre de rango (finito)norte{\displaystyle n}sobre un anillo conmutativoR{\displaystyle R}puede formularse de manera independiente de las coordenadas considerando lanorte{\displaystyle n}-otro poder exteriornorteV{\displaystyle \bigwedge ^{n}V}deV{\displaystyle V}. [ 46 ] El mapaT{\displaystyle T}induce un mapa lineal

norteT:norteVnorteVv1v2vnorteTv1Tv2Tvnorte.{\displaystyle {\begin{aligned}\bigwedge ^{n}T:\bigwedge ^{n}V&\rightarrow \bigwedge ^{n}V\\v_{1}\wedge v_{2}\wedge \dots \wedge v_{n}&\mapsto Tv_{1}\wedge Tv_{2}\wedge \dots \wedge Tv_{n}.\end{aligned}}}

ComonorteV{\displaystyle \bigwedge ^{n}V}es unidimensional, el mapanorteT{\displaystyle \bigwedge ^{n}T}se obtiene multiplicando por algún escalar, es decir, un elemento enR{\displaystyle R}. Algunos autores como ( Bourbaki 1998 ) utilizan este hecho para definir el determinante como el elemento enR{\displaystyle R}que satisfacen la siguiente identidad (para todoviV{\displaystyle v_{i}\in V}):

(norteT)(v1vnorte)=det(T)v1vnorte.{\displaystyle \left(\bigwedge ^{n}T\right)\left(v_{1}\wedge \dots \wedge v_{n}\right)=\det(T)\cdot v_{1}\wedge \dots \wedge v_{n}.}

Esta definición concuerda con la definición más concreta dependiente de coordenadas. Esto se puede demostrar utilizando la unicidad de una forma alternada multilineal ennorte{\displaystyle n}-tuplas de vectores enRnorte{\displaystyle R^{n}}Por esta razón, la potencia exterior no nula más altanorteV{\displaystyle \bigwedge ^{n}V}(a diferencia del determinante asociado a un endomorfismo) a veces también se denomina determinante deV{\displaystyle V}y de forma similar para objetos más complejos como haces vectoriales o complejos de cadenas de espacios vectoriales. Los menores de una matriz también se pueden expresar en este contexto, considerando formas alternas inferiores.kV{\displaystyle \bigwedge ^{k}V}conk<norte{\displaystyle k<n}. [ 47 ]

Integral de Berezin

La definición convencional del determinante, como una suma sobre permutaciones sobre un producto de elementos de matriz, puede escribirse utilizando la notación algo sorprendente de la integral de Berezin . En esta notación, el determinante puede escribirse como

exp[θTAη]dθdη=detA{\displaystyle \int \exp \left[-\theta ^{T}A\eta \right]\,d\theta \,d\eta =\det A}

Esto se aplica a cualquiernorte×norte{\displaystyle n\times n}matriz dimensionalA.{\displaystyle A.}Los símbolosθ,η{\displaystyle \theta ,\eta }son dosnorte{\displaystyle n}Vectores -dimensionales de números de Grassmann anticonmutativos (también conocidos como " supernúmeros "), tomados del álgebra de Grassmann .exp{\displaystyle \exp }Aquí está la función exponencial . El símbolo de integral debe entenderse como la integral de Berezin. A pesar del uso del símbolo de integral, esta expresión es, de hecho, una suma completamente finita.

Esta expresión de aspecto inusual puede entenderse como un truco de notación que reescribe la expresión convencional para el determinante.

detA=σSnortesgn(σ)a1,σ(1)anorte,σ(norte).{\displaystyle \det A=\sum _{\sigma \in S_{n}}\operatorname {sgn}(\sigma )a_{1,\sigma (1)}\cdots a_{n,\sigma (n)}.}

mediante el uso de alguna notación novedosa. La propiedad anticonmutativa de los números de Grassmann captura el signo (signatura) de la permutación, mientras que la integral combinada con laexp{\displaystyle \exp }asegura que se exploren todas las permutaciones. Es decir, la serie de Taylor paraexp{\displaystyle \exp }finaliza después exactamentenorte{\displaystyle n}términos, porque el cuadrado de un número de Grassmann es cero, y hay exactamentenorte{\displaystyle n}distintas variables de Grassmann. Mientras tanto, la integral se define para que se anule si el número de Grassmann correspondiente no aparece en el integrando. Por lo tanto, la integral selecciona solo aquellos términos en elexp{\displaystyle \exp }series que tienen exactamentenorte{\displaystyle n}variables distintas; todos los términos de orden inferior se desvanecen. Por lo tanto, la combinación un tanto mágica del signo integral, el uso de variables anticonmutativas y la serie de Taylor paraexp{\displaystyle \exp }Simplemente codifica una suma finita, idéntica a la sumatoria convencional.

Esta forma es popular en física, donde se usa frecuentemente como sustituto del determinante jacobiano. Su atractivo reside en que, en notación, la integral adopta la forma de una integral de trayectoria , como en la formulación de la integral de trayectoria para la mecánica hamiltoniana cuantizada . Un ejemplo se encuentra en la teoría de los fantasmas de Fadeev-Popov ; si bien esta teoría puede parecer bastante abstrusa, conviene recordar que el uso de los campos fantasma no es más que un truco notacional para expresar un determinante jacobiano.

El PfaffianoPAGFA{\displaystyle \mathrm {Pf} \,A}de una matriz antisimétricaA{\displaystyle A}es la raíz cuadrada del determinante: es decir,(PAGFA)2=detA.{\displaystyle \left(\mathrm {Pf} \,A\right)^{2}=\det A.}La forma integral de Berezin para el pfaffiano es aún más sugerente; es

exp[12θTAθ]dθ=PAGFA{\displaystyle \int \exp \left[-{\tfrac {1}{2}}\theta ^{T}A\theta \right]\,d\theta =\mathrm {Pf} \,A}

El integrando tiene exactamente la misma estructura formal que una distribución gaussiana normal , aunque con números de Grassmann en lugar de números reales. Esta semejanza formal explica la aparición ocasional de supernúmeros en la teoría de la dinámica estocástica y las ecuaciones diferenciales estocásticas .

Los determinantes, tal como se trataron anteriormente, admiten varias variantes: el permanente de una matriz se define como el determinante, excepto que los factoressgn(σ){\displaystyle \operatorname {sgn}(\sigma )}Se omiten los elementos que aparecen en la regla de Leibniz. El inmanente generaliza ambos al introducir un carácter del grupo simétrico.Snorte{\displaystyle S_{n}}en el gobierno de Leibniz.

Determinantes para álgebras de dimensión finita

Para cualquier álgebra asociativaA{\displaystyle A}que es de dimensión finita como un espacio vectorial sobre un campoF{\displaystyle F}, hay un mapa determinante [ 48 ]

det:AF.{\displaystyle \det :A\to F.}

Esta definición procede estableciendo el polinomio característico independientemente del determinante, y definiendo el determinante como el término de orden más bajo de este polinomio. Esta definición general recupera el determinante para el álgebra matricial.A=Esteranorte×norte(F){\displaystyle A=\operatorname {Mat} _{n\times n}(F)}, pero también incluye varios casos adicionales, incluido el determinante de un cuaternión ,

det(a+ib+jdo+kd)=a2+b2+do2+d2{\displaystyle \det(a+ib+jc+kd)=a^{2}+b^{2}+c^{2}+d^{2}},

la normanorteL/F:LF{\displaystyle N_{L/F}:L\to F}de una extensión de campo , así como el pfaffiano de una matriz antisimétrica y la norma reducida de un álgebra simple central .

Matrices infinitas

Para matrices con un número infinito de filas y columnas, las definiciones anteriores del determinante no son directamente aplicables. Por ejemplo, en la fórmula de Leibniz, habría que calcular una suma infinita (cuyos términos son productos infinitos). El análisis funcional proporciona diferentes extensiones del determinante para estas situaciones de dimensión infinita, que, sin embargo, solo funcionan para tipos específicos de operadores.

El determinante de Fredholm define el determinante para operadores conocidos como operadores de clase traza mediante una generalización apropiada de la fórmula.

det(I+A)=exp(tr(registro(I+A))).{\displaystyle \det(I+A)=\exp(\operatorname {tr} (\log(I+A))).}

Otra noción de determinante de dimensión infinita es el determinante funcional .

Operadores en álgebras de von Neumann

Para operadores en un factor finito , se puede definir un determinante real positivo llamado determinante de Fuglede-Kadison utilizando la traza canónica. De hecho, para cada estado de traza en un álgebra de von Neumann, existe una noción de determinante de Fuglede-Kadison.

Para matrices sobre anillos no conmutativos, la multilinealidad y las propiedades alternantes son incompatibles para n ≥ 2 , [ 49 ] por lo que no hay una buena definición del determinante en este contexto.

Para matrices cuadradas con entradas en un anillo no conmutativo, existen diversas dificultades para definir determinantes de forma análoga a como se hace para anillos conmutativos. Se puede dar un significado a la fórmula de Leibniz siempre que se especifique el orden del producto, y de manera similar para otras definiciones del determinante, pero la no conmutatividad conlleva la pérdida de muchas propiedades fundamentales del determinante, como la propiedad multiplicativa o que el determinante permanece inalterado bajo la transposición de la matriz. Sobre anillos no conmutativos, no existe una noción razonable de forma multilineal (la existencia de una forma bilineal no nula con un elemento regular de R como valor en algún par de argumentos implica que R es conmutativo). Sin embargo, se han formulado diversas nociones de determinante no conmutativo que preservan algunas de las propiedades de los determinantes, en particular los cuasideterminantes y el determinante de Dieudonné . Para algunas clases de matrices con elementos no conmutativos, se puede definir el determinante y demostrar teoremas de álgebra lineal que son muy similares a sus análogos conmutativos. Ejemplos incluyen el q -determinante en grupos cuánticos, el determinante de Capelli en matrices de Capelli y el bereziniano en supermatrices (es decir, matrices cuyas entradas son elementos deZ2{\displaystyle \mathbb {Z} _{2}}- anillos graduados ). [ 50 ] Las matrices de Manin forman la clase más cercana a las matrices con elementos conmutativos.

Cálculo

Los determinantes se utilizan principalmente como herramienta teórica. Rara vez se calculan explícitamente en álgebra lineal numérica , donde para aplicaciones como la comprobación de invertibilidad y el cálculo de autovalores, el determinante ha sido ampliamente sustituido por otras técnicas. [ 51 ] Sin embargo, la geometría computacional sí utiliza con frecuencia cálculos relacionados con determinantes. [ 52 ]

Aunque el determinante se puede calcular directamente utilizando la regla de Leibniz, este enfoque es extremadamente ineficiente para matrices grandes, ya que esa fórmula requiere calcularnorte¡{\displaystyle n!}(norte{\displaystyle n}factorial ) productos para unnorte×norte{\displaystyle n\times n}matriz. Por lo tanto, el número de operaciones requeridas crece muy rápidamente: es de ordennorte¡{\displaystyle n!}El desarrollo en serie de Laplace es igualmente ineficiente. Por lo tanto, se han desarrollado técnicas más complejas para calcular determinantes.

eliminación gaussiana

La eliminación gaussiana consiste en multiplicar por la izquierda una matriz por matrices elementales para obtener una matriz en forma escalonada por filas . El cálculo se puede restringir a matrices elementales de determinante 1. En este caso, el determinante de la matriz resultante en forma escalonada por filas es igual al determinante de la matriz inicial. Dado que una matriz en forma escalonada por filas es una matriz triangular , su determinante es el producto de los elementos de su diagonal.

Por lo tanto, el determinante se puede calcular casi sin coste alguno a partir del resultado de una eliminación gaussiana.

Métodos de descomposición

Algunos métodos de cálculodet(A){\displaystyle \det(A)}escribiendo la matriz como un producto de matrices cuyos determinantes se pueden calcular más fácilmente. Estas técnicas se denominan métodos de descomposición. Algunos ejemplos son la descomposición LU , la descomposición QR o la descomposición de Cholesky (para matrices definidas positivas ). Estos métodos son de ordenO(norte3){\displaystyle \operatorname {O} (n^{3})}, lo que supone una mejora significativa con respecto aO(norte¡){\displaystyle \operatorname {O} (n!)}. [ 53 ]

Por ejemplo, la descomposición LU expresaA{\displaystyle A}como producto

A=PAGLU.{\displaystyle A=PLU.}

de una matriz de permutaciónPAG{\displaystyle P}(que tiene exactamente un solo1{\displaystyle 1}en cada columna, y en caso contrario ceros), una matriz triangular inferiorL{\displaystyle L}y una matriz triangular superiorU{\displaystyle U}Los determinantes de las dos matrices triangularesL{\displaystyle L}yU{\displaystyle U}se pueden calcular rápidamente, ya que son los productos de las entradas diagonales respectivas. El determinante dePAG{\displaystyle P}es solo la señalε{\displaystyle \varepsilon }de la permutación correspondiente (que es+1{\displaystyle +1}para un número par de permutaciones y es1{\displaystyle -1}para un número impar de permutaciones). Una vez que se conoce dicha descomposición LU paraA{\displaystyle A}, su determinante se calcula fácilmente como

det(A)=εdet(L)det(U).{\displaystyle \det(A)=\varepsilon \det(L)\cdot \det(U).}

Métodos adicionales

El pedidoO(norte3){\displaystyle \operatorname {O} (n^{3})}El resultado obtenido mediante métodos de descomposición se ha mejorado con diferentes métodos. Si dos matrices de ordennorte{\displaystyle n}puede multiplicarse en el tiempoMETRO(norte){\displaystyle M(n)}, dóndeMETRO(norte)nortea{\displaystyle M(n)\geq n^{a}}para algunosa>2{\displaystyle a>2}, entonces hay un algoritmo que calcula el determinante en tiempoO(METRO(norte)){\displaystyle O(M(n))}. [ 54 ] Esto significa, por ejemplo, que unO(norte2.376){\displaystyle \operatorname {O} (n^{2.376})}Existe un algoritmo para calcular el determinante basado en el algoritmo de Coppersmith-Winograd . Este exponente se redujo aún más, a partir de 2016, a 2,373. [ 55 ]

Además de la complejidad del algoritmo, se pueden utilizar otros criterios para comparar algoritmos. Especialmente para aplicaciones relacionadas con matrices sobre anillos, existen algoritmos que calculan el determinante sin divisiones. (Por el contrario, la eliminación de Gauss requiere divisiones). Uno de estos algoritmos, con complejidadO(norte4){\displaystyle \operatorname {O} (n^{4})}Se basa en la siguiente idea: se reemplazan las permutaciones (como en la regla de Leibniz) por los llamados caminos ordenados cerrados , en los que se pueden repetir varios elementos. La suma resultante tiene más términos que en la regla de Leibniz, pero en el proceso se pueden reutilizar varios de estos productos, lo que la hace más eficiente que el cálculo ingenuo con la regla de Leibniz. [ 56 ] Los algoritmos también se pueden evaluar según su complejidad de bits , es decir, cuántos bits de precisión se necesitan para almacenar los valores intermedios que ocurren en el cálculo. Por ejemplo, el método de eliminación gaussiana (o descomposición LU) es de ordenO(norte3){\displaystyle \operatorname {O} (n^{3})}, pero la longitud en bits de los valores intermedios puede volverse exponencialmente larga. [ 57 ] En comparación, el algoritmo de Bareiss , que es un método de división exacta (por lo que sí usa la división, pero solo en los casos en que estas divisiones se pueden realizar sin resto), es del mismo orden, pero la complejidad en bits es aproximadamente el tamaño en bits de las entradas originales en la matriz multiplicado pornorte{\displaystyle n}. [ 58 ]

Si ya se han calculado el determinante de A y el inverso de A , el lema del determinante de matrices permite calcular rápidamente el determinante de A + uv T , donde u y v son vectores columna.

Charles Dodgson (es decir, Lewis Carroll, el autor de Alicia en el País de las Maravillas) inventó un método para calcular determinantes llamado condensación de Dodgson . Este método no siempre funciona en su forma original. [ 59 ]

Véase también

Notas

  1. Lang 1985 , §VII.1
  2. "Determinantes y volúmenes" . textbooks.math.gatech.edu . Consultado el 16 de marzo de 2018 .
  3. McConnell (1957). Aplicaciones del análisis tensorial . Dover Publications. págs. 10–17 . 
  4. Harris 2014 , §4.7
  5. Serge Lang , Álgebra lineal , 2.ª edición, Addison-Wesley, 1971, págs. 173, 191.
  6. Lang 1987 , §VI.7, Teorema 7.5
  7. Alternativamente, Bourbaki 1998 , §III.8, Proposición 1 demuestra este resultado utilizando la funtorialidad de la potencia exterior.
  8. Horn y Johnson 2018 , §0.8.7
  9. Kung, Rota y Yan 2009 , pág. 306 
  10. Horn & Johnson 2018 , §0.8.2 .
  11. Silvester, JR (2000). "Determinantes de matrices de bloques" . Math . Gaz . 84 (501): 460– 467. doi : 10.2307/3620776 . JSTOR 3620776. S2CID 41879675 .  
  12. Sothanaphan, Nat (enero de 2017). "Determinantes de matrices de bloques con bloques no conmutativos". Álgebra lineal y sus aplicaciones . 512 : 202–218 . arXiv : 1805.06027 . doi : 10.1016/j.laa.2016.10.004 . S2CID 119272194 . 
  13. Bunch y Hopcroft 1974 .
  14. Las demostraciones se pueden encontrar en http://www.ee.ic.ac.uk/hp/staff/dmb/matrix/proof003.html
  15. Lin, Minghua; Sra, Suvrit (2014). "Superaditividad completamente fuerte de funciones matriciales generalizadas". arXiv : 1410.1958 [ math.FA ].
  16. Paksoy; Turkmen; Zhang (2014). "Desigualdades de funciones matriciales generalizadas mediante productos tensoriales" . Revista electrónica de álgebra lineal . 27 : 332–341 . doi : 10.13001/1081-3810.1622 .
  17. Serre, Denis (18 de octubre de 2010). "Concavidad de det 1 / n sobre HPD n " . MathOverflow .
  18. Lang 1985 , §VIII.2 , Horn & Johnson 2018 , Def. 1.2.3
  19. Horn y Johnson 2018 , Observación 7.1.2, Teorema 7.2.5
  20. Se puede encontrar una prueba en el Apéndice B de Kondratyuk, LA; Krivoruchenko, MI (1992). "Materia de quarks superconductores en el grupo de colores SU (2)". Zeitschrift für Physik A. 344 (1): 99– 115. Código Bib : 1992ZPhyA.344...99K . doi : 10.1007/BF01291027 . S2CID 120467300 . 
  21. Horn & Johnson 2018 , § 0.8.10 
  22. Grattan-Guinness 2003 , §6.6
  23. Cajori, F. Historia de las matemáticas, pág.  80
  24. 1 2 3 Campbell, H: "Álgebra lineal con aplicaciones", páginas 111–112. Appleton Century Crofts, 1971
  25. Eves 1990 , pág. 405 
  26. Breve historia del álgebra lineal y la teoría de matrices en: "Breve historia del álgebra lineal y la teoría de matrices" . Archivado del original el 10 de septiembre de 2012. Consultado el 24 de enero de 2012 .
  27. Kleiner 2007 , pág. 80 
  28. Bourbaki (1994 , p. 59) 
  29. Muir, Sir Thomas, La teoría de los determinantes en el orden histórico del desarrollo [Londres, Inglaterra: Macmillan and Co., Ltd., 1906]. JFM 37.0181.02 
  30. Kleiner 2007 , §5.2
  31. ^ El primer uso de la palabra "determinante" en el sentido moderno apareció en: Cauchy, Augustin-Louis "Memoire sur les fonctions qui ne peuvent obtenir que deux valeurs égales et des signes contraires par suite des transpositions operées entre les variables qu'elles renferment", que se leyó por primera vez en el Instituto de Francia en París el 30 de noviembre de 1812 y que se publicó posteriormente en el Journal de l'Ecole Polytechnique , Cahier 17, Tomo 10, páginas 29-112 (1815).
  32. Orígenes de los términos matemáticos: http://jeff560.tripod.com/d.html
  33. Historia de las matrices y los determinantes: http://www-history.mcs.st-and.ac.uk/history/HistTopics/Matrices_and_determinants.html
  34. Eves 1990 , pág. 494 
  35. Cajori 1993 , vol. II, pág. 92, núm. 462
  36. Historia de la notación matricial: http://jeff560.tripod.com/matrices.html
  37. Habgood y Arel 2012
  38. Lang 1985 , §VII.3
  39. Lang 2002 , §IV.8
  40. Lang 1985 , §VII.6, Teorema 6.10
  41. Lay, David (2021). Álgebra lineal y sus aplicaciones, 6.ª edición . Pearson. pág. 172. 
  42. ^ Dr. Jeff Cruza. «Área del triángulo» .
  43. Dummit & Foote 2004 , §11.4
  44. Dummit y Foote 2004 , §11.4, Teorema 30
  45. Mac Lane 1998 , §I.4 . Véase también Transformación natural §  Determinante .
  46. Bourbaki 1998 , §III.8
  47. Lombardi & Quitté 2015 , §5.2 , Bourbaki 1998 , §III.5
  48. Garibaldi 2004
  49. En un contexto no conmutativo, la linealidad izquierda (compatibilidad con la multiplicación izquierda por escalares) debe distinguirse de la linealidad derecha. Suponiendo que la linealidad en las columnas se considera linealidad izquierda, se tendría, para escalares no conmutativos a , b : ab=ab|1001|=a|100b|=|a00b|=b|a001|=ba|1001|=ba,{\displaystyle {\begin{aligned}ab&=ab{\begin{vmatrix}1&0\\0&1\end{vmatrix}}=a{\begin{vmatrix}1&0\\0&b\end{vmatrix}}\\[5mu]&={\begin{vmatrix}a&0\\0&b\end{vmatrix}}=b{\begin{vmatrix}a&0\\0&1\end{vmatrix}}=ba{\begin{vmatrix}1&0\\0&1\end{vmatrix}}=ba,\end{aligned}}} una contradicción. No existe una noción útil de funciones multilineales sobre un anillo no conmutativo.
  50. Varadarajan, V. S (2004), Supersimetría para matemáticos: Una introducción , American Mathematical Soc., ISBN 978-0-8218-3574-6.
  51. "... mencionamos que el determinante, aunque es una noción conveniente desde el punto de vista teórico, rara vez encuentra un papel útil en los algoritmos numéricos.", véase Trefethen y Bau III 1997 , Lección 1 .
  52. Fisikopoulos y Peñaranda 2016 , §1.1, §4.3
  53. Camarero, Cristóbal (2018-12-05). "Algoritmos simples, rápidos y prácticos para la descomposición de Cholesky, LU y QR mediante la multiplicación rápida de matrices rectangulares". arXiv : 1812.02056 [ cs.NA ].
  54. Bunch y Hopcroft 1974
  55. Fisikopoulos y Peñaranda 2016 , §1.1
  56. Rote 2001
  57. Fang, Xin Gui; Havas, George (1997). "Sobre la complejidad del peor caso de la eliminación gaussiana entera" (PDF) . Actas del simposio internacional de 1997 sobre computación simbólica y algebraica . ISSAC '97. Kihei, Maui, Hawái, Estados Unidos: ACM. págs. 28–31 . doi : 10.1145/258726.258740 . ISBN  0-89791-875-4Archivado del original (PDF) el 7 de agosto de 2011. Consultado el 22 de enero de 2011 .
  58. Fisikopoulos & Peñaranda 2016 , §1.1 , Bareiss 1968
  59. Abeles, Francine F. (2008). "Condensación de Dodgson: El desarrollo histórico y matemático de un método experimental" . Álgebra lineal y sus aplicaciones . 429 ( 2– 3): 429– 438. doi : 10.1016/j.laa.2007.11.022 .

Referencias

  • Anton, Howard (2005), Álgebra lineal elemental (versión de aplicaciones) (9.ª  ed.), Wiley International
  • Axler, Sheldon Jay (2015). Álgebra lineal bien hecha (3.ª  ed.). Springer . ISBN 978-3-319-11079-0.
  • Bareiss, Erwin (1968), "Identidad de Sylvester y eliminación gaussiana multietapa que preserva los enteros" (PDF) , Mathematics of Computation , 22 (102): 565–578 , doi : 10.2307/2004533 , JSTOR 2004533 , archivado (PDF) del original el 25 de octubre de 2012. 
  • de Boor, Carl (1990), "Un ejercicio vacío" (PDF) , ACM SIGNUM Newsletter , 25 (2): 3–7 , doi : 10.1145/122272.122273 , S2CID 62780452 , archivado (PDF) del original el 1 de septiembre de 2006 
  • Bourbaki, Nicolas (1998), Álgebra I, Capítulos 1-3 , Springer, ISBN 9783540642435
  • Bunch, James R.; Hopcroft, John E. (1974). "Factorización triangular e inversión mediante multiplicación rápida de matrices" . Matemáticas de la computación . 28 (125): 231– 236. doi : 10.1090/S0025-5718-1974-0331751-8 . hdl : 1813/6003 .
  • Dummit, David S.; Foote, Richard M. (2004), Álgebra abstracta (3.ª  ed.), Hoboken, NJ: Wiley, ISBN 9780471452348, OCLC 248917264 
  • Fisikopoulos, Vissarion; Peñaranda, Luis (2016), "Algoritmos geométricos más rápidos mediante el cálculo dinámico de determinantes", Geometría Computacional , 54 : 1–16 , arXiv : 1206.7067 , doi : 10.1016/j.comgeo.2015.12.001
  • Garibaldi, Skip (2004), "El polinomio característico y el determinante no son construcciones ad hoc", American Mathematical Monthly , 111 (9): 761–778 , arXiv : math/0203276 , doi : 10.2307/4145188 , JSTOR 4145188 , MR 2104048  
  • Habgood, Ken; Arel, Itamar (2012). "Una aplicación de la regla de Cramer basada en la condensación para resolver sistemas lineales a gran escala" (PDF) . Journal of Discrete Algorithms . 10 : 98–109 . doi : 10.1016/j.jda.2011.06.007 . Archivado (PDF) del original el 5 de mayo de 2019.
  • Harris, Frank E. (2014), Matemáticas para las ciencias físicas y la ingeniería , Elsevier, ISBN 9780128010495
  • Kleiner, Israel (2007), Kleiner, Israel (ed.), Una historia del álgebra abstracta , Birkhäuser, doi : 10.1007/978-0-8176-4685-1 , ISBN 978-0-8176-4684-4, MR 2347309 
  • Kung, Joseph PS; Rota, Gian-Carlo; Yan, Catherine (2009), Combinatoria: The Rota Way , Cambridge University Press, ISBN 9780521883894
  • Lay, David C. (22 de agosto de 2005), Álgebra lineal y sus aplicaciones (3.ª  ed.), Addison Wesley, ISBN 978-0-321-28713-7
  • Lombardi, Enrique; Quitté, Claude (2015), Álgebra conmutativa: métodos constructivos , Springer, ISBN 9789401799447
  • Mac Lane, Saunders (1998), Categorías para el matemático en activo , Textos de posgrado en matemáticas 5 (2.ª  ed.), Springer-Verlag, ISBN 0-387-98403-8
  • Meyer, Carl D. (15 de febrero de 2001), Análisis matricial y álgebra lineal aplicada , Sociedad de Matemáticas Industriales y Aplicadas (SIAM), ISBN 978-0-89871-454-8Archivado del original el 31 de octubre de 2009.
  • Muir, Thomas (1960) [1933], Tratado sobre la teoría de los determinantes , Revisado y ampliado por William H. Metzler, Nueva York, NY: Dover
  • Poole, David (2006), Álgebra lineal: una introducción moderna (2.ª  ed.), Brooks/Cole, ISBN 0-534-99845-3
  • G. Baley Price (1947) "Algunas identidades en la teoría de determinantes", American Mathematical Monthly 54:75–90 MR 0019078 
  • Horn, Roger Alan ; Johnson, Charles Royal (2018) [1985]. Análisis matricial (2.ª  ed.). Cambridge University Press . ISBN 978-0-521-54823-6.
  • Lang, Serge (1985), Introducción al álgebra lineal , Textos de matemáticas para estudiantes de pregrado (2.ª  ed.), Springer, ISBN 9780387962054
  • Lang, Serge (1987), Álgebra lineal , Textos de matemáticas para estudiantes de pregrado (3.ª  ed.), Springer, ISBN 9780387964126
  • Lang, Serge (2002). Álgebra . Textos de posgrado en matemáticas. Nueva York, NY: Springer. ISBN 978-0-387-95385-4.
  • Leon, Steven J. (2006), Álgebra lineal con aplicaciones (7.ª  ed.), Pearson Prentice Hall
  • Rote, Günter (2001), "Algoritmos sin división para el determinante y el pfaffiano: enfoques algebraicos y combinatorios" (PDF) , Matemáticas discretas computacionales , Lecture Notes in Comput. Sci., vol.  2122, Springer, pp. 119–135 , doi : 10.1007/3-540-45506-X_9 , ISBN  978-3-540-42775-9MR 1911585 , archivado del original (PDF) el 1 de febrero de 2007 , consultado el 4 de junio de 2020. 
  • Trefethen, Lloyd ; Bau III, David (1997), Álgebra lineal numérica (1.ª  ed.), Filadelfia: SIAM, ISBN 978-0-89871-361-9

Referencias históricas

  • Bourbaki, Nicolas (1994), Elementos de la historia de las matemáticas , traducido por Meldrum, John , Springer, doi : 10.1007/978-3-642-61693-8 , ISBN 3-540-19376-6
  • Cajori, Florian (1993), Historia de las notaciones matemáticas: Incluye Vol. I. Notaciones en matemáticas elementales; Vol. II. Notaciones principalmente en matemáticas superiores, Reimpresión de los originales de 1928 y 1929 , Dover, ISBN 0-486-67766-4, MR 3363427 
  • Bézout, Étienne (1779), Théorie générale des ecuaciones algébriques , París
  • Cayley, Arthur (1841), "Sobre un teorema en la geometría de la posición", Cambridge Mathematical Journal , 2 : 267–271
  • Cramer, Gabriel (1750), Introducción al análisis de las líneas courbes algébriques , Ginebra: Frères Cramer & Cl. Filiberto, doi : 10.3931/e-rara-4048
  • Eves, Howard (1990), Introducción a la historia de las matemáticas (6.ª  ed.), Saunders College Publishing, ISBN 0-03-029558-0, MR 1104435 
  • Grattan-Guinness, I., ed. (2003), Companion Encyclopedia of the History and Philosophy of the Mathematical Sciences , vol.  1, Johns Hopkins University Press , ISBN 9780801873966
  • Jacobi, Carl Gustav Jakob (1841), "De Determinantibus funcionalibus" , Journal für die reine und angewandte Mathematik , 1841 (22): 320– 359, doi : 10.1515/crll.1841.22.319 , S2CID 123637858 
  • Laplace, Pierre-Simon, de (1772), "Recherches sur le calcul intégral et sur le systéme du monde" , Histoire de l'Académie Royale des Sciences (segunda parte), París: 267– 376{{citation}}: CS1 maint: varios nombres: lista de autores ( enlace )
  • Robert Forsyth Scott (1880): Tratado sobre la teoría de los determinantes y sus aplicaciones en análisis y geometría , Cambridge University Press
  • ER Hedrick: Sobre determinantes tridimensionales , Anales de Matemáticas, vol. 1, n.º 1/4 (1899-1900), págs. 49-67 (19 páginas). https://doi.org/10.2307/1967268 # Nota: Este no es el determinante ordinal.