Dátum: 2000. december 7., 19:39
Feladó: Külső hozzászóló --
Tárgy: > eleg lenne nekem egy ilyen kis gagyi keplet is, es
> eleg lenne nekem egy ilyen kis gagyi keplet is, es
> esetleg ha nagyon nagy pontossag kell, akkor bele lehet rakni egy olyan
> reszt is a prg-be, ami lassan, de pontosan szamol, es ezzel lehetne
> korrigalni az idokozben felgyulemlo pontatlansagokat. de most meg boven
> eleg lenne egy nagyon kis pontossaggal szamolo keplet is. (ilyen gondolom
> van, vagy lehet irni, csak az a baj, hogy fizikaval mar eleg regen
> foglalkoztam.)
t: az idopont JD-ben, amire szamitanod kell a poziciokat.
(F,f: fokot jelent)
T=(t-2 415 020)/36525
Merkur:
lambda = 178f 10' 44,68" + 538 106 654,80" * T + 1.084 *T*T
kpuf 53' 58.91" + 5599.76" *T + 1.061 *T *T
OM = 47F 08' 45.40" + 4266.75 *T +0.626 *T*T
e = 0.20561421+0.00002046*T-0.000000030*T*T
i = 7F 0' 10.37" + 6.699*T-0.066*T*T
A=0.3870984
VENUSZ
lambda = 342f 46' 1.39"+210669162.88"*T+1.1148"*T*T
kp0F 09' 49.8" + 5068.93"*T-3.515"*T*T
OM = 75F 46' 46.73" + 3239.46"*T + 1.476"*T*T
e = 0.00682069-0.00004774*T+0.000000091*T*T
i = 3F 23' 37.07" + 3.621"*T-0.0035"*T*T
A=0.72333015
FOLD
lambda = 99F 41' 48.04" + 129602768.13"*T+1.089"*T*T
kp1F 13' 15" + 6189.03"*T+1.63"*T*T+0.012*T*T*T
OM = 0
e = 0.01675104-0.00004180*T-0.000000126*T*T
i = 0
A=1.00000013
MARS
lambda = 293F 44' 51.46" + 68910117.33"*T+1.1184"*T*T
kp34F 13' 5.53" + 6626.73"*T+0.4675"*T*T-0.0043"*T*T*T
OM = 48F 47' 11.19" + 2775.57" *T -0.005"*T*T-0.0192"*T*T*T
e = 0.09331290+0.000092064*T-0.000000077*T*T
i = 1F 51' 1.20" - 2.430*T+0.0454*T*T
A=1.52368839
JUPITER
lambda = 238F 2' 57.32" + 10 930687.148"*T+1.20486"*T*T-0.005936"*T*T*T
kpF 43' 15.34" + 5795.862"*T+3.80258"*T*T-0.01236"*T*T*T
OM = 99F 26' 36.19"+3637.908"*T+1.2680"*T*T-0.03064"*T*T*T
e = 0.04833475+0.000164180*T-0.0000004676*T*T-0.0000000017*T*T*T
i = 1F 18' 31.45" - 20.506"*T+0.014"*T*T
A=5.202561
SZATURNUSZ
lambda = 266F 33' 51.76" + 4404635.5810"*T+1.16835"*T*T-0.021"*T*T*T
kp‘F 5' 53.38"+7050.297"*T+2.9749"*T*T+0.0166"*T*T*T
OM = 112F 47' 25.40"+3143.5025"*T-0.54785"*T*T-0.0191"*T*T*T
e = 0.05589232-0.00034550*T-0.000000728*T*T+0.00000000074*T*T*T
i = 2F 29' 33.07"-14.108"*T-0.05576"*T*T+0.00016*T*T*T
A=9.554747
URANUSZ
lambda = 244F 11' 50.89"+1547508.765"*T+1.13714"*T*T-0.002176"*T*T*T
kp1F 32' 55.14" + 5343.958"*T+0.8539"*T*T-0.00218*T*T*T
OM = 73F 28' 37.55" + 1795.204"*T+4.722"*T*T
e = 0.0463444-0.00002658*T+0.000000077*T*T
i = 0F 46' 20.87" + 2.251"*T+0.1422"*T*T
A.21814
NEPTUNUSZ
lambda = 178f 10' 44,68" + 538 106 654,80" * T + 1.084 *T*T
kpFF 43' 38.37" + 5128.468"*T+1.40694"*T*T-0.002176*T*T*T
OM = 130F 40' 52.89" + 3956.166"*T+0.89952"*T*T-0.002176"*T*T*T
e = 0.00899704+0.000006330*T-0.000000002*T*T
i = 1F 46' 45.27" - 34.357"*T-0.0328"*T*T
A0.10957
(A Plutora e modszer nem jo.)
M = lamba - kp
o = kp - OM
valahogy old meg az E - e*sinE = M egyenletet E-re.
v=2*arctg( ((1+e)/(1-e))^(1/2) tan(E/2) )
r= A*(1 - e *cosE)
k=r*cosv
z=r*sinv
Px=coso*cosOM-sino*sinOM*cosi
Py=coso*sinOM+sino*cosOM*cosi
Pz=sino*sini
Qx=-sino*cosOM-coso*sinOM*cosi
Qy=-sino*sinOM+coso*cosOM*cosi
Qz=coso*sini
x=Px*k+Qx*z
y=Py*k+Qy*z
z=Pz*k+Qz*z
x,y,z a bolygo heliocentrikus derekszogu koordinatai. Ebbol mar tudsz barmilyen
mas koordinatat szmaolni, az rad van bizva. A MODSZER NEHANY SZAZ EVIG JO
1900-TOL SZAMITVA, 1 IVPERC EREJEIG!!!!
Szilard
P.S.: Gepelesi hibakert felelosseget nem vallalok. Irodalom: Csillagaszat,
(1989), szerk. Marik Miklos, Akademiai Kiado.