Julia LinearAlgebra.LDLt用法及代码示例

LDLt <: Factorization

实数 SymTridiagonal 矩阵S 的LDLt 矩阵分解类型，使得S = L*Diagonal(d)*L'，其中L 是 UnitLowerTriangular 矩阵，d 是向量。 LDLt 因式分解 F = ldlt(S) 的主要用途是用 F\b 求解方程的线性系统 Sx = b。这是 ldlt 的返回类型，对应的矩阵分解函数。

分解 F::LDLt 的各个组件可以通过 getproperty 访问：

零件	说明
F.L	L(单位下三角)LDLt的一部分
F.D	D(对角线)LDLt 的一部分
F.Lt	Lt(单位上三角)LDLt的一部分
F.d	D 的对角线值作为 Vector

例子

julia> S = SymTridiagonal([3., 4., 5.], [1., 2.])
3×3 SymTridiagonal{Float64, Vector{Float64}}:
 3.0  1.0   ⋅
 1.0  4.0  2.0
  ⋅   2.0  5.0

julia> F = ldlt(S)
LDLt{Float64, SymTridiagonal{Float64, Vector{Float64}}}
L factor:
3×3 UnitLowerTriangular{Float64, SymTridiagonal{Float64, Vector{Float64}}}:
 1.0        ⋅         ⋅
 0.333333  1.0        ⋅
 0.0       0.545455  1.0
D factor:
3×3 Diagonal{Float64, Vector{Float64}}:
 3.0   ⋅        ⋅
  ⋅   3.66667   ⋅
  ⋅    ⋅       3.90909