Полярное разложение

редактировать
представление обратимых матриц в виде унитарных умножений на эрмитов оператор

В математике символ полярное разложение квадрата вещественной или сложной матрицы A {\ displaystyle A}A - это факторизация формы A = UP {\ displaystyle A = UP}{\ displaystyle A = UP} , где U {\ displaystyle U}U - это унитарная матрица и P {\ displaystyle P}P представляет собой позитивно-полуопределённую эрмитову матрицу, квадратную и одинакового размера.

Интуитивно, если действительная n × n {\ displaystyle n \ times n}n \ times n матрица A {\ displaystyle A}A интерпретируется как линейное преобразование из n {\ displaystyle n}n -размерное пространство R n {\ displaystyle \ mathbb {R} ^ {n} }\ mathbb {R} ^ {n} , полярное разложение разделяет его на вращение или отражение U {\ displaystyle U}U из R n {\ displaystyle \ mathbb {R} ^ {n}}\ mathbb {R} ^ {n} , и масштабирование пространства вдоль набор ортогональных осей n {\ displaystyle n}n .

Полярное разложение квадратной матрицы A {\ displaystyle A}A всегда существует. Если A {\ displaystyle A}A является обратимым, разложение будет уникальным, и коэффициент P {\ displaystyle P}P будет положительно-определенный. В этом случае A {\ displaystyle A}A можно записать однозначно в форме A = U e X {\ displaystyle A = Ue ^ {X}}{ \ Displaystyle A = Ue ^ {X}} , где U {\ displaystyle U}U - унитарный, а X {\ displaystyle X}X - уникальный самосопряженный логарифм матрица P {\ displaystyle P}P . Это разложение полезно при вычислении фундаментальной группы из (матрицы) групп Ли.

Полярное разложение также можно определить как A = PU {\ displaystyle A = PU}{\ displaystyle A = PU } где P {\ displaystyle P}P - симметричная положительно-определенная, но в целом другая матрица, а U {\ displaystyle U}U - та же матрица, что и выше.

Полярное разложение матрицы можно рассматривать как матричный аналог полярной формы комплексного числа z {\ displaystyle z}z as z = ur {\ displaystyle z = ur}{\ displaystyle z = ur} , где r {\ displaystyle r}r - его абсолютное значение (неотрицательное действительное число ), а u {\ displaystyle u}u - комплексное число с единичной нормой (элемент группы кругов ).

Содержание
  • 1 Свойства
    • 1.1 Интуитивная интерпретация
    • 1.2 Связь с SVD
    • 1.3 Конструкция и доказательства существования
      • 1.3.1 Случай A {\ displaystyle A}A нормальный
      • 1.3.2 Случай A {\ displaystyle A}A обратимый
      • 1.3.3 Общий случай
  • 2 Ограниченные операторы в гильбертовом пространстве
  • 3 Неограниченные операторы
  • 4 Кватернионное полярное разложение
  • 5 Альтернативные планарные разложения
  • 6 Численное определение матричного полярного разложения
  • 7 См. Также
  • 8 Ссылки
Свойства

полярное разложение комплексно-сопряженного элемента A {\ displaystyle A}A дается выражением A ¯ = U ¯ P ¯. {\ displaystyle {\ overline {A}} = {\ overline {U}} {\ overline {P}}.}{\ displaystyle {\ overline {A}} = {\ overline {U}} {\ overline {P }}.} Обратите внимание, что

det A = det U det P = ei θ ⋅ r {\ displaystyle \ det A = \ det U \ det P = e ^ {i \ theta} \ cdot r}{\ displaystyle \ det A = \ det U \ det P = e ^ {i \ theta} \ cdot r}

дает соответствующее полярное разложение детерминанта числа A, поскольку det U = ei θ {\ displaystyle \ det U = e ^ {i \ theta}}\ Det U = е ^ {{я \ theta}} и det P = r = | det A | {\ Displaystyle \ Det P = R = | \ Det A |}\ det P = r = | \ det A | . В частности, если A {\ displaystyle A}A имеет определитель 1, то оба U {\ displaystyle U}U и P {\ displaystyle P}P имеют определитель 1.

Положительно-полуопределенная матрица P всегда уникальна, даже если A сингулярна, и обозначается как

