Next Article in Journal
Artificial Intelligence Applied to Drone Control: A State of the Art
Next Article in Special Issue
Decentralized UAV Swarm Control: A Multi-Layered Architecture for Integrated Flight Mode Management and Dynamic Target Interception
Previous Article in Journal
Lightweight Oriented Detector for Insulators in Drone Aerial Images
Previous Article in Special Issue
Multi-UAV Formation Path Planning Based on Compensation Look-Ahead Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fault-Tolerant Tracking Control of Hypersonic Vehicle Based on a Universal Prescribe Time Architecture

1
Air Traffic Control and Navigation School, Air Force Engineering University, Xi’an 710038, China
2
Equipment Management and Unmanned Aerial Vehicle Engineering School, Air Force Engineering University, Xi’an 710038, China
3
Modern Aviation College, Guangzhou Institute of Science and Technology, Guangzhou 510540, China
*
Author to whom correspondence should be addressed.
Drones 2024, 8(7), 295; https://doi.org/10.3390/drones8070295
Submission received: 10 May 2024 / Revised: 10 June 2024 / Accepted: 16 June 2024 / Published: 2 July 2024
(This article belongs to the Special Issue Distributed Control, Optimization, and Game of UAV Swarm Systems)

Abstract

:
An adaptive tracking control strategy with a prescribe tracking error and the convergence time is proposed for hypersonic vehicles with state constraints and actuator failures. The peculiarity is that constructing a new time scale coordinate translation mapping method, which maps the prescribe time on the finite field to the time variable on the infinite field, and the convergence problem of the prescribe time is transformed into the conventional system convergence problem. The improved Lyapunov function, the improved tuning function, and the adaptive fault-tolerant mechanism are further constructed. Combined with the neural network, the prescribe time tracking control of the speed subsystem and the height subsystem are realized respectively. Combined with the Barbalat lemma and Lyapunov stability theory, the boundedness of the closed-loop system is proved. The simulation results have proven that, compared with other control strategies, it can ensure that the tracking error converges to the prescribe interval in the prescribe time and meets the constraints of the whole state of the system.

1. Introduction

In recent years, the research focus on hypersonic vehicles (HSVs) [1,2,3] has gradually shifted from general controller design to solving control design problems in practical scenarios, such as executing evaluation and reconnaissance, long-range transportation and delivery, and strategic strike missions in [4,5,6]. For the research on the HSV control problem, early works studied many linear feedback control methods, and then proposed nonlinear control methods such as sliding mode control in [7,8,9,10], adaptive backstepping control in [11,12], and fault-tolerant control in [13,14,15,16]. Subsequently, intelligent control methods based on general approximators such as neural networks in [17,18,19] or fuzzy logic systems in [20,21,22,23] were widely studied to solve the nonlinear control problem of HSV with completely unknown dynamics. Furthermore, adaptive control algorithms based on disturbance observers are used to handle unknown system uncertainties in [24,25,26]. In addition to requiring the robustness of the control system, good tracking performance is also the main goal of HSV controller design, and the convergence speed and convergence accuracy are the two important indicators for evaluating the quality of tracking performance. In order to improve the convergence speed of HSV, a finite-time control algorithms was designed to adaptive anti-saturation robust for flexible aspirated HSV under the actuator saturation state in [27]; this paper designed a self-adaptive fixed-time anti-saturation compensator by introducing auxiliary variables and adjusting gain, which not only avoids the influence of tracking error on convergence characteristics, but also reduces the complex calculation burden of inversion control. The finite-time deterministic learning control problem of HSVs with model uncertainty has been studied in [28], and achieved better learning and tracking performance through two stages: offline training and online control. However, the stability time in finite time control depends on the initial conditions of the system, which limits its practical application range. To overcome this limitation, a fixed-time adaptation strategy for HSV actuator faults has been proposed in [29], and both actual constraints and model uncertainty have been taken into account, which obtained a fault-tolerant attitude control rate that allows for fault compensation to be completed within a fixed-time. A fuzzy adaptive fault-tolerant control strategy for the fast fixed time-constrained tracking problem of HSV has been proposed in [30]; by estimating the upper and lower bounds of the actuator parameters for adaptive compensation, a smaller fixed convergence time was derived, and a segmented differentiable switching control rate was introduced to avoid singularity problems. An adaptive neural network control scheme based on the integral barrier Lyapunov function has been proposed in [31] for the fixed-time tracking control problem of HSV under asymmetric time-varying angle of attack constraints, which ensures the error always converges to a bounded compact set by directly designing the asymmetric time-varying angle of attack constraint scheme. The issue of resource consumption in flight control systems has been researched in [32], and the implementation of event triggered fixed-time control for HSVs by switching dynamic event triggering mechanisms has been investigated.
Although the above significant work makes the HSV controller more adapt to practical application needs, there is still room for improvement in the following areas. Firstly, most of the current control schemes are related to finite-time or fixed-time control, and there are relatively few research results that require achieving stability within a pre-set time in specific flight missions. Secondly, for most existing control schemes that address state constraints, prior knowledge of the initial tracking conditions is required, that is, the initial tracking error must meet certain predetermined ranges. Thirdly, when the dynamic equations of hypersonic vehicle are unknown, for established finite time and fixed-time control algorithms, the tracking error can only ensure convergence to an unknown residual set, but cannot guarantee convergence to a predetermined range. Therefore, in the case where the initial tracking conditions are completely unknown, it is worthwhile to conduct in-depth research on how to design a satisfactory controller to ensure that the tracking error asymptotically converges to the predetermined range within a predetermined time. In addition, due to the complex and ever-changing external environment, it may not be possible to ensure accurate tracking of the aircraft without considering the faults of the actuator.
Based on the above discussion, this paper proposes an adaptive asymptotic tracking control scheme with prescribe convergence time and tracking error for HSVs with actuator faults. The main innovative work is as follows:
  • A time scale coordinate mapping function was introduced to ensure the asymptotic convergence of the tracking error after transformation, achieving the convergence for the original tracking error within a predetermined time.
  • An improved Lyapunov function and a class of improved tuning functions were constructed, while incorporating the Barbalat lemma (if a continuous differentiable function has a finite limit value as time approaches infinity, and its derivative is uniformly continuous, then it tends towards stabilizing at infinity) to ensure that even if the initial tracking conditions are completely unknown, the tracking error can converge to the predetermined range.
  • Two adaptive parameters were designed to estimate unknown faults, achieving the adaptive fault-tolerant control of the system and enhancing robustness during actual flight processes.
The rest of this paper is organized as follows: Section 2 establishes the longitudinal motion model of the HSV, Section 3 designs the controllers for the speed subsystem and altitude subsystem of the HSV, Section 4 provides the closed-loop stability analysis of the HSV system, Section 5 and Section 6 list and analyzes the flight simulation results, and Section 7 draws the final conclusion.

2. Longitudinal Motion Model of Hypersonic Aircraft

