scipy.linalg.

pinv#

scipy.linalg.pinv(a, *, atol=None, rtol=None, return_rank=False, check_finite=True)[source]#

计算矩阵的（Moore-Penrose）伪逆。

使用矩阵的奇异值分解 U @ S @ V（在经济模式下）并仅选取与重要奇异值相关联的列/行来计算矩阵的广义逆。

如果 s 是 a 的最大奇异值，则重要性截止值由 atol + rtol * s 确定。任何低于此值的奇异值都被认为是不重要的。

参数：
a(M, N) array_like

要进行伪逆的矩阵。

atolfloat, 可选

绝对阈值项，默认值为 0。

在版本 1.7.0 中添加。

rtolfloat, 可选

相对阈值项，默认值为 max(M, N) * eps，其中 eps 是 a 数据类型机器精度值。

在版本 1.7.0 中添加。

return_rankbool, 可选

如果为 True，则返回矩阵的有效秩。

check_finitebool, 可选

是否检查输入矩阵是否仅包含有限数字。禁用可能带来性能提升，但如果输入确实包含无穷大或 NaN，则可能会导致问题（崩溃、不终止）。

返回值：
B(N, M) ndarray

矩阵 a 的伪逆。

rankint

矩阵的有效秩。如果 return_rank 为 True，则返回。

异常：
LinAlgError

如果 SVD 计算不收敛。

参见

pinvh

厄米矩阵的 Moore-Penrose 伪逆。

备注

如果 A 可逆，则 Moore-Penrose 伪逆恰好是 A 的逆 [1]。如果 A 不可逆，则 Moore-Penrose 伪逆计算 Ax = b 的 x 解，使得 ||Ax - b|| 最小化 [1]。

参考文献

[1] (1,2,3)

Penrose, R. (1956). On best approximate solutions of linear matrix equations. Mathematical Proceedings of the Cambridge Philosophical Society, 52(1), 17-19. doi:10.1017/S0305004100030929

示例

给定一个 m x n 矩阵 A 和一个 n x m 矩阵 B，四个 Moore-Penrose 条件是

  1. ABA = A (B 是 A 的广义逆)，

  2. BAB = B (A 是 B 的广义逆)，

  3. (AB)* = AB (AB 是厄米特)，

  4. (BA)* = BA (BA 是厄米特) [1]。

这里，A* 表示共轭转置。Moore-Penrose 伪逆是满足所有这四个条件的唯一 B，并且对于任何 A 都存在。请注意，与标准矩阵逆不同，A 不需要是方阵或具有线性无关的列/行。

例如，我们可以计算随机非方阵的 Moore-Penrose 伪逆，并验证它是否满足四个条件。

>>> import numpy as np
>>> from scipy import linalg
>>> rng = np.random.default_rng()
>>> A = rng.standard_normal((9, 6))
>>> B = linalg.pinv(A)
>>> np.allclose(A @ B @ A, A)  # Condition 1
True
>>> np.allclose(B @ A @ B, B)  # Condition 2
True
>>> np.allclose((A @ B).conj().T, A @ B)  # Condition 3
True
>>> np.allclose((B @ A).conj().T, B @ A)  # Condition 4
True
本页内容