En álgebra lineal, la fórmula de Leibniz expresa el determinante de una matriz cuadrada en términos de permutaciones de los elementos de la matriz. La fórmula debe su nombre a Gottfried Leibniz, la fórmula para una matriz de orden
es:
![{\displaystyle \det(A)=\sum _{\sigma \in S_{n}}\operatorname {sgn}(\sigma )\prod _{i=1}^{n}a_{\sigma (i),i},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f01618d6347d0680413998a0295a7b55aaa9e28)
donde
![{\displaystyle A=(a_{ij})_{i,j=1,\dots ,n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d83da6fd10359ef6906613b27f6c19f27dbb30fd)
y donde sgn es la función signo de permutaciones en el grupo de permutación Sn, que devuelve +1 si la permutación es par y −1 si es impar.
Otra notación común usada para la fórmula utiliza símbolos de Levi-Civita y la notación de Einstein, quedando:
![{\displaystyle \det(A)=\epsilon ^{i_{1}\cdots i_{n}}{a}_{1i_{1}}\cdots {a}_{ni_{n}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/25a332b7883c7f44f878247c17a26e870516d265)
que puede ser más familiar para los físicos.
Evaluar directamente la fórmula de Leibniz requiere
operaciones en general. Esto es, necesita un número de operaciones asintóticamente proporcional a n factorial, ya que n! es el número de permutaciones de orden n, lo que resulta aparatoso para valores grandes de n. En su lugar, el determinante se puede evaluar en O(n3) operaciones mediante la descomposición LU de la matriz (
, normalmente a través de la eliminación gaussiana o métodos similares). En ese caso,
, y los determinantes de las matrices triangulares L y U serán los productos de las entradas de sus respectivas diagonales principales. En la práctica de álgebra lineal, sin embargo, rara vez se requiere el cálculo explícito del determinante. Ver, por ejemplo, Trefethen y Bau (1997).
Teorema. Existe exactamente una función multilineal alternada:
![{\displaystyle F:{\mathfrak {M}}_{n}(\mathbb {K} )\rightarrow \mathbb {K} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6559e4230eb9ddb355446bd4487fb323fc98d98e)
tal que
.
Prueba.
Existencia: La función
, definida por la fórmula de Leibniz, debe cumplir estas tres propiedades:
- Multilinealidad (esto es, conserva la suma y el producto por escalar en cada componente):
.
.
- Función alternada (o sea, que si dos componentes son iguales, la función se anula).
.
- Para cualquier
, sea
la permutación igual a
salvo porque intercambia las imágenes de
y
. ![{\displaystyle {\begin{aligned}F(A)&=\sum _{\sigma \in S_{n},\sigma (j_{1})<\sigma (j_{2})}\left[\operatorname {sgn}(\sigma )\left(\prod _{i=1,i\neq j_{1},i\neq j_{2}}^{n}a_{\sigma (i)}^{i}\right)a_{\sigma (j_{1})}^{j_{1}}a_{\sigma (j_{2})}^{j_{2}}+\operatorname {sgn}(\sigma ')\left(\prod _{i=1,i\neq j_{1},i\neq j_{2}}^{n}a_{\sigma '(i)}^{i}\right)a_{\sigma '(j_{1})}^{j_{1}}a_{\sigma '(j_{2})}^{j_{2}}\right]\\&=\sum _{\sigma \in S_{n},\sigma (j_{1})<\sigma (j_{2})}\left[\operatorname {sgn}(\sigma )\left(\prod _{i=1,i\neq j_{1},i\neq j_{2}}^{n}a_{\sigma (i)}^{i}\right)a_{\sigma (j_{1})}^{j_{1}}a_{\sigma (j_{2})}^{j_{2}}-\operatorname {sgn}(\sigma )\left(\prod _{i=1,i\neq j_{1},i\neq j_{2}}^{n}a_{\sigma (i)}^{i}\right)a_{\sigma (j_{2})}^{j_{1}}a_{\sigma (j_{1})}^{j_{2}}\right]\\&=\sum _{\sigma \in S_{n},\sigma (j_{1})<\sigma (j_{2})}\operatorname {sgn}(\sigma )\left(\prod _{i=1,i\neq j_{1},i\neq j_{2}}^{n}a_{\sigma (i)}^{i}\right)\left(a_{\sigma (j_{1})}^{j_{1}}a_{\sigma (j_{2})}^{j_{2}}-a_{\sigma (j_{1})}^{j_{2}}a_{\sigma (j_{2})}^{j_{_{1}}}\right)\\\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b83899b02158f013dde5b363b838368208beb43)
- Por lo tanto, si
entonces
.
- Finalmente,
:
.
Singularidad: Sea
una función de ese tipo, y sea
una matriz
. Llámese
a la
-ésima columna de
, i.e.
de modo que
También, sea
la
-ésima columna de la matriz de identidad, en forma de vector.
Ahora se escribe cada uno de los vectores
en términos de los vectores
, por ejemplo:
.
Como
es multilineal, se tiene
![{\displaystyle {\begin{aligned}F(A)&=F\left(\sum _{k_{1}=1}^{n}a_{k_{1}}^{1}E^{k_{1}},\dots ,\sum _{k_{n}=1}^{n}a_{k_{n}}^{n}E^{k_{n}}\right)\\&=\sum _{k_{1},\dots ,k_{n}=1}^{n}\left(\prod _{i=1}^{n}a_{k_{i}}^{i}\right)F\left(E^{k_{1}},\dots ,E^{k_{n}}\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f208fe6e2dd47425a643af9db402e6dd4a380ba6)
Ahora bien, como
es alternada, cualquier combinación con índices repetidos es cero; por lo que el sumatorio puede reducirse a las tuplas con índices no repetidos, es decir, a únicamente las permutaciones:
![{\displaystyle F(A)=\sum _{\sigma \in S_{n}}\left(\prod _{i=1}^{n}a_{\sigma (i)}^{i}\right)F(E^{\sigma (1)},\dots ,E^{\sigma (n)}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f050fe9e3712b950e0b09c0be08b08d047f4e7a)
donde
es el grupo simétrico de los primeros n enteros (es decir, el conjunto de todas las permutaciones de los primeros n enteros). Como la función
es alternada, el orden de los vectores
sólo afecta para el signo del resultado, de forma que se puede extraer la permutación por la función signo. Como
:
![{\displaystyle {\begin{aligned}F(A)&=\sum _{\sigma \in S_{n}}\operatorname {sgn}(\sigma )\left(\prod _{i=1}^{n}a_{\sigma (i)}^{i}\right)F(I)\\&=\sum _{\sigma \in S_{n}}\operatorname {sgn}(\sigma )\prod _{i=1}^{n}a_{\sigma (i)}^{i},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a8cdcb815eb8b595634ab99158c96ce6a81a0ee)
que es precisamente la función definida por la fórmula de Leibniz.
Véase también
- Matriz.
- Regla de Sarrus.
- Teorema de Laplace para el cálculo de determinantes.
- Regla de Cramer.
Referencias
- Lloyd N. Trefethen y David Bau, Numerical Linear Algebra (SIAM, 1997) ISBN 978-0898713619
Control de autoridades | - Proyectos Wikimedia
Datos: Q538557 |
---|
Datos: Q538557