American scholars Bolender and Doman [33] focused on aspirated hypersonic aircraft and considered the highly integrated propulsion system in the aircraft in the Air Force Research Laboratory. They used oblique shock waves and Prandtl–Mayer expansion theory to solve the impact of oscillating bow shock waves on the propulsion system performance, and derived the motion equation of the flexible aircraft using the Lagrange equation, which captured the inertial coupling effect between rigid body acceleration and flexible body dynamics in structural dynamics, and established a representative longitudinal dynamic nonlinear physical model.
V ˙ = T cos α D m g sin γ h ˙ = V sin γ γ ˙ = L + T sin α m V g cos γ V α ˙ = Q γ ˙ Q ˙ = M I y y η ¨ i = 2 ξ i ω i η ˙ i ω i 2 η i + N i , i = 1 , , n
In the longitudinal motion model of HSV, the rigid body states V, h, γ , α and Q, respectively, represent velocity, altitude, trajectory angle (ballistic angle), angle of attack, and pitch angular velocity. The elastic state η i represents the amplitude of the i order bending mode of the fuselage, while m, g, I y y , ξ i and ω i , respectively, represent the mass, gravity acceleration, rotational inertia, damping ratio, and flexible mode frequency of the fuselage.
T q ¯ S C T , Φ ( α ) Φ + C T ( α ) + C T η η D q ¯ S C D α , δ e , δ c , η L q ¯ S C L α , δ e , δ c , η M z T T + q ¯ S c ¯ C M α , δ e , δ c , η N i q ¯ S N i α 2 α 2 + N i α α + N i δ e δ e + N i δ c δ c + N i 0 + N i η η
where T, D, L, M, and N i , respectively, represent thrust, drag, lift, pitch moment, and generalized elastic force; the parameter fitting values for aerodynamic force and moment are shown above. Among them, η = [ η 1 , η ˙ 1 , , η n , η ˙ n ] , n N + , due to the unmeasurable elastic state, it is considered as an unknown disturbance in the control law design. In addition, if the state and control input of the rigid body in the longitudinal motion model are bounded, then the elastic state is also bounded. q ¯ , S, z T , and c ¯ , respectively, represents flight dynamic pressure, reference area, thrust arm, and reference length. The approximate coefficients of the curve fitting model are expressed as
C T , Φ ( · ) = C T , Φ α 3 α 3 + C T , Φ α 2 α 2 + C T , Φ α α + C T , Φ 0 C T ( · ) = C T α 3 α 3 + C T α 2 α 2 + C T α α + C T 0 C D ( · ) = C D α 2 α 2 + C D α α + C D δ e 2 δ e 2 + C D δ e δ e + C D δ c 2 δ c 2 + C D δ c δ c + C D 0 + C D η η C L ( · ) = C L α α + C L δ e δ e + C L δ c δ c + C L 0 + C L η η C M ( · ) = C M α 2 α 2 + C M α α + C M δ e δ e + C M δ c δ c + C M 0 + C M η η C j η = C j η 1 , 0 , , C j η n , 0 , j = T , M , L , D N i η = N i η 1 , 0 , , N i η n , 0 , i = 1 , , n
Control input Φ , δ e and δ c , respectively, represents the fuel equivalence ratio, elevator deflection angle and canard deflection angle of HSV, which are implied in aerodynamic forces (moments). It is worth noting that the model adopts a duck layout to eliminate the coupling effect of lifting, so there is a relationship δ c = k e c δ e and k e c = C L δ e / C L δ c between the canard deflection angle and the elevator deflection angle, so HSV actually becomes two control inputs: the fuel equivalence ratio Φ and elevator deflection angle δ e .
The control objective of HSV is by designing and controlling the input fuel equivalence ratio Φ and elevator deflection angle δ e , the output signal speed V and height h can accurately track their respective reference commands in the longitudinal motion plane; at the same time, it ensures that, even in the event of actuator failure, the prescribe tracking accuracy can be achieved within the prescribe time set by the designer.
Because the control input fuel equivalence ratio Φ is the decisive factor affecting the thrust, the speed V in the control output signal changes according to the effect of Φ . In addition, since the elevator deflection angle δ e in the control input changes the pitch angle and track angle, the height h in the control output signal is mainly controlled by δ e . At the same time, the flexible dynamics is ignored, and for the elastic state, because it is unmeasurable, it is considered an unknown disturbance for the convenience of modeling and subsequent control law design. Thus, an uncertain simplified HSV model is obtained, which is mainly composed of five rigid body dynamic equations:
V ˙ = f V + g V Φ + d V
h ˙ = f h + g h γ + d h γ ˙ = f γ + g γ α + d γ α ˙ = f α + g α Q + d α Q ˙ = f Q + g Q δ e + d Q
where (4) is related to speed V and (5) is related to height h.
In the dynamic Equation (4) about velocity V:
f V = q ¯ S m α 3 cos α C T α 3 + α 2 cos α C T α 2 + α cos α C T α + cos α C T 0 α 2 C D α 2 α C D α δ e 2 C D δ D 2 + k e , c 2 δ D δ L 2 δ e C D e δ e + k e , c C D δ c C D 0 m g sin γ q ¯ S g V = q ¯ S m cos α α 3 C T , Φ α 3 + α 2 C T , Φ α 2 + α C T , Φ α + C T , Φ 0 d V = q ¯ S m C T η η cos α q ¯ S m C D η η + Δ V
The composite disturbance d V includes external disturbances such as gust, turbulence, and atmospheric disturbance, as well as structural flexibility caused by aerothermoelasticity, and Δ V represents the uncertainty and external disturbance in velocity dynamics.
In the four dynamic Equations (5) about height h:
f h = 0 , g h = V , d h = Δ h
f γ = q ¯ S m V C L 0 g V cos γ , g γ = q ¯ S m V C L α , d γ = q ¯ S m V C L η η + Δ γ
f α = q ¯ V α S m C L α S m C L 0 + g q ¯ cos γ , g α = 1 , d α = q ¯ S m V C L η η + Δ α
f Q = q ¯ S I y y α 3 Φ z T C T , Φ α 3 + α 2 Φ z T C T , Φ α 2 + α Φ z T C T , Φ α + Φ z T C T , Φ 0 + α 3 z T C T α 3 + α 2 z T C T α 2 + c ¯ C M α 2 + α z T C T α + c ¯ C M α + z T C T 0 + c ¯ C M 0 g Q = q ¯ S c ¯ I y y C M δ e + k e , c C M δ c d Q = z T q ¯ S I y y C T η η + q ¯ S c ¯ I y y C M η η + Δ Q
where Δ h , Δ γ , Δ α , and Δ Q represent uncertainty and external disturbance in high dynamics. Since γ of the actual cruise phase is very small, for simplicity, take sin γ = γ . In addition, the thrust control term T sin α is usually much smaller than the lift term L, so it can be ignored in the control design process.

3. Controller Design of Hypersonic Vehicle

The simplified HSV model is decomposed into velocity subsystem (including one dynamic equation of velocity V) and altitude subsystem (including four dynamic equations of altitude h, track angle γ , angle of attack α , and pitch angular velocity Q), and the control laws are designed, respectively.
At the same time, the adaptive controller is designed considering actuator failures including control failure and jamming. The form of actuator failure is expressed as follows:
Φ j ( t ) = ρ j h 1 Φ j ( t ) + ψ j h 1 , t t j h , k 1 , t j h , e 1
δ e j ( t ) = ρ j h 2 δ e j ( t ) + ψ j h 2 , t t j h , k 2 , t j h , e 2
where ρ j h 1 , ρ j h 2 0 , 1 , h represents the hth failure model of the system and ψ j h 1 , ψ j h 2 is an unknown number. t j h , k 1 and t j h , e 1 represent the occurrence and end time of the jth actuator failure in speed dynamics, and t j h , k 2 and t j h , e 2 represent the occurrence and end time of the jth actuator failure in high dynamics. Taking the speed subsystem as an example, note that (11) includes the three following cases:
  • When ρ j h 1 = 1 and ψ j h 1 = 0 , no fault occurred.
  • When 0 < ρ ̲ j h 1 ρ j h 1 ρ ¯ j h 1 < 1 and ψ j h 1 = 0 , partial actuator failure occurs.
  • When ρ j h 1 = 0 and ψ j h 1 0 , the actuator will no longer be affected by the control input, which means that the actuator fails completely.
Let us recall the following lemmas [34].
Lemma 1.
Let S ( Z ) be any nonlinear continuous function defined on the compact set Ω Z R n , and use the radial basis function neural network to approximate function S ( Z ) , then for any given ε > 0 , select a sufficiently large positive integer l to satisfy
S ( Z ) = Θ T ψ ( Z ) + ε ( Z ) , Z Ω Z
where ε ( Z ) is the approximation error of neural network and ε ( Z ) ε M , ε M i > 0 is the unknown normal number, Θ is the Θ value that makes ε ( Z ) the minimum among all Z Ω Z , i.e., Θ = arg min W R l sup Z Ω Z | S ( Z ) Θ T ψ ( Z ) | . In addition, ψ ( Z ) is selected as the common Gaussian function form
ψ ( Z ) = exp ( Z i ϕ i ) T ( Z i ϕ i ) ϖ i 2 , i = 1 , 2 , , w
where ϕ i , ϖ i , and w represent the center, width, and number of Gaussian functions, respectively.
Lemma 2.
For x and any constant ε > 0 , the following inequality
0 x < ε + x 2 x 2 + ε 2

3.1. Time Scale Coordinate Mapping

In this section, a time scale coordinate mapping is proposed, which transforms the convergence problem of prescribe time into a general asymptotic convergence problem.
The prescribe convergence time is recorded as T P , and it is required to achieve the convergence effect within the specified time interval, i.e., t 0 , T P . Using the following time scale coordinate mapping method, the specified time t 0 , T P in the finite field is mapped to the time τ [ 0 , + ) in the infinite field.
t = T P e τ 1 e τ + 1 t ( τ )
Matching the above expression to the transformed time τ leads to
d t d τ = 2 T P e τ ( e τ + 1 ) 2 λ ( τ )
It is easy to draw the following conclusion: λ ( τ ) is a monotonically decreasing bounded function, satisfying λ ( τ ) λ ¯ , where λ ¯ is a normal number.
Remark 1.
The design of time scale coordinate mapping function needs to meet the following properties: (1) t ( 0 ) = 0 ; (2) lim τ + t ( τ ) = T P ; (3) The function is differentiable and monotonically increasing, and the derivative is always positive, so as to ensure that the sign after combining with the gain function does not change. It is worth noting that the design of functions is not strictly limited to this form, as long as the above conditions are met, such as exp type, tan type, log type, etc.

3.2. Controller Design of Speed Subsystem

