Matrix 拉帕克'；s行缩减

我正在尝试编写一个函数，该函数为一个未充分表示的方程组生成一个解（例如，描述该系统的矩阵比它高的矩阵宽）。为了做到这一点，我在LAPACK文档中寻找一种将矩阵行缩减为其缩减梯队形式的方法，类似于Mathematica和TI计算器中的函数rref（）。我遇到的最接近的是这根细线。然而，这条线索似乎暗示，简单地取“U”上三角矩阵（并将每行除以对角线）与矩阵的简化梯队形式相同，我不认为是这样。我可以自己编写rref（），但我不相信我能达到LAPACK著名的性能

1） 对于代表性不足的系统，是否有更好的方法可以简单地获得任何一种特定的解决方案？
2） 如果没有，LAPACK是否有办法减少矩阵的行数

---

谢谢

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

LU分解中的U矩阵是a

ref

。将每行除以前导项将把

ref

转换为

rref

。