Inverse Problem of Astrodynamics

We consider the problem of determining the center of mass of an unknown gravitational body, using the disturbances in the motion of observed celestial bodies. In this paper an universal approach to obtain the approximate and stable estimate of problem solution is suggested. This approach can be used in other fields of Science. For example, it can be applied for investigation of interactions between fields of forces and elementary particles using known trajectories of elementary particles motions.


Introduction
In 1843-1845 famous astronomers and mathematicians Urbain Jean Le Verrier (1811-1877) and John Couch Adams (1819-1892) independently of one another performed the mathematical research and came to the conclusion that the Solar system includes a celestial body (at least one) which has not observed earlier.
In fact not only the existence of a previously unknown planet has been proven, but also its orbit has been determined with an accuracy, which was sufficient for its detection and surveillance.The planet Neptune has been discovered as result.Mathematically these problems belong to the category of inverse problems of mathematical physics, i.e., to ill-posed problems.The solution of this problem was executed by the method of least squares using some hypothesizes.After discovering Neptune, Le Verrier started the recalculation of theory motion of Uranus by taking into account the motion of Neptune.After finishing his investigation, Le Verrier was able to achieve results with high accuracy, which unfortunately still disagreed with results obtained by observation.This difference was not due to an error in theory or observation [1]- [3].Theoretical success of Adams and Le Verrier was attracted thousands of enthusiasts and professional astronomers and mathematicians.Hundreds of scientific calculations were published based on the work of Le Verrier and Adams, but the results yielded nothing.There were studies, which simultaneously proved the existence of up to two dozen new celestial bodies [3].
Much later it was discovered a natural property of inverse problems their instability.Of course, our knowledge of the Solar system has not been and will never be final and the level of our knowledge is entirely determined by the level of theoretical and observational studies.However, a theoretical analysis of the constructed motion of large (and primarily external) planets indicates, that there are yet unexplained discrepancies between theory and observation.Despite the fact that the theoretical parameters of motion were refined with results made from observation, which were made over a long period of time.For example, there are latitudinal variations in the motion of Uranus and Neptune and the deviation in the movement perihelion of Halley's comet, that cannot be explained by gravitational forces of known solar system bodies.These circumstances have led to the fact that in the 60s of the last century a hypothesis for the existence of a tenth planet emerged.This 10th planet should have a mass equal to the mass of Jupiter, with an approximate distance to the sun of 60 AU and an orbital tilt of 120 i =  .The joint solution of the equations of motion of known planets and of a hypothetical planet and subsequent thorough review of photographic plates of "suspicious" parts of the sky, have not yielded any positive results.Although, according to preliminary estimates, the hypothetical planet was supposed to be 13 -14th magnitude, and on photographic plates were considered objects to 16.5 magnitude, but the tenth planet was not found.
Analysis of solution methods of Le Verrier and Adams shown these methods did not take into account the inaccuracy of mathematical model of planets motion.The success of their investigations was guaranteed with help of right hypothesis about tilt of unknown planet orbit to plane of the ecliptic and orbit eccentricity.
This fact is explained failure of big numbers of investigations after Le Verrier for searches of planet Pluto.Thus the development of stable methods of approximate solutions of the inverse problem astrodynamics in more general statement remains relevant.