This section aims to design an adaptive tracking control scheme for the speed subsystem under actuator failure, and ensure that the speed tracking error converges to the range prescribed by the designer within the prescribe time. Firstly, the nonlinear dynamic equation of the speed subsystem of HSV is described as
V ˙ ( t ) = f V ( t ) + g V ( t ) j = 1 m 1 ρ j h 1 Φ j ( t ) + ψ j h 1 + d V ( t )
where the fuel equivalence ratio Φ j R and speed V R represent the input and output of the speed subsystem, respectively, f V ( t ) represents the unknown differentiable nonlinear system function, g V ( t ) represents the known differentiable control gain function, and d V ( t ) represents the compound disturbance.
The control objective of the HSV speed subsystem is to ensure that the output signal V can stably track the reference instruction V r e f with the prescribe tracking accuracy within the prescribe time of the designer, and the speed will not exceed the constrained set Ω V : = V R n : V < k c V , even if the actuator fault occurs, by designing the control input Φ .
Assumption 1.
The reference instruction V r e f and its derivatives of order n are smooth and bounded, that is, for any t > 0 , there is a normal number Y V such that the reference instruction V r e f satisfies V r e f ( t ) Y V < k c V , and all derivatives of the reference instruction V r e f satisfy V r e f ( i ) ( t ) Y V i , i = 1 , 2 , , n .
Using the time scale coordinate mapping method in Section 3.1, the velocity subsystem (18) is rewritten as
d V d τ = λ ( τ ) f V ( τ ) + g V ( τ ) j = 1 m 1 ρ j h 1 Φ j ( τ ) + ψ j h 1 + d V ( τ )
Remark 2.
When the original system is transformed into a new system (19), the prescribe time control problem is transformed into a general asymptotic convergence problem. As long as stable convergence can be achieved in the infinite time domain of the new system, stable convergence can be achieved in the original image in the mapping relationship, that is, within a predetermined time period.
Define speed tracking error
e V = V V r e f
Before performing the backstepping process design, first introduce a reduced order function s g V ·
s g V e V = e V e V , e V π V e V π V 2 e V 2 2 + e V , e V < π V
and a switching function
ϕ V e V = 1 , e V π V 0 , e V < π V
where π V > 0 is a constant that the designer can design in advance.
According to (21) and (22), the following conclusion can be obtained
s g V e V ϕ V e V = e V e V , e V π V 0 , e V < π V
and
ϕ V e V n = ϕ V e V
Recalling the speed subsystem (19) and (20), we can obtain
d e V d τ = λ ( τ ) f V ( τ ) + g V ( τ ) j = 1 m 1 ρ j h 1 Φ j ( τ ) + ψ j h 1 + d V ( τ ) V ˙ r e f
Choose an improved quadratic Lyapunov function
L V = 1 2 e V π V 2 ϕ V
Remark 3.
It is worth noting that the key difference between the proposed improved Lyapunov function and the traditional Lyapunov function used in similar studies is the addition of a switching function. In the traditional Lyapunov function, the constraint problem is only handled with an error constraint, but there is a drawback that it is not necessary to deal with the situation that the error itself is within the set interval. Therefore, an improved Lyapunov function is proposed with a switching function, take the following two conditions of the switching function into account: if the error is within the set interval, the function value is 0; if the error exceeds the set interval, the function value is 1. Meanwhile, because the switching function is a piecewise function, singularity may occur during the derivation process. Therefore, the proposed switching function ensures that its n-th power remains its own, which can convert discreteness into continuity and effectively avoid such a singularity problems.
Then, the time derivative of L V is expressed as
d L V d τ = e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) j = 1 m 1 ρ j h 1 Φ j ( τ ) + ψ j h 1 + λ ( τ ) f V ( τ ) + λ ( τ ) d V ( τ ) V ˙ r e f
Define an unknown nonlinear function S V ( Z V ) , where Z V = V , V ˙ r e f T R 2
S V ( Z V ) = λ ( τ ) f V ( τ ) + λ ( τ ) d V ( τ ) V ˙ r e f
Using the general form of neural network (13) to approximate the unknown function S V ( Z V ) , then
d L V d τ = e V π V ϕ V s g V e V Θ V T ψ V ( Z V ) + ε V ( Z V ) + λ ( τ ) g V ( τ ) j = 1 m 1 ρ j h 1 Φ j ( τ ) + ψ j h 1 e V π V ϕ V s g V e V θ V ρ V + λ ( τ ) g V ( τ ) j = 1 m ρ j h 1 Φ j ( τ ) + ψ j h 1
Define
θ V = Θ V T Θ V + ε ¯ V 2
ρ V = s g V e V φ V T Z V φ V Z V + l 0
where l 0 is a positive design constant.
Next, define
s 1 = inf t 0 j = 1 m 1 ρ j h 1 , ϑ 1 = 1 s 1 ζ 1 = sup t 0 j = 1 m 1 ψ j h 1
The unknown parameters ϑ 1 and ζ 1 will be estimated by designing an adaptive law.
Next, rewrite the time derivative of L V as
d L V d τ e V π V ϕ V s g V e V θ V ρ V + λ ( τ ) g V ( τ ) s 1 Φ j ( τ ) + ζ 1 + v 1 ( τ ) v 1 ( τ )
A new switching regulation function V is proposed
V = e V π V s g V e V ρ V , e V π V 0 , e V < π V
Construct an auxiliary control law as follows
v 1 ( τ ) = θ ^ V ρ V + c V s g V e V + k V e V π V 3 s g V e V + 0.275 λ ( τ ) g V ( τ ) ϵ ζ ^ 1 e V π V ϕ V s g V e V + γ 2 ρ V V + λ ( τ ) g V ( τ ) ζ ^ 1 tanh e V π V ϕ V s g V e V / ϵ
where θ ^ V is the estimate of θ V , ζ ^ 1 is the estimate of ζ 1 , and there is ζ ˜ 1 = ζ 1 ζ ^ 1 . c V , k V and γ are positive design constants.
Substituting (35) into (33) yields
d L V d τ θ ˜ V V c V e V π V ϕ V k V e V π V 4 ϕ V 0.275 λ ( τ ) g V ( τ ) ϵ ζ ^ 1 + e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) s 1 Φ j ( τ ) + λ ( τ ) g V ( τ ) ζ 1 + v 1 ( τ ) λ ( τ ) g V ( τ ) ζ ^ 1 tanh e V π V ϕ V s g V e V / ϵ γ 2 V 2
where θ ˜ V = θ V θ ^ V is the weight estimation error of the generalized neural network.
Then, construct the actual smoothing control input Φ j ( τ ) as
Φ j ( τ ) = e V π V ϕ V s g V e V ϑ ^ 1 2 v 1 2 ( τ ) λ ( τ ) g V ( τ ) e V π V 2 ϕ V 2 s g V 2 e V ϑ ^ 1 2 v 1 2 ( τ ) + σ 1 2 σ 1 e V π V ϕ V s g V e V λ ( τ ) g V ( τ )
where σ 1 is a positive design constant.
The following adaptive update laws are designed as
θ ^ ˙ V = Proj γ V , θ ^ V ( 0 ) Ω θ V
ζ ^ ˙ 1 = r e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) tanh e V π V ϕ V s g V e V / ϵ + 0.275 r λ ( τ ) g V ( τ ) ϵ
ϑ ^ ˙ 1 = Γ e V π V ϕ V s g V e V v 1 ( τ )
where ϑ ^ 1 is the estimation of ϑ 1 and there is an estimation error ϑ ˜ 1 = ϑ 1 ϑ ^ 1 . Proj · projection operator, Ω θ V is a known compact set satisfying θ V Ω θ V .
So far, the design process of adaptive tracking controller for speed subsystem in infinite time domain has been completed.

3.3. Controller Design of Height Subsystem

This section aims to design an adaptive tracking control scheme for the altitude subsystem under actuator failure, and ensure that the altitude tracking error converges to the range prescribed by the designer within the prescribe time. Firstly, the nonlinear dynamic equation of the height subsystem of HSV is described as
h ˙ ( t ) = f h ( t ) + g h ( t ) γ ( t ) + d h ( t ) γ ˙ ( t ) = f γ ( t ) + g γ ( t ) α ( t ) + d γ ( t ) α ˙ ( t ) = f α ( t ) + g α ( t ) Q ( t ) + d α ( t ) Q ˙ ( t ) = f Q ( t ) + g Q ( t ) j = 1 m 2 ρ j h 2 δ e j ( t ) + ψ j h 2 + d Q ( t )
where altitude h, track angle γ , angle of attack α , and pitch angular velocity Q represent the four state variables of the system respectively, and elevator deflection angle δ e j R and altitude h R represent the input and output of the altitude subsystem, respectively. For the convenience of description, the state variables of the system are recorded as χ η , η = 1 , 2 , 3 , 4 , where χ 1 = h , χ 2 = γ , χ 3 = α , χ 4 = Q , then f χ η ( t ) represents the unknown differentiable nonlinear system function, g χ η ( t ) represents the known differentiable control gain function, and d χ η ( t ) represents the compound disturbance.
The control objective of the HSV height subsystem is to ensure that the output signal h can stably track the reference instruction h r e f with the prescribe tracking accuracy within the prescribe time by the designer, and the state quantity χ η will not exceed the constrained set Ω χ η : = χ η R n : χ η < k c η , even if the actuator fault occurs, by designing the control input δ e .
Assumption 2.
The reference instruction h r e f and its n-order derivatives are smooth and bounded, that is, for any t > 0 , there is a normal number Y h such that the reference instruction h r e f satisfies h r e f ( t ) Y 0 < k c 1 , and the derivatives of the reference instruction h r e f satisfy h r e f ( η ) ( t ) Y η , η = 1 , 2 , 3 , 4 .
Using the time scale coordinate mapping method in Section 3.1, the height subsystem Equation (41) is rewritten as
d h d τ = λ ( τ ) f h ( τ ) + g h ( τ ) γ ( τ ) + d h ( τ ) d γ d τ = λ ( τ ) f γ ( τ ) + g γ ( τ ) α ( τ ) + d γ ( τ ) d α d τ = λ ( τ ) f α ( τ ) + g α ( τ ) Q ( τ ) + d α ( τ ) d Q d τ = λ ( τ ) f Q ( τ ) + g Q ( τ ) j = 1 m 2 ρ j h 2 δ e j ( τ ) + ψ j h 2 + d Q ( τ )
It should be mentioned that backstepping technique is used to construct an adaptive controller for nonlinear system, the recursive design procedure contains four steps. To facilitate the readers’ comprehension, the general block diagram of the proposed control scheme is given in Figure 1.
The design of adaptive control law is firstly based on the definition of tracking error
e h = h ( τ ) h r e f ( τ ) e γ = γ ( τ ) α 1 ( τ ) e α = α ( τ ) α 2 ( τ ) e Q = Q ( τ ) α 3 ( τ )
A reduced order function s g χ η · , η = 1 , 2 , 3 , 4 is introduced before performing the backstepping process design
s g χ η e χ η = e χ η e χ η , e χ η π χ η e χ η π η 2 e η 2 2 + e χ η , e χ η < π χ η
and a switching function
ϕ χ η e χ η = 1 , e χ η π χ η 0 , e χ η < π χ η
where a is π χ η > 0 constant that the designer can design in advance.
According to (21) and (22), the following conclusion can be obtained
s g χ η e χ η ϕ χ η e χ η = e χ η e χ η , e χ η π χ η 0 , e χ η < π χ η
and
ϕ χ η e χ η n = ϕ χ η e χ η
where n is positive integer.
Next, the detailed controller design process based on backstepping technology is as follows.

3.3.1. Step 1

