La matriz puede ser tan grande como $2500\times 2500$, ¿cuál es el mejor algoritmo para hacer eso, ¿existe algún algoritmo que es fácil escribir un programa, ¿hay algún paquete conveniente para?
Respuestas
¿Demasiados anuncios?
KP.
Puntos
1177
Echa un vistazo en A encuesta de Software para problemas de valor propio escaso por Hernández et al.
Mohit Jain
Puntos
412
No sé mucho sobre valores propios o lo que son aplicables a y R parece que disponen de una función para ello llamado eigen()
. Calcular los autovalores y autovectores de una matriz de 2500 * 2500 tomó ~ 1 minuto en mi máquina.
> sampData <- runif(6250000, 0, 2)
> x <- matrix(sampData, ncol = 2500, byrow = TRUE)
> system.time(eigen(x))
user system elapsed
79.74 2.90 65.69
Esta pregunta ha surgido también en Stack Overflow.