Trapezoidna metoda u integralima
U primijenjenoj matematici često se susrećemo s problemom izračunavanja površine ispod krive. Teoretski, ovo se može riješiti korištenjem određenih integrala. Međutim, u praksi se ne mogu sve funkcije lako analitički integrirati. Neke funkcije su složene, podaci su dostupni samo u tabelarnom obliku ili modelima nedostaju jednostavni antiderivativni elementi. U takvim situacijama, numeričke metode postaju neophodne. Jedna od najpoznatijih i najčešće korištenih numeričkih metoda je trapezoidno pravilo, tehnika aproksimacije za integrale koja površinu ispod krive dijeli na trapezoidne oblike.
Osnovni koncept određenog integrala
Geometrijski, određeni integral \(\int_a^bf(x)\,dx\) može se shvatiti kao površina ograničena krivuljom \(y=f(x)\), osom \(x\) i vertikalnim linijama \(x=a\) i \(x=b\). Ako je \(f(x)\ge 0\), integral je pozitivna površina. Ako funkcija ima negativne vrijednosti na određenom intervalu, integral daje „predznačenu površinu“.
Glavni problem nastaje kada:
1. Funkcija \(f(x)\) nema antiderivaciju koja se može izraziti elementarnim funkcijama.
2. Vrijednost \(f(x)\) je poznata samo u nekim tačkama (eksperimentalni podaci).
3. Simboličko računanje je neefikasno ili nemoguće.
Ovdje se koristi trapezoidna metoda za aproksimaciju integralne vrijednosti korištenjem vrijednosti funkcije u određenim tačkama.
Ideja trapezoidne metode
Trapezoidna metoda počinje jednostavnom idejom: aproksimirati krivulju pravom linijom na svakom podintervalu. Dok Riemann-ova metoda sume koristi okomite stranice i pravougaonu površinu, trapezoidna metoda koristi dvije krajnje tačke na podintervalu i povezuje ih pravom linijom. Površina ispod prave linije formira trapezoid, a ne pravougaonik.
Pretpostavimo da želimo izračunati:
\[
\int_a^bf(x)\,dx
\]
Interval \([a,b]\) dijelimo na \(n\) podintervala jednake dužine. Dužina svakog podintervala je:
\[
h=\frac{ba}{n}
\]
Tačke podjele:
\[
x_0=a,\;x_1=a+h,\;x_2=a+2h,\;tačke,\;x_n=b
\]
Sa vrijednostima funkcije:
\[
f(x_0), f(x_1), \tačke, f(x_n)
\]
Na svakom podintervalu \([x_{i-1}, x_i]\), površina ispod krive je aproksimirana površinom trapeza:
\[
A_i \approx \frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Zatim se ukupna površina (integralna aproksimacija) dobija sabiranjem svih trapezoida:
\[
\int_a^bf(x)\,dx \approx \sum_{i=1}^{n}\frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Ako se sredi, formula složenog trapezoidnog pravila postaje:
\[
\int_a^b f(x)\,dx \approx \frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+\cdots+2f(x_{n-1})+f(x_n)\right]
\]
Obratite pažnju na obrazac koeficijenata: krajnje tačke \(x_0\) i \(x_n\) se množe sa 1, dok se tačke u sredini množe sa 2.
Jednostavan primjer izračuna
Pretpostavimo da želimo pristupiti:
\[
\int_0^2 x^2\,dx
\]
Analitički, integral je \(\left.\frac{x^3}{3}\right|_0^2=\frac{8}{3}\approx 2{,}6667\). Međutim, koristit ćemo trapezoidnu metodu sa \(n=4\).
1. Interval \([0,2]\), tada:
\[
h=\frac{2-0}{4}=0{,}5
\]
2. Titik-titik: \(x_0=0\), \(x_1=0{,}5\), \(x_2=1\), \(x_3=1{,}5\), \(x_4=2\)
3. Vrijednost funkcije:
\(f(0)=0\)
(f(0{,}5)=0{,}25)
\(f(1)=1\)
(f(1{,}5)=2{,}25)
\(f(2)=4\)
Koristite formulu složenog trapeza:
\[
T=\frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+2f(x_3)+f(x_4)\right]
\]
\[
T=\frac{0{,}5}{2}\left[0+2(0{,}25)+2(1)+2(2{,}25)+4\right]
\]
\[
T=0{,}25\left[0{,}5+2+4{,}5+4\right]=0{,}25(11)=2{,}75
\]
Približan rezultat od 2,75 je prilično blizu tačnoj vrijednosti od 2,6667, sa greškom od oko 0,0833.
Greška u trapezoidnoj metodi
Trapezoidna metoda je aproksimacija, što znači da uvijek postoji razlika između približne vrijednosti i stvarne integralne vrijednosti. Na veličinu greške utiču:
1. Broj podintervala \(n\): što je veći \(n\), to je manji \(h\) i generalno su rezultati tačniji.
2. Zakrivljenost funkcije: ako je funkcija vrlo zakrivljena (ima veliki drugi izvod), trapezoidna metoda može biti manje tačna od metoda višeg reda.
Za dovoljno glatke funkcije (koje imaju kontinuirane druge derivacije), greška složene trapezoidne metode može se procijeniti:
\[
E_T = -\frac{(ba)}{12}h^2 f”(\xi)
\]
za neko \(\xi\) između \(a\) i \(b\). Iz ove formule se može vidjeti da je greška proporcionalna \(h^2\). To jest, ako prepolovimo \(h\) (tj. udvostručimo broj podintervala), greška se smanjuje za približno jednu četvrtinu.
Prednosti trapezoidne metode
Trapezoidna metoda je popularna iz nekoliko razloga:
1. Jednostavno i brzo
Izračun je relativno jednostavan i može se čak izvršiti ručno za manje slučajeve.
2. Pogodno za tabelarne podatke
Ako je vrijednost \(f(x)\) poznata samo u diskretnim tačkama (npr. rezultati mjerenja), ova metoda se može direktno primijeniti.
3. Bolja tačnost od pravougaone metode
Budući da koristi dvije krajnje tačke i aproksimira funkciju pravom linijom, generalno je tačnija od pravougaone aproksimacije za isti broj particija.
4. Stabilan za mnoge inženjerske primjene
Mnogi fizički i inženjerski problemi (npr. izračunavanje rada, energije, površine poprečnog presjeka, pražnjenja itd.) mogu se riješiti ovom metodom.
Ograničenja trapezoidne metode
Iako korisna, ova metoda ima ograničenja:
1. Manje precizno za visoko zakrivljene funkcije
Budući da koristi samo linearnu aproksimaciju na svakom podintervalu, funkcije s brzim promjenama ili oštrim krivuljama zahtijevaju veliko \(n\).
2. Nije metoda višeg reda
Greške reda veličine \(h^2\); druge metode poput Simpsonove mogu dati manje greške na istom broju podintervala (pod uslovom da je funkcija dovoljno glatka).
3. Osjetljivo na odabir intervala
Ako je interval preširok i \(n\) mali, rezultati mogu značajno odstupati.
Zatvaranje
Trapezoidna metoda integracije jedna je od najosnovnijih tehnika u numeričkoj integraciji. Dijeljenjem intervala na male dijelove i zamjenom krivulje pravom linijom u svakom dijelu, možemo aproksimirati integral bez potrebe za analitičkim pronalaženjem antiderivacije. Formula je jednostavna, lako se implementira i posebno korisna kada su podaci diskretni ili je funkciju teško integrirati.
Međutim, korisnici trapezoidne metode moraju obratiti pažnju na broj podintervala i ponašanje funkcije koja se integrira. Za visoko zakrivljene funkcije ili kada je potrebna visoka preciznost, trapezoidna metoda se i dalje može koristiti povećanjem \(n\), ili se mogu razmotriti druge numeričke metode poput Simpsonove ili Rombergove metode. Uz dobro razumijevanje, trapezoidna metoda postaje efikasan i praktičan alat za rješavanje različitih integralnih problema u nauci, inženjerstvu i ekonomiji.
