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
结论
如果你对教程有任何疑问,请在评论部分提到。
广告