Next Article in Journal
Investigation of the Impact of a Particle Foam Insulation on Airflow, Temperature Distribution, Pressure Profile and Frost Buildup on the Aircraft Structure
Previous Article in Journal
A Text-Driven Aircraft Fault Diagnosis Model Based on Word2vec and Stacking Ensemble Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Guidance Law with Impact-Angle Constraint and Acceleration Limit for Exo-Atmospheric Interception

School of Astronautics, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Aerospace 2021, 8(12), 358; https://doi.org/10.3390/aerospace8120358
Submission received: 6 October 2021 / Revised: 18 November 2021 / Accepted: 19 November 2021 / Published: 23 November 2021

Abstract

:
This paper aims to develop an optimal guidance law for exo-atmospheric interception, in which impact-angle constraint and acceleration limit are considered. Firstly, an optimal control problem with constraints on terminal miss and impact-angle is formulated, in which the control energy performance index is weighted by a power function of the time-to-go. The closed-loop guidance command, which is expressed as a linear combination of zero-effort miss distance and the zero-effort angle error, is derived using a traditional order reduction transformation. Then, an analytical solution to the maximal acceleration during the flight is obtained by analyzing the boundary points and critical points of the guidance command curve. It is found that the maximal acceleration is a function of the weighted gain in the performance index. Therefore, the maximal acceleration can be efficiently limited by using the variable weighted gain. Furthermore, the relationship between the total control energy and the weighted gain is studied. As a result, a systematic method is proposed for selecting the weighted gain so as to meet the constraint of the acceleration while the total control energy is minimal. Nonlinear simulations have been carried out to test the performance of the proposed method. The results show that this method performs well in intercepting the maneuvering target with a negligible miss distance and intercept angle error. And it can tolerate a stricter acceleration limit in comparison with the typical method.

1. Introduction

The terminal guidance law, which has direct influence on kill efficiency, is immensely important for the weapon system. During terminal attack, the achievement of a proper impact angle is an important goal of homing systems. Using impact-angle-control, the missile is capable of attacking a weak spot of the target to maximize warhead effectiveness and ensure a high kill probability [1,2]. In addition, the acceleration limit is another crucial constraint on homing missiles, which will lead to large miss-distance and impact angle error if the acceleration saturation occurs near the target.
Guidance law design with specific impact direction has been a hot topic during the past four decades. As an application of optimal control theory, Bryson and Ho [3] have proposed the optimal solution to the simple rendezvous problem with terminal constraints on the lateral position and velocity. In [4], a similar formulation was used to derive a linear quadratic optimal guidance law and linear quadratic differential game-based laws for maneuvering target scenarios. These two formulations can be utilized to impose a terminal intercept angle by selecting the ratio between the relative terminal velocity and the closing speed. Kim and Grider [5] first proposed an optimal guidance law to ensure the impact angle for a reentry vehicle, which pursues a fixed or slowly moving ground target. In [6], York and Pastrick further improved the guidance law presented in [5], in which the lag in the autopilot is considered. The resulting guidance law was applied to an air defense interceptor in [7]. In [8], a generalized form of energy minimization optimal guidance law was proposed for satisfying the impact-angle constraints. Song et al. [9] have studied an optimal impact angle guidance law for varying velocity missiles against a maneuvering ship. Shaferman [10] derived two impact angle control guidance laws based on the optimal control and differential game theory for an evasive maneuvering target. An optimal control planar interception law against maneuvering targets with known trajectories was devised in [11]. A terminal body angle control guidance law considering the angle-of-attack dynamics was developed in [12]. It is noted that the abovementioned optimal guidance laws are all based on a constant weighted performance index, which results in the weak ability of the trajectory shaping. And the terminal acceleration cannot be guaranteed to be zero, which possibly increases the miss-distance.
In order to address this drawback, Ryoo [13] proposed a time-to-go weighted optimal guidance with zero terminal acceleration because the weighting of cost infinitely increases as the time-to-go approaches zero. Ohlmeyer [14,15] proposed an optimal guidance law with a weighted gain specified by the user, called Generalized Vector Explicit Guidance, which improves the ability to shape the trajectory and achieves the impact-angle constraint. He [16] investigated a similar guidance law which considers the influence of the gravity and thrust in the exo-atmospheric interception engagement. Lee [17] further developed a generalized optimal guidance law with arbitrary weighted function to satisfy the impact-angle constraint. By considering the penalty term regarding the miss distance and intercept angle weighted by hyperbolic tangent functions in the performance index, Xiong [18] proposed a new optimal guidance law, of which the guidance performance was less affected by the time-varying speed. Li [19] further designed a unified optimal impact angle guidance law with the energy consumption weighted by an arbitrary positive function, which can adapt to different types of interceptors and targets. In addition to the guidance law based on the optimal control theory, biased proportional navigation [20,21,22,23,24] and sliding-mode guidance law [25,26,27,28,29] can also realize the impact with a specified angle. However, all these guidance laws do not consider the acceleration limit, which possibly leads to large miss-distance and impact angle error if the acceleration saturation occurs near the target.
In order to limit acceleration, Rusnak [30] and Hexner [31] directly solved the optimal control problem with control constraint, which leads to a significantly complicated solution. Weiss [32] introduced a new class of optimal guidance law in which a special performance index penalizing the variability of the guidance command is used, and therefore, an additional design parameter can be profitably used to avoid acceleration saturation. Although the acceleration limit is considered in this research, the impact-angle cannot be guaranteed. This is due to the fact that if the constraint of impact-angle is considered simultaneously, the solved guidance command is extremely complex and even cannot be solved analytically.
In this paper, an optimal guidance law with constraint on the impact-angle and acceleration limit is proposed. Firstly, an impact-angle constrained optimal guidance law (IACOGL) with a variable weighted gain is formulated. The closed-loop guidance command can be expressed as a linear combination of zero-effort miss distance (ZEM) and the zero-effort angle error (ZEAE). Then, we analyze the extreme value of the guidance command curve and derive the closed-form function between the maximum of guidance command with the weighted gain. Therefore, the maximal acceleration can be efficiently limited by using the variable weighted gain. Furthermore, the relationship between the total control energy and the weighted gain is studied. Finally, a systematic method to select the optimal weighted gain, which minimizes the total control energy and keeps the acceleration within the constraint, is proposed. The greatest difference from the previous works is that the proposed method is able to avoid acceleration saturation by choosing the appropriate weighted gain in an analytical manner, while having optimal performance for the guidance command. Several numerical simulations are carried out to evaluate the performance of the proposed method. The results show that this method performs well in intercepting the maneuvering target with the impact angle constraint and acceleration limit.
The remainder of this paper is organized as follows: In Section 2, the mathematical models used for the guidance law derivation and simulation are presented. Section 3 derives the time-to-go weighted optimal guidance law with constraint on impact-angle. Section 4 provides the derivation process on the maximum of guidance command and the total maneuvering energy. After which, a method used to select the weighted gain to avoid the acceleration saturation is also presented in Section 4. Finally, nonlinear simulations are carried out to evaluate the performance of the proposed law.

2. Model Formulation

Consider a planar engagement geometry for a maneuvering target shown in Figure 1. Here, X I O I Y I is a Cartesian inertial reference frame. We denote the missile and target by the subscripts M and T, respectively. The velocity, acceleration, and flight-path angles are denoted by V , a and γ , respectively. It is assumed that the velocity of missile and target for guidance law derivation is constant. The range between the adversaries is r , and λ is the angle between the line of sight (LOS) and the X I axis. The X R axis, aligned with the initial LOS, is denoted for linearization. The relative displacement between the target and the missile normal to this axis is ξ . The intercept angle between missile and angle is defined as γ I = γ T + γ M , which is required to be a given value γ f c at the intercept time.
From Figure 1, the equations of relative motion between missile and target are given by
r ˙ = V T cos γ T + λ + V M cos γ M λ
λ ˙ = V M sin γ M λ + V T sin γ T + λ / r
And, due to the existence of acceleration and gravitational effect, the flight-path angle and velocity of the missile as well as the target evolve according to
γ ˙ M = a M cos γ M λ g M cos γ M / V M
V ˙ M = a M sin γ M λ g M sin γ M
where, g M stands for the gravitational acceleration of the missile. Similarly, due to the gravitational effect, the flight-path angle and velocity of target evolve according to
γ ˙ T = g T cos γ T / V T
V ˙ T = g T sin γ T
where, g T is the gravitational acceleration of the target. g M and g T are assumed as a constant g because the duration of the terminal guidance phase is typically very short.
It is assumed that the missile’s closed-loop dynamics can be represented by the following first-order system:
a ˙ M = u a M / T
where, u is the guidance command, T is the time constant of the missile’s dynamics.
As was done in [10], the derivation of the guidance law in this paper will be performed based on a linearized model. In the linear model, the kinematics equations of ξ and γ I are
ξ ¨ = a M
γ ˙ I = k M a M + Δ g γ / V M
where, k M = cos γ M λ and Δ g γ = g cos γ M + cos γ T V M / V T . In the derivation of the guidance law, V M , k M , and Δ g γ are regarded as the constants and will be updated in each guidance step. The state vector of the linearized problem is defined as
x = ξ ξ ˙ γ I T
The matrix form of the kinematic equations is
x ˙ = A x + B u + C
where
A = 0 1 0 0 0 0 0 0 0 , B = 0 1 k M / V M ,   C = 0 0 Δ g γ / V M

