Сучасна молодь: Основні формули обчислення інтеграла

Формулою Сімпсона називається інтеграл від інтерполяційного многочлена другого степеня на відрізку $[a,b]$ :

${\int\limits_a^b f(x) dx} \approx {\int\limits_{a}^{b} {p_2(x)} dx} = \frac{b-a}{6}{ \left( f(a) + 4 f\left(\frac{a+b}{2}\right) + f(b) \right)},$

де $f(a)$ , $f((a+b)/2)$ і $f(b)$ — значення функції у відповідних точках .

[ред.] Похибка

При умові, що функція $f(x)$ на відрізку $[a,b]$ має похідну четвертого порядку, похибка $E(f)$ , дорівнює:

$E(f) = - \frac{(b-a)^5}{2880}{{f^{(4)}(\zeta)}}, \ \ \ \zeta \in [a,b].$

Зважаючи, що значення $\zeta$ переважно не є відомим, для оцінки похибки використовується нерівність:

$\left| E(f) \right| \le \frac{(b-a)^5}{2880} \max\limits_{x\in[a,b]} {\left| f^{(4)}(x) \right|}.$

[ред.] Виведення формули

Формула Сімпсона може бути виведена за допомогою багатьох різних способів.

[ред.] Квадратична інтерполяція

Якщо замінити функцію $f(x)$ квадратичним поліномом $P(x)$ що приймає ті ж значення що й $f(x)$ у точках a,b і m = (a+b) / 2. використавши інтерполяційну формулу Лагранжа, то одержимо формулу:

$P(x) = f(a) \frac{(x-m)(x-b)}{(a-m)(a-b)} + f(m) \frac{(x-a)(x-b)}{(m-a)(m-b)} + f(b) \frac{(x-a)(x-m)}{(b-a)(b-m)}.$

Після необхідних обчислень одержуємо:

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

[ред.] Використання методів прямокутників і трапецій

У цьому способі виведення використовуються метод прямокутників:

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

і метод трапецій:

$T = \tfrac12 (b-a) (f(a)+f(b)).$

Похибки цих наближень дорівнюють

$-\tfrac1{24} (b-a)^3 f''(a) + O((b-a)^4) \quad$ і $\quad \tfrac1{12} (b-a)^3 f''(a) + O((b-a)^4),$

відповідно. Звідси випливає, що аби позбутися третього степеня слід взяти для наближення величину

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

Однак таким чином одержується формула Сімпсона.

[ред.] Метод невизначених коефіцієнтів

Запишемо в загальному виді:

$\frac{1}{b-a} \int_{a}^{b} f(x) \, dx \approx \alpha f(a) + \beta f\left(\frac{a+b}{2}\right) + \gamma f(b).$

Коефіцієнти α, β і γ можуть бути знайдені з вимоги, що дане наближення є точним для всіх многочленів другого степеня. Таким чином знову ж одержується метод Сімпсона.

[ред.] Ітераційна формула

Для точнішого обчислення інтеграла проміжок $[a,b]$ розбивають на $N$ відрізків однакової довжини і застосовують формулу Сімпсона на кожному з них. Значення інтеграла є сумою для всіх відрізків.

${\int\limits_a^b f(x) dx} \approx \frac h3 \cdot \left( \frac 12 f(x_0)+\sum_{k=1}^{N-1}f(x_k)+2\sum_{k=1}^{N}f \left( \frac{x_{k-1}+x_k}2 \right)+\frac 12 f(x_N) \right)$