3

我正在尝试编写一个函数,该函数可以为代表性不足的方程组生成单一解(例如,描述系统的矩阵比它高)。为了做到这一点,我一直在 LAPACK 文档中寻找一种将矩阵行缩减为缩减梯队形式的方法,类似于 Mathematica 和 TI 计算器中的函数 rref()。我遇到的最接近的是http://software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/这个细线。然而,这个线程似乎暗示简单地采用“U”上三角矩阵(并将每一行除以对角线)与矩阵的简化梯形形式相同,我不认为是这种情况。我可以自己编写 rref() 代码,但我不相信我能达到 LAPACK 著名的性能。

1)有没有更好的方法来简单地为代表性不足的系统提供任何一种特定的解决方案?
2)如果没有,LAPACK有没有办法对矩阵进行行缩减?

谢谢!

4

1 回答 1

3

一种常用的方法是最小二乘法,参见 lapack 的 sgelsx

于 2009-07-21T15:28:21.207 回答