P = (A ∗ A) 1 2, {\ displaystyle P = \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}{\ displaystyle P = \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}

где A обозначает сопряженное транспонирование из A. Уникальность P гарантирует, что это выражение хорошо определено. Уникальность гарантируется тем, что A ∗ A {\ displaystyle A ^ {*} A}A ^ {*} A является положительно-полуопределенной эрмитовой матрицей и, следовательно, имеет единственную положительно-полуопределенную эрмитову матрицу квадратный корень. Если A обратимо, то P положительно определено, следовательно, также обратимо, и матрица U однозначно определяется с помощью

U = A P - 1. {\ displaystyle U = AP ^ {- 1}.}{\ displaystyle U = AP ^ {- 1}.}

Интуитивная интерпретация

Действительный квадрат m × m {\ displaystyle m \ times m}м \ раз м матрица A {\ displaystyle A}A можно представить как линейное преобразование для R m {\ displaystyle \ mathbb {R} ^ {m}}\ mathbb {R} ^ {m} , который принимает вектор-столбец от x {\ displaystyle x}x до A x {\ displaystyle Ax}Ax . Тогда в полярном разложении A = RP {\ displaystyle A = RP}{\ displaystyle A = RP } коэффициент R {\ displaystyle R}Р равен m × m {\ displaystyle m \ times m}м \ раз м вещественная ортонормированная матрица. Полярное разложение в таком случае можно рассматривать как выражение линейного преобразования, определенного A {\ displaystyle A}A , в масштабирование пространства R m {\ displaystyle \ mathbb {R} ^ {m}}\ mathbb {R} ^ {m} вдоль каждого собственного вектора ei {\ displaystyle e_ {i}}e_ {i} of A {\ displaystyle A}A на коэффициент масштабирования σ я {\ displaystyle \ sigma _ {i}}\ sigma _ {i} (действие P {\ displaystyle P}P ), за которым следует однократное вращение или отражение R m {\ displaystyle \ mathbb {R} ^ {m}}\ mathbb {R} ^ {m} (действие R {\ displaystyle R}Р ).

В качестве альтернативы, разложение A = PR {\ displaystyle A = PR}{\ displaystyle A = PR} выражает преобразование, определенное A {\ displaystyle A}A как поворот (R {\ displaystyle R}Р ) с последующим масштабированием (P {\ displaystyle P}P ) в определенных ортогональных направлениях. Коэффициенты масштабирования такие же, но направления разные.

Связь с SVD

В терминах разложения по сингулярным значениям (SVD) A {\ displaystyle A}A , A = W Σ V ∗ {\ displaystyle A = W \ Sigma V ^ {*}}{\ displaystyle A = W \ Sigma V ^ {*}} , один имеет

P = V Σ V ∗ U = WV ∗ {\ displaystyle {\ begin {выровнено} P = V \ Сигма V ^ {*} \\ U = WV ^ {*} \ end {align}}}{\ displaystyle {\ begin {align} P = V \ Sigma V ^ {*} \\ U = WV ^ {*} \ end {align}}}

где U {\ displaystyle U}U , V {\ displaystyle V}V , и W {\ displaystyle W}W - унитарные матрицы (называемые ортогональными матрицами, если поле является вещественным R {\ displaystyle \ mathbb {R}}\ mathbb {R} ). Это подтверждает, что P {\ displaystyle P}P является положительно определенным, а U {\ displaystyle U}U унитарным. Таким образом, существование SVD эквивалентно существованию полярного разложения.

Можно также разложить A {\ displaystyle A}A в форме

A = P ′ U {\ displaystyle A = P'U}{\displaystyle A=P'U}

