Джордан в нормальной форме

Пример матрицы в жордановой нормальной форме. Все не показанные элементы матрицы равны нулю. Очерченные квадраты известны как «иорданские блоки». Каждый блок Джордана содержит одно число лямбда на главной диагонали и одно число над главной диагональю. Лямбды — это собственные значения матрицы; они не обязательно должны быть различимы.

В линейной алгебре — йордановая нормальная форма , также известная как йордановая каноническая форма ( JCF ), — верхнетреугольная матрица определенного вида, называемая йордановой матрицей, представляющая линейный оператор в конечномерном векторном пространстве относительно некоторого базиса . Такая матрица имеет каждый ненулевой внедиагональный элемент, равный 1, непосредственно над главной диагональю (на супердиагонали ) и с одинаковыми диагональными элементами слева и под ними.

Пусть V векторное пространство над полем K. — Тогда базис, относительно которого матрица имеет требуемый вид, существует тогда и только тогда, когда все собственные значения матрицы лежат в K или, что то же самое, если характеристический полином оператора распадается на линейные множители над K . Это условие всегда выполняется, если К ( алгебраически замкнуто например, если оно является полем комплексных чисел ). Диагональные элементы нормальной формы являются собственными значениями (оператора), а количество раз, которое встречается каждое собственное значение, называется алгебраической кратностью собственного значения.

Если оператор изначально задан квадратной матрицей M , то его жорданова нормальная форма также называется жордановой нормальной M. формой Любая квадратная матрица имеет жордановую нормальную форму, если поле коэффициентов расширить до поля, содержащего все собственные значения матрицы. Несмотря на свое название, нормальная форма для данного M не совсем уникальна, поскольку это блочная диагональная матрица, образованная из жордановых блоков , порядок которых не фиксирован; принято группировать блоки для одного и того же собственного значения вместе, но никакой порядок не налагается ни между собственными значениями, ни между блоками для данного собственного значения, хотя последние можно, например, упорядочить путем слабо уменьшающегося размера.

Разложение Жордана –Шевалле особенно просто в отношении базиса, для которого оператор принимает свою жордановую нормальную форму. Диагональная форма для диагонализируемых матриц, например нормальных матриц , является частным случаем жордановой нормальной формы.

Нормальная форма Жордана названа в честь Камиля Джордана , который впервые сформулировал теорему Жордана о разложении в 1870 году.

Обзор

Обозначения

В некоторых учебниках есть те, что расположены на субдиагонали ; то есть сразу под главной диагональю, а не на супердиагонали. Собственные значения все еще находятся на главной диагонали.

Мотивация

Матрица n × n размера A диагонализуема тогда и только тогда, когда сумма размерностей собственных пространств равна n . Или, что то же самое, тогда и только тогда, когда A имеет n линейно независимых собственных векторов . Не все матрицы диагонализуемы; Матрицы, недиагонализуемые, называются дефектными . Рассмотрим следующую матрицу:

$A=\left[{\begin{array}{*{20}{r}}5&4&2&1\\[2pt]0&1&-1&-1\\[2pt]-1&-1&3&0\\[2pt]1&1&-1&2\end{array}}\right].$

С учетом кратности собственные значения A равны λ = 1, 2, 4, 4. Размерность собственного пространства, соответствующего собственному значению 4, равна 1 (а не 2), поэтому A не является диагонализируемым. Однако существует обратимая матрица P такая, что J = P ⁻¹ АП , где

$J={\begin{bmatrix}1&0&0&0\\[2pt]0&2&0&0\\[2pt]0&0&4&1\\[2pt]0&0&0&4\end{bmatrix}}.$

Матрица $J$ почти диагональная. Это жордановая нормальная форма A . В разделе «Пример» ниже приведены подробности вычислений.

Комплексные матрицы

В общем, квадратная комплексная матрица A похожа . на блочную диагональную матрицу

$J={\begin{bmatrix}J_{1}&\;&\;\\\;&\ddots &\;\\\;&\;&J_{p}\end{bmatrix}}$

где каждый блок J _i представляет собой квадратную матрицу вида

$J_{i}={\begin{bmatrix}\lambda _{i}&1&\;&\;\\\;&\lambda _{i}&\ddots &\;\\\;&\;&\ddots &1\\\;&\;&\;&\lambda _{i}\end{bmatrix}}.$

Итак, существует обратимая матрица P такая, что P ⁻¹ AP = J таков, что единственные ненулевые элементы J находятся на диагонали и супердиагонали. J называется нормальной A. формой жордановой Каждый J _i называется блоком A . жордановым В данном жордановом блоке каждая запись на супердиагонали равна 1.

Предполагая этот результат, мы можем вывести следующие свойства:

