Simpsonin sääntö

Simpsonin sääntö on numeerinen menetelmä, jolla voidaan approksimoida määrättyä integraalia, kun integroitavan funktion integraalifunktiota ei tunneta tai haluta käyttää. Jos väli [a,b] jaetaan kahteen yhtä pitkään osaväliin, niin tällöin pätee likimääräisesti:

.

Usein merkitään osavälien pituutta kirjaimella h = (b-a)/2, jolloin yllä oleva kaava voidaan kirjoittaa muodossa:
.

Simpsonin sääntö on nimetty englantilaisen matemaatikon Thomas Simpsonin (1710–1761) mukaan.

Johtaminen 
Simpsonin sääntö voidaan johtaa useammalla eri tavalla, joilla kaikilla saadaan luonnollisesti sama tulos. Tässä artikkelissa paneudutaan kuitenkin vain yhteen tapaan. Muista tavoista kiinnostunut lukija voi katsoa lisätietoa viitteistä.

Interpolointi toisen asteen polynomifunktiolla 
Lähtökohtana on korvata funktio f toisen asteen polynomifunktiolla g(x)=px2+qx+r, joka saa samat arvot kuin funktio f välin [a,b] päätepisteissä ja keskipisteessä. Lasketaan nämä arvot, sillä niitä tarvitaan lopuksi:

Tällöin integraali voidaan kirjoittaa polynomin g avulla muodossa:

Hyödyntämällä kaavoja 

saadaan edelleen

 

Edellinen lasku on hieman helpompi laskea loppuun jos tehdään aluksi muuttujan vaihto:

jolloin

ja integraali tulee muotoon

Simpsonin säännön virhe 
Jaettaessa väli [a,b] kahteen yhtä suureen osaan eli n=2 jolloin osavälien pituus on h=(b-a)/2 saadaan Simpsonin säännön virheeksi 

missä on jokin piste väliltä (a,b) ja on funktion f neljäs derivaatta kyseisessä pisteessä.
 
Simpsonin sääntö antaa tarkan vastauksen kaikille polynomeille, joiden asteluku on pienempi tai yhtä suuri kuin kolme, koska virhetermi sisältää funktion neljännen derivaatan, mikä on tietysti = 0 kun funktio f on korkeintaan kolmatta astetta.

Yleistetty Simpsonin sääntö 
Koska Simpsonin säännön virhe on suoraan verrannollinen laskettavan integroimisvälin pituuteen, saadaan virhe mahdollisimman pieneksi kun välin pituutta pienennetään. Käytännössä tämä tehdään jakamalla väli [a,b] useaan pienempään osaväliin [xi,xi+h], joihin jokaiseen käytetään Simpsonin sääntöä. Jos väli [a,b] jaetaan yhtä pitkiin väleihin, joita on N=2n kappaletta, (välejä tulee olla parillinen määrä, sillä Simpsonin sääntöä sovelletaan aina kahteen osaväliin) tulee jokaisen osavälin pituudeksi h=(b-a)/N. Tällöin saadaan yleistetty Simpsonin sääntö:

Yleistetyn Simpsonin säännön virhe  
Jaettaessa väli [a,b] yhtä suuriin osaväleihin, joita on N=2n kappaletta, jolloin osavälien pituus on h=(b-a)/N, saadaan yleistetyn Simpsonin säännön virheeksi:

missä on jokin välin (a,b) piste

Simpsonin 3/8-sääntö 
Jos funktiota f ei approksimoidakaan toisen asteen polynomifunktiolla g (kahdella osavälillä) vaan kolmannen asteen polynomifunktiolla (kolmella osavälillä), jolta vaaditaan, että g saa samat arvot kuin funktio f osavälien päätepisteissä, saadaan likimääräiseksi arvioksi:

Simpsonin 3/8-säännön virhe 

missä on jokin välin (a,b) piste

Lähteet

Viitteet

Kirjallisuutta 
 
 

Numeerinen analyysi