Lapack的行減少

我正在嘗試編寫一個函數，爲一個代表性不足的系統的方程（例如描述系統的矩陣比它的高度更寬）產生單一的解決方案。爲了做到這一點，我一直在尋找LAPACK文檔，以便將矩陣行降低爲梯形圖，類似於Mathematica和TI計算器中的函數rref（）。我遇到的最接近的是http://software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/這個微小的線程。然而，這個線程似乎意味着簡單地採用「U」上三角矩陣（並將每一行除以對角線）與矩陣的簡化階梯形式相同，我不相信這是這種情況。我可以自己編寫rref（），但我不相信我可以實現LAPACK的出色表現。Lapack的行減少

1）有沒有更好的方法來簡單地得到任何一個具體解決方案的代表性不足的系統？
2）如果沒有，有沒有辦法讓LAPACK行減少矩陣？

謝謝！

來源

2009-07-20 Litherum

LU分解中的U矩陣是`ref`。用引導條目劃分每一行將把`ref`轉換爲`rref`。 – Sharpie 2010-06-29 16:26:53

一個常用的方法是最小二乘方法，參見lapack的sgelsx

來源

2009-07-21 15:28:21 second

我試過了，它似乎工作。謝謝。 – Litherum 2009-07-24 15:34:02

Lapack的行減少

回答

相關問題