我有一个方阵(nxn),我知道它的等级< ñ。我还知道其列之一可以表示为少数其他列的线性组合(< n-1)。例如如果我有一个5x5矩阵,那么它的第2列可以表示为第1列和第3列的组合。 我想找到一种方法来检测这些线性相关的列。 https://stackoverflow.com/a/24548118/8245075上的答案确实可以完成我想要的工作。但是,我使用C编程并使用lapack例程进行线性代数。
如何在lapack中执行此操作?
答案 0 :(得分:0)
进行特征值分解(在此选择合适的:http://www.netlib.org/lapack/explore-html/d9/d8e/group__double_g_eeigen.html - 这些是真正的双精度驱动程序)。具有特征值0的cols(转置=>行)线性地依赖于其他。