Tengo una matriz cuadrada simétrica real densa. La dimensión es de aproximadamente 1000x1000. Necesito calcular el primer componente principal y preguntarme cuál sería el mejor algoritmo para hacerlo.
Parece que MATLAB usa los algoritmos Arnoldi / Lanczos (para eigs
). Pero al leer sobre ellos, no estoy seguro de si tienen alguna ventaja sobre la iteración de potencia simple , ya que mi matriz no es escasa y solo estoy interesado en el primer vector propio.
¿Alguna recomendación de cuál es el algoritmo más rápido en este caso?