Microsoft MVP성태의 닷넷 이야기
Math: 20. Matlab/Octave로 Gram-Schmidt 정규 직교 집합 구하는 방법 [링크 복사], [링크+제목 복사],
조회: 20071
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
 
(연관된 글이 1개 있습니다.)

Matlab/Octave로 Gram-Schmidt 정규 직교 집합 구하는 방법

a 행렬을 정의하고,

a = [1,1,0; 0,1,2; 1,2,1]'
a =

   1   0   1
   1   1   2
   0   2   1

matlab/octave로 QR 분해하면 다음과 같은 결과가 나옵니다.

[q,r] = qr(a)
q =

  -0.70711   0.23570   0.66667
  -0.70711  -0.23570  -0.66667
  -0.00000  -0.94281   0.33333

r =

  -1.41421  -0.70711  -2.12132
   0.00000  -2.12132  -1.17851
   0.00000   0.00000  -0.33333

참고로, wolframalpha 웹 사이트에서 확인하면 실수 계산 값이 아닌 분수 형태의 값으로 구할 수 있습니다. (단지, 부호가 반대입니다.)

gram schmidt {{1,0,1},{1,1,2},{0,2,1}}
; http://m.wolframalpha.com/input/?i=gram+schmidt+%7B%7B1%2C1%2C0%7D%2C%7B0%2C1%2C2%7D%2C%7B1%2C2%2C1%7D%7D&x=0&y=0

1/sqrt(2), -1/(3*sqrt(2)), -2/3
1/sqrt(2), 1/(3*sqrt(2)) , 2/3
0        , 2*sqrt(2) / 3 , -1/3

여기서, Q 결괏값은 직교 행렬로 Xi * Xj = 0 (i != j)인 조건을 만족합니다. 검증은 다음과 같이 쉽게 할 수 있습니다.

# a1 = [-0.70711, -0.70711, 0]
# a2 = [0.23570, -0.23570, -0.94281]
# a3 = [0.66667, -0.66667, 0.33333]

a1 = q(:,1)'
a2 = q(:,2)'
a3 = q(:,3)'

>> dot(a1, a2)
ans = 0
>> dot(a2, a3)
ans = 1.3807e-006  // 실수 계산 값의 오차 누적으로 인한 것일 뿐 0 값임!
>> dot(a3, a1)
ans = 0

또한, Q 결괏값은 정규 직교 행렬이기 때문에 Xi * Xj = 1 (i == j)인 조건을 만족하고, 이는 norm의 값이 1이므로 다음과 같이 검증할 수 있습니다.

>> dot(a1,a1)
ans =  1.0000
>> norm(a1)
ans =  1.0000

// 이하 a2, a3에 대해서도 동일한 결괏값

QR 분해의 Q 행렬은 Gram-Schmidt 정규 직교화 과정에 해당하기 때문에 다음과 같이 수작업으로도 구할 수 있습니다.

a1 = [1,1,0]

na1 = norm(a1)
ea1 = a1 / na1 // 계산 값: 0.70711, 0.70711, 0

a2 = [0,1,2]
_a2 = a2 - dot(a2,ea1) * ea1

n_a2 = norm(_a2)
ea2 = _a2 / n_a2 // 계산 값: -0.23570   0.23570   0.94281

a3 = [1,2,1]
_a3 = a3 - dot(a3,ea1) * ea1 - dot(a3,ea2) * ea2

n_a3 = norm(_a3)
ea3 = _a3 / n_a3 // 계산 값: -0.66667   0.66667  -0.33333

재미있는 것은, 직접 계산한 값의 경우에도 부호는 wolframalpha의 것과 같은 반면 matlab/octave의 QR 함수와는 반대입니다. (혹시 이 이유를 알고 계신 분은 덧글 부탁드립니다. ^^)

실제로 이런 수작업을 별도의 함수로 구현하고 그것이 matlab/octave의 qr 함수와 같다는 것을 다음의 문서에서 정리하고 있습니다.

The Gram-Schmidt process in Matlab
; https://www.math.purdue.edu/~wang838/teaching/GramSchmidt.pdf




기왕 해본 김에 octave 명령어도 알아볼 겸, 정규 직교 집합을 시각화해보겠습니다.

a1 = [1,1,0]'
a2 = [0,1,2]'
a3 = [1,2,1]'

a = [a1 a2 a3]
[q,r] = qr(a)

