TOPICS
Search

Kahan Matrix


A Kahan matrix is an upper triangular matrix used to test matrix rank estimation from a QR decomposition (Kahan 1966, pp. 791-792). An unperturbed n×n version has entries

 a_(ij)={s^(i-1)   if i=j; -cs^(i-1)   if i<j; 0   if i>j,
(1)

where s=sintheta, c=costheta, and 0<theta<pi. Each column has unit Euclidean norm. Thus column choices in a pivoted QR decomposition can leave the columns in their original order, with Q=I and R=A as a factorization. The last diagonal entry can substantially overestimate the smallest singular value. Implementations may add a small diagonal matrix to preserve the intended column order in floating-point arithmetic (Higham 1995, pp. 30-31 and 49).

TestMatrixTableKahan

The plot shows |a_(nn)|/sigma_(min)(A) for the unperturbed family with theta=pi/3. Its growth illustrates why the size of a diagonal entry alone need not reveal the smallest singular value.


See also

Matrix Rank, QR Decomposition, Singular Value, Test Matrix

Explore with Wolfram|Alpha

References

Higham, N. J. The Test Matrix Toolbox for Matlab (Version 3.0). Numerical Analysis Report 276. Manchester, England: Manchester Centre for Computational Mathematics, 1995. https://nhigham.com/wp-content/uploads/2023/10/high95m.pdf.Kahan, W. "Numerical Linear Algebra." Canad. Math. Bull. 9, 757-801, 1966. https://doi.org/10.4153/CMB-1966-083-2.

Cite this as:

Weisstein, Eric W. "Kahan Matrix." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/KahanMatrix.html

Subject classifications