On Principal Angles between Subspaces of Euclidean Space

The cosines of the principal angles between the column spaces of full column rank matrices $X\in{\hbox{{\bf R}$^{m\times p}$}}$ and $Y\in{\hbox{{\bf R}$^{m\times q}$}}$ are efficiently computed, using the Bjorck--Golub algorithm, as the singular values of $Q_x^{T}Q_y,$ where Qx and Qy are orthonorma...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:SIAM journal on matrix analysis and applications 2000-01, Vol.22 (1), p.173-194
1. Verfasser: Drmac, Zlatko
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:The cosines of the principal angles between the column spaces of full column rank matrices $X\in{\hbox{{\bf R}$^{m\times p}$}}$ and $Y\in{\hbox{{\bf R}$^{m\times q}$}}$ are efficiently computed, using the Bjorck--Golub algorithm, as the singular values of $Q_x^{T}Q_y,$ where Qx and Qy are orthonormal matrices computed by the QR factorizations of X and Y, respectively. This paper shows that the Bjorck--Golub algorithm is mixed stable in the following sense: the computed singular values approximate with small relative error the exact cosines of the principal angles between the column spaces of $X+\Delta X$ and $Y+\Delta Y,$ where $\Delta X,$ $\Delta Y$ are small backward errors. Further, theoretical analysis and numerical evidence show that the algorithm becomes more robust if the QR factorizations are computed with the complete pivoting scheme of Powell and Reid. Moreover, it is shown that Gaussian elimination with complete pivoting can be used as an efficient preconditioner in computation and as a useful tool in analysis of the sensitivity of the QR factorization.
ISSN:0895-4798
1095-7162
DOI:10.1137/S0895479897320824