3. Derivation of IACOGL

In this section, an impact-angle constraint guidance law with a variable weighted gain is derived on the foundation of optimal control theory. Now, let us consider the following optimal control problem: Find u t that minimizes defined by
J = 1 2 t 0 t f u 2 / t g o N d t
subject to (11) and terminal constraints given by
ξ t f = 0 , γ I t f γ f c = 0
where, t g o = t f t is the time-to-go; N is weighted gain, which is the shaping parameter of the guidance command.

3.1. Order Reduction Transformation

Bryson and Ho [3] introduced a transformation enabling reducing an optimal control problem’s order, which is denoted as terminal projection and can be used to simplify the derivation. This transformation can be easily applied to the resulting optimal control problem in previous section.
As was done in [10], let us define a new state vector Z t that satisfies
Z t = D Φ t f , t x t + D t t f Φ t f , τ C d τ D x f c
where Φ t f , t is the transition matrix associated with Equation (11), x f c indicates the terminal constraint of x t , and D is a constant matrix
D = 1 0 0 0 1 × n 0 0 1 0 1 × n
Since
Φ ˙ t f , t = Φ t f , t A
the time derivative of the new state vector Z t is
Z ˙ = D Φ ˙ t f , t x + Φ t f , t x ˙ D Φ t f , t C = D Φ t f , t B u = φ ξ t g o φ γ t g o u
which is state independent and only the function of the guidance command u . The terminal value of Z t can be expressed using (15) as
Z t f = D x t f x f c = ξ t f , γ I t f γ f c T
It is obvious that Z t f should be zeros. Then, using these new variables, the augmented cost function can be written as
J ¯ = υ 1 Z 1 t f + υ 2 Z 2 t f + 1 2 t 0 t f u 2 / t g o N d t
where, υ 1 and υ 2 are Lagrange multipliers; Z 1 t f and Z 2 t f are the first and second variables of Z t f .
Note that, besides reducing the order of the problem, the two variables of the new state vector Z t have an important physical meaning. Z 1 t is zero-effort miss distance (ZEM), and Z 2 t is zero-effort angle error (ZEAE), respectively referring to the miss distance and the impact-angle error if, from the current time onward, the missile will not apply any control.

3.2. Closed-Loop Optimal Guidance Command

According to the theory of optimal control, the Hamiltonian of the problem is
H = 1 2 u 2 / t g o N + λ 1 φ ξ t g o u + λ 2 φ γ t g o u
where λ 1 and λ 2 are co-state variables, respectively. Due to the order reduction transformation, the time derivative of the zero-effort variables is state independent, and therefore the adjoint equations can be presented in a simplified manner.
λ ˙ 1 = H / Z 1 = 0 λ ˙ 2 = H / Z 2 = 0
Then we can obtain the co-state variables:
λ 1 t = λ 1 t f = υ 1 λ 2 t = λ 2 t f = υ 2
For obtaining an analytic solution for the guidance command, it is assumed that the dynamics of missile are lag-free. Therefore, the time derivatives of the Z 1 and Z 2 are simplified as
Z ˙ 1 = φ ξ t g o u = t g o u Z ˙ 2 = φ γ t g o u = u / V M
where V M = V M / k M . Therefore, the optimal guidance command for the missiles, which satisfy
u = arg u min H
are
H / u = 0 u t = λ 1 t g o N + 1 λ 2 t g o N / V M
Substituting (23) into (26) and integrating from t to t f yields the following two coupled algebraic equations:
Z 1 t f = Z 1 t υ 1 t g o N + 3 N + 3 + υ 2 t g o N + 2 N + 2 V M = 0 Z 2 t f = Z 2 t + υ 1 t g o N + 2 N + 2 V M υ 2 t g o N + 1 N + 1 V M 2 = 0
Solving υ 1 , υ 2 and substituting the solution into (26), the optimal guidance command [10] is
u t = N Z E M t g o 2 Z 1 t + N Z E A E V M t g o Z 2 t
where, N Z E M and N Z E A E are the navigation gains, which are expressed as
N Z E M = N + 2 N + 3 N Z E A E = N + 1 N + 2
And the new states Z 1 and Z 2 can be expressed using state x as
Z 1 t = ξ + ξ ˙ t g o Z 2 t = Δ g γ t g o / V M + γ I γ f c
For the guidance law implementation, the time-to-go must be estimated. It can be approximated simply by
t g o = r / r ˙
If the more accurate t g o estimation is needed, the method presented in [8] can be adopted.
It is found that the command of the proposed guidance law is expressed as a linear combination of zero-effort miss distance and the zero-effort angle error, which is similar to that of the OGL-CTIA proposed in [10]. However, the navigation gains are absolutely different from that of the OGL-CTIA. They possess an additional parameter to shape the intercepting trajectory, which can be used to limit the maximal acceleration. Next, we will focus on how to select the most appropriate weighted gain to achieve missile acceleration within the limit.

4. IACOGL with Acceleration Limit

In the previous section, a universal impact-angle constrained optimal guidance law with great trajectory shaping capability is obtained, but acceleration limit is not considered. In this section, the maximum acceleration is limited by adjusting the shaping parameter, which is the weighted gain, N , to avoid acceleration saturation. Firstly, the closed-form solution to the guidance command is proposed, which is used to derive the function of maximum of guidance command and total maneuvering energy. Then a systematic method is proposed for selecting the weighted gain meeting the constraint of the acceleration while the total maneuvering energy is minimal.

4.1. Analytical Solution to Guidance Command

In order to derive an analytical solution to the guidance command, let us focus on the linear dynamics equation presented in (18). Substituting (28) into (18) results in the following optimal states dynamics:
Z ˙ 1 = N Z E M Z 1 / t g o N Z E A E V M Z 2 Z ˙ 2 = N Z E M Z 1 / V M t g o 2 + N Z E A E Z 2 / t g o
It is obvious that the resulting equation set is a first-order linear time-invariant ordinary differential equation (ODE). Its solution can be analytically formulated as
Z 1 t = N + 2 C 1 t g o N + 3 N + 1 C 2 t g o N + 2 Z 2 t = N + 3 C 1 t g o N + 2 N + 2 C 2 t g o N + 1 / V M
where
C 1 = N + 2 Z 10 / t f N + 3 + N + 1 V M Z 20 / t f N + 2 C 2 = N + 3 Z 10 / t f N + 2 + N + 2 V M Z 20 / t f N + 1
and Z 10 , Z 20 are initial values, respectively.
For simplification, a non-dimensional parameter, which is a function of Z 1 , Z 2 , V M , and t g o is introduced.
κ = Z 1 / V M Z 2 t g o
Therefore, Equation (33) can be rewritten as
Z 1 t = V M Z 20 t f N + 2 C ¯ 1 t ¯ g o N + 1 C ¯ 2 t ¯ g o N + 2 Z 2 t = Z 20 N + 3 C ¯ 1 t ¯ g o N + 2 C ¯ 2 t ¯ g o N + 1
where, C ¯ 1 and C ¯ 2 are non-dimensional coefficients which expressed as
C ¯ 1 = N + 2 κ 0 + N + 1 C ¯ 2 = N + 3 κ 0 + N + 2
where, κ 0 = Z 10 / V M Z 20 t f is the initial value of κ . Substituting (36) into (28), the expression to guidance command is derived as
u t = V M Z 20 t f η N C ¯ 1 N + 1 t ¯ g o N + 1 C ¯ 2 N + 3 t ¯ g o N
where, η N = N + 1 N + 2 N + 3 and t ¯ g o = t g o / t f is the dimensionless time-to-go. Therefore, the expression to guidance command can be derived directly if the initial conditions and terminal constraints are provided.

4.2. Maximum of Guidance Command

From a practical point of view, the available acceleration that the missile can provide is limited, so it is necessary to analyze the maximum acceleration during the homing flight. According to the analytical solution to guidance command presented in (38), the maximum of guidance command for the given initial condition and the selected weighted gain N can be derived. Therefore, the designer can judge whether command saturations occur when the acceleration limit is imposed. In this subsection, the maximum of guidance command for different κ 0 is derived in an analytical manner, which is a function of the weighted gain N . Furthermore, the weighted gain which minimizes the maximum of the guidance command is derived, which can be used to limit the maximum of the guidance command during the engagement.

4.2.1. Expression to Maximum of Guidance Command

