[논문 리뷰] Calculation of Coefficients of the Optimal Quadrature Formulas in the $W_2^{m,m-1}(0,1)$ Space

Kh. M. Shadimetov, A.R. Hayotov|ArXiv.org|2008. 10. 30.

Material Science and Thermodynamics인용 수 8

한 줄 요약

이 논문은 수치적 적분을 위한 소볼레프 알고리즘을 사용하여 Sobolev 공간 $W_2^{(m,m-1)}(0,1)$ 내에서 최적의 구적 공식을 구성한다. $m=1,2$인 경우, 구적 공식 $\int_0^1 \varphi(x)dx \approx \sum_{\beta=0}^N C_\beta \varphi(x_\beta)$ 내 최적의 계수 $C_\beta$에 대한 명시적 공식을 유도하며, 이는 차수 $m-2$ 이하의 다항함수와 함수 $e^{-x}$에 대해 정확하게 작동한다. $h=1/N$인 간격과 특성 방정식의 근을 이용해 $C_0$ 및 $C_N$에 대한 닫힌 형태의 표현식을 제공한다.

ABSTRACT

In this paper problem of construction of optimal quadrature formulas in $W_2^{(m,m-1)}(0,1)$ space is considered. Here by using Sobolev's algorithm when $m=1,2$ we find the optimal coefficients of the quadrature formulas of the form $$ \int\limits_0^1ϕ(x)dx\cong \sum\limits_{β=0}^NC_βϕ(x_β). $$

연구 동기 및 목표

• 고정된 정밀도로 함수를 수치적으로 적분하기 위한 힐버트 공간 $W_2^{(m,m-1)}(0,1)$ 내 최적의 구적 공식을 구성하는 것.
• 고정된 노드 $x_\beta$에 대해 오차 함수의 노름을 최소화하기 위해 최적의 계수 $C_\beta$를 결정하는 것.
• 모든 차수 $m-2$ 이하의 다항함수와 함수 $e^{-x}$에 대해 구적 공식이 정확하게 작동하도록 보장하는 것.
• 소볼레프 알고리즘을 사용하여 라그랑주 승수와 컨볼루션 기반 공식화로부터 유도된 연립방정식을 해결하는 것.

제안 방법

• 오차 함수 $\ell_N(x)$를 구간 $[0,1]$ 위의 디랙 델타 함수와 지시 함수를 포함하는 분포로 기술한다.
• 최소화할 오차 측도로 이중 공간 $W_2^{(m,m-1)*}(0,1)$ 내 오차 함수의 노름을 사용한다.
• 조건 $(\ell_N, x^\alpha) = 0$ ($\alpha = 0, \dots, m-2$) 및 $(\ell_N, e^{-x}) = 0$ 을 만족시키기 위해 라그랑주 승수를 적용한다.
• 핵함수 $\psi_m(x)$를 포함하는 컨볼루션 방정식으로 시스템을 단순화하며, 계수 $C_\beta$는 균일 메쉬 $x_\beta = h\beta$, $h = 1/N$ 상에서 정의된다.
• 해의 형태를 $C_\beta = h + a_1 \lambda_1^\beta + b_1 \lambda_1^{N-\beta}$로 가정하고, 재귀관계에서 계수를 일치시켜 $a_1$, $b_1$을 결정한다.
• 결정된 매개변수 $a_1$, $b_1$, 및 특성 근 $\lambda_1$를 이용해 첫 번째 및 마지막 계수 $C_0$ 및 $C_N$에 대한 닫힌 형태의 표현식을 도출한다.

실험 결과

연구 질문

• RQ1 고정된 $m=1,2$에 대해 $W_2^{(m,m-1)}(0,1)$ 공간 내 오차 노름을 최소화하는 최적의 계수 $C_\beta$는 무엇인가요?
• RQ2 다항함수의 차수 $m-2$와 함수 $e^{-x}$에 대해 동시에 정확한 구적 공식을 어떻게 확보할 수 있나요?
• RQ3 노드가 균일하게 배치된 경우 최적의 계수에 대한 해석적 구조는 어떻게 되나요?
• RQ4 소볼레프의 알고리즘이 계수에 대한 유도된 연립방정식을 해결하는 데 어떻게 적용될 수 있나요?

주요 결과

• $m=1,2$인 경우, 소볼레프의 알고리즘을 사용하여 최적의 계수 $C_\beta$를 명시적으로 도출하였으며, 이는 $W_2^{(m,m-1)}(0,1)$ 내에서 최소 오차 노름을 보장한다.
• 계수는 선형 재귀를 만족하며, $C_\beta = h + a_1 \lambda_1^\beta + b_1 \lambda_1^{N-\beta}$ 형태로 표현되며, 여기서 $\lambda_1$는 특성 방정식의 근이다.
• 첫 번째 계수는 $C_0 = \frac{e^h - 1 - h}{e^h - 1} - K(h)(\lambda_1^2 + \lambda_1^N - \lambda_1 - \lambda_1^{N+1})$로 표현되며, $K(h)$는 $h$와 $\lambda_1$로 정의된다.
• 마지막 계수는 $C_N = \frac{he^h - e^h + 1}{e^h - 1} - e^h K(h)(\lambda_1^2 + \lambda_1^N - \lambda_1 - \lambda_1^{N+1})$로 표현된다.
• 컨볼루션 방정식에서 $e^{h\beta}$ 및 $h\beta$의 계수를 일치시켜 $a_1$ 및 $b_1$에 대한 선형 시스템을 유도한다.
• 해는 제약 조건에 따라 모든 $x^\alpha$ ($\alpha = 0, \dots, m-2$) 및 $e^{-x}$에 대해 정확하게 작동함을 보장한다.