Предположим, у меня есть три матрицы A, P и Q, где известны только Q и A. Как я могу решить для P, используя numpy / scipy?
Это выглядит как непрерывное уравнение Ляпунова . Вы можете использовать функцию scipy.linalg.solve_continuous_lyapunov.
scipy.linalg.solve_continuous_lyapunov