For simplification, the dimensionless guidance command is defined as
u ¯ t ¯ u t V M Z 20 / t f = η N C ¯ 1 N + 1 t ¯ g o N + 1 C ¯ 2 N + 3 t ¯ g o N
For a continuous function on a closed and bounded interval, a maximum value can be observed at a critical point or boundary point by the extreme value theorem. The critical point of a function is where its time-derivative is zero. It is obvious that the dimensionless guidance command profile is continuously differentiable in the bounded time interval [0,1] and its derivative to dimensionless time is
d u ¯ / d t ¯ = η N C ¯ 1 t ¯ g o N C ¯ 2 / N + 3 t ¯ g o N 1
Then, considering the critical point condition d u ¯ / d t ¯ = 0 , two critical points can be obtained:
t ¯ g o 1 = 0 , t ¯ g o 1 = N C ¯ 2 / N + 3 C ¯ 1
It is obvious that the first solution is one of boundary points, which should be removed. For the second solution, if it is within 0 , 1 , it can be regarded as the critical point. Now, let us define a set N e x t of weighted gain, which makes the second solution within 0 , 1 . It is noted that the set N e x t is dependent on κ 0 . If κ 0 = 1 , 0 < t ¯ g o 1 < 1 degenerates to 0 < N / N + 3 < 1 . Thus, it is apparent that N e x t = 0 , + . And if κ 0 1 , 0 < t ¯ g o 1 < 1 is equivalent to
N + 3 N 2 N 2 / 3 > N N N 1 N N 1 > 0   or   N + 3 N 2 N 2 / 3 < N N N 1 N N 1 < 0
where
N 1 = 3 κ 0 + 2 / κ 0 + 1 ,   N 2 = 3 κ 0 + 3 / 2 / κ 0 + 1
Then we can obtain
N e x t = 0 , min 0 , N 1 , N 2 max 0 , N 1 , N 2 , +
Now, let us consider the special case of N = 0 . In this case, it can be obtained from (39) that the terminal dimensionless guidance command is not zero and its value is
u ¯ f N = 0 = 2 3 κ 0 + 2
when N N e x t , Substituting (41) into (39), the dimensionless guidance command at critical point is
u ¯ 1 N = N + 2 C ¯ 2 t ¯ g o 1 N
It is noted that, if N 0 , u ¯ 1 N is approach to u ¯ f N = 0 . Thus, if N = 0 , the critical point is at the terminal point. Therefore, the set N e x t should include the case, N = 0 . As a result, N e x t for different κ 0 can be expressed as
N e x t = 0 , min 0 , N 1 , N 2 max 0 , N 1 , N 2 , +                         = 0 , N 1 N 2 , + , 1 < κ 0 < 2 / 3 N 2 , + , 2 / 3 κ 0 < 1 / 2 0 , + , κ 0 1 , κ 0 1 / 2
Therefore, the maximum of dimensionless guidance command, u ¯ M a x , can be obtained by the extreme value theorem as
u ¯ M a x N = max u ¯ 0 , u ¯ 1 , N N e x t u ¯ 0 , N N e x t
where u ¯ 0 is the value of u ¯ t ¯ at the boundary point, 0, which is the function of N :
u ¯ 0 N = N + 2 C ¯ 2 1
All in all, the maximum value can be expressed by a uniform formula which is function of N . In order to analyze the specific expression of u ¯ M a x , a function of N is defined to denote the difference between the magnitudes of u ¯ 0 and u ¯ 1 .
Ξ N = u ¯ 0 u ¯ 1 , N N e x t
Obviously, t ¯ g o 1 0 for arbitrary κ 0 and N N e x t . Then, substitute (46) and (48) to (50), Ξ N can be expressed as
Ξ N = N + 2 C ¯ 2 1 C ¯ 2 t ¯ g o 1 N                             = N + 2 f a d d N , C ¯ 2 0 N + 2 f s u b N , C ¯ 2 0 , 1 N + 2 f a d d N , C ¯ 2 1
where, the functions f a d d N and f s u b N are defined as
f a d d N = C ¯ 2 1 C ¯ 2 t ¯ g o 1 N , N N add
f s u b N = C ¯ 2 1 + C ¯ 2 t ¯ g o 1 N , N N sub
where the domains of N add and N sub are
N add = N N e x t | C 2 , 0 1 , ,   N sub = N N e x t | C 2 0 , 1 .
It is obvious that the magnitude relationship between C ¯ 2 and 0,1 is different for various κ 0 . And therefore, the sets N add and N sub have different interval for various κ 0 , which is present in Figure 2. In the figure, N 3 is the special weighted gain which makes C ¯ 2 = 1 . It is
N 3 = 3 κ 0 + 1 / κ 0 + 1
It can be seen that the sign of Ξ N is depended on f a d d N and f s u b N . The Propositions A1 and A2 are presented in Appendix A to analyze the sign of f a d d N and f s u b N . Actually, the sign of Ξ N can be determined according to the magnitude of κ 0 . And then the expression of u ¯ M a x is derived. u ¯ M a x is a piecewise function of N , the piecewise point of which is dependent on κ 0 . Based on the set N e x t , the specific expression of u ¯ M a x is divided into the following four cases.
(1) If κ 0 1 , it can be obtained from the Proposition A1 that C ¯ 2 < 0 and f a d d N < 0 for any N N e x t . And therefore, it is obvious that Ξ N > 0 from (51), i.e., u ¯ 0 is larger than u ¯ 1 for N N e x t . Note that the sign of u ¯ 0 is negative according to (49). Therefore, the maximum of dimensionless guidance command is
u ¯ M a x = u ¯ 0 = u ¯ 0
(2) If 1 < κ 0 < 2 / 3 , the set N e x t = 0 , N 1 N 2 , + . According to Figure 2, Ξ N divides into three parts according to different weighted gain, which is
Ξ N = N + 2 f a d d N , N 0 , N 1 N + 2 f s u b N , N N 2 , N 3 N + 2 f a d d N , N N 3 , +
Then, the sign of Ξ N can be derived according to the Propositions A1 and A2, which is
Ξ N = u ¯ 0 u ¯ 1 > 0 , N 0 , N 1 u ¯ 0 u ¯ 1 < 0 , N N 2 , N z e r o a d d u ¯ 0 u ¯ 1 > 0 , N N z e r o a d d , +
where N z e r o a d d is the zero point of function f a d d N . It should be noted that there is a special case, N N 1 , N 2 , in which the maximum of dimensionless guidance command is u ¯ 0 according to (48). Additionally, the sign of u ¯ 0 and u ¯ 1 is dependent on the value of C ¯ 2 which can be easily conclude that
C ¯ 2 κ 0 + 1 N 2 + 3 κ 0 + 2 < 1 , N 0 , N 2 C ¯ 2 > κ 0 + 1 N 2 + 3 κ 0 + 2 > 0 , N N 2 , N z e r o a d d C ¯ 2 > κ 0 + 1 N 3 + 3 κ 0 + 2 = 1 , N N z e r o a d d , +
As a result, the maximum of dimensionless guidance command is
u ¯ M a x = u ¯ 0 = u ¯ 0 , N 0 , N 2 u ¯ 1 = u ¯ 1 , N N 2 , N z e r o a d d u ¯ 0 = u ¯ 0 , N N z e r o a d d , +
(3) If 2 / 3 κ 0 < 1 / 2 , the set N e x t is N 2 , + . Because the derivation process is similar to the case of 1 < κ 0 < 2 / 3 , the sign of Ξ N is
Ξ N = u ¯ 0 u ¯ 1 < 0 , N N 2 , N z e r o a d d u ¯ 0 u ¯ 1 > 0 , N N z e r o a d d , +
It should be noted that there is also a special case, N 0 , N 2 , in which the maximum of the dimensionless guidance command is u ¯ 0 according to (48). Thus, the maximum of the dimensionless guidance command is the same as that presented in (60).
(4) If κ 0 1 / 2 , the set N e x t = 0 , + . It can be easily obtained that for any N 0 , + ,
C ¯ 2 = κ 0 + 1 N + 3 κ 0 + 2 3 κ 0 + 2 > 0
Furthermore, it is clear that if κ 0 1 / 3 , C ¯ 2 > 1 for any N . As for the case 1 / 2 κ 0 < 1 / 3 , it can be obtained from Figure 2 that
C ¯ 2 < 1 , N 0 , N 3 C ¯ 2 > 1 , N N 3 , +
where N 3 is a function of κ 0 , which is presented in (55). Thus, Ξ N can be expressed as
Ξ N = N + 2 f s u b N , N 0 , m a x 0 , N 3 N + 2 f a d d N , N m a x 0 , N 3 , +
According to the Propositions A1 and A2, the sign of Ξ N can be derived.
Ξ N = u ¯ 0 u ¯ 1 < 0 , N 0 , N z e r o a d d u ¯ 0 u ¯ 1 > 0 , N N z e r o a d d , +
Thus, the maximum of dimensionless guidance command is
u ¯ M a x = u ¯ 1 = u ¯ 1 , N 0 , N z e r o a d d u ¯ 0 = u ¯ 0 , N N z e r o a d d , +
In conclusion, the closed-form solution of u ¯ M a x for various initial states are determined by (56), (60), and (66). Noted that, though u ¯ M a x is the piecewise function of N , it is still continuous.

4.2.2. Weighted Gain Minimal Maximum of Guidance Command

