Предположим, что - симметричная положительно определенная матрица. достаточно большой, чтобы решить дорого .
Существует ли итерационный алгоритм для наименьшего собственного значения , который не включает инвертирование в каждой итерации?
То есть мне нужно было бы использовать итерационный алгоритм, такой как сопряженные градиенты, для решения , поэтому многократное применение кажется дорогим «внутренним циклом». Мне нужен только один собственный вектор.
Благодаря!
linear-algebra
eigensystem
eigenvalues
iterative-method
Джастин Соломон
источник
источник
eigs
-routine. Это итеративный метод. Есть варианты, чтобы указать, какое собственное значение вы хотите, например, наименьшее действительное .Ответы:
Вычислите собственное значение наибольшей величины из A (скажем, с ).λmax A
eigs('lm')
Затем вычислить величину наибольшей (отрицательный) собственное значение λ т а х из М = А - λ м х I (опять же , с помощью стандартного вызова ).λ^max M=A−λmaxI
eigs('lm')
Заметим , что λ т х + λ макс = λ м I п ( А ) . Причина, по которой это верно, объясняется здесь .λ^max+λmax=λmin(A)
Найдите свой собственный вектор , решив .( A - λ m i n I ) v = 0v (A−λminI)v=0
источник