C++矩阵运算2

        矩阵的行列式计算，很多朋友们在学习的时候基本上只记2阶或者3阶行列式如何进行运算，因为考试的时候基本上就是考这些，3阶以上的行列式的计算的量确实太大了，不过这个事情我们可以交给计算机来做，只要将公式转换为代码那么无论是几阶的行列式都可以顺利的进行计算了。

1、行列是计算公式

2、公式说明

是关于，(n是行列式的列数或者式行数)的全排序的其中一种排序，这个求和公式需要遍历所有的排列情况，计算的是排序的逆序数个数。

生成全排序

以下图为列

        我们可以将某个数放在队列的队首或者是队尾，然后将其余的数视为一个完整的队列重复的进行这个操作，知道队列当中只剩下一个数为止，这样我们就得到了全排序，并且全排序的数量

具体代码如下:

std::vector<std::vector<int>> MathUntil::MakeTotalSort(const std::vector<int>& iNumbers) {
	std::vector<std::vector<int>> oTotalSort;
	oTotalSort.reserve(Factorial(iNumbers.size()));

	if (iNumbers.size() == 1)
		oTotalSort.emplace_back(iNumbers);
	
	for (int i = 0; i < iNumbers.size(); i++) {
		int lastNumber = *(iNumbers.end() - 1 - i);
		std::vector<int> Numbers = iNumbers;
		Numbers.erase(Numbers.end() - 1 - i);
		std::vector<std::vector<int>> TotalSort = MakeTotalSort(Numbers);
		for (int k = 0; k < TotalSort.size(); k++)
			TotalSort[k].emplace_back(lastNumber);
		
		oTotalSort.insert(oTotalSort.end(), TotalSort.begin(), TotalSort.end());
	}

	return oTotalSort;
}

计算逆序数

        在一个数的队列当中如果某个数前面存在有比它大的数，那么这个队列就存在逆序数。如何计算逆序数：从前到后遍历整个数列，某个数前面存在n个比它大的数，那么此时逆序数的个数为n，下一个数前面存在m个比它大的数那么此时的逆序数个数为n+m，以此类推直到遍历完所有的数为止。比如说：

序列：5 1 5 2

逆序数：0 1 0 2

逆序数的个数：1 + 2 = 3

代码如下：

int  MathUntil::ComputeInverseNumber(const std::vector<int>& iNumbers) {
	if (iNumbers.size() < 2)
		return 0;
	
	int InverseQuantities = 0;
	for (int i = 1; i < iNumbers.size(); i++) {
		int number = iNumbers[i];

		for (int j = i - 1; j >= 0; j--)
			if (number < iNumbers[j])
				InverseQuantities += 1;
	}

	return InverseQuantities;
}

行列式计算

具体代码如下

template<typename T>
T Matrix<T>::Det(const Matrix<T>& iMatrix) {
	MATRIXISNULL(iMatrix);MATRIX_ROW_EQUAL_COL(iMatrix);

	std::vector<int> indexNumbers;
	indexNumbers.reserve(iMatrix.mRow);

	for (int i = 0; i < iMatrix.mRow; i++)
		indexNumbers.emplace_back(i);

	std::vector<std::vector<int>> indexesNumbers = MathUntil::MakeTotalSort(indexNumbers);

	T result = static_cast<T>(0.0f);

	for (auto& indexes : indexesNumbers) {
		T mult = static_cast<T>(1.0f);

		for (int i = 0; i < iMatrix.mRow; i++) 
			mult *= iMatrix.GetElement(i, indexes[i]);
		
		mult *= pow(-1.0, MathUntil::ComputeInverseNumber(indexes));

		result += mult;
	}

	return result;
}

如果想要这些代码请移步到此网站https://gitee.com/HonyOrange_227/matrix-compute

直接复制下面三个文件即可

3、使用行列式计算

代码如下

#include"Matrix.h"
#include"MatrixUtil.h"

int main() {
	int arr[] = { 23,54,2,	22,23,3,	56,11,45};
	Matrix<int> m1(arr, 3, 3);

	ThrowIfFailed(std::cout << Matrix<int>::Det(m1) << std::endl);
}

计算结果，-23434

4、验证一下

第三方验证没有问题