Abstract

This paper presents a noncertainty equivalent adaptive backstepping control scheme for advanced fighter attitude tracking, in which unsteady effects, parameter uncertainties, and input constraints are all considered which increase the design difficulty to a large extent. Based on unsteady attitude dynamics and the noncertainty equivalent principle, a new observer is first developed to reconstruct the immeasurable and time-varying unsteady states. Afterwards, the unsteady aerodynamics is compensated in the backstepping controller where the command filter is introduced to impose physical constraints on actuators. In order to further enhance the robustness, the noncertainty equivalent adaptive approach is again used to estimate the uncertain constant parameters. Moreover, stability of the closed-loop system that includes the state observer, parameter estimator, and backstepping controller is proven by the Lyapunov theorem in a unified architecture. Finally, simulation results show that performance of the deterministic control system can be captured when attractive manifolds are achieved. The effectiveness and robustness of the proposed control scheme are verified by the Herbst maneuver.

1. Introduction

Flight envelopes are substantially extended by the fourth generation fighter which directly leads to urgent demands for higher performance control laws. Since full envelope dynamics of aircraft is nonlinear and involves wide variations in aerodynamic parameters, traditional linear control methods are incapable of addressing such challenges [1]. Therefore, to ensure adequate stability and tracking performance in extreme flight regimes, a number of nonlinear control techniques have been extensively investigated, such as dynamic inversion [2, 3], fuzzy logic [4], neural network [5], backstepping [6], and sliding mode control [7]. In particular, Lyapunov-based backstepping control is among the most widely studied of these methods. Due to the cascaded structure of aircraft dynamics, many efforts have been made to develop a flight control system via the combination of backstepping theory and other control technologies, such as disturbance observer [8], radial basis function neural network [9], and adaptive control method [10, 11].

Generally, adaptive backstepping provides a systematic approach to solve the tracking or regulation problems of uncertain nonlinear systems. However, it must be noted that the adaptive backstepping control methods mentioned above have not considered input constraints which cannot be ignored in practice due to the existence of physical limits on actuators [12]. When aircraft maneuvers at a high angle of attack (AOA), drastic decreases of control authority often causes aerodynamic surface saturation which may severely degrade control performance by giving undesirable inaccuracy or leading to instability [13]. In some applications, this problem is crucial, especially in combination with online approximation-based control, which tends to be aggressive in seeking desired tracking performance [14]. For the purpose of circumventing this dilemma, it is necessary to take the input constraints into account at the level of the control design. In constraint adaptive control, the key issue is how to analyze the constraint effects on a closed-loop system. Toward this end, an auxiliary system of which states are used for feedback control is proposed in [15, 16]. Recently, it was incorporated into backstepping architecture to solve trajectory tracking problems for reentry vehicles and air-breathing hypersonic vehicles in the presence of actuator constraints [1719]. In addition, the command filter adaptive backstepping control proposed in [14, 20] is also a wise choice to tackle input constraints, in which the command filter is used to calculate the derivatives of virtual controls and impose constraints on states and inputs. This control scheme is presented for the F-16/MATV fighter in [21], where tracking errors in adaptive laws are replaced by compensated errors; hence, the adaption process can be isolated, i.e., not affected by saturation of virtual controls or actuators.

Besides model uncertainties and input constraints, unsteady effects are also potential issues that should be considered for an advanced fighter. Unsteady aerodynamic states are immeasurable in practice and may change rapidly during maneuvers which pose additional challenges for control [22]. Due to the difficulties in unsteady aerodynamic modelling, unsteady effects are usually treated as a part of model uncertainties which increase robust requirements for control system design. Based on an unsteady aerodynamic model, an alternative way is to estimate the unsteady effects via state observation and then make compensation in the controller design. Generally, it is not difficult to merely handle state observation, but the complexity increases drastically when uncertainties and input constraints are all considered. In recent years, a novel approach for stabilization and adaptive control of uncertain nonlinear systems based on immersion and invariance (I&I) methodology has been proposed in [23] and then further developed in [24, 25]. This method gives a noncertainty equivalent adaptive (NCEA) law which is different from the traditional certainty equivalent adaptive (CEA) method. The main feature of this approach lies in the construction of an estimator, which is a sum of a partial estimate generated by an update law and a judiciously chosen nonlinear function. The essential idea is to create a manifold in the extended space of states and parameters to which trajectories are attracted. As the trajectory evolves on this manifold, the closed-loop system captures the behavior of a deterministic system [26]. Due to the advantages in prescribing estimate error dynamics and separately synthesizing controller and estimator, this method shows great potential for uncertain nonlinear systems with complex structures. However, it relies on solving a partial differential equation (PDE), which is difficult for multivariable systems [27]. To overcome this difficulty, auxiliary state filter and regressor matrix filter are introduced [2830]. Recently, its applications to the control of missiles, quad-rotor UAVs, and satellites are studied successively [26, 31, 32]. Apart from flight control, the I&I approach is also applied in aeroelastic control [33] and mechatronic systems [34].

