Content deleted Content added
incorrect extra period removed |
|||
(46 intermediate revisions by 34 users not shown) | |||
Line 1:
{{short description|Solution to the spacecraft attitude determination problem}}
The '''TRIAD method''' is the earliest published algorithm for [[Spacecraft attitude determination|determining spacecraft attitude]], which was first introduced by Harold Black in 1964.<ref>{{cite journal|last=Black|first=Harold|title=A Passive System for Determining the Attitude of a Satellite|journal=AIAA Journal|date=July 1964|volume=2|issue=7|pages=1350–1351|doi=10.2514/3.2555|bibcode = 1964AIAAJ...2.1350. }}</ref><ref>{{cite journal|last=Black|first=Harold|title=Early Developments of Transit, the Navy Navigation Satellite System|journal=Journal of Guidance, Control and Dynamics|date=July–August 1990|volume=13|issue=4|pages=577–585|doi=10.2514/3.25373|bibcode = 1990JGCD...13..577B }}</ref><ref>{{Cite journal |last=Markley |first=F. Landis |date=1999 |title=Attitude Determination Using Two Vector Measurements |url=https://www.researchgate.net/publication/4706531_Attitude_Determination_Using_Two_Vector_Measurements |journal=1999 Flight Mechanics Symposium |pages=2 |via=ResearchGate}}</ref> Given the knowledge of two vectors in the reference and body coordinates of a satellite, the TRIAD algorithm obtains the direction cosine matrix relating to both frames. Harold Black played a key role in the development of the guidance, navigation, and control of the U.S. Navy's Transit satellite system at Johns Hopkins Applied Physics Laboratories. TRIAD represented the state of practice in spacecraft attitude determination before the advent of [[Wahba's problem]]<ref>{{cite journal|last=Wahba|first=Grace|title=A Least Squares Estimate of Satellite Attitude, Problem 65.1|journal=SIAM Review|date=July 1966|pages=385–386|doi=10.1137/1008080|volume=8}}</ref> and its several optimal solutions. Covariance analysis for Black's solution was subsequently provided by Markley.<ref>{{cite journal|last=Markley|first=Landis|title=Attitude Determination Using Vector Observations: A Fast Optimal Matrix Algorithm|journal=The Journal of Astronautical Sciences|date=April–June 1993|volume=41|issue=2|pages=261–280|url=http://www.malcolmdshuster.com/FC_Markley_1993_J_FOAM_JAS_MDSscan.pdf|accessdate=April 18, 2012}}</ref>
{{NumBlk|:|
Line 10 ⟶ 11:
|{{EquationRef|1}}}}
for <math> i = 1,2 </math>, where <math> A </math> is
{{NumBlk|:|
Line 18 ⟶ 19:
|{{EquationRef|2}}}}
{{NumBlk|:|
Line 26 ⟶ 27:
|{{EquationRef|3}}}}
where <math> \vdots </math> have been used to separate different column vectors.
The solution presented above works well in the noise-free case. However, in practice, <math> \vec{r}_1, \vec{r}_2 </math> are noisy and the orthogonality condition of the attitude matrix (or the direction cosine matrix) is not preserved by the above procedure.
{{NumBlk|:|
Line 54 ⟶ 55:
|{{EquationRef|7}}}}
to be used in place of the first two columns of equation ({{EquationNote|3}}). Their cross product is used as the third column in the linear system of equations obtaining a proper orthogonal matrix for the spacecraft attitude given by the following:
{{NumBlk|:|
Line 62 ⟶ 63:
|{{EquationRef|8}}}}
While the normalizations of equations ({{EquationNote|4}}) - ({{EquationNote|7}}) are not necessary, they have been carried out to achieve a computational advantage in solving the linear system of equations in ({{EquationNote|8}}). Thus an estimate of the spacecraft attitude is given by the proper orthogonal matrix as
{{NumBlk|:|
<math>
\hat{A} = \left[ \hat{S} ~ \vdots ~ \hat{M} ~\vdots~ \hat{S} \times \hat{M} \right] . \left[ \hat{s} ~\vdots~ \hat{m} ~\vdots~ \hat{s} \times \hat{m} \right]^T
</math>
|{{EquationRef|9}}}}
Note that computational efficiency has been achieved in this procedure by replacing the matrix inverse with a transpose. This is possible because the matrices involved in computing attitude are each composed of a TRIAD of [[Orthonormality|orthonormal]] basis vectors. "TRIAD" derives its name from this observation.
==TRIAD Attitude Matrix and Handedness of Measurements==
It is of consequence to note that the
{{NumBlk|:|
Line 86 ⟶ 87:
\Delta = \left[ \hat{s} ~\vdots~ \hat{m} ~\vdots~ \hat{s} \times \hat{m} \right].
</math>
Note that if the columns of <math> \Gamma </math> form a left
{{NumBlk|:|
<math>
Line 94 ⟶ 95:
This is quite useful in practical applications since the analyst is always guaranteed a proper orthogonal matrix irrespective of the nature of the reference and measured vector quantities.
==See also==
Line 105 ⟶ 106:
{{Reflist}}
[[Category:
[[Category:Rotation in three dimensions]]
|