극성 분해

Polar decomposition

수학에서 정사각형 실행렬 또는 A A 극분해는 A P A=) 인수분해이다.UP U {\U}는 단일 이고 P {\ P 정사각형 및 [1]같은 크기의 양의 반정의 에르미트 입니다.

직관적으로 실제 ×n n\ n \ n \ n n \ displaystyle \ ^ { }의 선형 변환으로 해석하면 극분해는 회전 반사U \ U로 분할한다.{R}n개의 직교 축에 공간의 스케일링.

정사각형 극성 분해는 항상 존재합니다A A 반전 가능한 분해는 고유하며 P(\ P 의 확정입니다.이 경우 A})는 A X A 으로 고유하게 쓸 수 있습니다. U({ U ({ X행렬 P({[2]P})의 고유 자기점 로그입니다.이 분해는 (행렬) Lie [3]그룹의 기본 그룹을 계산하는 데 유용합니다.

극성 분해는 A U A=)로도 정의할 수 있습니다.PU는 P P 대칭 정의 정의이지만 일반적으로 다른 매트릭스이며 U U 위와 같은 매트릭스이다.

행렬의 극분해는 zz 극형식의 행렬 아날로그로 볼 수 서 r 절대값(부정수가 아닌 실수 u 단위규범(ci의 요소)을 갖는 복소수이다.rcle 그룹).

A A=)UP는 × \ U \ { ^ {\ n} 、 P ×n p\ n} C ×n C display C × n \ times n display C display c display C mathbb n \ display 로 확장할 수 .rix. 분해는 항상 하며 P P 항상 고유합니다. U디스플레이 스타일 A A 전체 순위를 가질 에만 고유합니다.[4]

직관적인 해석

실제 m × A(\times x(\ x Ax로 변환하는 m(\displaystyle ^{ 변환으로 해석할 수 있습니다.다음으로 극분해 (\ A에서 계수R (\ R m× (\ mm)의 직교 정규 행렬이다.극분해는 e i(\displaystyle })에 따라 공간 m^{ 축척으로 A A 의해 정의된 선형 변환을 스케일 계수 {})로 표현한 것으로 볼 수 있다.}( P의 작용), 그 \ {m})( R)이 1회 회전 또는 됩니다

분해 A R{\ A A{\ A 정의된 변환({\ R에 이어 특정 직교 방향에 따른 스케일링( P을 나타냅니다.축척 요인은 같지만 방향이 다릅니다.

특성.

의 극분해는 A P .{ { } = \ {U } } { \ overline { U } { \{ P } }

U \U = e θ \ U ^ { \ } a P r = A = r \ A } det 、 gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives 、 Uant gives then gives gives gives = e gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives gives행렬식 1이 있습니다.

양의-반정의 행렬 P는 A가 단수일지라도 항상 유일하며 다음과 같이 표시된다.

서 A { \ A^ { * }는공역 전치(\ A를 나타냅니다.P가 고유하기 때문에 이 표현은 명확하게 정의됩니다.A A{ A^ { * A }는 정의-반정의 에르미트 행렬이며, 따라서 고유한 양의-반정의 에르미트 제곱근을 [5]갖는다는 에 의해 고유성이 보장된다.A가 가역일 경우 P는 양의 유한이며, 따라서 가역 행렬 U는 다음과 같이 고유하게 결정된다.

SVD와의 관계

A A s V of{ \ A \ V { 특이치 분해(SVD)에 관해서는 다음과 같다.

U(\ U V W(\ W 단일 매트릭스입니다(필드가 Rdisplaystyle {인 경우 직교 매트릭스라고 합니다).이를 통해 P P 양의 정의이고 U 단일임을 할 수 있습니다.따라서 SVD의 존재는 극성 분해의 존재와 동등합니다.

A를수도 있습니다(\ A로 할 수도 있습니다.

서 U U 이전과 같으며 P P 다음과 같습니다.

이것을 왼쪽 극성 분해라고 하는 반면, 이전 분해를 오른쪽 극성 분해라고 합니다.좌극 분해는 역극 분해라고도 합니다.

정사각형 가역 극성 다음과 같다.

서 A ( T ) 2{\ A =\ 양의 정의 항목 이고 1 A^{-1는 직교 행렬입니다.

정규 행렬과의 관계

극분해 A A A=)UP는U U P P 출퇴근하는 경우 정상입니다. 또는 이와 동등하게 대각선화할 수 있습니다.

건축 및 존재 증명

극성 분해의 구성 뒤에 있는 핵심 아이디어는 특이값 분해 계산에 사용되는 것과 유사합니다.

의 A A에 대해 A A(\ A)는 에르미트이며 양의 반확정 대각 행렬과 단일적으로 동등합니다.으로 VVD({D}) 대각선 및 양의 A 행렬로 하자

정규 행렬에 대한 파생

A A 정상, (\ A) V(\ A V V V 및 대각행렬 V)에 대해 대각행렬 Adisplaystyle \Lambda V) )의 대각행렬)이 됩니다.s 그러면 쓸 수 있습니다.

여기서 \ \ _ { \ } \ style \ 요소의 위상을 포함하는 대각행렬, 즉 / / I\ Da } 0( (\ _})}( { _}= 경우).

극성 분해는 A A=)입니다. U U P Adisplaystyle A\ A\displaystyle A의 위상 및 절대값과 각각 을 가진다.

가역 행렬에 대한 유도

특이치 분해에서 A A A{ A { A { A { { *} A A ) { * } A {\ A ) { * } ( 동등하게 A { \ AA } } and 만 반전됨을 알 수 있다.또한 이는 AA a A 고유값이 모두 [6]0이 아닌 경우에만 해당된다.

이 경우 극성 분해는 다음과 같은 방법으로 직접 얻을 수 있다.

A( A) - 2 ( \ A \ ( ^ { * * A \ - { \ {2}}}}이 단일임을 관찰합니다.이를 위해 A A 분해를 이용하여 A - 1 - 2 { \ A \ ( A^ { * } \ - { \ { } { } }로 된다.

이 표현에서 V V V{\ V이므로 V displaystyle V^{*}는 유니터리입니다. - 2({ ^ { - { \ {1} {} } })도 유니터리임을 나타내려면 SVD를 사용하여 A D V ({ A=) 라고 쓸 수 있습니다. 다음과 같이 합니다.

서도 W W 구조적으로 통일되어 있습니다.

A A1 A\A\right)^{-{\2}}})의 단일성을 직접적으로 보여주는 또 다른 방법은A(\displaystyleA)의 를 A k)로 표기하는 이다.{\ 단수값입니다

행렬이 단수 절대값을 갖는 경우에만 단일이기 때문에 행렬은 A A) - 2({right2}}}})의 단일성을 직접 나타냅니다.

위의 구조에서 가역행렬의 극분해에서 유니터리행렬이 고유하게 정의되는 방법에 주목하십시오.

일반 파생

A})의 SVD는 A V로 됩니다 A=). W D 대각의 양의 반확정 행렬입니다.W W V V 추가 쌍을 삽입하는 것만으로 A A의 두 가지 극성 분해 형식을 얻을 수 있습니다.

힐베르트 공간의 유계 연산자

복소 힐베르트 공간 사이의 경계 선형 연산자 A의 극성 분해부분 등각학과 비음성 연산자의 산물로서 정규 인수 분해이다.

행렬의 극성 분해는 다음과 같이 일반화된다. A가 유계 선형 연산자일 경우, A의 고유한 인수분해는 A = UP로서 존재하며, 여기서 U는 부분 등각학이고, P는 음이 아닌 자기 점 연산자이고, U의 초기 공간은 P의 범위의 폐쇄이다.

연산자 U는 다음과 같은 문제로 인해 단일이 아닌 부분 등각계로 약화되어야 한다.A가 l(N)의2 단측 이동이면 A = {AA*}1/2 = I입니다.따라서 A = U A이면 U는 A여야 하며, 이는 단일이 아닙니다.

극성 분해의 존재는 더글러스의 보조개념의 결과이다.

Lemma — A, B가 힐베르트 공간 H에서 유계 연산자이고 AA* bb* BB이면 A = CB가 되는 수축 C가 존재합니다.또한 Ker(B*)가 Ker(C)일 경우 C는 고유합니다.

연산자 C는 H모든 h에 대해 C(Bh) := Ah정의할 수 있으며, Ran(B)의 폐쇄에 대한 연속성으로 확장되며, H 모두에 대한 직교 보형에서 0으로 정의할 수 있다.이어서 AA 'BB'가 Ker(B) 'Ker(A)'를* 의미하기 때문* lema가 계속됩니다.

특히.AA = BB* 경우* C는 부분 등각계로 Ker(B*) ) Ker(C)일 경우 고유합니다.일반적으로, 모든 유계 연산자 A에 대해,

여기서 (AA*)1/2는 일반적인 함수 미적분에 의해 주어진 AA의 고유* 양의 제곱근이다.그래서 보조군단에 의해, 우리는
일부 등각 U의 경우 Ker(A*) ker Ker(U)일 경우 고유합니다.P를 (AA*)1/2하면 극분해 A = UP를 얻는다.유사한 인수를 사용하여 A = P'U'나타낼 수 있습니다. 여기서 P'는 양수이고 U'는 부분 등각이다.

H가 유한 차원일 경우, U는 단일 연산자로 확장될 수 있습니다. 일반적으로 이것은 사실이 아닙니다(위의 예 참조).또는 극분해는 특이치 분해 연산자 버전을 사용하여 나타낼 수 있다.

연속함수 미적분의 성질에 의해 A는 A에 의해 생성된 C*대수에 있다.부분 등각법에 대해서도 비슷하지만 약한 문장이 존재한다: U는 A에 의해 생성된 폰 노이만 대수이다.A가 반전 가능한 경우 극성 부품 U도 C* 대수에 포함됩니다.

무제한 연산자

A가 복잡한 힐버트 공간 사이의 닫힌, 조밀하게 정의된 무한 연산자일 경우, 여전히 (독특한) 극성 분해를 갖는다.

여기서 A는 A와 같은 도메인을 가진 (아마도 무제한인) 비음성 자기 인접 연산자이고, U는 Ran(A) 범위의 직교 상보에서 소실되는 부분 등각계이다.

증명은 위와 동일한 보조항목을 사용하며, 이는 일반적으로 무제한 연산자에 적용됩니다.모든 h µ 돔(AA*)에 대해 돔(AA*) = (BB*) 및 AAh* = BBh이면*, A = UB되는 부분 등각 U가 존재한다.Ran(B) ker Ker(U)일 경우 U는 고유합니다.연산자 A는 폐쇄적이고 조밀하게 정의되어 연산자* AA가 (밀도의 도메인과) 자기접속되므로 (AA*)1/2를 정의할 수 있다.보조제를 바르는 것은 극성 분해를 일으킨다.

무한 연산자 A가 폰 노이만 대수 M에 관련되고, A = UP가 극분해라면, UM에 있고, [0, δ]의 보렐 집합 B에 대한B P, 1(P)의 스펙트럼 투영도 마찬가지이다.

사분극 분해

4원소 H의 극성 분해는2차원 구{ i + j + H : 2 + + }({ xi H}=)의 제곱근 빼기 에 의존합니다.이 구면상의 임의의 r 및 각도 < θ θ가 주어졌을 때, e cos () + θ ( }=\ H의 단위 3구면상있다. A = 0 또는 a = θ, 1에 관계없이, 또는 1이다.사분위수 q의 노름 t는 원점에서 q까지의 유클리드 거리이다.4분의 1이 단순히 실수가 아닌 경우, 고유한 극성 a. \ q= te}가 존재합니다.

대체 평면 분해

데카르트 평면에서 다음과 같은 대체 평면 링 분해가 발생합니다.

  • x 0 0일 경우 z = x(1 + y(y/x))는 이중수 z = x + 의 극성 분해이다. 여기서 θ2 = 0, 즉 θ0이다.이 극성 분해에서 단위 원은 x = 1, 극각은 기울기 y/x로 대체되었으며 반지름 x는 왼쪽 반평면에서 음수입니다.
  • x2 y2 y이면 단위 쌍곡선2 - y = 1과 그 켤레2 x - y22 = - 1을 사용하여 단위 쌍곡선을 통과하는 분기에 기초한 극성 분해를 형성할 수 있다.이 분지는 쌍곡선 앵글에 의해 매개 변수화되며 다음과 같이 쓰여집니다.

    여기2 j = +1분할 복소수 산술이[7] 사용됩니다.(-1, 0)까지의 브랜치는 -e로aj 추적됩니다.j를 곱하는 연산은 y = x 의 점을 반영하므로, 두 번째 쌍곡선은 je 또는 -jeaj 추적되는aj 분기를 가진다.따라서 사분원 중 하나의 점은 다음 중 하나의 형태로 극성 분해가 됩니다.

    세트 {1, -1, j, -j}에는 클라인 4그룹과 동형인 제품이 있습니다.분명히 이 경우 극성 분해는 그 그룹의 원소를 포함합니다.

매트릭스 극성 분해 수치 측정

극성 분해 A = UP의 근사치를 계산하기 위해 일반적으로 단일 인자 U는 [8][9]근사치입니다.반복은 1의 제곱근에 대한 Heron의 방법을 기반으로 하며 0 {\}=부터 시작하여 계산됩니다.

반전 및 헤르미트 활용의 조합은 특이값 분해에서 단일 인자가 동일하게 유지되고 반복이 특이값에 대한 헤론의 방법으로 감소하도록 선택된다.

이 기본적인 반복은 프로세스를 고속화하기 위해 조정할 수 있습니다.

  • 매 단계 또는 정기적인 으로Uk(\k})의 특이값 범위를 추정하여 행렬의 크기를 _})로 조정합니다. 단수값의 중심을 1에 맞춥니다.스케일링 계수 k _ 매트릭스와 그 역행렬의 매트릭스 규범을 사용하여 계산한다.이러한 규모 추정치의 예는 다음과 같다.

    행-합 및 열-합 매트릭스 규범을 사용하거나
    프로베니우스 규범을 사용해서요스케일 팩터를 포함하여 반복은 다음과 같습니다.

  • QR 분해는 준비 단계에서 단일 행렬 A를 더 작은 정규 행렬로 축소하고 모든 단계에서 역행렬의 계산을 가속화하기 위해 사용할 수 있습니다.
  • 2 - 루트를 계산하는 헤론의 은 예를 들어 핼리의 3차 방법한 고차법으로 대체될 수 있으며, 결과적으로 다음과 같은 결과가 초래됩니다
    이 반복은 다시 재스케일링과 결합할 수 있습니다.이 특정 공식은 단수 행렬 또는 직사각형 행렬 A에도 적용할 수 있다는 장점이 있다.

「 」를 참조해 주세요.

레퍼런스

  1. ^ Hall 2015 섹션 2.5
  2. ^ 2015 정리 2.17
  3. ^ Hall 2015 섹션 13.3
  4. ^ Higham, Nicholas J.; Schreiber, Robert S. (1990). "Fast polar decomposition of an arbitrary matrix". SIAM J. Sci. Stat. Comput. Philadelphia, PA, USA: Society for Industrial and Applied Mathematics. 11 (4): 648–655. CiteSeerX 10.1.1.111.9239. doi:10.1137/0911038. ISSN 0196-5204.
  5. ^ 2015 레마 2.18
  6. ^ A A {\ AA의 양성에 의해 고유값이 모두 실재하고 엄밀하게 양수임을 알 수 있습니다.
  7. ^ Sobczyk, G.(1995) "초과기수 평면", 대학수학저널 26:268-80
  8. ^ Higham, Nicholas J. (1986). "Computing the polar decomposition with applications". SIAM J. Sci. Stat. Comput. Philadelphia, PA, USA: Society for Industrial and Applied Mathematics. 7 (4): 1160–1174. CiteSeerX 10.1.1.137.7354. doi:10.1137/0907079. ISSN 0196-5204.
  9. ^ Byers, Ralph; Hongguo Xu (2008). "A New Scaling for Newton's Iteration for the Polar Decomposition and its Backward Stability". SIAM J. Matrix Anal. Appl. Philadelphia, PA, USA: Society for Industrial and Applied Mathematics. 30 (2): 822–843. CiteSeerX 10.1.1.378.6737. doi:10.1137/070699895. ISSN 0895-4798.