当前位置:首页 > 开发 > 行业应用 > 正文

矩阵求逆(JAVA)利用伴随矩阵

发表于: 2015-06-24   作者:qiuwanchi   来源:转载   浏览:
摘要: package gaodai.matrix; import gaodai.determinant.DeterminantCalculation; import java.util.ArrayList; import java.util.List; import java.util.Scanner; /** * 矩阵求逆(利用伴随矩阵) * @author 邱万迟
package gaodai.matrix;

import gaodai.determinant.DeterminantCalculation;

import java.util.ArrayList;
import java.util.List;
import java.util.Scanner;

/**
 * 矩阵求逆(利用伴随矩阵)
 * @author 邱万迟
 *
 */
public class InverseOfMatrix2 {
	
	private List<List<Double>> matrix;
	private int lineNum;
	private int columnNum;
	private double determinantValue;
	private List<List<Double>> adjointMatrix;

	public List<List<Double>> getMatrix() {
		return matrix;
	}

	public double getDeterminantValue() {
		return determinantValue;
	}

	public void setMatrix(List<List<Double>> matrix) {
		this.matrix = matrix;
	}

	public void setDeterminantValue(double determinantValue) {
		this.determinantValue = determinantValue;
	}
	
	public InverseOfMatrix2(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("此矩阵不能求逆>>>>>>>>>>>>>>>>>");
		}
		
		determinantValue = result;
		
		adjointMatrix = new ArrayList<List<Double>>();//伴随矩阵
		
		for(int i = 0; i < lineNum; i++){
			List<Double> line = new ArrayList<Double>();
			adjointMatrix.add(line);
			for(int j = 0; j < columnNum; j++){
				
				List<List<Double>> list = new ArrayList<List<Double>>();//余子式
				for(int t = 0; t < lineNum; t++){
					if(i == t){
						continue;
					}
					List<Double> newList = new ArrayList<Double>();
					list.add(newList);
					for(int k = 0; k < columnNum; k++){
						if(j == k){
							continue;
						}
						newList.add(matrix.get(t).get(k));
					}
				}
				
				DeterminantCalculation determinant = new DeterminantCalculation(list);
				determinant.chang2UpperTriangle();
				double tempValue = determinant.getValue();
				int sign = getSign(i, j);
				System.out.println( "(" + i + j + ")的余子式:" + tempValue);
				line.add(sign * tempValue);
			}
		}
		System.out.println("伴随矩阵的转置>>>>>>>>>>>>>>>>>>>>");
		print(adjointMatrix);
		int lineNum = adjointMatrix.size();
		int columnNum = adjointMatrix.get(0).size();
		for(int i = 0; i < lineNum; i++){
			for(int j = i; j < columnNum; j++){
				double t = adjointMatrix.get(i).get(j);
				adjointMatrix.get(i).set(j,adjointMatrix.get(j).get(i));
				adjointMatrix.get(j).set(i,t);
			}
		}
		
		System.out.println("伴随矩阵如下>>>>>>>>>>>>>>>>>>>>");
		print(adjointMatrix);
		System.out.println("逆矩阵如下>>>>>>>>>>>>>>>>>>>>");
		for(int i = 0; i < lineNum; i++){
			for(int j = i; j < columnNum; j++){
				adjointMatrix.get(i).set(j,adjointMatrix.get(i).get(j) / determinantValue);
			}
		}
		print(adjointMatrix);
	}
	
	private int getSign(int i,int j){
		if((i + j) % 2 == 0){
			return 1;
		}
		return -1;
	}
	
	/**
	 * 打印
	 */
	public void print(List<List<Double>> data) {
		int i = 0, j = 0;
		for (List<Double> line : data) {
			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();
	}
	
	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]));
			}
		}
		
		InverseOfMatrix2 m = new InverseOfMatrix2(matrix);
		m.print(m.getMatrix());
		try {
			m.caculate();
		} catch (Exception e) {
			e.printStackTrace();
		}
	}
}

 

 

矩阵求逆(JAVA)利用伴随矩阵

  • 0

    开心

    开心

  • 0

    板砖

    板砖

  • 0

    感动

    感动

  • 0

    有用

    有用

  • 0

    疑问

    疑问

  • 0

    难过

    难过

  • 0

    无聊

    无聊

  • 0

    震惊

    震惊

编辑推荐
设R是一个交换环,A是一个以R中元素为系数的 n×n 的矩阵。A的伴随矩阵可按如下步骤定义: 定义:A关
设R是一个交换环,A是一个以R中元素为系数的 n×n 的矩阵。A的伴随矩阵可按如下步骤定义: 定义:A关
作为一只数学基础一般般的程序猿,有时候连怎么求逆矩阵都不记得,之前在wikiHow上看了一篇不错的讲
一、矩阵乘法的五种表示方法 1、一般形式 2、矩阵与列向量相乘 <span style="font-family:'Micro
高斯消元法可以用来找出一个可逆矩阵的逆矩阵。设A 为一个N * N的矩阵,其逆矩阵可被两个分块矩阵表
为了实现Hill密码,我们必须先实现如何求逆矩阵。 A的逆矩阵公式为 A*/|A|,我们需要分别求出A*、|A
本节介绍伴随矩阵的概念,伴随矩阵可用于解超定方程 1 伴随 设X和U是两个欧几里得空间,若A是X到U的
求一个矩阵的逆(inverse or multiplicative inverse)使用矩阵的.i()方法或用inv()函数 m.i() //返
高斯约旦消去法是求矩阵逆的一种常用算法,但使用计算机来求解时,需要开辟另一个内存存放变化时的
1.代码思路 1)对矩阵进行合法性检查:矩阵必须为方阵 2)计算矩阵行列式的值(Determinant函数) 3
版权所有 IT知识库 CopyRight © 2009-2015 IT知识库 IT610.com , All Rights Reserved. 京ICP备09083238号