Kalman's minimal realization theory involves geometric objects (controllable, unobservable subspaces) which are subject to structural instability. Specifically, arbitrarily small perturbations in a model may cause a change in the dimensions of the associated subspaces. This situation is manifested in computational difficulties which arise in attempts to apply textbook algorithms for computing a minimal realization. Structural instability associated with geometric theories is not unique to control; it arises in the theory of linear equations as well. In this setting, the computational problems have been studied for decades and excellent tools have been developed for coping with the situation. One of the main goals of this paper is to call attention to principal component analysis (Hotelling, 1933), and an algorithm (Golub and Reinsch, 1970) for computing the singular value decompositon of a matrix. Together they form a powerful tool for coping with structural instability in dynamic systems. As developed in this paper, principal component analysis is a technique for analyzing signals. (Singular value decomposition provides the computational machinery.) For this reason, Kalman's minimal realization theory is recast in terms of responses to injected signals. Application of the signal analysis to controllability and observability leads to a coordinate system in which the "internally balanced" model has special properties. For asymptotically stable systems, this yields working approximations of X_{c}, X_{\bar{o}} , the controllable and unobservable subspaces. It is proposed that a natural first step in model reduction is to apply the mechanics of minimal realization using these working subspaces.

Principal component analysis in linear systems: Controllability, observability, and model reduction

A method, called the Eigensystem Realization Algorithm (ERA), is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system modes and noise modes. For illustration of the algorithm, examples are shown using simulation data and experimental data for a rectangular grid structure.

/pdf/an-eigensystem-realization-algorithm-for-modal-parameter-1d3nimncg6.pdf

An eigensystem realization algorithm for modal parameter identification and model reduction

Subspace identification for linear systems

This paper reviews stochastic system identification methods that have been used to estimate the modal parameters of vibrating structures in operational conditions. It is found that many classical input-output methods havean output-only counterpart. For instance, the Complex Mode Indication Function (CMIF) can be applied both to Frequency Response Functions and output power and cross spectra. The Polyreference Time Domain (PTD) method applied to impulse responses is similar to the Instrumental Variable (IV) method applied to output covariances. The Eigensystem Realization Algorithm (ERA) is equivalent to stochastic subspace identification.

Stochastic System Identification for Operational Modal Analysis: A Review

In this paper, we present two novel algorithms to realize a finite dimensional, linear time-invariant state-space model from input-output data. The algorithms have a number of common features. They are classified as one of the subspace model identification schemes, in that a major part of the identification problem consists of calculating specially structured subspaces of spaces defined by the input-output data. This structure is then exploited in the calculation of a realization. Another common feature is their algorithmic organization: an RQ factorization followed by a singular value decomposition and the solution of an overdetermined set (or sets) of equations. The schemes assume that the underlying system has an output-error structure and that a measurable input sequence is available. The latter characteristic indicates that both schemes are versions of the MIMO Output-Error State Space model identification (MOESP) approach. The first algorithm is denoted in particular as the (elementary MOESP scheme)...

Subspace model identification-Part 1. The output-error state-space model identification class of algorithms

Ho's Algorithm generates exact realizations of finite-dimensional linear systems, given exact data. We discuss here a variant of the algorithm that generates approximate realizations of specified dimension from approximate data.

A. McEwen

Papers

Approximate linear realizations of given dimension via Ho's algorithm