文档章节

矩阵的LU分解

千面人
 千面人
发布于 2016/12/08 22:27
字数 322
阅读 58
收藏 0

矩阵的高斯消元法,通过行变换化成上三角阵。

可以表示为...E32E31E21A = U, 其中E表示行变换初等矩阵或消元矩阵,系数ij代表让A的ij位置变成0(暂不考虑行交换)

于是A = LU 其中L等于(...E32E31E21)的逆 编辑器不能用上下标写公式真是费劲 

而(...E32E31E21)取逆后,行变换互不干扰,计算L只需要把消元矩阵中的值取反放到相应位置上即可。(Eij的逆就是ij位置上的数取负号)

于是分解的计算过程伪代码如下:(只表达计算思路,不优化存储,返回结果等):

#从第一行到最后一行循环
for k from 1 to n:
    for i from k+1 to n:
        L[i][k] = A[i][k]/A[k][k]
        #对应的行相减
        row[i] -= L[i][k]*row[k]

#最后L下三角的数据(不包括对角线,对角线上都是1),A的上三角数据为U 为最终分解结果。

网上很多直接套公式,求单个位置上的数的算法,但必须要知道公式怎么来的,否则不好记忆。

公式如下图【1】:

 

【1】LU分解 - 百度百科

http://baike.baidu.com/link?url=FgR_JLQXZpkS1CPtC-ChvnRh_RtS7H-8Kktev5YEjS9TuEitEE6ZyivQVF8ntTkHx_CdkJpAmAj6xfNlmwPv_nDfKAPnpr5W2KN2YhIkb13

© 著作权归作者所有

共有 人打赏支持
千面人
粉丝 16
博文 48
码字总数 24904
作品 0
杭州
高级程序员
私信 提问
Armadillo之LU分解(LU factorisation/LU decomposition)

在armadillo库中,矩阵的LU分解(LU factorisation or LU decomposition)使用lu函数,lu函数有两个版本 1 lu(L,U,P,X) 其中X是欲进行分解的矩阵,分解生成L,U,P满足 1)P是一个置换矩阵(p...

桑梓狼狼
2014/08/01
0
0
基于javascript的矩阵LU分解的实现

在线性代数中,LU分解是将一个矩阵分解为 L(单位下三角矩阵)和 U(上三角矩阵),可用于求解线性方程组、反矩阵和计算行列式。本文结合LU分解,用javascript实现线性方程组的求解; 假设存...

qq_37338983
04/12
0
0
你知道MATLAB ,但你知道 NMATH吗?

NMath是一个.NET的数学库,包含了NET平台上的面向对象数字计算的基础类。 产品特点如下: 单精度和双精度复数类 为以下四种数据类型提供全功能的向量和矩阵类:单精度浮点数,双精度浮点数,...

愤怒的小吉
2014/10/24
360
0
Eigen 3.2.0-beta1 发布,线性算术的C++模板库

这个beta版本引入了内置的稀疏矩阵,和真正的QZ分解和广义特征求解稠密矩阵的LU和QR因子分解,以及Ref<>参考类。同时修复了一些bug。 Eigen 是一个线性算术的C++模板库,包括:vectors, matr...

zino
2013/03/08
894
6
SP++3.0 发布,欢迎大家使用

消息来自 Jerry 的博客: SP++ (Signal Processing in C++) 是一个关于信号处理与数值计算的开源C++程序库,该库提供了信号处理与数值计算中常用算法的C++实现。SP++中所有算法都以C++类模板...

红薯
2011/02/12
4.5K
4

没有更多内容

加载失败,请刷新页面

加载更多

Spring核心概念--Spring01

小生初出程序茅庐,走上编程之路,还请各位大佬多多管照。 初学Spring框架: 创建HelloWorldSpring项目 用eclipse开发Web项目,添加jar包: commons-logging-1.2.jar log4j-1.2.17.jar spring-...

小橙子的曼曼
21分钟前
1
0
MySQL 运行状态监控方法

一、通过shell脚本监控mysql的qps,tps,io详细见附件脚本《mysqlgather.sh》 因在脚本中直接设置密码会提示告警信息,需要在my.cnf文件中的[client]标签下增加默认的用户名和密码,并重启mys...

PeakFang-BOK
25分钟前
2
0
ROS实操笔记 四 topic (主题)

主题 消息以一种发布/订阅的方式传递。一个节点可以在一个给定的主题中发布消息。一个节点针对某个主题关注与订阅特定类型的数据。可能同时有多个节点发布或者订阅同一个主题的消息。总体上,...

placido
26分钟前
1
0
【NLP】【六】gensim之doc2vec

【一】总述 doc2vec是指将句子、段落或者文章使用向量来表示,这样可以方便的计算句子、文章、段落的相似度。 【二】使用方法介绍 1. 预料准备 def read_corpus(fname, tokens_only=False):...

muqiusangyang
29分钟前
0
0
node中process.nextTick & promise & 异步IO & setTimeout & setImmediate 的优先级

process.nextTick > promise > setTimeout > 异步IO > setImmediate

小草先森
35分钟前
1
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部