Robust Backstepping Control of Wing Rock Using Disturbance Observer

Wing rock is a highly nonlinear phenomenon when the aircraft suffers undesired roll-dominated oscillatory at high angle of attack (AOA). Considering the strong nonlinear and unsteady aerodynamic characteristics, an uncertain multi-input and multi-output (MIMO) nonlinear wing rock model is studied, and system uncertainties, unsteady aerodynamic disturbances and external disturbances are considered in the design of wing rock control law. To handle the problem of multiple disturbances, a robust control scheme is proposed based on the extended state observer (ESO) and the radial basis function neural network (RBFNN) technique. Considering that the effectiveness of actuators are greatly decreased at high AOA, the input saturation problem is also handled by constructing a corresponding auxiliary system. Based on the improved ESO and the auxiliary system, a robust backstepping control law is proposed for the wing rock control. In addition, the dynamic surface control (DSC) technique is introduced to avoid the tedious computations of time derivatives for the virtual control laws in the backstepping method. The stability of the closed-loop system is guaranteed via rigorously Lyapunov analysis. Finally, simulation results are presented to illustrate the effectiveness of the ESO and the proposed wing rock control approach.


Introduction
Wing rock motion is an undesired roll-dominated oscillatory coupled with yaw oscillation for the aircraft flying at subsonic speed and high angle of attack (AOA).When an aircraft suffers the wing rock motion, the large amplitude and high frequency oscillatory of the rolling and yawing angle will severely limit the maneuverability, tracking accuracy, and operational safety of the aircraft.Hence, the wing rock phenomenon should be paid more attention, and the corresponding wing rock controller should be further studied.
In past decades, many researches have been done to analyze the mechanism of the wing rock.Mathematically, three possible nonlinear factors were considered as the inducements of the wing rock phenomenon, which are the variation of damping in roll with angle of sideslip, cubic variation of lateral derivatives with roll rate and sideslip, and the aerodynamic hysteresis in steady-state rolling moment [1,2].Meanwhile, different wing-rock models have been developed [3][4][5][6].Based on these simplified models, many control schemes were developed to suppress the influence of wing rock.Adaptive feedback linearization technique was employed for the wing rock in [7,8].An L 1 adaptive controller was designed for the wing-body rock motion [9].Fuzzy adaptive approaches have been employed in the wing rock during the last decade [10][11][12].The neural network technique and backstepping approach have been also extensively studied [13][14][15].In [16][17][18], simplified single degree-of-freedom wing-rock models were studied, and the disturbance observer or traditional extended state observer (ESO) was adopted in the design of robust wing rock control law.However, the yaw oscillation was often ignored in the design of wing rock controller for the convenience of research.Considering the serious coupling effect between the roll and yaw motion, multi-degree-of-freedom wing rock models should be further studied.
In this paper, an uncertain multi-input and multi-output (MIMO) nonlinear wing rock model is studied, and the coupling effect between the roll oscillation and yaw oscillation will be fully considered in the design of the wing rock controller.Meanwhile, to enhance the system robustness, the ESO and radial basis function neural network (RBFNN) techniques are combined to tackle the external disturbances and system uncertainties, respectively.Up to now, more and more researchers have paid attention to the study of the ESO or RBFNN based control techniques [19][20][21].In [22], an adaptive neural network based control law was developed for an uncertain MIMO system with unknown control coefficient matrices and input nonlinearities.In [23], the observer and adaptive neural network (NN) techniques were combined to tackle the uncertainties.To tackle the state constraint problem, the adaptive neural network technique was employed for an uncertain robot system [24].And the ESO technique was firstly proposed in the active disturbance rejection control (ADRC) [25].It has been successfully used in many challenging engineering problems [26,27].The main advantage of the traditional ESO is that it can handle various disturbances with little model information, and it takes both structured and unstructured uncertainties as generalized disturbances.Hence, considering the satisfactory uncertainty estimation performance, the ESO and RBFNN are fully combined in the design of the robust wing-rock control law.Meanwhile, except for uncertainties, the input saturation problem will also cause adverse effects on the performance of the wing rock control.
Input saturation is a challenging problem in the wing rock control.With the increase of AOA, the aircraft will lose the control efficiency gradually, which severely degrades the closed-loop system performance.During the past decades, there exist extensive researches on the control of mechanical systems with input saturation [28][29][30].In [31], an adaptive neural tracking control was considered for a class of stochastic nonlinear systems with input saturation, and a smooth nonaffine function of the control input signal was designed to approximate the input saturation function.In [32], a modified fault tolerant control law was designed to handle the input saturation problem.Constrained adaptive backstepping control was proposed based on the command filters in [33].Inspired by [34], the auxiliary system is designed to weaken the influence of the input saturation based on the backstepping technique in this paper.
The backstepping design method has been widely employed in the control of nonlinear systems.In recent years, many robust control methods have been introduced to combine with backstepping technique for the uncertain MIMO nonlinear systems [35][36][37][38][39].However, the ESO-based robust backstepping technique should be further developed for the wing-rock motion due to the relatively easy realization.Inspired by the above discussions, a robust attitude tracking control law is proposed for the wing-rock motion in the presence of unsteady aerodynamic disturbances, external disturbances, system uncertainties and input saturation.
In this paper, considering the characters of different disturbances, efficient processing methods are adopted.The system uncertainty is estimated online by adaptive RBFNN.By taking full advantage of the output of the RBFNN and the known model information of the unsteady aerodynamic disturbance, the ESO is employed to estimate the compounded disturbance, which consists of the unsteady aerodynamic disturbance, external disturbance, and the unknown neural network approximation error.By making full use of the known information of the system, the disturbance estimation performance will be greatly improved.In addition, it is meaningful in the wing rock control at high AOA.
The paper is organized as follows.Problem formulation and preliminaries are described in Section 2. Section 3 presents the detailed design process of the disturbance observers.In Section 4, the robust backstepping wing rock control law is developed considering the input saturation.In Section 5, simulation results are given to demonstrate the effectiveness of the proposed robust tracking control scheme, followed by concluding remarks in Section 6.