• При подсчете кратностей собственные значения J и, следовательно , A являются диагональными элементами.
• Учитывая собственное значение λ _i , его геометрическая кратность равна размерности ker( A − λ _i I ), где I - единичная матрица , и это количество жордановых блоков, соответствующих λ _i .
• Сумма размеров всех жордановых блоков, соответствующих собственному значению λi _, есть его алгебраическая кратность .
• A диагонализуема тогда и только тогда, когда для каждого собственного значения λ оператора A его геометрическая и алгебраическая кратности совпадают. В частности, жордановые блоки в этом случае представляют собой матрицы размера 1 × 1; то есть скаляры.
• Джордановый блок, соответствующий λ, имеет вид λI + N , где N — нильпотентная матрица, определяемая как N _ij = δ _{i , j −1} (где δ — дельта Кронекера ). Нильпотентность N можно использовать при вычислении f ( A ), где f — комплексная аналитическая функция. Например, в принципе жордановая форма может дать выражение в замкнутой форме для экспоненты exp( A ).
• Число жордановых блоков, соответствующих λ _i размера не менее j, равно dim ker( A − λ _i I ) ^дж - dim ker( А - λ _я я ) ^{j −1} . Таким образом, количество жордановых блоков размера j равно

  $2\dim \ker(A-\lambda _{i}I)^{j}-\dim \ker(A-\lambda _{i}I)^{j+1}-\dim \ker(A-\lambda _{i}I)^{j-1}$

• Учитывая собственное значение λ _i , его кратность в минимальном многочлене равна размеру его наибольшего жорданового блока.

Пример

Рассмотрим матрицу $A$ из примера в предыдущем разделе. Жорданова нормальная форма получается некоторым преобразованием подобия :

$P^{-1}AP=J;$ то есть, $AP=PJ.$

Позволять $P$ иметь векторы-столбцы $p_{i}$, $i=1,\ldots ,4$, затем

$A{\begin{bmatrix}p_{1}&p_{2}&p_{3}&p_{4}\end{bmatrix}}={\begin{bmatrix}p_{1}&p_{2}&p_{3}&p_{4}\end{bmatrix}}{\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}={\begin{bmatrix}p_{1}&2p_{2}&4p_{3}&p_{3}+4p_{4}\end{bmatrix}}.$

Мы видим это

$(A-1I)p_{1}=0$

$(A-2I)p_{2}=0$

$(A-4I)p_{3}=0$

$(A-4I)p_{4}=p_{3}.$

Для $i=1,2,3$ у нас есть $p_{i}\in \ker(A-\lambda _{i}I)$, то есть, $p_{i}$ является собственным вектором $A$ соответствующее собственному значению $\lambda _{i}$. Для $i=4$, умножив обе части на $(A-4I)$ дает

$(A-4I)^{2}p_{4}=(A-4I)p_{3}.$

Но $(A-4I)p_{3}=0$, так

$(A-4I)^{2}p_{4}=0.$

Таким образом, $p_{4}\in \ker(A-4I)^{2}.$

Векторы, такие как $p_{4}$ называются обобщенными собственными векторами оператора A .

Пример: Получение нормальной формы

В этом примере показано, как вычислить жорданову нормальную форму заданной матрицы.

Рассмотрим матрицу

$A=\left[{\begin{array}{rrrr}5&4&2&1\\0&1&-1&-1\\-1&-1&3&0\\1&1&-1&2\end{array}}\right]$

о котором сказано в начале статьи.

Характеристический полином A ​ равен

${\begin{aligned}\chi (\lambda )&=\det(\lambda I-A)\\&=\lambda ^{4}-11\lambda ^{3}+42\lambda ^{2}-64\lambda +32\\&=(\lambda -1)(\lambda -2)(\lambda -4)^{2}.\,\end{aligned}}$

Это показывает, что собственные значения равны 1, 2, 4 и 4 в соответствии с алгебраической кратностью. Собственное пространство, соответствующее собственному значению 1, можно найти, решив уравнение Av = λv . Он натянут на вектор-столбец v = (−1, 1, 0, 0) ^Т . Аналогично, собственное пространство, соответствующее собственному значению 2, охватывается w = (1, −1, 0, 1). ^Т . Наконец, собственное пространство, соответствующее собственному значению 4, также является одномерным (хотя это двойное собственное значение) и натянуто на x = (1, 0, −1, 1). ^Т . Итак, геометрическая кратность (т. е. размерность собственного пространства данного собственного значения) каждого из трёх собственных значений равна единице. Следовательно, два собственных значения, равные 4, соответствуют одному жордановому блоку, а жордановая нормальная форма матрицы A представляет собой прямую сумму

$J=J_{1}(1)\oplus J_{1}(2)\oplus J_{2}(4)={\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}.$

Есть три сети Jordan . Два имеют длину один: { v } и { w }, что соответствует собственным значениям 1 и 2 соответственно. Существует одна цепочка длины два, соответствующая собственному значению 4. Чтобы найти эту цепочку, вычислите

$\ker(A-4I)^{2}=\operatorname {span} \,\left\{{\begin{bmatrix}1\\0\\0\\0\end{bmatrix}},\left[{\begin{array}{r}1\\0\\-1\\1\end{array}}\right]\right\}$