Statement of the Inverse Problem Astrodynamics
We consider n interacting masses moving under the forces of mutual attraction in an inertial coordinate system.Masses i m are given the index ( ) , ik r denotes the vector connecting the mass i m and the mass k m .According to Newton's law the resulting force j F acting to the mass j m is equal [4] 3 1 where G is the gravitational constant.The mass j m under the influence of this force executes the motion which satisfies the differential equation ( ) where oj r is the radius vector connecting the origin of inertial coordinate system with the mass j m .Let us make the transition of the variables , It is assumed that among n gravitational masses the location of only mass n m is unknown.The last term in the sum on the right of Equation ( 3) is uncertain.
Equation (3) takes the form where ( ) ( ) is the function to be determined.
In terms of projections on the axis of the inertial coordinate system the Equation (4) can be written in the form: ( ) ( )  5)-( 7) twice from 0 t to t: Each equation of the system Equation ( 8) can be presented in the form where Equations ( 9) are known as Volterra integral equations of the first kind with respect to the unknown functions . By finding the solutions of the Equations ( 9) you can restore the force vector ( ) , exerted by the mass n m on the mass j m up to a constant factor.Performing similar calculations and solving equations of the type Equation ( 9) for the mass l m , we can de- termine (up to a constant factor) the force ( ) , which is acting on the mass l m from the mass n m .The in- tersection of the lines of action of the vectors ( ) gives the position of the mass n m in space (in a chosen inertial system).
As is easily seen, ( ) jk u t are defined by the functions , which are assumed to be known from astronomical observations of the motions of the masses  9) are an ill-posed problem [5].
In the Equations ( 5)-( 7) of motion the coefficients G are determined from astronomical observations and experimental investigations and so these values are known only approximately.Thus it is assumed that each coefficient in Equations ( 5)-( 7) is in some interval: Introduce into consideration the following notations where 1 The Inequalities (10) define a closed region of D in the ( ) , which can be equipped with following norm [6]: The Equation ( 9) can be written in the form where the operator A  is a compact operator, : Suppose instead ( ) 11) is given by the approximate function The approximate value of u  is the right side of the Equation (11).So Let us estimate the deviation of u  of ( ) ex u t , assuming that the exact operator ex B is linear: where [ ] ( ) Since real processes can be described by mathematical methods only approximately.It is assumed that the exact operator ex A in the Equation (11) differs from the approximation of the operator A  (if the exact opera- tor is linear) by a predetermined amount In this case it is possible to use the algorithm for solving the inverse problem with approximate operator A  proposed in [7].

Y. Menshikov
However, the assumption of linearity of the exact operator ex A and of information relatively of size of h in most cases does not correspond to reality.
Then the set of possible solutions 0 ,h Q δ of Equation ( 11) with account of the linearity of operators , p A B  will have the form: It is easy to show that if the operator A  is a compact operator, then the set 0 ,h Q δ is an unbounded, closed and convex set for any ( ) ( ) Then it is assumed that exact solution ex f of Equation (11) be- longs to functional space [ ] For solving ill-posed problem Equation (11) we use the Tikhonov regularization method with the stabilizing functional where 0 1 0, 0 q q ≥ > , 0 1 , q q are constants.Thus, the problem of finding an approximate solution of Equation ( 11) reduces to the solution of the extreme problem: It should be noted that there is no way to determine the size of 1 d , since the exact operator ex B is the un- known operator.Moreover, the exact operators , ex ex A B cannot be constructed, in principle, because the ma- thematical models, which describe the real processes, are always approximations.Of course, under certain assumptions about the exact operators , ex ex A B we can get some error estimates of 1 , d h , but such estimates are unrealistic.
Therefore, the approximate solution of inverse problems of measurement are not of interest for practical use due to instability of the solution.
The way out of this impasse exists, if by the investigation of inverse problems of measurement restrict only some estimates of exact solutions.

