Tengo un vector que puede descomponerse en el espacio propio del operador escaso hermitiano :
¿Hay una manera de encontrar el m i (el vector propio sí mismo) que se corresponden con la mayor v i (en magnitud)?
Básicamente, quiero los pocos términos más grandes de la suma, incluidos los vectores propios de , que no conozco de antemano.
Específicamente, quiero encontrar simultáneamente los vectores propios de que corresponden a los más grandes , junto con encontrar el mayor . Preferiblemente sin encontrar el espectro completo de primero.
Algunas posibilidades en las que he estado pensando:
Podemos "inflar" la matriz usando el opuesto de "Deflación de Wieldant":
Los valores propios de diferentes m i se desplazan λ i + sigma | v i | 2 . Creo que podemos extraer σ y v i porque los vectores propios no cambian. El problema es que el producto externo de V es denso.
Otra posibilidad:
¿Hay alguna forma de controlar esto para que solo converjamos en el componente más grande?
fuente
Respuestas:
Como la matriz es hermitiana, puede usarla como hamiltoniana para propagarla en tiempo imaginario. Es decir, resuelva el siguiente sistema de ecuaciones diferenciales:
La solución general a esto es:
Una vez que tenga los valores propios, encuentre los vectores propios con cualquier solucionador de valores propios específico que prefiera.
fuente
fuente