Python - 如何求解线性矩阵方程 AP + PA^T = Q for P (Q, A known)

Question

假设我有三个矩阵 A、P 和 Q，其中只有 Q 和 A 是已知的。如何使用 numpy/scipy?

求解 P

Answer 1

看起来像连续的Lyapunov equation. You can use the function scipy.linalg.solve_continuous_lyapunov。