Линейная алгебра (numpy.linalg)

Функции линейной алгебры NumPy опираются на BLAS и LAPACK для обеспечения эффективных низкоуровневых реализаций стандартных алгоритмов линейной алгебры. Эти библиотеки могут предоставляться самим NumPy с использованием C-версий подмножества их эталонных реализаций, но, когда это возможно, отдается предпочтение высокооптимизированным библиотекам, которые используют возможности специализированных процессоров. Примерами таких библиотек являются OpenBLAS, MKL (TM) и ATLAS. Поскольку эти библиотеки многопоточны и зависят от процессора, могут потребоваться переменные среды и внешние пакеты, такие как threadpoolctl, для управления количеством потоков или указания архитектуры процессора.

Библиотека SciPy также содержит подмодуль linalg, и существует перекрытие в функциональности, предоставляемой подмодулями SciPy и NumPy. SciPy содержит функции, отсутствующие в numpy.linalg, такие как функции, связанные с LU-разложением и разложением Шура, несколько способов вычисления псевдообратной матрицы и матричные трансцендентные функции, такие как матричный логарифм. Некоторые функции, которые существуют в обоих, имеют расширенную функциональность в scipy.linalg. Например, scipy.linalg.eig может принимать второй матричный аргумент для решения обобщенных задач нахождения собственных значений. Однако некоторые функции NumPy имеют более гибкие возможности вещания. Например, numpy.linalg.solve может обрабатывать «упакованные» массивы, в то время как scipy.linalg.solve принимает только один квадратный массив в качестве своего первого аргумента.

Матричные и векторные произведения

`dot`(a, b[, out])	Скалярное произведение двух массивов.
`linalg.multi_dot`(arrays, \*[, out])	Вычисление скалярного произведения двух или более массивов в одном вызове функции, автоматически выбирая наиболее быстрый порядок вычислений.
`vdot`(a, b)	Возвращает скалярное произведение двух векторов.
`inner`(a, b)	Внутреннее произведение двух массивов.
`outer`(a, b[, out])	Вычисляет внешнее произведение двух векторов.
`matmul`(x1, x2, /[, out, casting, order, …])	Матричное произведение двух массивов.
`tensordot`(a, b[, axes])	Вычисление тензорного скалярного произведения по указанным осям.
`einsum`(subscripts, *operands[, out, dtype, …])	Вычисляет сумму по индексам Эйнштейна на операндах.
`einsum_path`(subscripts, *operands[, optimize])	Оценивает порядок сокращения с наименьшей стоимостью для выражения einsum, рассматривая создание промежуточных массивов.
`linalg.matrix_power`(a, n)	Возводит квадратную матрицу в (целую) степень `n`.
`kron`(a, b)	Произведение Кронекера двух массивов.

Разложения

`linalg.cholesky`(a)	Разложение Холецкого.
`linalg.qr`(a[, mode])	Вычисляет QR-разложение матрицы.
`linalg.svd`(a[, full_matrices, compute_uv, …])	Сингулярное разложение.

Собственные значения матриц

`linalg.eig`(a)	Вычисляет собственные значения и правые собственные векторы квадратного массива.
`linalg.eigh`(a[, UPLO])	Возвращает собственные значения и собственные векторы комплексной эрмитовой (сопряженно-симметричной) или вещественной симметричной матрицы.
`linalg.eigvals`(a)	Вычисляет собственные значения произвольной матрицы.
`linalg.eigvalsh`(a[, UPLO])	Вычисляет собственные значения комплексной эрмитовой или вещественной симметричной матрицы.

Нормы и другие числа

`linalg.norm`(x[, ord, axis, keepdims])	Норма матрицы или вектора.
`linalg.cond`(x[, p])	Вычисляет число обусловленности матрицы.
`linalg.det`(a)	Вычисляет определитель массива.
`linalg.matrix_rank`(M[, tol, hermitian])	Возвращает ранг матрицы массива, используя метод SVD.
`linalg.slogdet`(a)	Вычисляет знак и (натуральный) логарифм определителя массива.
`trace`(a[, offset, axis1, axis2, dtype, out])	Возвращает сумму по диагоналям массива.

Решение уравнений и обращение матриц

`linalg.solve`(a, b)	Решает линейное матричное уравнение или систему линейных скалярных уравнений.
`linalg.tensorsolve`(a, b[, axes])	Решает тензорное уравнение `a x = b` для x.
`linalg.lstsq`(a, b[, rcond])	Возвращает наименьшее решение в квадратных скобках для линейного матричного уравнения.
`linalg.inv`(a)	Вычисляет (мультипликативное) обращение матрицы.
`linalg.pinv`(a[, rcond, hermitian])	Вычисляет (обращение по Муру-Пенрозу) псевдообратную матрицу.
`linalg.tensorinv`(a[, ind])	Вычисляет «обращение» N-мерного массива.

Исключения

linalg.LinAlgError

Объект исключения Python, производный от исключения, поднимаемый функциями linalg.

Линейная алгебра для нескольких матриц сразу

Новое в версии 1.8.0.

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

Это указано в документации посредством спецификаций параметров ввода, таких как a : (..., M, M) array_like. Это означает, что, например, если на вход подается массив a.shape == (N, M, M), он интерпретируется как «стек» из N матриц, каждая размером M на M. Аналогичное описание применяется к возвращаемым значениям, например, определитель имеет det : (...) и в этом случае вернет массив формы det(a).shape == (N,). Это обобщается на линейные алгебраические операции над многомерными массивами: последние 1 или 2 измерения многомерного массива интерпретируются как векторы или матрицы соответственно операции.

© 2005–2020 NumPy Developers
Licensed under the 3-clause BSD License.
https://numpy.org/doc/1.19/reference/routines.linalg.html