Abstract
Multiple Perron eigenvectors of non-negative matrices occur in applications, where they often become a source of trouble. A usual way to avoid it and to make the Perron eigenvector simple is a regularization of matrix: an initial non-negative matrix A is replaced by \(A + \varepsilon M\), where M is a strictly positive matrix and \(\varepsilon > 0\) is small. However, this operation is numerically unstable and may lead to a significant increase of the Perron eigenvalue, especially in high dimensions. We define a selected Perron eigenvector of A as the limit of normalized Perron eigenvectors of the regularizations \(A + \varepsilon M\) as \(\varepsilon \rightarrow 0\). It is shown that if the matrix M is rank-one, then the limit eigenvector can be found by an explicit formula and, moreover, is efficiently computed by the power method. The role of the rank-one condition is explained.
Similar content being viewed by others
References
Anderson, J.: Distance to the nearest stable Metzler matrix. In: IEEE 56th Annual Conference on Decision and Control (CDC), pp. 6567–6572 (2017)
Berman, A., Plemmons, R.J.: Nonnegative Matrices in the Mathematical Sciences. Acedemic Press, Now York (1979)
Byers, R.: A bisection method for measuring the distance of a stable to unstable matrices. SIAM J. Sci. Stat. Comput. 9, 875–881 (1988)
Cicone, A., Serra-Capizzano, S.: Google pageranking problem: the model and the analysis. J. Comput. Appl. Math. 234(11), 3140–3169 (2010)
Cvetkovic, A., Protasov, V.: The greedy strategy in optimizing the Perron eigenvalue. submitted to Math. Progr. arXiv:1807.05099
Gantmacher, F.R.: The Theory of Matrices. Chelsea, New York (2013)
Guglielmi, N., Manetta, M.: Approximating real stability radii. IMA J. Numer. Anal. 35(3), 1402–1425 (2015)
Guglielmi, N., Protasov, V.: On the closest stable/unstable nonnegative matrix and related stability radii. SIAM J. Matrix Anal. 39(4), 1642–1669 (2018)
Horn, R.A., Johnson, C.R.: Matrix Analysis. Cambridge University Press, Cambridge (1990)
Kato, T.: Perturbation Theory for Linear Operators, Classics in Mathematics. Springer, Berlin (2013)
Krause, U.: Swarm dynamics and positive dynamical systems. In: Differential and Difference Equations with Applications. Springer Proceedings in Mathematics and Statistics, vol. 47, pp. 69–81 (2013)
Nesterov, Y., Nemirovski, A.: Finding the stationary states of Markov chains by iterative methods. Appl. Math. Comput. 255, 58–65 (2015)
Nesterov, Y., Protasov, V.Y.: Optimizing the spectral radius. SIAM J. Matrix Anal. Appl. 34(3), 999–1013 (2013)
Protasov, V.Y.: The spectral simplex method. Math. Prog. 156, 485–511 (2016)
Serra-Capizzano, S.: Jordan Canonical form of the Google matrix: a potential contribution to the pagerank computation. SIAM J. Matrix Anal. Appl. 27(2), 305–312 (2006)
Stewart, G.W., Sun, J.G.: Matrix Perturbation Theory. Academic Press, New York (1990)
Valcher, M.E., Misra, P.: On the stabilizability and consensus of positive homogeneous multi-agent dynamical systems. IEEE Trans. Autom. Control 59, 1936–1941 (2014)
Acknowledgements
The author is grateful to the anonymous Referees for attentive reading and for many valuable comments and suggestions. The article was prepared within the framework of the HSE University Basic Research Program and funded by the Russian Academic Excellence Project ‘5-100’.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Scientific results of Sections 3 and 4 of this paper were obtained with support of RSF Grant 17-11-01927.
Rights and permissions
About this article
Cite this article
Protasov, V.Y. How to make the Perron eigenvector simple. Calcolo 56, 17 (2019). https://doi.org/10.1007/s10092-019-0314-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10092-019-0314-7
Keywords
- Non-negative matrix
- Spectral radius
- Perron eigenvector
- Regularization
- Uniqueness
- Geometric multiplicity
- Power method