190 research outputs found
The Anderson model of localization: a challenge for modern eigenvalue methods
We present a comparative study of the application of modern eigenvalue
algorithms to an eigenvalue problem arising in quantum physics, namely, the
computation of a few interior eigenvalues and their associated eigenvectors for
the large, sparse, real, symmetric, and indefinite matrices of the Anderson
model of localization. We compare the Lanczos algorithm in the 1987
implementation of Cullum and Willoughby with the implicitly restarted Arnoldi
method coupled with polynomial and several shift-and-invert convergence
accelerators as well as with a sparse hybrid tridiagonalization method. We
demonstrate that for our problem the Lanczos implementation is faster and more
memory efficient than the other approaches. This seemingly innocuous problem
presents a major challenge for all modern eigenvalue algorithms.Comment: 16 LaTeX pages with 3 figures include
Structured Eigenvalue Backward Errors of Matrix Pencils and Polynomials with Hermitian and Related Structures
Structure preservation: a challenge in computational control
Current and future directions in the development of numerical methods and numerical software for control problems are discussed. Major challenges include the demand for higher accuracy, robustness of the method with respect to uncertainties in the data or the model, and the need for methods to solve large scale problems. To address these demands it is essential to preserve any underlying physical structure of the problem. At the same time, to obtain the required accuracy it is necessary to avoid all inversions or unnecessary matrix products. We will demonstrate how these demands can be met to a great extent for some important tasks in control, the linear-quadratic optimal control problem for first and second order systems as well as stability radius and H-infinity norm computations. (C) 2003 Elsevier Science B.V All rights reserved
Lengths of quasi-commutative pairs of matrices
In this paper we discuss some partial solutions of the length conjecture which describes the length of a generating system for matrix algebras. We consider mainly the algebras generated by two matrices which are quasi-commuting. It is shown that in this case the length function is linearly bounded. We also analyze which particular natural numbers can be realized as the lengths of certain special generating sets and prove that for commuting or product-nilpotent pairs all possible numbers are realizable, however there are non-realizable values between lower and upper bounds for the other quasi-commuting pairs. In conclusion we also present several related open problems
Decay rate estimations for linear quadratic optimal regulators
Let be the optimal control of the open-loop system
in a linear quadratic optimization problem. By using
different complex variable arguments, we give several lower and upper estimates
of the exponential decay rate of the closed-loop system .
Main attention is given to the case of a skew-Hermitian matrix .
Given an operator , for a class of cases, we find a matrix that
provides an almost optimal decay rate.
We show how our results can be applied to the problem of optimizing the decay
rate for a large finite collection of control systems , , and illustrate this on an example of a concrete mechanical system. At the
end of the article, we pose several questions concerning the decay rates in the
context of linear quadratic optimization and in a more general context of the
pole placement problem.Comment: 25 pages, 1 figur
The level set method for the two-sided eigenproblem
We consider the max-plus analogue of the eigenproblem for matrix pencils
Ax=lambda Bx. We show that the spectrum of (A,B) (i.e., the set of possible
values of lambda), which is a finite union of intervals, can be computed in
pseudo-polynomial number of operations, by a (pseudo-polynomial) number of
calls to an oracle that computes the value of a mean payoff game. The proof
relies on the introduction of a spectral function, which we interpret in terms
of the least Chebyshev distance between Ax and lambda Bx. The spectrum is
obtained as the zero level set of this function.Comment: 34 pages, 4 figures. Changes with respect to the previous version: we
explain relation to mean-payoff games and discrete event systems, and show
that the reconstruction of spectrum is pseudopolynomia
- …