본문 바로가기
AI 수학/선형대수

Frobenius Norm 최소화 문제

by 깊은대학 2022. 11. 3.

행렬 ARn1×n2, BRn3×n4, YRn1×n4 가 주어졌을 때, 다음과 같은 프로베니우스 놈(Frobenius norm)을 최소화하는 행렬 XRn2×n3 를 구하는 문제를 프로베니우스 놈 최소화 문제라고 한다.

 

(1)Xopt=argminXAXBYF

 

 

 

참고로 어떤 행렬 M 의 프로베니우스 놈 MF 는 다음과 같이 정의된다.

 

(2)MF=tr(MTM)

 

여기서 trtrace 를 의미한다.

 

 

프로베니우스 놈 최소화 문제는 모델식별 문제나 모델축소 문제에서 자주 등장한다. 식 (1)의 일반해는 다음과 같이 주어진다.

 

(3)Xopt=A+YB++ZA+AZBB+

 

여기서 ZRn2×n3 는 임의의 행렬이고 A+ 는 무어-펜로즈(Moore-Penrose ) 유사 역행렬이다.

식 (3)이 식 (1)의 해가 되는 것을 증명하기 위해서 다음과 같은 함수 f(X) 를 정의하고 이 함수가 항상 f(X)0 이 됨을 보이도록 한다.

 

(4)f(X)=AXBYF2AXoptBYF2

 

그러면 AXBYF2AXoptBYF2 이므로 Xopt 가 식 (1)의 해가 되는 것이 증명된다.

프로베니우스 놈의 정의에 따라 식 (4)를 전개하면 다음과 같이 된다.

 

(5)f(X)=tr((AXBY)T(AXBY))         tr((AXoptBY)T(AXoptBY))=2tr(A(XXopt)B(AXoptBY)T)+A(XXopt)BF2

 

여기서 위 식의 첫번째 항을 더 전개하면,

 

(6)tr(A(XXopt) B (AXoptBY)T)     =tr((XXopt) B (AXoptBY)TA)     =tr((XXopt) B ((AXoptB)TYT)A)

 

위 식에 식 (3)을 대입하면 다음과 같이 된다.

 

(7)tr(A(XXopt) B ((A(A+YB++ZA+AZBB+)B)TY)TA)     =tr((XXopt) B ((AA+YB+B+AZBAA+AZBB+B)TYT)A)     =tr((XXopt) B ((AA+YB+B)TYT)A)     =tr((XXopt) (BB+BYTAA+BYTA))     =0

 

따라서 식 (5)는 다음과 같이 된다.

 

(8)f(X)=A(XXopt)BF20

 

따라서 식 (3)이 식 (1)의 해가 되는 것이 증명되었다.

 

 

'AI 수학 > 선형대수' 카테고리의 다른 글

심플렉틱 행렬 (Symplectic Matrix)  (0) 2023.07.01
좌(왼쪽) 고유벡터 (left eigenvector)  (0) 2023.02.10
행렬의 조건수 (Condition Number)  (0) 2021.03.02
놈 (norm)  (0) 2020.10.24
내적 (Inner Product)  (0) 2020.10.21

댓글