矩阵的逆

最新推荐文章于 2024-07-23 15:29:26 发布

水蛙菌

最新推荐文章于 2024-07-23 15:29:26 发布

阅读量1.5w

点赞数 7

分类专栏：其他文章标签：线性代数数据结构矩阵的逆 c++

本文链接：https://blog.csdn.net/qq_44826711/article/details/115294447

版权

矩阵运算逆矩阵初等行变换时间复杂度数值计算

关键词由CSDN通过智能技术生成

其他专栏收录该内容

9 篇文章 4 订阅

订阅专栏

给一个矩阵，求它的逆。

将原矩阵A和一个单位矩阵E作成大矩阵(A, E)，用初等行变换将大矩阵中的A变为E，则会得到(E, A^-1)的形式。

时间复杂度：O(n³)

inline vector<double> operator * (vector<double> a, double b)
{
	int N = a.size();
	vector<double> res(N, 0);
	for(int i = 0; i < N; ++i)
		res[i] = a[i] * b;
	return res;
}

inline vector<double> operator - (vector<double> a, vector<double> b)
{
	int N = a.size();
	vector<double> res(N, 0);
	for(int i = 0; i < N; ++i)
		res[i] = a[i] - b[i];
	return res;
}

inline void inverse(vector<double> A[], vector<double> C[], int N)
{
	for(int i = 0; i < N; ++i)
		C[i] = vector<double>(N, 0);
	for(int i = 0; i < N; ++i)
		C[i][j] = 1;
	for(int i = 0; i < N; ++i)
	{
		for(int j = i; j < N; ++j)
			if(fabs(A[i][j]) > 0)
			{
				swap(A[i], A[j]);
				swap(C[i], C[j]);
				break;
			}
		C[i] = C[i] * (1 / A[i][i]);
		A[i] = A[i] * (1 / A[i][i]);
		for(int j = 0; j < N; ++j)
			if(j != i && fabs(A[i][j] > 0))
			{
				C[j] = C[j] - C[i] * A[j][i];
				A[j] = A[j] - A[i] * A[j][i];
			}
	}
}