Firstly, the intermediate control α 1 is designed to make the corresponding subsystem (track angle γ ) toward the equilibrium position.
Consider the altitude subsystem and note the tracking error e h = h ( τ ) h r e f ( τ ) of altitude h
d e h d τ = λ ( τ ) f h ( τ ) + g h ( τ ) γ ( τ ) + d h ( τ ) h ˙ r e f
Choose an improved quadratic Lyapunov function
L 1 e = 1 2 e h π h 2 ϕ h
By introducing e γ = γ ( τ ) α 1 ( τ ) , the time derivative of L 1 e is described as
d L 1 e d τ = e h π h ϕ h s g h e h λ ( τ ) g h ( τ ) e γ + α 1 + λ ( τ ) f h ( τ ) + λ ( τ ) d h ( τ ) h ˙ r e f
Define an unknown nonlinear function S 1 ( Z 1 )
S 1 ( Z 1 ) = λ ( τ ) f h ( τ ) + λ ( τ ) d h ( τ ) h ˙ r e f
Using the general form of neural network (13) to approximate the unknown function S 1 ( Z 1 ) , then
d L 1 e d τ = e h π h ϕ h s g h e h Θ 1 T ψ 1 ( Z 1 ) + ε 1 ( Z 1 ) + λ ( τ ) g h ( τ ) e γ + α 1 e h π h ϕ h s g h e h θ 1 T φ 1 Z 1 + λ ( τ ) g h ( τ ) e γ + α 1 e h π h ϕ h s g h e h θ ρ 1 + λ ( τ ) g h ( τ ) e γ + α 1
where θ i = Θ i T , ε i ( Z i ) T , φ i Z i = ψ i T Z i , 1 T .
Define
θ = max Θ i T Θ i + ε ¯ i 2 , i = 1 , 2 , 3 , 4
ρ 1 = s g h e h φ 1 T Z 1 φ 1 Z 1 + l 0
where l 0 is a positive design constant.
A new switching regulation function 1 is proposed
1 = e h π h s g h e h ρ 1 , e h π h 0 , e h < π h
Design a virtual controller α 1 ( τ ) as follows
α 1 ( τ ) = 1 λ ( τ ) g h ( τ ) θ ^ ρ 1 c 1 s g h e h k 1 e h π h 3 s g h e h 1 4 e h π h s g h e h π γ + 1 s g h e h γ 2 ρ 1 1
where θ ^ is the estimate of θ , c 1 , k 1 , and γ are positive design constants.
Substituting (56) into (52) yields
d L 1 e d τ e h π h ϕ h s g h e h θ ˜ ρ 1 c 1 e h π h ϕ h s g h 2 e h k 1 e h π h 4 ϕ h s g h 2 e h 1 4 e h π h 2 ϕ h s g h 2 e h γ 2 ρ 1 e h π h ϕ h s g h e h e h π h s g h e h ρ 1 + e h π h ϕ h s g h e h λ ( τ ) g h ( τ ) e γ π γ + 1 e h π h ϕ h s g h 2 e h 1 θ ˜ c 1 e h π h ϕ h k 1 e h π h 4 ϕ h 1 4 e h π h 2 ϕ h γ 2 1 2 + e h π h ϕ h λ ( τ ) g h ( τ ) e γ π γ 1
where θ ˜ = θ θ ^ is the weight estimation error of generalized neural network.

3.3.2. Step 2

Secondly, the intermediate control α 2 is designed to make the corresponding subsystem (attack angle α ) tend towards equilibrium position.
Consider the tracking error e γ = γ ( τ ) α 1 ( τ ) and e α = α ( τ ) α 2 ( τ ) of track angle γ ( τ )
d e γ d τ = λ ( τ ) f γ ( τ ) + g γ ( τ ) α ( τ ) + d γ ( τ ) α ˙ 1 = λ ( τ ) f γ ( τ ) + g γ ( τ ) e α + α 2 ( τ ) + d γ ( τ ) α 1 h h ˙ α 1 γ ( τ ) γ ˙ ( τ ) m = 0 2 α 1 h r e f ( m ) h r e f ( m + 1 ) α 1 θ ^ θ ^ ˙
Choose an improved quadratic Lyapunov function
L 2 e = L 1 e + 1 2 e γ π γ 2 ϕ γ
Define an unknown nonlinear function S 2 ( Z 2 )
S 2 ( Z 2 ) = λ ( τ ) f γ ( τ ) + λ ( τ ) d γ ( τ ) α 1 h h ˙ α 1 γ ( τ ) γ ˙ ( τ ) m = 0 2 α 1 h r e f ( m ) h r e f ( m + 1 )
Using the general form of neural network (13) to approximate unknown functions S η ( Z η ) , η = 1 , 2 , 3 , 4 to obtain
S η ( Z η ) = Θ η T ψ η ( Z η ) + ε η ( Z η ) = θ η T φ η Z η
Invoking (58) and (61), the time derivative of L 2 e can be obtained by (59)
d L 2 e d τ = d L 1 e d τ + e γ π γ ϕ γ s g γ e γ d e γ d τ d L 1 e d τ + e γ π γ ϕ γ s g γ e γ λ ( τ ) g γ ( τ ) e α + α 2 ( τ ) + θ ρ 2 α 1 θ ^ θ ^ ˙
where ρ 2 is designed as
ρ 2 = s g γ e γ φ 2 T Z 2 φ 2 Z 2 + l 0
Similarly, a switching adjustment function 2 is proposed as
2 = e γ π γ s g γ e γ ρ 1 + 1 , e γ π γ 1 , e γ < π γ
Design a virtual controller α 2 τ as follows
α 2 τ = 1 λ ( τ ) g γ ( τ ) θ ^ ρ 2 c 2 s g γ e γ k 2 e γ π γ 3 s g γ e γ π α + 1 s g γ e γ 1 4 e γ π γ s g γ e γ γ 2 ρ 2 2 + 1 λ ( τ ) g h ( τ ) e γ π γ 2 e γ π γ s g γ e γ 1 4 e γ π γ s g γ e γ α 1 θ ^ 2
where c 2 and k 2 are positive design constants.
Substituting (65) into (62)
d L 2 e d τ η = 1 2 c η e χ η π χ η ϕ χ η η = 1 2 k η e χ η π χ η 4 ϕ χ η + θ ˜ τ 1 1 4 e γ π γ 2 ϕ γ + e γ π γ ϕ γ λ ( τ ) g γ ( τ ) e α π α 1 + θ ˜ e γ π γ ϕ γ s g γ e γ ρ 2 γ 2 ρ 2 2 + 1 e γ π γ ϕ γ s g γ e γ 1 4 e γ π γ 2 ϕ γ α 1 θ ^ 2 e γ π γ ϕ γ s g γ e γ α 1 θ ^ θ ^ ˙ γ 2 1 2 + Y 2
where
Y 2 = 1 4 e h π h 2 ϕ h λ ( τ ) g h ( τ ) e γ π γ 2 ϕ γ + e h π h ϕ h λ ( τ ) g h ( τ ) e γ π γ 1
It is worth mentioning that the result of Y 2 0 can be obtained through analysis. If λ ( τ ) g h ( τ ) e γ π γ + 1 , there is obviously e h π h ϕ h λ ( τ ) g h ( τ ) e γ π γ 1 0 , so Y 2 0 ; if λ ( τ ) g h ( τ ) e γ > π γ + 1 , then ϕ γ = 1 , using Young’s inequality, we obtain
Y 2 1 4 e h π h 2 ϕ h λ ( τ ) g h ( τ ) e γ π γ 2 + 1 4 e h π h 2 ϕ h + λ ( τ ) g h ( τ ) e γ π γ 1 2 0
According to the definition in (64) to obtain
γ 2 ρ 2 2 + 1 e γ π γ ϕ γ s g γ e γ γ 2 1 2 = γ 2 ρ 2 e γ π γ s g γ e γ ρ 2 + 2 1 e γ π γ ϕ γ s g γ e γ + 1 2 = γ 2 ρ 2 e γ π γ s g γ e γ 2 + 2 1 ρ 2 e γ π γ s g γ e γ + 1 2 = γ 2 2 2
Substituting (67) and (69) into (66) yields
d L 2 e d τ θ ˜ 2 η = 1 2 c η e χ η π χ η ϕ χ η η = 1 2 k η e χ η π χ η 4 ϕ χ η γ 2 2 2 1 4 e γ π γ 2 ϕ γ + e γ π γ ϕ γ λ ( τ ) g γ ( τ ) e α π α 1 1 4 e γ π γ 2 ϕ γ α 1 θ ^ 2 e γ π γ ϕ γ s g γ e γ α 1 θ ^ θ ^ ˙

3.3.3. Step 3

Thirdly, the intermediate control α 3 is designed to make the corresponding (pitch angular velocity Q) subsystem tend towards equilibrium position.
Consider the tracking error e α = α ( τ ) α 2 ( τ ) of the angle of attack α ( τ )
d e α d τ = λ ( τ ) f α ( τ ) + g α ( τ ) Q ( τ ) + d α ( τ ) α ˙ 2 = λ ( τ ) f α ( τ ) + g α ( τ ) e Q + α 3 ( τ ) + d α ( τ ) η = 1 3 α 2 χ η χ ˙ η m = 0 3 α 2 h r e f ( m ) h r e f ( m + 1 ) α 2 θ ^ θ ^ ˙
Choose an improved quadratic Lyapunov function
L 3 e = L 2 e + 1 2 e α π α 2 ϕ α
Define an unknown nonlinear function S 3 ( Z 3 )
S 3 ( Z 3 ) = λ ( τ ) f α ( τ ) + λ ( τ ) d α ( τ ) η = 1 3 α 2 χ η χ ˙ η m = 0 3 α 2 h r e f ( m ) h r e f ( m + 1 )
Using (58) and (61), the time derivative of L 3 e can be obtained by (59)
d L 3 e d τ = d L 2 e d τ + e α π α ϕ α s g α e α d e α d τ d L 1 e d τ + e α π α ϕ α s g α e α λ ( τ ) g α ( τ ) e Q + α 3 ( τ ) + θ ρ 3 α 2 θ ^ θ ^ ˙
where ρ 3 is designed as
ρ 3 = s g α e α φ 3 T Z 3 φ 3 Z 3 + l 0
Similarly, a switching adjustment function 3 is proposed as
3 = e α π α s g α e α ρ 2 + 2 , e α π α 2 , e α < π α
Design a virtual controller α 3 τ as follows
α 3 τ = 1 λ ( τ ) g α ( τ ) θ ^ ρ 3 c 3 s g α e α k 3 e α π α 3 s g α e α γ 2 ρ 3 3 + 2 1 4 e α π α s g α e α λ ( τ ) g γ ( τ ) e α π α 2 e α π α s g α e α π Q + 1 s g α e α 1 4 e α π α s g α e α α 2 θ ^ 2
where c 3 and k 3 are positive design constants.
Similar to Step 2
d L 3 e d τ θ ˜ 3 η = 1 3 c η e χ η π χ η ϕ χ η η = 1 3 k η e χ η π χ η 4 ϕ χ η 1 4 e α π α 2 ϕ γ η = 2 3 e χ η π χ η ϕ χ η s g χ η e χ η α η 1 θ ^ θ ^ ˙ 1 4 η = 2 3 e χ η π χ η 2 ϕ χ η α η 1 θ ^ 2 + e α π α ϕ α λ ( τ ) g α ( τ ) e Q π Q 1 γ 2 3 2

