Affordable Access

Model Reduction of Second-Order Systems.

Springer Berlin Heidelberg
Publication Date
  • 15 Linear And Multilinear Algebra
  • Matrix Theory
  • 65 Numerical Analysis
  • 93 Systems Theory
  • Control
  • Computer Science
  • Engineering


Model Reduction of Second-Order Systems. Younes Chahlaoui, Kyle A Gallivan, Antoine Vandendorpe and Paul Van Dooren 2006 MIMS EPrint: 2008.11 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available from: And by contacting: The MIMS Secretary School of Mathematics The University of Manchester Manchester, M13 9PL, UK ISSN 1749-9097 6 Model Reduction of Second-Order Systems Y. Chahlaoui1, K. A. Gallivan1, A. Vandendorpe2, P. Van Dooren2 1 School of Computational Science, Florida State University, Tallahassee, U.S.A., 2 CESAME, Universite´ catholique de Louvain, Louvain-la-Neuve, Belgium,, 6.1 Introduction In this chapter, the problem of constructing a reduced order system while preserving the second-order structure of the original system is discussed. Af- ter a brief introduction on second-order systems and a review of first order model reduction techniques, two classes of second-order structure preserving model reduction techniques – Krylov subspace-based and SVD-based – are presented. For the Krylov techniques, conditions on the projectors that guar- antee the reduced second-order system tangentially interpolates the original system at given frequencies are derived and an algorithm is described. For SVD-based techniques, a Second-Order Balanced Truncation method is de- rived from second-order Gramians. Second-order systems arise naturally in many areas of engineering (see, for example, [Pre97, WJJ87]) and have the following form:{ Mq¨(t) +Dq˙(t) + Sq(t) = F in u(t), y(t) = F out q(t). (6.1) We assume that u(t) ∈ Rm, y(t) ∈ Rp, q(t) ∈ RN , F in ∈ RN×m, F out ∈ R p×N , and M,D,S ∈ RN×N with M invertible. For mechanical systems the matrices M , D and S represent, respectively, the mass (or inertia), damping and stiffness matrices, u(t) corresponds to

There are no comments yet on this publication. Be the first to share your thoughts.