## Info

DECLINATION UNCLES «BEG«

Fig. II-5. Right Ascension Versus Declination Plot Including Both Solutions From Each Ambiguous Pair. Concentration of points near declination of -8 deg indicates correct solution should be in that region.

183. &, 183.0 182.S 162.0 181. S 181.0 180. S 180.0 179.s 179.0 176.S

0ECUNBTI0» UNCLES "DEC"

Fig. 11-6. Right Ascension Versus Declination Plot Keeping Only Selected Solution From Each Ambiguous Pair for the Same Data as in Fig. 11-5. Note the greatly reduced range on the axes relative to Fig. 11-5.

A third method is to use each of the solutions in each ambiguous set as trial solutions and choose the one which provides the best fit to the data. Thus, the first solution in the block is arbitrarily selected to provide the trial attitude and used to resolve the ambiguities. The set of remaining attitude solutions is averaged to obtain an initial attitude estimate. Because erroneous attitudes due to noisy or biased data or unrejected terminator crossings may be present, a residual edit process is performed. The residual, for each attitude is defined as the angle between the individual attitude vector and the average attitude vector. The Rt values are used to compute a standard deviation, o, for a set of selected attitudes, according to

v where Wf is the weight assigned to each selected attitude solution. The summation is over all the selected attitudes in the block. The selected attitudes are then compared with the average attitude, and any that differ in arc-length separation by more than No are rejected. (N is normally chosen in the range 3 to 5.) A new averaged attitude is calculated from the edited group of attitudes. This editing procedure is repeated until no additional attitudes are rejected. A goodness-of-fit parameter associated with the averaged spin vector and the M attitudes remaining is defined as o/A/. In a similar manner, other ambiguous solutions present in the original unresolved set of attitudes are selected as trial attitudes and the whole procedure is repeated. The averaged spin vector and the set of attitudes associated with the lowest o/M value are selected as the true set of attitude solutions.

After the resolution of ambiguous solutions, the processing is the same for all the methods. The a priori attitude, or the selected trial attitude, is used to process the unresolved block of attitudes. The chosen attitudes are then averaged to obtain a new attitude estimate. The chosen attitudes are residually edited in iterative fashion until a self-consistent set of solutions remains. In practice, the laborious search for a trial attitude is conducted over a small subset of data because processing time is proportional to the square of the number of solutions in the original set. In addition, for a trial search to be successful, it is necessary for at least one attitude in the block to lie close to the final averaged attitude and for the true attitude solution to vary more slowly than the false attitude solutions. A time-varying attitude or large systematic or random errors present in constant attitude data would result in wildly fluctuating attitude solutions and can cause the block averaging process to fail.

Reliability of the Averaged Attitude Solution. The quality of the computed average attitude solution may be evaluated using cither statistical measures or solution plots. As an example of the latter, the plots of right ascension versus declination shown in Figs. 11-5 and 11-6 immediately reveal if the chosen attitudes are clustering about a constant value. These plots display either the selected and the rejected attitude solutions or only the selected solutions. Plots of attitude solutions versus frame number may also be used to search for systematic variations or incorrect editing. For example, Fig. 9-18 in Section 9.4 shows a plot of declination versus frame number before residual editing for real SMS-2 data. The downward spikes are due to a systematic anomaly in the sensor performance; these data should be removed before further processing. A set of suspect attitude solutions would be revealed by a wide scatter in the plotted data. Similarly, poor quality solutions or a processing method giving inconsistent results can be eliminated. (See Section 11.5.)

A second possible method for evaluating the quality of solutions is to compute statistical indicators. For example, standard deviations can be calculated for the sets of solutions associated with each method separately, for an average of the single-sensor methods for each Earth sensor and for an average over all the attitude methods.

113 Single-Axis Attitude Determination Accuracy

### Lily C. Chen James R. Wertz

In this section, we calculate the uncertainty in deterministic single-axis solutions due to both the statistical noise on the data and estimates of any systematic errors which may be present. The purpose of this calculation is both to determine the attitude accuracy available from given measurements and to provide weights for the various data and measurement types as described in Section 11.2.

Attitude uncertainties can be obtained through two different approaches. In the direct calculation procedure discussed in Section 12.3, the attitude uncertainty is obtained directly from the uncertainties of the various observables via the partial derivatives of the attitude parameters with respect to the observables. Although this is adequate for determining the attitude uncertainty for specific values of the observables, it provides little insight into the underlying causes of the attitude uncertainty and does not lend itself to mission or maneuver planning or analytic attitude studies, where a wide range of alternatives are considered.

An alternative method, described in this section and in Section 11.4, is to express the attitude uncertainty in terms of three factors (involving the various partial derivatives) which have well-defined physical and geometrical interpretations. We then use this factorization to develop explicit analytic expressions for the attitude uncertainty and to discuss the geometrical causes of large uncertainties to provide the perspective necessary for prelaunch analysis and mission or maneuver planning.

In this section, the discussion is restricted to single-axis attitude determined from measurements taken at a single time. The attitude, which corresponds to a point on the celestial sphere as discussed in Chapter 10, is defined as the spin axis for a spinning spacecraft or as the direction of some convenient axis fixed in the body for a three-axis stabilized spacecraft.

