Solution for Nonlinear Three-Dimensional Intercept Problem with Minimum Energy

Abstract

Classical orbit intercept applications are commonly formulated and solved as Lambert-type problems, where the time-of-flight (TOF) is prescribed. For general three-dimensional intercept problems, selecting a meaningful TOF is often a difficult and an iterative process. This work overcomes this limitation of classical Lambert’s problem by reformulating the intercept problem in terms of a minimum-energy application, which then generates both the desired initial interceptor velocity and the TOF for the minimum-energy transfer. The optimization problem is formulated by using the classical Lagrangian 𝑓 and 𝑔 coefficients, which map initial position and velocity vectors to future times, and a universal time variable 𝑥. A Newton-Raphson iteration algorithm is introduced for iteratively solving the problem. A generalized problem formulation is introduced for minimizing the TOF as part of the optimization problem. Several examples are presented, and the results are compared with the Hohmann transfer solution approaches. The resulting minimum-energy intercept solution algorithm is expected to be broadly useful as a starting iterative for applications spanning: targeting, rendezvous, interplanetary trajectory design, and so on.

1. Introduction

A fundamental problem of astrodynamics is concerned with computing intercept trajectories or interplanetary mission orbit for objects in space [1, 2]. These calculations are often performed assuming a predetermined time-of-flight (TOF). This is the well-known Lambert’s problem [3–6]. Applications of Lambert’s problem are common in interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on. These problems are formulated by specifying the initial position vectors of an interceptor and target satellite. When the TOF is specified, the initial velocity vector for the interceptor is an unknown implicit function of the local gravity field and can only be recovered by a successive approximation strategy. Other authors have considered alternative problem formulations for a specified TOF that have included minimum-fuel and multiple-impulse strategies [7, 8]. This work overcomes the limitations of these approaches by introducing a minimum-energy approach that simultaneously generates both the TOF and the initial velocity vector for the interceptor.

For the special case of coplanar orbits, the Hohmann transfer algorithm generates a two-impulse minimum-energy orbit transfer by using tangential burns [3–5]. This technique provides a reference orbit transfer for various space applications. For direct applications of the Hohmann transfer to interplanetary orbit transfer, the position vectors of the target planet and initial departure planet are specified assuming a prescribed TOF. When the spacecraft is far from the initial position, one must be alert to the possibility that a multiorbit maneuver may be required.

Clearly, the TOF is a critical parameter for various applications. Once a TOF is determined, the rest of the procedure is solved readily by the solution of Lambert’s problem. This work addresses the problem that there are no adequate methods available for determining a TOF, especially, in general three-dimensional (3D) cases. The problem of finding an optimal TOF only becomes well defined when one specifies a minimization criteria. To this end, a minimum-energy version of classical interceptor problem is formulated for recovering the TOF for a 3D orbit transfer. The results of this calculation are useful as a reference value for interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on. Of course, one can also bound the range of achievable transfer trajectory times by solving for the minimum TOF consistent with the maximum energy that can be generated. Yielding a mission design space that spans the range of TOF consists in the range[TOF_min,TOF_max].

The design goal for the optimization problem is to simultaneously recover the required initial interceptor velocity and the TOF for the intercept. The mathematical advantage of this approach is that the problem has a unique optimal solution, rather than the family of solutions that characterize the classical Lambert’s problem. Mathematically, the problem is defined by a constrained optimization algorithm. Particular care is exercised in formulating the problem for handling the near-parabolic orbits that arise in intercept applications. Analytically, this is handled in a comprehensive way by introducing a universal variable that permits a single TOF equation to be developed that is valid for all conic orbits.

This work is organized in three sections. First, Kepler’s equation is used to define the TOF equation. This is followed by a description of the universal variable used for the problem formulation. For completeness, Lambert’s problem is briefly described. Second, the minimum-energy problem for the intercept problem is introduced and solved. Third, simulation results are presented which compare the TOF obtained for an interplanetary trajectory design with a trajectory developed using the Hohmann transfer methodology and interceptor design solution approaches.

2. Mathematical Review

A fundamental approach for determining the TOF for spacecraft starts with Kepler’s equation that is given by 𝑀=𝑛(𝑡−𝑇)=𝐸−𝑒sin𝐸, (1) where𝑀is the mean anomaly,𝐸and𝑒denote the eccentric anomaly and the eccentricity, respectively,𝑇is the time of periapsis passage,𝑡is the TOF,𝑛is the mean motion defined as𝜇/𝑎³,𝜇denotes the gravitational constant, and𝑎is the semimajor axis of orbit.