In order to find the weighted gain corresponding to the minimal of u ¯ M a x , the derivative of u ¯ M a x to N is obtained. There are still three types according to the different initial states.
(1) If κ 0 1 , it can be obtained from (56) that the derivative of u ¯ M a x to N is
u ¯ ˙ M a x = 2 κ 0 + 1 N 5 κ 0 + 3
It is obvious that u ¯ ˙ M a x > 0 for any N 0 , i.e., u ¯ M a x is the monotonic increasing function of N . Therefore, the weighted gain which minimizes the maximum of guidance command is
N min = 0
(2) If 1 < κ 0 < 1 / 2 , according to (60), the derivative of u ¯ M a x to N is
u ¯ ˙ M a x = 2 κ 0 + 1 N 5 κ 0 + 3 , N 0 , N 2 u ¯ 1 Γ N / N + 2 , N N 2 , N z e r o a d d 2 κ 0 + 1 N + 5 κ 0 + 3 , N N z e r o a d d , +
where
Γ N = 1 + C ¯ 1 / C ¯ 2 + N + 2 f N + 1
If N N z e r o a d d , + , it is obvious that u ¯ ˙ M a x is the monotonic increasing function of N . According to the Proposition A1, N 3 is larger than N z e r o a d d . It is easy to find that u ¯ ˙ M a x in interval N z e r o a d d , + is larger than zero.
u ¯ ˙ M a x > u ¯ ˙ M a x N z e r o a d d > u ¯ ˙ M a x N 3 = 1 κ 0 > 0
If N N 2 , N z e r o a d d , the sign of u ¯ ˙ M a x is depended on Γ N . It can be seen from the Proposition A3 presented in Appendix A, there is one and only one zero point of Γ N defined as N e x t u 1 in the interval N 2 , + for the case 1 < κ 0 < 1 / 2 . In order to analyze the sign of u ¯ ˙ M a x N in interval N N 2 , N z e r o a d d , it is necessary to study the relationship between N e x t u 1 and N z e r o a d d . Substituting N 3 into (70), the value of Γ N at N 3 is
Γ N 3 = 2 x ln 1 x + 2 7 κ 0 x > 2 x ln 1 x + 2 x > 0
where, x = κ 0 + 1 / 2 κ 0 0 , 1 / 2 . Note that, because Γ N is the monotonic increasing function and N e x t u 1 is the zero point of Γ N , N e x t u 1 < N 3 . Furthermore, it is clear that N 3 < N z e r o a d d because N z e r o a d d is the zero point of f add N in the interval max 0 , N 3 , + . It can be obtained that N e x t u 1 < N 3 < N z e r o a d d . Therefore, it is easy to find that u ¯ ˙ M a x < 0 when N N 2 , N e x t u 1 and u ¯ ˙ M a x > 0 when N N e x t u 1 , N z e r o a d d .
If N 0 , N 2 , it is clear that there is one and only one zero point of u ¯ ˙ M a x N .
N e x t u 0 = 5 κ 0 / 2 + 3 / 2 / κ 0 + 1
Moreover, u ¯ ˙ M a x > 0 , when N 0 , N e x t u 0 and u ¯ ˙ M a x < 0 , when N N e x t u 0 , N 2 . In conclusion, the sign of u ¯ ˙ M a x is derived as
u ¯ ˙ M a x > 0 , N 0 , N e x t u 0 u ¯ ˙ M a x < 0 , N N e x t u 0 , N e x t u 1 u ¯ ˙ M a x > 0 , N N e x t u 1 , +
Consequently, the minimum of u ¯ M a x is at N = 0 or N = N e x t u 1 . The weighted gain which minimizes the maximum of guidance command is
N min = 0 , u ¯ M a x 0 u ¯ M a x N e x t u 1 N e x t u 1 , u ¯ M a x 0 > u ¯ M a x N e x t u 1
It can be seen from Figure 3 that, there is only one point, κ 0 = 0.7745 , in the interval 1 , 1 / 2 such that u ¯ M a x 0 = u ¯ M a x N e x t u 1 . Therefore, (75) can be rewritten as
N min = 0 , 1 < κ 0 0.7745 N e x t u 1 , 0.7745 < κ 0 < 1 / 2
(3) If κ 0 1 / 2 , according to (66), the derivative of u ¯ M a x to N is
u ¯ ˙ M a x = u ¯ 1 Γ N / N + 2 , N 0 , N z e r o a d d 2 κ 0 + 1 N + 5 κ 0 + 3 , N N z e r o a d d , +
It is obvious that u ¯ ˙ M a x > 0 for any N N z e r o a d d , + . Furthermore, if κ 0 > 1 / 2 , it can be obtained from Proposition A3 that Γ N < 0 when N 0 , N e x t u 1 and Γ N > 0 when N N e x t u 1 , + . If N e x t u 1 < N z e r o a d d , the sign of u ¯ ˙ M a x is
u ¯ ˙ M a x < 0 , N 0 , N e x t u 1 u ¯ ˙ M a x > 0 , N N e x t u 1 , +
Thus, N e x t u 1 is the minimal point of u ¯ M a x . And if N e x t u 1 > N z e r o a d d , the sign of u ¯ ˙ M a x is
u ¯ ˙ M a x < 0 , N 0 , N z e r o a d d u ¯ ˙ M a x > 0 , N N z e r o a d d , +
Thus, N z e r o a d d is the minimal point of u ¯ M a x . Therefore, for the case κ 0 > 1 / 2 , the weighted gain which minimizes the maximum of guidance command is
N min = min N e x t u 1 , N z e r o a d d
It can be seen from Figure 4, there is only one point, κ 0 = 0.0634 , in the set 1 / 2 , + such that N e x t u 1 = N z e r o a d d . Therefore, (80) can be rewritten as
N min = N e x t u 1 , 1 / 2 < κ 0 0.06335 N z e r o a d d , κ 0 > 0.06335
If κ 0 = 1 / 2 , it can be obtained from the Proposition A3 that Γ N > 0 when N 0 , + . Thus, the sign of u ¯ ˙ M a x is positive for any N . The weighted gain which minimizes the maximum of guidance command is
N min = 0
In conclusion, the weighted gain which minimizes the maximum of guidance command is determined by (68), (76), and (81), which can be summarized as:
N min = 0 , κ 0 0.7745 , κ 0 = 0.5 N e x t u 1 , 0.7745 < κ 0 0.0634 & κ 0 0.5 N z e r o a d d , κ 0 > 0.06335

4.3. Total Maneuvering Energy

For the missile guidance and control system, it is desired to impact the target with less maneuvering in the course of guidance. Therefore, a special function E is defined to represent the energy consumption caused by maneuvering in the whole course of guidance.
E = 0 1 u ¯ 2 t ¯ d t ¯
Substituting (40) into (84), the expression of E can be obtained as
E N = N + 2 2 2 N + 1 2 N + 3 κ 0 + 1 2 N 3 + 7 κ 0 2 + 10 κ 0 + 3 N 2                             + 5 3 κ 0 2 + 3 κ 0 + 1 N + 3 3 κ 0 2 + 3 κ 0 + 1
It is obvious that E is related to the initial value κ 0 and the weighted gain N . Therefore, it is necessary to analyze the minimum of E .
The derivative of E to N is
E ˙ N = N N + 2 2 N + 1 2 2 N + 3 2 Θ a N , κ 0
where
Θ a N , κ 0 = a 2 N κ 0 2 + a 1 N κ 0 + a 0 N
And
a 2 N = 12 N 4 + 96 N 3 + 275 N 2 + 334 N + 147 a 1 N = 24 N 4 + 160 N 3 + 394 N 2 + 436 N + 183 a 0 N = 12 N 4 + 64 N 3 + 139 N 2 + 142 N + 57
It is obvious that the sign of E ˙ N is determined by the quadratic polynomial Θ a N , κ 0 . The trajectory of function Θ a N , κ 0 = 0 with different N and κ 0 is shown in Figure 5. It is obvious that it divides space N , κ 0 into two regions, S 1 and S 2 , the one is E ˙ N < 0 and the other is E ˙ N > 0 .
As can be seen from Figure 5, there is a maximum κ max of κ 0 in the trajectory of Θ a N , κ 0 = 0 which can be regarded as an implicit function of κ 0 with respect to N . Calculating the derivative of both sides of Θ a N , κ 0 = 0 to N , it yields
2 a 2 N κ 0 d κ 0 / d N + a 1 N d κ 0 / d N + a ˙ 2 N κ 0 2 + a ˙ 1 N κ 0 + a ˙ 0 N = 0
where, d κ 0 / d N , a ˙ 2 N , a ˙ 1 N , and a ˙ 0 N are the derivatives of κ 0 , a 2 , a 1 , and a 0 to N , respectively. It is obvious that κ max is at the point which makes d κ 0 / d N zero. Therefore, according to (89), κ max must meet the equation as
a ˙ 2 N κ max 2 + a ˙ 1 N κ max + a ˙ 0 N = 0
Because a 2 , a 1 , and a 0 are polynomial functions of N , those derivatives can be determined in analytical manner. Additionally, because κ max is at the trajectory of Θ a N , κ 0 = 0 , the following equation should be satisfied for κ max and N .
a 2 N κ max 2 + a 1 N κ max + a 0 N = 0
It is easy to solve the equation set (90) and (91) to obtain κ max = 0.6418 .
Obviously, if κ 0 κ max , E ˙ N 0 for any N , i.e., E N is the monotonic increasing function of N . And if κ 0 < κ max , two weighted gains satisfy Θ a N , κ 0 = 0 , which is defined as N E 1 and N E 2 (might as well assume that N E 1 < N E 2 ). Then, the sign of E ˙ N is derived as
E ˙ N > 0 , N 0 , N E 1 N E 2 , + E ˙ N < 0 , N N E 1 , N E 2
It is apparent that E N increases with the increase of N when N 0 , N E 1 , decreases with the increase of N when N N E 1 , N E 2 , and increases again when N N E 2 , + . Therefore, the minimum of E N is at N = 0 or N = N E 2 , where the values of E N are
E 0 = 4 3 κ 0 2 + 3 κ 0 + 1 , E 2 = E N E 2 ; κ 0
In order to determine the minimum of E 2 and E 0 , the difference between them is defined as
E 2 E 0 = N E 2 2 2 N E 2 + 1 2 N E 2 + 3 Θ b N E 2 , κ 0
where
Θ b N E 2 , κ 0 = b 2 N E 2 κ 0 2 + b 1 N E 2 κ 0 + b 0 N E 2
and
b 2 N = N 3 + 11 N 2 + 47 N + 49 b 1 N = 2 N 3 + 18 N 2 + 63 N + 61 b 0 N = N 3 + 7 N 2 + 21 N + 19
It is obvious that the sign of E 2 E 0 is determined by the quadratic polynomial Θ b N E 2 , κ 0 , the discriminant of which satisfies
Δ = N E 2 1 2 4 N E 2 2 + 8 N E 2 + 3 < 0
Thus, it leads to the fact that Θ b N E 2 , κ 0 > 0 because the highest degree coefficient of Θ b N E 2 , κ 0 is greater than zero. Furthermore, it can be obtained that E 2 > E 0 . It means that E 0 is the minimum value of E .

