Gaussova metoda (numerická integrace)

Gaussova  metoda je numerická integrační metoda , která umožňuje zvýšit algebraické pořadí přesnosti metod založených na interpolačních vzorcích speciální volbou integračních uzlů bez zvýšení počtu použitých hodnot integrandu. Gaussova metoda umožňuje dosáhnout maximální algebraické přesnosti pro daný počet integračních uzlů.

Například pro dva uzly můžete získat metodu přesnosti 3. řádu

,

zatímco pro ekvidistantní uzly metody nad 2. řádem nelze získat. Obecně platí, že pomocí bodů můžete získat metodu s řádem přesnosti . Hodnoty uzlů Gaussovy metody podle bodů jsou kořeny Legendreho polynomu stupně . Hodnoty vah se počítají podle vzorce , kde je první derivace Legendreho polynomu .

Pro uzly a váhy mají následující hodnoty: , váhy : .

(Polynom je definován na segmentu ).

Nejznámější je Gaussova pětibodová metoda.

Gauss-Kronrodova metoda

Nevýhodou Gaussovy metody je, že nemá jednoduchý (z výpočetního hlediska) způsob, jak odhadnout chybu získané hodnoty integrálu. Použití pravidla Runge při dělení integračního segmentu vyžaduje výpočet integrandu na přibližně stejném počtu bodů, přičemž nepřináší téměř žádné zvýšení přesnosti, na rozdíl od jednoduchých metod, kde se přesnost s každým novým rozdělením několikrát zvyšuje. Kronrod navrhl následující metodu pro odhad hodnoty integrálu

,

kde  jsou uzly Gaussovy metody podle bodů a parametry , , jsou zvoleny tak, že řád přesnosti metody je roven . K odhadu chyby pak můžete použít empirický vzorec :

,

kde  je přibližná hodnota integrálu získaného Gaussovou metodou nad body. Knihovny gsl a SLATEC pro výpočet určitých integrálů obsahují rutiny využívající Gauss-Kronrodovu metodu pro 15, 21, 31, 41, 51 a 61 bodů.

Viz také

Literatura

  1. Boltachev G.Sh. Numerické metody v tepelné fyzice. Přednáškový kurz Přednáška 3: Numerická integrace