侯体宗的博客
  • 首页
  • Hyperf版
  • beego仿版
  • 人生(杂谈)
  • 技术
  • 关于我
  • 更多分类
    • 文件下载
    • 文字修仙
    • 中国象棋ai
    • 群聊
    • 九宫格抽奖
    • 拼图
    • 消消乐
    • 相册

矩阵求逆(JAVA)初等行变换

Java  /  管理员 发布于 5年前   588
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();
}
}
}

 



  • 上一条:
    深入浅出Java Annotation(元注解和自定义注解)
    下一条:
    矩阵求逆(JAVA)利用伴随矩阵
  • 昵称:

    邮箱:

    0条评论 (评论内容有缓存机制,请悉知!)
    最新最热
    • 分类目录
    • 人生(杂谈)
    • 技术
    • linux
    • Java
    • php
    • 框架(架构)
    • 前端
    • ThinkPHP
    • 数据库
    • 微信(小程序)
    • Laravel
    • Redis
    • Docker
    • Go
    • swoole
    • Windows
    • Python
    • 苹果(mac/ios)
    • 相关文章
    • 在java中实现的脱敏工具类代码示例分享(0个评论)
    • zookeeper安装流程步骤(0个评论)
    • 在java中你背的“八股文”可能已经过时了(2个评论)
    • 在php8.0+版本中使用属性来增加值代码示例(3个评论)
    • java 正则表达式基础,实例学习资料收集大全 原创(0个评论)
    • 近期文章
    • 在go语言中使用api.geonames.org接口实现根据国际邮政编码获取地址信息功能(1个评论)
    • 在go语言中使用github.com/signintech/gopdf实现生成pdf分页文件功能(0个评论)
    • gmail发邮件报错:534 5.7.9 Application-specific password required...解决方案(0个评论)
    • 欧盟关于强迫劳动的规定的官方举报渠道及官方举报网站(0个评论)
    • 在go语言中使用github.com/signintech/gopdf实现生成pdf文件功能(0个评论)
    • Laravel从Accel获得5700万美元A轮融资(0个评论)
    • 在go + gin中gorm实现指定搜索/区间搜索分页列表功能接口实例(0个评论)
    • 在go语言中实现IP/CIDR的ip和netmask互转及IP段形式互转及ip是否存在IP/CIDR(0个评论)
    • PHP 8.4 Alpha 1现已发布!(0个评论)
    • Laravel 11.15版本发布 - Eloquent Builder中添加的泛型(0个评论)
    • 近期评论
    • 122 在

      学历:一种延缓就业设计,生活需求下的权衡之选中评论 工作几年后,报名考研了,到现在还没认真学习备考,迷茫中。作为一名北漂互联网打工人..
    • 123 在

      Clash for Windows作者删库跑路了,github已404中评论 按理说只要你在国内,所有的流量进出都在监控范围内,不管你怎么隐藏也没用,想搞你分..
    • 原梓番博客 在

      在Laravel框架中使用模型Model分表最简单的方法中评论 好久好久都没看友情链接申请了,今天刚看,已经添加。..
    • 博主 在

      佛跳墙vpn软件不会用?上不了网?佛跳墙vpn常见问题以及解决办法中评论 @1111老铁这个不行了,可以看看近期评论的其他文章..
    • 1111 在

      佛跳墙vpn软件不会用?上不了网?佛跳墙vpn常见问题以及解决办法中评论 网站不能打开,博主百忙中能否发个APP下载链接,佛跳墙或极光..
    • 2016-11
    • 2018-03
    • 2020-03
    • 2023-05
    • 2023-11
    • 2024-01
    Top

    Copyright·© 2019 侯体宗版权所有· 粤ICP备20027696号 PHP交流群

    侯体宗的博客