Loading [MathJax]/jax/output/CommonHTML/jax.js

선형머수학

5. 직교행렬과 부분공간

hwijin97 2022. 4. 9. 23:07

 

직교 orthogornal 

  1. 서로 직교인 벡터 x , yxTy=x1y1++xnyn=0 으로 판정한다. x, y 가 복소 성분을 포함하면 ˉxTy=0 으로 판정한다.
  2. 부분 공간의 직교 기저 : 모든 벡터 vTivj=0 을 만족한다. ij
    정규 직교 기저 : vTivi = 1 인 직교 기저, 직교 기저를 정규 직교 기저로 만들려면 각 벡터를 ||vi|| 으로 나눈다.
  3. 직교 부분공간 orthogornal subspace RN :  공간 R 의 모든 벡터는 공간 N 의 모든 벡터는 직교한다. Ax=0 에서 행공간과 영공간은 직교한다. A의 각행과 x 의 곱은 0 이다. 행공간의 기저가 영공간의 기저와 직교하면 모든 벡터는 직교한다.
  4. 정규 직교인 열을 가지는 행렬 Q에서 QTQ=I 를 만족한다. 이때 행렬 Q 는 길고 얇다. ( 행보다 열이 더많다.) 행렬 Q 에 어떤 벡터 x 를 곱하더라도 길이는 변하지 않는다.
  5. 직교 행렬 은 정규직교 열이 있는 정사각행렬 Q 는  QT=Q1 이 성립한다.
    QTQ=I 인 정사각 행렬은 QQT=I 이 성립한다.
    Q 의 좌역행렬 QT 는 우역행렬이다.
    n x n 행렬의 직교인 열들은 Rn 의 정규직교기저가 되고, 행들도 Rn 의 정규직교기저가 된다. 직교행렬은 실제로 정규직교행렬을 의미한다.


직교 벡터

 

판정식 xTy=0 은 피타고라스 정리와 관련이 있다. ||xy||2=||x||2+||y||2

||xy||2=(xy)T(xy)=xTx+yTyxTyyTx 으로 전개되고, xTy, yTx 는 항상 ||x||||y||cosθ 와 동일하다. 따라서 코사인 법칙을 만족하고, 직교하면 두 항이 0 이 된다.

 

직교 기저

 

표준 기저는 Rn 에서 정규 직교한다. 

아다마르 Hadamard 행렬은 R2, R4, R8 , ... 의 직교기저를 포함한다.

H2=[1111],

H4=[H2H2H2H2]

H8=[H4H4H4H4] , ...

위 행렬은 각 열의 길이가 1 이 아니기 때문에, 직교행렬은 아니다. 대신 각 길이로 나누면 간단히 정규직교기저를 찾을 수 있다.

H12 행렬은 존재할까? ( 1 과 -1 로 이루어진 직교기저를 가진 행렬 ) -> 가능!

H6 은? -> 불가능! 왜?

n 이 4의 배수인 것들만 가능함.

Rn 공간의 모든 부분공간은 직교 기저를 가진다.

R3 공간의 평면에 독립인 벡터 a b 에서 직교 기저 c 를 얻기위해, b 에서 a 의 방향요소를 뺀다.

직교 기저 ac , c=baTbaTaa 이고 aTc=0 이다.

이렇게 기저를 직교 기저로 바꾸는 방법이 그람-슈미트 아이디어 이다.

 

직교 부분공간

 

위에서 보았듯 Ax=0 에서 A 의 모든 행은 영공간의 벡터 x 에 곱해지고 그 값은 0 이다. 각 행 ( 각 행의 모든 결합 ) 은 N(A) 에서 x 와 직교한다. 즉 A 의 행공간과 A 의 영공간은 직교한다.

동일하게, ATy=0 에서 A 의 열공간과 AT 의 영공간은 직교한다.

 

여기서 행공간 차원 r 과 영공간 차원 nr 을 보면, 행공간 기저와 열공간 기저는 Rn 의 모든 벡터를 생성한다.

열공간 차원 r 과 좌영공간 차원 mr 도 두 공간의 기저는 Rm 에서 모든 벡터를 생성한다.

 

특이값 분해는 A 행공간에 대한 정규직교기저 v1, , vr 과, A 열공간에 대한 u1, , ur 을 찾는 것이다. 이는 그람-슈미트 방식으로도 가능하다.

특이값 분해로 얻은 특별한 기저는 각 쌍 (v, u) 이 A 에 의해 연결된다는 특징을 추가적으로 지닌다.

특이벡터 : Av1=σ1u1 , Av2=σ2u2, , Avr=σrur

 

정규직교인 열이 있는 길고 얇은 행렬 Q : QTQ=I

 

정규 직교인 행렬 3 x 1, 3 x 2, 3 x 3 행렬 Q1, Q2, Q3 행렬은 모두 QTQ=I 를 만족한다. 따라서 QTQ 의 좌역행렬이다. 마지막 3 x 3 행렬 Q3Q3QT3=I 를 만족해 QT3 Q3 의 는 우역행렬도 된다.