4.4. Weighted Gain for Limiting Maximum Guidance Command

In this subsection, a systematic method based on the characteristic presented in the previous two subsections is proposed for limiting maximum guidance command by adjusting the weighed gain N . Note that if the initial conditions are given, the maximum of guidance command can be expressed as a function of weighed gain N . Therefore, the set of weighted gain N that meets the acceleration limit can be determined inversely. Additionally, the total maneuvering energy, which is an important index especially for the missile’s maneuverability with direct force in exo-atmospheric interception, is also dependent on the weighted gain. Therefore, it is necessary to select the weighted gain which minimizes the total maneuvering energy while satisfying the acceleration constraint.
It can be seen from last subsection that when N = 0 , the value of E is at its minimum, and the navigation gains N Z E M , N Z E A E are the minimum. It is beneficial for the missile to increase the robustness to the error. Obviously, if the maximum of guidance command for the case N = 0 , u M a x 0 , is less than the acceleration limit, u lim , the weighted gain should be selected as
N = 0
If u M a x N min < u lim < u M a x 0 , there exists an interval N lim 1 , N lim 2 satisfying u M a x N < u lim . The weighted gain should be selected as the one makes E minimal in this interval. If κ 0 κ max , E is the monotonic increasing function of N according to the previous subsection, and therefore the weighted gain should be selected as
N = N lim 1
And if κ 0 < κ max , E is also a monotonic increasing function of N for N < N E 1 or N > N E 2 , and a monotonic decreasing function for N E 1 < N < N E 2 . Obviously, the minimum of E in interval N lim 1 , N lim 2 may occur at N lim 1 , N lim 2 , or N E 2 . Therefore, the weighted gain should be selected as the minimal values of E at these three points.
N = N j min E N j , j lim 1 , lim 2 , E 2
Let us consider a stricter case, u M a x N min u lim . It is well known that acceleration saturation will generally occur during the initial period. Therefore, the weighted gain N min should be updated at each guidance step until the current guidance command is within the acceleration limit. However, If the weighted gain is selected as N min , the acceleration will still be saturated near the interception. A varying weighted gain N should be employed to avoid acceleration saturation near the interception so as to improve guidance performance. For this case, the weighted gain should be re-selected to zero at the point κ = 1 / 2 , which can significantly reduce the maximum of guidance command according to the Proposition A4 presented in Appendix A.
The weighted gain satisfying the acceleration limit is obtained. Applying it to the guidance command in the form of state feedback presented in (28), a complete guidance law with impact-angle constraint and acceleration limit is successfully derived. It is noted that although this method has been developed for lag-free systems, it may be applied to lag systems as well.

4.5. Implementation of Proposed Guidance Law

In this paper, the time-to-go weighted optimal guidance law with impact-angle constraint and acceleration limit is developed. The procedure for implementing the proposed guidance law is included in the flow chart in Figure 6 and summarized in detail as follows.
(1)
Initialization: Set initial simulation parameters.
(2)
Linearization of relative motion model: Construct the linear model for guidance law derivation based on the current value of V M , V T , γ M , γ T , and λ .
(3)
Determine whether to update the weighted gain according to the following conditions:
  • The current guidance step is the first one.
  • There is the case mentioned in the previous subsection where the weighted gain needs to be updated:
    i:
    The current guidance command generated by the current weighted gain is not within the acceleration limit, update the weighted gain as N min .
    ii:
    If u max N min > u lim , re-select the weighted gain to zero at the point κ = 1 / 2 .
If yes, go to next step, if not, go to step 7.
(4)
Derive the closed-form expression to guidance command of the lag-free dynamic.
(5)
Calculate the maximum of guidance command and the total maneuvering energy based on the closed-form expression to guidance command.
(6)
Select the optimal weighted gain that minimal the total maneuvering energy while satisfying the acceleration constraint.
(7)
Using (28), generate the current guidance command with the selected weighted gain.
(8)
Apply the guidance command to the real nonlinear relative motion model and update the initial state, then return step 2.

5. Simulation Results

In this section, numerical simulations using the nonlinear engagement kinematics are carried out to evaluate the performance of the proposed guidance law. Three cases with different initial flight-path angle and impact-angle constraint are considered, which are: (a) case 1: γ M 0 = 0 , γ f c = 0 ; (b) case 2: γ M 0 = 20 , γ f c = 0 ; (c) case 1: γ M 0 = 20 , γ f c = 20 . The initial parameters κ 0 corresponding to these three cases are 0, −0.6, and −1.7, respectively, which are typical cases for the category of κ 0 . And the other required initial conditions are listed in Table 1. A first-order lag with a time constant of 0.3 s is considered for the interceptor. All programs are implemented on a computer with a 2.6 GHz processor and 2016a MATLAB environment.

5.1. Normal Simulation of Proposed Guidance Law

Numerical simulation is carried out in this subsection in which different initial conditions and weighted gains are considered to verify the previous theoretical results. Figure 7 provide the simulation results of instantaneous ZEM and ZEAE obtained by the proposed guidance law with N = 0 , 0.5 , 1 . It is obvious that both instantaneous ZEM and ZEAE converge to zeros at the time of impact, which means that the proposed guidance law is capable of guiding the missile to intercept the target with desired terminal angle for N 0 .
Figure 8 shows the histories of the missile accelerations. Obviously, for κ 0 = 0 and κ 0 = 1.7 , there exists a critical point for the cases with N = 0.5 and N = 1 ; for κ 0 = 0.6 there exists a critical point for the case with N = 1 . The simulation results are consistent with the theory presented in Section 4.2.1, which implies that, if N N e x t , there is a critical point in the acceleration profile.
In order to further verify the validity of the proposed method for limiting maximum acceleration, some comparisons with linear lag-free system, nonlinear lag-free system, as well as nonlinear first-order lag system have been done. The simulation results are demonstrated in Figure 9, Figure 10 and Figure 11. The subfigure (a), (b), and (c) respectively show the maximum accelerations of the three system with various N , the dimensionless time-to-go corresponding to the maximum acceleration, and the histories of the real acceleration for three systems, in which the weighted gain is selected as N = N min .
Obviously, the maximum acceleration of nonlinear simulation is approximate to, in most cases even slightly less than, that of linear lag-free system except for the case 2. The reason is that the first-order lag system is considered in the simulation, which leads to a lag for the response of command tracking. Therefore, it can be seen from Figure 9c and Figure 11c, the maximum acceleration at the beginning, is less than that in linear lag-free system. Accordingly, for the case 1, the acceleration at the critical point t 1 is also less than that in linear lag-free system. As for the case 2, the maximum acceleration is a little bigger than its theoretical value due to nonlinear and first-order lag. This may result in acceleration saturation. But the maximum for the acceleration is happened at the beginning of the flight. And then, it shapely reduces to the non-saturation. Therefore, the acceleration saturation will not have much influence on the terminal miss. The phenomenon can be inferred from Figure 10b,c. Conclusively, the method derived from the linear lag-free systems is still valid for the nonlinear lag systems.

5.2. Comparison with Other Guidance Law