If the uncertainties in two measurements are due to independent error sources such as random noise or unrelated systematic errors, then the uncertainties are uncorrelated. In other words, an error in one measurement does not imply any error in the other measurement, and vice versa. Alternatively, part of the uncertainties in two measurements may come from the same error source. For example, if the attitude is determined by measuring the Earth width with a horizon telescope at two different times, then a misalignment in the sensor mounting angle would cause an error in both measurements and the uncertainties weT correlated. Measurement uncertainties are also correlated whenever there exists a systematic error which can introduce uncertainties in both measurements.

113.1 Attitude Accuracy for Measurements With Uncorrelated Uncertainties

Quantized Measurements. The easiest measurements to interpret physically are quantized measurements, i.e., those for which the measurement uncertainty is the result of the step size or bucket size in which the measurements are made or transmitted.* As discussed in Section 10.4, two measurements, m and n, each imply that the attitude lies within a band on the celestial sphere as shown in Fig. 11-7. Here, a plane geometry approximation is made such that the constant measurement curves on the celestial sphere are approximated by straight lines. Lm is the locus of attitudes corresponding to measurement m; ALm = Um/dm is the change in the attitude perpendicular to Lm due to the uncertainty, Um (1/2 step size): dm is the measurement density for measurement m as defined in Section 10.1, and &m/n is the correlation angle between measurements m and n, as defined in Section 10.4.

The correct attitude solution may be anywhere inside the error parallelogram shown in Fig. 11-7. The probability of the correct attitude being in any small area of the parallelogram is the same regardless of the location within the parallelogram.

* In a strict sense, all attitude measurements are quantized by the process of transmitting them as binary numbers. However, we will regard the quantization as important only when the quantization step is sufficiently larger than the noise such that in a continuous string of measurements of an approximately constant observable, only one or two discrete values are reported.

For example, the probability of the attitude being in some small area at the tip of the parallelogram is the same as the probability of the attitude being in an equal area at the center of the parallelogram. The probability density, or probability per unit area on the celestial sphere, is constant inside the parallelogram and zero outside the parallelogram.

The attitude uncertainty for quantized measurements is the error parallelogram. To fully specify this uncertainly requires both the size and orientation of the parallelogram, which depend on four independent parameters. For example, we could give the width of both bands and the azimuthal orientation of each band relative to an arbitrary reference direction.

It is frequently convenient to characterize the attitude uncertainty by a single number. Clearly, this cannot be done in any precise sense because no one number completely defines the error parallelogram. We define three convenient error parameters which may be used depending on the nature of the uncertainty requirements. The component uncertainty is the distance from the center to the edge of the parallelogram along some specified direction, e.g., right ascension uncertainty. The maximum uncertainty, Umax, is the semilength of the longest diagonal, or, equivalently, the radius of a circle cimcumscribed about the parallelogram.

Vmax= |sin^/n| [(ALg,)2+(ALn)z + 2(ALw)(ALn)|cos0gl/n|]l/2

If we have not specified otherwise, U^ will be taken as the attitude uncertainty for quantized measurements. Finally, we define the mean uncertainty, Umean, as the radius of a circle with area equal to that of the parallelogram:

ALmOLn V" VmVn m n

Continuous Measurements. If the uncertainty in a transmitted measurement is due to either Gaussian-distributed random noise or any unknown systematic error which is assumed to have a Guassian probability distribution, then the attitude uncertainty corresponds to an error ellipse on the celestial sphere. For illustration, we first consider the simplest case in which the two independent measurements, m and n, correspond to attitude loci which are orthogonal on the celestial sphere, as shown in Fig. 11-8. Let x be the attitude component perpendicular to L„ and ax be the standard deviation in x resulting from the uncertainty in m\ i.e., ax= Vm/dm, where Um is now the standard deviation of the measurement m. By the definition of a Gaussian distribution, the probability of the jr-component of the attitude lying between x and x + Sx is given by p(x)6x

72 ol)8x

where p(x) is called the probability density for x. Similarly, if y is the attitude component perpendicular to Ln, then the /-component probability is p(y)8y f/lol

The probability both that the x-component lies between x and x + Sx and that the/-component lies between / and y+fy is

gv Fig. 11-8. Error Ellipse for Orthogonal Measurements

gv Fig. 11-8. Error Ellipse for Orthogonal Measurements

p(x,y) is the two-dimensional probability density on the celestial sphere. From Eq. (11-18) it is clear that the lines of constant probability density are ellipses defined by (x2/Oj)+(y2/a*)= K, where K is a constant. The standard deviations of x and y are just the semimajor and semiminor axes of the K= 1 ellipse. As shown in Fig. 11-8, the standard deviation of any arbitrary component, x', is the perpendicular projection of the K= 1 ellipse onto the x' axis. That is,

where v is the angle between the x' axis and the major axis of the ellipse. The same relationship holds for the Ka uncertainties in any attitude component.

In general, we would like to consider the independent measurements, m and n, corresponding to nonorthogonal loci, as shown in Fig. 11-9. For computation, we choose an orthogonal coordinate system, x and y, for which the /-component is perpendicular to L„. Thus, the standard deviation for the/-component is

The standard deviation for the x-component is now more complex. As we will show later (Eq. (ll-27b):

## Post a comment