LU分解

deepfuture

浏览: 4411093 次
性别:
来自: 湛江

最近访客更多访客>>

linxl2011

mars36

jccz_zys

zkm0309

博主相关

博客

微博

相册

留言

关于我

博客专栏

: SQLite源码剖析
浏览量：80123

: WIN32汇编语言学习应用...
浏览量：70323

: 神奇的perl
浏览量：103573

: lucene等搜索引擎解析...
浏览量：286520

: 深入lucene3.5源码...
浏览量：15039

: VB.NET并行与分布式编...
浏览量：67776

: silverlight 5...
浏览量：32280

: 算法下午茶系列
浏览量：46069

文章分类

社区版块

存档分类

博客分类：

数学与计算

设A是一个方块矩阵。A的LU分解是将它分解成如下形式：

$A = LU, \,$

其中L和U分别是下三角矩阵和上三角矩阵。

例如对于一个 $3 \times 3$ 的矩阵，就有

$\begin{bmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \\ \end{bmatrix} = \begin{bmatrix} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \\ \end{bmatrix} \begin{bmatrix} u_{11} & u_{12} & u_{13} \\ 0 & u_{22} & u_{23} \\ 0 & 0 & u_{33} \\ \end{bmatrix}$ 。

一个LDU分解是一个如下形式的分解：

$A = LDU$

其中D是对角矩阵，L和U是单位三角矩阵（对角线上全是1的三角矩阵）。

一个LUP分解是一个如下形式的分解：

$A = LUP$

其中L和U仍是三角矩阵，P是一个置换矩阵。

一个充分消元的LU分解为如下形式：

$PAQ = LU$

将以下矩阵进行LU分解：

$\begin{bmatrix} 4 & 3 \\ 6 & 3 \\ \end{bmatrix} = \begin{bmatrix} l_{11} & 0 \\ l_{21} & l_{22} \\ \end{bmatrix} \begin{bmatrix} u_{11} & u_{12} \\ 0 & u_{22} \\ \end{bmatrix}.$

由于矩阵阶数只是2，可以直接列方程解：

$l_{11}* u_{11} + 0 * 0 = 4$

$l_{11}* u_{12} + 0 * u_{22} = 3$

$l_{21}* u_{11} + l_{22} * 0 = 6$

$l_{21}* u_{12} + l_{22} * u_{22} = 3.$

这个线性方程组有无数多组解。因此，可以假设其中一个是单位三角矩阵，比如说L，也就是说其对角线上的两个系数都是1。这时可以解出：

$l_{21} = 1.5$

$u_{11} = 4$

$u_{12} = 3$

$u_{22} = -1.5$ 。

也就是说

$\begin{bmatrix} 4 & 3 \\ 6 & 3 \\ \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ 1.5 & 1 \\ \end{bmatrix} \begin{bmatrix} 4 & 3 \\ 0 & -1.5 \\ \end{bmatrix}.$

分享到：

Forth简介　 | matlab-线性网络训练过程

2012-08-08 20:31
浏览 1504
评论(0)
分类:编程语言
查看更多

发表评论

您还没有登录,请您登录后再发表评论

最近访客更多访客>>

博主相关

博客专栏

文章分类

社区版块

存档分类

最新评论

LU分解

评论

发表评论

相关推荐

最近访客 更多访客>>

博主相关

博客专栏

文章分类

社区版块

存档分类

最新评论

LU分解

评论

发表评论

相关推荐

R语言与数据分析

机器学习实践指南：案例应用解析

matlab-矩阵合并

人工智能与数据分析所需要的知识

麦哈普的AI乐园【myhaspl@qq.com】我的另一个博客（机器学习、数据分析、智能计算的原创）

R-并行计算

谱聚类

对变化建模-用差分方程-动力系统及常数解

逻辑斯蒂映射-伪随机数

matlab-多项式乘除法及式子和导数

matlab-数组-元胞数据与结构数组

矩阵-范数

向量-范数

矩阵-求逆

lisp-猜数字算法与全局函数、变量

开源 Lisp 相关项目

四分位数求法

matlab-神经网络-自定义多层感知器解决异或(2)

matlab-神经网络-自定义多层感知器解决异或(1)

matlab-模态对话框

最近访客更多访客>>