Здесь U {\ displaystyle U}U то же самое, что и раньше, а P ′ {\ displaystyle P '}P'задается как

P ′ = UPU - 1 = (AA ∗) 1 2 = W Σ W ∗. {\ displaystyle P '= UPU ^ {- 1} = \ left (AA ^ {*} \ right) ^ {\ frac {1} {2}} = W \ Sigma W ^ {*}.}{\displaystyle P'=UPU^{-1}=\left(AA^{*}\right)^{\frac {1}{2}}=W\Sigma W^{*}.}

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

Матрица A {\ displaystyle A}A является нормальной тогда и только тогда, когда P '= P {\ displaystyle P' = P}{\displaystyle P'=P}. Тогда U Σ = Σ U {\ displaystyle U \ Sigma = \ Sigma U}{\ displaystyle U \ Sigma = \ Sigma U} , и можно диагонализовать U {\ displaystyle U}U с помощью унитарная матрица сходства S {\ displaystyle S}S , которая коммутирует с Σ {\ displaystyle \ Sigma}\ Sigma , давая SUS ∗ = Φ - 1 {\ displaystyle SUS ^ {*} = \ Phi ^ {- 1}}{\ displaystyle SUS ^ {*} = \ Phi ^ {- 1}} , где Φ {\ displaystyle \ Phi}\ Phi - диагональная унитарная матрица фаз ei φ {\ Displaystyle е ^ {я \ varphi}}{\ displaystyle e ^ {i \ varphi}} . Положив Q = VS ∗ {\ displaystyle Q = VS ^ {*}}{\ displaystyle Q = VS ^ {*}} , можно переписать полярное разложение как

A = (Q Φ Q ∗) (Q Σ Q *), {\ displaystyle A = \ left (Q \ Phi Q ^ {*} \ right) \ left (Q \ Sigma Q ^ {*} \ right), \,}{\ displaystyle A = \ left (Q \ Phi Q ^ {*} \ right) \ left (Q \ Sigma Q ^ {*} \ right), \,}

так что A { \ displaystyle A}A , то, таким образом, также имеет спектральное разложение

A = Q Λ Q ∗ {\ displaystyle A = Q \ Lambda Q ^ {*}}{\ displaystyle A = Q \ Lambda Q ^ {*}}

с комплексными собственными значениями, такими как что Λ Λ ∗ = Σ 2 {\ displaystyle \ Lambda \ Lambda ^ {*} = \ Sigma ^ {2}}{\ displaystyle \ Lambda \ Lambda ^ {*} = \ Sigma ^ {2}} и унитарная матрица комплексных собственных векторов Q {\ displaystyle Q }Q .

Полярное разложение квадратной обратимой вещественной матрицы A {\ displaystyle A}A имеет форму

A = | А | R {\ displaystyle A = | A | R}{\ displaystyle A = | A | R}

, где | А | = (AAT) 1 2 {\ displaystyle | A | = \ left (AA ^ {\textf {T}} \ right) ^ {\ frac {1} {2}}}{\ displaystyle | A | = \ left (AA ^ {\ textf {T}} \ right) ^ {\ frac {1} {2}}} - это положительно определенная матрица и R = | А | - 1 A {\ displaystyle R = | A | ^ {- 1} A}{\ displaystyle R = | A | ^ {- 1} A} - ортогональная матрица.

Построение и доказательства существования

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

Для любого A {\ displaystyle A}A , матрица A ∗ A {\ displaystyle A ^ {*} A}A ^ {*} A эрмитова и положительно полуопределенная, и поэтому унитарно эквивалентна положительная полуопределенная диагональная матрица. Пусть тогда V {\ displaystyle V}V будет унитарным таким, что A ∗ A = VDV ∗ {\ displaystyle A ^ {*} A = VDV ^ {*}}{\ displaystyle A ^ {*} A = VDV ^ {*}} , с диагональю D {\ displaystyle D}D и положительным полуопределенным.

Случай A {\ displaystyle A}A нормальный

Если A {\ displaystyle A}A нормальный, то он унитарно эквивалентен диагональной матрице: A = V Λ V ∗ {\ displaystyle A = V \ Lambda V ^ {*}}{\ displaystyle A = V \ Lambda V ^ {*}} для некоторого унитарного V {\ displaystyle V}V и некоторая диагональная матрица Λ {\ displaystyle \ Lambda}\ Lambda . Тогда мы можем написать

A = V Φ ​​Λ | Λ | V ∗ знак равно (В Φ Λ V ∗) ⏟ ≡ U (V | Λ | V ∗) ⏟ ≡ P, {\ Displaystyle A = V \ Phi _ {\ Lambda} | \ Lambda | V ^ {*} = \ underbrace {\ left (V \ Phi _ {\ Lambda} V ^ {*} \ right)} _ {\ Equiv U} \ underbrace {\ left (V | \ Lambda | V ^ {*} \ right)} _ {\ Equiv P},}{\ displaystyle A = V \ Phi _ {\ Lambda} | \ Lambda | V ^ {*} = \ underbrace {\ left (V \ Phi _ {\ Lambda} V ^ {*} \ right)} _ {\ Equiv U} \ underbrace {\ left (V | \ Lambda | V ^ {*} \ right)} _ {\ Equiv P},} где Φ Λ {\ displaystyle \ Phi _ {\ Lambda}}\ Phi _ {\ Lambda} - диагональная матрица, содержащая фазы элементов Λ {\ displaystyle \ Lambda}\ Lambda , то есть (Φ Λ) ii ≡ Λ ii / | Λ i i | {\ displaystyle (\ Phi _ {\ Lambda}) _ {ii} \ Equiv \ Lambda _ {ii} / | \ Lambda _ {ii} |}{\ displaystyle (\ Phi _ {\ Lambda}) _ {ii} \ Equiv \ Lambda _ {ii} / | \ Lambda _ {ii} |} или (Φ Λ) ii { \ displaystyle (\ Phi _ {\ Lambda}) _ {ii}}{\ displaystyle (\ Phi _ {\ Lambda}) _ {ii}} произвольное комплексное число с единичной величиной, когда Λ ii = 0 {\ displaystyle \ Lambda _ {ii} = 0}{\ displaystyle \ Lambda _ {ii} = 0} .

Таким образом, полярное разложение имеет вид A = UP {\ displaystyle A = UP}A = UP , где U {\ displaystyle U}U и P {\ displaystyle P}P диагональ на основе собственных значений A {\ displaystyle A}A и с собственными значениями, равными фазам и абсолютным значениям собственных значений A {\ displaystyle A}A соответственно.

Случай A {\ displaystyle A}A обратимый

Из разложения по единственному числу можно показать, что A {\ displaystyle A}A обратимо тогда и только тогда, когда A ∗ A {\ displaystyle A ^ {*} A}A ^ {*} A (эквивалентно, AA ∗ {\ displaystyle AA ^ {*}}AA^{*}) есть. Более того, это верно тогда и только тогда, когда все собственные значения A * A {\ displaystyle A ^ {*} A}A ^ {*} A не равны нулю.

В этом случае полярное разложение непосредственно получается записью

A = A (A ∗ A) - 1 2 (A ∗ A) 1 2, {\ displaystyle A = A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}} \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}{\ displaystyle A = A \ left (A ^ {*} A \ справа) ^ {- {\ frac {1} {2}}} \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}

