矩阵的LU分解
矩阵的LU分解
千面人 发表于1年前
矩阵的LU分解
  • 发表于 1年前
  • 阅读 41
  • 收藏 0
  • 点赞 0
  • 评论 0

腾讯云 新注册用户 域名抢购1元起>>>   

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

可以表示为...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

共有 人打赏支持
粉丝 13
博文 37
码字总数 20569
×
千面人
如果觉得我的文章对您有用,请随意打赏。您的支持将鼓励我继续创作!
* 金额(元)
¥1 ¥5 ¥10 ¥20 其他金额
打赏人
留言
* 支付类型
微信扫码支付
打赏金额:
已支付成功
打赏金额: