Цитата:
Решение несимметричной задачи собственных значений осуществляется в несколько этапов. На первом этапе матрица приводится ортогональным преобразованием к верхней форме Хессенберга. На втором этапе, занимающем больше всего времени, матрица приводится ортогональным преобразованием к верхней форме Шура. Если требуются только собственные значения, то этого достаточно, т.к. собственные числа матрицы располагаются в диагональных блоках квазитреугольной матрицы из каконической формы Шура. Если же требуются собственные векторы, то они могут быть получены по векторам Шура и квазитреугольной матрице путем обратной подстановки (фактически - решения системы линейных уравнений; сам процесс обратной подстановки занимает незначительную часть времени работы алгоритма, но необходимость накапливать проводимые над матрицей преобразования для применения их к собственным векторам замедляет алгоритм более чем в два раза).
|
Информация взята
отсюда, там есть и примеры реализации алгоритма на разных языках программирования