如何解决numpy中的冗余线性系统？

Question

我试图解决的问题如下。我有一些带冗余的线性方程 - Ax = s其中A不是平方但有多行而不是列。任何建议将不胜感激。

以下是我想尽可能优雅地做的事情。我有一些矩阵A，有n行和m列; n> = m。

我想要生成的是来自A的方阵B，使得从A中选择m行，并且B是规则矩阵。

这样我就可以获得常规矩阵，我可以使用numpy.linalg.solve来计算解决方案。

Answer 1

您可以使用np.linalg.lstsq来计算问题的最小二乘解决方案（即，您不需要像在提问时那样抛出任何数据）。

这将找到向量 x ，它最小化距离 | s - Axe |