где I — единичная матрица 4 × 4. Выберите вектор в указанном выше диапазоне, который не входит в ядро ​​A − 4 I ; например, y = (1,0,0,0) ^Т . Теперь ( A − 4 I ) y = x и ( A − 4 I ) x = 0, поэтому { y , x } — это цепочка длины два, соответствующая собственному значению 4.

Матрица перехода P такая, что P ⁻¹ AP = J формируется путем размещения этих векторов рядом друг с другом следующим образом

$P=\left[{\begin{array}{c|c|c|c}v&w&x&y\end{array}}\right]=\left[{\begin{array}{rrrr}-1&1&1&1\\1&-1&0&0\\0&0&-1&0\\0&1&1&0\end{array}}\right].$

Расчет показывает, что уравнение P ⁻¹ AP = J действительно имеет место.

$P^{-1}AP=J={\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}.$

Если бы мы поменяли порядок появления цепных векторов, то есть изменили порядок v , w и { x , y } вместе, жордановые блоки поменялись бы местами. Однако жордановые формы являются эквивалентными жордановыми формами.

Обобщенные

Основная статья: Обобщенный собственный вектор

Учитывая собственное значение λ , каждый соответствующий жордановый блок порождает жордановую цепочку линейно независимых векторов p _i , i = 1 , ..., b , где b — размер жорданового блока. Генератор ( или ведущий вектор p b _цепи является обобщенным собственным вектором таким, что A − λ I ) ^б p _b = 0. Вектор p ₁ = ( A − λ I ) ^{б -1} p _b — обычный собственный вектор, соответствующий λ . В общем, p _i является прообразом p _{i −1} при A − λ I . Таким образом, ведущий вектор генерирует цепочку путем умножения на A − λ I . Следовательно, утверждение о том, что каждую квадратную матрицу A можно привести к жордановой нормальной форме, эквивалентно утверждению, что базис основного векторного пространства состоит из жордановых цепей.

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

Мы докажем по индукции , что любую комплексную квадратную матрицу A можно привести к жордановой нормальной форме. Поскольку базовое векторное пространство можно показать Будучи прямой суммой инвариантных подпространств, связанных с собственными значениями, A можно предположить, что имеет только одно собственное значение λ . Случай 1 × 1 тривиален. Пусть A — матрица размера n × n . Диапазон − A , − λ I обозначаемый Ran( A . λ I является инвариантным подпространством A ) , Кроме того, поскольку λ является собственным значением A , размерность Ran( A − λ I ), r , строго меньше n , поэтому по индуктивному предположению Ran ( A − λ I ) имеет базис { p ₁ , ..., p _r } составленная из жордановых цепей.

Далее рассмотрим ядро ​​, то есть подпространство ker( A − λ I ). Если

$\operatorname {Ran} (A-\lambda I)\cap \ker(A-\lambda I)=\{0\},$

искомый результат непосредственно следует из теоремы о ранге-нулевости . (Так было бы, например, если бы А было эрмитовым .)

В противном случае, если

$Q=\operatorname {Ran} (A-\lambda I)\cap \ker(A-\lambda I)\neq \{0\},$

пусть размерность Q равна s ≤ r . Каждый вектор из Q является собственным вектором, поэтому Ran( A − λ I ) должен содержать s жордановых цепей, соответствующих s линейно независимым собственным векторам. Поэтому базис { p ₁ , ..., p _r } должен содержать s векторов, скажем { p ₁ , ..., p _s }, которые являются ведущими векторами этих жордановых цепей. Мы можем «продлить цепочки», взяв прообразы этих ведущих векторов. (Это ключевой шаг.) Пусть q _i таково, что

$\;(A-\lambda I)q_{i}=p_{i}{\mbox{ for }}i=1,\ldots ,s.$

Наконец, мы можем выбрать любую основу для

$\ker(A-\lambda I)/Q$

а затем поднимемся до векторов { z ₁ , ..., z _t } в ker( A − λI ). Каждый z _i образует жорданову цепь длины 1. Нам просто нужно показать, что объединение { p ₁ , ..., p _r }, { z ₁ , ..., z _t } и { q ₁ , . .., q _s } образует основу векторного пространства.

По теореме о нулевом ранге dim(ker( A − λ I ))= nr , поэтому t=nrs , и поэтому количество векторов в потенциальном базисе равно n. Чтобы показать линейную независимость, предположим, что некоторая линейная комбинация векторов равна 0. Применяя A − λ I, мы получаем некоторую линейную комбинацию pi _, где q _i становится ведущим вектором среди pi _. Из линейной независимости pi _{следует ,} что коэффициенты векторов q _i должны быть равны нулю. Более того, никакая нетривиальная линейная комбинация z i _не может равняться линейной комбинации p _i , потому что тогда она принадлежала бы Ran( A − λ I ) и, следовательно, Q, что невозможно по конструкции z _i . Следовательно, коэффициенты z i _также будут равны 0. Остается только p _i членов, которые считаются линейно независимыми, и поэтому эти коэффициенты также должны быть равны нулю. Мы нашли базис, состоящий из жордановых цепей, и это показывает, что A можно привести к жордановой нормальной форме.

