用xA=b解决scipy.linalg.solve_triangular问题

Question

我想使用scipy.linalg.solve_triangular()来解决表单xA=b (而不是Ax=b)的系统。有什么简单的方法吗？我想在使用这个函数之前，我可以把所有的东西都转过来，但似乎不起作用。任何帮助感谢，是一个初学者在这方面！

Answer 1

向量x必须是(1xn)；矩阵A必须是( n )；向量b必须是(1xm)。

如果你接受两边的转位，你会得到：

(xA)^T = b^T

重新安排LHS：

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

现在A^T是(m )矩阵；x是(Nx1)向量；b是(Mx1)向量。

如果A是正方形和对称的，那么根据定义A^T = A.不需要做任何工作。

您可以使用通常的技术来解决x^T = (A^T)^-1 (b^T)问题。

我不建议计算矩阵逆。如果你的矩阵是正方形的，你最好用LU分解和前向替换代替。它要稳定得多。