представление обратимых матриц в виде унитарных умножений на эрмитов оператор
В математике символ полярное разложение квадрата вещественной или сложной матрицы - это факторизация формы , где - это унитарная матрица и представляет собой позитивно-полуопределённую эрмитову матрицу, квадратную и одинакового размера.
Интуитивно, если действительная матрица интерпретируется как линейное преобразование из -размерное пространство , полярное разложение разделяет его на вращение или отражение из , и масштабирование пространства вдоль набор ортогональных осей .
Полярное разложение квадратной матрицы всегда существует. Если является обратимым, разложение будет уникальным, и коэффициент будет положительно-определенный. В этом случае можно записать однозначно в форме , где - унитарный, а - уникальный самосопряженный логарифм матрица . Это разложение полезно при вычислении фундаментальной группы из (матрицы) групп Ли.
Полярное разложение также можно определить как где - симметричная положительно-определенная, но в целом другая матрица, а - та же матрица, что и выше.
Полярное разложение матрицы можно рассматривать как матричный аналог полярной формы комплексного числа as , где - его абсолютное значение (неотрицательное действительное число ), а - комплексное число с единичной нормой (элемент группы кругов ).
Содержание
- 1 Свойства
- 1.1 Интуитивная интерпретация
- 1.2 Связь с SVD
- 1.3 Конструкция и доказательства существования
- 1.3.1 Случай нормальный
- 1.3.2 Случай обратимый
- 1.3.3 Общий случай
- 2 Ограниченные операторы в гильбертовом пространстве
- 3 Неограниченные операторы
- 4 Кватернионное полярное разложение
- 5 Альтернативные планарные разложения
- 6 Численное определение матричного полярного разложения
- 7 См. Также
- 8 Ссылки
Свойства
полярное разложение комплексно-сопряженного элемента дается выражением Обратите внимание, что
дает соответствующее полярное разложение детерминанта числа A, поскольку и . В частности, если имеет определитель 1, то оба и имеют определитель 1.
Положительно-полуопределенная матрица P всегда уникальна, даже если A сингулярна, и обозначается как
где A обозначает сопряженное транспонирование из A. Уникальность P гарантирует, что это выражение хорошо определено. Уникальность гарантируется тем, что является положительно-полуопределенной эрмитовой матрицей и, следовательно, имеет единственную положительно-полуопределенную эрмитову матрицу квадратный корень. Если A обратимо, то P положительно определено, следовательно, также обратимо, и матрица U однозначно определяется с помощью
Интуитивная интерпретация
Действительный квадрат матрица можно представить как линейное преобразование для , который принимает вектор-столбец от до . Тогда в полярном разложении коэффициент равен вещественная ортонормированная матрица. Полярное разложение в таком случае можно рассматривать как выражение линейного преобразования, определенного , в масштабирование пространства вдоль каждого собственного вектора of на коэффициент масштабирования (действие ), за которым следует однократное вращение или отражение (действие ).
В качестве альтернативы, разложение выражает преобразование, определенное как поворот () с последующим масштабированием () в определенных ортогональных направлениях. Коэффициенты масштабирования такие же, но направления разные.
Связь с SVD
В терминах разложения по сингулярным значениям (SVD) , , один имеет
где , , и - унитарные матрицы (называемые ортогональными матрицами, если поле является вещественным ). Это подтверждает, что является положительно определенным, а унитарным. Таким образом, существование SVD эквивалентно существованию полярного разложения.
Можно также разложить в форме
Здесь то же самое, что и раньше, а задается как
Это называется левым полярным разложением, тогда как предыдущее разложение известно как правое полярное разложение. Левополярное разложение также известно как обратное полярное разложение.
Матрица является нормальной тогда и только тогда, когда . Тогда , и можно диагонализовать с помощью унитарная матрица сходства , которая коммутирует с , давая , где - диагональная унитарная матрица фаз . Положив , можно переписать полярное разложение как
так что , то, таким образом, также имеет спектральное разложение
с комплексными собственными значениями, такими как что и унитарная матрица комплексных собственных векторов .
Полярное разложение квадратной обратимой вещественной матрицы имеет форму
, где - это положительно определенная матрица и - ортогональная матрица.
Построение и доказательства существования
Основная идея построения полярного разложения аналогична той, которая используется для вычисления разложения по сингулярным числам.
Для любого , матрица эрмитова и положительно полуопределенная, и поэтому унитарно эквивалентна положительная полуопределенная диагональная матрица. Пусть тогда будет унитарным таким, что , с диагональю и положительным полуопределенным.
Случай нормальный
Если нормальный, то он унитарно эквивалентен диагональной матрице: для некоторого унитарного и некоторая диагональная матрица . Тогда мы можем написать
где
- диагональная матрица, содержащая фазы элементов
, то есть
или
произвольное комплексное число с единичной величиной, когда
.
Таким образом, полярное разложение имеет вид , где и диагональ на основе собственных значений и с собственными значениями, равными фазам и абсолютным значениям собственных значений соответственно.
Случай обратимый
Из разложения по единственному числу можно показать, что обратимо тогда и только тогда, когда (эквивалентно, ) есть. Более того, это верно тогда и только тогда, когда все собственные значения не равны нулю.
В этом случае полярное разложение непосредственно получается записью
и наблюдая, что унитарен. Чтобы убедиться в этом, мы можем использовать спектральное разложение для записи .
В этом выражении унитарен, потому что есть. Чтобы показать, что является унитарным, мы можем использовать SVD написать , так что
где снова унитарен по построению.
Еще один способ напрямую показать унитарность означает, что запись SVD из с точки зрения ранга-1 матрицы как , где - сингулярные значения , мы имеем
что прямо подразумевает унитарность
, потому что матрица унитарна тогда и только тогда, когда ее сингулярные значения имеют унитарное абсолютное значение.
Обратите внимание, как из приведенной выше конструкции следует, что унитарная матрица в полярном разложении обратимой матрицы определена однозначно.
Общий случай
SVD читает , с унитарными матрицами и диагональная положительно полуопределенная матрица. Просто вставив дополнительную пару s или s, мы получим две формы полярного разложения из :
Ограниченные операторы в гильбертовом пространстве
полярное разложение любого линейного ограниченного оператора A между комплексными гильбертовыми пространствами - это каноническая факторизация как произведение частичной изометрии и неотрицательного оператора.
Полярное разложение для матриц обобщается следующим образом: если A - ограниченный линейный оператор, то существует уникальная факторизация A как произведения A = UP, где U - частичная изометрия, P - неотрицательная самость. -сопряженный оператор и начальное пространство U является закрытием диапазона P.
Оператор U должен быть ослаблен до частичной изометрии, а не унитарной, из-за следующих проблем. Если A - односторонний сдвиг на l (N ), то | A | = {AA} = I. Итак, если A = U | A |, U должно быть A, которое не является унитарным.
Существование полярного разложения является следствием леммы Дугласа :
- Лемма Если A, B - ограниченные операторы в гильбертовом пространстве H и AA ≤ BB, то существует сжатие C такое, что A = CB. Кроме того, C уникален, если Ker (B) ⊂ Ker (C).
Оператор C может быть определен как C (Bh): = Ah для всех h в H, продолженный по непрерывности до замыкания Ran (B), и нулем на ортогональном дополнении ко всему H. Из этого следует, что из AA ≤ BB следует Ker (B) ⊂ Ker (A).
В частности. Если AA = BB, то C - частичная изометрия, которая единственна, если Ker (B) ⊂ Ker (C). В общем, для любого ограниченного оператора A
где (AA) - единственный положительный квадратный корень из AA, полученный с помощью обычного функционального исчисления. Итак, по лемме
для некоторой частичной изометрии U, которая единственна, если Ker (A) ⊂ Ker (U). Возьмем P равным (AA), и получим полярное разложение A = UP. Обратите внимание, что аналогичный аргумент можно использовать, чтобы показать A = P'U ', где P' положительно, а U '- частичная изометрия.
Когда H конечномерно, U может быть расширен до унитарного оператора; в целом это неверно (см. пример выше). В качестве альтернативы полярное разложение может быть показано с использованием операторной версии разложения по сингулярным значениям.
По свойству непрерывного функционального исчисления, | A | находится в C * -алгебре, порожденной A. Аналогичное, но более слабое утверждение справедливо для частичной изометрии: U находится в алгебре фон Неймана, порожденной A. Если A обратима, полярная часть U также будет в C * -алгебре.
Неограниченные операторы
Если A - замкнутый, плотно определенный неограниченный оператор между комплексными гильбертовыми пространствами, то он все еще имеет (уникальное) полярное разложение
где | A | является (возможно, неограниченным) неотрицательным самосопряженным оператором с той же областью определения, что и A, а U - частичная изометрия, исчезающая на ортогональном дополнении диапазона Ran (| A |).
Доказательство использует ту же лемму, что и выше, которая проходит для неограниченных операторов в целом. Если Dom (AA) = Dom (BB) и AAh = BBh для всех h ∈ Dom (AA), то существует частичная изометрия U такая, что A = UB. U единственно, если Ran (B) ⊂ Ker (U). Замкнутый и плотно определенный оператор A гарантирует, что оператор AA является самосопряженным (с плотной областью определения) и, следовательно, позволяет определить (AA). Применение леммы дает полярное разложение.
Если неограниченный оператор A связан с алгеброй фон Неймана M, и A = UP является его полярным разложением, то U находится в M, а также спектральная проекция P, 1 B (P), для любого борелевского множества B в [0, ∞).
Полярное разложение кватернионов
Полярное разложение кватернионов Hзависит от единичной 2-мерной сферы из квадратный корень из минус единицы. Для любого r на этой сфере и угла −π < a ≤ π, the versor находится на блоке 3-сфера из H . Для a = 0 и a = π версор равен 1 или -1 независимо от того, какой r выбран. norm t кватерниона q - это евклидово расстояние от начала координат до q. Когда кватернион - это не просто действительное число, тогда существует уникальное полярное разложение
Альтернативные плоские разложения
В декартовой плоскости альтернативные плоские кольцевые разложения возникают следующим образом:
- Если x ≠ 0, z = x (1 + ε (y / x)) - полярное разложение двойственного числа z = x + yε, где ε = 0; то есть ε нильпотентен. В этом полярном разложении единичный круг заменен линией x = 1, полярный угол - наклоном y / x, а радиус x отрицателен в левой полуплоскости.
- Если x ≠ y, то единичная гипербола x - y = 1 и ее сопряженная x - y = −1 могут быть использованы для формирования полярного разложения на основе ветви единичной гиперболы через ( 1, 0). Эта ветвь параметризуется гиперболическим углом a и записывается как
где j = +1 и используется арифметика комплексных чисел с разбиением. Ветвь через (−1, 0) отслеживается −e. Поскольку операция умножения на j отражает точку на прямой y = x, вторая гипербола имеет ветви, обозначенные je или −je. Следовательно, точка в одном из квадрантов имеет полярное разложение в одной из форм:
В наборе {1, -1, j, -j} есть продукты, которые делают его изоморфным четверке Клейна -group. Очевидно, полярное разложение в этом случае включает элемент из этой группы.
Численное определение полярного разложения матрицы
Для вычисления приближения полярного разложения A = UP, обычно унитарного коэффициент U приближается. Итерация основана на методе Герона для квадратного корня из 1 и вычисляет, начиная с , последовательность
Комбинация инверсии а сопряжение Эрмита выбрано таким образом, чтобы при разложении по сингулярным значениям унитарные множители оставались неизменными, а итерация сводилась к методу Герона по сингулярным значениям.
Эта базовая итерация может быть уточнена для ускорения процесса:
- Каждый шаг или через равные промежутки времени, диапазон сингулярных значений оценивается, а затем матрица масштабируется до , чтобы центрировать сингулярные значения вокруг 1. Масштабирование Фактор вычисляется с использованием норм матрицы и обратной матрицы. Примеры таких масштабных оценок:
с использованием суммы строки и сумма столбцов матричные нормы или
с использованием Фробениуса норма. Включая масштабный коэффициент, итерация теперь имеет вид
- QR-разложение можно использовать на этапе подготовки уменьшить сингулярную матрицу A до регулярной матрицы меньшего размера и внутри каждого шага ускорить вычисление обратной.
- Метод Герона для вычисления корней можно заменить методами более высокого порядка, например, на основе метода Галлея третьего порядка, в результате чего получится Эту итерацию снова можно комбинировать с изменением масштаба. Эта конкретная формула имеет то преимущество, что она также применима к сингулярным или прямоугольным матрицам A.
См. Также
Ссылки
- Конвей, Дж. Б. : курс функционального анализа. Тексты для выпускников по математике. Нью-Йорк: Springer 1990
- Дуглас, Р.Г. : О мажоризации, факторизации и включении диапазонов операторов в гильбертовом пространстве. Proc. Амер. Математика. Soc. 17, 413-415 (1966)
- Холл, Брайан К. (2015), Группы Ли, алгебры Ли и представления: элементарное введение, Тексты для выпускников по математике, 222 ( 2-е изд.), Springer, ISBN 978-3319134666.
- Helgason, Sigurdur (1978), Дифференциальная геометрия, группы Ли и симметричные пространства, Academic Press, ISBN 0-8218-2848-7