C++ 程序中的矩阵行列式


本教程中,我们将学习如何求矩阵行列式。

让我们来看看求矩阵行列式的步骤。

  • 初始化矩阵。

  • 编写一个函数以求矩阵行列式。

    • 如果矩阵大小为 1 或 2,则求矩阵行列式。这是一件很简单的事情。

    • 初始化行列式、子矩阵、符号的变量。

    • 从 1 迭代到矩阵大小 N。

    • 找到当前矩阵元素的子矩阵。

      • 不在当前元素行和列中的所有元素

    • 将当前元素及其余因子的乘积加到行列式中。

    • 更改符号。

  • 打印矩阵行列式。

示例

我们来看一下代码。

 演示

#include <bits/stdc++.h>
using namespace std;
#define N 3
void subMatrix(int mat[N][N], int temp[N][N], int p, int q, int n) {
   int i = 0, j = 0;
   // filling the sub matrix
   for (int row = 0; row < n; row++) {
      for (int col = 0; col < n; col++) {
         // skipping if the current row or column is not equal to the current
         // element row and column
         if (row != p && col != q) {
            temp[i][j++] = mat[row][col];
            if (j == n - 1) {
               j = 0;
               i++;
            }
         }
      }
   }
}
int determinantOfMatrix(int matrix[N][N], int n) {
   int determinant = 0;
   if (n == 1) {
      return matrix[0][0];
   }
   if (n == 2) {
      return (matrix[0][0] * matrix[1][1]) - (matrix[0][1] * matrix[1][0]);
   }
   int temp[N][N], sign = 1;
   for (int i = 0; i < n; i++) {
      subMatrix(matrix, temp, 0, i, n);
      determinant += sign * matrix[0][i] * determinantOfMatrix(temp, n - 1);
      sign = -sign;
   }
   return determinant;
}
int main() {
   int mat[N][N] = {{2, 1, 3}, {6, 5, 7}, {4, 9, 8}};
   cout << "Determinant: " << determinantOfMatrix(mat, N) << endl;
   return 0;
}

输出

如果你执行以上程序,你将得到以下结果。

Determinant: 36

结论

如果你对教程有任何疑问,请在评论部分提到。

更新于: 2021-01-27

6K+ 浏览

开启你的职业生涯

完成课程并获得认证

开始学习
广告