Tähtitiede:interpolointi
interpolointi
interpolointi |
Esimerkiksi tähtitieteellisissä vuosikirjoissa annetut planeettojan koordinaatit on ilmoitettu vain yhden tai useamman päivän välein. Jos muutokset ovat hitaita, voit arvioida koordinaatit mielivaltaiselle hetkelle olettamalla, että noiden kahden hetken välillä koordinaatti muuttuu tasaisella nopeudella. Tätä taulukkoarvojen väliin jäävien lukujen arviointia tarkoittaa interpolointi. Jos oletat muutoksen tasaiseksi, kyseessä on lineaarinen interpolointi.
Oletetaan, että taulukossa on ilmoitettu jonkin suureen arvot hetkille t0 ja t1, ja haluat sen hetkelle t, joka on näiden välissä. Jos suureen arvo hetkellä t0 on f0 ja hetkellä t1 f0, sen arvo hetkellä t on likimain
f = f0 + [ (t - t0) / (t1 - t0) ] (f1 - f0).
Esimerkki: Tähdet 1993 -kirjan taulukon mukaan Marsin rektaskensio 23.4. tunneiksi muutettuna oli 8.038 ja 9.5. 8.575 tuntia. Mikä oli sen rektaskensio 2.5?
Aikavälit vuorokausina ovat t1 - t0 = 16 ja t - t0 = 9. Saamme siis rektaskensioksi
8.038 + (9 / 16) (8.575 - 8.038) = 8.340.
Vastaava menetelmä soveltuu toki myös taulukoihin, joissa muuttujana ei ola aika, vaan jokin muu suure. Mikäli muutos ei ole tasaista tai tarvitaan suurempaa tarkkuutta, voidaan käyttää useamman pisteen interpolaatiota. Mistä tahansa numeerisen analyysin oppikirjasta löytyy suuri määrä erilaisiin tilanteisiin soveltuvia menetelmiä.
Seuraava kolmeen pisteeseen perustuva kvadraattinen interpolointi on lineaarisen interpoloinnin yksinkertaisin täsmennys. Siinä laskettavaa suuretta kuvataan toisen asteen käyrällä. Esimerkiksi planeettojen koordinaateille tällainen menetelmä antaa jo hyvin tarkat arvot.
Oletetaan, että hetkiä t0, t1 = t1 + h ja t2 = t1 + 2h vastaavat muuttujan arvot ovat f0, f1 ja f2. Nämä taulukkoarvot on valittu siten, että t1 on lähinnä haluttua hetkeä t.
Lasketaan aluksi interpolointivälin suhde taulukointiväliin:
d = (t - t1) / h.
Muuttujan arvo hetkellä t on silloin
f = f1 + (f2 - f0) d / 2 + (f2 - 2f1 + f0) d2 / 2
Esimerkki: Lasketaan edellisen esimerkin mukainen Marsin paikka tarkemmin, kun tiedämme, että 25.5 sen rektaskensio oli 9.135 tuntia.
Ensinnäkin saamme
d = (2 - 9) / 16 = -0.4375.
Rektaskensio on siten
8.575 - (9.135 - 8.038)(0.4375 / 2) + (9.135 - 2×8.575 + 8.038) 0.43752/2 = 8.337.
Kolmen pisteen interpolointi antaa esimerkiksi planeettojen koordinaatit riittävän tarkasti useimpiin tarkoituksiin, jos taulukoidut arvot on annettu vuorokauden välein. Käyttämällä vielä useampia pisteitä, saadaan yhä tarkempia interpolointikaavoja, mutta samalla ne tulevat mutkikkaammiksi ja niiden laskeminen käy hitaammaksi.
Erikieliset vastineet
interpolation | englanti (English) |
Käytetyt lähteet
Alaviitteet
Lähdeviittaus tähän sivuun:
Tieteen termipankki 23.12.2024: Tähtitiede:interpolointi. (Tarkka osoite: https://tieteentermipankki.fi/wiki/Tähtitiede:interpolointi.)