So ich hab jetzt zwecks mangelndem mathematischem Geschick einfach zwei Ortskoordinaten in kurzem zeitlichen Abstand genommen und daraus die Geschwindigkeitsvektoren errechnet. Hat auch funktioniert. Hier mal ein Ergebniss meiner Arbeit:
Der Jupiter Swing-By von Cassini auf dem Weg zum Saturn:
Beiträge von sonium im Thema „Karthesische Kooridinaten gesucht“
-
-
wie ich die Ortsvektoren berechne weiß ich schon mal, mein Problem ist jetzt nur, wie ich aus diesen Angaben die Vektoren der Geschwindigkeit bekomme. Ne passende Formel finde ich nirgends.
-
Code
Alles anzeigen******************************************************************************* Ephemeris / WWW_USER Thu Jan 20 06:21:39 2005 Pasadena, USA / Horizons ******************************************************************************* Target body name: Mars (499) {source: DE-0406LE-0406} Center body name: Earth (399) {source: DE-0406LE-0406} Center-site name: (User Defined Site) ******************************************************************************* Start time : A.D. 2005-Jan-20 00:00:00.0000 UT Stop time : A.D. 2005-Feb-04 00:00:00.0000 UT Step-size : 1440 minutes ******************************************************************************* Center geodetic : 241.759000, 34.0542, -0.00{E-lon(deg),Lat(deg),Alt(km)} Center cylindric: 241.759000, 5289.8941, 3551.42{E-lon(deg),Dxy(km),Dz(km)} Center pole/equ : High-precision EOP model {East-longitude +} Center radii : 6378.1 x 6378.1 x 6356.8 km {Equator, meridian, pole} Target pole/equ : IAU_MARS {East-longitude -} Target radii : 3396.2 x 3396.2 x 3376.2 km {Equator, meridian, pole} Target primary : Sun {source: DE-0406LE-0406} Interfering body: MOON (Req= 1737.400) km {source: DE-0406LE-0406} Deflecting body : Sun, EARTH {source: DE-0406LE-0406} Deflecting GMs : 1.3271E+11, 3.9860E+05 km^3/s^2 Atmos refraction: NO (AIRLESS) RA format : DEG Time format : CAL EOP file : eop.050119.p050412 EOP coverage : DATA-BASED 1962-JAN-20 TO 2005-JAN-19. PREDICTS-> 2005-APR-11 Units conversion: 1 AU= 149597870.691 km, c= 299792.458 km/s, 1 day= 86400.0 s Table cut-offs 1: Elevation (-90.0deg=NO ),Airmass (>38.000=NO), Daylight (NO ) Table cut-offs 2: Solar Elongation ( 0.0,180.0=NO ) ************************************************************************************* Date__(UT)__HR:MN R.A._(J2000.0)_DEC. dRA*cosD d(DEC)/dt delta deldot ************************************************************************************* 2005-Jan-20 00:00 *m 256.36474 -22.90539 104.24 -11.22 2.1208696875 -11.81075 2005-Jan-21 00:00 *m 257.11820 -22.98109 104.37 -10.71 2.1138452924 -11.84784 2005-Jan-22 00:00 *m 257.87299 -23.05334 104.50 -10.19 2.1068003514 -11.88354 2005-Jan-23 00:00 *m 258.62907 -23.12214 104.63 -9.66 2.0997354212 -11.91792 2005-Jan-24 00:00 *m 259.38642 -23.18746 104.75 -9.14 2.0926510163 -11.95109 2005-Jan-25 00:00 * 260.14502 -23.24928 104.88 -8.61 2.0855476105 -11.98322 2005-Jan-26 00:00 * 260.90482 -23.30758 105.00 -8.08 2.0784256385 -12.01449 2005-Jan-27 00:00 * 261.66580 -23.36235 105.12 -7.55 2.0712854983 -12.04511 2005-Jan-28 00:00 * 262.42792 -23.41356 105.24 -7.01 2.0641275544 -12.07531 2005-Jan-29 00:00 * 263.19115 -23.46120 105.36 -6.47 2.0569521423 -12.10529 2005-Jan-30 00:00 * 263.95545 -23.50525 105.47 -5.93 2.0497595728 -12.13518 2005-Jan-31 00:00 * 264.72078 -23.54569 105.59 -5.39 2.0425501388 -12.16510 2005-Feb-01 00:00 * 265.48710 -23.58251 105.69 -4.84 2.0353241219 -12.19502 2005-Feb-02 00:00 * 266.25437 -23.61569 105.79 -4.29 2.0280818011 -12.22485 2005-Feb-03 00:00 * 267.02253 -23.64522 105.89 -3.74 2.0208234650 -12.25434 2005-Feb-04 00:00 * 267.79155 -23.67107 105.99 -3.19 2.0135494278 -12.28309 ************************************************************************************* Column meaning: TIME Prior to 1962, times are UT1. Dates thereafter are UTC. Any 'b' symbol in the 1st-column denotes a B.C. date. First-column blank (" ") denotes an A.D. date. Calendar dates prior to 1582-Oct-15 are in the Julian calendar system. Later calendar dates are in the Gregorian system. The uniform Coordinate Time scale is used internally. Conversion between CT and the selected non-uniform UT output scale has not been determined for UTC times after the next July or January 1st. The last known leap-second is used over any future interval. NOTE: "n.a." in output means quantity "not available" at the print-time. SOLAR PRESENCE Time tag is followed by a blank, then a solar-presence symbol: '*' Daylight (refracted solar upper-limb on or above apparent horizon) 'C' Civil twilight/dawn 'N' Nautical twilight/dawn 'A' Astronomical twilight/dawn ' ' Night OR geocentric ephemeris LUNAR PRESENCE The solar-presence symbol is immediately followed by a lunar-presence symbol: 'm' Refracted upper-limb of Moon on or above apparent horizon ' ' Refracted upper-limb of Moon below apparent horizon OR geocentric ephemeris R.A._(J2000.0)_DEC. = J2000.0 astrometric right ascension and declination of target. Corrected for light-time. Units: DEGREES dRA*cosD d(DEC)/dt = The rate of change of target apparent RA and DEC (airless). d(RA)/dt is multiplied by the cosine of the declination. Units: ARCSECONDS PER HOUR delta deldot = Target apparent range ("delta") and range-rate ("delta-dot") relative to observer. Units: AU and KM/S ************************************************************************************* -
thx, das ist was ich brauche
Eine Frage habe ich noch, nämlich zur Geschwindigkeit. Warum verwendet man diese komische Einheit:
" dRA*cosD d(DEC)/dt =
The rate of change of target apparent RA and DEC (airless). d(RA)/dt is
multiplied by the cosine of the declination. Units: ARCSECONDS PER HOUR"welchen sinn hat es die Winkelgeschwindigkeit der Deklination mit dem mit dem cosinus der Deklination zu multiplizieren?
-
ok, eliptikale länge und breite habe ich jetzt gefunden sogar den abstand und die geschwindigkeit der Abstandsänderung. Allerdings ist das nur eine Komponente des Geschwindigkeitsvektors.
-
der Bezugspunkt für eine Achse muss wohl der Frühlingspunkt sein. Alle anderen Ergeben sich ja durch rechte Winkel. Ich denke dass ich die Koordinaten aus ekleotikaler breiter, ekliptikaler länge und Enfternung von der Sonne errechnen muss. Allerdings weiß ich nicht, welche Angaben das von dem Zeug das es auf der Seite gibt ist.
-
Hi, ich schreibe gerade eine Software zur nummerischen Berechnung von Ephemerieden. Allerdings brauche ich die Ausgangsdaten (Position und Geschwindigkeitsvektoren) in Form Karthesischer Koordinaten mit der Sonne als Ursprung (ok, theoretisch ist es egal, man kann auch die Erde als Bezugspunkt wählen, allerdings kommen dann wohl keine Kreisbahnen raus). Als Datenquelle würde ich http://ssd.jpl.nasa.gov/horizons.html benutzen. Allerdings weiß ich nicht wie ich daraus Karthesische Koordinaten bekomme/rauslese
