- MATH.APP.120
- 6. Integraalikäsitteen laajennuksia ja sovelluksia
- 6.4 Symbolinen ja numeerinen integrointi
\({}^*\)Symbolinen ja numeerinen integrointi¶
Symbolinen integrointi¶
Käytännössä integraalifunktioiden ja määrättyjen integraalien laskeminen perustuen integroimissääntöihin voi olla työlästä ja virhealtista. Onneksi on olemassa runsaasti symbolisen laskennan ohjelmistoja ja sovelluksia, joilla rutiini-integroinnit voidaan suorittaa. Seuraavassa esimerkissä tarkastellaan kahta opintojen kannalta keskeistä sovellusta.
Esimerkki 6.4.1
Lasketaan esimerkin 5.5.17 kohdan 4 määrätty integraali kahdella eri sovelluksella.
Matlab ja sen Symbolic Math Toolbox:
syms x;
int(x/(3*x+2)^(1/3),x,-1/3,2)
antaa tulokseksi ans = 16/15
.
int(x/(3*x+2)^(1/3),x,-1/3,2)
antaa tulokseksi
Numeerinen integrointi¶
Sovelluksissa törmätään usein tilanteisiin, joissa
- integrointi alkeisfunktioiden avulla ei onnistu (esimerkiksi \(f(x)=e^{x^2}\)) tai on vaikeaa,
- funktion \(f\) lauseketta ei tunneta, vaan tiedetään vain sen arvoja tietyissä pisteissä esimerkiksi mittaustuloksina.
Tällöin funktion \(f\) integraalia voidaa arvioida numeerisella integroinnilla käyttäen funktion \(f\) arvoja äärellisen monessa integroimisvälin pisteessä.
Riemannin summa¶
Jos \(P=\{x_0,x_1,\ldots,x_n\}\) on välin \([a,b]\) jako, niin mikä tahansa Riemannin summa antaa funktion \(f\) integraalille välillä \([a,b]\) arvion
Jos valitaan tasavälinen jako, jossa kunkin osavälin pituus on \(h\), sievenee arvio muotoon
Jos \(f\) on ei-negatiivinen, niin geometrinen tulkinta arviolle on se, että jokaisella välillä \([x_{i-1},x_i]\) funktion \(f\) kuvaajan ja \(x\)-akselin väliin jäävän alueen pinta-alaa arvioidaan suorakulmion pinta-alalla (vertaa kuvaan Riemannin summasta).
Esimerkki 6.4.2
Arvioi integraalia
Riemannin summalla, kun käytetään tasavälistä jakoa, jolle \(n=6\) ja \(x_i^*\) on osavälin keskipiste.
Nyt \(h=\frac{b - a}{n}=\frac26=\frac13\) ja välien keskipisteet ovat \(\frac76, \frac96,\ldots, \frac{17}{6}\), joten
Vertaa tarkkaan arvoon \(\ln(3)=1{,}098~612~288\cdots\).
Käytännössä Riemannin summaa ei juurikaan käytetä integraalin arvioimiseen, sillä voidaan kehittää huomattavasti tehokkaampia menetelmiä, joissa samalla määrällä jakopisteitä (eli samalla vaivalla tai tietokoneajalla) päästään huomattavasti parempaan tarkkuuteen. Käsitellään seuraavassa kahta yksinkertaista menetelmää.
Puolisuunnikassääntö¶
Puolisuunnikassäännön (trapezoid rule) ideana on (kun \(f\) on ei-negatiivinen) käyttää funktion \(f\) kuvaajan ja \(x\)-akselin väliin jäävän alueen pinta-alan arvioinnissa suorakulmioiden sijasta puolisuunnikkaita. Ne saadaan aikaan korvaamalla funktion \(f\) kuvaaja pisteiden \((x_i,f(x_i))\) kautta kulkevalla murtoviivalla. Käytetään tasavälistä jakoa, jossa osavälin pituus on \(h\). Tällöin puolisuunnikkaan \(i\) pinta-ala on
ja pinta-alojen summa on
Siis funktion \(f\) integraalille saadaan arvio
Arvio (1) on voimassa myös yleiselle \(f\) (eli vaikka \(f\) ei olisi ei-negatiivinen). Jos funktion \(g(x)\) kuvaaja on pisteiden \((x_i, f(x_i))\) kautta kulkeva murtoviiva, niin välillä \([x_{i-1},x_i]\) on
Integroimalla saadaan
ja summaamalla yli kaikkien osavälien
Esimerkki 6.4.3
Arvioi puolisuunnikassäännöllä samaa integraalia kuin Riemannin summilla arvioitiin esimerkissä 6.4.2, kun käytetään tasavälistä jakoa, jolle \(n=6\).
Nyt \(h=\frac13\) ja jakopisteet ovat \(1, \frac43, \frac53,\ldots, 3\), joten
Simpsonin sääntö¶
Yleensä vielä parempaan arvioon päädytään, jos suorien sijaan korvataan funktion \(f\) kuvaaja paraabelin kaarilla. Simpsonin kaavassa käytetään kolmen peräkkäisen pisteen \((x_i, f(x_i))\), \((x_{i + 1}, f(x_{i + 1}))\) ja \((x_{i + 2}, f(x_{i + 2}))\) kautta kulkevaa paraabelia. Jako valitaan tasaväliseksi, jossa osavälin pituus on \(h\) ja jossa on parillinen määrä osavälejä.
Oletetaan ensin, että \(x_0=-h\), \(x_1=0\) ja \(x_2=h\). Olkoon \(y(x)=Ax^2+Bx+C\) se toisen asteen polynomi, jonka kuvaaja kulkee pisteiden \((x_0,f(x_0))\), \((x_1,f(x_1))\) ja \((x_2,f(x_2))\) kautta. Nyt
Kauttakulkuehdot ovat
joten \(f(x_0)+4f(x_1)+f(x_2)=2Ah^2+6C\). Saatiin siis
Tämä kaava on voimassa myös ilman oletusta \(x_1=0\). Erityisesti
ja vastaavalla tavoin kaikilla väleillä \([x_{2i},x_{2(i+1)}]\), joten päädytään arvioon
Esimerkki 6.4.4
Arvioi esimerkin 6.4.2 integraalia Simpsonin säännöllä kuin Riemannin summilla ja puolisuunnikassäännöllä, kun käytetään tasavälistä jakoa, jolle \(n=6\).
Nyt \(h=\frac13\) ja jakopisteet ovat \(1, \frac43, \frac53,\ldots, 3\), joten
Huomautus 6.4.5
Riemannin summaa, jossa on \(n\) osaväliä ja \(x_i^*\) on osavälin keskipiste, kutsutaan keskipisteapproksimaatioksi \(M_n\). Vastaavasti \(n\) osavälin puolisuunnikassäännön antamaa arviota merkitään \(T_n\). Tällöin Simpsonin kaavan antama arvio yhteensä \(2n\) osavälillä on
Simpson-arvio saadaan siis keskipiste- ja puolisuunnikasapproksimaatioiden sopivasti painotettuna keskiarvona.
Gaussin numeerinen integrointi¶
Erityisesti fysikaalisten insinööritieteiden numeriikassa joudutaan ratkaisemaan erilaisia kenttäongelmia, joiden systemaattinen käsittely perustuun yleensä niin sanotun elementtimenetelmän (FEM) käyttöön. Menetelmässä on tarve laskea suuri määrä integraaleja mahdollisimman nopeasti kuitenkin niin, ettei laskenta-aika pitenisi kohtuuttomasti. Yleensä tässä menetelmässä integrointi perustuu Gaussin numeeriseen integrointiin tai kvadratuuriin.
Yleisesti Gaussin integroinnissa arvioidaan integraalia kaavalla
missä
- \(x_i\in[-1,1]\) ovat integrointipisteet,
- \(w_i\) pisteisiin liittyvät painokertoimet,
- \(n\) on kvadratuurin kertaluku.
Huomautus 6.4.6
Jos halutaan arvioida integroituvan funktion \(f:[a,b]\to\mathbb{R}\) integraalia yli välin \([a,b]\) Gaussin kvadratuurilla, voidaan käyttää muuttujanvaihtokaavaa (sopivan sijoituksen löytäminen jätetään lukijalle)
Integrointipisteet ja painokertoimet ovat annettuja parametreja ja niiden numeerisia arvoja löytyy kirjallisuudesta ja internetistä. Integrointipisteet ja painokertoimet on määrätty siten, että parittoman asteen (\(2n - 1\)) polynomit integroituvat tarkasti. Tarkastellaan vain kahta alimman kertaluvun kvadratuuria. Ensimmäisen kertaluvun Gaussin interointikaava (\(x_1=0\) ja \(w_1=2\)) on
ja toisen kertaluvun (\(x_1=-\frac{1}{ \sqrt{3}}\), \(x_2=\frac{1}{ \sqrt{3}}\), \(w_1=w_2=1\))
Esimerkki 6.4.7
Arvioidaan integraalia
- ja 2. kertaluvun Gaussin kvadratuureilla. Ensimmäisen kertaluvun kvadratuuri antaa
Toisen kertaluvun kvadratuuri puolestaan
Huomataan, että toisen kertaluvun kvadratuuri antaa jo suhteellisen tarkan tuloksen.