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