package gaodai.matrix; import gaodai.determinant.DeterminantCalculation; import java.util.ArrayList; import java.util.List; import java.util.Scanner; /** * 矩阵求逆(初等行变换) * @author 邱万迟 * */ public class InverseOfMatrix { private List<List<Double>> matrix; private int lineNum; private int columnNum; public InverseOfMatrix(List<List<Double>> data){ matrix = data; lineNum = data.size(); columnNum = data.get(0).size(); } public void caculate() throws Exception{ //1.非方正不能求逆 //2.方正的行列式值为零不能求逆 if( lineNum != columnNum){ throw new Exception("此矩阵不能求逆>>>>>>>>>>>>>>>>>"); } List<List<Double>> tempList = new ArrayList<List<Double>>(); for(List<Double> l : matrix){ List<Double> newList = new ArrayList<Double>(); newList.addAll(l); tempList.add(newList); } DeterminantCalculation d = new DeterminantCalculation(tempList); d.chang2UpperTriangle(); double result = d.getValue(); if(result == 0){ throw new Exception("此矩阵不能求逆>>>>>>>>>>>>>>>>>"); } //增加单位矩阵 for(int i = 0; i < lineNum; i++){ List<Double> list = matrix.get(i); for(int j = 0; j < columnNum; j++){ if(i == j){ list.add(1.0); }else{ list.add(0.0); } } } print(); chang2UpperTriangle();//化为上三角 changeReducedMatrix();//化为约化矩阵 print(); } public void getValue(){ boolean flag = true; for(int i = 0; i < lineNum; i++){ if(matrix.get(i).get(i) == 0){ flag = false; } if(!flag){ break; } } if(!flag){ System.out.println("此矩阵不可逆>>>>>>>>>>>>>>"); }else{ for(int i = 0; i < lineNum; i++){ List<Double> list = matrix.get(i); for(int j = 0; j < columnNum; j++){ list.remove(0); } } System.out.println("逆矩阵为>>>>>>>>>>>>>>>>>"); print(); } } /** * 打印 */ public void print() { int i = 0, j = 0; for (List<Double> line : matrix) { for (double element : line) { System.out.print(element); System.out.print("(" + i + "," + j + ") "); System.out.print(" "); j++; } System.out.println(); i++; j = 0; } System.out.println(); } /** * 校验是否是上三角,不是就的继续计算 * * @return */ public boolean isCaculate() { boolean hasCaculate = false; for (int i = 0; i < matrix.size(); i++) { for (int j = 0; j < i; j++) { if (matrix.get(i).get(j) != 0.0) { System.out.println("(" + (i + 1) + "," + (j + 1) + ")元素值不为零"); hasCaculate = true; break; } } if (hasCaculate) { break; } } return hasCaculate; } private int caculateTimes; /** * 化为上三角 * @throws Exception */ public void chang2UpperTriangle() throws Exception { if (!isCaculate()) { return; } int min = lineNum; caculateTimes++; System.out.println("--------------第" + caculateTimes + "次计算--------------"); for (int i = 0; i < min; i++) { for (int j = i + 1; j < min; j++) { double multiplyNum = -1 * matrix.get(j).get(i) / matrix.get(i).get(i); if (multiplyNum == 0) { continue; } this.lineMultiplyNumAdd2OtherLine(multiplyNum, (i + 1), (j + 1)); print(); } } print(); chang2UpperTriangle(); } /** * 变为约化矩阵 */ public void changeReducedMatrix() throws Exception{ for(int i = 0; i < lineNum; i++){//行 if(i == 0){ //continue; } List<Double> temp = matrix.get(i); for(Double d : temp){ if(d == 0){ continue; } double multiplyNum = 1.0 / d; for(int a = 0; a < temp.size(); a++){ temp.set(a, temp.get(a) * multiplyNum); } break; } print(); for(int j = 0; j <= columnNum; j++){//列 if(temp.get(j) != 0){//这个数不为零 ,此数为第 i行第j列 for(int t = 0; t < lineNum; t++){//行 if(t == i || matrix.get(t).get(j) == 0){//本列的其他行 continue; } double multiplyNum = -1 * matrix.get(t).get(j) / temp.get(j); this.lineMultiplyNumAdd2OtherLine(multiplyNum, (i + 1), (t + 1)); print(); } break; } } } } /** * 第a行乘以number 加到第b行上 * @param number 乘以的数 * @param a行号 * @param b行号 * @throws Exception */ public void lineMultiplyNumAdd2OtherLine(double number, int a, int b) throws Exception { if (a < 1 || a > matrix.size() || b < 1 || b > matrix.size()) { throw new Exception("输入的行号不合法"); } List<Double> aLine = matrix.get(a - 1); List<Double> bLine = matrix.get(b - 1); for (int i = 0; i < bLine.size(); i++) { double temp = bLine.get(i) + aLine.get(i) * number; bLine.set(i, temp); } System.out.println("第" + a + "行乘以" + number + "加到第" + b + "行:"); } public static void main(String[] args) { Scanner scanner = new Scanner(System.in); System.out.println("请输入矩阵的行数与列数,用逗号分隔:"); String sn = scanner.next(); String[] snArr = sn.split(","); int lineNum = Integer.valueOf(snArr[0]); int columnNum = Integer.valueOf(snArr[1]); List<List<Double>> matrix = new ArrayList<List<Double>>(); for(int i = 0; i < lineNum; i++){ System.out.println("请输入第" + (i + 1) + "行的数,用逗号分隔:"); String lineData = scanner.next(); String[] lineDataArr = lineData.split(","); List<Double> line = new ArrayList<Double>(); matrix.add(line); for(int j = 0; j < columnNum; j++){ line.add(Double.valueOf(lineDataArr[j])); } } InverseOfMatrix m = new InverseOfMatrix(matrix); m.print(); try { m.caculate(); m.getValue(); } catch (Exception e) { e.printStackTrace(); } } }
相关推荐
求逆矩阵的方法与矩阵的秩 一、矩阵的初等行变换 (由定理2.4给出的求逆矩阵的伴随矩阵法,要求计算矩阵A的行列式 ...下面介绍利用矩阵初等行变换求逆矩阵的方法.在介绍这种方法之前,先给出矩阵初等行变换的定义.)
复数矩阵求逆,用了STL中的vector表示矩阵,用了complex类表示复数,求逆矩阵使用了初等变换求逆矩阵的方法
逆矩阵的求法行列式初等变换 在前面我们以学习了用公式求逆矩阵,但当矩阵A的阶数较大时,求A*很繁琐,此方法不实用,因此必须找一种更简单的方法求逆矩阵,那么如何找到一种简单的方法呢?
这是矩阵运算的代码希望对大家有用,这个代码简单易懂,希望初学者可以下载后查看,相关解释我也比较详细
本人使用C语言编写使用初等行变换的方法,求出矩阵的逆矩阵。
一个计算矩阵求逆、矩阵转置和矩阵相乘的类。矩阵求逆包括初等变换法和伴随矩阵法。
复数矩阵求逆
下面是通过初等行变换求解逆矩阵的,非用伴随矩阵! 比如求: 1 2 3 2 2 1 3 4 3 则把下面矩阵: 1 2 3 1 0 0 2 2 1 0 1 0 3 4 3 0 0 1 通过初等行变换变为: 1 0 0 x x x 0 1 0 x x x 0 0 1 x x x 则x矩阵就是要求的...
该资源是对于多入多出矩阵求逆的C++编程语言实现,用于多入多出的矩阵求逆算法,本人是用机器人系统设计控制算法,亲自验证代码的可行性
求N阶矩阵的逆矩阵,由java编写,经测试完全正确,很好用。
在进行矩阵求逆等计算时遇到矩阵条件数较大导致病态时,常用的各种解决病态方程方法
分块矩阵的初等变换
实数矩阵求逆,用了STL中的vector表示矩阵,求逆矩阵使用了初等变换求逆矩阵的方法
矩阵算法中基础的初等行变换代码,详细有注释,易于理解
利用矩阵lu分解的优秀特性进行矩阵求逆的代码,减少求逆计算量,大约200行 求逆矩阵思路: 1.求矩阵的crout(LU)分解,其中L为下三角矩阵,U为上三角矩阵 2.求L,U矩阵的伴随阵,参考文献:三角形矩阵求伴随矩阵...
矩阵求逆在编程过程中会经常用到,所以一个简单而实用的矩阵求逆对于顺利的编好程序是十分有用处的
C++写的矩阵求逆的代码,输入维数和矩阵,利用矩阵初等变换求逆矩阵,代码比较精简。
矩阵求逆矩阵求逆矩阵求逆矩阵求逆矩阵求逆矩阵求逆
Gauss-Jordan法实矩阵求逆 求4阶矩阵的逆矩阵 算法很实用啊
用Mathematica处理矩阵是很方便的。 而用矩阵描述几何变换,又是非常给力的。