Разложение Шура

Материал из testwiki
Перейти к навигации Перейти к поиску

Разложение Шура — разложение матрицы на унитарную, верхнюю треугольную и обратную унитарную матрицы, названное именем Исая Шура.

Утверждение

Если A является квадратной матрицей порядка n с комплексными элементами, то её можно представить в виде[1][2]:

A=QUQ1

где Q — унитарная матрица (так что её обратная Q1 является эрмитово-сопряжённой Q* матрицы Q), а U — верхняя треугольная матрица, которая называется формой Шура матрицы A. Поскольку U подобна матрице A, она имеет то же мультимножество собственных значений, а поскольку она треугольна, эти собственные значения совпадают с диагональными элементами матрицы U.

Из разложения Шура следует, что существует вложенная последовательность A-инвариантных подпространств {0}=V0V1Vn=n и упорядоченный ортогональный базис, такие что линейная комбинация первых i базисных векторов даёт Vi для всех i в последовательности. Иными словами, первая часть говорит, что линейное отображение J на комплексном конечномерном векторном пространстве стабилизирует весь флаг V1,Vn.

Доказательство

Конструктивное доказательство разложения Шура следующее: любой оператор A в комплексном конечномерном векторном пространстве имеет собственное значение λ, соответствующее собственному пространству Vλ. Пусть Vλ — ортонормальное дополнение. При таком ортогональном разложении A имеет матричное представление (можно выбрать любые ортонормальные базисы 1 и 2 для натянутых на них пространств Vλ и Vλ соответственно):

[Z1Z2]*A[Z1Z2]=[λIλA120A22]:VλVλVλVλ,

где Iλ — тождественный оператор на Vλ. Полученная матрица треугольна за исключением блока A22. Но точно ту же процедуру можно совершить для подматрицы A22, которая рассматривается как оператор на Vλ и её подматрицы. Продолжив процедуру n раз, пространство n будет исчерпано и построение даст желаемый результат.

Особенности

Хотя любая квадратная матрица имеет разложение Шура, в общем случае такое разложение не единственно. Например, собственное пространство Vλ может иметь размерность более 1, и в этом случае любой ортонормальный базис для Vλ даст желаемый результат.

Треугольная матрица U может быть представлена в виде суммы диагональной D и строго верхней треугольной N: U=D+N. Строго верхняя треугольная матрица нильпотентна. Диагональная матрица D содержит собственные значения матрицы A в случайном порядке. Нильпотентная часть N в общем случае также не уникальна, но её норма Фробениуса единственным образом определяется матрицей A, так как норма Фробениуса матрицы A равна норме Фробениуса матрицы U=D+N.

Если A является нормальной, то её форма Шура U диагональна, а столбцы матрицы Q разложения QUQ1 будут собственными векторами матрицы A. Таким образом, разложение Шура обобщает спектральное разложение. В частности, если A является положительно определённой, её разложение Шура, её спектральное разложение и её сингулярное разложение совпадают.

Коммутативное семейство матриц {Ai} может быть приведено к треугольному виду одновременно, то есть существует унитарная матрица Q, такая что для любой Ai из данного семейства выполнено QAiQ* является верхней треугольной. Конечное утверждение доказывается индукцией. Как следствие, любое коммутативное семейство нормальных матриц может быть приведено к диагональному виду[3].

В бесконечномерном случае не всякий ограниченный оператор в банаховом пространстве имеет инвариантное подпространство. Однако приведение к треугольному виду произвольной квадратной матрицы обобщается для компактных операторов. Любой компактный оператор в банаховом пространстве имеет гнездо замкнутых инвариантных подпространств.

Вычисление

Декомпозиция Шура заданной матрицы выполняется QR-алгоритмом или его вариантами. С использованием таких алгоритмов для разложения Шура нет необходимости заранее вычислять корни характеристического многочлена, соответствующего матрице. И наоборот, QR-алгоритм можно использовать для вычисления корней любого заданного характеристического многочлена путём нахождения разложения Шура его сопровождающей матрицы. Таким же образом QR-алгоритм используется для вычисления собственных значений любой заданной матрицы, которые являются диагональными элементами верхней треугольной матрицы разложения Шура. Все необходимые алгоритмы реализованы, в частности, в библиотеке Lapack[4].

Приложения

Из разложения Шура следуют некоторые важные результаты Шаблон:Не переведено 5, в частности:

Обобщённое разложение Шура

Обобщённое разложение Шура двух квадратных матриц A и B — согласованная пара разложений обеих матриц A=QSZ* и B=QTZ*, где Q и Z — унитарны, а S и T — треугольные. Обобщённое разложение Шура иногда называется также QZ-разложением.

Обобщённые собственные значения λ, решающие задачу обобщённых значений Ax=λBx (где x — неизвестный ненулевой вектор), могут быть вычислены как отношение диагональных элементов S к соответствующит элементам T. То есть, i-е обобщённое собственное значение λi удовлетворяет равенству λi=Sii/Tii.

Примечания

Шаблон:Примечания

Литература

Шаблон:Rq