In this research, a novel adaptive backstepping control scheme based on noncertainty equivalent principle is proposed for advanced fighter attitude tracking. Compared with our preceding study in [35], the prior hypothesis about the upper bounds of the uncertainties is removed, and the input constraints are dealt with in a simpler way. Moreover, the coupling relationship of the observer and estimator is considered fully in the closed-loop design. Therefore, the main contributions of this paper can be summarized as follows: (1)The I&I method is used to estimate not only immeasurable states but also unknown parameters; it is further modified with different forms in these two situations and is closely associated with the characteristics of an advanced fighter(2)The mutual effects of unsteady states and parameter uncertainties make the closed-loop control design complex. The observer and estimator are coupled together, but the coupling effects can be weakened by selecting appropriate parameters; in other words, the observer and estimator can be designed to satisfy the spectrum separation principle(3)The proposed control scheme involves a state observer, a constraint backstepping controller, and two parameter estimators which are designed separately, but the stability of the closed-loop system is proven in a unified architecture

The rest of this paper is organized as follows. In the next section, unsteady attitude dynamics is modeled. The third section details the derivation of the state observer. In the fourth section, the constraint adaptive backstepping control design and stability analysis are presented. Simulations and conclusions are given in fifth and sixth sections, respectively.

2. Unsteady Attitude Dynamics

The model adopted in this research is a certain advanced fighter which is developed to investigate the flight dynamics during high AOA maneuvers. The fighter attitude dynamics is formulated as where , represent the rotation matrices from body axis system to wind axis system and stability axis system, respectively, and are inertia parameters defined in [21].

The fighter is equipped with two kinds of actuators, i.e., the aerodynamic surfaces and thrust vector engines, all of which are modeled as a first-order filter with limits in magnitude and rate (see Table 1). The sketch of the thrust vector control is depicted in Figure 1.

The thrust forces and moments are defined as

The aerodynamic model is extracted from wind tunnel tests which are conducted on sufficiently close points to capture the nonlinear behavior of the aerodynamics. Based on linear superposition principle, the aerodynamic forces and moments are expressed by where represent the unsteady states, and the expressions of are given in [36]. Unsteady effects on the axial force are neglected since the axial force measured in the test is small in contrast to other aerodynamic forces [36].

Usually, the quasisteady model in which unsteady states are equal to zero can be employed to describe the aerodynamics in normal flight condition. However, owing to notable effects caused by separated and vortical flow, it becomes inadequate to describe the aircraft dynamics at high AOA. Thus, unsteady aerodynamics are introduced and modeled in the following form [37, 38]: where is the immeasurable unsteady state; ; , ; and are fitting polynomials given in [36]. Note that the unsteady effects will converge to zero exponentially as the maneuver ends since (4) is stable.

3. Unsteady State Observation

The main objective of this section is to develop a nonlinear state observer which precisely reconstructs the unsteady states in the presence of parameter uncertainties. To this end, the attitude dynamics used for state observation and control design are rewritten as where is the incidence angle, represents the stability axis angular rate, denotes the measurable state, , and are constant parameter uncertainties, represents the intermediate control variable, denotes the control input, and are vectors or matrices with appropriate dimensions (see appendix).

First of all, related assumptions used in the subsequent developments are given below.

Assumption 1. The desired trajectories and their derivatives are bounded. The compact set is defined as follows: where is a known positive constant and stands for the absolute value of a scalar.

Assumption 2. For matrices , there exist positive constants such that , with compact subset containing the origin, and means the two norms of a vector or matrix.

To achieve precise reconstruction of unsteady states, a state filter is first introduced: where are positive definite matrices to be designed, denotes the estimate of , are filter errors, and are estimates of , respectively. According to (5) and (7), dynamics of take the form as where denotes the reconstruction error, and are estimate errors of uncertain parameters.

The state observer based on I&I theory is proposed as follows: where denotes the observer state, is a smooth nonlinear function, , . Taking the time derivative of and using (8), we can get

