Abstract:In this paper, we focus on using optimization methods to solve matrix equations by transforming the problem of solving the Sylvester matrix equation or continuous algebraic Riccati equation into an optimization problem. Initially, we use a constrained convex optimization method (CCOM) to solve the Sylvester matrix equation with $\ell_{2,1}$-norm, where we provide a convergence analysis and numerical examples of CCOM; however, the results show that the algorithm is not efficient. To address this issue, we employ classical quasi-Newton methods such as DFP and BFGS algorithms to solve the Sylvester matrix equation and present the convergence and numerical results of the algorithm. Additionally, we compare these algorithms with the CG algorithm and AR algorithm, and our results demonstrate that the presented algorithms are effective. Furthermore, we propose a unified framework of the alternating direction multiplier method (ADMM) for directly solving the continuous algebraic Riccati equation (CARE), and we provide the convergence and numerical results of ADMM. Our experimental results indicate that ADMM is an effective optimization algorithm for solving CARE. Finally, to improve the effectiveness of the optimization method for solving Riccati equation, we propose the Newton-ADMM algorithm framework, where the outer iteration of this method is the classical Newton method, and the inner iteration involves using ADMM to solve Lyapunov matrix equations inexactly. We also provide the convergence and numerical results of this algorithm, which our results demonstrate are more efficient than ADMM for solving CARE.

Low-rank generalized alternating direction implicit iteration method for solving matrix equations

Low-rank generalized alternating direction implicit iteration method for solving matrix equations

A general alternating direction implicit iteration method for solving complex symmetric linear systems

Low-rank alternating direction doubling algorithm for solving large-scale continuous time algebraic Riccati equations

An extension of the low-rank Lyapunov ADI to non-zero initial values and its applications

An implicit iterative algorithm for generalized algebraic Riccati matrix equations

An RADI-type method for stochastic continuous-time algebraic Riccati equations

An Implicit Iterative Algorithm for Algebraic Riccati Matrix Equations

Low-Rank Methods for Solving Discrete-Time Projected Lyapunov Equations

An increasing rank Riemannian method for generalized Lyapunov equations

ISC: an RADI-type Method for Stochastic Continuous-Time Algebraic Riccati Equations.

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

On a family of low-rank algorithms for large-scale algebraic Riccati equations

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

Alternately linearized implicit iteration methods for the minimal nonnegative solutions of the nonsymmetric algebraic Riccati equations

The relaxed gradient based iterative algorithm for solving the generalized coupled complex conjugate and transpose Sylvester matrix equations

Optimization methods for solving matrix equations

Preconditioned Low-Rank Riemannian Optimization for Symmetric Positive Definite Linear Matrix Equations

Gradient Based Iterative Algorithm for Solving Coupled Matrix Equations

A class of ADI method for the nonsymmetric coupled algebraic Riccati equations