Lunar Ephemeris Polynomials - Revision history

imported>Irnenginer: /* Apollo 7 Example */

2008-03-14T16:03:43Z

Apollo 7 Example

@@ Line 76: / Line 76: @@
 <math>TIMEMO(T_m) = 2440147 JD</math>
-First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.
+First a set of lunar state vectors for the time period during the mission is obtained from JPL’s horizions system.
 [[Category:Apollo Guidance Computer]]

imported>Irnenginer at 15:45, 14 March 2008

2008-03-14T15:45:17Z

← Older revision		Revision as of 15:45, 14 March 2008
Line 77:		Line 77:

	First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.		First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.

			[[Category:Apollo Guidance Computer]]

imported>Irnenginer: /* Apollo 7 Example */

2008-03-07T23:05:13Z

Apollo 7 Example

← Older revision		Revision as of 23:05, 7 March 2008
Line 73:		Line 73:

	<math>TEPHEM(T_0) = 2440140 JD</math>		<math>TEPHEM(T_0) = 2440140 JD</math>

	<math>TIMEMO(T_m) = 2440147 JD</math>		<math>TIMEMO(T_m) = 2440147 JD</math>

	First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.		First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.

imported>Irnenginer: /* Apollo 7 Example */

2008-03-07T23:04:48Z

Apollo 7 Example

← Older revision		Revision as of 23:04, 7 March 2008
Line 72:		Line 72:
	Based on the mission parameters, the following time parameters are chosen.		Based on the mission parameters, the following time parameters are chosen.

	<math>TEPHEM(T_0) = 2440140 JD		<math>TEPHEM(T_0) = 2440140 JD</math>
	TIMEMO(T_m) = 2440147 JD</math>		<math>TIMEMO(T_m) = 2440147 JD</math>

	First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.		First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.

imported>Irnenginer: /* Modern Method of Computing Polynomials */

2008-03-07T23:03:16Z

Modern Method of Computing Polynomials

← Older revision		Revision as of 23:03, 7 March 2008
Line 62:		Line 62:
	In addition the posistion vector data is put in the form		In addition the posistion vector data is put in the form

	<math>X=~~X_d100)~~2^-^3^1</math>		<math>X=X_d1000*2^-^3^1</math>

	Where Xd is the position data in all three axis.		Where Xd is the position data in all three axis.

imported>Irnenginer: /* Modern Method of Computing Polynomials */

2008-03-07T23:01:48Z

Modern Method of Computing Polynomials

← Older revision		Revision as of 23:01, 7 March 2008
Line 54:		Line 54:

	For the paremters selected the one that is required is to use Earth Mean Equator and Equinox of Reference Epoch coordinate system. This cooresonds to the universal coordinate system of the AGC.		For the paremters selected the one that is required is to use Earth Mean Equator and Equinox of Reference Epoch coordinate system. This cooresonds to the universal coordinate system of the AGC.

			The Data will then need to be scaled to the AGC so that the resultant polynomials are valid. First time is scaled and putinto the reference based on TIMEMO such that:

			<math>t=(t_d-t_m)8.64E6*2^-^2^6</math>

			where td is the Julian Date Data and tm is the Julian Date of the mission center point, also to be used for TIMEMO.
			In addition the posistion vector data is put in the form

			<math>X=X_d100)2^-^3^1</math>

			Where Xd is the position data in all three axis.

			Then from the data a set of polynomials can be determined used the least squres method. The resultant values are then stored in twos complement octal in earsable memory, without any additional scaling.

			== Apollo 7 Example==

			Based on the mission parameters, the following time parameters are chosen.

			<math>TEPHEM(T_0) = 2440140 JD
			TIMEMO(T_m) = 2440147 JD</math>

			First a set of lunar state vectors for the time period during the mission was obtained from JPL’s horizions system.

imported>Irnenginer: /* Modern Method of Computing Polynomials */

2008-03-07T22:05:42Z

Modern Method of Computing Polynomials