The OGL-CTIA proposed in [10] is a widely used and effective guidance law to impose a specified impact-angle. In this subsection, a comparison between the proposed guidance law and OGL-CTIA is provided to demonstrate the superior performance of the proposed guidance law in limiting maximum acceleration. In the simulations, the engagement condition is set to case 1 and the various acceleration limits are involved.
Figure 12 shows the interception trajectories of both guidance laws with u lim = 65 , 75 , 85   m / s 2 . As can be seen from those figures, both guidance laws perform well when the acceleration limit is 85   m / s 2 . If the acceleration limit is stricter to 65   m / s 2 or 75   m / s 2 , the proposed method also performs well. However, the OGL-CTIA fails to intercept the target. The statistics of the miss distances and the impact-angle errors are listed in Table 2. The computational time at every guidance step of the proposed method and OGL-CTIA is presented in Figure 13. It is obvious that the computational time of the proposed method is almost same as that of OGL-CTIA. Therefore, the increase of computational burden due to solving the optimal weighted gain is less and negligible, so the real-time computational feasibility of the proposed method can be guaranteed.
The missile’s accelerations for different guidance laws are given in Figure 14. It is obvious that the accelerations of all cases for the proposed guidance law are within the limiting value throughout the homing phase. However, the acceleration for OGL-CTIA saturates during the last seconds when the acceleration limits are restricted to 75   m / s 2 and 65   m / s 2 . That leads to the failure to intercept the target. It should be noted that only one case is successful because the acceleration limiting value is only 85   m / s 2 , which is much looser than that in other cases. Conclusively, the proposed guidance law improves the effectiveness of direct hit in comparison with other guidance laws. Note that a saltation in the acceleration profile of the proposed guidance law occurred for the case with u lim = 65   m / s 2 . The reason is that, at this point, N suddenly becomes zero in response to the coming acceleration saturation, which is shown in Figure 15.
The miss distance, impact-angle error, and the total maneuvering energy for different acceleration limiting values are displayed in Figure 16. Obviously, the proposed guidance law can achieve successful interception with desired impact-angle if the acceleration limiting value is as low as 62 m/s2. However, OGL-CTIA will fail if the acceleration limit is below 82 m/s2. This is because the acceleration saturates during the last seconds when the acceleration limits are stricter, causing large miss distance and terminal angle error. Although the total maneuvering energy has increased, it is not too much, by only 11.5%. Therefore, it is easy to conclude that the proposed guidance law has a superior performance in dealing with the stricter acceleration limit in comparison with typical OGL-CTIA.

6. Conclusions

A guidance law is proposed for imposing a predetermined impact-angle with the acceleration limit using linear optimal control methodology. This guidance law introduces a time-varying weighted function at the cost of control energy so that the acceleration profile can be regulated by adjusting the weighted gain. Furthermore, the relationship between the maximum of guidance command, the total maneuvering energy and the weighted gain is derived in the framework of linear lag-free system. And a method to select the suitable weighted gain to limit the maximum acceleration is developed. Several numerical simulations with linear lag-free system, nonlinear lag-free system, and nonlinear lag system are carried out to verify the performance of the proposed guidance law. The results demonstrate that, even when the scenario is initiated with largely different impact-angle requirement, the guidance laws exhibited excellent performance by providing near-zero miss distance and intercept angle error. Additionally, the proposed guidance law also has superior performance in dealing with the stricter acceleration limit in comparison with typical OGL-CTIA. Benefiting from this character, the proposed method greatly improves the interceptor’s warhead lethality, resulting in possible warhead size reduction.

Author Contributions

Conceptualization, S.Z. and W.C.; methodology, S.Z. and L.Y.; software, S.Z.; validation, S.Z., W.C. and L.Y.; investigation, S.Z. and L.Y.; writing—original draft preparation, S.Z.; writing—review and editing, L.Y.; funding acquisition, L.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the National Natural Science Foundation of China (NSFC) under Grant 62003019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proposition A1.
If   κ 0 1 , f a d d N < 0 in set  N add ; if  κ 0 > 1 , there is one and only one weighted gain  N z e r o a d d in set  N add such that  f a d d N = 0 and
f a d d N 0 , N N add | N N z e r o a d d f a d d N > 0 , N N add | N > N z e r o a d d  
Proof of Proposition A1.
(1) If κ 0 1 , the set N add = N e x t = 0 , + . It can be easily obtained that for any N 0 , +
C ¯ 2 = N + 3 κ 0 + N + 2 N + 3 + N + 2 = 1
Then, the sign of f a d d N can be derived as
f a d d N = C ¯ 2 1 C ¯ 2 t ¯ g o 1 N < C ¯ 2 1 C ¯ 2 < 0
(2) If κ 0 > 1 , let us firstly consider the case that N max 0 , N 3 , . For simplification, a function of N is defined as
g a d d N f a d d / C ¯ 2 = 1 1 / C ¯ 2 t ¯ g o 1 N , N max 0 , N 3 ,
Its derivative to N is
g ˙ a d d = κ 0 + 1 / C ¯ 2 2 t ¯ g o 1 N f N
where
f N = N 1 N + 1 N N 0 + 1 1 N + 3 1 N N 0 + ln N N N 0 + 1 N + 3 N N 0 , N 0 = 2 κ 0 + 1 / κ 0 + 1
If κ 0 > 1 , it is obvious that N 0 2 , + . In order to get the sign of f N , let us introduce a function
M a = ln N N a N N a ; N a > 0
The derivative of M a to a is
d M d a = a N a 2 = > 0 , a < 0 < 0 , a > 0
Using M a , the function f N can be rewritten as
f N = M 3 + M N 0 M N 0 1 M 0
If N 0 2 , 0 , 3 < N 0 1 < N 0 < 0 , then
f N = M 3 M N 0 1 + M N 0 M 0 < 0
If N 0 0 , 1 , 3 < N 0 1 < 0 < N 0 , then
f N = M 3 M N 0 1 + M N 0 M 0 < 0
And if N 0 1 , + , 3 < 0 < N 0 1 < N 0 , then
f N = M N 0 M N 0 1 + M 3 M 0 < 0
Therefore, it is obvious that f N < 0 for any κ 0 > 1 and N N add . Then, it is easy to find that g a d d is a monotonic function and its derivative is larger than zero.
g ˙ a d d = κ 0 + 1 / C ¯ 2 2 t ¯ g o 1 N f N > 0
when N + , the limit value of g a d d can be found, which is larger than zero.
g a d d + = lim N + 1 1 / C ¯ 2 t ¯ g o 1 N = 1 e 2 > 0
If κ 0 1 / 3 , N 3 0 , the left boundary value of g a d d N is
g a d d 0 = lim N 0 1 1 / C ¯ 2 t ¯ g o 1 N = 1 / 3 κ 0 + 2 < 0
And if 1 < κ 0 < 1 / 3 , N 3 > 0 , the left boundary value of g a d d N is
g a d d N 3 = t ¯ g o 1 N < 0
Conclusively, there is one and only one N = N z e r o a d d in interval max 0 , N 3 , such that f a d d N = 0 according to Zero Point Theorem. Moreover, the sign of f a d d N can be easily derived as
f a d d N 0 , N max 0 , N 3 , N z e r o a d d f a d d N > 0 , N N z e r o a d d , +
It should be noted that if 1 < κ 0 < 2 / 3 , N add include an interval 0 , N 1 apart from N 3 , + . For that case, C ¯ 2 < 0 . Therefore,
f a d d N = C ¯ 2 1 t ¯ g o 1 N 1 < 0
because t ¯ g o 1 < 1 , and therefore, there is no zero point of f a d d N in interval 0 , N 1 . □
Proposition A2.
If  κ 0 > 1 , for any  N N sub , f s u b N is always larger than zero.
Proof of Proposition A2.
For simplifying the derivation, κ 0 is regard as the independent variable. The derivative of f s u b to κ 0 is
d f s u b d κ 0 = N + 3 + t ¯ g o 1 N N + 3 N C ¯ 1
For any N N sub = max 0 , N 2 , N 3 , it is easy to find the expression of the minimum of κ 0 according to (43).
κ min = N + 3 / 2 / N + 3
Substitute (A20) into (37), the minimal value of C ¯ 1 is obtained as
C ¯ 1 , min = N + 2 κ min + N + 1 = N / 2 N + 3
Next, substitute (A21) into (A19), we can get the lower boundary of the derivative of f s u b to κ 0 , which is large than zero since t ¯ g o 1 N is within 0 , 1 .
d f s u b d κ 0 > N + 3 + t ¯ g o 1 N N + 3 N C ¯ 1 , min = N + 3 1 t ¯ g o 1 N > 0
It is obvious that f s u b is a monotonic increasing function of κ 0 . Furthermore, substituting (A20) into (53), the value of f s u b for the case κ 0 = κ min is
f s u b κ min = 0
Therefore, it is obvious that
f s u b N > f s u b κ min = 0
 □
