Julia LinearAlgebra.nullspace用法及代碼示例

nullspace(M; atol::Real=0, rtol::Real=atol>0 ? 0 : n*ϵ)
nullspace(M, rtol::Real) = nullspace(M; rtol=rtol) # to be deprecated in Julia 2.0

通過包含 M 的奇異向量，其奇異值的幅度大於 max(atol, rtol*σ₁) ，其中 σ₁ 是 M 的最大奇異值，計算 M 的零空間的基礎。

默認情況下，相對公差 rtol 是 n*ϵ ，其中 n 是 M 的最小尺寸的大小，而 ϵ 是 M 元素類型的 eps 。

例子

julia> M = [1 0 0; 0 1 0; 0 0 0]
3×3 Matrix{Int64}:
 1  0  0
 0  1  0
 0  0  0

julia> nullspace(M)
3×1 Matrix{Float64}:
 0.0
 0.0
 1.0

julia> nullspace(M, rtol=3)
3×3 Matrix{Float64}:
 0.0  1.0  0.0
 1.0  0.0  0.0
 0.0  0.0  1.0

julia> nullspace(M, atol=0.95)
3×1 Matrix{Float64}:
 0.0
 0.0
 1.0