Julia LinearAlgebra.lyap用法及代码示例

lyap(A, C)

计算连续 Lyapunov 方程 AX + XA' + C = 0 的解 X ，其中 A 的任何特征值都没有零实部，并且没有两个特征值是彼此的负复共轭。

例子

julia> A = [3. 4.; 5. 6]
2×2 Matrix{Float64}:
 3.0  4.0
 5.0  6.0

julia> B = [1. 1.; 1. 2.]
2×2 Matrix{Float64}:
 1.0  1.0
 1.0  2.0

julia> X = lyap(A, B)
2×2 Matrix{Float64}:
  0.5  -0.5
 -0.5   0.25

julia> A*X + X*A' + B
2×2 Matrix{Float64}:
 0.0          6.66134e-16
 6.66134e-16  8.88178e-16