Abstract
In this paper, we propose an extended block Krylov process to construct two biorthogonal bases for the extended Krylov subspaces \(\mathbb {K}_{m}^e(A,V)\) and \(\mathbb {K}_{m}^e(A^{T},W)\), where \(A \in \mathbb {R}^{n \times n}\) and \(V,~W \in \mathbb {R}^{n \times p}\). After deriving some new theoretical results and algebraic properties, we apply the proposed algorithm with moment matching techniques for model reduction in large scale dynamical systems. Numerical experiments for large and sparse problems are given to show the efficiency of the proposed method.
Similar content being viewed by others
Notes
Oberwolfach model reduction benchmark collection, 2003. http://www.imtek.de/simulation/benchmark.
References
Abidi, O., Heyouni, M., Jbilou, K.: On some properties of the extended block and global Arnoldi methods with applications to model reduction. Numer. Algorithms 75(1), 285–304 (2017)
Aliaga, J.I., Boley, D.L., Freund, R.W., Hernandez, V.: A Lanczos type method for multiple starting vectors. Math. Comput. 69(232), 1577–1601 (1999)
Antoulas, A.C., Sorensen, D.C., Gugercin, S.: A survey of model reduction methods for large scale systems. Contemp. Math. 280, 193–219 (2001)
Bai, Z., Day, D., Ye, Q.: ABLE: an adaptive block Lanczos method for non-Hermitian eigenvalue problems. SIAM J. Math. Anal. Appl. 20, 1060–1082 (1999)
Bai, Z.: Krylov subspace techniques for reduced-order modeling of large scale dynamical systems. Appl. Numer. Math. 43, 9–44 (2002)
Barkouki, H., Bentbib, A.H., Jbilou, K.: An adaptive rational block Lanczos-type algorithm for model reduction of large scale dynamical systems. J. Sci. Comput. 67, 221–236 (2016)
Beattie, C.A., Gugercin, S.: Krylov-based minimization for optimal \(\fancyscript {H}_{2}\) model reduction. In: Proceedings of the 46th IEEE Conference on Decision and Control, pp. 4385–4390 (2007)
Brezinski, C., Zaglia, M.Redivo, Sadok, H.: A breakdown-free Lanczos type algorithm for solving linear systems. Numer. Math. 63, 29–38 (1992)
Brezinski, C.: The block Lanczos and Vorobyev methods. C.R. Acad. Sci. Paris, Série I 331, 137–142 (2000)
Chan, T.F., de Pillis, L., van der Vorst, H.: Transpose-free formulations of Lanczos-type methods for nonsymmetric linear systems. Numer. Algorithms 17, 51–66 (1998)
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)
Frangos, M., Jaimoukha, I.M.: Adaptive rational interpolation: Arnoldi and Lanczos-like equations. Eur. J. Control 14(4), 342–354 (2008)
Freund, R.W., Gutknecht, M.H., Nachtigal, N.M.: An implementation of the look-ahead Lanczos algorithm for non-Hermitian matrices. SIAM J. Sci. Comput. 14, 137–158 (1993)
Gallivan, K., Grimme, E., Van Dooren, P.: Asymptotic waveform evaluation via a Lanczos method. Appl. Math. Lett. 7, 75–80 (1994)
Glover, K.: All optimal Hankel-norm approximation of linear multivariable systems and their \(L^{\infty }\)-error bounds. Int. J. Control 39(6), 1115–1193 (1984)
Grimme, E., Gallivan, K.: Rational Lanczos algorithm for model reduction II: interpolation point selection. Technical report, University of Illinois at Urbana Champaign (1998)
Grimme, E.: Krylov Projection methods for model reduction. Ph.D. thesis, ECE Department, University of Illinois, Urbana-Champaign (1997)
Gugercin, S., Antoulas, A.C.: Model reduction of large scale systems by least squares. Linear Algebra. Appl. 415(2–3), 290–321 (2006)
Gugercin, S., Antoulas, A.C., Beattie, C.A.: A Rational Krylov iteration for optimal \({\fancyscript {H}}_{2}\) model reduction. In: Proceedings of Mathematical Theory of Networks and Systems (2006)
Heyouni, M., Jbilou, K., Messaoudi, A., Tabaa, K.: Model in large scale MIMO dynamical systems via the block Lanczos method. Comput. Appl. Math. 27(2), 211–236 (2008)
Heyouni, M., Jbilou, K.: Matrix Krylov subspace methods for large scale model reduction problems. Appl. Math. Comput. 181, 1215–1228 (2006)
Heyouni, M., Jbilou, K.: An extended block Arnoldi algorithm for large-scale solutions of the continuous-time algebraic Riccati equation. Electron. Trans. Numer. Anal. 33, 53–62 (2009)
Lanczos, C.: An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. J. Res. Natl. Bur. Stand. 45, 225–280 (1950)
Mehrmann, V., Penzl, T.: Benchmark collections in SLICOT, Technical Report SLWN1998-5, SLICOT Working Note, ESAT, KU Leuven, K. Mercierlaan 94, Leuven-Heverlee 3100, Belgium (1998). http://www.win.tue.nl/niconet/NIC2/reports.html
Knizhnerman, L., Simoncini, V.: Convergence analysis of the extended Krylov subspace method for the Lyapunov equation. Numer. Math. 118(3), 567–586 (2011)
Parlett, B.N., Taylor, D.R., Li, Z.A.: A look-ahead Lanczos algorithm for unsymmetric matrices. Math. Comput. 44(169), 105–124 (1985)
Penzl, T.: LYAPACK-A MATLAB toolbox for Large Lyapunov and Riccati Equations, Model Reduction Problems and Linear-quadratic Optimal Control Problems. http://www.tu-chemnitz.de/sfb393/lyapack
Simoncini, V.: A new iterative method for solving large-scale Lyapunov matrix equations. SIAM J. Sci. Comput. 29, 1268–1288 (2007)
Wilkinson, J.H.: The Algebraic Eigenvalue Problem. Oxford Univ. Press, Oxford (1965)
Zhou, Y.: Numerical methods for larger scale matrix equations with applications in LTI system model reduction. Ph.D. thesis, CAAM Department, Rice University (2002)
Acknowledgements
We would like to thank the referees for valuable remarks and helpful suggestions.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Barkouki, H., Bentbib, A.H., Heyouni, M. et al. An extended nonsymmetric block Lanczos method for model reduction in large scale dynamical systems. Calcolo 55, 13 (2018). https://doi.org/10.1007/s10092-018-0248-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10092-018-0248-5