The main Hypothesis and Results
The set of possible solutions , p Q δ of Equation (11) for fixed operators , p A B  has the form: The regularization parameter α can be find from condition: To obtain useful information on the exact solution of the inverse problem the use of the following hypothesis (main hypothesis) is suggested: the inequality is valid where the function ex f is an exact solution of the inverse problem of measurement Equation ( ), or its absence but without guarantee (in case 0 , it is also possible to have a real celestial body n m as this does not take into account all the oper- ators , To study the effect of the parameters of mathematical model , we need to examine all possible estimates p f  .Note that the preparation of estimates does not use properties of the exact operators , ex ex A B .To assess the existence of an exact solution, it is possible to use the stabilizing functional [ ] f Ω .We give sufficient conditions for the existence of an element p f  . be continuous, non-negative and strictly convex on 1 F .Then the solution p f  of the extreme problem (15) exists, is unique for any 0 δ > and any In order to study the influence of the process parameters on the estimation of the exact solution, we consider the union of the sets , p Q δ for all vectors The set un Q is unbounded as a union of unbounded sets.Let un f be a solution to the extreme problem: [ ] The regularization parameter α was determined by the method of discrepancy: ( ) 15) will yield more accurate estimates for the exact solution.This approach is based on the work of [8].In addition, it is obvious that the inequality is valid.Evaluation of the exact solution of un f allows us to conclude on the existence of an unknown planet with a guarantee in case 0 un f ≠ or its absence, but without guarantee in case 0 un f ≡ .To study the influence of the parameters p on the estimation of the exact solution of the inverse problem it is necessary to have a possibility to select an operator 0 p B from operators 0 p B which satisfies the following condition: for any possible R and any Subsequently, the operator 0 p B in the right-hand side of Equation (11) will be called as "special minimal operator" in the sense of the Inequality (25).
If the operator 0 p B exists and is unique, then the problem of finding the greatest lower bound of the func- tional [ ] f Ω on the set un Q will have a solution, which coincides with the solution of the more simpler ex- treme problem [9]: find an element The problem (26) has a solution for every D ∈ p and δ , since the conditions of the Theorem 1 are satis- fied.
In this case the following inequality holds for any Evaluation of the exact solution 0 f allows us to conclude the existence of an unknown planet with a guar- antee in case 0 un f ≠ or its absence, but without guarantee in case 0 un f ≡ .Theorem 2. Special minimal operator 0 p B in the Equation (11) exists, is unique and corresponds to the vector ( ) be a realization of an astronomical observations.Consider the problem of determining the exact lower bound the functional [ ] where C is a real symmetric matrix Coefficients of the matrix C are given by: , then the Silvester's inequalities gives us:

[ ]
Ω p is strongly convex on D .Similarly, [10] we prove that [ ] Ω p achieves the greatest lower bound in a single point in the domain ( ) for any ( ) t R .□ Suppose that among the operators p B there is an operator Thus the following inequality is valid: Theorem 3. Special maximal operator 1 p B in the Equation (11) exists, is unique and corresponds to the vector ( )

Conclusion
In this paper we proposed an algorithm for finding the coordinates of an unknown gravitational mass as a result of astronomical observations.This problem was solved first by Urbain Jean Le Verrier and John Couch Adams.Here a more universal approach was suggested.The proposed hypothesis allows us to exclude the error of the approximate operator from the calculations.Also conditions for the existence of an approximate solution were obtained and several non-standard formulations of inverse problems were considered.Suggested approach can be used in other fields of Science.For example, it can be applied for investigation of interactions between fields of forces and elementary particles by help of known trajectories of elementary particles motions.
a linear restricted operator, which transforms elements of the space [ ] L t T .Let us assumed that the functional space [ ] 0 , F C t T = as according of the physical sense the searching functions must be the continuous functions.The operator A  in Equation (11) doesn't depend on parameters of the mathematical model of the motion of body n m .It is easy to see that the operator A  in Equation (11) is a completely continuous operator [6].The operator p B depends on the specific values of the parameters of the mathematical model of the process, which is p. the exact operators in Equation (11), the exact vector-valued function ( ) t R and the exact function of u in the right-hand side of Equation (11).

1 pBoperator 1 pB
is the inverse operator to A  .Subsequently, the operator in the right-hand side of Equation (11) will be called "special maximum operator" in the sense of satisfying InequalityIf the exists and is uniquely determined then we can consider the solution of the following extreme problem: find an element problem (31) has a solution for any D ∈ p and δ , since the conditions of Theorem 1 are satisfied.
proof is similar to Theorem 2. □