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

最新推荐文章于 2026-03-02 09:15:45 发布

原创最新推荐文章于 2026-03-02 09:15:45 发布 · 8.4k 阅读

5 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#矩阵求逆

▶ 数学基础温习专栏收录该内容

6 篇文章

订阅专栏

博客围绕给定矩阵A求其逆矩阵展开，介绍了两种方案。方案一是直接求逆，以简单二维矩阵和特殊上三角矩阵为例说明，并提及编程实现时的注意事项；方案二则是LU分解。

目标：给定任意一个矩阵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分解