This dissertation focuses on efficiently forming reduced-order models for large, linear dynamic systems. Projections onto unions of Krylov subspaces lead to a class of reducedorder models known as rational interpolants. The cornerstone of this dissertation is a collection of theory relating Krylov projection to rational interpolation. Based on this theoretical framework, three algorithms for model reduction are proposed. The first algorithm, dual rational Arnoldi, is a numerically reliable approach involving orthogonal projection matrices. The second, rational Lanczos, is an efficient generalization of existing Lanczos-based methods. The third, rational power Krylov, avoids orthogonalization and is suited for parallel or approximate computations. The performance of the three algorithms is compared via a combination of theory and examples. Independent of the precise algorithm, a host of supporting tools are also developed to form a complete model-reduction package. Techniques for choosing the matching frequencies, estimating the modeling error, insuring the model's stability, treating multiple-input multiple-output systems, implementing parallelism, and avoiding a need for exact factors of large matrix pencils are all examined to various degrees.

The utility of Lanczos methods for the approximation of large-scale dynamical systems is considered. In particular, it is shown that the Lanczos method is a technique for yielding Pad'e approximants which has several advantages over more traditional explicit moment matching approaches. An extension of the Lanczos algorithm is developed for computing multi-point Pad'e approximations of descriptor systems. Keywords: Dynamic system, Pad'e approximation, Lanczos algorithm, model reduction. 1. Introduction This paper explores the use of Lanczos techniques for the reduced-order modeling and simulation of large-scale, SISO dynamical systems. One can define such a system through the set of state space equations ae E x(t) = Ax(t) + bu(t) y(t) = cx(t) + du(t): (1) The scalar functions u(t) and y(t) are the system's input and output while x(t) is the state vector of dimension n. For simplicity, the direct-coupling term, d, will be assumed to be zero. The system matrix, A 2 R n\Thetan ...

In this paper, we describe some recent developments in the use of projection methods to produce reduced-order models for linear time-invariant dynamic systems. Previous related efforts in model reduction problems from various applications are also discussed. An overview is given of the theory governing the definition of the family of Rational Krylov methods, the practical heuristics involved and the important future research directions.

In this paper we try to show the relations between the Lanczos algorithm and Pad'e approximations as used e.g. in identification and model reduction of dynamical systems. 1 1 Introduction For simplicity we assume here that all systems are SISO, although some results do extend to the MIMO case. Let a n-th order dynamical system be described by x = Ax + bu (1.1) y = cx + du (1.2) where A is a square, b is a column vector, c is a row vector, and d is a scalar. It is well-known that the transfer function of this system : h(s) = c(sI \Gamma A) \Gamma1 b + d has a Taylor expansion around s = 1 that looks like : h(s) = d + cbs \Gamma1 + cAbs \Gamma2 + cA 2 bs \Gamma3 + cA 3 bs \Gamma4 + : : : The coefficients m \Gammai of the powers of s \Gammai satisfy thus m 0 = d ; m \Gammai = cA i\Gamma1 b ; i 1: For i 1 these are also called moments or Markov parameters of the system fA; b; cg. It follows already from the work of Hankel that the first 2n moments 1 To appea...

Introduction In these two lectures we try to show the relations between the Lanczos algorithm and Pad'e approximations as used e.g. in identification and model reduction of dynamical systems. These notes are based on material in the papers [10, 17, 11, 12] for which a lot of credit ought to be given to the respective coauthors. For simplicity we assume here that all systems are SISO, although some results do extend to the MIMO case. Let a n-th order dynamical system be described by x = Ax + bu (1) y = cx + du (2) where A is a square, b is a column vector, c is a row vector, and d is a scalar. It is well-known that the transfer function of this system : h(s) = c(sI \Gamma A) \Gamma1 b +<F29