и наблюдая, что A (A ∗ A) - 1 2 {\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}}}{\ displaystyle A \ left (A ^ { *} A \ right) ^ {- {\ frac {1} {2}}}} унитарен. Чтобы убедиться в этом, мы можем использовать спектральное разложение A ∗ A {\ displaystyle A ^ {*} A}A ^ {*} A для записи A (A ∗ A) - 1 2 = AVD - 1 2 V ∗ {\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}} = AVD ^ {- {\ frac {1} {2}}} V ^ {*}}{\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}} = AVD ^ {- {\ frac {1} {2}}} V ^ {*}} .

В этом выражении V ∗ {\ displaystyle V ^ {*}}V ^ {*} унитарен, потому что V {\ displaystyle V}V есть. Чтобы показать, что AVD - 1 2 {\ displaystyle AVD ^ {- {\ frac {1} {2}}}}{\ displaystyle AVD ^ {- {\ frac {1 } {2}}}} является унитарным, мы можем использовать SVD написать A = WD 1 2 V ∗ {\ displaystyle A = WD ^ {\ frac {1} {2}} V ^ {*}}{\ displaystyle A = WD ^ {\ frac {1} {2}} V ^ {*}} , так что

AVD - 1 2 = WD 1 2 V * VD - 1 2 = W, {\ displaystyle AVD ^ {- {\ frac {1} {2}}} = WD ^ {\ frac {1} {2}} V ^ {*} VD ^ {- {\ frac {1} {2}}} = W,}{\ displaystyle AVD ^ {- {\ frac {1} {2}}} = WD ^ {\ frac {1} {2}} V ^ {*} VD ^ {- {\ frac {1} {2}}} = W,}

