[논문 리뷰] Buhmann covariance functions, their compact supports, and their smoothness

Emilio Porcu, P. Zastavyi외 1명|arXiv (Cornell University)|2016. 06. 30.

Advanced Numerical Analysis Techniques참고 문헌 12인용 수 12

한 줄 요약

이 논문은 m차원 유클리드 공간에서 원점에서의 부하만 코바리언스 함수의 스무쓰니스를 높이면서도 양의 정부호성과 컴팩트 지지성을 유지하는 새로운 차분 연산자를 제안한다. 고전적인 몬테 연산자와는 달리, 낮은 차원에서의 양의 정부호성을 잃지 않으며, 공간 통계 및 수치 해석 분야에서 스무쓰하고 컴팩트 지지성을 갖는 라디얼 기저 함수를 구성하는 데 있어 더 우수한 대안을 제공한다.

ABSTRACT

We consider the Buhmann class of compactly supported radial basis functions, whih includes a wealth of special cases that have been studied in both numerical analysis and spatial statistics literatures. In particular, the celebrated Wu, Wendland and Missing Wendland functions are notable special cases of this class. We propose a very simple difference operator and show the conditions for which the application of it to Buhmann functions preserves positive definiteness on $m$-dimensional Euclidean spaces. We also show that the application of the difference operator increases smoothness at the origin, whilst keeping positive definiteness in the same $m$-dimensional Euclidean space, as well as compact support. Thus, our operator is a competitor of the celebrated Mont{é}e operator, which allows to increase the smoothness at the origin, at the expense of losing positive definiteness in the space where the radial basis function is originally defined. The proofs of our results highlight surprising connections with past literatures on celebrated class of functions. Amongst them, absolute and completely monotone functions.

연구 동기 및 목표

• 원점에서 부하만 유형의 컴팩트 지지성 라디얼 기저 함수의 스무쓰니스를 높이되, 양의 정부호성을 유지하는 방법을 개발하는 것.
• 스무쓰니스를 높일 때 낮은 차원에서 양의 정부호성을 상실하는 고전적 몬테 연산자의 한계를 극복하는 것.
• 제안된 연산자와 고전적인 완전 단조 함수 및 로그 단조 함수 클래스 사이의 관계를 설정하는 것.
• 웬랜드, 우, 그리고 누락된 웬랜드 함수와 같은 알려진 특수 케이스를 일반화하는 통합 프레임워크를 제공하는 것.
• 임의의 차원에서 부하만 함수에 대해 제안된 연산자를 적용했을 때 얻을 수 있는 정확한 미분 가능성 수준을 규명하는 것.

제안 방법

• Buhmann 함수를 매개변수 $ k $, $ eta_1 $, $ eta_2 $, $ u $ 로 매개변수화한 간단한 차분 연산자를 제안하며, 이는 이동 및 스케일링된 Buhmann 함수들의 특정 선형 조합을 통해 정의된다.
• Buhmann 함수 $ ilde{ ho}_{ u,k,eta_1,eta_2} $ 에 이 연산자를 적용하여, $ bR^m $ 에서 컴팩트 지지성과 양의 정부호성을 유지하는 새로운 함수 $ f_{ u,k+1,2k+1,eta_1,eta_2} $ 를 도출한다.
• Schoenberg의 정리를 활용하여, 양의 정부호성을 유일한 확률 측도 $ F $ 가 존재함으로써 특징짓는다. 이때 $ ho(t) = bE[bOmega_m(rt)] $ 를 만족하며, $ bOmega_m $ 는 $ bR^m $ 에서 단위 구면 위의 균일 분포를 가진 랜덤 벡터의 특성 함수이다.
• 완전 단조 함수 및 로그 단조 함수 클래스와의 연결 고리를 설정하며, 특히 Berg, Porcu, Mateu 의 보조 함수 $ f_{eta, u}(x) = x^{-eta}(1+x^{ u})^{-1} $ 에 기반한 결과를 활용한다.
• 라플라스 및 푸리에 변환 기법을 활용하여 결과 함수의 스무쓰니스와 양의 정부호성을 분석한다.
• 이 연산자가 원점에서의 미분 가능성의 차수를 두 단계 증가시킴을 입증한다 (예: $ C^{2k} $ 에서 $ C^{2k+2} $ 로). 이와 동시에 $ bR^m $ 에서 컴팩트 지지성과 양의 정부호성을 유지한다.

실험 결과

연구 질문

• RQ1 원점에서의 스무쓰니스를 높이되 컴팩트 지지성과 양의 정부호성을 손상시키지 않는 차분 연산자를 구성할 수 있는가?
• RQ2 제안된 연산자는 다양한 차원에서 양의 정부호성을 유지하는 데 있어 고전적 몬테 연산자와 어떻게 비교되는가?
• RQ3 Buhmann 함수에 차분 연산자를 적용한 후 원점에서 얻을 수 있는 정확한 미분 가능성 수준은 무엇인가?
• RQ4 제안된 연산자와 고전적인 완전 단조 함수 또는 로그 단조 함수 클래스 사이에 어떤 관계가 존재하는가?
• RQ5 제안된 방법을 사용하여 웬랜드, 우, 그리고 누락된 웬랜드 함수와 같은 알려진 특수 케이스를 통합된 프레임워크 내에서 일반화할 수 있는가?

주요 결과

• 제안된 차분 연산자는 원점에서 Buhmann 함수의 스무쓰니스를 두 단계 높여, $ C^{2k} $ 함수를 $ C^{2k+2} $ 함수로 변환한다.
• 결과 함수는 몬테 연산자와 달리 $ bR^m $ 에서도 양의 정부호성을 유지한다. 몬테 연산자는 $ m o bR^m $ 에서 $ bR^{m-2} $ 에서만 양의 정부호성을 유지한다.
• 연산자는 컴팩트 지지성을 유지하므로, 결과 라디얼 기저 함수는 진정으로 밴드형이며 대규모 공간 예측에 적합하다.
• 특히 $ k=0,1,2 $ 일 때, 각각 $ C^2 $, $ C^4 $, $ C^6 $ 에서 $ C^4 $, $ C^6 $, $ C^8 $ 로 미분 가능성 수준이 향상되며, 이는 표 1에 나타나 있다.
• 이 구성은 일반적이며, 웬랜드, 우, 그리고 누락된 웬랜드 함수와 같은 잘 알려진 특수 케이스를 특정 사례로 포함한다.
• 증명 과정은 완전 단조 함수 및 로그 단조 함수 이론과 깊은 연결을 드러내며, 특히 보조 함수 $ f_{eta, u} $ 가 특정 매개변수 제약 조건 하에서 $ bL(0,bR) $ 에 속한다는 것이 입증된다.