numpy.linalg.solve

numpy.linalg.solve(a, b) [source]

Решает линейное матричное уравнение или систему линейных скалярных уравнений.

Вычисляет «точное» решение, x, хорошо определённого, т.е. с полным рангом, линейного матричного уравнения ax = b.

Параметры

a(…, M, M) array_like

Матрица коэффициентов.

b{(…, M,), (…, M, K)}, array_like

Значения ординаты или «зависимой переменной».

Возвращает

x{(…, M,), (…, M, K)} ndarray

Решение системы a x = b. Форма возвращаемого значения идентична форме b.

Возбуждает

LinAlgError

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

Примечания

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

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

Решения вычисляются с помощью процедуры LAPACK _gesv.

a должна быть квадратной и с полным рангом, т.е. все строки (или, что равносильно, столбцы) должны быть линейно независимыми; если это не так, используйте lstsq для получения наилучшего «решения» системы/уравнения по методу наименьших квадратов.

Ссылки

1

Г. Странг, Линейная алгебра и её приложения, 2-е изд., Орландо, Флорида, Издательство Academic Press, Inc., 1980, стр. 22.

Примеры

Решите систему уравнений 3 * x0 + x1 = 9 и x0 + 2 * x1 = 8:

>>> a = np.array([[3,1], [1,2]])
>>> b = np.array([9,8])
>>> x = np.linalg.solve(a, b)
>>> x
array([2.,  3.])

Проверьте, что решение правильное:

>>> np.allclose(np.dot(a, x), b)
True