P=QQT 인 모든 행렬은 P2=P 를 만족한다. P2=(QQT)(QQT)=Q(QTQ)QT)=QQT=P

P2=P 는 행렬 P 가 사영행렬임을 의미한다. P2=P=PT 이면, PbP 의 열공간위로의 b 의 직교사영이다.

P 의 열공간은 Q 의 열공간과 동일하다.

 

EX) Q1=13[221], Q2=13[222112] 위로 b=[333] 를 사영해라.

 

P1b=19[221][221][333]=[221] 이 값이 선위로의 정사형이다.

행렬 P1b 를 두 개의 수직성분 e , P1b 으로 나눈다. (사영 P1b 와의 오차 e=(IP1)b )

 

P2b=Q2QT2b=[441] 이 값이 평면 위로의 정사형이다.

여기서 Q2Q1 이 포함되기 때문에, 오차 bP2bbP1b 보다 작다.

 

Q3 에서의 정사형은 전체공간 R3 위로의 정사형이기 때문에, b 그대로이다.

 

 

직교행렬 : Q 는 정사각행렬이다. 따라서 QT=Q1 이다.

 

이러한 Q 는 아주 중요하다. Q 가 2 x 2 이면, 평면의 회전 또는 반사를 의미한다. 

Qrotate=[cosθsinθsinθcosθ],

Qreflect=[cosθsinθsinθcosθ]

각각 행렬은 θ 회전, θ2 선을 경계로 반사

Q1, Q2 는 서로 직교한다. (Q1Q2)T(Q1Q2)=QT2QT1Q1Q2=QT2Q2=I

 


직교 기저 는 Rn 에서 직교 축

 

n x n 직교행렬 Q 의 열 단위벡터 q1qnRn 의 모든 벡터 v=c1q1++cnqn 를 나타낼수 있다.

여기서 ciqii 축위에 있는 v 의 성분이다. 이 말은 i 축위로의 v 의 정사형이라는 의미를 갖는다.

정규직교기저의 계수    c1=qT1v        cn=qTnv 

qTkv=c1qTkq1++ckqTkqk++cnqTkqn=ck, 모든 qTkv=ck 이다.

행렬 방정식 v=Qc 양변에 QT 를 곱하면, QTv=QTQc 가 되고 이를 통해 한번에 정규직교계수를 구할 수 있다.

이것이 직교 기저 적용의 핵심이다. 기저벡터가 직교하면, 각 계수 c1, , cn 을 구할 수 있다.

 


하우스 홀더 반사

 

반사행렬 Q=Hn 의 깔끔한 예로,  단위벡터 u 에 대해서, I2uuT 는 하우스홀더 행렬이 된다.

하우스 홀더의 예    u=[11],    Hn=I2uuT=I2nones(n,n)

Hn 은 대칭행렬이고, H2=HTH=(I2uuT)(I2uuT)=I4uuT+4uuTuuT=I), 직교행렬이다.

직교행렬을 만들기위해 그람-슈미트보다 더간단한 방법으로 자주사용된다.

반사행렬 Hnu=(I2uuT)u=u, u 와 직교하는 벡터 w 에 대해서, Hnw=+w 이다.

이 반사행렬을 포함한 모든 반사행렬은 -1 과 1 을 고유값으로 한다. 여기서 H 의 고유값은 -1 ( 1개 ), 1 ( n-1 ) 개 이다.


웨이블릿 wavelet 행렬

 

직교 기저를 가지는 4x4 행렬을 표현할때, 디지털 신호같이

  • 1 - 1 - 1 - 1
  • 1 - 1 - (-1) - (-1)
  • 1 - (-1) - 0 - 0
  • 0 - 0 - 1 - (-1)

로 표현되어 행렬 W4=[1110111011011101] 로 표현하는 행렬이 하르 웨이블릿 Haar Wavelet 행렬이다. 1 열은 평균, 2 열은 차이, 3열은 작은부분의 차이 등등을 의미할 수 있다.


푸리에 행렬

 

4 x 4 푸리에 행렬 F 에 대해서

F4=[11111ii2i31i2i4i61i3i6i9] 이고, 각 col 은 직교한다.

허수가 포함된 벡터의 직교를 확인하는 법은 내적을 할때, 켤레를 취해야한다. ¯cicj


https://jamboard.google.com/d/13jiMd_aRPe_r5S95hDpdP0JAJCGAe_ms4lpAzdgWn_4

 

로그인 - Google 계정

하나의 계정으로 모든 Google 서비스를 Google 계정으로 로그인

accounts.google.com

 

'선형머수학' 카테고리의 다른 글

7. 대칭인 양의 정부호 행렬 Symmetric Positive Definite  (0) 2022.04.11
6. 고유값과 고유벡터  (0) 2022.04.11
4. 소거법과 A=LU  (0) 2022.04.09
3. 네 가지 기본 부분공간  (0) 2022.04.06
2. 행렬 곱셈 AB  (0) 2022.04.05