Moody T. Chu

Inverse Eigenvalue Problems

A collection of inverse eigenvalue problems are identified and classified according to their characteristics. Current developments in both the theoretic and the algorithmic aspects are summarized and reviewed in this paper. This exposition also reveals many open questions that deserve further study. An extensive bibliography of pertinent literature is attached.

Structured inverse eigenvalue problems

An inverse eigenvalue problem concerns the reconstruction of a structured matrix from prescribed spectral data. Such an inverse problem arises in many applications where parameters of a certain physical system are to be determined from the knowledge or expectation of its dynamical behaviour. Spectral information is entailed because the dynamical behaviour is often governed by the underlying natural frequencies and normal modes. Structural stipulation is designated because the physical system is often subject to some feasibility constraints. The spectral data involved may consist of complete or only partial information on eigenvalues or eigenvectors. The structure embodied by the matrices can take many forms. The objective of an inverse eigenvalue problem is to construct a matrix that maintains both the specific structure as well as the given spectral property. In this expository paper the emphasis is to provide an overview of the vast scope of this intriguing problem, treating some of its many applications, its mathematical properties, and a variety of numerical techniques.

The projected gradient methods for least squares matrix approximations with spectral constraints

The problems of computing least squares approximations for various types of real and symmetric matrices subject to spectral constraints share a common structure. This paper describes a general procedure in using the projected gradient method. It is shown that the projected gradient of the objective function on the manifold of constraints usually can be formulated explicitly. This gives rise to the construction of a descent flow that can be followed numerically. The explicit form also facilitates the computation of the second-order optimality conditions. Examples of applications are discussed. With slight modifications, the procedure can be extended to solve least squares problems for general matrices subject to singular-value constraints.

Structured low rank approximation

This paper concerns the construction of a structured low rank matrix that is nearest to a given matrix. The notion of structured low rank approximation arises in various applications, ranging from signal enhancement to protein folding to computer algebra, where the empirical data collected in a matrix do not maintain either the specified structure or the desirable rank as is expected in the original system. The task to retrieve useful information while maintaining the underlying physical feasibility often necessitates the search for a good structured lower rank approximation of the data matrix. This paper addresses some of the theoretical and numerical issues involved in the problem. Two procedures for constructing the nearest structured low rank matrix are proposed. The procedures are flexible enough that they can be applied to any lower rank, any linear structure, and any matrix norm in the measurement of nearness. The techniques can also be easily implemented by utilizing available optimization packages. The special case of symmetric Toeplitz structure using the Frobenius matrix norm is used to exemplify the ideas throughout the discussion. The concept, rather than the implementation details, is the main emphasis of the paper.

A rank-one reduction formula and its applications to matrix factorizations

Let

Parallel solution of ODE's by multi-block methods

A \in R^{m \times n}

ON THE CONTINUOUS REALIZATION OF ITERATIVE PROCESSES

denote an arbitrary matrix. If

Numerical methods for inverse singular value problems3

x \in R^n

On a multivariate eigenvalue problem, part I: algebraic theory and a power method

and

Linear algebra algorithms as dynamical systems

y \in R^m