使用 GNU 科学图书馆寻找矩阵的内核

Using the GNU Scientific Library to Find the Kernel of a Matrix

给定一个 gsl_matrix * 对象,它是一个 M-by-N 矩阵,最简单的方法是什么找到A的内核?我尝试使用奇异值分解(特别是 gsl_linalg_SV_decomp 方法),但是当 M<N 时它不起作用,这是有时我的矩阵就是这种情况。

最简单的方法(如果您受限于该库)是用零行填充您的矩阵到 N x N,然后您可以将其传递给您的函数。填充后的矩阵将具有相同的空 space.