[논문 리뷰] On an optimal quadrature formula in Sobolev space $L_2^{(m)} (0,1)$

Kh. M. Shadimetov, A.R. Hayotov외 1명|ArXiv.org|2008. 12. 11.

Differential Equations and Boundary Problems인용 수 4

한 줄 요약

이 논문은 함수 값이 $N+1$개의 노드와 끝점에서의 1차 도함수 값을 사용하는 적분을 위한 $L_2^{(m)}(0,1)$ 소볼레프 공간에서 최적의 구적 공식을 구성한다. 임의의 $N$과 $m \geq 2$에 대해 최적의 계수를 명시적으로 유도하고 오차 기능의 노름을 계산하며, $m=2$ 및 $m=3$일 때 오일러-마클라우린 공식이 최적임을 보여준다. 핵심 결과는 오차 노름의 제곱에 대한 닫힌 표현식이며, $m=2$와 $m=3$에 대해 검증되었으며, 각각 $\|\ell\|^2 = \frac{h^4}{720}$ 및 $\frac{h^6}{30240}$이다.

ABSTRACT

In this paper in the space $L_2^{(m)}(0,1)$ the problem of construction of optimal quadrature formulas is considered. Here the quadrature sum consists on values of integrand at nodes and values of first derivative of integrand at the end points of integration interval. The optimal coefficients are found and norm of the error functional is calculated for arbitrary fixed $N$ and for any $m\geq 2$. It is shown that when $m=2$ and $m=3$ the Euler-Maclaurin quadrature formula is optimal.

연구 동기 및 목표

• 함수 값이 $N+1$개의 등간격 노드와 끝점에서의 1차 도함수 값을 사용하는 $L_2^{(m)}(0,1)$ 소볼레프 공간에서의 최적 구적 공식을 구성하는 것.
• 오차 기능의 노름을 이중 공간 $L_2^{(m)*}(0,1)$에서 최소화하는 계수 $C[\beta]$, $A$, $B$를 결정하는 것.
• 임의의 $N$과 $m \geq 2$에 대해 오차 기능의 제곱 노름에 대한 닫힌 표현식을 유도하여 사르드(Sard)와 소볼레프(Sobolev)의 작업을 일반화하는 것.
• 이 설정에서 오일러-마클라우린 구적 공식이 $m=2$ 및 $m=3$일 때 최적임을 증명하는 것.

제안 방법

• 구적 공식을 $h = 1/N$을 사용하여 $\int_0^1 \varphi(x)dx \approx \sum_{\beta=0}^N C[\beta]\varphi(\beta h) + A\varphi'(0) + B\varphi'(N)$로 설정한다.
• 오차 기능 $\ell(x)$를 델타 함수와 그 도함수를 포함하는 분포로 정의하고, 이중 소볼레프 공간에서의 노름을 계산한다.
• 소볼레프 공간 $L_2^{(m)}(0,1)$에서의 극값 함수 이론과 재생 커널 이론을 활용하여, $\|\ell\|_{L_2^{(m)*}}$를 최소화하는 최적의 계수를 도출한다.
• 오차 노름 표현식을 단순화하기 위해 오일러 다항식, 베르누이 수, 유한 차분의 성질을 적용한다.
• 생성 함수와 잔여물 계산법을 사용하여 일반적인 $\|\ell\|^2$ 공식을 유도하고, 특성 방정식의 근과 베르누이 수를 포함하는 간결한 표현식을 도출한다.
• 일반 결과를 $m=2$ 및 $m=3$에 대해 검증하기 위해 일반 공식에 대입하고 알려진 항등식을 사용해 간단화한다.

실험 결과

연구 질문

• RQ1 함수 값이 $N+1$개의 등간격 노드와 끝점에서의 1차 도함수 값을 사용하는 $L_2^{(m)}(0,1)$에서의 최적 계수 $C[\beta]$, $A$, $B$는 무엇인가?
• RQ2 임의의 $N$과 $m \geq 2$에 대해 이중 공간 $L_2^{(m)*}(0,1)$에서 오차 기능의 제곱 노름에 대한 정확한 표현식은 무엇인가?
• RQ3 오일러-마클라우린 구적 공식은 $m=2$ 및 $m=3$일 때 $L_2^{(m)}(0,1)$에서 최적인가?
• RQ4 베르누이 수와 특성 방정식의 근은 최적 공식에서 오차 노름에 어떻게 기여하는가?
• RQ5 저차수 소볼레프 공간($m=2,3$)에 대해 일반 오차 노름 공식을 단순화하여 명시적인 닫힌 형태를 도출할 수 있는가?

주요 결과

• 모든 $m=2$일 때 최적 계수는 $C[\beta] = h/2$ ($\beta=0,N$), $C[\beta] = h$ ($1 \leq \beta \leq N-1$), $A = h^2/12$, $B = -h^2/12$이며, 이에 따라 $\|\ell\|^2 = h^4/720$이다.
• 모든 $m=3$일 때 동일한 계수 $C[\beta]$, $A$, $B$가 최적이며, $\|\ell\|^2 = h^6/30240$이다.
• 오일러-마클라우린 구적 공식은 주어진 함수 형태에서 $L_2^{(2)}(0,1)$ 및 $L_2^{(3)}(0,1)$에서 최적이다.
• 일반 오차 노름 공식은 베르누이 수 $B_k$, 유한 차분 $\Delta^i 0^k$, 기본 함수에서 유도된 특성 방정식의 근 $q_k$를 포함한다.
• 일반적인 $m$에 대해 오차의 제곱 노름은 $h^{2m}$, $h^{2m+1}$ 및 $q_k$에 의존하는 계수를 포함하는 유한 합으로 표현되며, $\alpha > 1$이 홀수일 때 $B_{\alpha} = 0$이므로 홀수 베르누이 항이 상쇄된다.
• 유도 과정은 계수가 소볼레프 공간에서의 극값 함수 조건에서 유도된 연립방정식을 만족할 때 오차 노름이 최소화됨을 확인한다.