где снова W {\ displaystyle W}W унитарен по построению.

Еще один способ напрямую показать унитарность A (A ∗ A) - 1 2 {\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac { 1} {2}}}}{\ displaystyle A \ left (A ^ { *} A \ right) ^ {- {\ frac {1} {2}}}} означает, что запись SVD из A {\ displaystyle A}A с точки зрения ранга-1 матрицы как A = ∑ kskvkwk ∗ {\ displaystyle A = \ sum _ {k} s_ {k} v_ {k} w_ {k} ^ {*}}{\ displaystyle A = \ sum _ {k} s_ {k} v_ {k} w_ {k} ^ {*}} , где sk {\ displaystyle s_ {k}}s_ {k} - сингулярные значения A {\ displaystyle A}A , мы имеем

A (A ∗ A) - 1 2 = (∑ j λ jvjwj ∗) (∑ k | λ k | - 1 wkwk ∗) = ∑ k λ k | λ k | vkwk *, {\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}} = \ left (\ sum _ {j} \ lambda _ {j} v_ {j} w_ {j} ^ {*} \ right) \ left (\ sum _ {k} | \ lambda _ {k} | ^ {- 1} w_ {k} w_ {k} ^ {*} \ right) = \ sum _ {k} {\ frac {\ lambda _ {k}} {| \ lambda _ {k} |}} v_ {k} w_ {k} ^ {*},}{\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}} = \ left (\ sum _ {j} \ lambda _ {j} v_ {j} w_ {j} ^ {*} \ right) \ left (\ sum _ {k} | \ lambda _ {k} | ^ {- 1} w_ {k} w_ {k} ^ {*} \ right) = \ sum _ {k} {\ гидроразрыва {\ lambda _ {k}} {| \ lambda _ {k} |}} v_ {k} w_ {k} ^ {*},} что прямо подразумевает унитарность A (A ∗ A) - 1 2 {\ displaystyle A \ left (A ^ {*} A \ right) ^ {- {\ frac {1} {2}}}}{\ displaystyle A \ left (A ^ { *} A \ right) ^ {- {\ frac {1} {2}}}} , потому что матрица унитарна тогда и только тогда, когда ее сингулярные значения имеют унитарное абсолютное значение.

Обратите внимание, как из приведенной выше конструкции следует, что унитарная матрица в полярном разложении обратимой матрицы определена однозначно.

Общий случай

SVD A {\ displaystyle A}A читает A = WD 1 2 V ∗ {\ displaystyle A = WD ^ {\ frac {1} {2}} V ^ {*}}{\ displaystyle A = WD ^ {\ frac {1} {2}} V ^ {*}} , с W, V {\ displaystyle W, V}{\ displaystyle W, V} унитарными матрицами и D {\ displaystyle D}D диагональная положительно полуопределенная матрица. Просто вставив дополнительную пару W {\ displaystyle W}W s или V {\ displaystyle V}V s, мы получим две формы полярного разложения из A {\ displaystyle A}A :

A = WD 1 2 V ∗ = (WD 1 2 W ∗) ⏟ P (WV ∗) ⏟ U = (WV ∗) ⏟ U (VD 1 2 V ∗) ⏟ P ′. {\ Displaystyle A = WD ^ {\ frac {1} {2}} V ^ {*} = \ underbrace {\ left (WD ^ {\ frac {1} {2}} W ^ {*} \ right)} _ {P} \ underbrace {\ left (WV ^ {*} \ right)} _ {U} = \ underbrace {\ left (WV ^ {*} \ right)} _ {U} \ underbrace {\ left (VD ^ {\ frac {1} {2}} V ^ {*} \ right)} _ {P '}.}{\displaystyle A=WD^{\frac {1}{2}}V^{*}=\underbrace {\left(WD^{\frac {1}{2}}W^{*}\right)} _{P}\underbrace {\left(WV^{*}\right)} _{U}=\underbrace {\left(WV^{*}\right)} _{U}\underbrace {\left(VD^{\frac {1}{2}}V^{*}\right)} _{P'}.}
Ограниченные операторы в гильбертовом пространстве

