If the matrices C1,C2 and C0=C1+C2 are positive-definite it is possible to compute a linear proposed identification; but is the linear formalism a reasonable approximation of the full nonlinear problem?
This difficulty is not a separate one from the bad conditioning one:
if the conditioning number of C is very large, there is a very large
eigenvalue
and a very small one .
Then the
inverse matrix
must have a small eigenvalue
and
a large eigenvalue
,
and a corresponding
eigenspace along which the RMS uncertainty is ;
this is also
the length of the longest semi-axis of the confidence ellipsoid
defined by the inequality
The problem of nonlinearity does not have an easy solution, applicable to the orbit identification problem; ``brute force'' (that is computationally intensive) solutions can be applied to analyze a couple of orbits for an already suspected identification, but cannot be used to propose identifications, that is to select among hundred of millions couples of orbits.
One possibility is to use a reduced covariance matrix, as in our tests in Section 4. Another solution is to discard the couples of orbits X1 and X2 such that the difference is too large, by using a simple metric , even when the linear identification algorithm would indicate the couple as a possible identification. This solution is quite unsatisfactory, because one of the goals of the algorithm we are proposing is to allow identification of couples even with large ; the cases with small could be found, and indeed are at present found, with much simpler algorithms, such as the one used in [Sansaturio et al. 1996].
An additional problem of nonlinearity is introduced if a singular set of elements is used. For example, if the least squares fit is performed with the state vector X expressed in the usual keplerian elements , the linear approximation breaks down at the e=0 and/or I=0 singularities of the keplerian elements. If the confidence region contains e=0 and/or I=0, then the linear approximation necessarily fails within the confidence region, even when the latter is small. This problem is easily solved by using nonsingular elements such as the equinoctial ones, for which the covariance propagation is more regular; see Paper I, Section 4.1.