← Older revision		Revision as of 22:05, 7 March 2008
Line 30:		Line 30:

	The first step is to obtain data on the location of the moon over the period of time of a lunar mission. This data is avalible via the JPL Horizons system [http://ssd.jpl.nasa.gov/?horizons]. The following is an example header of the output.		The first step is to obtain data on the location of the moon over the period of time of a lunar mission. This data is avalible via the JPL Horizons system [http://ssd.jpl.nasa.gov/?horizons]. The following is an example header of the output.
			<pre>Ephemeris / WWW_USER Thu Jun 14 08:22:42 2007 Pasadena, USA / Horizons
			*******************************************************************************
			Target body name: Moon (301) {source: DE405}
			Center body name: Earth (399) {source: DE405}
			Center-site name: BODY CENTER
			*******************************************************************************
			Start time : A.D. 1968-Oct-01 00:00:00.0000 CT
			Stop time : A.D. 1968-Oct-30 00:00:00.0000 CT
			Step-size : 30 minutes
			*******************************************************************************
			Center geodetic : .000000000,.000000000,.00000000 {E-lon(deg),Lat(deg),Alt(km)}
			Center cylindric: .000000000,.000000000,.00000000 {E-lon(deg),Dxy(km),Dz(km)}
			Center radii : 6378.1 x 6378.1 x 6356.8 km {Equator, meridian, pole}
			Output units : KM-S
			Output format : 01
			Reference frame : ICRF/J2000.0
			Output type : LT+S CORRECTED cartesian states
			Coordinate systm: Earth Mean Equator and Equinox of Reference Epoch
			*******************************************************************************
			JDCT , , X, Y, Z,
			*******************************************************************************
			</pre>

			For the paremters selected the one that is required is to use Earth Mean Equator and Equinox of Reference Epoch coordinate system. This cooresonds to the universal coordinate system of the AGC.

imported>Irnenginer at 05:38, 7 March 2008

2008-03-07T05:38:51Z

← Older revision		Revision as of 05:38, 7 March 2008
Line 22:		Line 22:
	Scaling the time and distance differently maximizes the accuracy of the resultant polynomials. The time is scaled to cs2^-26 and the distance is scaled to meters2^-31. This scaling produces polynomials that maintain a +/-1 result constraint. The change in scaling does however, reduces the accuracy of the distance, though but does not effect it significantly. The time scaling does however limit the effective time range to approximately 15 days.		Scaling the time and distance differently maximizes the accuracy of the resultant polynomials. The time is scaled to cs2^-26 and the distance is scaled to meters2^-31. This scaling produces polynomials that maintain a +/-1 result constraint. The change in scaling does however, reduces the accuracy of the distance, though but does not effect it significantly. The time scaling does however limit the effective time range to approximately 15 days.

	Given the time range constrains and the computational constants at the time, the ~~Ephemeredes~~ Data was centered to a time independent of the mission and time of launch. Assuming that the computation of the polynomials was not trivial at the time, and the that the polynomials were computed far ahead of mission with this constraint, the polynomials were based on point in time independent of other time related mission parameters. While most time commutations used the mission time based off the time of launch. The lunar ~~ephemeredes~~ were based on a point during the mission. This point was selected in such a way that placed the negative extreme –(2^-26)/2 cs at a point slightly prior to the launch time. This method maximizes the time range that the polynomials are valid, and allows for changes in the launch time without the need to recompute the polynomials.		Given the time range constrains and the computational constants at the time, the Ephemerids Data was centered to a time independent of the mission and time of launch. Assuming that the computation of the polynomials was not trivial at the time, and the that the polynomials were computed far ahead of mission with this constraint, the polynomials were based on point in time independent of other time related mission parameters. While most time commutations used the mission time based off the time of launch. The lunar ephemeris were based on a point during the mission. This point was selected in such a way that placed the negative extreme –(2^-26)/2 cs at a point slightly prior to the launch time. This method maximizes the time range that the polynomials are valid, and allows for changes in the launch time without the need to recompute the polynomials.


			== Modern Method of Computing Polynomials ==

			Computing the polynomials can be accomplished using freely available data and any number of mathematical computational software.

			The first step is to obtain data on the location of the moon over the period of time of a lunar mission. This data is avalible via the JPL Horizons system [http://ssd.jpl.nasa.gov/?horizons]. The following is an example header of the output.

imported>Irnenginer at 05:18, 7 March 2008

2008-03-07T05:18:49Z

← Older revision		Revision as of 05:18, 7 March 2008
Line 10:		Line 10:
	<math>X(t) &= x_0+x_1t+x_2t^2+\cdots x_9t^9 </math>		<math>X(t) &= x_0+x_1t+x_2t^2+\cdots x_9t^9 </math>

			Where X is the position of the moon with respect to the earth at a given time. There is a complete set of polynomials for each axis X,Y, and Z.

	== ~~Determination~~ of ~~Polynomials~~ ==		Time for the polynomials are expressed as the difference between the mission time and a predetermined time that was defined as a point midway through the lunar mission. This time, called TIMEMO is the zero intercept the polynomial expression.

			== Constraints of the AGC ==

			Aspects of the AGC presented an interesting challenge to programmers. Given that the polynomials would be expressed as double precision numbers in the AGC that had to fall within the range of -1 to +1. In addition the input time and output distance would need to be scaled to the AGC defaults. And finally the calculations would need to be performed to maximize the accuracy of the Ephemeredes.

			The scaling of the parameters is the most important aspect of the computation. The AGC measures distance primary in meters2^-29, and time in cs(centiseconds)2^-28. In order to produce accurate polynomials, the data that is used to produce them must be scaled the same as the AGC. Unfortunately, producing polynomials with this scaling does not produce polynomials within the +/-1 range.

			Scaling the time and distance differently maximizes the accuracy of the resultant polynomials. The time is scaled to cs2^-26 and the distance is scaled to meters2^-31. This scaling produces polynomials that maintain a +/-1 result constraint. The change in scaling does however, reduces the accuracy of the distance, though but does not effect it significantly. The time scaling does however limit the effective time range to approximately 15 days.

			Given the time range constrains and the computational constants at the time, the Ephemeredes Data was centered to a time independent of the mission and time of launch. Assuming that the computation of the polynomials was not trivial at the time, and the that the polynomials were computed far ahead of mission with this constraint, the polynomials were based on point in time independent of other time related mission parameters. While most time commutations used the mission time based off the time of launch. The lunar ephemeredes were based on a point during the mission. This point was selected in such a way that placed the negative extreme –(2^-26)/2 cs at a point slightly prior to the launch time. This method maximizes the time range that the polynomials are valid, and allows for changes in the launch time without the need to recompute the polynomials.

imported>Irnenginer at 05:03, 7 March 2008

2008-03-07T05:03:42Z

← Older revision		Revision as of 05:03, 7 March 2008
Line 8:		Line 8:
	The lunar ephemeris are stored in the AGC using a set of 3 9th order Taylor series of the form:		The lunar ephemeris are stored in the AGC using a set of 3 9th order Taylor series of the form:

	<math> X(t) &= x_0+x_1t+x_2t^2+\cdots x_9t^9 </math>		<math>X(t) &= x_0+x_1t+x_2t^2+\cdots x_9t^9 </math>


	== Determination of Polynomials ==		== Determination of Polynomials ==