重温数学基础——矩阵求逆

changshen_xu

墨衍会员 · AI 创作全网分发码龄9年

关注

原创最新推荐文章于 2026-07-30 13:04:42 发布 8.5k 阅读 AI 写同款· GEO 优化 ›

线代——求逆矩阵的快捷方法分块矩阵求逆矩阵阅读详情

目标：给定任意一个矩阵A，求它的逆矩阵。

方案一：直接求逆

【例1】：简单的二维矩阵

设 $A=[abcd]A=\begin{bmatrix} a & b \\ c & d\end{bmatrix}$ ，则 $A−1=1det(A)[d−b−ca]A^{-1}=\frac{1}{det(A)} \begin{bmatrix} d & -b \\ -c & a\end{bmatrix}$ 。
【编程实现】：
注意1.det(A)=ad-bc的值等于0时，矩阵A没有逆；
注意2.det(A)=ad-bc的值非常接近于0时，求得的逆 $A^{-1}$ 可能是不准确的（由于一些数值舍入误差？）。

int QRInvM2(double A[2][2]) {
	double detA = A[0][0] * A[1][1] - A[0][1] * A[1][0];
	if ((detA < 1.0e-20) && (detA > -1.0e-20))
		return(0);
	double t = 1/detA ;
	A[0][0] = A[0][0] * t;
	A[1][1] = A[1][1] * t;
	A[0][1] = -A[0][1] * t;
	A[1][0] = -A[1][0] * t;
	t = A[0][0];
	A[0][0] = A[1][1];
	A[1][1] = t;
	printf("%lf\n", detA);
	printf("%lf, %lf\n", A[0][0], A[0][1]);
	printf("%lf, %lf\n", A[1][0], A[1][1]);
	return(1);
}

【例2】：特殊的上三角矩阵
设矩阵 $A=[a00a01a02a030a11a12a1300a22a23000a33]A=\begin{bmatrix} a_{00} & a_{01} & a_{02} & a_{03} \\ 0 & a_{11} & a_{12} & a_{13} \\ 0 & 0 & a_{22} & a_{23} \\ 0 & 0 & 0 & a_{33} \\ \end{bmatrix}$ ，同时设它的逆矩阵 $A−1=[a00′a01′a02′a03′a10′a11′a12′a13′a20′a21′a22′a23′a30′a31′a32′a33′]A^{-1}=\begin{bmatrix} a'_{00} & a'_{01} & a'_{02} & a'_{03} \\ a'_{10} & a'_{11} & a'_{12} & a'_{13} \\ a'_{20} & a'_{21} & a'_{22} & a'_{23} \\ a'_{30} & a'_{31} & a'_{32} & a'_{33} \\ \end{bmatrix}$ ，根据 $AA^{-1}=I$ ，将 $I$ 第一列对应的四个方程取出有：
$a00∗a00′+a01∗a10′+a02∗a20′+a03∗a30′=1a11∗a10′+a12∗a20′+a13∗a30′=0a22∗a20′+a23∗a30′=0a33∗a30′=0\begin{aligned} a_{00}*a'_{00}+a_{01}*a'_{10}+a_{02}*a'_{20}+a_{03}*a'_{30} & =1 \\ a_{11}*a'_{10}+a_{12}*a'_{20}+a_{13}*a'_{30} &=0 \\ a_{22}*a'_{20}+a_{23}*a'_{30} &=0 \\ a_{33}*a'_{30} &=0 \end{aligned}$
则，（四个方程四个未知量）我们可以求出 $A^{-1}$ 的第一列。后面三列求取的思路与此类似。
【编程实现】：

int UInvM4(double U[4][4]) {
	double detU = U[0][0] * U[1][1] * U[2][2] * U[3][3];
	if ((detU < 1.0e-20) && (detU > -1.0e-20))
		return(0);
	double invU[4][4];
	for (int i = 0; i < 4; i++) {
		for (int j = 0; j < 4; j++) {
			invU[i][j] = 0.0;
		}
	}
	double b[4];
	for (int k = 0; k < 4; k++) {
		for (int i = 0; i < 4; i++) {
			b[i] = 0.0;
		}
		b[k] = 1.0;
		invU[3][k] = b[3] / U[3][3];
		invU[2][k] = (b[2] - U[2][3] * invU[3][k]) / U[2][2];
		invU[1][k] = (b[1] - U[1][3] * invU[3][k] - U[1][2] * invU[2][k]) / U[1][1];
		invU[0][k] = (b[0] - U[0][3] * invU[3][k] - U[0][2] * invU[2][k] - U[0][1] * invU[1][k]) / U[0][0];
	}
	return(1);
}

方案二：LU分解

收录于