ad1 = quiver3(0,0,0,-q(1),-q(2),-q(3))
hold on
ad2 = quiver3(0,0,0,-q(4),-q(5),-q(6))
hold on
ad3 = quiver3(0,0,0,-q(7),-q(8),-q(9))
hold on

axis equal

text(-q(1),-q(2),-q(3), 'ad1')
text(-q(4),-q(5),-q(6), 'ad2')
text(-q(7),-q(8),-q(9), 'ad3')

그럼, 다음과 같은 결과를 얻을 수 있고 "Rotate" 기능을 이용해 마우스로 이리저리 돌려 보면 3개의 벡터가 직교하고 있음을 눈으로 쉽게 확인할 수 있습니다.

gramschmidt_1.png




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]

[연관 글]






[최초 등록일: ]
[최종 수정일: 4/6/2023]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer at outlook.com

비밀번호

댓글 작성자
 




... 151  152  153  154  155  [156]  157  158  159  160  161  162  163  164  165  ...
NoWriterDateCnt.TitleFile(s)
1151정성태10/18/201122706Java: 13. 자바도 64비트에서 (2GB) OutOfMemoryException 예외가 발생할까?
1150정성태10/18/201129886.NET Framework: 250. WPF - ComboBox의 SelectionChagned 이벤트파일 다운로드1
1149정성태10/16/201125601.NET Framework: 249. WPF - d:DesignHeight 값을 구할 수 있을까?
1148정성태10/14/201131585Java: 12. 자바에서 LINQ 사용? [7]
1147정성태10/13/201127509.NET Framework: 248. 닷넷에서 지원되는 문자열 인코딩 이름 목록
1146정성태10/12/201133232.NET Framework: 247. LINQ에서의 Max 기능 구현 [10]파일 다운로드1
1144정성태10/10/201128914.NET Framework: 246. WCF - 서버 측에서의 유효한 Timeout 설정파일 다운로드1
1143정성태10/9/201134483.NET Framework: 245. ASP.NET 서버 측 코드에서 페이스북 계정 연동하는 방법
1142정성태10/8/201134983.NET Framework: 244. 윈도우 폼을 열고 닫는 것만으로 메모리 leak이 발생할까? [2]파일 다운로드1
1141정성태10/7/201133610.NET Framework: 243. DataTable에 대해서 Dispose 메서드를 호출할 필요가 있을까? [4]파일 다운로드1
1140정성태10/6/201126775.NET Framework: 242. 닷넷 개발자 입장에서 이해해 보는 자바의 서블릿, JSP
1138정성태10/1/201144441Java: 11. 웹 로직에서 MS-SQL 서버 연결 [2]
1137정성태9/30/201129997Java: 10. 닷넷 개발자가 설치해 본 Oracle WebLogic Server - 설치 및 기본 도메인 구성
1136정성태9/29/201125396개발 환경 구성: 131. Visual Studio - ASP.NET의 Code-behind처럼 cs 파일을 그룹핑하는 매크로 함수 [2]파일 다운로드1
1135정성태9/29/201122852오류 유형: 138. TF10216: Team Foundation services are currently unavailable
1134정성태9/27/201130301.NET Framework: 241. C# 5.0에 새로 추가된 Caller Info 특성 [5]
1133정성태9/25/201133581VC++: 54. C++로 만든 WinRT 프로그램 [2]
1132정성태9/24/201173056Java: 9. 자바의 keytool.exe 사용법과 Tomcat의 SSL 통신 설정
1131정성태9/23/201129261Java: 8. 닷넷 개발자가 구현해 본 자바 웹 서비스 (2)
1130정성태9/23/201137254Java: 7. 닷넷 개발자가 구현해 본 자바 웹 서비스 (1)파일 다운로드2
1129정성태9/22/201128895개발 환경 구성: 130. Hyper-V에 MS-DOS VM 만드는 방법 - MSDN 구독자 대상 [3]
1128정성태9/20/201129053오류 유형: 137. KB2449742 보안 업데이트로 인한 충돌 문제 해결 - 두 번째 이야기
1127정성태9/19/201133144Java: 6. Java에서 MySQL 사용 [2]
1126정성태9/18/201128267Math: 3. "유클리드 호제법"과 "Bezout's identity" 구현 코드(C#)파일 다운로드1
1125정성태9/17/201126170Windows: 54. Windows 8 개발자 Preview를 사용해 보고... [2]
1124정성태9/17/201126423.NET Framework: 240. System.Collections.ArrayList가 .NET 4.5에서 지원이 안된다??? [2]
... 151  152  153  154  155  [156]  157  158  159  160  161  162  163  164  165  ...