Contact info
Word
Research
   Publications
Studies
Free Software
Hobbies
Articles
Photography
About me
   Curriculum Vitae

©Marko Grönroos, 1998

USENET News sfnet.harrastus.astronomia

Säie: Planeettojen rataparametrit

Timo Lindfors <timo PISTE lindfors AT kolumbus PISTE fi> writes: > Nyt tarvitsisin kuitenkin aurinkokunnan tarkan tilan jollakin > ajanhetkellä. Käynnistin vsop87-ohjelman (fortrania, löytyy wepistä jostain, kääntyy f2c:llä ja gcc:llä) ja copy-pastesin siitä noita koordinaatteja vuoden 2001 uudenvuodenhetkelle. Vasemmalla XYZ ovat koordinaatit ja luulisin että X'Y'Z' taitaa olla liikevektori AU:ta per päivä. VSOP87 on ns. hillittömän tarkka (tarkkuus jotain 0.001-0.01 kaarisekunttia, riippuen planeetasta ja aikavälistä). Siinä ei kuitenkaan ole Pluton rataa. Huomaa, että aika on TDB:tä.       PLANETARY SOLUTION VSOP87A BODY : MERCURY       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 0.1687266681 au X' : 0.0203799656 ua/d       Y : -0.4096394268 au Y' : 0.0121338833 ua/d       Z : -0.0489499153 au Z' : -0.0008793178 ua/d       PLANETARY SOLUTION VSOP87A BODY : VENUS       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 0.4943091928 au X' : -0.0148141797 ua/d       Y : 0.5270081943 au Y' : 0.0137504798 ua/d       Z : -0.0213286722 au Z' : 0.0010430238 ua/d       PLANETARY SOLUTION VSOP87A BODY : EARTH       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : -0.1813707407 au X' : -0.0171908893 ua/d       Y : 0.9664440449 au Y' : -0.0032434938 ua/d       Z : 0.0000004907 au Z' : 0.0000002783 ua/d       PLANETARY SOLUTION VSOP87A BODY : MARS       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : -1.6473415133 au X' : 0.0010132438 ua/d       Y : -0.0575231270 au Y' : -0.0127910557 ua/d       Z : 0.0392847846 au Z' : -0.0002928412 ua/d       PLANETARY SOLUTION VSOP87A BODY : JUPITER       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 1.7967535956 au X' : -0.0071540176 ua/d       Y : 4.7149317120 au Y' : 0.0030468382 ua/d       Z : -0.0597981925 au Z' : 0.0001475004 ua/d       PLANETARY SOLUTION VSOP87A BODY : SATURN       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 4.6857555359 au X' : -0.0050854089 ua/d       Y : 7.8099720508 au Y' : 0.0028614295 ua/d       Z : -0.3222348696 au Z' : 0.0001524501 ua/d       PLANETARY SOLUTION VSOP87A BODY : URANUS       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 15.3733746418 au X' : 0.0024718881 ua/d       Y :-12.7245485198 au Y' : 0.0028499740 ua/d       Z : -0.2465780580 au Z' : -0.0000213776 ua/d       PLANETARY SOLUTION VSOP87A BODY : NEPTUNE       -----------------------------------------------------------------------------       Heliocentric rectangular coordinates. Dynamical equinox and ecliptic J2000.       Date : 01 / 01 / 2001 Julian date : 2451910.500000       Time : 00h 00m 00s TDB Precision of the complete solution       X : 17.7418959191 au X' : 0.0025080201 ua/d       Y :-24.3248267429 au Y' : 0.0018712935 ua/d       Z : 0.0920971334 au Z' : -0.0000962588 ua/d Oma rinnakkaistettu N:n kappaleen simulointiohjelmani (C++) löytyy osoitteesta: http://magi.yok.utu.fi/~magi/opinnot/mpi/ Siinä on esimerkkinä myös planeettojen ratakierto, vaikkakin melko likimääräisillä arvoilla. -- -- Marko Grönroos, magi AT iki PISTE fi (http://www.iki.fi/magi/) -- Paradoxes are the source of truth and the end of wisdom

Timo Lindfors <timo PISTE lindfors AT kolumbus PISTE fi> writes: > Mieleeni tuli seuraavia asioita; > a) Luultavasti nuo koordinaatit esittävät siis sitä paikka jossa > planeetta *näkyy*. Käsittääkseni valon äärellisen nopeuden takia menetän > noita koordinaatteja käyttämällä paljon tarkkuudessa. Nuo ovat käsittääkseni ihan "oikeat" heliosentriset koordinaatit tiettynä ajanhetkenä. Jos halutaan laskea näennäiset geosentriset koordinaatit, täytyy valonnopeuden aiheuttamat virheet ottaa huomioon erikseen iteroimalla kerran tai kaksi. En nyt muista tarkkaan mitä kaikkea noissa laskuissa piti ottaa huomioon. Vähän kaikenlaista. Suosittelen Jean Meeuksen kirjaa Astronomical Algorithms. Siinä on kaikki nämä jutut suht "yksinkertaisesti" ilmaistuna.Tosin niiden juttujen koodaukseen voi mennä toinenkin viikko. Jos haluaa päästä helpolla, suosittelen valmiita ohjelmia, joita löytyy ainakin kaupallisena. Itselläni on joitain noista laskuista toteutettuna, mutta ei aivan kaikkia. > b) Oma mallini olettaa, että painovoiman vaikutus etenee äärettömän > suurella nopeudella. Tuleekohan tästä virheitä? En tiedä miten tuo vaikuttaa. Newtonin kaavoillahan ei pärjää kuin tiettyyn tarkkuuteen asti, sehän se suhteellisuusteorian suuri juttu juuri oli. Merkuriuksen radassa virhe on jopa ihan mitattava. En kyllä yhtään tiedä miten nuo laskut täytyisi silloin tehdä. Voi myös olla että iterointi muutenkin aiheuttaa sen verran virhettä ettei riittävään tarkkuuteen pääse muutenkaan. En tiedä, en ole kokeillut. > c) Miksiköhän Maa ei ollut mukana? Olihan se siinä kolmantena. Tietty maan liikkeessä täytyy ottaa huomioon kuu, eli sinun pitäisi tietää senkin paikka. Se on kuitenkin suht helppo laskea riittävällä tarkkuudella. En tosin tiedä missä olisi valmis implementaatio, itselläni ei muistaakseni vielä ole. Saattaisit päästä hieman helpommalla jos et ota kuuta malliisi. Eli tekisit laskut maa-kuu-barysentrin suhteen, jolloin tosin geosentristen koordinaattien lasku ei jälkeenpäin taida olla enää mahdollista. Tarvitset lisäksi myös Pluton, koska on tunnettua että se häiritsee Neptunuksen rataa, varsinkin nyt. Ei tosiaankaan mitenkään helppo projekti, jos ihan oikeasti haluat iteroida mahdollisimman tarkoilla koordinaateilla... -- -- Marko Grönroos, magi AT iki PISTE fi (http://www.iki.fi/magi/) -- Paradoxes are the source of truth and the end of wisdom

Edellinen säie: Mustista aukoista
Seuraava säie: Valoa nopeammin ??
[Muut säikeet] [Muut uutisryhmät]