Julia LinearAlgebra.rank用法及代码示例

rank(A::AbstractMatrix; atol::Real=0, rtol::Real=atol>0 ? 0 : n*ϵ)
rank(A::AbstractMatrix, rtol::Real)

通过计算A 有多少奇异值的幅度大于max(atol, rtol*σ₁) 来计算矩阵的秩，其中σ₁ 是A 的最大奇异值。 atol 和 rtol 分别是绝对公差和相对公差。默认的相对公差是 n*ϵ ，其中 n 是 A 的最小尺寸的大小，而 ϵ 是 A 的元素类型的 eps 。

例子

julia> rank(Matrix(I, 3, 3))
3

julia> rank(diagm(0 => [1, 0, 2]))
2

julia> rank(diagm(0 => [1, 0.001, 2]), rtol=0.1)
2

julia> rank(diagm(0 => [1, 0.001, 2]), rtol=0.00001)
3

julia> rank(diagm(0 => [1, 0.001, 2]), atol=1.5)
1