Problem Statement
According to the typical characteristics of wing rock motion, the following nonlinear MIMO lateral attitude motion model is considered [40]: where is an attitude angle vector which includes sideslip angle and flight-path roll angle, respectively; x 2 = [p, r] T ∈ R 2 is the body-axis angular rate vector; y is the output; x 2 ) are the system uncertainties which include the modeling errors and the perturbation of aerodynamic coefficients; δ a , δ y T denotes the practical control input; δ a is the aileron deflection angle and δ y is the lateral thrust vectoring angle; is the unsteady moment aerodynamic coefficient vector, and h l , h n are unsteady rolling moment aerodynamic coefficient and unsteady yawing moment aerodynamic coefficient, respectively.On the basis of work in [41], the unsteady aerodynamic coefficient can be modeled by first-order linear differential equation which is given by where At high AOA, the lateral attitude motion of (1) becomes unstable.When system (1) suffers no disturbances, simulation results for an angle of attack of 55 1 and 2. From Figures 1 and 2, it can be observed that the roll angle exhibits large amplitude oscillatory motion, and the sideslip angle motion also exhibits considerable oscillation.Meanwhile, the maximum amplitude of the sideslip motion is closed to 20 • , which is very dangerous for aircrafts at high AOA.When system (1) suffers the disturbances, simulation results for an angle of attack of 55 3  and 4.And the equations of the system uncertainty, unsteady aerodynamic disturbance and external disturbance will be given in details in the simulation section.In Figure 3 and 4, the states of the aircraft diverge quickly, then stable wing rock appears after 50 s.Comparing Figures 3 and 4 with Figures 1  and 2, it can be concluded that the disturbances will further aggravate the oscillation.Thus, it is very important to design a robust controller to suppress the wing-rock motion, and render the roll angle and sideslip angle converge to zero quickly.

t(s)
Considering that the wing rock motion is characterized by uncommanded roll-dominated oscillatory coupled with yaw oscillation, the aileron deflection angle δ a and the lateral thrust vectoring angle δ y are chosen as the practical control input.However, with the increase of AOA, the aileron control efficiency is reduced drastically [40].Hence, the input saturation must be taken into consideration in the wing rock control.And the input saturation model can be written as where u M i is the bounds of u i (t), i = 1, 2.