Proposition A3.
If  κ 0 > 1 and  κ 0 1 / 2 , there is one and only one zero point of  Γ N in the interval  max 0 , N 2 , + . And, if  κ 0 = 1 / 2 , there is no zero point of  Γ N in the interval  max 0 , N 2 , + .
Proof of Proposition A3
. If κ 0 > 1 , according to (70), the derivative of Γ N to N is
Γ ˙ N = ln N N N 0 + 1 N + 3 N N 0 + 2 N N 0 + 2 N N 0 + N 0 3 N 0 + 2 N N 0 2 + N 0 + 1 N N 0 + 1 N 0 N 0 + 1 N N 0 + 1 2 + 1 N + 3 + 3 N + 3 2
where
N 0 = 2 κ 0 + 1 / κ 0 + 1
For simplifying the derivation, N 0 is also regard as the independent variable and Γ ˙ N is defined as Λ N , N 0 . According to (A26), it is obvious that N 0 is within 2 , + if κ 0 > 1 . Furthermore, in order to satisfy N max 0 , N 2 , the variables N and N 0 must meet the inequality, 2 N 3 N 0 0 . The partial derivative of Λ N , N 0 to N 0 is
Λ N 0 = N 0 5 N + N 0 + 4 N N 0 3 N + N 0 + 3 N N 0 + 1 3
In order to get the sign of Λ / N 0 , let define
Ω N , N 0 = 5 N + N 0 + 4 N N 0 3 N + N 0 + 3 N N 0 + 1 3
It is obvious that 0 < N N 0 < N N 0 + 1 if 2 N 3 N 0 0 . And then, Ω N , N 0 is large than zero.
Ω N , N 0 > 5 N + N 0 + 4 N N 0 + 1 3 N + N 0 + 3 N N 0 + 1 3 = 4 N + 1 N N 0 + 1 3 > 0
Therefore, if N 0 2 , 0 , Λ N , N 0 is decreasing in the direction of N 0 , and if N 0 0 , + , Λ N , N 0 is increasing, i.e., N 0 = 0 is the minimum point of Λ N , N 0 .
Λ N , N 0 Λ N , N 0 = 0
Subsisting N 0 = 0 into (A25), the value of Λ N , N 0 = 0 is
Λ N , N 0 = 0 = ln N + 1 N + 3 + 1 N + 1 + 1 N + 3 + 3 N + 3 2
The derivative of Λ N , N 0 = 0 to N is
d Λ N , N 0 = 0 d N = 3 N + 1 2 N + 3 2 6 N + 3 3 < 0
Therefore, it is obvious that
Λ N , N 0 = 0 > lim N Λ N , N 0 = 0 = 0
It can be easily obtained from (A30) and (A33) that Γ ˙ N = Λ N , N 0 > 0 , i.e., Γ N is the monotonic increasing function. When N + , the limit value of Γ N can be found, which is bounded by two.
Γ + = lim N + Γ N = 2
According to (70), if κ 0 > 1 and κ 0 1 / 2 , the left boundary value of Γ N is
Γ max 0 , N 2 = lim N 0 Γ N = , κ 0 > 1 / 2 lim N N 2 Γ N = 2 κ 0 < 0 , 1 < κ 0 < 1 / 2
Conclusively, there is one and only one zero point, N ext u 1 , of Γ N in interval max 0 , N 2 , according to Zero Point Theorem. Moreover, the sign of Γ N can be easily derived as
Γ N 0 , N max 0 , N 2 , N ext u 1 Γ N > 0 , N N ext u 1 , +
If κ 0 = 1 / 2 , the left boundary value of Γ N is
Γ max 0 , N 2 = lim N 0 Γ N = 3 2 ln 3 > 0
Therefore, there is no zero point of Γ N in interval max 0 , N 2 , because Γ N > Γ max 0 , N 2 > 0 . □
Proposition A4.
The maximum of guidance command can be reduced by reselecting  N = 0 when κ = 1 / 2 .
Proof of Proposition A4.
Substituting (36) into (35), the parameter κ relate with the optimal trajectory is derived as
κ = 1 N + 3 N + 2 C ¯ 2 N + 2 C ¯ 2 N + 3 C ¯ 1 t ¯ g o
κ = 1 / 2 happens at t 2 where the value of t ¯ g o is
t ¯ g o 2 = N C ¯ 2 / N + 1 C ¯ 1 > t ¯ g o 1
Substituting (A39) into (36), it yields
Z 2 t 2 = 2 C ¯ 1 t ¯ g o 2 N + 1 Z 20 / N + 1
According to (39), if the weighted gain is reselected to zero at time t 2 , it is easy to find that
u ¯ t = 1
Then, the subsequent guidance command is a constant.
u c = V M Z 2 t 2 / t ¯ g o 2
Also, a variable η u u 1 / u c is defined to denote the difference between u 1 and u c , which can be expressed as
η u = 1 2 N + 1 N + 2 N + 1 N + 3 N
and its derivative to N is
d η u / d N = η u Π N / 2
where
Π N = ln N + 1 ln N + 3 + 1 / N + 2 + 3 / N + 3
It is obvious that η u > 0 from (A43), thus the sign d η u / d N is determined by Π N . The derivative of Π N to N is
d Π d N = 2 N 3 + 8 N 2 + 7 N 3 N + 1 N + 2 2 N + 3 2
As analyzing from the numerator, it is obvious that there is one and only on zero point of d Π / d N , and Π N will increase at the begging and then decrease as N increases. Therefore, the minimum of Π N is at boundary points, N = 0 or N = + .
Π N min Π 0 , Π + = min 3 / 2 ln 3 , 0 = 0
Finally, it is concluded that
η u η u 0 = 1 u c * < u 1 *
It also means that the maximum of guidance command is reduced. □

References

  1. Han, T.; Hu, Q.; Xin, M. Analytical solution of field-of-view limited guidance with constrained impact and capturability analysis. Aerosp. Sci. Technol. 2020, 97, 105586. [Google Scholar] [CrossRef]
  2. Liu, B.; Hou, M.; Feng, D. Nonlinear mapping based impact angle control guidance with seeker’s field-of-view constraint. Aerosp. Sci. Technol. 2019, 86, 724–736. [Google Scholar] [CrossRef]
  3. Bryson, A.E.; Ho, Y.-C. Applied Optimal Control; John Wiley & Sons: New York, NY, USA, 1975; pp. 154–155. [Google Scholar]
  4. Ben-Asher, J.Z.; Yaesh, I. Advances in Missile Guidance Theory; AIAA, Inc.: Reston, VA, USA, 1998; pp. 123–130. [Google Scholar]
  5. Kim, M.; Grider, K.V. Terminal Guidance for Impact Attitude Angle Constrained Flight Trajectories. IEEE Trans. Aerosp. Electron. Syst. 1973, AES-9, 852–859. [Google Scholar] [CrossRef]
  6. York, R.J.; Pastrick, H.L. Optimal Terminal Guidance with Constraints at Final Time. J. Spacecr. Rocket. 1977, 14, 381–383. [Google Scholar] [CrossRef]
  7. Pastrick, H.; York, R. Optimal control for an air defense interceptor: Part I. In Proceedings of the 1980 IEEE Southeastcon, Knoxville, TN, USA, 10–13 April 1980; pp. 229–232. [Google Scholar]
  8. Ryoo, C.-K.; Cho, H.; Tahk, M.-J. Optimal Guidance Laws with Terminal Impact Angle Constraint. J. Guid. Control Dyn. 2005, 28, 724–732. [Google Scholar] [CrossRef]
  9. Song, T.L.; Shin, S.J.; Cho, H. Impact Angle Control for Planar Engagements. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 1439–1444. [Google Scholar] [CrossRef]
  10. Shaferman, V.; Shima, T. Linear Quadratic Guidance Laws for Imposing a Terminal Intercept Angle. J. Guid. Control Dyn. 2008, 31, 1400–1412. [Google Scholar] [CrossRef]
  11. Rusnak, I.; Weiss, H.; Eliav, R. Missile guidance with constrained terminal body angle. In Proceedings of the 2010 IEEE 26th Convention of Electrical and Electronics Engineers in Israel, Eilat, Israel, 17–20 November 2010; pp. 45–49. [Google Scholar]
  12. Idan, M.; Golan, O.M.; Guelman, M. Optimal planar interception with terminal constraints. J. Guid. Control Dyn. 1995, 18, 1273–1279. [Google Scholar] [CrossRef]
  13. Chang-Kyung, R.; Hangju, C.; Min-Jea, T. Time-to-go weighted optimal guidance with impact angle constraints. IEEE Trans. Control Syst. Technol. 2006, 14, 483–492. [Google Scholar] [CrossRef]
  14. Ohlmeyer, E.J. Control of terminal engagement geometry using generalized vector explicit guidance. In Proceedings of the 2003 American Control Conference, Denver, CO, USA, 4–6 June 2003; pp. 396–401. [Google Scholar]
  15. Ohlmeyer, E.J.; Phillips, C.A. Generalized Vector Explicit Guidance. J. Guid. Control Dyn. 2006, 29, 261–268. [Google Scholar] [CrossRef]
  16. He, S.; Lee, C.-H. Optimal Impact Angle Guidance for Exoatmospheric Interception Utilizing Gravitational Effect. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1–10. [Google Scholar] [CrossRef] [Green Version]
  17. Lee, C.-H.; Tahk, M.-J.; Lee, J.-I. Generalized Formulation of Weighted Optimal Guidance Laws with Impact Angle Constraint. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1317–1322. [Google Scholar] [CrossRef] [Green Version]
  18. Xiong, S.; Wei, M.; Zhao, M.; Xiong, H.; Wang, W.; Zhou, B. Hyperbolic tangent function weighted optimal intercept angle guidance law. Aerosp. Sci. Technol. 2018, 78, 604–619. [Google Scholar] [CrossRef]
  19. Li, C.; Wang, J.; He, S.; Lee, C.-H. Collision-geometry-based generalized optimal impact angle guidance for various missile and target motions. Aerosp. Sci. Technol. 2020, 106, 106204. [Google Scholar] [CrossRef]
  20. Kim, B.S.; Lee, J.G.; Han, H.S. Biased PNG law for impact with angular constraint. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 277–288. [Google Scholar] [CrossRef]
  21. Ratnoo, A.; Ghose, D. Impact Angle Constrained Guidance Against Nonstationary Nonmaneuvering Targets. J. Guid. Control Dyn. 2010, 33, 269–275. [Google Scholar] [CrossRef]
  22. Erer, K.S.; Merttopçuoglu, O. Indirect Impact-Angle-Control Against Stationary Targets Using Biased Pure Proportional Navigation. J. Guid. Control Dyn. 2012, 35, 700–704. [Google Scholar] [CrossRef]
  23. Kim, T.-H.; Park, B.-G.; Tahk, M.-J. Bias-Shaping Method for Biased Proportional Navigation with Terminal-Angle Constraint. J. Guid. Control Dyn. 2013, 36, 1810–1816. [Google Scholar] [CrossRef]
  24. Lee, C.-H.; Kim, T.-H.; Tahk, M.-J. Interception Angle Control Guidance Using Proportional Navigation with Error Feedback. J. Guid. Control Dyn. 2013, 36, 1556–1561. [Google Scholar] [CrossRef]
  25. Kumar, S.R.; Rao, S.; Ghose, D. Sliding-Mode Guidance and Control for All-Aspect Interceptors with Terminal Angle Constraints. J. Guid. Control Dyn. 2012, 35, 1230–1246. [Google Scholar] [CrossRef]
  26. Cho, D.; Kim, H.J.; Tahk, M.-J. Impact angle constrained sliding mode guidance against maneuvering target with unknown acceleration. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1310–1323. [Google Scholar] [CrossRef]
  27. Ji, Y.; Lin, D.; Wang, W.; Hu, S.; Pei, P. Three-dimensional terminal angle constrained robust guidance law with autopilot lag consideration. Aerosp. Sci. Technol. 2019, 86, 160–176. [Google Scholar] [CrossRef]
  28. Chen, Z.; Chen, W.; Liu, X.; Cheng, J. Three-dimensional fixed-time robust cooperative guidance law for simultaneous attack with impact angle constraint. Aerosp. Sci. Technol. 2021, 110, 106523. [Google Scholar] [CrossRef]
  29. Hu, Q.; Cao, R.; Han, T.; Xin, M. Field-of-view limited guidance with impact angle constraint and feasibility analysis. Aerosp. Sci. Technol. 2021, 114, 106753. [Google Scholar] [CrossRef]
  30. Rusnak, I.; Meir, L. Optimal guidance for acceleration constrained missile and maneuvering target. IEEE Trans. Aerosp. Electron. Syst. 1990, 26, 618–624. [Google Scholar] [CrossRef]
  31. Hexner, G.; Shima, T.; Weiss, H. LQG Guidance Law with Bounded Acceleration Command. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 77–86. [Google Scholar] [CrossRef]
  32. Weiss, M.; Shima, T. Optimal Linear-Quadratic Missile Guidance Laws with Penalty on Command Variability. J. Guid. Control Dyn. 2015, 38, 226–237. [Google Scholar] [CrossRef]
Figure 1. The Planar engagement geometry.
Figure 1. The Planar engagement geometry.
Aerospace 08 00358 g001
Figure 2. The value of C ¯ 2 for the various κ 0 . (a) κ 0 , 1 , N add = 0 , + , N sub = ; (b) κ 0 1 , 2 / 3 , N add = 0 , N 1 N 3 , + , N sub = N 2 , N 3 ; (c) κ 0 2 / 3 , 1 / 3 , N add = N 3 , + , N sub = max 0 , N 2 , N 3 ; (d) κ 0 1 / 3 , + , N add = 0 , + , N sub = .
Figure 2. The value of C ¯ 2 for the various κ 0 . (a) κ 0 , 1 , N add = 0 , + , N sub = ; (b) κ 0 1 , 2 / 3 , N add = 0 , N 1 N 3 , + , N sub = N 2 , N 3 ; (c) κ 0 2 / 3 , 1 / 3 , N add = N 3 , + , N sub = max 0 , N 2 , N 3 ; (d) κ 0 1 / 3 , + , N add = 0 , + , N sub = .
Aerospace 08 00358 g002
Figure 3. u ¯ M a x 0 and u ¯ M a x N e x t u 1 for various κ 0 .
Figure 3. u ¯ M a x 0 and u ¯ M a x N e x t u 1 for various κ 0 .
Aerospace 08 00358 g003
Figure 4. N e x t u 1 and N z e r o a d d for various κ 0 .
Figure 4. N e x t u 1 and N z e r o a d d for various κ 0 .
Aerospace 08 00358 g004
Figure 5. The trajectory of Θ a N , κ 0 = 0 .
Figure 5. The trajectory of Θ a N , κ 0 = 0 .
Aerospace 08 00358 g005
Figure 6. The flowchart of implementation.
Figure 6. The flowchart of implementation.
Aerospace 08 00358 g006
Figure 7. ZEM and ZEAE. (a) case 1; (b) case 2; (c) case 3.
Figure 7. ZEM and ZEAE. (a) case 1; (b) case 2; (c) case 3.
Aerospace 08 00358 g007
Figure 8. Missile acceleration. (a) case 1; (b) case 2; (c) case 3.
Figure 8. Missile acceleration. (a) case 1; (b) case 2; (c) case 3.
Aerospace 08 00358 g008
Figure 9. Some comparisons for three systems in case 1. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Figure 9. Some comparisons for three systems in case 1. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Aerospace 08 00358 g009
Figure 10. Some comparisons for three systems in case 2. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Figure 10. Some comparisons for three systems in case 2. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Aerospace 08 00358 g010
Figure 11. Some comparisons for three systems in case 3. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Figure 11. Some comparisons for three systems in case 3. (a) Maximum acceleration; (b) Dimensionless time-to-go corresponding to the maximum acceleration; (c) History of the missile acceleration.
Aerospace 08 00358 g011
Figure 12. Interception trajectories. (a) u lim = 65   m / s 2 ; (b) u lim = 75   m / s 2 ; (c) u lim = 85   m / s 2 .
Figure 12. Interception trajectories. (a) u lim = 65   m / s 2 ; (b) u lim = 75   m / s 2 ; (c) u lim = 85   m / s 2 .
Aerospace 08 00358 g012
Figure 13. Computational time. (a) u lim = 65   m / s 2 ; (b) u lim = 75   m / s 2 ; (c) u lim = 85   m / s 2 .
Figure 13. Computational time. (a) u lim = 65   m / s 2 ; (b) u lim = 75   m / s 2 ; (c) u lim = 85   m / s 2 .
Aerospace 08 00358 g013
Figure 14. Missile acceleration. (a) u lim = 65 m / s 2 ; (b) u lim = 75 m / s 2 ; (c) u lim = 85 m / s 2 .
Figure 14. Missile acceleration. (a) u lim = 65 m / s 2 ; (b) u lim = 75 m / s 2 ; (c) u lim = 85 m / s 2 .
Aerospace 08 00358 g014
Figure 15. κ and N with u lim = 65   m / s 2 .
Figure 15. κ and N with u lim = 65   m / s 2 .
Aerospace 08 00358 g015
Figure 16. Interception results for various acceleration limiting values. (a) Miss distance; (b) Impact-angle error; (c) Total maneuvering energy.
Figure 16. Interception results for various acceleration limiting values. (a) Miss distance; (b) Impact-angle error; (c) Total maneuvering energy.
Aerospace 08 00358 g016
Table 1. Initial Conditions for Engagement.
Table 1. Initial Conditions for Engagement.
ParametersValues
Missile-target initial relative range, r 0 230 km
Initial LOS angle, λ 0
Missile initial velocity, V M 0 4000 m/s
Target initial velocity, V T 0 3000 m/s
Target initial flight-path angle, γ T 0
Gravitational acceleration, g 9.8 m/s2
Table 2. Miss Distances and Impact-angle Errors of Two Guidance Laws.
Table 2. Miss Distances and Impact-angle Errors of Two Guidance Laws.
Acceleration
Limit (m/s2)
Proposed GuidanceOGL-CTIA
Miss (m)Angle Error (deg)Miss (m)Angle Error (deg)
652.04 × 10-31.52 × 10-3257.541.26
757.42×10-38.08×10-332.880.26
853.68 × 10-34.11 × 10-33.68 × 10-34.11 × 10-3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, S.; Chen, W.; Yang, L. Optimal Guidance Law with Impact-Angle Constraint and Acceleration Limit for Exo-Atmospheric Interception. Aerospace 2021, 8, 358. https://doi.org/10.3390/aerospace8120358

AMA Style

Zhao S, Chen W, Yang L. Optimal Guidance Law with Impact-Angle Constraint and Acceleration Limit for Exo-Atmospheric Interception. Aerospace. 2021; 8(12):358. https://doi.org/10.3390/aerospace8120358

Chicago/Turabian Style

Zhao, Shilei, Wanchun Chen, and Liang Yang. 2021. "Optimal Guidance Law with Impact-Angle Constraint and Acceleration Limit for Exo-Atmospheric Interception" Aerospace 8, no. 12: 358. https://doi.org/10.3390/aerospace8120358

APA Style

Zhao, S., Chen, W., & Yang, L. (2021). Optimal Guidance Law with Impact-Angle Constraint and Acceleration Limit for Exo-Atmospheric Interception. Aerospace, 8(12), 358. https://doi.org/10.3390/aerospace8120358

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop