Assessment of Non-Coplanar Maneuver Parameters and Perturbing Accelerations Using the Minimal Number of Observations

: The algorithms for determining the active space object maneuver parameters in the conditions of near-circular orbits are presented in the paper. The right ascension and declination angles are used to determine the parameters of a single maneuver with transversal and lateral components (the application moment and the velocity impulse magnitude). Two pairs of angles are used to determine the parameters of the maneuver with only a lateral component. Two pairs of angles are needed for the determination of the parameters of the long-duration maneuver performed by a low-thrust engine (moments of the maneuver start and finish, and components of the acceleration delivered by the engine). The essential detail that makes it possible to determine the parameters of long-duration maneuvers is their symmetry relative to the center of the engine operating interval. Considerable perturbing accelerations, which are not accounted for by traditional perturbing models, affect passive objects, especially membranous objects with a big and variable area-to-mass ratio. This decreases dramatically the accuracy of these objects’ motion propagation. In this paper, the magnitude of perturbing acceleration is determined with the assumption that it is constant and is active throughout the whole time interval from the moment of the last orbit determination to the moment of the new measurement used. Examples of the assessment of maneuvers performed by an object in the geostationary orbit are presented.


Introduction
Nowadays, the catalogs of different organizations contain information for more than 25,000 space objects (SOs) in near-Earth orbits.The majority of them move along the slightly perturbed Keplerian orbits, and their motion propagation presents no considerable difficulties.However, there are objects for which motion can be quite challenging to propagate.First of all, this concerns maneuvering space objects (MSOs).More than 5000 of such objects are already in orbit and more are going to come due to the deployment of constellations with hundreds of active satellites.These objects maneuver to maintain their orbits or with the purpose of changing their orbits in order to fulfill new mission tasks.As a rule, only MSO operators have real-time data containing information about alterations of orbit parameters.Moreover, the maneuvers of a number of satellites are being calculated on board.Hence, currently, we have to accumulate enough information for the accurate determination of the new orbit of these objects after the fulfilled maneuvers.The rapid propagation of motion of the maneuvering SO, with consideration of the performed maneuver, will allow us to calculate the collision probability between the protected and the maneuvering objects in a prompt and accurate fashion, operatively determine the MSO orbit changing goal, and assess the accuracy of the fulfillment of the maneuvers performed by our own spacecraft, i.e., assess the correctness of the engine functioning.
Due to its importance, the problem of assessing the maneuver performed by an MSO has been in the process of being solved for years by specialists of different space organizations.Many works on this topic have been published.The works differ by the measurements used (optical or radio locational) and by the types and heights of the MSO orbits.The parameters of the orbit after a maneuver are considered to be known in some papers, and the orbit is used for the assessment of the maneuver; in some works, the orbit is determined simultaneously with the maneuver parameters.Usually, only the fact that the maneuver took place is established in the papers.Paper [1] contains 55 items in bibliography list considering the maneuver assessment, paper [2]-49, and our work [3] contains 26 items.These lists differ from each other almost completely.Despite the fact that not all cited works in these three articles relate to the maneuver assessment topic, their total number exceeds 100.This list is far from being full.A review of the results from more than 100 works is worth conducting in a separate paper, which will be far more substantial than this work, not mentioning that the contents of the aforementioned papers with high quality reviews can be easily accessed.In this light, we will not repeat ourselves and consider study [1], which is one of the most recent works on the maneuver assessment topic.The choice of this work can be also explained by the fact that the problem statement in [1] practically coincides with the problem statement in this work.Hence, it can be agreed that study [1] contains the most recent accomplishments relating to the problem considered in this paper.It is not necessary to review the works considering low Earth orbits and radiolocational measurements in this paper due to the total difference between the problem of interest and these problems.
Paper [1] and this paper have the same basis.Optical measurements are used in both papers.It is assumed that there are series of close optical measurements with the distance of several seconds to several dozens of seconds.The total duration of the series is far less than the orbit period.Such series are denoted as tracklets in [1].The basic problem of the maneuvering objects catalogue maintenance is being solved in the present paper, and it is necessary to determine the MSO orbit after the maneuver and assess the parameters of the maneuver itself with the help of the minimal number of tracklets.Although this problem is not solved in [1] (it is mentioned that it is to be solved in the future), another important problem is solved: the measurement identification problem.Since the geostationary area is heavily inhabited, several objects may occur in one limited area.It is not an easy matter to tell which new measurement corresponds to which object in the area.The identification problem is not solved in this work; it is assumed that measurements are already assigned to the specific objects.But the simplicity of the suggested maneuver assessment and hence the simplicity of the orbit determination after the maneuver allows for easy consideration of the possible variants (there are not many as only several objects may be in one GSO point) and use of the results obtained in the paper for the identification of measurements.It is mentioned several times in paper [1] that is not possible to determine the orbit by one short tracklet.This is correct if there is no information on the orbit before the maneuver.For example, Laplace and Gauss methods are now used to determine the orbit of a newly discovered object.Three short measurements on a relatively short curve are needed for these methods.However, it is shown in the present work that one measurement (two at most) is needed for the determination of the orbit after a maneuver for the maneuvering objects for which the orbits are known.
As it has been already mentioned in a number of works (including the author's works), the orbit after a maneuver is considered to be known and it is used during the maneuver assessment.However, in other works, the orbit after a maneuver is determined simultaneously with the parameters of the maneuver [4][5][6][7].These works are noted due to the use of optical measurements in them.The next step is taken in this paper.The maneuver parameters are assessed by the minimal number of observations without an accumulation of the measurements necessary for the initial orbit determination.Then, the parameters of an orbit shaped by a maneuver are determined with the help of information about the maneuver.The methodology for the assessment of the maneuvers of active SOs in a geostationary orbit (GSO) developed in this paper will allow us to sped up the accurate modeling of the motion of these objects with the use of the minimal number of pairs α, δ (the angles of the right ascension and declination) which set the direction from the observer on the Earth's surface to the points in which MSOs are detected.
The suggested methodology to assess active SO maneuvers will allow us to determine the parameters of the single-impulse maneuvers with different engine orientations, and the maneuvers can be performed both by high-and low-thrust engines.These objects comprise approximately 20% of all objects registered in geostationary orbits.The developed approach to the solution of the problem allows us to also assess the non-modeled perturbations, which are currently not accounted for during the assessment of the orbit of some passive SOs.The motion propagation of approximately 30% of the passive objects of space debris registered in the geostationary orbit is not accurate enough for the calculation of the tolerable possibility of collision with these objects.It happens due to the fact that traditional perturbing models used for the determination of orbits of these objects do not correspond to real perturbations.An algorithm to account for the additional non-modeled perturbations is suggested in this paper.This will allow us to improve the accuracy of the motion propagation of the objects of this type.This paper contains the development of the methodology published in [3], in which only coplanar maneuvers were considered and the perturbing accelerations were not determined.

General Problem Statement
The MSO state vector X = X(t 0 ) is known at time t 0 .The observer's coordinates on the Earth's surface Y = Y(t i ) for t i (t i > t 0 ) and the angles α i , δ i , which point out the directions from the observer to MSO in this moment of time, are also known.We determine the times at which the engine starts and stops (the engine started after t 0 but before t i ), its orientation, the MSO acceleration and the maneuver total delta-v.

Equations Describing the Influence of the Velocity Impulses
In case the burn duration is small enough in regard to the MSO orbit period, the following assumption can be made: the orbit parameters' alteration due to the influence of the burn was instant.
Each impulse, applied at points with the angles φ i (i = 1, . . ., N), causes deviations in the orbital elements at the specified point with the angle φ f , and the sums of these deviations after N impulses can be written as follows [8]: where r 0 , V 0 = the reference circular orbit radius and velocity,∆V r i , ∆V t i , ∆V z i = the radial, transversal, and lateral components of the i-th impulse, respectively, and ∆r, ∆V r , ∆V t , ∆n = deviations by radius, radial and transversal components of velocity and along the orbit deviation caused by the velocity impulses.The angles φ i , φ f are being counted from the MSO position at t 0 to the direction to the MSO movement.The motion occurs in the vicinity of the reference circular orbit.
The linearized equations of unperturbed motion in the vicinity of the circular orbit were used to obtain Equations ( 1)-( 6) [9].

Determination of the Along-the-Orbit Distance and Lateral Deviations and Deviation by Radius, Using the Known Observation Angles α, δ
In order to solve the stated problem one should determine the along-the-orbit distance ∆n i , perpendicular ∆z i deviations and the deviation by radius ∆r i caused by the maneuver using the known angles α i , δ i in the moment t i first.For this sake, the initial orbit X(t 0 ), which is known for t 0 , is propagated from t 0 to t i with the help of numerical integration.The orbit orientation of the uncorrected MSO orbit for t i is determined.
Then, the point (x i , y i , z i ) of intersection between the beam, the direction of which is set by the position of the observer on the Earth's surface, the angles α i , δ i and the plane of the uncorrected orbit of the MSO, is determined.The vector r p (t i ) pointed from the center of Earth to the point of intersection of the beam with the plane of the initial orbit (this vector belongs to the plane of the MSO uncorrected orbit) is determined.Then, ∆r i = r p (t i ) − r orbi is calculated, where r orbi = the magnitude of the radius vector of the uncorrected orbit, pointed along the vector r p (t i ), and r p (t i ) = vector r p (t i ) magnitude.The along-the-orbit distance ∆n i between the position of an MSO in the uncorrected orbit in the moment t i (set by vector r(t i )) and its position in the orbit shaped by the maneuver in this moment (set by r p (t i )) is calculated analogously.The angle between vectors r(t i ) and r p (t i ) is multiplied by r 0 in order to calculate ∆n i .
In order to obtain the lateral deviation, we use the fact that the lateral velocity impulse does not change the orbit geometry.It changes only the orbital plane orientation.Hence, the distance from the Earth to the point of intersection of the beam with the new orbit plane after the maneuver coincides with the distance from the Earth to the MSO in the uncorrected orbit in this moment.Thus, the magnitude r(t i ) of the vector r(t i ) is used for the calculation of the lateral deviation.Then, the distance L from the position of the observer on the Earth's surface to the end of the vector r(t i ) is found.Vector R(t i ) corresponds to the point on the beam, the distance from which to the center of the Earth is r(t i ).This beam is set by the position of the observer and the angles α, δ.The distance from the observer to the abovementioned point is approximately L. We obtain the lateral deviation ∆z(t i ) by calculating the deviation vector ∆R = R(t i ) − r(t i ) and by projecting ∆R on the perpendicular to the initial orbit plane (effective for the moment t i ).

Possible Variants of the Problems Solved
Due to the existence of several possible variants regarding the realization of the maneuvers with lateral components, several problems need to be considered and solved.The methods for solving the problems and the number of the used sets of deviations ∆r, ∆n, ∆z(α, δ) differ too.Let us enumerate these problems.

1.
The velocity impulse with only the lateral component has been performed.We determine its application angle φ and the lateral component magnitude ∆V z .There are two unknowns in the problem.Two lateral deviations ∆z 1 , ∆z 2 (two tracklets) are needed for its solution.

2.
The velocity impulse with transversal and lateral components has been performed.We determine its application angle φ and the magnitude of the transversal and lateral components ∆V t , ∆V z .There are three unknowns in the problem, and the deviations ∆r, ∆n, ∆z are used for its solution (only one tracklet is used).

3.
The velocity impulse with transversal, radial, and lateral components has been performed.We determine its application angle φ and the magnitude of the transversal, radial and lateral constituents ∆V t , ∆V r , ∆V z .There are four unknowns in the problem.Hence, two sets of deviations ∆r, ∆n, ∆z are needed for its solution.4.
The long-duration transversal and lateral maneuver has been performed.We determine the maneuver angular duration ∆φ, its middle φ m , the angle between the engine's axis and the orbit plane γ and the acceleration w delivered to the MSO.There are four unknowns.Hence, two pairs of deviations ∆z, ∆n are needed.

5.
The long-duration maneuver with transversal and radial components has been performed.We determine the maneuver angular duration ∆φ, its middle φ m , the angle of rotation of the engine's axis in the orbit plane θ and the acceleration w delivered to the MSO.There are four unknowns in the problem.Hence, two pairs of deviations ∆r, ∆n(α, δ) are needed for its solution.Though the accelerations in the plane of the initial orbit are calculated in this problem, the solution to the problem can be used to determine the perturbing accelerations.This problem has no direct relation to non-coplanar maneuvers, but it will be used while obtaining the non-modeled accelerations.6.
The long-duration maneuver with transversal, radial and lateral components has been performed.It is necessary to determine the maneuver angular duration ∆φ, its middle φ m , the angle of the engine axis's turn in the orbit plane θ, the angle between the engine's axis and the orbit plane γ and the acceleration w delivered to the MSO.There are five unknowns in the problem.Hence, two sets of deviations ∆r, ∆n, ∆z(α, δ) are needed to solve it.
It is worth mentioning that maneuvers in the geostationary area are usually performed for the purpose of correcting the east-west longitudinal drift and the north-south inclination drift.They are performed to maintain the set longitudinal boundaries with the minimal inclination and to resist the perturbing factors of natural Earth gravity and a third body.

Lateral Velocity Impulse Assessment
The application angle φ z and the magnitude of the lateral component ∆V z of the velocity impulse are to be determined.There are two unknowns in the problem, and two deviations ∆z 1 , ∆z 2 calculated for the moments t 1 , t 2 are needed for its solution.
Equation ( 5) is used twice: By excluding ∆V z , we obtain which gives us φ z , which in turn gives us ∆V z .
It is a widespread but simple problem with a guaranteed single solution which can be found fast.It goes without saying that the quality of the solution is affected by the accuracy of calculation of the deviations ∆z 1 , ∆z 2 .

Assessment of the Transversal and Lateral Components of the Velocity Impulse
The equations for the velocity impulse application angle φ and for the determination of the magnitudes of the transversal and lateral components ∆V t , ∆V z can be written as The angle φ f corresponds to the point in which the deviations ∆r, ∆n, ∆z were calculated.
One can obtain ∆V t from the first equation and apply it for the second one.Hence, φ is found and can be used to find ∆V t and ∆V z .

Assessment of the Impulse with Transversal, Radial and Lateral Components
There are four unknowns in the problem: the application angle φ and the magnitudes of the transversal, radial and lateral components ∆V t , ∆V r , ∆V z .Hence, two sets of the deviations ∆r, ∆n, ∆z calculated for the moments of time with which the angles φ 1 and φ 2 correspond are needed for its solution.
The equations for the velocity impulse application angle φ and the magnitudes of the transversal, radial and lateral components of the velocity impulse ∆V t , ∆V r , ∆V z are as follows: The ∆V t and ∆V r are obtained from Equations ( 7) and ( 8) and are later used in (10).We obtain the Equation (13) and can determine φ.
With the known φ, we can obtain ∆V t and ∆V r .Equation (13) may have several solutions for this problem.For this case, the velocity impulse parameters for every solution are used in Equation (9).The final solution for the problem will be the solution for which Equation ( 9) is fulfilled with the best accuracy.
The velocity impulse of the lateral component is determined after the velocity impulse application angle.One should choose between Equation (11) or Equation (12) depending on a bigger value of ∆z.

Shaped Orbit Determination
By altering the velocity in the initial orbit in the moment of the velocity impulse application by the magnitude of the velocity impulse components, we obtain the orbit after the maneuver.

Long-Duration Maneuver Parameters Assessment
The impulsive and the long-duration maneuvers are performed in the GSO.The type of maneuver can be distinguished by its duration in comparison to the period of the orbit in which the SO maneuvers.The maneuver will be impulsive if its duration is far less than the orbit period.When the duration of the maneuver is comparable to the orbit period, a low-thrust engine maneuver can be considered.However, the difference between two types of the maneuvers is nominal and depends on the complexity of the problem being solved.The more complex the problem is, the earlier we need to account for the maneuver duration.It was shown in [8] that a maneuver with the duration of 40 • or less produces an alteration of the same orbit elements as the impulsive maneuver.The effectiveness of the long-duration maneuver decreases (especially for the eccentricity correction [8]) with a longer duration, and the maneuver duration should be accounted for.
In this problem, for the assessment of long-duration maneuvers, there is a need to determine the maneuver's start and finish and the components of the acceleration produced by the engine w t , w r , w z .The symmetry relative to the center of the engine operating interval allows us to use the middle of the interval φ m and its duration ∆φ instead of the beginning and end of the engine operating interval.

Assessment of the Influence of the Transversal and Radial Components of the Thrust Vector on the Deviations Caused by Radius and along the Orbit
Let us first determine the effect from the low-thrust maneuver on the deviation caused by radius in the moment set by the angle φ f .The orientation of the engine is fixed by the transversal component in the orbital coordinate frame.
It is supposed that the maneuver has an angular duration, ∆φ, the middle of the maneuver is defined by the angle φ m , and the acceleration is delivered to the MSO w t .By using (1), one can find the deviation caused by the radius, caused by the velocity alteration ∆V t , equally distributed throughout the whole interval of the argument of latitude ∆φ: By using the equation ∆φ Pr 0 = w c w , w c = centripetal acceleration of the reference circular orbit (w c = V 2 0 r 0 ), and w = engine acceleration (w = P m ), we obtain the final formula for the deviation as follows: Similarly, the influence of the radial component distributed on the interval of the argument of latitude ∆φ on the by-radius deviation is determined as follows: With the help of ( 4), one can also find the effect of the maneuver performed by the low-thrust engine, the orientation of which is fixed by transversal component in the orbital coordinate frame, on the deviation along the orbit.Keeping in mind that the maneuver has an angular duration, ∆φ, the middle of the maneuver is defined by the angle φ m : The influence of the radial component, distributed on the interval of the latitude argument ∆φ, on the deviation along the orbit can be found in the same fashion: where φ f = the point in which the deviations are calculated.

Long-Duration Coplanar Maneuver Parameters Assessment
There are four unknowns in the problem of the coplanar long-duration maneuver φ m , ∆φ, w t , w r .Hence, two pairs of deviations ∆r, ∆n(α, δ) should be used for its solution.By solving the equation system in which the influence of the burns is described by the Formulas ( 14)-(17), we obtain all the unknowns.If it is known that the object is using only the transversal orientation, only three unknowns remain, φ m , ∆φ, w t .Only three deviations are needed: ∆n for the first point and ∆r, ∆n for the second.It is more preferable to use the deviations along the orbit, as they can be calculated more accurately, and the deviations which correspond to the more distant point from the maneuver.

Assessment of the Non-Modeled Accelerations in the Orbit Plane
The Formulas ( 14)-(17) allow us to determine the non-modeled perturbations affecting the non-maneuvering objects, the orbits of which suffer from being determined not accurately enough.It is also supposed that the non-modeled accelerations are caused not by natural causes, which is, in fact, so, but by the work of some engine.It is assumed that the perturbing acceleration is effective throughout the whole interval from the moment of the last SO orbit determination (φ 0 ) till the moment φ f , for which the deviations ∆r, ∆n(α, δ) are being calculated.Thus, we know the angular duration of the maneuver ∆φ = φ f − φ 0 and its middle φ m = φ f − 0.5∆φ.It can be seen that φ f is not used anymore after the substitution of φ m in Equations ( 14)-(17).Only w t , w r are left to determine.It is also agreed that the perturbing acceleration is constant regarding its value, does not change the orientation in the orbital coordinate frame, and is effective in the orbital plane in this problem.Thus, its action can be fully described by Equations ( 14)-(17).It is necessary to determine the transversal w t and radial w r components of the perturbing acceleration.There are two unknowns in the problem.Hence, only one pair ∆r, ∆n(α, δ) should be used for its solution.
By using Equations ( 14)-( 17), we obtain the system of two linear equations and can obtain transversal w t and radial w r components of the perturbing acceleration.The influence of the perturbations constant in the inertial coordinate frame can be calculated in a similar way.These formulae were omitted in order not to overload the understanding of the paper.

Assessment of the Influence of the Lateral Component of the Thrust Vector
In general, the perturbing acceleration can have the lateral component w z , which needs to be obtained too.
Let us determine the influence of the lateral component of the perturbing acceleration w z on the turn of the orbit plane.
Since the maneuvers in the plane and the plane-turning maneuvers do not correlate with each other, they can be considered separately.
The angular alteration of the orbit plane orientation ∆γ, caused by the lateral velocity ∆V z in the angular interval ∆φ, can be calculated using the following formula: where ∆φ = the angular magnitude of the active part, and φ = the angle between the middle of the active part and the current point.The middle of the active part is situated on the line of intersection of the orbit planes (Figure 1).The influence of the perturbations constant in the inertial coordinate frame can be calculated in a similar way.These formulae were omitted in order not to overload the understanding of the paper.

Assessment of the Influence of the Lateral Component of the Thrust Vector
In general, the perturbing acceleration can have the lateral component z w , which needs to be obtained too.
Let us determine the influence of the lateral component of the perturbing acceleration z w on the turn of the orbit plane.
Since the maneuvers in the plane and the plane-turning maneuvers do not correlate with each other, they can be considered separately.
The angular alteration of the orbit plane orientation γ Δ , caused by the lateral velocity z V Δ in the angular interval ϕ Δ , can be calculated using the following formula: where ϕ Δ = the angular magnitude of the active part, and ϕ = the angle between the middle of the active part and the current point.The middle of the active part is situated on the line of intersection of the orbit planes (Figure 1).).
The line of intersection of the orbit planes, around which the turn caused by the angle ∆γ is performed, is set by the angle φ m .By using Equation (5) in the form of r 0 ∆V z V 0 sin( ∆φ 2 ) = z, where ∆V z V 0 ≈ ∆γ, and the found equation for ∆γ, we obtain Equation ( 19) allows us to find the normal acceleration for the orbit plane.

