The Algorithms logo
The Algorithms
À proposFaire un don

Determinant

P
H
/**
 * Given a square matrix, find its determinant using Laplace Expansion.
 * Time Complexity : O(n!)
 *
 * For more info: https://en.wikipedia.org/wiki/Determinant
 *
 * @param {number[[]]} matrix - Two dimensional array of integers.
 * @returns {number} - An integer equal to the determinant.
 *
 * @example
 * const squareMatrix = [
 *                          [2,3,4,6],
 *                          [5,8,9,0],
 *                          [7,4,3,9],
 *                          [4,0,2,1]
 *                      ];
 *
 * const result = determinant(squareMatrix);
 * // The function should return 858 as the resultant determinant.
 */

const subMatrix = (matrix, i, j) => {
  let matrixSize = matrix[0].length
  if (matrixSize === 1) {
    return matrix[0][0]
  }
  let subMatrix = []
  for (let x = 0; x < matrixSize; x++) {
    if (x === i) {
      continue
    }
    subMatrix.push([])
    for (let y = 0; y < matrixSize; y++) {
      if (y === j) {
        continue
      }
      subMatrix[subMatrix.length - 1].push(matrix[x][y])
    }
  }
  return subMatrix
}

const isMatrixSquare = (matrix) => {
  let numRows = matrix.length
  for (let i = 0; i < numRows; i++) {
    if (numRows !== matrix[i].length) {
      return false
    }
  }
  return true
}

const determinant = (matrix) => {
  if (
    !Array.isArray(matrix) ||
    matrix.length === 0 ||
    !Array.isArray(matrix[0])
  ) {
    throw new Error('Input is not a valid 2D matrix.')
  }
  if (!isMatrixSquare(matrix)) {
    throw new Error('Square matrix is required.')
  }
  let numCols = matrix[0].length
  if (numCols === 1) {
    return matrix[0][0]
  }
  let result = 0
  let setIndex = 0
  for (let i = 0; i < numCols; i++) {
    result +=
      Math.pow(-1, i) *
      matrix[setIndex][i] *
      determinant(subMatrix(matrix, setIndex, i))
  }
  return result
}
export { determinant }
À propos de cet Algorithme

Description

Principe

  • Calculer le déterminant d'une matrice permet de savoir si elle est inversible ou non.
  • Avoir cette information en elle-même ne sert pas à grand-chose, mais si cette matrice représente (par exemple) un système d'équation, cela permet de savoir si ce système possède une solution.

Complexités

  • L'algorithme possède une complexité temporelle de l'ordre de $n!$ (À vérifier).
  • On utilise une fonction récursive, ce qui est gourmand en termes de mémoire.

Applications

  • Trouver le déterminant d'une matrice permet de savoir si on peut utiliser la méthode du pivot de Gauss ou pas.

Étapes

  • Puisque le déterminant n'est défini que pour les matrices carrées, il faut d'abord vérifier que la matrice est carré.
  • Ensuite, nous développons la matrice par rapport à la première ligne.
  • On répète ces étapes à la matrice ainsi développer jusqu'à tomber sur une matrice 2×2, où on applique la formule connue $\bigg(\begin{vmatrix}a&b\c&d\end{vmatrix}=ad-cb\bigg)$.
  • De là, on remonte la cascade d'appel en effectuant les calculs.

Exemple

Calculons le déterminant de la matrice $\begin{bmatrix}1&2&3\4&5&6\7&8&9\end{bmatrix}$.

Exemple en vidéo

Conclusion

Le déterminant de $\begin{bmatrix}1&2&3\4&5&6\7&8&9\end{bmatrix}$ est $0$.