[논문 리뷰] High-dimensional periodic sampling on Smolyak grids based on B-spline quasi-interpolation

Ðinh Dũng|arXiv (Cornell University)|2015. 02. 05.

Mathematical Approximation and Integration참고 문헌 26인용 수 3

한 줄 요약

이 논문은 B-spline 쿼asi-보간을 사용하여 고차원 주기적 함수의 혼합 스무스니스에 대해 선형 샘플링 복구 및 적분 공식을 Smolyak 격자에서 제안한다. 이는 차원 $d$, 정밀도 수준 $m$, 활성 변수 수 $u$에 대해 명시적인 오차 한계를 확립하며, 희소 Smolyak 격자가 훨씬 적은 노드 수로도 전체 격자와 동일한 수렴 속도를 달성할 수 있음을 보여준다.

ABSTRACT

We constructed linear algorithms of sampling recovery and cubature formulas on Smolyak grids parametrized by $m \in \mathbb{N}$ of periodic $d$-variate functions having Lipschitz-Hölder mixed smoothness $α> 0$ based on B-spline quasi-interpolation, and studied their optimality. We established lower estimates (for $α\le 2$) and upper bounds of the error of the optimal sampling recovery and the optimal integration on Smolyak grids, explicit in $d$, $m$ and the number $ν$ of active variables of functions when $d$ and $m$ may be large.

연구 동기 및 목표

• 혼합 스무스니스 $\alpha>0$를 가진 주기적 $d$변수 함수에 대해 Smolyak 격자에서 선형 샘플링 복구 및 적분 알고리즘을 개발한다.
• 이러한 알고리즘의 최적성은 $d$, $m$, 및 $\nu$에 대해 명시적인 오차 한계의 관점에서 분석한다.
• 스무스니스 $\alpha \leq 2$에 대해 하한을 확립하고, Smolyak 격자에서 최적의 샘플링 복구 및 적분의 오차에 상한을 설정한다.
• 희소 Smolyak 격자가 전체 격자와 동일한 수렴 속도를 달성하면서도 훨씬 적은 샘플링 포인트를 사용함을 보여준다.

제안 방법

• B-spline 쿼اسي-보간 연산자를 사용하여 Smolyak 격자 $G^d(m)$에서 선형 샘플링 복구 및 적분 공식을 구성한다.
• 일변수 대칭 라우렌트 다항식 $P_\Lambda$를 사용하여 근사 오차를 명시적으로 제어할 수 있는 쿼اسي-보간 연산자 $Q$를 정의한다.
• 혼합 차분 연산자 $\Delta_h^{r,u}$를 사용하여 공간 $H^\alpha_\infty(\mathbb{T}^d)$에서 리프시츠-홀더 반노름을 정의한다.
• $\|P_{\text{even}}^*\|$, $\|P_{\text{odd}}^*\|$, 및 $\|L_\Lambda\|_{C(\mathbb{T})}$의 분석을 통해 오차 한계를 유도하며, 이는 안정성과 근사 품질을 제어한다.
• 조각별 선형 B-spline($r=1$)을 통한 페벨 급수 표현을 적용하고, 이는 삼차($r=2$) 및 오차($r=3$) B-spline 쿼اسي-보간으로 확장한다.
• 계층적 분해와 $|k|_1 = m$ 조건을 사용한 희소 격자 구축을 통해 희소성을 확보하면서도 최적의 수렴 속도를 유지한다.

실험 결과

연구 질문

• RQ1 B-spline 쿼اسي-보간이 혼합 스무스니스 $\alpha>0$를 가진 고차원 주기적 함수에 대해 Smolyak 격자에서 최적의 수렴 속도를 달성할 수 있는가?
• RQ2 Smolyak 격자에서 샘플링 복구 및 적분의 명시적 오차 한계는 $d$, $m$, 및 $\nu$에 대해 어떻게 표현되는가?
• RQ3 고차원에서 B-spline 쿼اسي-보간의 근사 성질은 표준 전체 격자와 비교해 어떻게 되는가?
• RQ4 라우렌트 다항식 $P_\Lambda$는 쿼اسي-보간 연산자의 안정성과 정확도를 제어하는 데 어떤 역할을 하는가?
• RQ5 스무스니스 $\alpha \leq 2$에 대해 오차 한계가 날카로운가? 그리고 하한과 비교해보면 어떻게 되는가?

주요 결과

• Smolyak 격자에서 최적의 샘플링 복구 오차는 $C \cdot 2^{-m\alpha}$로 유계이며, 이때 $C$는 $d$, $\nu$, 및 $\alpha$에 대해 명시적으로 의존한다. 이 한계는 $\alpha \leq 2$일 때 최적이며, 이는 최적성 보장한다.
• 조각별 선형 B-spline 쿼اسي-보간($r=1$)의 경우, $\|P_{\text{even}}^*\| = 1/2$, $\|P_{\text{odd}}^*\| = 0$, 및 $\|L_\Lambda\| = 1$이며, 이는 안정성과 최적의 수렴 속도 보장한다.
• 삼차 B-spline 쿼اسي-보간($r=2$)의 경우, $\|P_{\text{even}}^*\| = 3/8$, $\|P_{\text{odd}}^*\| = 1/2$, 및 $\|L_\Lambda\| = 11/9$이며, 이는 매우 높은 근사 정확도를 제공한다.
• 오차 B-spline 쿼اسي-보간($r=3$)의 경우, $\|P_{\text{even}}^*\| \approx 0.97$, $\|P_{\text{odd}}^*\| \approx 1.00$, 및 $\|P_\Lambda\| \approx 2.63$이며, 이는 거의 최적의 안정성과 고차 정확도를 보여준다.
• 격자 $G^d(m)$의 격자 포인트 수는 $\frac{1}{(d-1)!}2^m m^{d-1}$보다 작으며, 이는 $2^{dm}$보다 크게 감소한 것으로, Smolyak 격자의 희소성을 확인한다.
• 제안된 알고리즘은 전체 격자에서의 알고리즘과 동일한 수렴 속도를 달성하며, 이는 희소 격지가 혼합 스무스니스를 가진 고차원 주기적 함수에 대해 최적이며, 최적의 성능을 보임을 증명한다.