Results 1 
6 of
6
Krylov Projection Methods For Model Reduction
, 1997
"... This dissertation focuses on efficiently forming reducedorder 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 p ..."
Abstract

Cited by 118 (3 self)
 Add to MetaCart
This dissertation focuses on efficiently forming reducedorder 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 Lanczosbased 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 modelreduction package. Techniques for choosing the matching frequencies, estimating the modeling error, insuring the model's stability, treating multipleinput multipleoutput systems, implementing parallelism, and avoiding a need for exact factors of large matrix pencils are all examined to various degrees.
Model reduction of state space systems via an Implicitly Restarted Lanczos method
 Numer. Algorithms
, 1996
"... The nonsymmetric Lanczos method has recently received significant attention as a model reduction technique for largescale systems. Unfortunately, the Lanczos method may produce an unstable partial realization for a given, stable system. To remedy this situation, inexpensive implicit restarts are de ..."
Abstract

Cited by 56 (8 self)
 Add to MetaCart
The nonsymmetric Lanczos method has recently received significant attention as a model reduction technique for largescale systems. Unfortunately, the Lanczos method may produce an unstable partial realization for a given, stable system. To remedy this situation, inexpensive implicit restarts are developed which can be employed to stabilize the Lanczos generated model.
Padé Approximation Of LargeScale Dynamic Systems With Lanczos Methods
, 1994
"... The utility of Lanczos methods for the approximation of largescale 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 ..."
Abstract

Cited by 18 (1 self)
 Add to MetaCart
The utility of Lanczos methods for the approximation of largescale 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 multipoint 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 reducedorder modeling and simulation of largescale, 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 directcoupling term, d, will be assumed to be zero. The system matrix, A 2 R n\Thetan ...
On some modifications of the Lanczos algorithm and the relation with Padé approximations
, 1995
"... 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 ..."
Abstract

Cited by 6 (0 self)
 Add to MetaCart
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 nth 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 wellknown 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...
The Lanczos algorithm and Padé approximations
 Short Course, Benelux Meeting on Systems and Control
, 1995
"... 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 give ..."
Abstract

Cited by 4 (0 self)
 Add to MetaCart
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 nth 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 wellknown that the transfer function of this system : h(s) = c(sI \Gamma A) \Gamma1 b +<F29
Computing Minimal Partial Realizations Via a LanczosType Algorithm for Multiple Starting Vectors
 Proc. 36th IEEE CDC
, 1997
"... We describe a Lanczostype procedure that reduces a given realization of a finite sequence of (moment) matrices to a minimal partial realization. A key feature of this procedure is that the underlying Lanczostype algorithm is directly applied to the matrix triplet describing the given realization, ..."
Abstract

Cited by 3 (2 self)
 Add to MetaCart
We describe a Lanczostype procedure that reduces a given realization of a finite sequence of (moment) matrices to a minimal partial realization. A key feature of this procedure is that the underlying Lanczostype algorithm is directly applied to the matrix triplet describing the given realization, rather than to the moment matrices. It thus avoids explicit formulation of and the usually unstable computation with the moment matrices. 1 Introduction Given a finite sequence of matrices, M 0 ; M 1 ; : : : ; M r\Gamma1 2 C p\Thetam ; (1) the partial realization problem is to find a triplet of matrices A n 2 C n\Thetan ; B n 2 C n\Thetam ; C n 2 C n\Thetap ; (2) such that M i = C T n A i n B n for all i = 0; 1; : : : ; r \Gamma 1: (3) A partial realization (2) is called minimal if the statespace dimension n in (2) and (3) is as small as possible. The minimal partial realization problem has been studied extensively; see, e.g., [2, 3, 5, 6, 7, 8, 15, 18, 22, 23] and the refere...