[Java] 详解用java描述矩阵求逆的算法

/**<br>   * 求矩阵在i,j处余子式<br>   * @param mat<br>   * @param i<br>   * @param j<br>   * @return<br>   */<br>  public static matrix getcomplementminor(matrix mat, int i, int j) {<br>    //创建一个新的矩阵用于接收表示该余子式，需删除本行本列的数值<br>    matrix m = new matrix(mat.getrow()-1,mat.getcol()-1); <br>    //用于遍历新矩阵m的变量<br>    int row =0 ,col=0;<br>    /*<br>     * 遍历原矩阵的数据，j2表示行,k表示列<br>     */<br>    for (int j2 = 0; j2 < mat.getrow(); j2++) {<br>      //在第i行除的数据省略<br>      if(j2 == i) continue; <br>      for (int k = 0; k < mat.getcol(); k++) {<br>        //在第j列的数据省略<br>        if(k == j) continue;<br>        //赋值<br>        m.setvalue(row, col,mat.getvalue(j2, k));<br>        //遍历新矩阵的变量<br>        col++;<br>        if(col >= m.getcol() ) {<br>          col = 0;<br>          row++;<br>        }<br>      }<br>    }<br>    return m;<br>  }

/**
  * 求矩阵的行列式的值
  * @param mat
  * @return
  */
 public static double getmatrixvalue(matrix mat) {
  if(mat.getrow() != mat.getcol()) {
   system.out.println("该矩阵不是方阵，没有行列式");
   return double.min_value;
  }
  //若为1*1矩阵则直接返回
  if(mat.getrow() == 1) return mat.getvalue(0, 0); 
  //若为2*2矩阵则直接计算返回结果
  if(mat.getrow() == 2) {    
   return mat.getvalue(0, 0)*mat.getvalue(1, 1) - mat.getvalue(0, 1)*mat.getvalue(1, 0);
  }
  //行列式的值
  double matrixvalue = 0; 
  for (int i = 0; i < mat.getcol(); i++) {
   //获取0，i位置的余子式，即第一行的余子式
   matrix m = getcomplementminor(mat, 0, i);
   //将第一行的余子式相加 ，递归下去
   matrixvalue += math.pow(-1, i) * getmatrixvalue(m);
   
  }
  return matrixvalue;
 }

/**<br>   * 求矩阵的伴随矩阵<br>   * @param mat<br>   * @return<br>   */<br>  public static matrix getwithmatrix(matrix mat) {<br>    //创建一个矩阵存放伴随矩阵的值<br>    matrix withmatrix = new matrix(mat.getrow(),mat.getcol());<br>    //遍历withmatrix存放对应的mat的值<br>    for (int i = 0; i < withmatrix.getrow(); i++) {<br>      for (int j = 0; j < withmatrix.getcol(); j++) {<br>        double temp = math.pow(-1, i+j) * matrixutil.getmatrixvalue(matrixutil.getcomplementminor(mat, j, i));<br>        if(math.abs(temp) <= 10e-6) temp = 0;<br>        withmatrix.setvalue(i, j,temp);<br>      }<br>    }<br>    //返回结果<br>    return withmatrix;  <br>  }

/**
  * 求逆矩阵
  * @param mat
  * @return
  */
 public static matrix getrematrix(matrix mat) {
  //创建一个矩阵接收逆矩阵数据
  matrix rematrix = new matrix(mat.getrow(),mat.getcol()); 
  //得到原矩阵行列式的值
  double value = matrixutil.getmatrixvalue(mat); 
  //判断矩阵行列式的值是否为零
  if(math.abs(value) <= 10e-6) {
   system.out.println("该矩阵不可逆！");
   return null;
  }
  //将原矩阵mat赋值除以原行列式的值value给逆矩阵
  for (int i = 0; i < rematrix.getrow(); i++) {
   for (int j = 0; j < rematrix.getcol(); j++) {
  rematrix.setvalue(i, j, matrixutil.getwithmatrix(mat).getvalue(i, j) / value);
   }
  }
  return rematrix;
  
 }