For a symmetric matrix $A$, the rationale for the success of applying QR to compute the spectral decomposition of $A = UDU^T$ is, for large $k$, the QR factorization of $A^k = Q_kR_k$ obeys, \begin{equation} Q_k\approx U. \end{equation} Could anyone help to justify the approximation?
Thanks, Jack