如何解決numpy中的冗餘線性系統？

Question

我想解決的問題如下。我有一些具有冗餘性的線性方程 - Ax = s，其中A不是平方的，但有比列更多的行。任何建議，將不勝感激。

以下是我希望儘可能優雅地做的事情。我有一些有n行m列的矩陣A， n> = m。

我想生成的是來自A的方矩陣B，使得m行從A中選擇，而B是正則矩陣。

這樣我就可以得到規則矩陣，我可以使用numpy.linalg.solve來計算解決方案。

Answer 1

您可以使用np.linalg.lstsq來計算問題的最小二乘解（即，您不需要拋出任何數據，因爲您在提問時提出）。

這將發現矢量x，它最小化的距離| s - Ax |