矩阵计算-线性系统和 LU 分解

2024-03-23 13:04

本文主要是介绍矩阵计算-线性系统和 LU 分解,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、三角系统

……

二、高斯消元法

……

三、LU分解--直接三角分解法 求解线性方程Ax=b

参考视频:【数值分析】矩阵LU三角分解| 速成讲解 考试宝典_哔哩哔哩_bilibili

令A=L·U,其中L是单位下三角矩阵(对角线上元素都是1),U是上三角矩阵

无论矩阵A是多少,都有唯一的单位下三角矩阵L和上三角矩阵U,使得A可以分解为L·U。

由矩阵乘法可推出:

1.矩阵U的第一行元素等于矩阵A的第一行元素,即a1i=u1i,i=1,2……,n

2.矩阵L的第一列元素等于矩阵A的第一列元素除以u11,即li1=ai1/u11,i=1,2……,n

L是单位下三角矩阵 对角线上元素都是1,即lii=1 ,i=1,2……,n 

U是上三角矩阵。所以矩阵中的未知元素其实就不多了。

由于Ax=b即L·Ux=b

令Ly=b,Ux=y

第一步:首先根据A=L·U,可以很快推导出唯一的单位下三角矩阵L和唯一的上三角矩阵U

第二步:已知单位下三角矩阵L和向量b,根据Ly=b,很方便就求出y=(y1,y2,y3)^T

第三步:已知上三角矩阵U和向量y,根据Ux=y,很方便就求出x=(x1,x2,x3)

至此,线性方程求解完毕。

动手实操一下,就能清晰掌握了:

 第一步:根据A=L·U,推出矩阵L和U

第二步:根据Ly=b 求出y

第三步:根据Ux=求出x

向量x求解完毕。

四、部分选主元的 LU 分解

……

五、LDM 分解, 对称矩阵的 LDL 分解,正定矩阵的 Choleskey 分解

……

六、其他线性系统,如带状系统

……

这篇关于矩阵计算-线性系统和 LU 分解的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/838390

相关文章

Python如何计算两个不同类型列表的相似度

《Python如何计算两个不同类型列表的相似度》在编程中,经常需要比较两个列表的相似度,尤其是当这两个列表包含不同类型的元素时,下面小编就来讲讲如何使用Python计算两个不同类型列表的相似度吧... 目录摘要引言数字类型相似度欧几里得距离曼哈顿距离字符串类型相似度Levenshtein距离Jaccard相

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显