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