심프슨 공식

"오늘의AI위키"는 AI 기술로 일관성 있고 체계적인 최신 지식을 제공하는 혁신 플랫폼입니다.
"오늘의AI위키"의 AI를 통해 더욱 풍부하고 폭넓은 지식 경험을 누리세요.

1. 개요
2. 기본
3. 심프슨 공식의 확장
4. 다른 수치 적분 방법과의 관계
5. 좁은 피크를 갖는 함수의 경우
6. 관련 항목
참조

1. 개요

심프슨 공식은 주어진 구간의 세 점(양 끝점과 중간점)을 지나는 이차 함수를 이용하여 함수를 근사하는 수치 적분 방법이다. 복합 심프슨 공식과 심프슨 3/8 법칙과 같은 확장된 형태가 있으며, 적분 구간을 여러 부분으로 나누어 정확도를 높인다. 심프슨 공식은 중점 규칙과 사다리꼴 공식을 결합하여 유도할 수 있으며, 좁은 피크를 갖는 함수에 대해서는 사다리꼴 공식보다 효율성이 떨어진다. 뉴턴-코츠 공식, 수치 적분, 사다리꼴 공식, 롬베르크 적분법 등이 관련 항목이다.

2. 기본

심프슨 공식은 주어진 구간의 양 끝점과 중간점, 총 세 점을 지나는 이차 함수를 이용하여 원래 함수를 근사하는 방식이다. 이 이차 함수는 라그랑주 다항식 보간법을 통해 구할 수 있다.^[1]

기본 공식은 다음과 같다.

: $\int_{a}^{b} f(x) \, dx\approx \frac{b-a}{6}\left[f(a) + 4f\left(\frac{a+b}{2}\right)+f(b)\right]$

이 공식으로 적분을 구할 때 생기는 오차는 다음과 같다. 여기서 $h = \textstyle \frac{b-a}2$ 이고, $\xi$ 는 ''a''와 ''b'' 사이에 있는 임의의 숫자이다.

: $-\frac{h^5}{90}f^{(4)}(\xi)$

3. 심프슨 공식의 확장

심프슨 공식은 적분 구간이 작을 때 정확하지만, 구간이 넓어지면 오차가 커진다. 이 문제를 해결하기 위해 적분 구간을 여러 개의 작은 부분 구간으로 나누어 각 구간에 심프슨 공식을 적용하고 그 결과들을 합하는 방식을 사용한다. 이를 통해 정확도를 높일 수 있다.^[2]

확장된 공식은 측량학에서 면체적 측량 시 쓰이며 심프슨 1법칙이라고도 부른다.^[2]

이 외에도 겹치는 구간에 심프슨 공식을 적용하는 다른 형태의 복합 심프슨 공식도 있으며, 임의 차수 다항식 보간에 대한 일반화는 뉴턴-코츠 공식으로 이어진다.

