Skip to main content
Log in

A preconditioned block Arnoldi method for large scale Lyapunov and algebraic Riccati equations

  • Published:
Journal of Global Optimization Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1

Similar content being viewed by others

References

  1. Anderson, B.D.O., Moore, J.B.: Linear Optimal Control. Prentice-Hall, Englewood Cliffs, NJ (1971)

    MATH  Google Scholar 

  2. Arnold III, W.F., Laub, A.J.: Generalized eigenproblem algorithms and software for algebraic Riccati equations. Proc. IEEE 72, 1746–1754 (1984)

    Article  Google Scholar 

  3. Bartels, R.H., Stewart, G.W.: Algorithm 432: solution of the matrix equation AX + XB \(=\) C. Circ. Syst. Signal Proc. 13, 820–826 (1994)

    Google Scholar 

  4. 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)

    MathSciNet  MATH  Google Scholar 

  5. Datta, B.N.: Numerical Methods for Linear Control Systems. Elsevier, Amsterdam (2004)

    MATH  Google Scholar 

  6. Druskin, V., Simoncini, V.: Adaptive rational Krylov subspaces for large-scale dynamical systems. Syst. Control Lett. 60, 546–560 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  7. 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)

    Article  MathSciNet  MATH  Google Scholar 

  8. Feitzinger, F., Hylla, T., Sachs, E.W.: It inexact Kleinman–Newton method for Riccati equations. SIAM J. Matrix Anal. Appl. 31, 272–288 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  9. 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)

    MathSciNet  MATH  Google Scholar 

  10. Jaimoukha, I.M., Kasenally, E.M.: Krylov subspace methods for solving large Lyapunov equations. SIAM J. Matrix Anal. Appl. 31(1), 227–251 (1994)

    MathSciNet  MATH  Google Scholar 

  11. Jbilou, K.: Block Krylov subspace methods for large continuous-time algebraic Riccati equations. Numer. Algorithms 34, 339–353 (2003)

    Article  MathSciNet  MATH  Google Scholar 

  12. Jbilou, K.: An Arnoldi based algorithm for large algebraic Riccati equations. Appl. Math. Lett. 19, 437–444 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  13. Jbilou, K., Riquet, A.J.: Projection methods for large Lyapunov matrix equations. Linear Algebra Appl. 415(2), 344–358 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  14. Jbilou, K.: ADI preconditioned Krylov methods for large Lyapunov matrix equations. Linear Algebra Appl. 432, 2473–2485 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  15. Kelly, C.: Iterative Methods for Linear and Nonlinear Equations. SIAM, Philadelphia (1995)

    Book  Google Scholar 

  16. Kleinman, D.L.: On an iterative technique for Riccati equation computations. IEEC Trans. Autom. Control 13, 114–115 (1968)

    Article  Google Scholar 

  17. Lancaster, P., Rodman, L.: The Algebraic Riccati Equations. Clarendon Press, Oxford (1995)

    MATH  Google Scholar 

  18. Laub, A.J.: A Schur method for solving algebraic Riccati equations. IEEE Trans. Autom. Control AC 24, 913–921 (1979)

    Article  MathSciNet  MATH  Google Scholar 

  19. Li, J.-R., White, J.: Low-rank solutions of Lyapunov equations. SIAM J. Matrix Anal. Appl. 24(1), 60–280 (2002)

    MathSciNet  Google Scholar 

  20. Lu, A., Wachspress, E.L.: Solution of Lyapunov equations by alternating direction implicit iteration. Comput. Math. Appl. 21, 43–58 (1991)

    Article  MathSciNet  MATH  Google Scholar 

  21. 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)

  22. Peaceman, D., Rachford, H.: The numerical solution of elliptic and parabolic differential equations. J. Soc. Ind. Appl. Math. 3, 28–41 (1955)

    Article  MathSciNet  MATH  Google Scholar 

  23. Penzl, T.: A cyclic low-rank Smith method for large Lyapunov equations. SIAM J. Sci. Comput. 21(4), 1401–1418 (2000)

    Article  MathSciNet  MATH  Google Scholar 

  24. 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

  25. Ruhe, A.: Rational Krylov sequence methods for eigenvalue computations. Linear Algebra Appl. 58, 391–405 (1984)

    Article  MathSciNet  MATH  Google Scholar 

  26. 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)

    Google Scholar 

  27. Simoncini, V.: A new iterative method for solving large-scale Lyapunov matrix equations. SIAM J. Sci. Comput. 29, 1268–1288 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  28. 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)

  29. Wachspress, E.L.: Iterative solution of the Lyapunov matrix equation. Appl. Math. Lett. 1, 87–90 (1988)

    Article  MathSciNet  MATH  Google Scholar 

  30. Wonham, W.M.: On a matrix Riccati equation of Stochastic control. SIAM J. Control 6, 681–697 (1968)

    Article  MathSciNet  MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to K. Jbilou.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10898-015-0317-0

Keywords

Navigation