zbMATH — the first resource for mathematics

Numerical solutions to large-scale differential Lyapunov matrix equations. (English) Zbl 1416.65116
Summary: In the present paper, we consider large-scale differential Lyapunov matrix equations having a low rank constant term. We present two new approaches for the numerical resolution of such differential matrix equations. The first approach is based on the integral expression of the exact solution and an approximation method for the computation of the exponential of a matrix times a block of vectors. In the second approach, we first project the initial problem onto a block (or extended block) Krylov subspace and get a low-dimensional differential Lyapunov matrix equation. The latter differential matrix problem is then solved by the Backward Differentiation Formula method (BDF) and the obtained solution is used to build a low rank approximate solution of the original problem. The process is being repeated, increasing the dimension of the projection space until some prescribed accuracy is achieved. We give some new theoretical results and present numerical experiments.

65F30 Other matrix algorithms (MSC2010)
15A24 Matrix equations and identities
Full Text: DOI arXiv
[1] Abou-Kandil, Hisham; Freiling, Gerhard; Ionescu, Vlad; Jank, Gerhard, Applications to Robust Control, 467-493, (2003), Basel · Zbl 1027.93001
[2] Anderson, B.D.O., Moore, J.B.: Linear Optimal Control. Prentice-Hall, Englewood Cliffs (1971) · Zbl 0321.49001
[3] Bartels, RH; Stewart, GW, Algorithm 432: solution of the matrix equation AX + XB = c, Circ. Syst. Signal Proc., 13, 820-826, (1972) · Zbl 1372.65121
[4] Brezinski, C.: Computational Aspects of Linear Control. Kluwer, Dordrecht (2002)
[5] Corless, M.J., Frazho, A.E.: Linear Systems and Control—An Operator Perspective. Pure and Applied Mathematics. Marcel Dekker, New York-Basel (2003)
[6] Datta, B.N.: Numerical Methods for Linear Control Systems Design and Analysis. Elsevier Academic Press, Amsterdam (2003)
[7] Druskin, V.; Knizhnerman, L., Extended Krylov subspaces: approximation of the matrix square root and related functions, SIAM J. Matrix Anal. Appl., 19, 755-771, (1998) · Zbl 0912.65022
[8] El Guennouni, A.; Jbilou, K.; Riquet, AJ, Block Krylov subspace methods for solving large Sylvester equations, Numer. Alg., 29, 75-96, (2002) · Zbl 0992.65040
[9] Gallopoulos, E.; Saad, Y., Efficient solution of parabolic equations by Krylov approximation methods, SIAM J. Sci. Statist. Comput., 13, 1236-1264, (1992) · Zbl 0757.65101
[10] Glover, K., All optimal Hankel-norm approximations of linear multivariable systems and their L-infinity error bounds, Int. J. Control., 39, 1115-1193, (1984) · Zbl 0543.93036
[11] Golub, GH; Nash, S.; Loan, C., A Hessenberg Schur method for the problem AX + XB = C, IEEE Trans. Automat. Contr., 24, 909-913, (1979) · Zbl 0421.65022
[12] Higham, NJ, The scaling and squaring method for the matrix exponential revisited, SIAM J. Matrix Anal. Appl., 26, 1179-1193, (2005) · Zbl 1081.65037
[13] Higham, NJ; Al-Mohy, A., A new scaling and squaring algorithm for the matrix exponential, SIAM J. Matrix Anal. Appl., 31, 970-989, (2009) · Zbl 1194.15021
[14] Hu, DY; Reichel, L., Krylov-subspace methods for the Sylvester equation, Lin. Alg. Appl., 172, 283-313, (1992) · Zbl 0777.65028
[15] Jaimoukha, IM; Kasenally, EM, Krylov subspace methods for solving large Lyapunov equations, SIAM J. Numer. Anal., 31, 227-251, (1994) · Zbl 0798.65060
[16] Jbilou, K., Low-rank approximate solution to large Sylvester matrix equations, App. Math. Comput., 177, 365-376, (2006) · Zbl 1095.65041
[17] Jbilou, K.; Riquet, AJ, Projection methods for large Lyapunov matrix equations, Lin. Alg. Appl., 415, 344-358, (2006) · Zbl 1094.65039
[18] Moler, CB; Loan, CF, Nineteen dubious ways to compute the exponential of a matrix, SIAM Rev., 20, 801836, (1978) · Zbl 0395.65012
[19] Moore, BC, Principal component analysis in linear systems: controllability, observability and model reduction, IEEE Trans. Autom. Control, AC-26, 17-32, (1981) · Zbl 0464.93022
[20] Mullis, CT; Roberts, RA, Synthesis of minimum roundoff noise fixed point digital filters., IEEE Trans. Circ. Syst., CAS-23, 551-562, (1996) · Zbl 0342.93066
[21] Saad, Y.: Numerical solution of large Lyapunov equations, in Signal Processing, Scattering, Operator Theory and Numerical Methods. In: Kaashoek, M.A., van Schuppen, J.H., Ran, A.C. (eds.) Proceedings of the international symposium MTNS-89, vol. 3, pp 503-511. Birkhauser, Boston (1990)
[22] Saad, Y., Analysis of some Krylov subspace approximations to the matrix exponential operator, SIAM J. Numer. Anal., 29, 209-228, (1992) · Zbl 0749.65030
[23] Sandberg, H.: Linear Time-Varying Systems: Modeling and Reduction. Ph.D. thesis, Department of Automatic Control, Lund Institute of Technology, Lund (2002)
[24] Shokoohi, S.; Silverman, L.; Dooren, P., Linear time-variable systems: balancing and model reduction, IEEE Trans. Autom. Control, 28, 810-822, (1983) · Zbl 0517.93014
[25] Simoncini, V., A new iterative method for solving large-scale Lyapunov matrix equations, SIAM J. Sci. Comp., 29, 1268-1288, (2007) · Zbl 1146.65038
[26] Verriest, EI; Kailath, T., On generalized balanced realizations, IEEE Trans. Autom. Control, 28, 833-844, (1983) · Zbl 0528.93021
[27] van der Vorst, H.: Iterative Krylov Methods for Large Linear Systems. Cambridge University Press, Cambridge (2003) · Zbl 1023.65027
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.