k =rank(A) %a is matrix s = svd(A); tol = max(size(A))*eps(max(s)); r = sum(s > tol); 转载于:https://www.cnblogs.com/vigorz/p/10499128.html