Уникальность

Можно показать, что жордановая нормальная форма данной матрицы A единственна с точностью до порядка жордановых блоков.

Знания алгебраической и геометрической кратностей собственных значений недостаточно для определения жордановой нормальной формы A . Предполагая, что алгебраическая кратность m ( λ ) собственного значения λ известна, структуру жордановой формы можно выяснить, анализируя ранги степеней ( A − λI ) ^{м ( λ )} . Чтобы убедиться в этом, предположим, что размера n × n матрица A имеет только одно собственное значение λ . Итак м ( λ ) знак равно п . Наименьшее целое число k ₁ такое, что

$(A-\lambda I)^{k_{1}}=0$

размер самого большого жорданового блока в жордановой форме A. — (Это число k ₁ также называется индексом λ . См . обсуждение в следующем разделе.) Ранг

$(A-\lambda I)^{k_{1}-1}$

— количество жордановых блоков размера k ₁ . Аналогично, ранг

$(A-\lambda I)^{k_{1}-2}$

в два раза больше количества жордановых блоков размера k ₁ плюс количество жордановых блоков размера k ₁ − 1. Общий случай аналогичен.

Это можно использовать, чтобы показать уникальность формы Джордана. Пусть J ₁ и J ₂ — две жордановые нормальные формы A . Тогда J ₁ и J ₂ подобны и имеют одинаковый спектр, включая алгебраические кратности собственных значений. Для определения структуры этих матриц можно использовать процедуру, описанную в предыдущем параграфе. Поскольку ранг матрицы сохраняется преобразованием подобия, существует биекция между жордановыми блоками J ₁ и J ₂ . Это доказывает единственность утверждения.

Реальные матрицы

Если A — действительная матрица, ее жорданова форма все равно может быть невещественной. Вместо представления его комплексными собственными значениями и значениями на супердиагонали, как обсуждалось выше, существует действительная обратимая матрица P такая, что P ⁻¹ AP = J — вещественная блочная диагональная матрица , каждый блок которой является действительным жордановым блоком. Действительный жордановый блок либо идентичен комплексному жордановскому блоку (если соответствующее собственное значение $\lambda _{i}$ является вещественным) или представляет собой саму блочную матрицу, состоящую из блоков 2×2 (для недействительных собственных значений $\lambda _{i}=a_{i}+ib_{i}$ с заданной алгебраической кратностью) вида

$C_{i}=\left[{\begin{array}{rr}a_{i}&-b_{i}\\b_{i}&a_{i}\\\end{array}}\right]$

и описать умножение на $\lambda _{i}$ в комплексной плоскости. Супердиагональные блоки представляют собой единичные матрицы размером 2 × 2, и, следовательно, в этом представлении размеры матрицы больше, чем комплексная жорданова форма. Полный реальный блок Джордана определяется выражением

$J_{i}={\begin{bmatrix}C_{i}&I&&\\&C_{i}&\ddots &\\&&\ddots &I\\&&&C_{i}\end{bmatrix}}.$

Эта действительная жорданова форма является следствием комплексной жордановой формы. Для вещественной матрицы всегда можно выбрать невещественные собственные векторы и обобщенные собственные векторы, образующие комплексно-сопряженные пары. Учитывая действительную и мнимую часть (линейную комбинацию вектора и сопряженного ему), матрица имеет такой вид по отношению к новому базису.

Матрицы с записями в поле

Жордановую редукцию можно распространить на любую квадратную матрицу M, которой лежат в поле K. элементы Результат утверждает, что любое можно записать в виде суммы D + N где D полупросто , , N нильпотентно M и DN = ND . Это называется разложением Жордана – Шевалле . Всякий раз, когда K содержит собственные значения M , в частности, когда K , алгебраически замкнут нормальная форма может быть явно выражена как прямая сумма жордановых блоков.

Аналогично случаю, когда K — комплексные числа, зная размерности ядер ( M − λI ) ^к для 1 ⩽ k ⩽ m , где — алгебраическая кратность собственного значения λ , позволяет определить жорданову форму M. m Мы можем рассматривать лежащее в основе векторное пространство V как K [ x ] -модуль , рассматривая действие x на V как применение M и расширяя его посредством K -линейности. Тогда полиномы ( x − λ ) ^к являются элементарными делителями M , а нормальная форма Жордана связана с представлением M в терминах блоков, связанных с элементарными делителями.

Доказательство жордановой нормальной формы обычно проводится как приложение к кольцу K [ x ] структурной теоремы для конечно порожденных модулей над областью главных идеалов , следствием которой оно является.

Последствия

Можно видеть, что жордановая нормальная форма по сути является результатом классификации квадратных матриц, и поэтому несколько важных результатов из линейной алгебры можно рассматривать как ее следствия.

