Cesko.wiki Simpsonovo pravidlo
Author
Albert FloresSimpsonovu metodu lze odvodit aproximací integrandu f (x) (modře) kvadratickou interpolační funkcí P(x) (červeně). Animace ukazující jak Simpsonova metoda aproximuje funkci parabolou, a jak se zmenšuje chyba při zmenšení velikosti kroku Animace ukazující, jak se aproximace Simpsonovou metodou zlepšuje při jemnějším dělení intervalu. Simpsonova pravidla jsou vzorce používané pro numerický výpočet určitých integrálů, které jsou pojmenovány po Thomasu Simpsonovi (1710-1761).
Nejjednodušší Simpsonovo pravidlo nazývané také Simpsonovo 1/3 pravidlo je
: \int_{a}^{b} f(x) \, dx \approx\tfrac{b-a}{6}\left[f(a) + 4f\left(\tfrac{a+b}{2}\right)+f(b)\right].
V němčině a některých dalších jazycích je pojmenované po Janu Keplerovi, který jej odvodil již v roce 1615 pro výpočet objemu vinných sudů . Přibližná rovnost v pravidle přechází na přesnou rovnost, pokud f je polynom nejvýše třetího stupně.
Pokud se interval \langle a, b\rangle rozdělí na n stejných dílů a na každý se aplikuje Simpsonovo pravidlo, dostaneme #Složené Simpsonovo pravidlo|složené Simpsonovo pravidlo, u něhož se pro výpočet integrálu počítá vážený průměr hodnot funkce se střídajícími se vahami 4/3 a 2/3.
Simpsonovo druhé pravidlo, také nazývané #Simpsonovo 3/8 pravidlo|Simpsonovo 3/8 pravidlo vyžaduje o jedno vyhodnocení funkce v integračním rozsahu více, ale poskytuje přesný výsledek pro polynomy do třetího stupně.
1/3 a 3/8 Simpsonova pravidla jsou dvěma speciálními případy uzavřených Newtonových-Cotesových vzorců.
Při konstrukci lodí a při odhadech stability lodí se používá také Simpsonovo třetí pravidlo, které pro numerickou matematiku nemá žádný speciální význam, viz Simpsonova pravidla (stabilita lodí).
Simpsonovo 1/3 pravidlo
Odvození Simpsonova pravidla
Kvadratická interpolace
Simpsonovo pravidlo lze odvodit tak, že nahradíme integrovanou funkci f(x) kvadratickým polynomem P(x) (tj. parabolou), který v koncových bodech integračního intervalu a a b a v jeho středu m = (a+b)/2 nabývá stejných hodnot jako f(x). +more Pro odvození vzorce tohoto polynomu můžeme použít Lagrangeovu polynomiální interpolaci:.
: P(x) = f(a) \tfrac{(x-m)(x-b)}{(a-m)(a-b)} + f(m) \tfrac{(x-a)(x-b)}{(m-a)(m-b)} + f(b) \tfrac{(x-a)(x-m)}{(b-a)(b-m)}. Pomocí substituční metody můžeme ukázat, že
: \int_{a}^{b} P(x) \, dx =\tfrac{b-a}{6}\left[f(a) + 4f\left(\tfrac{a+b}{2}\right)+f(b)\right].
Pokud velikost kroků, ve kterých počítáme hodnotu funkce f(x), vyjádříme jako h=(b-a)/2, dostaneme vzorec, ze kterého je patrné, proč je v názvu této metody numerické integrace zlomek 1/3:
: \int_{a}^{b} P(x) \, dx =\tfrac{h}{3}\left[f(a) + 4f\left(\tfrac{a+b}{2}\right)+f(b)\right]
(zobecnění je popsáno níže).
Průměrování obdélníkové a lichoběžníkové metody
Simpsonovo pravidlo lze odvodit také ze dvou jednodušších aproximací: obdélníkové metody (též Riemannův středový součet):
: M = (b-a) f \left( \tfrac{a+b}{2} \right)
: T = \tfrac12 (b-a) (f(a)+f(b)).
Chyby těchto aproximací jsou
: \tfrac1{24} (b-a)^3 f(a) + O((b-a)^4) \quad\text{a}\quad -\tfrac1{12} (b-a)^3 f(a) + O((b-a)^4),
kde O((b-a)^4) označuje člen asymptoticky úměrný (b-a)^4 . Dva členy O((b-a)^4) si nejsou rovné; pro více detailů viz Landauova notace. +more Z výše uvedeného vzorce pro chyby obdélníkové a lichoběžníkové metody vyplývá, že pokud použijeme vážený průměr výsledků obou metod a bude chtít, aby měl první chybový člen nulovou hodnotu, musíme použít váhu 2/3 pro výsledek obdélníkové metody a 1/3 pro výsledek lichoběžníkové metody:.
: \tfrac{2M+T}{3}.
což je právě Simpsonovo pravidlo.
Myšlenku použití váženého průměru můžeme aplikovat na kombinaci výsledků jiných aproximačních metod výpočtu určitého integrálu (například lichoběžníkové metody s dvojnásobným počtem bodů), a odstranit jiný chybový člen. To je princip Rombergovy metody.
Neurčité koeficienty
Třetí odvození Simpsonova pravidla vychází z myšlenky vyjádřit hodnotu integrálu pomocí kombinace hodnot integrované funkce v krajních bodech integračního intervalu a v jeho středu:
: \tfrac{1}{b-a} \int_{a}^{b} f(x) \, dx \approx \alpha f(a) + \beta f\left(\tfrac{a+b}{2}\right) + \gamma f(b).
Koeficienty α, β a γ lze zafixovat tím, že budeme vyžadovat, aby tato aproximace byla přesná pro všechny kvadratické polynomy. Výsledkem je Simpsonovo pravidlo.
Chyba
Chyba při aproximaci integrálu Simpsonovým pravidlem pro n=2 je
: -\tfrac{1}{90} \left(\tfrac{b-a}{2}\right)^5 f^{(4)}(\xi),
kde \xi (řecké písmeno ksí) je nějaké číslo z intervalu \langle a,\ b\rangle .
Chyba je asymptoticky úměrná (b-a)^5. Podle výše uvedených odvození by však chyba měla být úměrná (b-a)^4. +more Simpsonovo pravidlo však má ve vzorci pro chybu o řád vyšší mocninu, protože body, v nichž se integrand vyčísluje, jsou v intervalu \langle a,\ b\rangle rozložené symetricky.
Když je chybový člen úměrný čtvrté derivaci funkce f v bodě \xi, Simpsonovo pravidlo dává přesné výsledky, pokud je funkce f polynomem nejvýše třetího stupně, protože čtvrtá derivace takových polynomů je ve všech bodech nulová.
Pokud existuje druhá derivace f a v intervalu (a,\ b) je konvexní, platí nerovnost: : (b-a)f\left(\tfrac{a+b}{2}\right) + \tfrac{1}{3}\left(\tfrac{b-a}{2}\right)^3 f\left(\tfrac{a+b}{2}\right) \leq \int_{a}^{b} f(x) \, dx \leq \tfrac{b-a}{6}\left[f(a)+4f\left(\tfrac{a+b}{2}\right)+f(b)\right].
Složené Simpsonovo pravidlo
Pokud interval integrace \langle a, b\rangle je v nějakém smyslu „malý“, pak Simpsonovo pravidlo s n=2 podintervaly bude dávat dostatečně přesnou aproximaci hodnoty integrálu. Slovem „malý“ máme na mysli to, že integrovaná funkce je na tomto intervalu \langle a, b\rangle relativně hladká. +more Pro takové funkce bude hladký kvadratický interpolant, který se používá v Simpsonově pravidle, dávat dobré výsledky.
Často však funkce, kterou chceme integrovat, není na uvažovaném intervalu hladká. Obvykle to znamená, že buď silně osciluje nebo nemá v určitých bodech derivace. +more V těchto případech může Simpsonovo pravidlo dávat velmi špatné výsledky. Obvyklým způsobem, jak se vypořádat s tímto problémem je rozdělení intervalu \langle a, b\rangle na několik menších podintervalů. Simpsonovo pravidlo se pak aplikuje na každý podinterval a výsledky se sečtou, aby se získala aproximace integrálu na celém intervalu. Tento druh přístupu se nazývá složené Simpsonovo pravidlo.
Předpokládejme, že, že interval \langle a, b\rangle je rozdělují na n podintervalů, kde n je sudé číslo. Složené Simpsonovo pravidlo je pak dáno vzorcem
: \begin{align} \int_a^b f(x) \, dx & \approx \frac{h}{3} \sum_{j=1}^{n/2}\bigg[f(x_{2j-2})+4f(x_{2j-1})+f(x_{2j})\bigg]\\ {} & = \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], \end{align}
kde x_j=a+jh pro j=0, 1, ..., n-1, n s h=(b-a)/n; konkrétně, x_0=a a x_n=b. Toto složené pravidlo je pro n = 2 totožné s normálním Simpsonovým pravidlem popsaným v předchozí části.
Chyba složeného Simpsonova pravidla je
:-\frac{h^4}{180}(b-a)f^{(4)}(\xi),
kde \xi je nějaké číslo mezi a a b a h=(b-a)/n je „délka kroku“. Absolutní hodnota chyby je omezena výrazem
:\tfrac{h^4}{180}(b-a)\max_{\xi\in\langle a,b\rangle } |f^{(4)}(\xi)|.
Toto formulace rozděluje interval \langle a,b\rangle na podintervaly stejné délky. V praxi je často výhodné používat podintervaly různých délek, a podrobněji se zabývat místy, kde je integrand méně rozumný. +more To vede k adaptivní Simpsonově metodě.
Simpsonovo 3/8 pravidlo
Simpsonovo 3/8 pravidlo, také nazývaný Simpsonovo druhé pravidlo, je druhá metoda pro numerickou integraci, kterou navrhl Thomas Simpson. Používá kubickou interpolaci místo kvadratické. +more Simpsonovo 3/8 pravidlo popisuje vzorec: : \int_{a}^{b} f(x) \, dx \approx \tfrac{3h}{8}\left[f(a) + 3f\left(\tfrac{2a+b}{3}\right) + 3f\left(\tfrac{a+2b}{3}\right) + f(b)\right] = \tfrac{(b-a)}{8}\left[f(a) + 3f\left(\tfrac{2a+b}{3}\right) + 3f\left(\tfrac{a+2b}{3}\right) + f(b)\right], kde b − a = 3h. Chyba této metody je: :-\tfrac{(b-a)^5}{6480} f^{(4)}(\xi), kde \xi je nějaké číslo mezi a a b. 3/8 pravidlo má tedy dvojnásobnou přesnost než standardní metoda, ale používá o jednu funkční hodnotu více. Existuje také složené 3/8 pravidlo, obdobně jako výše.
Dalším zobecněním tohoto konceptu na interpolaci pomocí polynomů libovolného stupně jsou Newtonovy-Cotesovy vzorce.
Složené Simpsonovo 3/8 pravidlo
Pokud interval \langle a,b\rangle rozdělíme na n podintervalů délky h=(b-a)/n a budeme počítat hodnoty funkce f v bodech x_i = a + ih dostáváme vzorec
: \begin{align} \int_a^b f(x) \, dx & \approx \tfrac{3h}{8}\left[f(x_0) + 3f(x_1) + 3f(x_2) + 2f(x_3) + 3f(x_4) + 3f(x_5) + 2f(x_6) + \cdots + 3f(x_{n-2}) + 3f(x_{n-1}) + f(x_n)\right]. \\ & = \frac{3h}{8} \left[ f(x_0) + 3 \sum_{i \ne 3k}^{n-1} f(x_i) + 2 \sum_{j=1}^{n/3 - 1} f(x_{3j}) + f(x_n) \right] \qquad \text{For: } k \in \mathbb{N}_0 \end{align}
přičemž zbytek je:
:-\frac{h^4}{80}(b-a)f^{(4)}(\xi),
Pro použití tohoto pravidla musí být n dělitelné třemi.
Important
Alternativní rozšíření Simpsonova pravidla
Složené Simpsonovo pravidlo lze také aplikovat na překrývající se segmenty, což dává:
: \int_a^b f(x) \, dx\approx \tfrac{h}{48}\bigg[17f(x_0)+59f(x_1)+43f(x_2)+49f(x_3)+48 \sum_{i=4}^{n-4} f(x_i)+49f(x_{n-3})+43f(x_{n-2})+59f(x_{n-1})+17f(x_n)\bigg].
Tento vzorec lze získat zkombinováním původního složeného Simpsonova pravidla s pravidlem, kde se použije Simpsonovo 3/8 pravidlo v krajních podintervalech a standardní tříbodové pravidlo ve zbývajících podintervalech. Výsledek získáme použitím průměru obou vzorců.
Simpsonova pravidla v případě úzké špičky
Při zjišťování plochy ohraničené funkcí, která má úzké špičky jsou Simpsonova pravidla mnohem méně efektivní než lichoběžníková metoda. Jmenovitě složené Simpsonovo 1/3 pravidlo vyžaduje 1,8krát více bodů, aby se dosáhlo stejné přesnosti jako při použití lichoběžníkové metody. +more Složené Simpsonovo 3/8 pravidlo je v tomto případě ještě méně přesné. Integraci Simpsonovým 1/3 pravidlem lze považovat za součet 2/3 výsledku získaného lichoběžníkovou metodou s krokem h a 1/3 aplikace obdélníkového pravidla s krokem 2h. Zprůměrování složených součtů Simpsonova 1/3 pravidla se správně posunutými rámci dává následující vzorce:.
\int_a^b f(x) \, dx\approx \tfrac{h}{24}\bigg[-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})\bigg]
kde se používají dva body mimo integrační oblast a
\int_a^b f(x) \, dx\approx \tfrac{h}{24}\bigg[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})\bigg]
kde se používají pouze body uvnitř integrační oblasti. Aplikace druhého pravidla na oblast se 3 body dává Simpsonovo 1/3 pravidlo, na oblast se 4 body Simpsonovo 3/8 pravidlo.
Tato pravidla se velmi podobají Pressovu alternativnímu rozšíření Simpsonova pravidla. Koeficienty ve velké části oblasti, kde se počítá integrál, jsou rovny jedné, rozdíly jsou pouze na jejích okrajích. +more Tato tři pravidla lze spojit s Eulerovým-MacLaurinovým vzorcem s první derivací, který se nazývá Eulerova-MacLaurinova integrační pravidla prvního řádu. Dvě výše uvedená pravidla se liší pouze způsobem, jak se počítá první derivace na konci oblasti. Člen s první derivací v Eulerově-MacLaurinově integračním pravidle zodpovídá za integrál druhé derivace, který se rovná rozdílu prvních derivací na okrajích integrační oblasti. Eulerova-Maclaurinova pravidla vyššího řádu je možné generovat přidáním rozdílu třetí, páté, atd. derivace s koeficienty, které jsou definovány Eulerovým-MacLaurinovým vzorcem.
Složené Simpsonovo pravidlo pro nepravidelně vzdálené body
V případech, kdy je potřeba, aby byl integrační interval I=\langle a,b\rangle rozdělen na nestejné intervaly - například kvůli nerovnoměrnému vzorkování dat nebo chybějícím nebo poškozeným datům. Předpokládejme, že rozdělíme interval I na sudé počet N podintervalů se šířkami h_k. +more Pak složené Simpsonovo pravidlo popisuje vztah : \int_a^b f(x)\,\mathrm{d}x = \sum_{i=0}^{N/2-1} \left( \alpha_i f_{2i+2} + \beta_i f_{2i+1} + \eta_i f_{2i} \right),.
kde f_k=f\left(a+\sum_{i=0}^{k-1} h_i \right) jsou funkce hodnoty v k-tém bodě na intervalu I, a koeficienty \alpha_i,\,\beta_i, a \eta_i jsou
: \alpha_i = \frac{2 h_{2i+1}^3-h_{2i}^3+3h_{2i}h_{2i+1}^2}{6h_{2i+1}(h_{2i+1}+h_{2i})},
: \beta_i = \frac{h_{2i+1}^3+h_{2i}^3+3h_{2i+1}h_{2i}(h_{2i+1}+h_{2i})}{6h_{2i+1}h_{2i}}, \text{and}
: \eta_i = \frac{2h_{2i}^3-h_{2i+1}^3+3h_{2i+1}h_{2i}^2}{6 h_{2i}(h_{2i+1}+h_{2i})}.
Pro lichý počet N podintervalů se výše uvedený vzorec používá až na druhý až poslední interval, a poslední interval je ošetřen odděleně přidáním následujícího členu k výsledku: : \alpha f_N + \beta f_{N-1} - \eta f_{N-2},
kde
: \alpha = \frac{2h_{N-1}^2+3h_{N-1}h_{N-2}}{6(h_{N-2}+h_{N-1})},
: \beta = \frac{h_{N-1}^2+3h_{N-1}h_{N-2}}{6h_{N-2}},\,\text{and}
: \eta = \frac{h_{N-1}^3}{6 h_{N-2}(h_{N-2}+h_{N-1})}.
| Ukázka implementace Simpsonovy metody v Pythonu | |
|---|---|
| from collections. abc import Sequence def simpson_nonuniform(x: Sequence[float], f: Sequence[float]) -> float: """ Integrace funkce Simpsonovým pravidlem pro nestajně vzdálené body. +more :param x: Seznam bodů, v nichž jsou známy hodnoty integrované funkce :param f: Seznam hodnot integrované funkce v bodech zadaných v x :return: approximation for the integral Viz ``scipy. integrate. simpson`` a ``_basic_simpson`` pro efektivnější implementaci používající broadcast z modulu numpy. """ N = len(x) - 1 h = [x[i + 1] - x[i] for i in range(0, N)] assert N > 0 result = 0. 0 for i in range(1, N, 2): h0, h1 = h[i - 1], h[i] hph, hdh, hmh = h1 + h0, h1 / h0, h1 * h0 result += (hph / 6) * ( (2 - hdh) * f[i - 1] + (hph**2 / hmh) * f[i] + (2 - 1 / hdh) * f[i + 1] ) if N % 2 == 1: h0, h1 = h[N - 2], h[N - 1] result += f[N] * (2 * h1 ** 2 + 3 * h0 * h1) / (6 * (h0 + h1)) result += f[N - 1] * (h1 ** 2 + 3 * h1 * h0) / (6 * h0) result -= f[N - 2] * h1 ** 3 / (6 * h0 * (h0 + h1)) return result |
Odkazy
Poznámky
Reference
Související články
Externí odkazy
[url=https://encyclopediaofmath. org/index. +morephp. title=Simpson_formula]Simpson formula[/url] na Encyclopedia of Mathematics, EMS Press, 2001 [1994] * [url=https://mathworld. wolfram. com/SimpsonsRule. html]"Simpson's Rule"[/url] MathWorld. Weisstein, Eric W. * [url=http://www. howtoexcel. info/Civil_Engineering/Earthwork_Volume. htm]Aplikace Simpsonova pravidla - Earthwork Excavation[/url] (Note: Vzorec popsaný v toto stránka je správný ale existují chyby v výpočet který musí dávají výsledek 569m3 a not 623m3 jako uvedeny) * [url=http://numericalmethods. eng. usf. edu/topics/simpsons_13rd_rule. html]Simpson's 1/3. pravidlo integrace - Poznámky/zmiňuje, PPT, Mathcad, Matlab, Mathematica, Maple[/url] v [url=http://numericalmethods. eng. usf. edu]Numerický Metody pro KMEN univerzitní[/url] * Podrobný popis počítačové implementace popsal Dorai Sitaram v Teach Yourself Scheme v Fixnum Dnů, [url=http://www. ccs. neu. edu/home/dorai/t-y-schéma/t-y-schéma-Z-H-22. html#node_chap_C]Appendix C[/url].
Kategorie:Integrální počet Kategorie:Numerická integrace Kategorie:Numerická matematika