numpy.linalg.inv

linalg.inv(a)

Вычислить обратную матрицу.

Для квадратной матрицы a, вернуть матрицу ainv удовлетворяющий a @ ainv = ainv @ a = eye(a.shape[0]).

Параметры:

a(…, M, M) array_like

Матрица для обращения.

Возвращает:

ainv(…, M, M) ndarray или matrix

Обратная матрица a.

Вызывает:

LinAlgError

Если a не является квадратной или инверсия не удалась.

Смотрите также

scipy.linalg.inv

Похожая функция в SciPy.

numpy.linalg.cond

Вычислить число обусловленности матрицы.

numpy.linalg.svd

Вычислить сингулярное разложение матрицы.

Примечания

Применяются правила трансляции, см. numpy.linalg документации для подробностей.

Если a обнаружена как вырожденная, выдается LinAlgError вызывается. Если a является плохо обусловленной, LinAlgError может быть или не быть вызвано, и результаты могут быть неточными из-за ошибок с плавающей запятой.

Ссылки

[1]

Википедия, «Число обусловленности», https://en.wikipedia.org/wiki/Condition_number

Примеры

Попробуйте в вашем браузере!

>>> import numpy as np
>>> from numpy.linalg import inv
>>> a = np.array([[1., 2.], [3., 4.]])
>>> ainv = inv(a)
>>> np.allclose(a @ ainv, np.eye(2))
True
>>> np.allclose(ainv @ a, np.eye(2))
True

Если a является объектом матрицы, то возвращаемое значение также является матрицей:

>>> ainv = inv(np.matrix(a))
>>> ainv
matrix([[-2. ,  1. ],
        [ 1.5, -0.5]])

Обратные матрицы нескольких матриц могут быть вычислены одновременно:

>>> a = np.array([[[1., 2.], [3., 4.]], [[1, 3], [3, 5]]])
>>> inv(a)
array([[[-2.  ,  1.  ],
        [ 1.5 , -0.5 ]],
       [[-1.25,  0.75],
        [ 0.75, -0.25]]])

Если матрица близка к вырожденной, вычисленная обратная матрица может не удовлетворять a @ ainv = ainv @ a = eye(a.shape[0]) даже если LinAlgError не вызывается:

>>> a = np.array([[2,4,6],[2,0,2],[6,8,14]])
>>> inv(a)  # No errors raised
array([[-1.12589991e+15, -5.62949953e+14,  5.62949953e+14],
   [-1.12589991e+15, -5.62949953e+14,  5.62949953e+14],
   [ 1.12589991e+15,  5.62949953e+14, -5.62949953e+14]])
>>> a @ inv(a)
array([[ 0.   , -0.5  ,  0.   ],  # may vary
       [-0.5  ,  0.625,  0.25 ],
       [ 0.   ,  0.   ,  1.   ]])

Для обнаружения плохо обусловленных матриц можно использовать numpy.linalg.cond для вычисления его число обусловленности [1]. Чем больше число обусловленности, тем более плохо обусловлена матрица. Как правило, если число обусловленности cond(a) = 10**k, тогда вы можете потерять до k цифр точности сверх того, что было бы потеряно из-за численного метода из-за потери точности от арифметических методов.

>>> from numpy.linalg import cond
>>> cond(a)
np.float64(8.659885634118668e+17)  # may vary

Также возможно обнаружить плохую обусловленность, проверяя сингулярные значения матрицы напрямую. Отношение между наибольшим и наименьшим сингулярными значениями является числом обусловленности:

>>> from numpy.linalg import svd
>>> sigma = svd(a, compute_uv=False)  # Do not compute singular vectors
>>> sigma.max()/sigma.min()
8.659885634118668e+17  # may vary

НазадОткрыть во вкладке