полярное разложение любого линейного ограниченного оператора 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

A ∗ A = (A ∗ A) 1 2 (A ∗ A) 1 2, {\ displaystyle A ^ {*} A = \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}} \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}{\ displaystyle A ^ {*} A = \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}} \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}},}

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

A = U (A ∗ A) 1 2 {\ displaystyle A = U \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}}}{\ displaystyle A = U \ left (A ^ {*} A \ right) ^ {\ frac {1} {2}}}

для некоторой частичной изометрии 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 = U | А | {\ displaystyle A = U | A | \,}A = U | 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-мерной сферы {xi + yj + zk ∈ H: x 2 + y 2 + z 2 = 1} {\ displaystyle \ lbrace xi + yj + zk \ in H: x ^ {2} + y ^ {2} + z ^ {2} = 1 \ rbrace}{\ displaystyle \ lbrace xi + yj + zk \ in H: x ^ {2} + y ^ {2} + z ^ {2} = 1 \ rbrace} из квадратный корень из минус единицы. Для любого r на этой сфере и угла −π < a ≤ π, the versor ear = cos ⁡ (a) + r sin ⁡ (a) {\ displaystyle e ^ {ar} = \ cos (a) + r \ \ sin (a)}e ^ {{ar}} = \ cos (a) + r \ \ sin (a) находится на блоке 3-сфера из H . Для a = 0 и a = π версор равен 1 или -1 независимо от того, какой r выбран. norm t кватерниона q - это евклидово расстояние от начала координат до q. Когда кватернион - это не просто действительное число, тогда существует уникальное полярное разложение q = t e a r. {\ displaystyle q = te ^ {ar}.}{\ displaystyle q = te ^ {ar}.}

Альтернативные плоские разложения

В декартовой плоскости альтернативные плоские кольцевые разложения возникают следующим образом:

  • Если 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 и записывается как
    cosh ⁡ (a) + j sinh ⁡ (a) = exp ⁡ (aj) = eaj {\ displaystyle \ cosh (a) + j \ \ sinh (a) = \ exp (aj) = e ^ {aj}}{\ displaystyle \ cosh (a) + j \ \ sinh (a) = \ exp (aj) = e ^ {aj}}

    где j = +1 и используется арифметика комплексных чисел с разбиением. Ветвь через (−1, 0) отслеживается −e. Поскольку операция умножения на j отражает точку на прямой y = x, вторая гипербола имеет ветви, обозначенные je или −je. Следовательно, точка в одном из квадрантов имеет полярное разложение в одной из форм:

    reaj, - reaj, rjeaj, - rjeaj, r>0 {\ displaystyle re ^ {aj}, - re ^ {aj}, rje ^ {aj}, - rje ^ {aj}, \ quad r>0}{\displaystyle re^{aj},-re^{aj},rje^{aj},-rje^{aj},\quad r>0}
    В наборе {1, -1, j, -j} есть продукты, которые делают его изоморфным четверке Клейна -group. Очевидно, полярное разложение в этом случае включает элемент из этой группы.
Численное определение полярного разложения матрицы

Для вычисления приближения полярного разложения A = UP, обычно унитарного коэффициент U приближается. Итерация основана на методе Герона для квадратного корня из 1 и вычисляет, начиная с U 0 = A {\ displaystyle U_ {0} = A}{\ displaystyle U_ {0} = A} , последовательность

U k + 1 = 1 2 (U k + (U k ∗) - 1), k = 0, 1, 2,… {\ displaystyle U_ {k + 1} = {\ frac {1} {2}} \ left (U_ {k} + \ left (U_ {k} ^ {*} \ right) ^ {- 1} \ right), \ qquad k = 0,1,2, \ ldots}{\ displaystyle U_ {k + 1} = {\ frac {1} {2}} \ left (U_ {k} + \ left (U_ {k} ^ {*} \ right) ^ {- 1} \ right), \ qquad k = 0,1,2, \ ldots}

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

