Abstract
In this paper, we propose a structure-preserving doubling algorithm (SDA) for the computation of the minimal nonnegative solution to the nonsymmetric algebraic Riccati equation (NARE), based on the techniques developed for the symmetric cases. This method allows the simultaneous approximation to the minimal nonnegative solutions of the NARE and its dual equation, requiring only the solutions to two linear systems and several matrix multiplications per iteration. Similar to Newton's method and the fixed-point iteration methods for solving NAREs, we also establish global convergence for SDA under suitable conditions, using only elementary matrix theory. We show that sequences of matrices generated by SDA are monotonically increasing and quadratically convergent to the minimal nonnegative solutions of the NARE and its dual equation. Numerical experiments show that the SDA algorithm is feasible and effective, and outperforms Newton's iteration and the fixed-point iteration methods.
Similar content being viewed by others
References
Anderson, B.: Second-order convergent algorithm for the steady-state Riccati equation. Int. J. Control 28, 295–306 (1978)
Bai, Z., Demmel, J., Gu, M.: An inverse free parallel spectral divide and conquer algorithm for nonsymmetric eigenproblems. Numer. Math. 76, 279–308 (1997)
Barlow, M., Rogers, L., Williams, D.: Wiener-Hopf factorization of matrices. In: Séminaire de Probabilités XIV, no. 784, Lecture Notes in Math., Springer-Verlag, 1980, pp. 324–331
Bartels, R., Stewart, G.: Solution of the matrix equation AX+XB=C. Comm. ACM 15, 820–826 (1972)
Bazarra, M., Sheraii, H., Shetty, C.: Nonlinear programming. John Wiley & Sons, 1993
Benner, P.: Contributions to the numerical solutions of algebraic Riccati equations and related eigenvalue problems. PhD Dissertation, Fakultät für Mathematik, TU Chemnitz-Zwickau, Chemnitz, Germany, 1997
Benner, P., Byers, R.: Disk functions and their relationship to the matrix sign function. In: Proc. European Control Conf. ECC 97, Paper 936, (CD-ROM). Waterloo, Belgium, 1997, BELWARE Information Technology
Benner, P., Byers, R.: Evaluating products of matrix pencils and collapsing matrix products. Num. Lin. Alg. Appl. 8, 357–380 (2001)
Berman, A., Plemmons, R.: Nonnegative Matrices in the Mathematical Sciences. Academic Press, 1979
Chu, E.-W., Fan, H.-Y., Lin, W.-W.: structure-preserving doubling algorithm for cotinuous-time algebraic Riccati equations. Lin. Alg. Appl. 396, 55–80 (2005)
Chu, E.-W., Fan, H.-Y., Lin, W.-W., Wang, C.-S.: Structure-preserving doubling algorithms for periodic discrete-time algebraic Riccati equations. Int. J. Control 77, 767–788 (2004)
Golub, G., Nash, S., Van Loan, C.: A Hessenberg-Schur method for the problem AX+XB=C. IEEE Trans. Auto. Control 24, 909–913 (1979)
Golub, G.H., Van Loan, C.F.: Matrix Computations, 3rd ed. The Johns Hopkins University Press, 1996
Guo, C.-H.: Nonsymmetric algebraic Riccati equations and Wiener-Hopf factorization for M-matrices. SIAM J. Matrix Analy. Appl. 23, 225–242 (2001)
Guo, C.-H., Laub, A.: On the iterative solution of a class of nonsymmetric algebraic Riccati equations. SIAM J. Matrix Analy. Appl. 22, 376–391 (2000)
Juang, J.: Existence of algebraic matrix Riccati equations arising in transport theory. Lin. Alg. Appl. 230, 89–100 (1995)
Juang, J., Lin, W.-W.: Nonsymmetric algebraic Riccati equations and Hamiltonian-like matrices. SIAM J. Matrix Analy. Appl. 20, 228–243 (1999)
Lin, W.-W., Xu, S.-F.: Analysis of structure-preserving doubling algorithms for Riccati-type matrix equations. SIAM J. Matrix Analy. Appl. 2005, to appear
London, R., Mckean, H., Rogers, L., Williams, D.: A martingale approach to some Wiener-Hopf problems II. In: Séminaire de Probabilités XVI. no. 920, Lecture Notes in Math., Springer-Verlag, 1982, pp. 68–90
Rogers, L.: Fluid models in queueing theory and Wiener-Hopf factorization of Markov chains. Ann. Appl. Probab. 4, 390–413 (1994)
Rogers, L., Shi, Z.: Computing the invariant law of a fluid model. J. Appl. Probab. 31, 885–896 (1994)
Sun, X., Quintana-Ortí, E.: Spectral division methods for block generalized Schur decompositions. Math. Comp. 73, 1827–1847 (2004)
Varga, R.: Matrix iterative analysis. Prentice-Hall, 1962
Williams, D.: A ``potential-theoretic'' note on the quadratic Wiener-Hopf equation for Q-matrices. In: Séminaire de Probabilités XVI. no. 920, Lecture Notes in Math., Springer-Verlag, 1982, pp. 91–94
Author information
Authors and Affiliations
Corresponding author
Additional information
This research was supported in part by RFDP (20030001103) & NSFC (10571007) of China and the National Center for Theoretical Sciences in Taiwan.
This author's research was supported by NSFC grant 1057 1007 and RFDP grant 200300001103 of China.
Rights and permissions
About this article
Cite this article
Guo, XX., Lin, WW. & Xu, SF. A structure-preserving doubling algorithm for nonsymmetric algebraic Riccati equation. Numer. Math. 103, 393–412 (2006). https://doi.org/10.1007/s00211-005-0673-7
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-005-0673-7