np.linalg.qr(A) 或 scipy.linalg.orth(A) 用于寻找正交基

Question

如果我有一个由五个向量 v1....v5 跨越的向量空间，找到 A 的正交基，其中 A=[v1,v2...v5] 并且 A 是 5Xn我应该使用np.linalg.qr(A) 还是 scipy.linalg.orth(A)？？提前致谢

Answer 1

你会想要使用：

scipy.linalg.orth(A)

普遍接受的规则是使用 scipy.linalg - 因为它比 np.linalg 涵盖更多的功能。希望有帮助！

Answer 2

请注意，sp.linalg.orth使用 SVD 时np.linalg.qr使用 QR 分解。两种分解都是通过 LAPACK 函数的包装器获得的。

我不认为对一个比另一个有强烈的偏好。SVD 会稍微稳定一些，但计算速度也会慢一些。在实践中，我认为您不会真正看到太大的不同。