Abstract
In the present paper, we propose a preconditioned Newton–Block Arnoldi method for solving large continuous time algebraic Riccati equations. Such equations appear in control theory, model reduction, circuit simulation amongst other problems. At each step of the Newton process, we solve a large Lyapunov matrix equation with a low rank right hand side. These equations are solved by using the block Arnoldi process associated with a preconditioner based on the alternating direction implicit iteration method. We give some theoretical results and report numerical tests to show the effectiveness of the proposed approach.
Similar content being viewed by others
References
Anderson, B.D.O., Moore, J.B.: Linear Optimal Control. Prentice-Hall, Englewood Cliffs, NJ (1971)
Arnold III, W.F., Laub, A.J.: Generalized eigenproblem algorithms and software for algebraic Riccati equations. Proc. IEEE 72, 1746–1754 (1984)
Bartels, R.H., Stewart, G.W.: Algorithm 432: solution of the matrix equation AX + XB \(=\) C. Circ. Syst. Signal Proc. 13, 820–826 (1994)
Bouhamidi, A., Hached, M., Heyouni, M., Jbilou, K.: A preconditioned block Arnoldi method for large Sylvester matrix equations. Numer. Linear Algebra Appl. 236(6), 1531–1542 (2011)
Datta, B.N.: Numerical Methods for Linear Control Systems. Elsevier, Amsterdam (2004)
Druskin, V., Simoncini, V.: Adaptive rational Krylov subspaces for large-scale dynamical systems. Syst. Control Lett. 60, 546–560 (2011)
Druskin, V., Knizhnerman, L., Simoncini, V.: Analysis of the rational Krylov subspace and ADI methods for solving the Lyapunov equation. SIAM J. Numer. Anal. 49, 1875–1898 (2011)
Feitzinger, F., Hylla, T., Sachs, E.W.: It inexact Kleinman–Newton method for Riccati equations. SIAM J. Matrix Anal. Appl. 31, 272–288 (2009)
Heyouni, M., Jbilou, K.: An extended block Arnoldi algorithm for large-scale solutions of the continuous-time algebraic Riccati equation. Elect. Trans. Numer. Anal. 33, 53–62 (2009)
Jaimoukha, I.M., Kasenally, E.M.: Krylov subspace methods for solving large Lyapunov equations. SIAM J. Matrix Anal. Appl. 31(1), 227–251 (1994)
Jbilou, K.: Block Krylov subspace methods for large continuous-time algebraic Riccati equations. Numer. Algorithms 34, 339–353 (2003)
Jbilou, K.: An Arnoldi based algorithm for large algebraic Riccati equations. Appl. Math. Lett. 19, 437–444 (2006)
Jbilou, K., Riquet, A.J.: Projection methods for large Lyapunov matrix equations. Linear Algebra Appl. 415(2), 344–358 (2006)
Jbilou, K.: ADI preconditioned Krylov methods for large Lyapunov matrix equations. Linear Algebra Appl. 432, 2473–2485 (2010)
Kelly, C.: Iterative Methods for Linear and Nonlinear Equations. SIAM, Philadelphia (1995)
Kleinman, D.L.: On an iterative technique for Riccati equation computations. IEEC Trans. Autom. Control 13, 114–115 (1968)
Lancaster, P., Rodman, L.: The Algebraic Riccati Equations. Clarendon Press, Oxford (1995)
Laub, A.J.: A Schur method for solving algebraic Riccati equations. IEEE Trans. Autom. Control AC 24, 913–921 (1979)
Li, J.-R., White, J.: Low-rank solutions of Lyapunov equations. SIAM J. Matrix Anal. Appl. 24(1), 60–280 (2002)
Lu, A., Wachspress, E.L.: Solution of Lyapunov equations by alternating direction implicit iteration. Comput. Math. Appl. 21, 43–58 (1991)
Mehrmann, V.: The autonomous linear quadratic control problem, theory and numerical solution. In: Number in Lecture Notes in Control and Information Sciences. Springer, Heidelberg (1991)
Peaceman, D., Rachford, H.: The numerical solution of elliptic and parabolic differential equations. J. Soc. Ind. Appl. Math. 3, 28–41 (1955)
Penzl, T.: A cyclic low-rank Smith method for large Lyapunov equations. SIAM J. Sci. Comput. 21(4), 1401–1418 (2000)
Penzl, T.: LYAPACK A MATLAB toolbox for Large Lyapunov and Riccati Equations, Model Reduction Problems, and Linear-quadratic Optimal Control Problems. http://www.tu-chemintz.de/sfb393/lyapack
Ruhe, A.: Rational Krylov sequence methods for eigenvalue computations. Linear Algebra Appl. 58, 391–405 (1984)
Saad, Y.: Numerical solution of large Lyapunov equations. In: Kaashoek, M.A., Van Shuppen, J.H., Ran, A.C. (eds.) Signal Processing, Scattering, Operator Theory and Numerical Methods, pp. 503–511. Birkhaser, Boston (1990)
Simoncini, V.: A new iterative method for solving large-scale Lyapunov matrix equations. SIAM J. Sci. Comput. 29, 1268–1288 (2007)
Van Dooren, P.: Gramian based model reduction of large-scale dynamical systems. In: Numerical Analysis, pp. 231–247, Chapman and Hall/CRC Press, London (2000)
Wachspress, E.L.: Iterative solution of the Lyapunov matrix equation. Appl. Math. Lett. 1, 87–90 (1988)
Wonham, W.M.: On a matrix Riccati equation of Stochastic control. SIAM J. Control 6, 681–697 (1968)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Bouhamidi, A., Hached, M. & Jbilou, K. A preconditioned block Arnoldi method for large scale Lyapunov and algebraic Riccati equations. J Glob Optim 65, 19–32 (2016). https://doi.org/10.1007/s10898-015-0317-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10898-015-0317-0