Metoda obdélníků je metoda numerické integrace funkce jedné proměnné, která spočívá v nahrazení integrandu polynomem nultého stupně, tedy konstantou, na každém elementárním segmentu. Pokud budeme uvažovat graf integrandu, pak bude metoda spočívat v přibližném výpočtu plochy pod grafem sečtením ploch konečného počtu obdélníků, jejichž šířka bude určena vzdáleností mezi odpovídajícími sousedními integracemi. uzly a výšku podle hodnoty integrandu v těchto uzlech. Algebraický řád přesnosti je 0. (Pro vzorec středních obdélníků je to 1).
Pokud je segment elementární a nepodléhá dalšímu dělení, lze hodnotu integrálu zjistit z
V případě rozdělení integračního segmentu na elementární segmenty jsou výše uvedené vzorce aplikovány na každý z těchto elementárních segmentů mezi dvěma sousedními uzly. Výsledkem jsou složené kvadraturní vzorce
Vzorec s výpočtem hodnoty ve středu mezi dvěma uzly lze použít pouze tehdy, když je integrand specifikován analyticky nebo jiným způsobem, který umožňuje výpočet hodnoty v libovolném bodě. V úlohách, kde je funkce dána tabulkou hodnot, zbývá pouze vypočítat průměrnou hodnotu mezi integrály vypočítanou vzorcem levého a pravého obdélníku, což vede ke složenému kvadraturnímu lichoběžníkovému vzorci .
Protože složené kvadraturní vzorce nejsou nic jiného než součty zahrnuté v definici Riemannova integrálu , konvergují k přesné hodnotě integrálu. V souladu s tím se s rostoucí přesností zvyšuje výsledek získaný pomocí přibližných vzorců.
Jednotnou mřížku lze popsat pomocí následující sady vzorců:
kde je mřížkový krok.
Pro jednotné mřížky lze vzorce obdélníku zapsat jako následující vzorce Cotes :
U vzorců pravého a levého obdélníku je chyba
Pro vzorec obdélníků (střední)
Pro složené vzorce pravého a levého obdélníku na jednotné mřížce:
Pro složený vzorec obdélníků:
Vzorec středních obdélníků pro analyticky danou funkci napsaný v C
double InFunction ( double x ) { //Integrovaná funkce vrátit hřích ( x ); } double CalcIntegral ( double a , double b , int n ) { dvojitý výsledek = 0 , h = ( b - a ) / n ; for ( int i = 0 ; i < n ; i ++ ) { výsledek += InFunkce ( a + h / 2 + i * h ); } výsledek *= h ; vrátit výsledek ; }Integrální počet | ||
---|---|---|
Hlavní | ||
Zobecnění Riemannova integrálu | ||
Integrální transformace |
| |
Numerická integrace | ||
teorie míry | ||
související témata | ||
Seznamy integrálů |