3.3.4. Step 4

Fourthly, the stabilization of system can be achieved with the actual control input u (elevator deflection angle δ e j ) being designed.
Consider the tracking error e Q = Q ( τ ) α 3 ( τ ) of the pitch angular velocity Q
d L 3 e d τ θ ˜ 3 η = 1 3 c η e χ η π χ η ϕ χ η η = 1 3 k η e χ η π χ η 4 ϕ χ η 1 4 e α π α 2 ϕ γ η = 2 3 e χ η π χ η ϕ χ η s g χ η e χ η α η 1 θ ^ θ ^ ˙ 1 4 η = 2 3 e χ η π χ η 2 ϕ χ η α η 1 θ ^ 2 + e α π α ϕ α λ ( τ ) g α ( τ ) e Q π Q 1 γ 2 3 2
Choose an improved quadratic Lyapunov function
L 4 e = L 3 e + 1 2 e Q π Q 2 ϕ Q
Define an unknown nonlinear function S 4 ( Z 4 )
S 4 ( Z 4 ) = λ ( τ ) f Q ( τ ) + λ ( τ ) d Q ( τ ) η = 1 4 α 3 χ η χ ˙ η m = 0 4 α 3 h r e f ( m ) h r e f ( m + 1 )
And define
s 2 = inf t 0 j = 1 m 2 ρ j h 2 , ϑ 2 = 1 s 2 ζ 2 = sup t 0 j = 1 m 2 ψ j h 2
where the unknown parameters ϑ 2 and ζ 2 will be estimated by designing an adaptive law.
With the help of (79) and (81), the time derivative of L 4 e can be obtained by (80)
d L 4 e d τ = d L 3 e d τ + e Q π Q ϕ Q s g Q e Q d e Q d τ d L 3 e d τ + e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) s 2 δ e j ( τ ) + ζ 2 + θ ρ 4 α 3 θ ^ θ ^ ˙
where ρ 4 is designed as
ρ 4 = s g Q e Q φ 4 T Z 4 φ 4 Z 4 + l 0
Next, an auxiliary control law v 2 ( τ ) is constructed as follows
v 2 ( τ ) = θ ^ ρ 4 k 4 e Q π Q 3 s g Q e Q γ 2 ρ 4 4 + 3 + 0.275 λ ( τ ) g Q ( τ ) ϵ ζ ^ 2 e Q π Q ϕ Q s g Q e Q 1 4 e Q π Q s g Q e Q α 2 θ ^ 2 λ ( τ ) g α ( τ ) e Q π Q 2 e Q π Q s g Q e Q + λ ( τ ) g Q ( τ ) ζ ^ 2 tanh e Q π Q ϕ Q s g Q e Q / ϵ c 4 s g Q e Q
where c 3 and k 4 are positive design constants, ζ ^ 2 is the estimate of ζ 2 and ζ ˜ 2 = ζ 2 ζ ^ 2 .
Similarly, the switching adjustment function 4 is constructed as follows
4 = e Q π Q s g Q e Q ρ 4 + 3 , e Q π Q 3 , e Q < π Q
In combination with (85), add and subtract the term v 2 ( τ ) of (83) to obtain
d L 4 e d τ θ ˜ 4 η = 1 4 c η e χ η π χ η ϕ χ η η = 1 4 k η e χ η π χ η 4 ϕ χ η 0.275 λ ( τ ) g Q ( τ ) ϵ ζ ^ 2 1 4 η = 2 4 e χ η π χ η 2 ϕ χ η α η 1 θ ^ 2 η = 2 4 e χ η π χ η ϕ χ η s g χ η e χ η α η 1 θ ^ θ ^ ˙ + e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) ζ ^ 2 tanh e Q π Q ϕ Q s g Q e Q / ϵ + λ ( τ ) g Q ( τ ) s 2 δ e j ( τ ) + ζ 2 + v 2 ( τ ) γ 2 4 2
Next, construct the actual smoothing control input δ e j ( τ ) as
δ e j ( τ ) = e Q π Q ϕ Q s g Q e Q ϑ ^ 2 2 v 2 2 ( τ ) λ ( τ ) g Q ( τ ) e Q π Q 2 ϕ Q 2 s g Q 2 e Q ϑ ^ 2 2 v 2 2 ( τ ) + σ 2 2 σ 2 e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ )
where σ 2 is a positive design constant.
Design the following adaptive update laws are as
θ ^ ˙ = Proj γ 4 , θ ^ ( 0 ) Ω θ
ζ ^ ˙ 2 = r e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) tanh e Q π Q ϕ Q s g Q e Q / ϵ + 0.275 r λ ( τ ) g Q ( τ ) ϵ
ϑ ^ ˙ 2 = Γ e Q π Q ϕ Q s g Q e Q v 2 ( τ )
where ϑ ^ 2 is the estimation of ϑ 2 with estimation error ϑ ˜ 2 = ϑ 2 ϑ ^ 2 , Ω θ is a known compact set satisfying θ Ω θ .
So far, the design process of the adaptive tracking controller for the height subsystem in the infinite time domain has been completed.

4. Stability Analysis

4.1. Stability Analysis of Speed Subsystem

Theorem 1.
Considering the nonlinear system (18) under Assumption 1, the actual control input (103) is designed by introducing the auxiliary control law (102) and the parameter adaptive update law (104)–(106). The control scheme can ensure that: (1) All the signals of the speed subsystem are bounded; (2) The tracking error can converge to the interval defined by the designer within the prescribe fixed time T P , i.e., lim t T P e V ( t ) π V , in which the parameters T P and π V can be designed in advance; (3) The speed will not exceed the set Ω V : = V R n : V < k c V of constraints.
Proof. 
In order to analyze the stability of the velocity subsystem, the following Lyapunov function is considered
L 1 = L V + L θ + L ϑ + L ζ
where L θ = 1 2 γ θ ˜ V 2 , L ϑ = s 1 2 Γ ϑ ˜ 1 2 , L ζ = 1 2 r ζ ˜ 1 2 . □
Combining with θ ˜ ˙ V = θ ˙ V θ ^ ˙ V = θ ^ ˙ V , substitute (38) into (36) to obtain
d L 1 d τ θ ˜ V V c V e V π V ϕ V k V e V π V 4 ϕ V 0.275 λ ( τ ) g V ( τ ) ϵ ζ ^ 1 s 1 Γ ϑ ˜ 1 ϑ ^ ˙ 1 1 r ζ ˜ 1 ζ ^ ˙ 1 + e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) s 1 Φ j ( τ ) + λ ( τ ) g V ( τ ) ζ 1 + v 1 ( τ ) λ ( τ ) g V ( τ ) ζ ^ 1 tanh e V π V ϕ V s g V e V / ϵ 1 γ θ ˜ V Proj γ V
We know the property Proj V V of the inner projection operator of compact set Ω θ V , so we obtain
d L 1 d τ 1 r ζ ˜ 1 r e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) tanh e V π V ϕ V s g V e V / ϵ ζ ^ ˙ 1 + e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) s 1 Φ j ( τ ) + e V π V ϕ V s g V e V v 1 ( τ ) c V e V π V ϕ V k V e V π V 4 ϕ V 0.275 λ ( τ ) g V ( τ ) ϵ ζ ^ 1 s 1 Γ ϑ ˜ 1 ϑ ^ ˙ 1 + λ ( τ ) g V ( τ ) ζ 1 e V π V ϕ V s g V e V e V π V ϕ V s g V e V × tanh e V π V ϕ V s g V e V / ϵ
Recalling (37) and Lemma 2 to obtain
e V π V ϕ V s g V e V λ ( τ ) g V ( τ ) s 1 Φ j ( τ ) s 1 e V π V 2 ϕ V 2 s g V 2 e V ϑ ^ 1 2 v 1 2 ( τ ) e V π V 2 ϕ V 2 s g V 2 e V ϑ ^ 1 2 v 1 2 ( τ ) + σ 1 2 s 1 σ 1 s 1 e V π V ϕ V s g V e V ϑ ^ 1 v 1 ( τ ) s 1 e V π V ϕ V s g V e V ϑ ^ 1 v 1 ( τ )
With inequality 0 q q tanh q / ϵ 0.2785 ϵ , it holds that
λ ( τ ) g V ( τ ) ζ 1 e V π V ϕ V s g V e V e V π V ϕ V s g V e V × tanh e V π V ϕ V s g V e V / ϵ 0.275 λ ( τ ) g V ( τ ) ζ 1 ϵ
Substituting (39), (40), (95), and (96) into (94), we obtain
d L 1 d τ c V e V π V ϕ V k V e V π V 4 ϕ V 0.275 λ ( τ ) g V ( τ ) ϵ ζ ^ 1 s 1 e V π V ϕ V s g V e V ϑ ^ 1 v 1 ( τ ) + e V π V ϕ V s g V e V v 1 ( τ ) + 0.275 λ ( τ ) g V ( τ ) ϵ ζ 1 0.275 λ ( τ ) g V ( τ ) ϵ ζ ˜ 1 s 1 ϑ ˜ 1 e V π V ϕ V s g V e V v 1 ( τ ) c V e V π V 2 ϕ V 1 / 2 k V e V π V 2 ϕ V 2
According to the result of (97), it is obtained that L 1 has the property of non-increasing, so the boundedness of e V , θ ˜ V , ϑ ˜ 1 and ζ ˜ 1 is guaranteed. Since θ ˜ V = θ V θ ^ V , and θ V are bounded constants, θ ^ V L . Because V r e f and its derivatives are bounded, V L . Through similar analysis, all closed-loop signals are bounded.
It is now proven that the proposed control method can allow the tracking error to converge to a range that can be defined by the designer within a predetermined time.
With the help of (97), the following inequality holds
k V e V π V 2 ϕ V 2 d L 1 d τ c V e V π V 2 ϕ V 1 / 2 d L 1 d τ
Integrate both sides of (98) to obtain
0 e V π V 2 ϕ V 2 d τ 1 k V L 1 + L 1 0 1 k V L 1 0
where indicates e V π V 2 ϕ V L 2 .
By applying the Barbalat lemma yields
lim τ e V τ π V 2 ϕ V = 0
Therefore, the tracking error e V τ asymptotically approaches a prescribe interval π V in the infinite domain time τ .
In turn, the adaptive controller in the original nonlinear velocity subsystem (18) is designed. According to the inverse transformation of time scale, τ in the control law in Section 3.2 is replaced by t, so the prescribe time control is realized through the mapping function from the perspective of time scale.
The auxiliary control law v 1 ( t ) and the actual control input Φ j ( t ) are designed as follows:
v 1 ( t ) = θ ^ V ρ V + c V s g V e V + k V e V π V 3 s g V e V + 0.275 λ ( τ ) g V ( τ ( t ) ) ϵ ζ ^ 1 e V π V ϕ V s g V e V + γ 2 ρ V τ V + λ ( τ ) g V ( τ ( t ) ) ζ ^ 1 tanh e V π V ϕ V s g V e V / ϵ
Φ j ( t ) = e V π V ϕ V s g V e V ϑ ^ 1 2 v 1 2 ( τ ( t ) ) λ ( τ ) g V ( τ ( t ) ) e V π V 2 ϕ V 2 s g V 2 e V ϑ ^ 1 2 v 1 2 ( τ ( t ) ) + σ 1 2 σ 1 e V π V ϕ V s g V e V λ ( τ ) g V ( τ ( t ) )
The adaptive updating law of parameters is
θ ^ ˙ V t = Proj γ τ V τ ( t ) , θ ^ V ( 0 ) Ω θ V
ζ ^ ˙ 1 t = r e V π V ϕ V s g V e V λ ( τ ) g V τ ( t ) tanh e V π V ϕ V s g V e V / ϵ + 0.275 r λ ( τ ) g V τ ( t ) ϵ
ϑ ^ ˙ 1 t = Γ e V π V ϕ V s g V e V v 1 τ ( t )
Similar to the above proof process, consider the following Lyapunov function
V 1 t = V V t + V θ t + V ϑ t + V ζ t
where V θ t = 1 2 γ θ ˜ V 2 t , V ϑ t = s 1 2 Γ ϑ ˜ 1 2 t , V ζ t = 1 2 r ζ ˜ 1 2 t .
Along the same lines, we can draw a conclusion
d V 1 d t c V e V t π V 2 ϕ V 1 / 2 k V e V t π V 2 ϕ V 2
Because the time scale coordinate mapping method is adopted, the prescribe time t [ 0 , T P ) of the finite field is mapped to the time τ [ 0 , + ) of the infinite field, which is obtained from (98)–(100)
lim t T P e V ( t ) π V
Therefore, the actual speed tracking error e V ( t ) can converge to the interval π V defined by the designer within the prescribe time T P . Next, it is proved that the speed will not exceed the set Ω V : = V R n : V < k c V of constraints. Because lim t T P e V ( t ) π V , and π V can be designed in advance, if π V = k b V k c V Y V is defined, it can be inferred that there is V t e V t + V r e f t π V + Y V k c V for t T P , so the conclusion is proved.
So far, the proof of Theorem 1 has been completed.