Эта базовая итерация может быть уточнена для ускорения процесса:

  • Каждый шаг или через равные промежутки времени, диапазон сингулярных значений U k {\ displaystyle U_ {k}}U_ {k} оценивается, а затем матрица масштабируется до γ k U k {\ displaystyle \ gamma _ {k} U_ {k}}{\ displaystyle \ gamma _ {k} U_ {k}} , чтобы центрировать сингулярные значения вокруг 1. Масштабирование Фактор γ k {\ displaystyle \ gamma _ {k}}\ gamma _ {k} вычисляется с использованием норм матрицы и обратной матрицы. Примеры таких масштабных оценок:
    γ k = ‖ U k - 1 ‖ 1 ‖ U k - 1 ‖ ∞ ‖ U k ‖ 1 ‖ U k ‖ ∞ 4 {\ displaystyle \ gamma _ {k} = {\ sqrt [{4}] {\ frac {\ left \ | U_ {k} ^ {- 1} \ right \ | _ {1} \ left \ | U_ {k} ^ {- 1} \ right \ | _ { \ infty}} {\ left \ | U_ {k} \ right \ | _ {1} \ left \ | U_ {k} \ right \ | _ {\ infty}}}}}{\ displaystyle \ gamma _ {k} = {\ sqrt [{4}] {\ frac {\ l eft \ | U_ {k} ^ {- 1} \ right \ | _ {1} \ left \ | U_ {k} ^ {- 1} \ right \ | _ {\ infty}} {\ left \ | U_ { k} \ right \ | _ {1} \ left \ | U_ {k} \ right \ | _ {\ infty}}}}}

    с использованием суммы строки и сумма столбцов матричные нормы или

    γ k = ‖ U k - 1 ‖ F ‖ U k ‖ F {\ displaystyle \ gamma _ {k} = {\ sqrt {\ frac {\ left \ | U_ {k} ^ {- 1} \ right \ | _ {F}} {\ left \ | U_ {k} \ right \ | _ {F}}}}}{\ displaystyle \ gamma _ {k} = {\ sqrt {\ frac {\ left \ | U_ {k} ^ {- 1 } \ right \ | _ {F}} {\ left \ | U_ {k} \ right \ | _ {F}}}}}

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

    U k + 1 = 1 2 (γ k U k + 1 γ k (U k ∗) - 1), k = 0, 1, 2,… {\ displaystyle U_ {k + 1} = {\ frac {1} {2}} \ left (\ gamma _ {k} U_ {k} + {\ frac {1} {\ gamma _ {k}}} \ left (U_ { k} ^ {*} \ right) ^ {- 1} \ right), \ qquad k = 0,1,2, \ ldots}{\ displaystyle U_ {k + 1} = {\ frac {1} {2}} \ left (\ gamma _ {k} U_ {k} + {\ frac {1} {\ gamma _ {k}}} \ left (U_ {k} ^ {*} \ right) ^ {- 1} \ right), \ qquad k = 0,1,2, \ ldots}
  • QR-разложение можно использовать на этапе подготовки уменьшить сингулярную матрицу A до регулярной матрицы меньшего размера и внутри каждого шага ускорить вычисление обратной.
  • Метод Герона для вычисления корней x 2 - 1 = 0 {\ displaystyle x ^ {2} -1 = 0}{\ displaystyle x ^ {2} -1 = 0} можно заменить методами более высокого порядка, например, на основе метода Галлея третьего порядка, в результате чего получится
    U k + 1 знак равно U К (I + 3 U К * U К) - 1 (3 I + U К * U К), К = 0, 1, 2,… {\ Displaystyle U_ {k + 1} = U_ {k} \ left (I + 3U_ {k} ^ {*} U_ {k} \ right) ^ {- 1} \ left (3I + U_ {k} ^ {*} U_ {k} \ right), \ qquad k = 0,1,2, \ ldots}{\ displaystyle U_ {k + 1} = U_ {k} \ left (I + 3U_ {k} ^ {*} U_ {k} \ right) ^ {- 1} \ слева (3I + U_ {k} ^ {*} U_ {k} \ right), \ qquad k = 0,1,2, \ ldots}
    Эту итерацию снова можно комбинировать с изменением масштаба. Эта конкретная формула имеет то преимущество, что она также применима к сингулярным или прямоугольным матрицам 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
Последняя правка сделана 2021-06-02 09:36:30
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте