Если матрицы горбатые или ленточные можно попробовать простой метод Якоби, или для комплекнозначных матриц метод Якоби-Эберляйн (Справочник алгоритмов на языке Алгол, Уилкинсон, Райнш). Дело в том, что внедиагональные элементы очень быстро убывают, причем кажется, что даже экспоненциально, и спектр определяется достаточно хорошо после нескольких больших итераций. Изобрел это С.К. Годунов в 50-ые, им с товарищем пришлось решать задачу переноса нейтронов, и быстро-быстро на машине Стрела определять весь спектр.
|