4.2. Stability Analysis of Altitude Subsystem

Theorem 2.
Considering the nonlinear system (41) under Assumption 2, the actual control input (122) is designed by introducing the auxiliary control law (121) and the parameter adaptive update law(123)–(126). The control scheme can ensure that: (1) all signals of the altitude subsystem are bounded; (2) The tracking error can converge to the interval defined by the designer within the prescribe fixed time T P , that is, lim t T P e χ η ( t ) π χ η , η = 1 , 2 , 3 , 4 , in which the parameters T P and π χ η can be designed in advance; (3) All state quantities χ η do not exceed the set Ω χ η : = χ η R n : χ η < k c η of constraints.
Proof. 
In order to analyze the stability of the velocity subsystem, the following Lyapunov function is considered
L 2 = L 4 e + L 2 θ + L 2 ϑ + L 2 ζ
where L 2 θ = 1 2 γ θ ˜ 2 , L 2 ϑ = s 2 2 Γ ϑ ˜ 2 2 , L 2 ζ = 1 2 r ζ ˜ 2 2 . □
Combining with θ ˜ ˙ = θ ˙ θ ^ ˙ = θ ^ ˙ , substituting (89) and (87) into (109) yields
d L 2 d τ θ ˜ γ γ 4 Proj γ 4 η = 1 4 c η e χ η π χ η ϕ χ η η = 1 4 k η e χ η π χ η 4 ϕ χ η s 2 Γ ϑ ˜ 2 ϑ ^ ˙ 2 1 r ζ ˜ 2 ζ ^ ˙ 2 + + e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) s 2 δ e j ( τ ) + ζ 2 + v 2 ( τ ) λ ( τ ) g Q ( τ ) ζ ^ 2 tanh e Q π Q ϕ Q s g Q e Q / ϵ 0.275 λ ( τ ) g Q ( τ ) ϵ ζ ^ 2
where
= η = 2 4 e χ η π χ η ϕ χ η s g χ η e χ η α η 1 θ ^ Proj γ 4 γ 2 4 2 1 4 η = 2 4 e χ η π χ η 2 ϕ χ η α η 1 θ ^ 2
We know the property Proj 4 2 4 2 and Proj 4 4 of the inner projection operator of compact set Ω θ , so we obtain
η = 2 4 e χ η π χ η ϕ χ η s g χ η e χ η α η 1 θ ^ Proj γ 4 1 4 η = 2 4 e χ η π χ η 2 ϕ χ η α η 1 θ ^ 2 + γ 2 4 2
According to (112), 0 is obtained. Therefore, d L 2 / d τ can be further expressed as
d L 2 d τ e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) s 2 δ e j ( τ ) η = 1 4 c η e χ η π χ η ϕ χ η s 2 Γ ϑ ˜ 2 ϑ ^ ˙ 2 + e Q π Q ϕ Q s g Q e Q v 2 ( τ ) + 1 r ζ ˜ 2 r e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) × tanh e Q π Q ϕ Q s g Q e Q / ϵ ζ ^ ˙ 2 0.275 λ ( τ ) g Q ( τ ) ϵ ζ ^ 2 + λ ( τ ) g Q ( τ ) ζ 2 e Q π Q ϕ Q s g Q e Q e Q π Q ϕ Q s g Q e Q × tanh e Q π Q ϕ Q s g Q e Q / ϵ η = 1 4 k η e χ η π χ η 4 ϕ χ η
According to (88) and Lemma 2, we can obtain
e Q π Q ϕ Q s g Q e Q λ ( τ ) g Q ( τ ) s 2 δ e j ( τ ) s 2 e Q π Q 2 ϕ Q 2 s g Q 2 e Q ϑ ^ 2 2 v 2 2 ( τ ) e Q π Q 2 ϕ Q 2 s g Q 2 e Q ϑ ^ 2 2 v 2 2 ( τ ) + σ 2 2 s 2 σ 2 s 2 e Q π Q ϕ Q s g Q e Q ϑ ^ 2 v 2 ( τ )
With the inequality 0 q q tanh q ϵ 0.2785 ϵ yields
λ ( τ ) g Q ( τ ) ζ 2 e Q π Q ϕ Q s g Q e Q e Q π Q ϕ Q s g Q e Q × tanh e Q π Q ϕ Q s g Q e Q / ϵ 0.275 λ ( τ ) g Q ( τ ) ζ 2 ϵ
Similarly to the previous section, substitute (90), (91), (114), and (115) into (113) to obtain
d L 2 d τ η = 1 4 c η e χ η π χ η 2 ϕ χ η 1 / 2 η = 1 4 k η e χ η π χ η 2 ϕ χ η 2
According to the result of (116), it is obtained that L 2 has the property of non-increasing, so the boundedness of e χ η , θ ˜ , ϑ ˜ 2 and ζ ˜ 2 is guaranteed. Since θ ˜ = θ θ ^ , and θ are bounded constants, θ ^ L . Because h r e f and its derivatives are bounded, h L . Since α 1 is a function of bounded signals, there is α 1 L , so γ = e γ + α 1 L . Through similar analysis, the boundedness of α and Q can be deduced in turn, so all closed-loop signals are bounded.
It is now proven that the proposed control method can allow the height tracking error to converge to a range that can be defined by the designer within a predetermined time.
With the help of (116), the following inequality holds
η = 1 4 k η e χ η π χ η 2 ϕ χ η 2 d L 2 d τ η = 1 4 c η e χ η π χ η 2 ϕ χ η 1 / 2 d L 2 d τ
Take η = 1 and integrate both sides of (117) to obtain
0 e h π h 2 ϕ h 2 d τ 1 k 1 L 2 + L 2 0 1 k 1 L 2 0
which indicates that e h π h 2 ϕ h L 2 .
Applying the Barbalat lemma yields
lim τ e h τ π h 2 ϕ h = 0
Therefore, the tracking error e h τ asymptotically approaches a prescribe interval π h within the time τ of the infinite domain.
In turn, the adaptive controller in the original nonlinear height subsystem (41) is designed. According to the inverse transformation of time scale, τ in the control law in Section 3.3 is replaced by t, so the prescribe time control is realized through the mapping function from the perspective of time scale. Similarly, the auxiliary control law v 2 ( t ) and the actual control input δ e j ( t ) are designed as follows:
v 2 ( t ) = θ ^ ρ 4 c 4 s g Q e Q k 4 e Q π Q 3 s g Q e Q + 0.275 λ ( τ ( t ) ) g Q ( τ ( t ) ) ϵ ζ ^ 2 e Q π Q ϕ Q s g Q e Q 1 4 e Q π Q s g Q e Q α 2 θ ^ 2 λ ( τ ( t ) ) g α ( τ ( t ) ) e Q π Q 2 e Q π Q s g Q e Q + λ ( τ ( t ) ) g Q ( τ ( t ) ) ζ ^ 2 tanh e Q π Q ϕ Q s g Q e Q / ϵ γ 2 ρ 4 4 + 3
δ e j ( t ) = e Q π Q ϕ Q s g Q e Q ϑ ^ 2 2 v 2 2 ( τ ( t ) ) λ ( τ ( t ) ) g Q ( τ ( t ) ) e Q π Q 2 ϕ Q 2 s g Q 2 e Q ϑ ^ 2 2 v 2 2 ( τ ( t ) ) + σ 2 2 σ 2 e Q π Q ϕ Q s g Q e Q λ ( τ ( t ) ) g Q ( τ ( t ) )
where the parameter adaptive updating law is
θ ^ ˙ t = Proj γ 4 ( t ) , θ ^ ( 0 ) Ω θ
ζ ^ ˙ 2 t = r e Q π Q ϕ Q s g Q e Q λ ( τ ( t ) ) g Q ( τ ( t ) ) tanh e Q π Q ϕ Q s g Q e Q / ϵ + 0.275 r λ ( τ ( t ) ) g Q ( τ ( t ) ) ϵ
ϑ ^ ˙ 2 t = Γ e Q π Q ϕ Q s g Q e Q v 2 ( τ ( t ) )
Similarly to the former proof process, consider the following Lyapunov function
V 2 t = V 4 e t + V 2 θ t + V 2 ϑ t + V 2 ζ t
where V 2 θ t = 1 2 γ θ ˜ 2 t , V 2 ϑ t = s 2 2 Γ ϑ ˜ 2 2 t , V 2 ζ t = 1 2 r ζ ˜ 2 2 t .
Along the same lines, we can draw a conclusion
d V 2 d t η = 1 4 c η e χ η π χ η 2 ϕ χ η 1 / 2 η = 1 4 k η e χ η π χ η 2 ϕ χ η 2
Due to the use of the time scale coordinate mapping method, the prescribe time t [ 0 , T P ) of the finite field is mapped to the time τ [ 0 , + ) of the infinite field, as obtained from Equations (117)–(119)
lim t T P e h ( t ) π h
Therefore, the actual height error e h ( t ) can converge to the interval π h defined by the designer within a prescribe time T P .
Next, it is proven that all state variables χ η will not exceed the set of constraints Ω χ η : = χ η R n : χ η < k c η .
Since lim t T P e h ( t ) π h and π h can be designed in advance, if π h = k b 1 k c 1 Y 0 is defined, it can be inferred that t T P . Due to the boundedness of the Lyapunov function V 2 t , lim t T P e χ η ( t ) π χ η can be obtained. Based on the conclusion of α 1 L , there must be a normal number α ¯ 1 such that α 1 α ¯ 1 , and if π γ = k b 2 k c 2 α ¯ 1 is defined, there is γ t e γ t + α 1 π γ + α ¯ 1 k c 2 . Using the similar approach, predefining the value of π χ η by the designer, the conclusion of χ η < k c η , η = 1 , 2 , 3 , 4 can be fully guaranteed, that is, none of the state variables χ η will exceed the set of constraints.
At this point, the proof of Theorem 2 was completed.