Neural Networks
The RBFNN has been widely employed to tackle unknown nonlinear functions due to its satisfactory approximation ability.Considering an unknown function f (x) : R m → R, the RBFNN approximation output over a compact set Ω m ⊆ R m can be written as [22] where W * ∈ R p is the optimal weight value and ι * is the smallest approximation error; is a basis function vector, which can be written in the form of where c i and b i are the center and width of the neural cell of the ith hidden layer.
The optimal weight value of the RBFNN is given by [42] W * = arg min where Ŵ is the estimation of the optimal weight value W * ; Ω f = Ŵ : Ŵ ≤ M is a valid field of the parameter and M is a design parameter; Ω m is an allowable set of the state vector.Using the optimal weight value yields where ῑ > 0.

Review of the Extended State Observer
An uncertain nonlinear single-input and single-output dynamic system is given by [43] where d(t) is the external disturbance; b = 0 is the control gain; f y (n−1) (t) , • • • , y (t) , d (t) is the unknown nonlinear function.
Treat f as an extended system state and define χ n+1 = f .Then, the system (8) can be written in the following augmented state space form [44]: where h represents the time derivative of , u ∈ R and y ∈ R are the state, input and output of the system, respectively.Hence, an ESO can be constructed as [44] ż1  [45,46], we can obtain In this paper, the RBFNN and ESO technique are combined to improve the robustness of the wing rock controller.Based on the output of the RBFNN and ESO, a robust backsteping attitude control law will be developed.To promote the control law design, the following lemma and assumptions are required.Lemma 1. [22]: For bounded initial conditions, if there exists a C 1 continuous and positive-definite Lyapunov function Assumption 4. For the uncertain nonlinear MIMO lateral model ( 1), there exist unknown positive constants λi such that λ max (g i g T i ) ≤ λi (i = 1, 2) where λ max (•) denotes the maximum eigenvalue.
Remark 1. Considering that the wing rock motion is a typical lateral attitude motion, the flight speed can be assumed to be constant.In practice, the wing rock motion is a limit-cycle oscillation, and the aerodynamic angle β is bounded in the wing rock motion.Hence, invoking the expression of the gain matrices g 1 (x 1 ), and g 2 (x 1 , x 2 ) in [40], Assumption 4 is reasonable.
Remark 2. The unsteady aerodynamic effect denoted by h 2 is induced by the periodical separation and attachment of the vortical flow.Meanwhile, the unsteady aerodynamic effect is also limited by the breakdown of the vortical flow.Hence, h 2 is bounded in practice.

Design of Disturbance Observer for the First Subsystem
According to (1), the first subsystem can be written as In ( 12), the system uncertainty ∆ f 1 (x 1 ) and external disturbance d 1 (t) are considered simultaneously.In traditional ESO method, ∆ f 1 (x 1 ) and d 1 (t) are merged into a single equivalent uncertainty.To improve the disturbance estimation performance, it is necessary to introduce different mechanisms to handle the system uncertainty and the external disturbance separately.
To efficiently tackle the unknown function vector ∆ f 1 (x 1 ), the RBFNN is employed to approximate it.The approximation output of the RBFNN can be written as where 13) into (12), we obtain where Extending D 1 as an additional state variable, i.e., defining x 3 = D 1 , ( 14) can be rewritten as where ] T is an unknown bounded time derivative function vector of D 1 .
According to the structure of extended system (15), a linear ESO combined with the RBFNN can be constructed as where x1 , x3 are the states of the observer (16); x1 is the estimation of Invoking (15) and ( 16), we obtain . Then ( 17) can be written as where Hence, it can be easily proved that A si is Hurwitz, if , and k si is a positive constant.Thus, there exists a positive definite matrix P si , i = 1, 2 satisfying the following Lyapunov equation: In this subsection, the disturbance observer (16) for the first subsystem ( 12) has been proposed.Selecting the parameters of ( 16) by the given method, A si has been proved to be Hurwitz.And the positive definite matrix P si , i = 1, 2 will be used in the design of the Lyapunov function (38), which is very important in the stability analysis of the proposed observer ( 16) and the closed-loop system.

