北京pk10最新历史记录:详解用java描述矩阵求逆的算法

北京快三走势图 www.mxld9.cn  更新时间:2019年03月12日 15:48:02   作者:Thomas_Wang   我要评论

这篇文章主要介绍了用java描述矩阵求逆的算法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

今天很开心把困扰几天的问题解决了,在学习线性代数这门课程的时候。想通过程序实现里面的计算方法,比如矩阵求逆,用java代码该如何描述呢?

首先,咱们先用我们所交流语言描述一下算法思路:

1.求出一个矩阵A对应的行列式在第i,j(i表示行,j表示列)位置的余子式(余子式前面乘以-1^(i+j)即得代数余子式);

2.根据代数余子式求得矩阵A行列式的值。(行列式展开法);

3.根据代数余子式和行列式的值求出伴随矩阵;

4.由伴随矩阵和矩阵行列式值求逆矩阵。(A^-1 = A* / |A|)。

了解上述算法思路后,废话少说,上代码。

1.求出一个矩阵A对应的行列式在第i,j(i表示行,j表示列)位置的余子式(余子式前面乘以-1^(i+j)即得代数余子式);

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

A行列式的值。(行列式展开法);

/**
  * 求矩阵的行列式的值
  * @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;
 }

3.根据代数余子式和行列式的值求出伴随矩阵;

/**
     * 求矩阵的伴随矩阵
     * @param mat
     * @return
     */
    public static Matrix getWithMatrix(Matrix mat) {
        //创建一个矩阵存放伴随矩阵的值
        Matrix withMatrix = new Matrix(mat.getRow(),mat.getCol());
        //遍历withMatrix存放对应的mat的值
        for (int i = 0; i < withMatrix.getRow(); i++) {
            for (int j = 0; j < withMatrix.getCol(); j++) {
                double temp = Math.pow(-1, i+j) * MatrixUtil.getMatrixValue(MatrixUtil.getComplementMinor(mat, j, i));
                if(Math.abs(temp) <= 10e-6) temp = 0;
                withMatrix.setValue(i, j,temp);
            }
        }
        //返回结果
        return withMatrix;   
    }

4.由伴随矩阵和矩阵行列式值求逆矩阵。(A^-1 = A* / |A|)。

/**
  * 求逆矩阵
  * @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;
  
 }

以上所述是小编给大家介绍的用java描述矩阵求逆的算法详解整合,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对脚本之家网站的支持!

相关文章

最新评论

  • 星战外传再起——关于《游侠索罗》你应该知道的N件事 2019-05-26
  • “威马逊”在海南翁田镇沿海登陆 最大风力达17级 2019-05-25
  • 新版外商投资负面清单即将发布 2019-05-24
  • 解说米勒谈WE:LPL夏季赛夺冠可能性不大 2019-05-23
  • 陈坤:虽然会被困扰,但感谢情绪无法删除 2019-05-22
  • 问题来了!那些年你遇到的噪声扰民怎么办? 2019-05-21
  • 俞正声:广泛开展对外友好交往 2019-05-21
  • 运用互联网思维介入城市化会如何? 2019-05-20
  • 【加拿大房产网加拿大新房加拿大房产信息网】 2019-05-20
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-19
  • [网连中国]5地社保费率降低后工资条看变化 职工:只少交几块钱 2019-05-19
  • 全国政协委员侯一筠:发挥山东海洋优势,打造国家级科技成果转化平台 2019-05-18
  • 广西河池消防支队召开新兵第二阶段集训工作部署会 2019-05-18
  • 龙峰:帮助更多企业用好互联网—上游新闻对话重庆经济 2019-05-17
  • 其实,生产力发展了,社会财富丰富了,把小萌们养起来也不是什么问题……但你们不能被养着还养出脾气来还妄图对真正的劳动者指手画脚! 2019-05-16
  • 545| 209| 111| 553| 265| 916| 620| 207| 588| 617|