5. Simulation Results and Analysis

In this section, the simulation results of hypersonic vehicle are given to verify the effectiveness of the proposed controller.
For the longitudinal dynamic model, the initial state is set as velocity V 0 = 7700 ft/s, altitude h 0 = 88,000 ft, initial velocity error e V 0 = 1 ft/s, altitude error e h 0 = 2 ft/s, track angle γ 0 = 0 rad, angle of attack α 0 = 0.0284926 rad, pitch angular velocity Q 0 = 0 rad/s, and elastic state quantities η 1 0 = 0.97008 , η 2 0 = 0.796696 and η ˙ 1 0 = η ˙ 2 0 = 0 . The reference trajectories V r e f and h r e f are generated by a second-order filter V r e f s V c s = h r e f s h c s = 0.03 2 s 2 + 2 × 0.95 × 0.03 s + 0.03 2 , where V c = 1000 ft/s and h c = 3000 ft. The parameters of the controller are set to c V = 2 , k V = 5 , γ = 1 , σ 1 = 0.1 ; c 4 = 2 , k h = 55 , σ 2 = 0.1 .
The set scene requires that the tracking error converge to the prescribe range e V ( t ) 0.02 at the prescribe time T P = 20 s, and the height tracking error converge to the prescribe range e h ( t ) 0.1 at the prescribe time T P = 25 s. The analog actuator starts to fail in the 100th second, and the failure parameter is set to ρ j h 1 = ρ j h 2 = 0.8 , ψ j h 1 = 0.1 , ψ j h 2 = 2 / 57.3 rad. In order to prove the superiority of the proposed method, it is compared with the traditional control scheme, where PCs represents the proposed control scheme and CCS represents the traditional control scheme.
The simulation results are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. From Figure 2 and Figure 3, it can be seen that, even considering the actuator fault, the speed and height tracking errors converge to the prescribe range within the prescribe time. In addition to the prescribe time tracking performance, it can also be verified from Figure 2 and Figure 3 that, even if an unknown actuator fault occurs, the speed and height have strong robustness to the tracking of its instructions, will not change greatly and the tracking error can quickly converge to the range prescribed by the designer. Figure 4 shows the three rigid body state quantities γ , α , and Q of HSV. The simulation results show that the state signals are bounded. Figure 5 and Figure 6, respectively, depict the corresponding curves of the adaptive parameters of the HSV speed subsystem and height subsystem. It can be seen that θ ^ , θ ^ V , ζ ^ 1 , ζ ^ 2 , ϑ ^ 1 and ϑ ^ 2 are bounded in the process of system control. Figure 7 describes the control input signals Φ and δ e of the designed HSV and the flexible state variables η i and i = 1 , 2 , it can be concluded that the signals are bounded.

6. Discussion

Compared with the traditional control scheme, the proposed control scheme can meet the different tracking performance requirements, has faster convergence speed and higher convergence accuracy, and better solves the adaptive prescribe time tracking problem with actuator failure. In addition, it is noted that there is no constraint on the tracking error at the beginning of the system performance. This means that, compared with the existing error constraint implementation schemes, the proposed control algorithm does not need to know the initial tracking conditions in advance, and can also achieve the prescribe time tracking performance. Additionally, the selection of actuator fault parameters is based on an appropriate value. If the fault parameters are too large, greater control energy is required and even system stability cannot be ensured. If the fault parameters are too small, the impact of the fault can be ignored, which cannot reflect the progressiveness of the proposed method. Due to the design of two adaptive parameters to estimate unknown faults, the parameters can be adjusted online to achieve adaptive fault-tolerant control when faults occur, thereby enhancing the robustness of the actual flight control system. According to the tracking curve of the speed subsystem, it can be seen that the analog actuator starts to fail in the 100th second. At this time, traditional control methods cannot guarantee that the tracking error converges to the prescribed range of 0.02, and the maximum tracking error even exceeds 0.25 with a chattering phenomenon that occurs. However, the control method proposed in the article can ensure that the tracking error is within the required range.

7. Conclusions

In this paper, the adaptive fault-tolerant tracking control of hypersonic vehicle with prescribe time and tracking error is studied. Firstly, a time scale coordinate mapping function is introduced and the controller is designed in the transformed infinite time domain. An improved Lyapunov function and an improved tuning function are constructed. When the initial tracking conditions are completely unknown, the tracking error converges in the range that can be defined by the designer combined with the Barbalat lemma. At the same time, the online parameter estimation is designed to ensure the adaptive fault-tolerant control of the system. On this basis, it is ensured that the actual tracking errors of speed, altitude, track angle, angle of attack, and pitch rate can converge to the range that can be defined in advance by the designer within the prescribe time. Simulation results verify the effectiveness of the method. The control scheme effectively solves the limitation of the unknown prior knowledge of the initial error, and meets the specific needs of the actual scene. In future research directions, as the system is based on a universal prescribe time architecture, this control strategy is also applicable to other vehicle models, including wing flying aircraft. In addition, this paper focuses on the control of a single aircraft, further research will expand it to the formation control of multiple aircraft.

Author Contributions

Methodology, F.G.; Conceptualization and funding acquisition, W.Z.; data analysis, M.L.; writing, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by China Postdoctoral Science Foundation (Grant no. 2023M744292).

Data Availability Statement

The data are available from the corresponding author on reasonable request.

Acknowledgments