The nonlinear function is selected so that has a stable behavior. So we choose with . Therefore, the nonlinear function is selected as

Then, substituting (12) into (10), we obtain

In (13), contains the immeasurable state and uncertain parameters which cannot be used directly. To overcome these difficulties, a command filter is employed to provide filtered derivative of ; the filter dynamics is expressed by in which the initial condition satisfies . Based on (13) and (14), the adaptation law is designed as

Substituting (15) into (13), we can get in which . According to Assumption 2, there exists a positive constant which depends on such that , [39].

Consider the following Lyapunov function candidate

Using (8) and (16), the time derivative of yields

Using Young’s inequality, we can get where . Then, (18) is modified as with and is the minimum eigenvalue of matrix. Note if and the following inequalities hold then , which implies that the state observation process is asymptotically stable. In the following parameter estimation, will be determined to regulate .

Remark 3. The differences between the preceding studies [26, 27, 2933] and this research lie in the following two aspects: (1)The I&I approach adopted in this section is used to reconstruct the immeasurable and time-varying states(2)The filtered derivative of the regressor matrix, instead of filtered regressor matrix itself, is employed to construct the observer to avoid solving PDE

4. Constraint Adaptive Backstepping Control

In this section, a constraint adaptive backstepping controller based on I&I method is developed. The control objective is to track the predefined trajectories in the presence of input constraints and parameter uncertainties.

4.1. Constraint Backstepping Control

The constraint backstepping control design procedure is initiated by defining the following tracking errors: where is the reference trajectory generated by the desired command and is the virtual control law produced by .

To eliminate the “explosion of term” problem, is filtered via a command filter to provide and . The difference between and is evaluated by an auxiliary filter where , , is the filter state, and will be defined later. The compensated errors are defined by where , and denotes the reference trajectory of stability axis roll rate generated by . Taking time derivative of and using (5), we obtain

Therefore, the nominal virtual control law is given by

Combining (23), (24), (25), and (26), the dynamics of takes the form as

Similar to (26), the nominal control law for angular rate subsystem is designed as where is pseudoinverse of , , can be obtained by filtering .

To impose physical constraints on actuators, is filtered via a command filter to provide the physical limited . The command filter integrated with magnitude and rate limits is described as [14] where are magnitude and rate saturation function, respectively. The impact of (29) can be evaluated by where . Combining (5), (24), (28), and (30), the dynamics of is expressed by

Then, the following Lyapunov function candidate is considered:

Taking the time derivative of along (27) and (31) and using Young’s inequality, we have where .

4.2. Uncertain Parameter Estimation

Now the design of the parameter estimator is considered. This process is completed by deriving the nonlinear functions and adaptation laws for partial estimates.

The estimator of based on the I&I approach is constructed as where is the partial estimate of and denotes the nonlinear function of and . Taking the time derivative of and combining (27), we can get

To ensure that has a stable behavior, the nonlinear function is chosen as where . In view of (35) and (36), if is updated by then substituting (37) into (35), we can get

For angular rate subsystem, the estimator of is given as follows with where . Combining (39), (40), and (41), the dynamics of can be expressed by

Define the Lyapunov function candidate

Differentiating and invoking (38) and (42), we have

Remark 4. Assume that the state observer performs perfectly in reconstructing the unsteady states, then can be achieved. Hence, according to (38), the dynamics of can be treated as a linear time-varying system given by . Note that if at some instant , , then for . That is, the manifold defined by is an invariant manifold. Therefore, the parameter estimate remains frozen at its actual value for . This invariant characteristic also holds for dynamics.

Remark 5. Note that availability of is assumed in the estimator design; otherwise, we will need to solve PDE. Although cannot be measured directly, we could attempt to compute via velocity components and incidence angles.

Remark 6. According to (38) and (42), we can find that the state reconstruction error may affect the estimator performance, that is to say, the observer and estimator are coupled together. To weaken the effects between each other, the design parameters of the observer and estimator need to be chosen appropriately. In practice, the reconstruction error should converge much more quickly than the estimate errors to satisfy the spectrum separation principle.

4.3. Stability Analysis

The stability analysis of closed-loop system is initiated by the following theorem.

Theorem 7. Consider the closed-loop system under the foregoing assumptions, with application of the state filter (7), state observer (9), constraint backstepping controller (26) and (28), and the parameter estimator (34) and (39), , can be guaranteed to converge to zero.