Теорема

Используя нормальную форму Жордана, прямые вычисления дают теорему о спектральном отображении для полиномиального функционального исчисления : пусть A — матрица размера n × n с собственными значениями λ ₁ , ..., λ _n , тогда для любого многочлена p , p ( A ) имеет собственные значения p ( λ ₁ ), ..., p ( λ _n ).

Характеристический полином

Характеристический полином A ​ равен $p_{A}(\lambda )=\det(\lambda I-A)$. Подобные матрицы имеют одинаковый характеристический полином.Поэтому, ${\textstyle p_{A}(\lambda )=p_{J}(\lambda )=\prod _{i}(\lambda -\lambda _{i})^{m_{i}}}$,где $\lambda _{i}$ является i- м корнем ${\textstyle p_{J}}$ и $m_{i}$ является его кратностью, поскольку это, очевидно, характеристический многочлен жордановой формы A .

Теорема Кэли–Гамильтона

Теорема Кэли-Гамильтона что каждая матрица A удовлетворяет своему характеристическому уравнению: если p является характеристическим полиномом A утверждает , , то $p_{A}(A)=0$. Это можно показать прямым расчетом в жордановой форме, так как если $\lambda _{i}$ является собственным значением кратности $m$,тогда это блок Джордана $J_{i}$ явно удовлетворяет $(J_{i}-\lambda _{i}I)^{m_{i}}=0$.Поскольку диагональные блоки не влияют друг на друга, i -й диагональный блок из $(A-\lambda _{i}I)^{m_{i}}$ является $(J_{i}-\lambda _{i}I)^{m_{i}}=0$; следовательно ${\textstyle p_{A}(A)=\prod _{i}(A-\lambda _{i}I)^{m_{i}}=0}$.

расширяющим базовое поле матрицы, например, над полем расщепления p Можно предположить, что жордановая форма существует над полем , ; не меняет матрицу p ( A ) это расширение поля никак .

Минимальный полином

Минимальный многочлен P квадратной матрицы A — это уникальный монический многочлен наименьшей степени m , такой, что P ( A ) = 0. Альтернативно, набор многочленов, которые аннулируют данное A, образуют идеал I в C [ x ], область главных идеалов многочленов с комплексными коэффициентами. порождающий I, — это именно P. Монический элемент ,

Пусть λ ₁ , ..., λ _q — различные собственные значения A , а s _i — размер наибольшего жорданового блока, соответствующего λ _i . формы ясно, что минимальный полином A имеет степень Σs i _. Из жордановой нормальной

Хотя нормальная форма Жордана определяет минимальный полином, обратное неверно. Это приводит к понятию элементарных делителей . Элементарные делители квадратной матрицы A — это характеристические многочлены ее жордановых блоков. Факторами минимального многочлена m являются элементарные делители наибольшей степени, соответствующие различным собственным значениям.

Степень элементарного дивизора — это размер соответствующего жорданова блока, а следовательно, и размерность соответствующего инвариантного подпространства. Если все элементарные делители линейны, A диагонализуемо.

Инвариантные разложения подпространств

Йордановая форма размера n × n матрицы A является блочно-диагональной и, следовательно, дает разложение n- мерного евклидова пространства на инвариантные подпространства A . Каждый жордановый блок J _i соответствует инвариантному подпространству X _i . Символически мы положили

$\mathbb {C} ^{n}=\bigoplus _{i=1}^{k}X_{i}$

где каждый X _i — это длина соответствующей жордановой цепи, а k — количество жордановых цепей.

Можно также получить несколько иное разложение через жорданову форму. собственное значение λ _i , размер его наибольшего соответствующего жорданового блока s _i называется индексом λ i _{Учитывая} и обозначается v ( λ _i ). (Следовательно, степень минимального многочлена равна сумме всех индексов.) Определим подпространство Y _i формулой

$Y_{i}=\ker(\lambda _{i}I-A)^{v(\lambda _{i})}.$

Это дает разложение

$\mathbb {C} ^{n}=\bigoplus _{i=1}^{l}Y_{i}$

где l — количество различных собственных значений A . Интуитивно мы объединяем вместе жордановы блочно-инвариантные подпространства, соответствующие одному и тому же собственному значению. В крайнем случае, когда A кратно единичной матрице, мы имеем k = n и l = 1.

Проекция на Y _i и вдоль всех остальных Y _j ( j ≠ i ) называется спектральной проекцией A в точке _vi и обычно обозначается P ( λ _i ; A ) . Спектральные проекции взаимно ортогональны в том смысле, что P ( λ _i ; A ) P (v _j ; A ) = 0, если я ≠ j . Кроме того, они коммутируют с A , и их сумма является единичной матрицей. Замена каждого v _i в йордановой матрице J на ​​единицу и обнуление всех остальных элементов дает P (vi _; J ) , причем, если UJU ⁻¹ — преобразование подобия такое, что A = UJU ⁻¹ тогда п ( λ _я ; А ) знак равно UP ( λ _я ; J ) U ⁻¹ . Они не ограничиваются конечными размерами. См. ниже их применение к компактным операторам и в голоморфном функциональном исчислении более общее обсуждение .