As𝑒∼1,the solution for Kepler’s equation becomes more difficult to obtain. This problem is overcome by introducing the universal variable given by [3] ̇√𝑥=𝜇𝑟, (2) where𝑟is the position of spacecraft. As shown in [3], by introducing the universal variable defined by (2), one can express Kepler’s equation and the radial spacecraft coordinate in the following form: √𝜇𝑡=𝑎𝑥−√𝑥𝑎sin√𝑎𝐫+𝑎⁰⋅𝐯₀√𝜇𝑥1−cos√𝑎+𝑟₀√𝑥𝑎sin√𝑎,𝐫𝑟=𝑎+𝑎⁰⋅𝐯₀√𝑥𝜇𝑎sin√𝑎+𝑟₀𝑎𝑥−1cos√𝑎, (3)(4) where𝑇is assumed to be zero without loss of generality and𝐫₀and𝐯₀are the initial position and velocity vectors of spacecraft, respectively. These necessary conditions describe the position and velocity of an orbiting object as a function of time. If the value of the universal variable from (3) is known, the position of the spacecraft at that time is evaluated. Even though (3) is transcendental in𝑥, a Newton’s iteration technique is used to successfully solve for𝑥when the TOF,𝑡, is given.

Assuming that there are no external forces, then the four vectors𝐫₀,𝐯₀,𝐫, and𝐯are assumed to be governed by Keplerian motion. To compute𝐯and𝐫in terms of𝐯₀,𝐫₀, and𝑥, the position and velocity vectors of spacecraft at time𝑡are described as [9] 𝐫=𝑓𝐫₀+𝑔𝐯₀,𝐯=𝑓𝐫̇₀+𝑔𝐯̇₀, (5) where𝑓,𝑔,𝑓̇, and𝑔̇are scalar time-dependent constants, which are subject to the following constraint: 𝑓𝑔̇−𝑓𝑔̇=1, (6) where 𝑎𝑓=1−𝑟₀𝑥1−cos√𝑎,𝑎𝑔=𝑡−√𝜇𝑥−√𝑥𝑎sin√𝑎,𝑓̇=−√𝜇𝑎𝑟𝑟₀𝑥sin√𝑎,̇𝑎𝑔=1−𝑟𝑥1−cos√𝑎. (7)

The energy minimum form of Lambert’s problem is solved by introducing the classical Lagrangian coefficients and universal variable in the problem necessary conditions.

3. Time-of-Flight for Minimum-Energy Orbit Transfer

The major objective in this paper is to compute (i) the TOF and (ii) the initial velocity for an interceptor object for two arbitrary given position vectors so that the transfer orbit energy is a minimum. This problem differs from the classical Lambert’s problem, which fixes a TOF and only recovers a solution for the initial velocity for the interceptor object, given initial and final position vectors of spacecraft.

The optimization problem is formulated by assuming that a target and an interceptor exist in arbitrary orbits, respectively. The problem geometry is illustrated in Figure 1, where𝐫₀and𝐯₀denote the initial position and velocity vectors of the target, respectively, and𝐫̂₀and𝐯₀represent the initial position and velocity vectors of interceptor, respectively.

Figure 1 

Geometry of the minimum-energy problem for a TOF.

The unknowns for the problem are the TOF and initial velocity correction for the interceptor. The goal of the trajectory optimization is to reduce the displacement position vector locating the interceptor relative to the target to zero values at the TOF, while minimizing the orbit energy of the interceptor. The problem is formulated as a nonlinear optimization problem.

4. Constrained Optimization Problem

For given𝐫₀,𝐯₀, and𝐫̂₀, find𝑡and𝐯₀by minimizing the performance index defined as the interceptor’s orbit energy,𝒥, defined as 𝑣̂𝒥=²⁰2−𝜇𝑟₀̂ (8) subject to 𝜼𝑥,𝑥,𝐯₀=𝜂()𝜂𝐯,𝑡𝑥,𝑡𝑥,₀̂,𝑡=0,𝐫−𝐫=0, (9)where𝑥and𝑥(9) denote the universal variables for the target and the interceptor, respectively, and𝑡is the TOF to be determined.

The displacement vectors for the target and interceptor are expressed using𝑓and𝑔as follows: 𝐫−𝐫̂=𝑓𝐫₀+𝑔𝐯₀−𝑓𝐫̂₀+𝑔𝐯₀. (10)

As a constraint vector,𝜼∈ℛ², (3) for𝑥and𝑥is rewritten as 𝜂(𝑥,𝑡)=𝑎𝑥−√𝑥𝑎sin√𝑎+𝑟₀√𝑥𝑎sin√𝑎𝐫+𝑎⁰⋅𝐯₀√𝜇𝑥1−cos√𝑎−√𝜂𝐯𝜇𝑡,𝑥,₀,𝑡=𝑎𝑥−√𝑥𝑎sin√𝑎+𝑟₀̂√𝑥𝑎sin√𝑎+𝑎𝐫̂₀⋅𝐯₀√𝜇𝑥1−cos√𝑎−√𝜇𝑡. (11) Note that the augmented variables to be obtained are𝐯₀,𝑥,𝑥, and𝑡.

5. Optimal Necessary Conditions

Since the second term of the energy is constant, it does not affect the performance index so that the index is redefined, without loss of generality, as [10] 𝒥𝐯₀=12𝐯^𝑇₀𝐯₀. (12)

The Hamiltonian is formed by appending the constraints of (9) with Lagrange multipliers as follows: 𝐻=𝒥𝐯₀+𝝀^𝑇𝜼𝑥,𝑥,𝐯₀,𝑡+𝝓^𝑇(𝐫−𝐫̂), (13) where𝝀=[𝜆₁,𝜆₂]^𝑇and𝝓=[𝜙₁,𝜙₂,𝜙₃]^𝑇. To minimize the performance index with respect to the augmented variables, the necessary conditions provide the following [11]: 𝜕𝐻𝜕𝑥=𝜆₁𝑟+𝝓^𝑇𝜕𝐫𝜕𝑥=0,𝜕𝐻𝜕𝑥=𝜆₂𝑟̂−𝝓^𝑇𝜕𝐫̂𝜕𝑥=0,𝜕𝐻√𝜕𝑡=−𝜇𝜆₁−√𝜇𝜆₂+𝝓^𝑇𝐯₀−𝐯₀=0,𝜕𝐻𝜕𝐯₀=𝐯^𝑇₀+𝜆₂𝜕𝜂𝜕𝐯₀−𝝓^𝑇𝜕𝐫̂𝜕𝐯₀=0, (14)(15)(16)(17) where𝜕𝜂/𝜕𝐯₀∈ℛ^1×3represents a row vector,𝜕𝐫̂/𝜕𝐯₀∈ℛ^3×3is a matrix (refer to the Appendix for detail derivation), and 𝜕𝜂𝜕𝜂𝜕𝑥=𝑟,𝜕𝑥=𝑟.̂ (18)

The necessary conditions of (14)–(17) are simplified by the following manipulations. First, from (15), the Lagrange multiplier𝜆₂is obtained as 𝜆₂=1𝑟̂𝝓^𝑇𝜕𝐫̂𝜕𝑥. (19)

Second, substituting (19) into (17) yields𝐯^𝑇₀+𝝓^𝑇1𝑟̂𝜕𝐫̂𝜕𝑥𝜕𝜂𝜕𝐯₀−𝜕𝐫̂𝜕𝐯₀=0 (20) which can be solved for𝝓, leading to 𝝓^𝑇=𝐯^𝑇₀𝜕𝐫̂𝜕𝐯₀−1𝑟̂𝜕𝐫̂𝜕𝑥𝜕𝜂𝜕𝐯₀⁻¹. (21)

Third, by using (14), the Lagrange multiplier𝜆₁follows as 𝜆₁1=−𝑟𝝓^𝑇𝜕𝐫.𝜕𝑥 (22)

Collecting the Lagrange multiplier solutions from (19) and (22), introducing the results into (16), one obtains √𝜇𝑟𝝓^𝑇𝜕𝐫−√𝜕𝑥𝜇𝑟̂𝝓^𝑇𝜕𝐫̂𝜕𝑥+𝝓^𝑇𝐯₀−𝐯₀=0. (23)

This equation is further simplified by recalling the terminal constraint𝑟=𝑟̂, leading to 𝝓^𝑇𝜕𝐫−𝜕𝐫̂𝜕𝑥𝜕𝑥+𝑟√𝜇𝐯₀−𝐯₀=0. (24)

