How long does it take to compute the eigenvalues of a random, symmetric matrix?
How long does it take to compute the eigenvalues of a random, symmetric matrix?
复制标题
计算随机对称矩阵的特征值需要多长时间?
DOI:
10.14288/1.0319078
复制
发表时间:
2012
期刊:
影响因子:
--
通讯作者:
Govind Menon
中科院分区:
文献类型:
--
作者:
C. Pfrang;P. Deift;Govind Menon
We present the results of an empirical study of the performance of the QR algorithm (with and without shifts) and the Toda algorithm on random symmetric matrices. The random matrices are chosen from six ensembles, four of which lie in the Wigner class. For all three algorithms, we observe a form of universality for the deflation time statistics for random matrices within the Wigner class. For these ensembles, the empirical distribution of a normalized deflation time is found to collapse onto a curve that depends only on the algorithm, but not on the matrix size or deflation tolerance provided the matrix size is large enough (see Figure 4, Figure 7 and Figure 10). For the QR algorithm with the Wilkinson shift, the observed universality is even stronger and includes certain non-Wigner ensembles. Our experiments also provide a quantitative statistical picture of the accelerated convergence with shifts.