Сравнивая два разложения, обратите внимание, что, вообще говоря, l ≤ k . Когда A нормально, подпространства X _i в первом разложении одномерны и взаимно ортогональны. Это спектральная теорема для нормальных операторов. Второе разложение легче обобщается на общие компактные операторы в банаховых пространствах.

Здесь может быть интересно отметить некоторые свойства индекса ν ( λ ). В более общем смысле, для комплексного числа λ его индекс можно определить как наименьшее неотрицательное целое число ν ( λ ), такое что

$\ker(A-\lambda I)^{\nu (\lambda )}=\ker(A-\lambda I)^{m},\;\forall m\geq \nu (\lambda ).$

Итак, ν (v) > 0 тогда и только тогда, когда λ является собственным значением A . В конечномерном случае ν (v) ⩽ алгебраическая кратность v.

Плоская (плоская) нормальная форма

Форма Жордана используется для нахождения нормальной формы матриц с точностью до сопряженности, такой, что нормальные матрицы составляют алгебраическое многообразие низкой фиксированной степени в объемлющем матричном пространстве.

Множества представителей классов матричной сопряженности для жордановой нормальной формы или рациональных канонических форм вообще не составляют линейных или аффинные подпространства в объемлющих матричных пространствах.

Владимир Арнольд позировал проблема:Найти канонический вид матриц над полем, для которого множество представителей классов сопряженности матриц представляет собой объединение аффинных линейных подпространств (квартир). Другими словами, инъективно отобразите набор классов сопряженности матриц обратно в исходный набор матриц так, чтобы образ этого вложения - набор всех нормальных матриц имел наименьшую возможную степень - это было объединение сдвинутых линейных подпространств.

Для алгебраически замкнутых полей ее решил Петерис Даугулис. Построение однозначно определенной плоской нормальной формы матрицы начинается с рассмотрения ее жордановой нормальной формы.

Матричные функции

Основная статья: Матричная функция

Итерация цепочки Джордана мотивирует различные расширения до более абстрактных настроек. Для конечных матриц получаются матричные функции; это можно распространить на компактные операторы и голоморфное функциональное исчисление, как описано ниже.

Жорданова нормальная форма является наиболее удобной для вычисления матричных функций (хотя она может быть не лучшим выбором для компьютерных вычислений). Пусть f ( z ) — аналитическая функция комплексного аргумента. Применение функции к размера n × n жордановому блоку J с собственным значением λ приводит к получению верхней треугольной матрицы:

