numpy.linalg.pinv

numpy.linalg.pinv(a, rcond=1e-15) [source]

Вычислить псевдообратное (по Муру-Пенрозу) матрицы.

Вычислить обобщённое обращение матрицы, используя её сингулярное разложение (SVD) и включение всех больших сингулярных значений.

Параметры:	a : (M, N) array_like Матрица, для которой нужно вычислить псевдообращение. rcond : float Порог для малых сингулярных значений. Сингулярные значения, меньшие (по модулю) чем rcond * наибольшее_сингулярное_значение (снова, по модулю) устанавливаются в ноль.
Возвращает:	B : (N, M) ndarray Псевдообращение матрицы a. Если a является экземпляром matrix, то и B является экземпляром этого типа.
Возможные исключения:	LinAlgError Если вычисление SVD не сходится.

Примечания

Псевдообратное матрицы A, обозначаемое , определяется как: «матрица, которая „решает“ [задачу наименьших квадратов] », т.е., если является этим решением, то — это матрица такая, что .

Можно показать, что если — сингулярное разложение матрицы A, то , где — ортогональные матрицы, — диагональная матрица, состоящая из сингулярных значений A (обычно, после которых следуют нули), и тогда — просто диагональная матрица, состоящая из обратных сингулярных значений A (опять же, после которых следуют нули). [R47]

Ссылки

[R47]

(1, 2) G. Strang, Линейная алгебра и её приложения, 2-е изд., Орландо, FL, Academic Press, Inc., 1980, стр. 139-142.

Примеры

Следующий пример проверяет, что a * a+ * a == a и a+ * a * a+ == a+:

>>> a = np.random.randn(9, 6)
>>> B = np.linalg.pinv(a)
>>> np.allclose(a, np.dot(a, np.dot(B, a)))
True
>>> np.allclose(B, np.dot(B, np.dot(a, B)))
True