2016年11月18日

矩阵的QR分解(三种方法)Python实现

摘要: 1.Gram-Schmidt正交化 假设原来的矩阵为[a,b],a,b为线性无关的二维向量,下面我们通过Gram-Schmidt正交化使得矩阵A为标准正交矩阵: 假设正交化后的矩阵为Q=[A,B],我们可以令A=a,那么我们的目的根据AB=I来求B,B可以表示为b向量与b向量在a上的投影的误差向量: 阅读全文

posted @ 2016-11-18 21:53 1357 阅读(19974) 评论(0) 推荐(0) 编辑

LU分解和求解线性方程组

摘要: 1 # coding:utf8 2 import numpy as np 3 4 def lu(mat): 5 r,c=np.shape(mat) 6 s=min(r,c) 7 for k in range(s): 8 x=1.0/mat[k][k] # 将后续除法变成乘法 9 for i in range(k+1,r... 阅读全文

posted @ 2016-11-18 17:36 1357 阅读(2303) 评论(0) 推荐(0) 编辑

导航