Proof. The Lyapunov function is chosen as Combining (20), (33), and (44), the time derivative of yields To ensure stability of the closed-loop system, we select appropriate matrices and adaptation gain to satisfy the following inequalities: and inequalities (21), then can be achieved. Thus, we can conclude that all closed-loop signals are bounded based on Assumptions 1 and 2. Further, because is lower-bounded and monotonic by the negative-semidefiniteness of , we know that exists and is finite, which in turn implies Since and their derivatives are all bounded, according to Barbalat’s lemma, it follows that Thus, can be guaranteed to converge to zero.

According to Theorem 7, asymptotically converges to zero. Thus, the manifold defined by with is an attractive manifold. Along any trajectories evolving on , the dynamics of can be described by

Apparently, (52) can be obtained if a controller based on attitude dynamics with measurable unsteady states and known aerodynamic parameters is implemented. Hence, the constraint backstepping controllers (26) and (28) capture the performance of the deterministic control system. From (52), it can be easily seen

which implies that will converge to zero as expected.

Remark 8. For the closed-loop stability analysis, two scenarios are considered: (1)When actuator constraints are not in effect, note the fact that and are bounded, then the auxiliary filter states converge to a sufficient small domain around origin by selecting appropriate natural angular frequency and control matrices . Therefore, from (24), it is obvious that the tracking errors can be guaranteed to be bounded. Consequently, the stability of the closed-loop system can be ensured(2)When input saturation occurs, the primary control object is to maintain the stability of online approximation process [14]. It is obvious that the tracking errors may increase while the compensated errors still converge to zero when saturation happens. Through replacing the tracking errors by compensated errors in the parameter estimation, the stability of the adaptation process can be ensured not affected by the saturation

In conclusion, the proposed control scheme consists of three modules: unsteady states observation, constraint backstepping control, and parameter estimation. The overview of the design architecture is depicted in Figure 2.

5. Numerical Simulations and Discussions

This section presents several numerical simulations which are carried out to investigate the performance of the proposed control scheme. The comparisons of the NCEA and CEA laws are completed in the first simulation. The second simulation is performed to examine the effectiveness of the proposed control scheme via the Herbst maneuver. Both the simulations are conducted in the MATLAB/Simulink environment, where the step size is set to 5 ms. The parameters of the fighter are given in [36].

5.1. Comparisons of NCEA and CEA Laws

In this subsection, comparisons of the constraint filter backstepping controller based on the NCEA and CEA laws are performed to evaluate the adaptation effects on the closed-loop system. The parameters of the proposed controller are chosen as , . The initial estimate values of the unsteady states and uncertain parameters are all set to zero. Besides the desired command of AOA, the commands of are set to zero to neglect the lateral motion for simplicity. Meanwhile, parameter uncertainties are added into to simulate model errors. The equilibrium states at are chosen as the initial condition.

To illustrate the essence of the NCEA law, adaptive backstepping control approach proposed in [14] is adopted for comparison, the CEA laws are given as where are adaptation gains to be assigned. Consider that the transient control performance is highly dependent on adaptation gains, four sets of gains with different adaptation rates are selected, that is,

Case 1. (slow): .

Case 2. (medium): .

Case 3. (fast): .

In order to compare the performance of different adaptive laws, the root-mean-square (RSME) of tracking errors and the energy consumption (EC) are chosen as the evaluation indexes which are defined as where the superscript means the sampling point, denotes the number of samples, denotes the final time of simulation, and are the commands of total control moment in each channel.

Simulation results are shown in Figures 36 and Table 2. It should be pointed out that the unsteady states are perfectly reconstructed with fast converge rates for both adaptive laws; therefore, the results are not presented herein. From Figures 3 and 4 and Table 2, it is apparent that the NCEA law outperforms the three CEA laws both in command tracking and energy consumption. The track performance of the CEA laws is sensitive to the adaptation gains. When compared with the case 1, the CEA law in case 3 achieves higher estimate precision and faster converging rate at the expense of sacrificing the energy consumption and transient control performance. It can be observed from Figure 5 that are converged to zero; thus, the closed-loop trajectory eventually confines to the manifold , and the NCEA law recovers the performance of the deterministic backstepping controller. Figure 6 shows the time history of parameter estimation; it is obvious that the parameters estimated by CEA laws drift from their true values, especially in case 1. For the NCEA law, can be precisely estimated while this is not the case for . The reason is that the rank of equals to , then means ; therefore, the invariant manifold is achieved and can be estimated accurately. However, the rank of is less than ; hence, cannot be guaranteed via . In addition, it can be seen that is frozen after at which instant the attractive manifold is achieved.

