求逆矩阵

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

有了行列式,可以求解矩阵的逆矩阵,逆矩阵可以表示为 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;
	}

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值