c++ 求任意阶矩阵行列式代码

最新推荐文章于 2024-02-04 21:42:24 发布

原创最新推荐文章于 2024-02-04 21:42:24 发布 · 1.7k 阅读

4 ·

CC 4.0 BY-SA版权

c++ 专栏收录该内容

72 篇文章

订阅专栏

本文介绍了如何使用定义法计算双精度矩阵的行列式，通过递归地构造子矩阵并应用代数余子式，详细展示了从单元素到复杂矩阵的求解过程。

double matrix::determinant() { // calculate the determinant of the matrix
	//使用定义法计算行列式
	//如果只有一个元素，则返回当前元素
	if (n == 1) {
		return values[0];
	}
	else {
		//对第一行使用代数余子式的方式计算
		double sub = 0;
		for (int i = 0; i < n; i++) {
			//初始化一个子矩阵
			int *subMatrix = new int[(n - 1)*(n - 1)];
			int index = 0;
			for (int j = n; j < n*n; j++) {
				int col = j % n;
				if (col != i % n) {
					subMatrix[index++] = values[j];
				}
			}
			matrix subm(subMatrix, n - 1);
			double result = subm.determinant();
			sub += pow(-1, 2 + i)*result*values[i];
		}
		return sub;
	}
}