5.2. Robust Analysis under Parameter Uncertainties and Unsteady Effects

To justify the effectiveness of the proposed control scheme, three different cases are considered for comparison. In the first case, the controller proposed in Section 4 is adopted in which the effects of unsteady aerodynamics and uncertain parameters are both compensated. For the second case, the effects of uncertain parameters are not considered in the controller design. The third case is the worst case in which both the unsteady effects and uncertain parameters are not compensated. The design parameters of the controller and initial condition are all kept the same as that in the previous subsection. The uncertainties of aerodynamic forces and moments in longitudinal and lateral channels are all considered in this simulation. Moreover, the Herbst maneuver is implemented to validate the control performance; the commands of are defined in [35].

The commands tracked in these three cases are compared in Figures 7 and 8, and the control performance of the proposed adaptive backstepping controller in the first case is presented in Figures 914. In Figure 14, are components of , respectively; are components of , respectively. It can be observed from Figures 7 and 8 that the tracking performance of the first case outperforms the other two cases. It is apparent that smaller steady tracking errors and better transient tracking performances can be achieved since the effects of unsteady aerodynamics and uncertain parameters are both compensated. It also can be concluded that the unsteady aerodynamics seriously affects the aircraft dynamics and the compensation of unsteady effects could significantly improve the control performance. Note that the thrust force presented in Figure 9 increases drastically in the maneuver; it is quite rational for the pilot to increase the throttle to provide a larger thrust force when the fighter operates at high AOA. In addition, the aerodynamic surfaces are all set to zero when and the attitude control was completed via thrust vector merely in this deep stall area. From Figures 1012, it can be seen that the tracking errors increase while the compensated errors still converge to zero even though thrust vector nozzles are saturated. Moreover, the proposed observer is effective in state reconstruction with strong robustness against parameter uncertainties, especially in the lateral channel. Meanwhile, the parameter estimator works fairly well in the presence of unsteady effects. The most essential reason can be concluded that the converging rate of state reconstruction error is much larger than the estimate error, thus the mutual effects nearly can be ignored. Figure 13 shows the uncertain parameters in aerodynamic moments are estimated exactly, but steady errors exist in the estimation of aerodynamic force parameters. This is consistent with the discussion in previous subsection. However, as presented in Figure 14, , and all converge to zero after ; the attractive manifold is achieved and the norms of the estimated states and parameters are maintained. Therefore, the performance of the deterministic control system is captured by the proposed constraint adaptive backstepping controller.

6. Conclusions

This paper is devoted to the constraint adaptive backstepping control design for the advanced fighter, where robust tracking of the predefined trajectories has been achieved. In particular, unsteady effects, parameter uncertainties, and input constraints are all considered which brings forward a great challenge for the control design. The proposed control scheme involves a state observer, a constraint backstepping controller, and two parameter estimators, which are designed separately. It should be pointed out that both the state observer and parameter estimator are designed based on the novel noncertainty equivalent principle which renders the closed-loop system with deterministic control performance. Moreover, stability of the closed-loop system can be guaranteed even though the saturation occurs. It can be concluded from simulation results that the backstepping controller integrated with the NCEA law outperforms the traditional adaptive backstepping controller. The effectiveness and robustness of the proposed control scheme are also verified by the Herbst maneuver. This paper mainly focuses on the theoretical aspect; our future work will be devoted to the practical applications of the flight control of scaled aircraft.

Appendix

The expressions of are given as follows:

Nomenclature

:Aircraft mass
:Airspeed
:Wing span and mean aerodynamic chord
:Dynamic pressure
:Reference area
:Product moment of inertia
:Total thrust force
:Angle of attack and sideslip angle
:Body-axis roll, pitch, and yaw rates
:Stability-axis roll, pitch, and yaw rates
:Control derivatives
:Roll, pitch, and yaw moments of inertia
:Aerodynamic moments along body axis system
:Thrust moments along body axis system
:Gravity components along body axis system
:Aileron, elevator, and rudder deflections
:Left nozzle deflections in yaw and pitch plane
:Right nozzle deflections in yaw and pitch plane
:Thrust point x-axis and z-axis locations
:Thrust forces along body axis system
:Aerodynamic forces along body axis system.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was cosupported in part by the National Natural Science Foundation of China under Grants 61673240 and 61873278.