Design of Disturbance Observer for the Second Subsystem
According to (1), the second subsystem can be written as It is more complicated that multiple disturbances are considered in (20).∆ f 2 (x 1 , x 2 ) is the unmodeled system uncertainty, which is an unknown function of the system states.h 2 is the unmeasurable unsteady moment aerodynamic coefficient vector whose model can be achieved by wind tunnel experiments.And d 2 (t) represents the external disturbance.Similarly, different mechanisms are adopted to handle these different disturbances in this section.Meanwhile, to improve the disturbance estimation performance, the model information of the unmeasurable unsteady moment aerodynamic coefficient vector h 2 is fully utilized.
Invoking ( 2) and ( 20), the augmented second subsystem state equations can be written as where G 2 = diag {G 21 , G 22 }, and G 21 > 0, G 22 > 0 are constants.Similar to (13), considering the neural network approximation of ∆ f 2 (x 1 , x 2 ), (21) can be written as where (14), the time derivative of x 4 is given by Then, invoking (23), the extended system of ( 22) can be written as where h2 = [h 21 , h22 ] T is an unknown bounded time derivative function vector of x 5 .
Then, a linear type of the ESO can be constructed as where x2 is the estimation of x 2 ; x4 is the estimation of x 4 ; x5 is the estimation of Invoking ( 24) and ( 25), we have h2i (t) (27)   where According to the Routh-Hurwitz criterion, it can be easily proved that if the parameters of ( 25) can be chosen as , and k f 1 > 0, k f 2 > 0 are design constants, A f i will be Hurwitz.Thus, there exists a positive definite matrix P f i , i = 1, 2 satisfying the following Lyapunov equation: In this subsection, the disturbance observer (25) for the second subsystem (20) has been proposed.Selecting the parameters of (25) by the given method, A f i has been proved to be Hurwitz.And the positive definite matrix P f i , i = 1, 2 will be used in the design of the Lyapunov function (44), which is very important in the stability analysis of the proposed observer (25) and the closed-loop system.Now, we have illustrated the design process of disturbance observers ( 16) and (25).However, according to (18) and ( 27), we still can't prove their stability due to the coupling between the observers and the RBFNN which will be design later.Hence, the stability of the disturbance observers will be proved in the stability analysis of the closed-loop system.

