Gram-Schmidt Orthogonalization, QR-factorization

/-@,.@-/·2023년 6월 19일
0

Linear Algebra

목록 보기
4/6

Gram-Schmidt Orthogonalization

그람슈미츠는 임의의 벡터들로 orthogonal한 벡터들을 구하는 과정입니다.

임의의 다른 벡터 4개가 존재한다고 가정하겠습니다.(Linear independent, non-orthogonal)

a1,a2,a3,a4a_1, a_2, a_3, a_4 전부 independent하다면 4차원을 전부 span하게 됩니다.

a1a_1,a2a_2 둘 중 하나를 잡고 projection하면 됩니다.
a1a_1v1v_1으로 잡고 a2a_2v1v_1에 proj시킨다면 먼저 v1v_1의 unit-vector를 구해야합니다.

q1=v1v12q_1= {v1\over {\Vert v_1 \Vert}_2}

지금 저희가 구하는 값은 v2v_2입니다. v2v_2v1v_1에 orthogonal한 벡터입니다.

v2=a2(a2Tv1v12)v1v12v_2 = a_2-(a^T_2 \cdot {v1\over {\Vert v_1 \Vert}_2}) \cdot {v1\over {\Vert v_1 \Vert}_2}

a3a_3는 least squares 생각하면 되겠네요!
한 평면 위에 평면 밖으로 튀어나온 벡터(a3a_3)를 직교하는 값을 찾는 것이니 생략하겠습니다.

oorthogonal=직교 90도의각을 이루는 것을 말하고
orthonormal= orthogonal + norm 이 합쳐진 것 -> 서로 직교하며 norm이 1인 것을 의미합니다

QR-factorization

A=QRA = QR

Q = orthogonal matrix
R = upper triangular matrix

[a1,a2,a3][a_1, a_2,a_3]=> rank=3, orthonormalize -> q1,q2,q3q_1,q_2,q_3(gram-schmidt)
a1=a1Tq1q1a_1= a^T_1q_1\cdot q_1
a2=a2Tq1q1+a2Tq2q2a_2= a^T_2q_1\cdot q_1 + a^T_2q_2\cdot q_2
a3=a3Tq1q1+a3Tq2q2+a3Tq3q3a_3= a^T_3q_1\cdot q_1 + a^T_3q_2\cdot q_2 + a^T_3q_3\cdot q_3

[a1a2a3]\begin{bmatrix} a_1 & a_2 & a_3 \end{bmatrix} = [q1q2q3]\begin{bmatrix} q_1 & q_2 & q_3 \end{bmatrix}[a1Tq1a2Tq1a3Tq10a2Tq2a3Tq200a3Tq3]\begin{bmatrix} a^T_1q_1 & a^T_2q_1 & a^T_3q_1\\0 & a^T_2q_2 & a^T_3q_2\\0 & 0 & a^T_3q_3 \end{bmatrix}

빨간색이 Q
보라색이 R

profile
공부한 것과 관심 있는 것을 정리합니다.

0개의 댓글