用numpy / scipy求解A x = b

Question

是否有可能使用numpy.linalg或scipy.linalg中的工具来求解方程式中的A

$$ A x = b $$

其中$ A $，$ x $和$ b $是$ n \乘以n $个矩阵（$ b $是对角矩阵）？本质上，我正在尝试“逆求解”一组标准的线性方程。

Answer 1

应用一些线性代数。

Ax = b
(Ax)^T = b^T
x^T A^T = b^T

用A^T解决numpy.linalg.solve，并转置：

A = numpy.linalg.solve(x.T, b.T).T