c++解算3*3矩阵的逆矩阵

LiChengwei47

于 2020-12-17 18:33:32 发布

阅读量3.1k

点赞数

CC 4.0 BY-SA版权

分类专栏： C++&Cmake个人学习记录文章标签： c++

本文链接：https://blog.youkuaiyun.com/it_LiChengwei/article/details/111321566

用伴随矩阵计算3行3列矩阵的逆矩阵（要求矩阵行列式不能为0）
原理参考：https://blog.youkuaiyun.com/feixia_24/article/details/41644335
矩阵子函数部分

double** MatrixInversion3_3(double A[3][3])
{
	double **Result = (double**)malloc(sizeof(double*) * 3);
	Result[0] = (double*)malloc(sizeof(double*) * 3);
	Result[1] = (double*)malloc(sizeof(double*) * 3);
	Result[2] = (double*)malloc(sizeof(double*) * 3);
	double a1 = *A[0]; double b1 = *(A[0]+1); double c1 = *(A[0]+2);
	double a2 = *A[1]; double b2 = *(A[1] + 1); double c2 = *(A[1] + 2);
	double a3 = *A[2]; double b3 = *(A[2] + 1); double c3 = *(A[2] + 2);
	double k = 1 / (a1*(b2*c3 - c2 * b3) - a2 * (b1*c3 - c1 * b3) + a3 * (b1*c2 - b2 * c1));
	*Result[0] = k * (b2*c3 - c2 * b3);
	*(Result[0]+1) = k * (c1*b3 - b1 * c3);
	*(Result[0] + 2) = k * (b1*c2 - c1 * b2)