Substituting (21) into (3) yields the final necessary condition required for finding the TOF for the intercept problem: 𝐯^𝑇₀𝐿𝜕𝐫−𝜕𝐫̂𝜕𝑥𝜕𝑥+𝑟√𝜇𝐯₀−𝐯₀=0, (25) where the new matrix is defined for simplicity as 𝜕𝐫̂𝐿=𝜕𝐯₀−1𝑟̂𝜕𝐫̂𝜕𝑥𝜕𝜂𝜕𝐯₀⁻¹. (26)

Satisfaction of the new equation implies that the interceptor can meet the target with minimum energy in a time provided by the computed TOF, not by a predetermined time.

6. Summary

The approach for obtaining the nonlinear 3D intercept problem by using the classical Lagrangian𝑓and𝑔coefficients is summarized. Computing the TOF and the interceptor’s initial velocity is the goal to meet the requirement that the final distance between the two spacecrafts becomes zero. Firstly, one can simply obtain𝑎and𝑟₀using the following: 𝑟₀=𝐫₀,𝜇𝑎=−,2ℰ (27) where the orbital energy is given by 𝑣ℰ=²⁰2−𝜇𝑟₀. (28)

Then, one can find the universal variables, the initial velocity for the interceptor, and the TOF using the following equations: 𝜼𝑥,𝑥,𝐯₀,𝑡=0,𝑓𝐫₀+𝑔𝐯₀−𝑓𝐫̂₀−𝑔𝐯₀𝐯=0,_𝑇0𝐿𝜕𝐫−𝜕𝐫̂𝜕𝑥𝜕𝑥+𝑟√𝜇𝐯₀−𝐯₀=0, (29) where the semimajor axis𝑎of the interceptor can be iteratively computed with estimated𝐯₀. The Newton-Raphson iteration algorithm is applied to solve the previous equations. Next, one can compute all of the𝑓and𝑔expressions using (7). Then, (5) is applied to obtain the final position and velocity vectors.

There are many feasible performance indices to specify a TOF. For example, consider the candidate performance index 𝒥𝐯₀=1,𝑡2𝐯^𝑇₀𝐯₀+𝛼𝑡, (30) where𝛼is nonnegative weight. By adding the time as one part of the performance index, the TOF is expected to be shortened with respect to the variation of𝛼. In a similar manner with the minimum-energy procedure in the previous section, the optimization solution to this problem is readily determined. The partial derivative of the Hamiltonian𝐻with respect to𝑡is given by 𝜕𝐻√𝜕𝑡=−𝜇𝜆₁−√𝜇𝜆₂+𝝓^𝑇𝐯₀−𝐯₀+𝛼=0. (31)

Finally, a cost-effective equation weighted to the time is obtained as 𝐯^𝑇₀𝐿𝜕𝐫−𝜕𝐫̂𝜕𝑥𝜕𝑥+𝑟√𝜇𝐯₀−𝐯₀+𝑟√𝜇𝛼=0. (32)

Numerical convergences based on different methods and their overall computational cost depend on the chosen parameterization, the initial guess, and the numerical technique used for solving the resulting equation. It is known that singularities exist when solving Lambert’s problem that prevent some algorithms from converging for particular cases or make convergence extremely slow. For example, Lambert’s method fails when the transfer angle is 180 degrees [9]. Therefore, the features of the suggested method must be analyzed. However, this is out of the scope of this paper, which is focused on approaches to determine the TOF and initial velocity of the interceptor.

7. Application Examples

The specification of a TOF for an intercept problem is generally not unique, and a family of solutions are possible when the initial trust level is variable. As a result, the process of determining a useful TOF requires experimentation and iteration. The minimum-energy optimization approach of this work finds a unique value for the TOF. The solution for the intercept problem simultaneously determines the initial interceptor velocity vector and TOF. Numerical examples are presented that compare and contrast the classical Hohmann transfer with the proposed method.

Let us briefly review the Hohmann transfer and compare the minimum-energy problem with it. The geometry of the Hohmann transfer is illustrated in Figure 2. The distances of the departure and arrival orbits are denoted as𝑟₁and𝑟₂, respectively.

Figure 2 

Geometry of the Hohmann transfer.

