1. 最佳求解方案
Most efficient way to solve a system of linear equations
求解形如 Ax=b 的最佳方式
将 A 分解为三角矩阵,A=M1⋅M2
通过 M1⋅y=b 求解 y,使用 back substitution
通过 M2⋅x=y 求解 x,使用 ba...
基于吉文斯变换的QR分解
# 基于吉文斯变换的QR分解 import numpy as np def givens(A): m = A.shape[0] n = A.shape[1] P = np.mat(np.eye(n)) p1 = P for i in range(0, n-1): for j in range(i+1, n): if ...