如何以更快的方式计算条件数？

Question

我正在 MATLAB 中进行优化以计算条件数，如下所示：

s = svd(H, 0);
cond = max(s) ./ min(s);

非常耗时，尤其是对于大型 H 矩阵。我正在寻找近似但快速的东西。

Answer 1

您正在计算所有奇异值并丢弃除其中两个以外的所有奇异值。你可以试试这个：

largestS = svds(H,1); %largest singular value
smallestS  = svds(H,1,0); %smallest singular value
cond = largestS/smallestS