The semimajor axis for the elliptic orbit and the energy are given by [3] 2𝑎_ℎ=𝑟₁+𝑟₂,ℰ_ℎ𝜇=−𝑟₁+𝑟₂, (33)(34) and, then, the departure velocity of the transfer orbit is readily obtained as 𝑣_ℎ0=2𝜇𝑟₁+ℰ_ℎ. (35)

Since the velocity of the departure orbit is given by 𝑣₁=𝜇𝑟₁, (36) the velocity change for the Hohmann transfer is calculated as Δ𝑣=𝑣_ℎ0−𝑣₁, (37) and the TOF of the Hohmann transfer is written as 𝑡=𝜋𝑎³_ℎ𝜇, (38) where𝑎_ℎis obtained from (33). Two circular orbits are assumed with the radii of𝑟₁=4000km and𝑟₂=6000km, respectively. Then, the velocity of the departure circular orbit is𝑣₁=9.9825km/s, and the remaining parameters for the Hohmann transfer are obtained as𝑡=1759.3sec,𝑣_ℎ0=10.9353km/s, andΔ𝑣=0.9528km/s. To navigate to the final position of the arrival orbit by the Hohmann transfer, the initial position and velocity vectors are assumed to be given by 𝐫₁=[0,−4000]^𝑇,𝐫₂=[4097.2993,−4383.1653]^𝑇,𝐯₁=[9.9825,0]^𝑇,𝐯₂=[5.9543,5.5660]^𝑇. (39)

Four cases including the Hohmann case are analyzed. Initial velocities, velocity changes, and TOF obtained by the solution of the proposed minimum-energy problem are arranged in Table 1. Also, the initial positions of target spacecraft and their resultant transfer trajectories are displayed in Figure 3.

Table 1 

Numerically computed transfer orbits.

Figure 3 

Trajectory generation by initial positions of target.

As shown in Table 1, the result for the Hohmann case is nearly identical to the output from the classical approach in (35)–(38) with a small numerical error. It proves that the proposed approach provides optimal solutions we are looking for. Moreover, it is obvious that it gives the minimum velocity change, which is tangential with the trajectory, compared with the other cases. If the target is positioned at case 1, 2, or 3, relative to the interceptor’s initial position and it is required to start the orbit transfer mission at this time, it would be a great advantage to have a reference minimum-energy trajectory to accomplish the mission. Fortunately, the results in Table 1 can be utilized, since they represent the minimum velocity in each case. This means that there are no more efficient trajectories in these cases than the transfer orbit listed in Table 1. When the target is positioned forward compared to the Hohmann transfer, the phase angle, sometimes called flight-path-angle, at departure should be negative to meet the optimal trajectory requirement. When the target, on the contrary, is positioned backward, the flight-path-angle should be positive.

Even if a circular orbit is selected for the comprehensive analysis by comparing with the Hohmann transfer, the application of the proposed approach is not limited. Therefore, an illustrative example in Figure 4 is conducted to demonstrate the performance of space maneuver of the interceptor. There are two arbitrary elliptic orbits, which are not coplanar. The initial positions of the target and interceptor orbit are depicted in Figure 4. By solving the nonlinear 3D intercept problem through minimizing energy, the TOF and initial velocity is computed. Then, applying minimum-energy velocity obtained for the interceptor confirms that the final distance between the two orbits is zero at the computed TOF. Note that the problems formulated by the universal variable and𝑓and𝑔expressions in this paper are solved in 3D space for supporting the design of arbitrary intercept problems with minimum energy.

Figure 4 

An illustrative example for the intercept problem.

The intercept problem is easily generalized by introducing a time weighting factor in the definition of the optimal control problem, which allows a systematic exploration of the intercept design space as a function of the allowed transfer energy. In this example, the outer orbit is selected as a target orbit to be captured as illustrated in Figure 5. The initial position of the interceptor is at the inner orbit. The TOF obtained from the optimization problem is shortened when𝛼increases, and the results are illustrated in Figure 5. As expected, the longest TOF is obtained when𝛼is zero. If𝛼approaches one in this simulation case, the interceptor can hit the target in a very short time. It means that shortening the TOF is the optimal solution to minimize the chosen cost function.

Figure 5 

TOF for intercept due to variation of𝛼.

8. Conclusions

