Разложение Шура
T = schur(A)
T = schur(A,flag)
[U,T] = schur(A,...)
schur функция вычисляет форму Шура матрицы.
T = schur(A) возвращает матрицу Шура T.
T = schur(A,flag) для вещественной матрицы A возвращает матрицу Шура T в одной из двух форм в зависимости от значения flag:
'complex' |
|
'real' |
|
Если A является сложным, schur возвращает комплексную форму Шура в матрице T и flag игнорируется. Комплексная форма Шура является верхней треугольной с собственными значениями A по диагонали.
Функция rsf2csf преобразует действительную форму Шура в сложную форму Шура.
[U,T] = schur(A,...) также возвращает унитарную матрицу U чтобы A = U*T*U' и U'*U = eye(size(A)).
H является матрицей теста 3 на 3 собственных значений:
H = [ -149 -50 -154
537 180 546
-27 -9 -25 ]Его форма Schur
schur(H)
ans =
1.0000 -7.1119 -815.8706
0 2.0000 -55.0236
0 0 3.0000Собственные значения, которые в данном случае 1, 2, и 3, находятся на диагонали. Тот факт, что внедиагональные элементы настолько велики, указывает на то, что эта матрица имеет плохо обусловленные собственные значения; небольшие изменения в элементах матрицы приводят к относительно большим изменениям ее собственных значений.