Design of Robust Backstepping Attitude Control Based on Disturbance Observer
At high AOA, the aerodynamic characteristics of the flow field around an aircraft is quite complicated.Any small changes of the attitude of an aircraft may cause drastic changes of the aerodynamic force and aerodynamic moment, which will severely degrade the aircraft flight performance.Hence, the robust backstepping attitude tracking control law should be further developed for the aircraft at high AOA.To compensate for the matched and unmatched disturbances, the ESOs proposed in Section 3 are employed to estimate them.And the disturbance estimations will be introduced into the design of the virtual control law and the practical control law.
Invoking ( 15) and ( 24), the system (1) can be written as Considering the input saturation problem, the following auxiliary system is constructed: where > 0 are design matrixes.Considering the backstepping control law design method, we define where x 2c is the virtual control law which will be designed later.
In the following, we consider the robust backstepping attitude control design for the wing-rock motion.In the attitude control design, the adaptive RBFNN and the ESO are combined, and the detailed design process is appended as follows.
Step 1: Considering the first equation in (29) and differentiating e 1 with respect to time yields The ideal virtual control law x2c for x 2 is designed as where c 1 = c T 1 > 0 is a matrix to be designed.And the weight value of RBFNN is updated by where The dynamic surface control (DSC) technique technique is adopted to obtain the approximation of the time derivative of x 2c .Introducing a first-order filter with time constant Γ 1 = diag {τ 11 , τ 12 }, and passing x2c through it yields Defining ε f 1 = x 2c − x2c and invoking Equation ( 35), we have where M 1 (•) is the sufficiently smooth function vector about Π 0 and Π 1 : (x 1 , Ŵ1 , x3 , ξ 1 , ξ 2 ).Since the set Π 0 and Π 1 are compacts, Π 0 × Π 1 is also compact.Therefore, M 1 (•) has a maximum M1 on Then, we have Choose the Lyapunov function candidate as Define ε 13i = [ε 1i , ε 3i ] T , then the time derivative of V 1 can be written as where P si C si = M PC si , h1i ≤ M h1i , P si B si S 1 (J 1 ) ≤ M PBS1 , and τ > 0 is a design constant.The detailed analysis of the time derivative of V 1 is shown in Appendix A.
Step 2: Invoking (29), differentiating e 2 yields ė2 The control law v is designed as where c 2 = c T 2 > 0 is a design matrix, and Ŵ2 is updated by where ρ 2 > 0, Λ 2 = Λ T 2 > 0 are design parameters.Substituting the controller v into Equation (40), we obtain where W2 = Ŵ2 − W * 2 , x4 = x4 − x 4 are the estimation errors.Define ε 245i = [ε 2i , ε 4i , ε 5i ] T , i = 1, 2, and consider the Lyapunov function candidate as The time derivative of V 2 is given by where The detailed analysis of the time derivative of V 2 is shown in Appendix B. Then, the above robust backstepping wing rock control design procedure can be summarized in the following theorem.Theorem 1. Considering the nonlinear MIMO lateral model (1) in the presence of unsteady aerodynamics disturbances, external disturbances, system uncertainties, and input saturation, the disturbance observers are designed as ( 16), (25), the auxiliary system is designed as (30), the filter is designed as (35), and the robust backstepping attitude control law is proposed as (41).Then, the attitude tracking error is convergent and bounded under the proposed adaptive robust backstepping control scheme.

Proof.
To analyze the stability of the closed-loop system, the following Lyapunov function candidate is chosen as Then, differentiating V and invoking Equations ( 39), ( 45), we have According to Equation ( 47) and Lemma 1, we obtain Considering Equations ( 31), ( 46), (49), we obtain Then, we obtain Hence, it is obvious that the attitude tracking error is convergent and bounded under the proposed adaptive robust backstepping control scheme.Similarly, we can conclude that the disturbance estimation errors ε 13i , ε 245i , i = 1, 2 of the disturbance observers ( 16) and ( 25) are also convergent and bounded.
This concludes the proof.

