求逆矩阵

上篇文章求解了矩阵的行列式，链接为：点击打开链接。

有了行列式，可以求解矩阵的逆矩阵，逆矩阵可以表示为 A^(-1) = A^* / |A|，其中  A^* 为伴随矩阵，|A| 为行列式。

	/**计算逆矩阵，求逆矩阵必须是n*n，因此放在行列式中求解
	 * @param A 要求解的矩阵
	 * @return 逆矩阵
	 */
	public double[][] inverse_Matrix(double A[][]){
		double result = get_Result(A);
		if(result == 0){
			System.out.println("行列式为0，不符合计算要求");
			return null;
		}
		double B[][] = new double[A.length][A.length];
		// 求伴随矩阵
		B = adjoint_Matrix(A.length);
		//求逆矩阵
		for(int i = 0; i < B.length; i++){
			for(int j = 0; j < B.length; j++){
				B[i][j] = B[i][j]/result;
			}
		}
		return B;
	}

在这里没有判断矩阵的行列是否相等，因为已经在行列式求解，也就是函数get_Result（A）中进行了判断。

现在需要求解伴随矩阵，本例通过 adjoint_Matrix()求解，传入参数为A数组的长度，伴随矩阵的求解方法可以参照该链接的内容：点击打开链接，相关实现代码如下：

	/**求伴随矩阵
	 * @param length 矩阵的行数
	 * @return 伴随矩阵
	 */
	private double[][] adjoint_Matrix(int length){
		double B[][] = new double[length][length];
		//如果是二阶
		if(length == 2){
			B[0][0] = A[1][1];
			B[0][1] = - A[1][0];
			B[1][0] = - A[0][1];
			B[1][1] = A[0][0];
			return B;
		}
		double sub_A[][] = new double[length - 1][length - 1];
		for(int i = 0; i < length; i++){
			for(int j = 0; j < length; j++){
				//获取当前行列式对应元素的子矩阵
				sub_A = get_SubMatrix(A, i, j, A.length);
				B[j][i] = Math.pow(-1, (i + j)) * get_Result(sub_A);
			}
		}
		return B;
	}

这里用到了get_SubMatrix（），该函数用于求解Aij对应的子矩阵，代码如下：

	/**获取当前矩阵某个元素的子矩阵
	 * @param A 当前矩阵
	 * @param m 第m行
	 * @param n 第n列
	 * @param length 矩阵的行
	 * @return 子矩阵
	 */
	public double[][] get_SubMatrix(double A[][], int m, int n, int length){
		double B[][] = new double[length - 1][length - 1];
		int i_b = 0;
		int j_b = 0;
		for(int i = 0; i < length; i++){
			//去掉第m行元素
			if(i == m)continue;
			for(int j = 0; j < length; j++){
				//去掉第n列元素
				if(j == n)continue;
				B[i_b][j_b] = A[i][j];
				j_b++;
			}
			i_b++;
			j_b = 0;
		}
		return B;
	}