Assessment of Influence of the Transversal and Radial Components of the Thrust Vector on Orbit Elements Deviations
Due to the burn with transversal orientation on the angular interval of the argument of latitude ∆φ, the eccentricity and the semimajor axis will be altered by ∆e and ∆a, correspondingly.They can be found by using equations [8] ∆e = 4 w t w c sin ∆φ 2 and ∆a = 2 w t w c ∆φ.The eccentricity alteration will be fulfilled in the direction described by the angle φ m .The radial component of acceleration will not alter the semimajor axis.The eccentricity will be altered in the direction perpendicular to the direction of alteration of the eccentricity by the transversal component according to the magnitude ∆e = 2 w r w c sin ∆φ 2 .The alteration of the orbit plane orientation is determined by the Formula (18).
By altering the orbit elements in this point by the magnitude of the calculated deviations, we get the orbit with the accounted influence of the non-modeled perturbations.By calculating of the non-modeled accelerations on the intervals between the observations, we can form the model of additional accelerations on the whole revolution and use it for subsequent orbit determinations and motion propagations of these objects.

Results
The following two examples (see Tables 1-4) demonstrate the capabilities of the aforementioned method.The examples with the known realized velocity impulses were taken to assess the accuracy of the found solution.

The First Tracklet
The Second Tracklet for the tolerable calculation of possibility of the collision with these objects and guaranteed collision avoidance.This can be explained by the fact that traditional perturbing models used for the motion propagation of these objects differ considerably from the real ones.Furthermore, the orbits of approximately 300 (10%) objects are being propagated so poorly that these objects are often lost and drop off from the catalogue.This mainly relates to the membranous objects with big varying area-to-mass ratios.This, and previously published works [3], allow us to increase the quality of motion propagation of the challenging objects, which comprise almost 50% of the objects in the GSO, registered in different catalogues.First of all, this refers to the maneuvering objects.In contrast to the traditional works in which the orbit is determined after the maneuver and the maneuver itself is assessed afterwards, or when the assessment of the maneuver occurs simultaneously with the orbit determination, in this work, the maneuvers are assessed by the measurements and the orbit of the maneuvering SO is assessed with their help.This allows us to determine the MSO orbit fast and with high accuracy as only one or two measurements are used.The given examples verify this.The simplicity of the described maneuver assessment method ensures its simple realization and reliability.
We also hope that this work will allow us to approach the solution for the problem of propagation of motion of the membranous objects with big and variable area-to-mass ratios.By substituting the non-modeled perturbation (for example, solar radiation) with the influence of the low-thrust engine work, one can make the model of these perturbations on the intervals between the measurements and account for the influence of these perturbations during the propagation.These models are updated on each interval between the measurements.
Thus, this work, with the work [3], will allow us to improve considerably the orbit determination technology of approximately half of the objects in the geostationary orbit, currently registered in the catalogues of different organizations.
The linearized equations were used while solving this problem.This restricts the use of the suggested method to the assessment of the maneuvers performed in near-circular orbits.Since the optical measurements are used, this method fits the assessment of the maneuvers of spacecraft in the GSO and the medium Earth orbits (the navigation system orbits).
The maneuvers of the maintenance of these orbits are small enough.Hence, the linearization conditions are not violated.The use of the iteration procedure described in [8] allows for the widening the constraints of the usage of this method for the magnitude of the assessed parameters (more than 100 m/s) and for the deviations from circular orbits (several hundreds of kilometers).However, this was not the case for the practical problems of this paper.
transversal t w and radial r w components of the perturbing acceleration.

Table 1 .
Real maneuver parameters for Example 1.

Table 2 .
The deviations for the tracklets for Example 1.

Table 3 .
Real maneuver parameters for Example 2.

Table 4 .
The deviations for the tracklets for Example 2.