Найдём собственные значения матриц
и
при помощи какой-нибудь программы, например, REDUCE.
Для этого определим матрицу
коммандой:
m3:=mat((0, 0, 0, 1, 0), (0, 0, 1, 0, 0), (0, 1, 0, 0, 0), (1, 0, 0, 0, 0), (0, 0, 0, 0, 4));
Определим единичную матрицу коммандой:
i5:=mat((1, 0, 0, 0, 0), (0, 1, 0, 0, 0), (0, 0, 1, 0, 0), (0, 0, 0, 1, 0), (0, 0, 0, 0, 1));
Вычислим характеристический полином матрицы
коммандой:
det(m4-i5 * x);
Получим:
Найдём собственные значения и собственные векторы матрицы
коммандой:
mateiden(m3, x);
Получим:
, то есть
,
и
являются собственными значениями матрицы
.
Проверяем коммандой:
(x-4) * (x-1)^2 * (x+1)^2;
Получим:
Комманда mateiden(m3, x); даёт нам также собственные векторы:
,
,
,
,
.
Собственные векторы, соответствующие разным собственным значениям ортогональны.
Это не случайно - есть такая теорема.
Нам повезло: собственные векторы одного и того же собственного значения ортогональны, так что нет необходимости задействовать алгоритм Грамм-Шмидта для нахождения ортогонального базиса.
Образуем матрицу из этих векторов коммандой:
u:=mat((0, 0, 0, 0, 1), (1, 0, 0, 1, 0), (0, 1, 1, 0, 0), (-1, 0, 0, 1, 0), (0, -1, 1, 0, 0));
Продолжение следует.