Решение собственных векторов / собственных значений разреженной матрицы в разных диапазонах последовательно - PullRequest
2 голосов
/ 24 апреля 2019

Я ищу собственные значения / собственные векторы огромной разреженной матрицы A (250k * 250k), которые лежат в разных диапазонах собственных значений без вычисления всех собственных значений между , поскольку просто займите слишком много времени.


Мое рассмотрение было

Шаг 1: Я пытаюсь найти k определенных собственных значений / собственных векторов матрицы A, близких к определенному значению sigma1. Для этого я использую [V1,D1] = eigs(A,k,sigma1). Пока это не имеет большого значения.

Step2: Я делаю то же самое для sigma2 [V2,D2] = eigs(A,k,sigma2).

К сожалению, собственные векторы V2 шага 2 не являются отогональными по отношению к собственным векторам V1 шага 1. Похоже, они не имеют той же основы, что и собственные векторы V1. Есть ли способ заставить Matlab использовать тот же базис или каким-либо образом передать команду eigs собственными векторами, найденными в шаге 1?


Использование одного и того же начального вектора для Step1 и Step2, к сожалению, не помогло. Заранее благодарю за помощь!

...