[논문 리뷰] Simple derivation of basic quadrature formulas

Erik Talvila, Matthew Wiersma|arXiv (Cornell University)|2012. 02. 01.

Iterative Methods for Nonlinear Equations참고 문헌 37인용 수 5

한 줄 요약

이 논문은 적절히 선택된 다항식 또는 조각다항식을 두 번째 또는 네 번째 도함수에 곱한 후 부분적분을 적용하여 중점법, 트라페지움 법칙, 그리고 심프슨의 적분 법칙을 단순하고 통합적으로 유도한다. 핵심 기여는 표준 중점법 및 트라페지움 법칙보다 더 작은 오차 상수를 갖는 보정된 트라페지움 법칙을 구성적으로 유도한 것으로, 최소한의 전제 조건으로 계산 및 수치해석 수업에 적합하다.

ABSTRACT

Simple proofs of the midpoint, trapezoidal and Simpson's rules are proved for numerical integration on a compact interval. The integrand is assumed to be twice continuously differentiable for the midpoint and trapezoidal rules, and to be four times continuously differentiable for Simpson's rule. Errors are estimated in terms of the uniform norm of second or fourth derivatives of the integrand. The proof uses only integration by parts, applied to the second or fourth derivative of the integrand, multiplied by an appropriate polynomial or piecewise polynomial function. A corrected trapezoidal rule that includes the first derivative of the integrand at the endpoints of the integration interval is also proved in this manner, the coefficient in the error estimate being smaller than for the midpoint and trapezoidal rules. The proofs are suitable for presentation in a calculus or elementary numerical analysis class. Several student projects are suggested.

연구 동기 및 목표

• 보간법이나 고급 평균값 정리에 의존하지 않고, 고전적인 수치적 적분 법칙—중점법, 트라페지움 법칙, 심프슨 법칙—에 대한 접근성 있고 원시적인 증명을 제공한다.
• 두 번째 또는 네 번째 도함수의 연속성만 가정하고, 부분적분과 원시적 최적화를 통해 이 법칙들의 오차 추정치를 도출한다.
• 끝점의 일阶 도함수를 포함하는 보정된 트라페지움 법칙을 제안하고 철저히 정당화하며, 표준 중점법 및 트라페지움 법칙보다 더 작은 오차 상수를 달성한다.
• 더 이상 스무스하지 않은 함수에 대해서는 헨스톡–쿠르즈와일 적분과 애크시에비치 노름을 사용하여, $ f'' $가 리만 또는 르베그 적분 가능하지 않은 경우에도 오차 추정치를 도출한다.
• 수업 활용을 지원하기 위해 대학 기초 수준과 고급 학부 수준의 학생 프로젝트 및 연습 문제를 제안한다.

제안 방법

• 이 방법은 적분된 함수의 두 번째 또는 네 번째 도함수에 모닉 다항식 또는 조각다항식 함수를 곱한 적분으로 시작하며, 이후 부분적분을 적용하여 적분 공식을 도출한다.
• 중점법과 트라페지움 법칙의 증명은 $ \int_a^b f''(x) p(x) \, dx $에 부분적분을 적용하며, 여기서 $ p(x) $는 경계항을 제거하기 위해 선택된 이차 또는 조각이차 함수이다.
• 심프슨 법칙의 경우, $ \int_a^b f^{(4)}(x) p(x) \, dx $를 사용하며, 모닉이면서 조각다항식인 4차 함수 $ p(x) $를 선택하여 표준 3점 공식을 도출한다.
• 보정된 트라페지움 법칙은 오차 항을 최소화하도록 $ p(x) $를 선택함으로써 도출되며, 이로 인해 오차 bound의 비례 상수는 표준 중점법 및 트라페지움 법칙보다 작아진다.
• 오차 추정치는 두 번째 또는 네 번째 도함수의 균일 노름을 기반으로 유도되며, $ |E^M(f)| \leq \frac{(b-a)^3 \|f''\|_\infty}{24} $ 와 같은 형태의 경계를 갖는다.
• 함수의 정규성에 떨어지는 경우, 헨스톡–쿠르즈와일 적분과 애크시에비치 노름을 사용하여, $ f'' $가 리만 또는 르베그 적분 가능하지 않은 경우에도 오차 추정치를 도출할 수 있다.

실험 결과

연구 질문

• RQ1 어떻게 부분적분과 원시적 다항식 최적화만을 사용하여 중점법, 트라페지움 법칙, 심프슨 법칙을 도출할 수 있는가?
• RQ2 표준 중점법 및 트라페지움 법칙보다 더 작은 오차 상수를 갖는 보정된 트라페지움 법칙을 구성할 수 있는가?
• RQ3 예를 들어 $ f''(x) = \sin(1/x) $ 또는 $ f(x) = \sqrt{x} $와 같이 유계 또는 적분 가능하지 않은 도함수를 갖는 경우, 어떤 오차 경계를 도출할 수 있는가?
• RQ4 함수의 두 번째 도함수가 리만 또는 르베그 적분 가능하지 않지만, 헨스톡–쿠르즈와일 적분 가능할 경우, 이 방법을 어느 정도 일반화할 수 있는가?
• RQ5 두 번째 도함수가 부분구간에서 크게 변할 경우, 오차 추정치를 어떻게 보다 정밀하게 개선할 수 있는가?

주요 결과

• 중점법의 오차는 $ \frac{(b-a)^3 \|f''\|_\infty}{24} $ 이하이며, 복합 법칙의 오차는 $ \leq \frac{(b-a)^3 \|f''\|_\infty}{24n^2} $ 이다.
• 트라페지움 법칙의 오차는 $ \frac{(b-a)^3 \|f''\|_\infty}{12} $ 이하이며, 복합 법칙의 오차는 $ \leq \frac{(b-a)^3 \|f''\|_\infty}{12n^2} $ 이다.
• 심프슨 법칙의 오차는 $ \frac{(b-a)^5 \|f^{(4)}\|_\infty}{2880} $ 이하이며, 복합 법칙의 오차는 $ \leq \frac{(b-a)^5 \|f^{(4)}\|_\infty}{2880n^2} $ 이다.
• 보정된 트라페지움 법칙은 오차 bound가 $ (b-a)^3 \|f''\|_\infty $에 비례하며, 이는 세 가지 표준 법칙 중에서 비례 상수가 가장 작은 것으로 확인된다.
• 예를 들어 $ f(x) = |\log x|^\alpha $ 와 같은 특이점을 갖는 함수의 경우, $ x \to 0^+ $ 일 때 $ f''(x) = O(1/x^2) $ 라면, 이 방법은 유효한 적분 법칙을 도출하며 오차는 $ \sup |f''(x)x^2| $ 이하로 경계된다.
• 만약 $ f'' $가 리만 또는 르베그 적분 가능하지 않더라도, 헨스톡–쿠르즈와일 적분과 애크시에비치 노름을 사용하면 여전히 오차를 추정할 수 있으며, 이는 이 방법의 적용 범위를 더 넓은 함수 클래스로 확장한다.