Razlika između verzija stranice "Numerička integracija"
[pregledana izmjena] | [pregledana izmjena] |
Uklonjeni sadržaj Dodani sadržaj
mNo edit summary |
m čišćenje, replaced: ) je → ) jeste, replaced: → (7) using AWB |
||
Red 104:
</source>
Neki detalji algoritma zahtjevaju pažljiv pristup. Za više slučajeva, pronalazak greške kvadrature na intervalu za funkciju ''f''(''x'') nije očita. Jedna poznata solucija je da se koriste dva pravila kvadratura i korištenjem razlike se može procijeniti greška kvadrature. Čest problem je odlučiti šta označava "premalo" ili "preveliko". Lokalni kriterij za "preveliko" je da kvadraturna greška ne smije biti veća od ''t''
=== Metoda ekstrapolacije ===
Red 112:
=== Konzervativna (a priori) procjena greške ===
Ako se pusti da ''f'' ima ograničen prvi izvod na zatvorenom intervalu [''a'',''b'']. Teorema srednje vrijednosti za ''f'', gdje je ''x''
: <math>(x - a) f'(v_x) = f(x) - f(a)\,</math>
Red 125:
: <math>\left| \int_a^b f(x)\,dx - (b - a) f(a) \right| \leq {(b - a)^2 \over 2} \sup_{a \leq x \leq b} \left| f'(x) \right|</math> (**)
Ako aproksimiramo integral '''∫<sub>''a''</sub><sup>''b''</sup>
: <math>{n^{-1} \over 2} \sup_{0 \leq x \leq 1} \left| f'(x) \right|</math>
Red 197:
* [[QUADPACK]] (dio SLATEC-a): opis [http://www.netlib.org/slatec/src/qpdoc.f], izvorni kod [http://www.netlib.org/slatec/src]. QUADPACK je kolekcija algoritama, u [[FORTRAN|Fortran]]<nowiki/>u, za numeričku integraciju baziranu na Gausovim kvadraturama.
* [http://openopt.org/interalg interalg]: softver od [[OpenOpt]]/[[FuncDesigner]] frameworksa, baziran na analizi intervala, garantovanoj preciznosti, licenci: BSD (besplatan za svaku svrhu)
* [http://www.gnu.org/software/gsl/ GSL]: The GNU Scientific Library (GSL)
* Algoritmi numeričke integracije se mogu naći u [[Guide to Available Mathematical Software|GAMS]] razredu [http://gams.nist.gov/serve.cgi/Class/H2 H2].
* [http://www.alglib.net/integral/ ALGLIB] je kolekcija algoritama napisanih u C# / C++ / Delphi / Visual Basic / itd. za numeričku integraciju (uključuje Bulirsch-Stoer i Runge-Kutta integratore).
|