Trapetsimeetod integraalides
Rakendusmatemaatikas puutume sageli kokku kõvera aluse pindala arvutamise probleemiga. Teoreetiliselt saab seda lahendada kindlate integraalide abil. Praktikas ei ole aga kõiki funktsioone analüütiliselt lihtne integreerida. Mõned funktsioonid on keerulised, andmed on saadaval ainult tabelina või mudelitel puuduvad lihtsad algfunktsioonid. Sellistes olukordades muutuvad numbrilised meetodid hädavajalikuks. Üks tuntumaid ja sagedamini kasutatavaid numbrilisi meetodeid on trapetsireegel, integraalide lähendustehnika, mis jagab kõvera aluse pindala trapetskujulisteks kujunditeks.
Määratud integraali põhimõiste
Geomeetriliselt võib kindlat integraali (\int_a^bf(x)\,dx\) mõista kui pindala, mida piiravad kõver (y=f(x)\), x-telg ja vertikaaljooned (x=a\) ja (x=b\). Kui (f(x)\ge 0\), on integraal positiivne pindala. Kui funktsioonil on teatud intervallis negatiivsed väärtused, annab integraal „märgilise pindala“.
Peamine probleem tekib siis, kui:
1. Funktsioonil \(f(x)\) puudub algfunktsioon, mida saaks väljendada elementaarfunktsioonidega.
2. Funktsiooni f(x) väärtus on teada ainult teatud punktides (eksperimentaalsed andmed).
3. Sümboolne arvutamine on ebaefektiivne või võimatu.
Siin kasutatakse trapetsikujulist meetodit integraali väärtuse ligikaudseks määramiseks, kasutades funktsiooni väärtusi teatud punktides.
Trapetsikujulise meetodi idee
Trapetsimeetod algab lihtsast ideest: lähendada kõver sirge joonega igas osaintervallis. Kui Riemanni summa meetod kasutab risti asetsevaid külgi ja ristkülikukujulist pindala, siis trapetsimeetod kasutab osaintervallil kahte lõpp-punkti ja ühendab need sirge joonega. Sirge all olev pindala moodustab trapetsi, mitte ristküliku.
Oletame, et tahame arvutada:
\[
_int_a^bf(x)_dx
\]
Jagame intervalli \([a,b]\) \(n\) võrdse pikkusega alamintervalliks. Iga alamintervalli pikkus on:
\[
h=\frac{ba}{n}
\]
Jagamispunktid:
\[
x_0=a,\;x_1=a+h,\;x_2=a+2h,\;\punktid,\;x_n=b
\]
Funktsioonide väärtustega:
\[
f(x_0), f(x_1), \dots, f(x_n)
\]
Igal alamintervallil \([x_{i-1}, x_i]\) on kõvera alune pindala ligikaudne trapetsi pindalaga:
\[
A_i \umx \frac{h}{2}\vasak[f(x_{i-1}) + f(x_i)\parem]
\]
Seejärel saadakse kogupindala (integraalne lähendus) kõigi trapetside liitmise teel:
\[
int_a^bf(x)\,dx \approx \sum_{i=1}^{n}\frac{h}{2}\vasak[f(x_{i-1}) + f(x_i)\parem]
\]
Korrastamisel saab trapetsikujulise joone liitvalemiks järgmine:
\[
\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]
\]
Pane tähele koefitsientide mustrit: lõpp-punktid \(x_0\) ja \(x_n\) korrutatakse 1-ga, samas kui keskmised punktid korrutatakse 2-ga.
Lihtne arvutusnäide
Oletame, et tahame läheneda:
\[
int_0^2 x^2 dx
\]
Analüütiliselt on integraal \(\left.\frac{x^3}{3}\right|_0^2=\frac{8}{3}\approx 2{,}6667\). Meie aga kasutame trapetsimeetodit, kus \(n=4\).
1. Intervall \([0,2]\), siis:
\[
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. Funktsiooni väärtus:
(f(0)=0)
\(f(0{,}5)=0{,}25\)
(f(1)=1)
\(f(1{,}5)=2{,}25\)
(f(2)=4)
Kasutage trapetsi liitvalemit:
\[
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
\]
Ligikaudne tulemus 2,75 on üsna lähedal täpsele väärtusele 2,6667, veaga umbes 0,0833.
Trapetsikujulise meetodi viga
Trapetsimeetod on lähendusmeetod, mis tähendab, et ligikaudse väärtuse ja tegeliku integraalväärtuse vahel on alati erinevus. Vea suurust mõjutavad:
1. Alamintervallide arv \(n\): mida suurem \(n\), seda väiksem \(h\) ja üldiselt täpsemad on tulemused.
2. Funktsiooni kõverus: kui funktsioon on väga kõver (suure teise tuletisega), võib trapetsimeetod olla vähem täpne kui kõrgema järgu meetodid.
Piisavalt siledate funktsioonide (millel on pidevad teised tuletised) korral saab liittrapetsi meetodi viga hinnata:
\[
E_T = -\frac{(ba)}{12}h^2 f”(\xi)
\]
mingi \(\xi\) väärtuse korral vahemikus \(a\) kuni \(b\). Sellest valemist on näha, et viga on võrdeline \(h^2\)-ga. See tähendab, et kui me vähendame \(h\) poole võrra (st kahekordistame alamintervallide arvu), väheneb viga ligikaudu veerandi võrra.
Trapetsikujulise meetodi eelised
Trapetsikujuline meetod on populaarne mitmel põhjusel:
1. Lihtne ja kiire
Arvutamine on suhteliselt lihtne ja väikeste juhtumite korral saab seda teha isegi käsitsi.
2. Sobib tabelina esitatud andmete jaoks
Kui f(x) väärtus on teada ainult diskreetsetes punktides (nt mõõtmistulemustes), saab seda meetodit otse rakendada.
3. Parem täpsus kui ristkülikukujulisel meetodil
Kuna see kasutab kahte lõpp-punkti ja lähendab funktsiooni sirgega, on see sama arvu partitsioonide korral üldiselt täpsem kui ristkülikukujuline lähendus.
4. Stabiilne paljudes insenerirakendustes
Selle meetodi abil saab lahendada paljusid füüsika- ja inseneriprobleeme (nt töö, energia, ristlõikepindala, vooluhulga jne arvutamine).
Trapetsikujulise meetodi piirangud
Kuigi see meetod on kasulik, on sellel piirangud:
1. Vähem täpne väga kõverate funktsioonide puhul
Kuna see kasutab iga alamintervalli puhul ainult lineaarset lähendust, vajavad kiirete muutuste või teravate kõveratega funktsioonid suurt \(n\).
2. Mitte kõrgetasemeline meetod
Vead suurusjärgus \(h^2\); teised meetodid, näiteks Simpsoni meetod, võivad sama arvu alamintervallide korral anda väiksemaid vigu (eeldusel, et funktsioon on piisavalt sujuv).
3. Tundlik intervalli valiku suhtes
Kui intervall on liiga lai ja n on väike, võivad tulemused oluliselt erineda.
Sulgemine
Trapetsmeetod integreerimiseks on numbrilise integreerimise üks põhilisemaid tehnikaid. Jagades intervalli väikesteks osadeks ja asendades kõvera igas osas sirgjoonega, saame integraali lähendada ilma algfunktsiooni analüütiliselt leidmata. Valem on lihtne, hõlpsasti rakendatav ja eriti kasulik, kui andmed on diskreetsed või funktsiooni on raske integreerida.
Trapetsmeetodi kasutajad peavad aga pöörama tähelepanu alamintervallide arvule ja integreeritava funktsiooni käitumisele. Väga kõverate funktsioonide või suure täpsuse vajaduse korral saab trapetsmeetodit endiselt kasutada, suurendades \(n\) väärtust, või kaaluda võib muid numbrilisi meetodeid, näiteks Simpsoni või Rombergi meetodeid. Hea arusaamise korral saab trapetsmeetodist tõhus ja praktiline tööriist mitmesuguste integraalsete probleemide lahendamiseks teaduses, inseneriteaduses ja majanduses.
