Вращений Метод

метод Якоби,-метод решения полной проблемы собственных значений эрмитовой матрицы, основанный на подобном преобразовании эрмитовой матрицы к диагональному виду с помощью последовательности плоских вращений. В. М. - итерационный метод, он имеет простую вычислительную схему и всегда сходится, причем скорость сходимости асимптотически квадратичная. Наличие кратных и близких собственных значений у матрицы не вызывает осложнений. В. М. Позволяет вычислить собственные значения как с нахождением собственных векторов, так и без них. Система собственных векторов, вычисленная по В. М., ортонормирована. Идеи В. М. Предложены в [1]. В современном виде это один из наиболее развитых и эффективных по реализации на ЭВМ методов решения полной проблемы собственных значений матрицы.

Классический В. М. Состоит в построении последовательности матриц - исходная матрица, - матрица плоского вращения, аннулирующего максимальный по модулю внедиагональный элемент матрицы . При этом, если то матрица отличается от единичной лишь элементами В действительном случае, когда А - симметрическая матрица, (*) В комплексном случае соотношения (*) незначительно усложняются. Последовательность матриц сходится к диагональной матрице , скорость сходимости асимптотически квадратичная. Диагональные элементы матрицы являются приближенными собственными значениями А, а столбцы матрицы - приближенными собственными векторами. Реализация описанного варианта В. М. Требует выбора максимального по модулю внедиагонального элемента матрицы на каждом шаге.

Для выполнения этой операции на ЭВМ требуется значительный объем вычислительной работы. Существуют другие варианты В. М., более эффективные в этом отношении. Циклич. В. М., В. М. С барьером, В. М. С выбором оптимального элемента. В циклическом В. М. Пары индексов аннулирующего элемента пробегают циклически все над-диагональные позиции. Недостатком этого процесса является возможность выполнения большого числа неэффективных вращений, аннулирующих малые внедиа-гональные элементы. Этот недостаток частично устраняется в барьерном В. М., в к-ром вводится монотонно убывающая к нулю последовательность чисел называемых барьерами, и при циклич. Просмотре индексов (р, q). Аннулируются лишь те внедиагональные элементы, к-рые по модулю меньше .

После того, как Все внедиагональные элементы стали меньше , барьер заменяется на барьер и процесс продолжается. Практич. Использование этого варианта В. М. Встречает ряд трудностей, связанных с выбором оптимального барьера. Наиболее эффективным для использования в вычислительной практике является В. М. С выбором оптимального элемента [4]. В этом методе пары индексов ( р, q).соответствуют почти максимальному элементу и выбираются так, что р - номер строки с максимальной евклидовой длиной, q - номер столбца максимального по модулю внедиагонального элемента р-й строки. Поскольку на каждом шаге процесса строки матрицы, за исключением р-й и q-й, не меняют длины, то выбор индексов ( р, q).существенно не увеличивает общего объема вычислений.

Вся теория классич. В. М. Полностью переносится на описанную модификацию [2]. Расчетные формулы В. М., реализующие вычисление (*), обеспечивают сходимость процесса В. М. В реальных условиях машинной арифметики и высокую точность как собственных значений, так и собственных векторов [5]. Лит..- [1] Jacobi С. G. J., "J. Reine imd anprew. Math.", 1846, Bd 30, S. 51-94. [2] Воеводин В. В., Численные методы алгебры, М., 1966. [3] Уилкинсон Д ж. X., Алгебраическая проблема собственных значений, пер. С англ., М., 1970. [4] Воеводин В. В., Ким Г. Д., Программа для нахождения собственных значений и собственных векторов симметрической матрицы методом вращений, в сб. Вычислительные методы и программирование, М., 1962, с. 269-77. [5] Ким Г. Д., в сб. Численный анализ на ФОРТРАН'е, в.

3, М., 1973, с. 97-113. Г. Д. Ким.