3. 1. 복합 심프슨 1/3 법칙 (Composite Simpson's 1/3 rule)

적분 구간

[a, b]

이 작지 않을 때는 먼저 구간을 몇 개의 작은 구간으로 나누고, 각 구간에 심프슨 법칙을 적용해 그 값들을 합해야 한다. 여기서 확장된 공식을 유도할 수 있다. 측량학에서는 면체적 측량 시 쓰이며 심프슨 1법칙이라고 부른다.^[2]

구간

[a, b]

를 짝수

n

개의 부분 구간으로 나누면, 복합 심프슨 공식은 다음과 같다.

:

\int_a^b f(x) \, dx\approx\frac{h}{3}\bigg[f(x_0)+2\sum_{j=1}^{n/2-1}f(x_{2j})+4\sum_{j=1}^{n/2}f(x_{2j-1})+f(x_n)\bigg]

이 식에서

n

은 구간

[a, b]

을 나눈 부분구간의 총 개수를 뜻하며 짝수여야 하고,

h = \textstyle \frac{b-a}n

은 각 부분구간의 길이이다. 면적측량 시 n이 홀수라면 남는 부분은 사다리꼴의 넓이로 계산하여 더해준다.^[2] 이 공식을 정리하면 다음과 같이 쓸 수도 있다.

:

\int_a^b f(x) \, dx\approx\frac{h}{3}\bigg[f(x_0)+4f(x_1)+2f(x_2)+4f(x_3)+...+4f(x_{n-1})+f(x_n)\bigg]

구간

[a, b]

를 길이

h = (b - a)/n

인

n

개의 부분 구간으로 나누고,

0 \leq i \leq n

에 대해 점

x_i = a + ih

을 도입하면 (특히,

x_0 = a

및

x_n = b

), 다음과 같이 나타낼 수 있다.

:

\int_a^b f(x)\, dx\approx \frac{1}{3} h\sum_{i = 1}^{n/2}\big[f(x_{2i - 2}) + 4f(x_{2i - 1}) + f(x_{2i})\big]= \frac{1}{3} h\big[f(x_0) + 4f(x_1) + 2f(x_2) + 4f(x_3) + 2f(x_4) + \dots + 2f(x_{n - 2}) + 4f(x_{n - 1}) + f(x_n)\big]= \frac{1}{3} h\left[f(x_0) + 4\sum_{i = 1}^{n/2} f(x_{2i - 1}) + 2\sum_{i = 1}^{n/2 - 1} f(x_{2i}) + f(x_n)\right].

복합 심프슨 공식에 의해 발생하는 오차는 다음과 같다.

:

-\frac{1}{180} h^4(b - a)f^{(4)}(\xi),

여기서

\xi

는

a

와

b

사이의 어떤 수이고,

h = (b - a)/n

은 "단계 길이"이다. 오차는 (절대값으로) 다음과 같이 제한된다.

:

\frac{1}{180} h^4(b - a) \max_{\xi \in [a, b]} \left|f^{(4)}(\xi)\right|.

3. 2. 심프슨 3/8 법칙 (Simpson's 3/8 rule)

n이 3의 배수일 때 3개의 h씩 묶어 면적을 계산하여 다음 식으로 전체 면적을 구할 수도 있다. n이 3의 배수가 아니면, 2법칙을 적용하고 남는 구간은 심프슨 1법칙으로 계산해서 더한다.^[3]

:

\frac{3}{8}h [f(x_0) + 2\Sigma f(x_{\text{3의  배 수 }}) + 3\Sigma f(x_{\text{남 은  수 }}) + f(x_n)]

토머스 심프슨이 제안한 수치 적분 방법 중 하나로, 3차 보간법을 기반으로 한다. 심프슨 3/8 법칙은 심프슨의 두 번째 공식이라고도 불린다. 공식은 다음과 같다.

\begin{align}\int_a^b f(x)\, dx&\approx \frac{b - a}{8} \left[f(a) + 3f\left(\frac{2a + b}{3}\right) + 3f\left(\frac{a + 2b}{3}\right) + f(b)\right]\\&= \frac{3}{8} h\left[f(a) + 3f\left(a + h\right) + 3f\left(a + 2h\right) + f(b)\right],\end{align}

여기서

h = (b - a)/3

는 단계 크기이다.

이 방법의 오차는 다음과 같다.

-\frac{3}{80} h^5f^{(4)}(\xi) = -\frac{(b - a)^5}{6480} f^{(4)}(\xi),

여기서

\xi

는

a

와

b

사이의 어떤 숫자이다. 따라서 3/8 공식은 표준 방법보다 약 두 배 더 정확하지만, 함수 값을 하나 더 사용한다.

구간

[a, b]

를 길이

h = (b - a)/n

인

n

개의 부분 구간으로 나누고, 점

x_i = a + ih

를

0 \leq i \leq n

에 대해 도입하면(특히,

x_0 = a

와

x_n = b

), 다음과 같이 나타낼 수 있다.

\begin{align}\int_a^b f(x)\, dx&\approx \frac{3}{8} h\sum_{i = 1}^{n/3} \big[f(x_{3i - 3}) + 3f(x_{3i - 2}) + 3f(x_{3i - 1}) + f(x_{3i})\big]\\&= \frac{3}{8} h\big[f(x_0) + 3f(x_1) + 3f(x_2) + 2f(x_3) + 3f(x_4) + 3f(x_5) + 2f(x_6) + \dots \\&\qquad+ 2f(x_{n - 3}) + 3f(x_{n - 2}) + 3f(x_{n - 1}) + f(x_n)\big]\\&= \frac{3}{8} h\left[f(x_0) + 3 \sum_{i = 1,\ 3\nmid i}^{n-1} f(x_i) + 2 \sum_{i=1}^{n/3 - 1} f(x_{3i}) + f(x_n)\right].\end{align}

이 공식에 대한 나머지 항은

-\frac{1}{80} h^4(b - a)f^{(4)}(\xi),

로 표시되지만, 이는

n

이 3의 배수일 때만 사용할 수 있다. 나머지 부분 구간에 대해서는 오차 항의 차수를 변경하지 않고 1/3 공식을 사용할 수 있다. (반대로, 홀수 번호의 부분 구간에 대해서는 3/8 공식과 합성 1/3 공식을 함께 사용할 수 있다.)

이 개념을 더 일반화 하면 뉴턴-코츠 공식을 얻을 수 있다.

3. 3. 불규칙 간격 데이터에 대한 복합 심프슨 공식

어떤 응용 분야에서는 적분 구간

I = [a, b]

를 불균등한 간격으로 나누어야 할 필요가 있다. 이는 데이터의 불균등한 샘플링, 또는 누락되거나 손상된 데이터 포인트 때문일 수 있다.

구간

I

를 너비

h_k

를 갖는 '''짝수 개의

N

개의 부분 구간'''으로 나눈다고 가정하면, 합성 심슨 공식은 다음과 같다.

:

\int_a^b f(x)\, dx \approx \sum_{i = 0}^{N/2 - 1} \frac{h_{2i} + h_{2i + 1}}{6} \left[\left(2 - \frac{h_{2i + 1}}{h_{2i}}\right) f_{2i} + \frac{(h_{2i} + h_{2i + 1})^2}{h_{2i} h_{2i + 1}} f_{2i + 1} + \left(2 - \frac{h_{2i}}{h_{2i + 1}}\right) f_{2i + 2}\right],

여기서

:

f_k = f\left(a + \sum_{i = 0}^{k - 1} h_i\right)

는 구간

I

의

k

번째 샘플링 지점에서의 함수 값이다.

'''홀수 개의

N

개의 부분 구간'''의 경우, 위 공식은 마지막에서 두 번째 구간까지 사용되며, 마지막 구간은 결과에 다음을 더하여 별도로 처리된다.^[1]

:

\alpha f_N + \beta f_{N - 1} - \eta f_{N - 2},

여기서

:

\begin{align}\alpha &= \frac{2h_{N - 1}^2 + 3h_{N - 1} h_{N - 2}}{6(h_{N - 2} + h_{N - 1})},\\[1ex]\beta &= \frac{h_{N - 1}^2 + 3h_{N - 1} h_{N - 2}}{6h_{N - 2}},\\[1ex]\eta &= \frac{h_{N - 1}^3}{6 h_{N - 2}(h_{N - 2} + h_{N - 1})}.\end{align}

4. 다른 수치 적분 방법과의 관계

심프슨 공식은 뉴턴-코츠 공식의 한 예로, 함수를 이차 다항식으로 근사하여 적분값을 계산한다.

심프슨 공식은 중점 규칙과 사다리꼴 공식을 조합하여 유도할 수도 있다. 중점 규칙은 다음과 같다.

: $M = (b - a)f\left(\frac{a + b}{2}\right)$

사다리꼴 공식은 다음과 같다.

: $T = \frac{1}{2} (b - a)\big(f(a) + f(b)\big).$

이 두 공식의 오차는 각각 다음과 같다.

: $\frac{1}{24} (b - a)^3 f''(a) + O\big((b - a)^4\big)$

: $-\frac{1}{12} (b - a)^3 f''(a) + O\big((b - a)^4\big)$

여기서 $O\big((b - a)^4\big)$ 는 $(b - a)^4$ 에 점근적으로 비례하는 항을 나타낸다. 중점 규칙과 사다리꼴 규칙의 오차에서 주요 항을 제거하기 위해 가중 평균을 취하면 심프슨 공식을 얻을 수 있다.

: $\frac{2M + T}{3}.$

롬베르크 적분법은 이와 같이 다른 근사값을 사용하여 오차 항을 제거하는 방법이다.

5. 좁은 피크를 갖는 함수의 경우

좁은 피크 형태 함수 전체 면적 추정 작업에서, 심프슨 공식은 사다리꼴 공식보다 효율성이 떨어진다. 특히, 복합 심프슨 1/3 공식은 사다리꼴 공식과 동일한 정확도를 얻기 위해 1.8배 더 많은 점이 필요하다.^[1] 복합 심프슨 3/8 공식은 정확도가 더 낮다. 심프슨 1/3 공식에 의한 적분은 단계 ''h''를 갖는 사다리꼴 공식에 의한 적분에서 2/3의 값과 단계 2''h''를 갖는 직사각형 공식에 의한 적분에서 1/3의 값이 나오는 가중 평균으로 표현될 수 있다. 정확도는 두 번째 (2''h'' 단계) 항에 의해 결정된다. 적절하게 이동된 프레임을 갖는 심프슨 1/3 공식 복합 합의 평균은 다음과 같은 공식을 생성한다.

: $\int_a^b f(x)\, dx \approx \frac{1}{24} h\left[-f(x_{-1}) + 12f(x_0) + 25f(x_1) + 24\sum_{i = 2}^{n - 2} f(x_i) + 25f(x_{n - 1}) + 12f(x_n) - f(x_{n + 1})\right],$

여기서 적분 영역 외부의 두 점이 사용되고,

: $\int_a^b f(x)\, dx \approx \frac{1}{24} h\left[9f(x_0) + 28f(x_1) + 23f(x_2) + 24\sum_{i = 3}^{n - 3} f(x_i) + 23f(x_{n - 2}) + 28f(x_{n - 1}) + 9f(x_n)\right],$

여기서 적분 영역 내의 점만 사용된다. 두 번째 공식을 3점 영역에 적용하면 1/3 심프슨 공식이 생성되고, 4점 영역에 적용하면 3/8 공식이 생성된다.

이러한 공식은 대안적인 확장된 심프슨 공식과 매우 유사하다. 적분되는 영역의 주요 부분 내의 계수는 가장자리에만 단위가 아닌 계수가 있는 1이다. 이 두 가지 공식은 첫 번째 미분 항을 갖는 오일러-매클로린 공식과 관련될 수 있으며, '''1차''' '''오일러-매클로린 적분 공식'''이라고 명명된다.^[1]

6. 관련 항목

뉴턴-코츠 공식

참조

_[1] 수식 m は“中点”
_[2] 서적 측량학2 형설출판사 2013
_[3] 서적 측량학2 형설출판사 2013

본 사이트는 AI가 위키백과와 뉴스 기사,정부 간행물,학술 논문등을 바탕으로 정보를 가공하여 제공하는 백과사전형 서비스입니다.
모든 문서는 AI에 의해 자동 생성되며, CC BY-SA 4.0 라이선스에 따라 이용할 수 있습니다.
하지만, 위키백과나 뉴스 기사 자체에 오류, 부정확한 정보, 또는 가짜 뉴스가 포함될 수 있으며, AI는 이러한 내용을 완벽하게 걸러내지 못할 수 있습니다.
따라서 제공되는 정보에 일부 오류나 편향이 있을 수 있으므로, 중요한 정보는 반드시 다른 출처를 통해 교차 검증하시기 바랍니다.

문의하기 : help@durumis.com