$f(J)={\begin{bmatrix}f(\lambda )&f'(\lambda )&{\tfrac {f''(\lambda )}{2}}&\cdots &{\tfrac {f^{(n-1)}(\lambda )}{(n-1)!}}\\0&f(\lambda )&f'(\lambda )&\cdots &{\tfrac {f^{(n-2)}(\lambda )}{(n-2)!}}\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&0&f(\lambda )&f'(\lambda )\\0&0&0&0&f(\lambda )\end{bmatrix}},$

так что элементы k -й супердиагонали полученной матрицы равны ${\tfrac {f^{(k)}(\lambda )}{k!}}$. Для матрицы общей жордановой нормальной формы приведенное выше выражение должно применяться к каждому жорданову блоку.

В следующем примере показано применение степенной функции f ( z ) = z ^н :

${\begin{bmatrix}\lambda _{1}&1&0&0&0\\0&\lambda _{1}&1&0&0\\0&0&\lambda _{1}&0&0\\0&0&0&\lambda _{2}&1\\0&0&0&0&\lambda _{2}\end{bmatrix}}^{n}={\begin{bmatrix}\lambda _{1}^{n}&{\tbinom {n}{1}}\lambda _{1}^{n-1}&{\tbinom {n}{2}}\lambda _{1}^{n-2}&0&0\\0&\lambda _{1}^{n}&{\tbinom {n}{1}}\lambda _{1}^{n-1}&0&0\\0&0&\lambda _{1}^{n}&0&0\\0&0&0&\lambda _{2}^{n}&{\tbinom {n}{1}}\lambda _{2}^{n-1}\\0&0&0&0&\lambda _{2}^{n}\end{bmatrix}},$

где биномиальные коэффициенты определяются как ${\textstyle {\binom {n}{k}}=\prod _{i=1}^{k}{\frac {n+1-i}{i}}}$. Для целого положительного n оно сводится к стандартному определению.коэффициентов. Для отрицательного n тождество ${\textstyle {\binom {-n}{k}}=(-1)^{k}{\binom {n+k-1}{k}}}$ может оказаться полезным.

Компактные операторы

Результат, аналогичный жордановой нормальной форме, справедлив для компактных операторов в банаховом пространстве . Ограничиваются компактными операторами, поскольку каждая точка x в спектре компактного оператора T является собственным значением; Единственным исключением является случай, когда x является предельной точкой спектра. Это неверно для ограниченных операторов в целом. Чтобы дать некоторое представление об этом обобщении, сначала переформулируем разложение Жордана на языке функционального анализа.

Голоморфное функциональное исчисление

Дополнительная информация: голоморфное функциональное исчисление.

Пусть X — банахово пространство, L ( X ) — ограниченные операторы на X , а ( T ) обозначает спектр T L ∈ σ ( X ). Голоморфное функциональное исчисление определяется следующим образом:

Зафиксируйте ограниченный оператор T . Рассмотрим семейство Hol( T комплексных функций ), голоморфное на некотором открытом множестве G, содержащем σ ( T ). Пусть Γ = { γi _} — конечный набор жордановых кривых такой, что σ ( T ) лежит внутри Γ , мы определяем f ( T ) как

$f(T)={\frac {1}{2\pi i}}\int _{\Gamma }f(z)(z-T)^{-1}\,dz.$

Открытое множество G может меняться в зависимости от f и не обязательно должно быть связным. Интеграл определяется как предел сумм Римана, как и в скалярном случае. Хотя интеграл имеет смысл для непрерывного f , мы ограничиваемся голоморфными функциями, чтобы применить аппарат классической теории функций (например, интегральную формулу Коши). Предположение о том, что σ ( T ) лежит внутри Γ, обеспечивает f ( T корректность определения ); оно не зависит от выбора Γ. Функциональное исчисление — это отображение Φ из Hol( T ) в L ( X ), заданное формулой

$\;\Phi (f)=f(T).$

Нам потребуются следующие свойства этого функционального исчисления:

1. Φ расширяет полиномиальное функциональное исчисление.
2. Верна теорема о спектральном отображении : σ ( f ( T )) = f ( σ ( T )).
3. Φ — гомоморфизм алгебры.

Конечномерный случай

В конечномерном случае σ ( T ) = { λ _i } — конечное дискретное множество в комплексной плоскости. Пусть e _i будет функцией, которая равна 1 в некоторой открытой окрестности λ _i и равна 0 в другом месте. По свойству 3 функционального исчисления оператор

$e_{i}(T)$

является проекцией. , пусть νi _— индекс λi _и Более того

$f(z)=(z-\lambda _{i})^{\nu _{i}}.$

Теорема о спектральном отображении говорит нам

$f(T)e_{i}(T)=(T-\lambda _{i})^{\nu _{i}}e_{i}(T)$

имеет спектр {0}. По свойству 1 f ( T ) может быть непосредственно вычислено в жордановой форме, и при проверке мы видим, что оператор f ( T ) e _i ( T ) является нулевой матрицей.

По свойству 3 f ( T ) e _i ( T ) = e _i ( T ) f ( T ). Итак, e _i ( T ) — это в точности проекция на подпространство

$\operatorname {Ran} e_{i}(T)=\ker(T-\lambda _{i})^{\nu _{i}}.$

Отношение

$\sum _{i}e_{i}=1$

подразумевает

$\mathbb {C} ^{n}=\bigoplus _{i}\;\operatorname {Ran} e_{i}(T)=\bigoplus _{i}\ker(T-\lambda _{i})^{\nu _{i}}$

где индекс i проходит через различные собственные значения T . Это инвариантное разложение подпространства

$\mathbb {C} ^{n}=\bigoplus _{i}Y_{i}$

приведено в предыдущем разделе. Каждое e _i ( T ) является проекцией на подпространство, натянутое жордановыми цепями, соответствующими λ _i, и вдоль подпространств, натянутых жордановыми цепями, соответствующими v _j, для j ≠ i . Другими словами, е _я ( Т ) = P ( λ _я ; Т ). Эта явная идентификация операторов ei _{) ,} ( T в свою очередь, дает явную форму голоморфного функционального исчисления для матриц:

Для всех f ∈ Hol( T ),

$f(T)=\sum _{\lambda _{i}\in \sigma (T)}\sum _{k=0}^{\nu _{i}-1}{\frac {f^{(k)}}{k!}}(T-\lambda _{i})^{k}e_{i}(T).$

Обратите внимание, что выражение f ( T ) является конечной суммой, потому что в каждой окрестности vi _мы выбрали в ряд Тейлора разложение f с центром в _vi .

Полюса оператора

Пусть T — ограниченный оператор λ — изолированная точка σ ( T ). (Как указано выше, когда T компактно, каждая точка его спектра является изолированной точкой, за исключением, возможно, предельной точки 0.)

Точка λ называется полюсом оператора T порядка ν , если резольвентная функция R _T определяется формулой

$R_{T}(\lambda )=(\lambda -T)^{-1}$

имеет полюс порядка ν в точке λ .

Покажем, что в конечномерном случае порядок собственного значения совпадает с его индексом. Этот результат справедлив и для компактных операторов.

Рассмотрим кольцевую область A с центром в собственном значении λ и достаточно малым радиусом ε такую, что пересечение открытого диска B _ε ( λ ) и σ ( T ) равно { λ }. Резольвентная функция R _T голоморфна на A .Расширяя результат классической теории функций, имеет _RT представление ряда Лорана на A :

$R_{T}(z)=\sum _{-\infty }^{\infty }a_{m}(\lambda -z)^{m}$

где

$a_{-m}=-{\frac {1}{2\pi i}}\int _{C}(\lambda -z)^{m-1}(z-T)^{-1}dz$ и C — небольшой круг с центром в точке λ .

Судя по предыдущему обсуждению функционального исчисления,

$a_{-m}=-(\lambda -T)^{m-1}e_{\lambda }(T)$ где $e_{\lambda }$ 1 на $B_{\varepsilon }(\lambda )$ и 0 в другом месте.

Но мы показали, что наименьшее натуральное число m такое, что

$a_{-m}\neq 0$ и $a_{-l}=0\;\;\forall \;l\geq m$

это в точности индекс λ , ν ( λ ). Другими словами, функция R _T имеет полюс порядка ν ( λ ) в точке λ .

Численный анализ

Если матрица A имеет несколько собственных значений или близка к матрице с несколькими собственными значениями, то ее жордановая нормальная форма очень чувствительна к возмущениям. Рассмотрим, например, матрицу

$A={\begin{bmatrix}1&1\\\varepsilon &1\end{bmatrix}}.$

Если ε = 0, то жордановая нормальная форма — это просто

${\begin{bmatrix}1&1\\0&1\end{bmatrix}}.$

Однако при ε ≠ 0 жордановая нормальная форма имеет вид

${\begin{bmatrix}1+{\sqrt {\varepsilon }}&0\\0&1-{\sqrt {\varepsilon }}\end{bmatrix}}.$

Из-за этой плохой обусловленности очень сложно разработать надежный численный алгоритм для нормальной формы Жордана, поскольку результат критически зависит от того, считаются ли два собственных значения равными. обычно избегают жордановой нормальной формы По этой причине при численном анализе ; стабильное разложение Шура или псевдоспектры являются лучшими альтернативами.

См. также

• Каноническая основа
• Каноническая форма
• Нормальная форма Фробениуса
• Жорданова матрица
• Разложение Жордана – Шевалле
• Разложение матрицы
• Модальная матрица
• Каноническая форма Вейра

Ссылки

• Борегар, Раймонд А.; Фрэли, Джон Б. (1973), Первый курс линейной алгебры: с дополнительным введением в группы, кольца и поля , Бостон: Houghton Mifflin Co. , ISBN  0-395-14017-Х
• Бронсон, Ричард (1970), Матричные методы: введение , Нью-Йорк: Academic Press , LCCN   70097490
• Каллен, Чарльз Г. (1966), Матрицы и линейные преобразования , Чтение: Аддисон-Уэсли , LCCN   66021267
• Данфорд, Н.; Шварц, Дж.Т. (1958), Линейные операторы, Часть I: Общая теория , Interscience
• Финкбайнер II, Дэниел Т. (1978), Введение в матрицы и линейные преобразования (3-е изд.), WH Freeman and Company
• Франклин, Джоэл Н. (1968), Теория матриц , Энглвуд Клиффс: Прентис-Холл , LCCN   68016345
• Голуб, Джин Х.; Ван Лоан, Чарльз Ф. (1996), Матричные вычисления (3-е изд.), Балтимор: Издательство Университета Джона Хопкинса , ISBN  0-8018-5414-8
• Голуб, Джин Х.; Уилкинсон, Дж. Х. (1976). «Плохо обусловленные собственные системы и вычисление жордановой нормальной формы». Обзор СИАМ . 18 (4): 578–619. дои : 10.1137/1018113 .
• Holt, Derek; Rumynin, Dmitriy (2009), Algebra I – Advanced Linear Algebra (MA251) Lecture Notes (PDF)
• Хорн, Роджер А.; Джонсон, Чарльз Р. (1985), Матричный анализ , Издательство Кембриджского университета , ISBN  978-0-521-38632-6
• Джеймс, Гленн; Джеймс, Роберт К. (1976), Математический словарь (2-е изд.), Ван Ностранд Рейнхольд
• Маклейн, Сондерс; Биркгоф, Гаррет (1967), Алгебра , Macmillan Publishers
• Мишель, Энтони Н.; Херже, Чарльз Дж. (1993), Прикладная алгебра и функциональный анализ , Dover Publications
• Неринг, Эвар Д. (1970), Линейная алгебра и теория матриц (2-е изд.), Нью-Йорк: Wiley , LCCN   76091646
• Шафаревич, ИР; Ремизов, АО (2012), Линейная алгебра и геометрия , Springer , ISBN  978-3-642-30993-9
• Shilov, Georgi E. (1977), Linear Algebra , Dover Publications
• Статья о канонической форме Джордана на mathworld.wolfram.com