The authors would like to greatly appreciate the support of National Key Laboratory of Unmanned Aerial Vehicle Technology and The Youth Innovation Team of Shaanxi University. Besides, thanks for the editor and all anonymous reviewers for their comments, which help to improve the quality of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, J.; Pu, Z.; Chang, Y.; Ding, S.; Yi, J. Appointed-Time Control for Flexible Hypersonic Vehicles with Conditional Disturbance Negation. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 6327–6345. [Google Scholar] [CrossRef]
  2. Wang, L.; Qi, R.; Wen, L.; Jiang, B. Adaptive Multiple-Model-Based Fault-Tolerant Control for Non-minimum Phase Hypersonic Vehicles with Input Saturations and Error Constraints. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 519–540. [Google Scholar] [CrossRef]
  3. Li, B.-Q.; Du, M.-Z.; Zhang, L.-H.; Zong, M.-Z. A comprehensive RFD-FTC-DCA system for hypersonic vehicles with saturation constraints. Int. J. Control 2024, 97, 871–883. [Google Scholar] [CrossRef]
  4. Lv, M.; De Schutter, B.; Baldi, S. Nonrecursive Control for Formation-Containment of HFV Swarms with Dynamic Event-Triggered Communication. IEEE Trans. Ind. Inform. 2023, 19, 3188–3197. [Google Scholar] [CrossRef]
  5. Guo, J.; Wang, J.; Guo, Z.; Su, Y. Attitude optimization control of hypersonic flight vehicle considering partially unknown control direction. Trans. Inst. Meas. Control 2023, 45, 1337–1350. [Google Scholar] [CrossRef]
  6. Lv, M.; De Schutter, B.; Wang, Y.; Shen, D. Fuzzy Adaptive Zero-Error-Constrained Tracking Control for HFVs in the Presence of Multiple Unknown Control Directions. IEEE Trans. Cybern. 2023, 53, 2779–2790. [Google Scholar] [CrossRef] [PubMed]
  7. Mu, C.; Ni, Z.; Sun, C.; He, H. Air-Breathing Hypersonic Vehicle Tracking Control Based on Adaptive Dynamic Programming. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 584–598. [Google Scholar] [CrossRef] [PubMed]
  8. Basin, M.V.; Yu, P.; Shtessel, Y.B. Hypersonic Missile Adaptive Sliding Mode Control Using Finite- and Fixed-Time Observers. IEEE Trans. Ind. Electron. 2018, 65, 930–941. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Shou, Y.; Zhang, P.; Han, W. Sliding mode based fault-tolerant control of hypersonic reentry vehicle using composite learning. Neurocomputing 2022, 484, 142–148. [Google Scholar] [CrossRef]
  10. Hu, X.; Guo, C.; Hu, C.; He, B. Sliding mode learning control for T-S fuzzy system and an application to hypersonic flight vehicle. Asian J. Control 2023, 25, 407–417. [Google Scholar] [CrossRef]
  11. Zhao, D.; Jiang, B.; Yang, H. Backstepping-Based Decentralized Fault-Tolerant Control of Hypersonic Vehicles in PDE-ODE Form. IEEE Trans. Autom. Control 2022, 67, 1210–1225. [Google Scholar] [CrossRef]
  12. Bu, X.; Wu, X.; Zhang, R.; Ma, Z.; Huang, J. Tracking differentiator design for the robust backstepping control of a flexible air-breathing hypersonic vehicle. J. Frankl. Inst.-Eng. Appl. Math. 2015, 352, 1739–1765. [Google Scholar] [CrossRef]
  13. Yu, X.; Li, P.; Zhang, Y. The Design of Fixed-Time Observer and Finite-Time Fault-Tolerant Control for Hypersonic Gliding Vehicles. IEEE Trans. Ind. Electron. 2018, 65, 4135–4144. [Google Scholar] [CrossRef]
  14. Chao, D.; Qi, R.; Jiang, B. Adaptive fault-tolerant attitude control for hypersonic reentry vehicle subject to complex uncertainties. J. Frankl. Inst.-Eng. Appl. Math. 2022, 359, 5458–5487. [Google Scholar] [CrossRef]
  15. Han, T.; Hu, Q.; Shin, H.-S.; Tsourdos, A.; Xin, M. Incremental Twisting Fault Tolerant Control for Hypersonic Vehicles with Partial Model Knowledge. IEEE Trans. Ind. Inform. 2022, 18, 1050–1060. [Google Scholar] [CrossRef]
  16. Wang, L.; Qi, R.; Jiang, B. Adaptive actuator fault-tolerant control for non-minimum phase air-breathing hypersonic vehicle model. ISA Trans. 2022, 126, 47–64. [Google Scholar] [CrossRef] [PubMed]
  17. Dai, P.; Feng, D.; Feng, W.; Cui, J.; Zhang, L. Entry trajectory optimization for hypersonic vehicles based on convex programming and neural network. Aerosp. Sci. Technol. 2023, 137, 108259. [Google Scholar] [CrossRef]
  18. Xu, B.; Yang, C.; Pan, Y. Global Neural Dynamic Surface Tracking Control of Strict-Feedback Systems with Application to Hypersonic Flight Vehicle. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2563–2575. [Google Scholar] [CrossRef]
  19. Bu, X.W.; Xiao, Y.; Lei, H.M. An Adaptive Critic Design-Based Fuzzy Neural Controller for Hypersonic Vehicles: Predefined Behavioral Nonaffine Control. IEEE-ASME Trans. Mechatron. 2019, 24, 1871–1881. [Google Scholar] [CrossRef]
  20. Wang, P.F.; Wang, J.; Bu, X.W.; Jia, Y.J. Adaptive fuzzy tracking control for a constrained flexible air-breathing hypersonic vehicle based on actuator compensation. Int. J. Adv. Robot. Syst. 2016, 13, 1–11. [Google Scholar] [CrossRef]
  21. Hu, X.X.; Xu, B.; Hu, C.H. Robust Adaptive Fuzzy Control for HFV with Parameter Uncertainty and Unmodeled Dynamics. IEEE Trans. Ind. Electron. 2018, 65, 8851–8860. [Google Scholar] [CrossRef]
  22. Chen, H.; Wang, P.; Tang, G. Fuzzy Disturbance Observer-Based Fixed-Time Sliding Mode Control for Hypersonic Morphing Vehicles with Uncertainties. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 3521–3530. [Google Scholar] [CrossRef]
  23. Bu, X.; Qi, Q. Fuzzy Optimal Tracking Control of Hypersonic Flight Vehicles via Single-Network Adaptive Critic Design. IEEE Trans. Fuzzy Syst. 2022, 30, 270–278. [Google Scholar] [CrossRef]
  24. Zhang, H.; Wang, P.; Tang, G.; Bao, W. Fuzzy disturbance observer-based dynamic sliding mode control for hypersonic morphing vehicles. Aerosp. Sci. Technol. 2023, 142, 108633. [Google Scholar] [CrossRef]
  25. Xu, B.; Wang, D.; Zhang, Y.; Shi, Z. DOB-Based Neural Control of Flexible Hypersonic Flight Vehicle Considering Wind Effects. IEEE Trans. Ind. Electron. 2017, 64, 8676–8685. [Google Scholar] [CrossRef]
  26. Dong, M.; Xu, X.; Xie, F. Constrained Integrated Guidance and Control Scheme for Strap-Down Hypersonic Flight Vehicles with Partial Measurement and Unmatched Uncertainties. Aerospace 2022, 9, 840. [Google Scholar] [CrossRef]
  27. Ding, Y.; Yue, X.; Liu, C.; Dai, H.; Chen, G. Finite-time controller design with adaptive fixed-time anti-saturation compensator for hypersonic vehicle. ISA Trans. 2022, 122, 96–113. [Google Scholar] [CrossRef] [PubMed]
  28. Guo, Y.; Xu, B. Finite-Time Deterministic Learning Command Filtered Control for Hypersonic Flight Vehicle. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 4214–4225. [Google Scholar] [CrossRef]
  29. Yu, X.; Li, P.; Zhang, Y. Fixed-Time Actuator Fault Accommodation Applied to Hypersonic Gliding Vehicles. IEEE Trans. Autom. Sci. Eng. 2021, 18, 1429–1440. [Google Scholar] [CrossRef]
  30. Lv, M.; Li, Y.; Wan, L.; Dai, J.; Chang, J. Fast Nonsingular Fixed-Time Fuzzy Fault-Tolerant Control for HFVs with Guaranteed Time-Varying Flight State Constraints. IEEE Trans. Fuzzy Syst. 2022, 30, 4555–4567. [Google Scholar] [CrossRef]
  31. Dong, Z.; Li, Y.; Lv, M. Adaptive nonsingular fixed-time control for hypersonic flight vehicle considering angle of attack constraints. Int. J. Robust Nonlinear Control 2023, 33, 6754–6777. [Google Scholar] [CrossRef]
  32. Zhang, H.; Wang, P.; Tang, G.; Bao, W. Fixed-time sliding mode control for hypersonic morphing vehicles via event-triggering mechanism. Aerosp. Sci. Technol. 2023, 140, 108458. [Google Scholar] [CrossRef]
  33. Bolender, M.A.; Doman, D.B. Nonlinear Longitudinal Dynamical Model of an Air-Breathing Hypersonic Vehicle. J. Spacecr. Rocket. 2007, 44, 374–387. [Google Scholar] [CrossRef]
  34. Zhang, W.; Dong, W.; Lv, M.; Liu, Z.; Zhou, Y.; Feng, H. Barrier Lypunov functions-based nonsingular fixed-time switching control for strict-feedback nonlinear dynamics with full state constraints. Int. J. Robust Nonlinear Control 2021, 31, 7862–7885. [Google Scholar] [CrossRef]
Figure 1. Block diagram of the height subsystem.
Figure 1. Block diagram of the height subsystem.
Drones 08 00295 g001
Figure 2. Speed tracking results with the actuator fault of HSV.
Figure 2. Speed tracking results with the actuator fault of HSV.
Drones 08 00295 g002
Figure 3. Height tracking results with the actuator fault of HSV.
Figure 3. Height tracking results with the actuator fault of HSV.
Drones 08 00295 g003
Figure 4. Response curve of HSV in rigid state γ , α , and Q.
Figure 4. Response curve of HSV in rigid state γ , α , and Q.
Drones 08 00295 g004
Figure 5. Adaptive parameters θ ^ V , ϑ ^ 1 , and ζ ^ 1 of HSV speed subsystem.
Figure 5. Adaptive parameters θ ^ V , ϑ ^ 1 , and ζ ^ 1 of HSV speed subsystem.
Drones 08 00295 g005
Figure 6. Adaptive parameters θ ^ , ϑ ^ 2 , and ζ ^ 2 of HSV height subsystem.
Figure 6. Adaptive parameters θ ^ , ϑ ^ 2 , and ζ ^ 2 of HSV height subsystem.
Drones 08 00295 g006
Figure 7. Actual fuel equivalent ratio Φ , elevator deflection angle δ e , and flexibility state quantity η i of HSV.
Figure 7. Actual fuel equivalent ratio Φ , elevator deflection angle δ e , and flexibility state quantity η i of HSV.
Drones 08 00295 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, F.; Zhang, W.; Lv, M.; Zhang, R. Fault-Tolerant Tracking Control of Hypersonic Vehicle Based on a Universal Prescribe Time Architecture. Drones 2024, 8, 295. https://doi.org/10.3390/drones8070295

AMA Style

Guo F, Zhang W, Lv M, Zhang R. Fault-Tolerant Tracking Control of Hypersonic Vehicle Based on a Universal Prescribe Time Architecture. Drones. 2024; 8(7):295. https://doi.org/10.3390/drones8070295

Chicago/Turabian Style

Guo, Fangyue, Wenqian Zhang, Maolong Lv, and Ruiqi Zhang. 2024. "Fault-Tolerant Tracking Control of Hypersonic Vehicle Based on a Universal Prescribe Time Architecture" Drones 8, no. 7: 295. https://doi.org/10.3390/drones8070295

APA Style

Guo, F., Zhang, W., Lv, M., & Zhang, R. (2024). Fault-Tolerant Tracking Control of Hypersonic Vehicle Based on a Universal Prescribe Time Architecture. Drones, 8(7), 295. https://doi.org/10.3390/drones8070295

Article Metrics

Back to TopTop