A general algorithm is presented for generalizing the classical Lambert’s transfer problem, where the determination of time-of-flight (TOF) for a spacecraft intercept, in arbitrary three-dimensional orbit, is addressed. A constrained optimization technique is introduced and iteratively solved to find both the TOF and the initial intercept velocity vector. The proposed algorithm provides a benchmark minimum-energy solution that provides an optimal reference trajectory. A significant advantage of this approach is that the TOF is unique when the energy is minimized. This implies that the interceptor with lower energy than the evaluated minimum energy cannot meet the target. Numerical results are presented, and they compare the intercept solutions with those obtained using the classical Hohmann transfer technique. The proposed algorithm is expected to be broadly useful for all classes of intercept problem that have a Lambert-like character.

Appendix

The Partial Derivatives

The partial derivatives𝐫and𝐫̂with respect to𝑥and𝑥, respectively, are given by 𝜕𝐫=𝜕𝑥𝜕𝑓𝐫𝜕𝑥⁰+𝜕𝑔𝐯𝜕𝑥⁰,𝜕𝐫̂𝜕𝑥=𝜕𝑓𝜕𝑥𝐫̂₀+𝜕𝑔𝜕𝑥𝐯₀. (A.1)

Applying (4) and using (7), the partial derivatives with respect to the universal variable are written as 𝜕𝑓√𝜕𝑥=−𝑎𝑟₀𝑥sin√𝑎,𝜕𝑔𝑎𝜕𝑥=−√𝜇𝑥1−cos√𝑎,𝜕𝑓𝜕𝑥√=−𝑎𝑟₀̂𝑥sin√𝑎,𝜕𝑔𝜕𝑥𝑎=−√𝜇𝑥1−cos√𝑎. (A.2)

The orbit energy has the following relationship: 𝐯^𝑇₀𝐯₀2−𝜇𝑟₀̂𝜇=−2𝑎. (A.3)

Since𝑟₀̂is a constant in this case, differentiating both sides with respect to𝐯₀yields 𝐯^𝑇₀=𝜇2𝑎²𝜕𝑎𝜕𝐯₀. (A.4)

The partial derivative of𝑎with respect to𝐯₀is readily written as 𝜕𝑎𝜕𝐯₀=2𝑎²𝜇𝐯^𝑇₀. (A.5)

The row vector, partial derivative of𝜂with respect to𝐯₀, is given by 𝜕𝜂𝜕𝐯₀=𝜕𝜂(𝑎)𝜕𝑎𝜕𝑎𝜕𝐯₀+𝜕𝜂(𝐯₀)𝜕𝐯₀𝑎, (A.6) where 𝜕𝜂(𝑎)𝜕𝑎=𝑥−√𝑥𝑎sin√𝑎+𝑎2𝑥𝑎𝑥cos√𝑎−1√𝑎𝑥sin√𝑎+𝐫̂^𝑇₀𝐯₀√𝜇𝑥1−cos√𝑎−𝑥2√𝑎𝑥sin√𝑎+𝑟₀̂21√𝑎𝑥sin√𝑎−𝑥𝑎𝑥cos√𝑎,𝜕𝜂(𝐯₀)𝜕𝐯₀𝑎=𝑎√𝜇𝑥1−cos√𝑎𝐫̂^𝑇₀. (A.7)

Next, the partial derivative of𝐫̂with respect to𝐯₀follows as 𝜕𝐫̂𝜕𝐯₀=𝜕𝑓𝜕𝑎𝐫̂₀+𝜕𝑔𝜕𝑎𝐯₀𝜕𝑎𝜕𝐯₀+𝜕𝐫̂𝜕𝐯₀𝑎, (A.8) where 𝜕𝑓𝜕𝑎1=−𝑟₀̂𝑥1−cos√𝑎−𝑥2√𝑎𝑥sin√𝑎,𝜕𝑔𝜕𝑎1=−√𝜇3√𝑥−𝑎2𝑥sin√𝑎+𝑥2𝑥cos√𝑎,𝜕𝐫̂𝜕𝐯₀𝑎=𝑔𝐼_3×3. (A.9)

The matrix,𝐿, in (26), consists of the combination of two vectors and one scaled identity matrix. If𝑔is not zero, the matrix would have a full rank of three. The position and velocity vectors in orbit are in general not parallel. Since the condition that𝑔is zero means that𝐫and𝐯are parallel, this is impossible in orbit. It could guarantee the existence of the inverse of the matrix.

Acknowledgment

This study was supported by research fund from Chosun University, 2013.