使用GNU Scientific Library查找矩阵的核心

Question

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

Answer 1

最简单的方法（如果你只限于那个库）将填充你的矩阵用零行到N x N，然后你可以传递给你的函数。填充矩阵将具有相同的零空间。