SummaryThe condition number of a problem measures the sensitivity of the answer to small changes in the input. We call the problem ill-posed if its condition number is infinite. It turns out that for many problems of numerical analysis, there is a simple relationship between the condition number of a problem and the shortest distance from that problem to an ill-posed one: the shortest distance is proportional to the reciprocal of the condition number (or bounded by the reciprocal of the condition number). This is true for matrix inversion, computing eigenvalues and eigenvectors, finding zeros of polynomials, and pole assignment in linear control systems. In this paper we explain this phenomenon by showing that in all these cases, the condition number κ satisfies one or both of the diffrential inequalitiesm·κ2≤∥Dκ∥≤M·κ2, where ‖Dκ‖ is the norm of the gradient of κ. The lower bound on ‖Dκ‖ leads to an upper bound 1/mκ(x) on the distance. fromx to the nearest ill-posed problem, and the upper bound on ‖Dκ‖ leads to a lower bound 1/(Mκ(X)) on the distance. The attraction of this approach is that it uses local information (the gradient of a condition number) to answer a global question: how far away is the nearest ill-posed problem? The above differential inequalities also have a simple interpretation: they imply that computing the condition number of a problem is approximately as hard as computing the solution of the problem itself. In addition to deriving many of the best known bounds for matrix inversion, eigendecompositions and polynomial zero finding, we derive new bounds on the distance to the nearest polynomial with multiple zeros and a new perturbation result on pole assignment.
[1]
C. Eckart,et al.
A principal axis transformation for non-hermitian matrices
,
1939
.
[2]
P. Hartman.
Ordinary Differential Equations
,
1965
.
[3]
J. H. Wilkinson.
The algebraic eigenvalue problem
,
1966
.
[4]
W. Kahan.
Numerical Linear Algebra
,
1966,
Canadian Mathematical Bulletin.
[5]
Tosio Kato.
Perturbation theory for linear operators
,
1966
.
[6]
Axel Ruhe.
Properties of a matrix with a very ill-conditioned eigenproblem
,
1970
.
[7]
G. Stewart.
Error Bounds for Approximate Invariant Subspaces of Closed Linear Operators
,
1971
.
[8]
W. Kahan.
Conserving Confluence Curbs Ill-Condition
,
1972
.
[9]
J. H. Wilkinson.
Note on matrices with a very ill-conditioned eigenproblem
,
1972
.
[10]
W. Wonham.
Linear Multivariable Control: A Geometric Approach
,
1974
.
[11]
D. Hough.
Explaining and ameliorating the ill condition of zeros of polynomials.
,
1977
.
[12]
J. Demmel.
The condition number of equivalence transformations that block diagonalize matrix pencils
,
1983
.
[13]
J. H. Wilkinson.
On neighbouring matrices with quadratic elementary divisors
,
1984
.
[14]
B N Parlett,et al.
Development of an Accurate Algorithm for EXP(Bt). Appendix
,
1985
.
[15]
N. Nichols,et al.
Robust pole assignment in linear state feedback
,
1985
.
[16]
J. Demmel.
Computing stable eigendecompositions of matrices
,
1986
.