您的位置:首页 > 编程语言

MIT线性代数课程笔记对应代码-【lecture 4】

2016-03-03 21:31 405 查看
矩阵的LU分解

from scipy.linalg import lu
import numpy as np

a = np.array([[2,1], [8,7]])
p, l, u = lu(a)


因为 PA = LU, 所以LU分解返回了三个参数,为了验证是否正确

np.mat(l) * np.mat(u)


发现就是a。(要转换为矩阵再操作,因为原始是array,array相乘是element的)
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: