杜利特尔分解(LU分解)的JAVA代码

2024-03-14 10:08
文章标签 java 代码 分解 lu 杜利特

本文主要是介绍杜利特尔分解(LU分解)的JAVA代码,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

杜利特尔分解(LU分解)的JAVA代码

数学基础:设A为n阶方阵,满足其各阶主子式:
det (Subscript[A, p]) != 0
A=LU

其代码如下:
public static double[][] LUDecomposition(double a[][],String X)//X代表返回矩阵的种类,若为U,返回u矩阵;若为l,返回l矩阵。若为其他,返回a;

    {int row=a.length;//row为矩阵行数int line=a[0].length;//line为矩阵列数double temp=0;if(row!=line) {System.out.println("矩阵行列数不等!!");return a;}//如果矩阵行列数不等,返回原矩阵。double[][] l=new double[row][line];//定义L矩阵double[][] u=new double[row][line];//定义U矩阵int k=0;for(int j=0;j<line;j++)//首先计算U矩阵第一行各元素u[k][j]=a[k][j];for(int i=0;i<line;i++)//计算L矩阵第一列各元素l[i][k]=a[i][k]/u[k][k];for(k=1;k<row;k++)//计算U矩阵第k行元素{for(int j=k;j<row;j++) {temp=0;for(int p=0;p<k;p++) {temp+=l[k][p]*u[p][j];}u[k][j]=a[k][j]-temp;}//求出第k行的ufor(int i=k;i<line;i++) //计算L矩阵第k列元素{temp=0;for(int p=0;p<k;p++) {temp+=l[i][p]*u[p][k];}l[i][k]=(a[i][k]-temp)/u[k][k];}//求出第k列的j}System.out.println("L矩阵为:");printMatrix(l,"0.00");//打印矩阵l,此为自定义函数System.out.println("U矩阵为:");printMatrix(u,"0.00");//打印矩阵l,此为自定义函数if(X=="U"||X=="u")return u;//返回U矩阵else if(X=="L"||X=="l")return l;//返回L矩阵else return a;//返回原矩阵}public static void printMatrix(double[][] a,String precision) {int row=a.length;int column=a[0].length;for(int i=0;i<row;i++) {for(int j=0;j<column;j++){System.out.print(new DecimalFormat(precision).format(a[i][j])+" ");}System.out.println();}System.out.println();
}

下面在主函数里调用这个函数:
public static void main(String[] args) {
double[][] a= {{6, 10, 7}, {1, 7, 8}, {10, 6, 3}};double[] b= {1,0,-3};
double[][] l=LUDecomposition(a,”l”);
double[][] u=LUDecomposition(a,”u”);
double[] y=SolveLy(l,b);
double[] x=SolveUx(u,y);
}

这篇关于杜利特尔分解(LU分解)的JAVA代码的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/808067

相关文章

SpringBoot集成EasyPoi实现Excel模板导出成PDF文件

《SpringBoot集成EasyPoi实现Excel模板导出成PDF文件》在日常工作中,我们经常需要将数据导出成Excel表格或PDF文件,本文将介绍如何在SpringBoot项目中集成EasyPo... 目录前言摘要简介源代码解析应用场景案例优缺点分析类代码方法介绍测试用例小结前言在日常工作中,我们经

Java实现预览与打印功能详解

《Java实现预览与打印功能详解》在Java中,打印功能主要依赖java.awt.print包,该包提供了与打印相关的一些关键类,比如PrinterJob和PageFormat,它们构成... 目录Java 打印系统概述打印预览与设置使用 PageFormat 和 PrinterJob 类设置页面格式与纸张

Spring Security常见问题及解决方案

《SpringSecurity常见问题及解决方案》SpringSecurity是Spring生态的安全框架,提供认证、授权及攻击防护,支持JWT、OAuth2集成,适用于保护Spring应用,需配置... 目录Spring Security 简介Spring Security 核心概念1. ​Securit

SpringBoot+EasyPOI轻松实现Excel和Word导出PDF

《SpringBoot+EasyPOI轻松实现Excel和Word导出PDF》在企业级开发中,将Excel和Word文档导出为PDF是常见需求,本文将结合​​EasyPOI和​​Aspose系列工具实... 目录一、环境准备与依赖配置1.1 方案选型1.2 依赖配置(商业库方案)二、Excel 导出 PDF

Python实现MQTT通信的示例代码

《Python实现MQTT通信的示例代码》本文主要介绍了Python实现MQTT通信的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 安装paho-mqtt库‌2. 搭建MQTT代理服务器(Broker)‌‌3. pytho

SpringBoot改造MCP服务器的详细说明(StreamableHTTP 类型)

《SpringBoot改造MCP服务器的详细说明(StreamableHTTP类型)》本文介绍了SpringBoot如何实现MCPStreamableHTTP服务器,并且使用CherryStudio... 目录SpringBoot改造MCP服务器(StreamableHTTP)1 项目说明2 使用说明2.1

spring中的@MapperScan注解属性解析

《spring中的@MapperScan注解属性解析》@MapperScan是Spring集成MyBatis时自动扫描Mapper接口的注解,简化配置并支持多数据源,通过属性控制扫描路径和过滤条件,利... 目录一、核心功能与作用二、注解属性解析三、底层实现原理四、使用场景与最佳实践五、注意事项与常见问题六

Spring的RedisTemplate的json反序列泛型丢失问题解决

《Spring的RedisTemplate的json反序列泛型丢失问题解决》本文主要介绍了SpringRedisTemplate中使用JSON序列化时泛型信息丢失的问题及其提出三种解决方案,可以根据性... 目录背景解决方案方案一方案二方案三总结背景在使用RedisTemplate操作redis时我们针对

Java中Arrays类和Collections类常用方法示例详解

《Java中Arrays类和Collections类常用方法示例详解》本文总结了Java中Arrays和Collections类的常用方法,涵盖数组填充、排序、搜索、复制、列表转换等操作,帮助开发者高... 目录Arrays.fill()相关用法Arrays.toString()Arrays.sort()A

Spring Boot Maven 插件如何构建可执行 JAR 的核心配置

《SpringBootMaven插件如何构建可执行JAR的核心配置》SpringBoot核心Maven插件,用于生成可执行JAR/WAR,内置服务器简化部署,支持热部署、多环境配置及依赖管理... 目录前言一、插件的核心功能与目标1.1 插件的定位1.2 插件的 Goals(目标)1.3 插件定位1.4 核