[논문 리뷰] General inner approximation of vector-valued functions

Mullier Olivier, Éric Goubault외 2명|arXiv (Cornell University)|2013. 10. 07.

Numerical Methods and Algorithms참고 문헌 6인용 수 8

한 줄 요약

이 논문은 도메인과 공역의 차원이 다를 경우에도 벡터 값 함수의 내부 근사화를 확장하는 방법을 제시한다. 자코비안 행렬의 질량과 구간 행렬에서 최대 질량을 갖는 부분행렬을 추출하기 위한 세 가지 기법을 활용함으로써, 자코비안의 질량 r에 따라 함수의 범위에서 최대 r개 성분에 대해 내부 근사화를 계산할 수 있으며, 이는 고차원 설정에서 더 날카운 범위 추정을 가능하게 한다.

ABSTRACT

This paper addresses the problem of evaluating a subset of the range of a vector-valued function. It is based on a work by Gold- sztejn and Jaulin which provides methods based on interval analysis to address this problem when the dimension of the domain and co-domain of the function are equal. This paper extends this result to vector-valued functions with domain and co-domain of different dimensions. This ex- tension requires the knowledge of the rank of the Jacobian function on the whole domain. This leads to the sub-problem of extracting an in- terval sub-matrix of maximum rank from a given interval matrix. Three different techniques leading to approximate solutions of this extraction are proposed and compared.

연구 동기 및 목표

• 기존의 내부 근사화 방법이 도메인과 공역 차원이 동일한 함수에만 적용된다는 한계를 해결하기 위해.
• m ≠ n 인 경우 벡터 값 함수 f: ℝ^m → ℝ^n 의 범위에서 최대 r개 성분에 대한 내부 근사화를 계산할 수 있도록 하기 위해.
• 일般화 과정에서 핵심 문제인 구간 행렬에서 최대 질량 부분행렬을 추출하기 위한 실용적인 기법을 개발하기 위해.
• 로버스트 제어, 글로벌 최적화, 추상 해석 등 응용 분야에서 더 날카운 범위 추정을 지원하는 프레임워크를 제공하기 위해.
• 자코비안의 질량이 일정한 r을 갖는 경우를 다루어, 방법의 타당성과 적용 가능성을 보장하기 위해.

제안 방법

• 자코비안 질량이 일정한 r인 함수 f: ℝ^m → ℝ^n 에 대해 Goldsztejn과 Jaulin의 접근 방식을 일반화한다. 여기서 r ≤ min(m,n) 이다.
• 암시함수 정리를 활용하여 자코비안 부분행렬의 질량이 r인 r개 성분으로 함수를 투영한다.
• 주어진 구간 행렬에서 최대 질량 부분행렬을 추출하기 위한 세 가지 다른 기법을 제안한다: (1) 질량을 드러내는 QR 분해, (2) SVD 기반 선택, (3) 구간 폭과 행렬식 한계를 기반으로 한 휴리스틱 그리디 선택.
• 시스템을 변환하여 역행렬 부분을 분리한 후, 선택된 r개 성분에 대해 구간 뉴턴 반복을 적용한다.
• 정밀도를 제어하는 매개변수 ε 가 존재하며, 더 작은 ε 는 더 날카운 경계를 제공하지만 계산 비용이 증가한다.
• 모든 계산된 내부 근사화가 수학적으로 타당하며 실제로 f에 의해 도달된 값들만 포함되도록 하기 위해 구간 분석에 의존한다.

실험 결과

연구 질문

• RQ1 도메인과 공역 차원이 동일하지 않은 경우를 초월해, 벡터 값 함수의 범위에 대한 내부 근사화를 어떻게 일반화할 수 있는가?
• RQ2 정확한 질량 추출이 NP-난해하므로, 구간 행렬에서 최대 질량 부분행렬을 추출하기 위한 기법은 무엇인가?
• RQ3 자코비안 행렬의 질량을 어떻게 활용하여 어떤 성분이 유효한 내부 근사화를 계산하는 데 사용될 수 있는지 식별할 수 있는가?
• RQ4 다른 부분행렬 추출 기법이 결과 내부 근사화의 날카움과 계산 비용에 어떤 영향을 미치는가?
• RQ5 알고리즘에서 ε 의 선택이 내부 근사화의 정확도와 성능에 어떤 영향을 미치는가?

주요 결과

• 자코비안의 일정한 질량 r에 따라, 도메인과 공역 차원이 다를 경우에도 함수의 범위에서 최대 r개 성분에 대해 내부 근사화를 성공적으로 계산한다.
• 침습 예제(f: ℝ² → ℝ³)에서 ε 가 감소함에 따라 내부 근사화가 점점 날카워지며, 계산 시간은 각각 ε=0.1일 때 0.026초, ε=0.06일 때 0.10초, ε=0.02일 때 0.64초였다.
• 포화 예제(f: ℝ³ → ℝ²)에서는 자코비안 질량이 최대가 아닌 영역이 존재해 더 긴 계산 시간이 소요되었으며, ε=0.1일 때 145.53초, ε=0.05일 때 838.67초였다.
• 세 가지 부분행렬 추출 기법을 비교한 결과, SVD 기반 및 QR 기반 방법이 테스트된 사례에서 휴리스틱 그리디 방법보다 더 우수한 질량 복원 성능을 보였다.
• 자코비안의 질량이 도메인 전역에서 일정한 r을 갖는 함수에 한해 적용 가능하며, 질량이 변하는 경우는 국소 역행렬성의 핵심 가정을 위반하므로 적용 불가능하다.
• 결과는 고차원 공역에서도 의미 있는 내부 근사화를 제공함을 보여주며, 침습 및 포화 예제에서 모두 시각적 검증이 이루어졌다.