Simulation Study
In this section, the simulation results are given to verify the effectiveness of the proposed robust control scheme in suppression of wing rock motion at high AOA.In this section, the aircraft model in [40] and the unsteady model in [41] are adopted to demonstrate the attitude control performance and the disturbance rejection performance of the above designed controller.
In simulations, initial conditions for the aircraft are take as α = 55 whereas initial conditions for the proposed ESO are taken as xT 1 (0) , xT 3 (0) = [5.7, 0, 0, 0],xT 2 (0) , xT 4 (0) , xT 5 (0) = [0, 0, 0, 0, 0, 0].And the initial condition for the auxiliary system is chosen as ξ T 1 , ξ T 2 = [0, 0, 0, 0].The time constant for ( 35) is designed as Γ 1 = diag {0.08, 0.08}.It is very important to choose suitable parameters for the ESO.And the observer gains of the error dynamics ( 18) and ( 27 [ω 21 , ω 22 ] = [6, 6].Large k si , k f i , ω 1i , ω 2i , i = 1, 2 will improve the disturbance performance, and make the disturbance estimation errors arbitrarily small.Meanwhile, large k si and k f i will reduce the conservatism of the controller proposed in this paper. The parameters for the controller ( 41) are chosen as c T 1 , c T 2 = [6, 6, 10, 10].Large control parameters c T 1 , c T 2 will improve the control performance and reduce steady-state tracking errors.However, large control parameters will also induce serious input saturation.The parameters for the auxiliary system ( 30) is designed as k T aux1 , k T aux2 = [6, 6, 60, 60].The parameters for the adaptive RBFNNs ( 34) and ( 42) are chosen as In the simulation study, the system uncertainties ∆ f 1 (x 1 ) and ∆ f 2 (x 1 , x 2 ) are assumed as +40% variation of the force and moment aerodynamic coefficients.In addition, we assume that the unknown time-varying disturbances imposed on the aircraft are ḣl ḣn = − 5 0 0 5 0.01 (sin (0.5πt) + sin (0.2πt) + 0.1) 0.02 (sin (0.5πt) + sin (0.2πt) + 0.1) (53) According to Figures 1-4, the states of the aircraft diverge rapidly when the wing rock motion appears.Hence, it is necessary to adopt appropriate control scheme immediately after the wing rock phenomenon is detected.In the following simulations, we assume that the wing rock motion is detected at 5 s, and then the controller ( 41) is executed immediately.
To facilitate the description, we define the compounded uncertainties Ψ s , Ψ f as following Meanwhile, we define the two estimation errors e Ψs = [e Ψs1 , e Ψs2 ] T , e

Wing Rock Control
In this subsection, the control objective is to track a desired attitude command while suppressing the wing-rock motion.The simulations are given by considering a reference attitude angle trajectory as The controller is executed at 5 s after the wing-rock motion appears, and the simulation results are presented in Figures 5-14.
Figures 5 and 6 show that the wing rock motion appears, and the aircraft starts to oscillate at the beginning.Then, after the controller ( 41) is executed in the fifth second, the states µ, β converge to the desired attitude command (57) quickly, which means that the controller can achieve satisfactory tracking performance while suppressing the wing rock motion.Figures 7-10 show that the ESO can estimate the uncertainty quite accurately.In Figure 7, the compounded disturbance Ψ s1 diverged quickly in fifth seconds, and the proposed disturbance observer can still estimate the disturbance well.The estimation errors are plotted in Figures 11 and 12, which show that the proposed disturbance observer not only has good steady-state performance, but also has good transient performance.Meanwhile, the corresponding control input histories are given in Figures 13 and 14, respectively.
From the simulation results, we know that the developed ESO based robust control law is valid for the nonlinear MIMO lateral model of an aircraft with system uncertainty, unsteady aerodynamic disturbance, and input saturation.One of the focuses of this paper is the robust control of wing rock motion suffering multiple disturbances.To enhance the system robustness, an improved ESO is proposed combined with the RBFNN.Next, the simulation will be carried out to compare the performance of the proposed ESO with existing traditional ESOs.

Comparison with Traditional ESO
In the traditional ESO design method, all the uncertainties are merged into a single compounded uncertainty.Hence, we rewrite (4) as follows Considering that a linear ESO is developed in this paper, the traditional linear ESO (LESO) is considered for comparison instead of traditional nonlinear ESO.Meanwhile, the use of linear ESO offers certain advantages.Firstly, the close loop stability for nonlinear ESO is hard to established conclusively.Secondly, the LESO is easy from hardware implementation point of view.
In [18], the authors have presented the wing rock control law based on the linear ESO.Hence, according to [18], the LESO for (58) is designed as where xLs is the estimation of x 1 ; ΨsL = ΨsL1 , ΨsL2 T is the estimation of Ψ s ; e Ls = xLs − x 1 ; xL f is the estimation of x 2 ; The estimation output of the LESO (59) is presented in Figures 15-18.In (58), all the uncertainties are merged into a single compounded uncertainty, which leads to a larger value of the time derivative of the compounded uncertainty.Hence, obvious estimation errors can be observed in the output of LESO.Comparing Figures 11 and 12 with Figures 19 and 20, it is obvious that the disturbance observer proposed in this paper can achieve better estimation performance under the same conditions.Meanwhile, we can conclude that the developed disturbance observer not only has good steady-state performance, but also has good transient performance.

Wing Rock Control under Narrow-Band Disturbances
In many classes of applications like active vibration control and active noise control, the broadband and band-limited disturbances have been considered, which can be of narrow band type (simple or multiple) or of broad band type.To make the simulation more realistic, a narrow-band disturbance (NBD) will be considered in the wing rock control.Usually, a common framework is the assumption that the broadband and band-limited disturbances are the results of a white noise or a Dirac impulse passing through the shaping filter.Similar to [47], the transfer function of the shaping filter is taken to be s s 2 +0.95409s+9103 .The mean, variance and sample time of the white noise are set as 0, 150, and 0.1 s.And the output of the filter is the NBD which is indicated by d nb .
All the design parameters, the considered disturbances, and the reference attitude angle trajectory are the same as stated in Section 5.1, except for the external disturbance d 2 which is written as (60).Comparing (54) and (60), it is obvious that the NBD d nb is considered in the external disturbance.

Conclusions
In this paper, a robust attitude control scheme has been developed for wing-rock motion control with the system uncertainties, unsteady aerodynamic disturbances, external disturbances, and input saturation.To be more consistent with the practical wing rock phenomenon, a uncertain MIMO nonlinear wing rock model was studied.To address the issues of multiple disturbances, an improved ESO was developed combined with the RBFNN technique, which relaxed the known boundary requirement of the uncertainties and disturbances.Furthermore, the input saturation problem was considered, and a corresponding auxiliary system was designed to weaken the effect of the input saturation.Finally, a robust backstepping control scheme has been designed.Exclusive simulations have proven that the improved ESO can achieve satisfactory disturbance estimation performance, and that the proposed control scheme is valid in the wing rock control.
Furthermore, considering the following facts: the time derivative of V 1 can be written as ≤e T 1 g 1 (x 1 )e 2 − c 1 − where P si C si = M PC si , h1i ≤ M h1i , P si B si S 1 (J 1 ) ≤ M PBS1 , and τ > 0 is a design constant.This concludes the analysis.

Appendix B
The detailed analysis of the time derivative of Lyapunov function V 2 is given as following: The Lyapunov function V 2 is given by Then, (A8) can be written as V2 ≤ − e T 2 g T 1 (x where This concludes the analysis. and b 21 (β) > 0, b 22 (β) > 0, a 21 (β) , a 22 (β) are known bounded function variables of β.

where π 1 ,Assumption 1 .Assumption 2 .Assumption 3 .
π 2 : R n → R are class K functions, κ and c are positive constants, then the solution x(t) is uniformly bounded.All the states of the uncertain nonlinear MIMO lateral model (1) are measurable.For the uncertain nonlinear MIMO lateral model (1), the time derivatives of d i , i = 1, 2 are bounded.For the uncertain nonlinear MIMO system (1), the desired system output y d and the derivatives ẏd , ÿd are bounded.Namely, the set Π 0 = (y d , ẏd , ÿd ) : y d 2 + ẏd 2 + ÿd 2 ≤ M 0 exists, where M 0 > 0 is an unknown positive constant.
are the observer gain parameters, and φ i , i = 1, 2, • • • , n + 1 are the linear functions or nonlinear functions of (z 1 − y ) to be chosen.
If the functions φ i , i = 1, 2, • • • , n + 1 and the observer gain parameters l i The simulation results are presented in Figures21-28.From Figure21and 22, it can be observed that the proposed controller can still achieve satisfactory tracking performance in the presence of the NBD.The disturbance estimation results are presented in Figures23-26.Especially Figures25 and 26indicate that the observer proposed in this paper can still achieve satisfactory disturbance estimation performance in the presence of the NBD.And the corresponding control input histories are given in Figures27 and 28 Then, the time derivative of V 2 can be written asV2 = − e T 2 g T 1 (x 1 ) e 1 − e T T 245i Iε 245i + ε T 245i P f i B 1 f i WT 2i S 2 (J 2 ) − 1 T 245i P f i B 2 f i G 2i a 2i (β) WT 11 S 1 (J 1 ) − ξ T 2 k aux2 ξ 2 + ξ T 2 g 2 (x 1 , x 2 ) ∆u 245i + max {ω 2i } e T 2 e 2 − c 2 e T 2 e 2 1 ) e 1 − e T 245i + max {ω 2i } e T 2 e 2 − c 2 e T 2 e 2 (x 1 ) e 1 − (c 2 − max {ω 2i }) e T 2 e 2 −