Neural Extended State Observer Based Intelligent Integrated Guidance and Control for Hypersonic Flight

Abstract: Near-pace hypersonic flight has great potential in civil and military use due to its high speed and low cost. To optimize the design and improve the robustness, this paper focuses on the integrated guidance and control (IGC) design with nonlinear actuator dynamics in the terminal phase of hypersonic flight. Firstly, a nonlinear integrated guidance and control model is developed with saturated control surface deflection, and third-order actuator dynamics is considered. Secondly, a neural network is introduced using an extended state observer (ESO) design to estimate the complex model uncertainty, nonlinearity and disturbance. Thirdly, a command-filtered back-stepping controller is designed with flexible designed sliding surfaces to improve the terminal performance. In this process, hybrid command filters are implemented to avoid the influences of disturbances and repetitive derivation, meanwhile solving the problem of unknown control direction caused by nonlinear saturation. The stability of the closed-loop system is proved by the Lyapunov theory, and the controller parameters can be set according to the relevant remarks. Finally, a series of numerical simulations are presented to show the feasibility and validity of the proposed IGC scheme.


Introduction
Since the 1980s, integrated guidance and control (IGC) has been the subject of lots of studies for tactical missiles, space shuttles and hypersonic flight vehicles.Firstly, the conventional design is not able to produce an optimal overall system.Because the guidance and control system is formed in two independent loops with different operation frequencies, the original performance objectives are lost when the two control loops are iteratively incorporated.Secondly, in order to improve the closed-loop performance, the modern guidance and control (G&C) design has involved various techniques such as sliding mode control, adaptive control, disturbance rejection control, anti-windup control, etc., because the system dynamics is obviously characterized by strong states coupling, imprecisely known aerodynamics, complex internal uncertainties and external disturbances, especially in hypersonic flight.Consequently, the IGC design has provided better solutions in terms of low cost, modular growth, design flexibility, simple logistics and has attracted great interest from researchers.
Schierman et al. [1] introduced an integrated adaptive guidance and control (IAG&C) program for Boeing's X-40A.The development of future autonomous launch platforms will benefit from integrating trajectory reshaping with the essential elements of control reconfiguration and guidance adaptation.Shima et al. [2] considered the actuator dynamics as an important part of the integrated missile guidance and control system in a low-order element.Furthermore, two first-order actuators are considered in the IGC model of dual-control missiles in [3,4]; however, the controller complexity becomes staggering.Nevertheless, the true missile dynamics has high order for second-or third-order actuator dynamics.In [5], the sliding surface of a sliding mode integrated missile guidance and control is formed by Predicted-Impact-Point (PIP) heading error, and then a terminal second-order sliding mode control is designed to achieve the convergence in a finite time without chattering.In [6][7][8][9], an impact angle constraint sliding surface is proposed to extend the application and improve the terminal performance of IGC design.Disturbance observer has been widely used in dynamic model compensation, and the design of ESO has good flexibility with respect to its general form.The nonlinear ESO-based controller in [10] has better performance under complex uncertainties and measurement noises than linear ESO in [11] for missile guidance and control systems, while the nonlinear ESO has more parameters to tune.Huang et al. [12] used a nonlinear disturbance observer to estimate the nonlinearities, uncertainties, disturbances and compensate the dynamic model; thus, the undesired chattering in control signals is depressed and the robustness of closed-loop system is enhanced.Guo et al. [13] employed the ESOs to estimate indirectly measured states and cancel out various parametric uncertainties.
Neural networks are capable of providing arbitrarily good approximations of prescribed functions of a finite number of real variables in [14].Application to dynamic system modeling, nonlinear complex valued signal processing associated with the radial basis function (RBF) neural network are described by Wu et al. [15].In [16], direct neural control is proposed to deal with the input nonlinearity in the model of a hypersonic vehicle.An adaptive neural network is employed to estimate the structure uncertainties, and then a back-stepping controller is proposed to guarantee the uniform asymptotic stability of the uncertain system in [17].In [18], RBF neural network approximation is combined with an adaptive back-stepping technique to achieve boundedness of all closed-loop system states.In [19], RBF neural networks are applied to approximate the lumped unknown nonlinearities to satisfy robustness against system uncertainties for a constrained flexible air-breathing hypersonic vehicle.Nonlinear saturation is modelled as a non-smooth function, and approximated with a smooth function with bounded error in backstepping controller design in [20].Nussbaum function is used when differentiating the approximate function to guarantee the closed-loop stability.Chen et al. [21,22] introduced Nussbaum functions to compensate the nonlinear term due to differentiating the saturated control input and solve the problem of unknown control direction.In [22], a novel Nussbaum function gain-based controller is proposed for systems with multiple unknown actuator directions and time-varying nonlinearities.
To simplify the implementation and satisfy the constraint in back-stepping control techniques, a command filter is introduced.Farrell et al. [23] first introduced the command-filtered back-stepping technique to get the time derivatives of the virtual control signals, and it is feasible for high-order dynamic system.In [24,25], a low-pass command filter is employed to construct the derivative of virtual control input in dynamic surface control.It solves the problem of "explosion of complexity" caused by the repeated differentiations of the virtual control signals in the missile IGC design.The performance of command filter-based back-stepping control scheme significantly improves in stability and steady-state tracking accuracy, while the analysis via singular perturbation theory is made in detail in [26].Xu et al. [27] proposed a second-order nonlinear command filter to directly impose magnitude and rate limitation on the control input.Wang et al. [28] proposed composite command-filtered back-stepping control for an ESO compensating IGC model with second-order actuator, and a Nussbaum function gain-based command filter was developed to deal with the control input constraints.
In brief, this work is motivated by (1) integrated guidance and control design, (2) complex uncertainty and disturbance estimate using an improved extended state observer based on neural network, (3) composite command-filtered back-stepping control with nonlinear constraint.This paper proposes a novel IGC approach to synthetically enhance the stability and robustness of an overall closed-loop system.Then, the application for hypersonic flight is able to be studied and discussed more practically and effectively.
The paper is organized as follows: in Section 1, the representative study in integrated guidance and control design are summarized, the control technique involved is introduced.In Section 2, an integrated guidance and control model in the terminal phase of hypersonic flight is developed with high-order actuator dynamics, and the control surface deflection is constrained by magnitude saturation.In Section 3, a neural extended state observer (NESO) is designed to synthetically estimate the complex uncertainty and time-varying disturbance, and three NESOs are used in different channels.In Section 4, a hybrid command-filtered back-stepping controller is designed.Command filters based on tracking differentiators are used to get the derivate of virtual inputs to avoid chattering, and a command filter based on the Nussbaum function is introduced when differentiating the saturated control input with unknown control direction.In Section 5, the stability of closed-loop system is analyzed by the Lyapunov theory, and the controller parameters can be set by relevant principles in the proof.In Section 6, numerical simulations in five scenarios, and Monte Carlo simulations of initial flight states, measurement noises, aerodynamics coefficients are presented to verify the proposed design in widely practical situation.Finally, in Section 7, all the above work is concluded.

IGC Model
In this section, the integrated guidance and control model in the terminal phase is developed with actuator dynamics.Figure 1 shows the planar geometry under inertial coordinate system without earth rotation.V is the flight velocity, R is the distance between the mass center and the landing point, q is the angle between line of sight (LOS) and local horizon, γ is the flight path angle, α is the angle of attack (AOA) and δ z is the control surface deflection.The angle directions of q and γ shown in Figure 1 are negative and α is positive.
Energies 2018, 11, x FOR PEER REVIEW 3 of 18 overall closed-loop system.Then, the application for hypersonic flight is able to be studied and discussed more practically and effectively.The paper is organized as follows: in Section 1, the representative study in integrated guidance and control design are summarized, the control technique involved is introduced.In Section 2, an integrated guidance and control model in the terminal phase of hypersonic flight is developed with high-order actuator dynamics, and the control surface deflection is constrained by magnitude saturation.In Section 3, a neural extended state observer (NESO) is designed to synthetically estimate the complex uncertainty and time-varying disturbance, and three NESOs are used in different channels.In Section 4, a hybrid command-filtered back-stepping controller is designed.Command filters based on tracking differentiators are used to get the derivate of virtual inputs to avoid chattering, and a command filter based on the Nussbaum function is introduced when differentiating the saturated control input with unknown control direction.In Section 5, the stability of closed-loop system is analyzed by the Lyapunov theory, and the controller parameters can be set by relevant principles in the proof.In Section 6, numerical simulations in five scenarios, and Monte Carlo simulations of initial flight states, measurement noises, aerodynamics coefficients are presented to verify the proposed design in widely practical situation.Finally, in Section 7, all the above work is concluded.

IGC Model
In this section, the integrated guidance and control model in the terminal phase is developed with actuator dynamics.Figure 1 shows the planar geometry under inertial coordinate system without earth rotation.V is the flight velocity, R is the distance between the mass center and the landing point, q is the angle between line of sight (LOS) and local horizon, γ is the flight path angle, α is the angle of attack (AOA) and z δ is the control surface deflection.The angle directions of q and γ shown in Figure 1 are negative and α is positive.( ) ( ) ( ) ( ) ( ) ( ) The state vector is chosen as T , and the control surface deflection is saturated as v = sat(x 4 ).Then the state space form of integrated guidance and control model can be built as follows: . . .
where d i (t), i = 1, 2, 3 are time-varying disturbances in different channels.The nonlinear functions f i (x), i = 1, 2, 3 and k, b are depicted as follows: In hypersonic flight, control-oriented forms of force and moment functions are introduced with curve-fitted approximations in Appendix A. The lift force can be simplified into two parts: 2a) and (2b) are the unified form of the parts, where m is the mass.The aerodynamics moment is also described as two parts: both are unified by the moment of inertia with respect to body z-axis I z : 2c).The high-order nonlinearity and states coupling mainly caused by aerodynamics coefficients uncertainties are synthesized as The nonlinear saturation characteristic of the actuator in Equation (1c) can be modeled as: where δ z,max is the maximum value of control surface deflection.The saturated x 4 can be approximated by a smooth function defined as: then sat(x 4 ) = g(x 4 )+, is the small approximation error.The actuator dynamics is considered to be a first order inertial element and a second-order element.It is expressed with the following transfer function: where τ s is time constant, ξ s is damping rate, ω s is natural frequency and K s is constant gain.

Neural Extended State Observer
The extended state observer is widely implemented, with a simple structure and good estimation accuracy.Considering the following system with uncertainty and time-varying disturbance: .
The basic idea of ESO is to estimate the synthetical disturbances ∆ p = ∆ f p x p + d p (t).As the introduction section states, neural networks have good approximation performance in dealing with complex system uncertainties and disturbances.To compensate the complex aerodynamics uncertainties and time-varying disturbances in the IGC model, the technique of the extended state observer is exploited with a neural network.Defining ∆ NESO , x NESO as the states of NESO to estimate, respectively, ∆ p and x p , and the estimate error of x p as e p = x NESO − x p , then the state space model of NESO is designed as follows: . where ) is a vector of Gaussian functions, O is the construction error of neural network observer.
T is the center of the Gaussian functions vector, d is the affect size.By employing the σ-modification type of adaptation law, the weight vector W can be functioned as: .
where η > 0 is constant, Γ > 0 is gain.Defining x pr as the reference signal of x p , then the tracking error can be expressed as e p = x p − x pr .
In the integrated guidance and control design, three NESOs are used to compensate the synthetical disturbances

Hybrid Command-Filtered Controller
In this section, a hybrid command-filtered back-stepping sliding mode controller is designed to achieve LOS rate convergence based on the estimate results of NESOs.Three low-pass filters are used to get one-order derivatives of virtual control input; the Nussbaum function is introduced to deal with the nonlinear saturation of control surface deflection with one-order auxiliary dynamics, and then a nonlinear third-order differentiator is implemented in the terminal sliding mode control.
Step 1: The first sliding surface is defined as s 1 = x 1 .The virtual control input x 2 is designed according to proportional reaching law with positive constant k 1 .
Step 2: The second sliding surface is chosen as s 2 = x 2 − x 2 , and the virtual control input x 3 is obtained by the following control law, like in Step 1: where .
x 2r is the derivative of x 2 .Directly differentiating x 2 with respect to time results in chatter and peak of the virtual control signal, so a low-pass command filter is introduced to get the derivative: where τ 1 is the positive time constant.The state error produced by Equation ( 12) is defined as which can be sufficiently small when τ 1 is appropriately set, thus replacing s 2 with s 2 = x 2 − x 2r .A low-pass filter in the same form as Equation ( 12) is used to avoid directly differentiating x 3 .
Then the error produced by the low-pass Equation ( 13) can be expressed as λ 2 = x 3r − x 3 .
Energies 2018, 11, 2605 6 of 17 Step 3: In this step, there are two sub-steps to deal with the nonlinear saturation of control surface deflection.
Firstly, the virtual control input x 4 is designed to achieve the convergence of s 3 = x 3 − x 3r under the following control law: Secondly, existing a sliding surface s 4 = g(x 4 ) − x 4r , its derivative with respect to time is depicted as: .
where x 4r is obtained by the following low-pass filter with an error λ 3 = x 4r − x 4 under time constant τ 3 : To get the first-order derivative .
x 4 , x 4 is passed through a low-pass filter with time constant τ s .
In the process of reaching the law design of s 4 , the control direction ξ is unknown because of saturation.A Nussbaum function N(χ) command filter is introduced to design the virtual control input.The Nussbaum function is defining with the following properties: Then the virtual control input x 4u is established as follows: According to the above properties in Equation ( 18), the following Nussbaum function is implemented: The parameter χ is set as an adaptive parameter according to the following principle: .
In the end, a third-order differentiator is used as a command filter to get first-and second-order derivatives of the virtual control input x 4u : .
where ρ = x CF1 − x 4u + 1 ε x CF2 Step 4: The task of tracking x 4u is completed by a terminal sliding mode controller.From the second order dynamics in Equation ( 5), the following model is derived: .
The model uncertainty ∆ f s (x) and the time-varying disturbance d s (t) are bounded by positive constants F, D: By defining error vector E = e s , .e s T with e s = x 4s − x 4u , and its derivative can be written as .e s = .
x 4s − x CF2 .Then the following nonlinear sliding mode surface is designed: In Equation ( 26) ..
where T is the convergence time of the terminal sliding mode controller.Finally, the control input is designed as follows:
The following Lyapunov function is chosen: Differentiating the actuator dynamics related part of Equation ( 24) with respect to time: Energies 2018, 11, 2605 8 of 17 The differentiator can be set appropriately such that x CF3 − ..
x 4u ≤ CF , CF is a small positive constant, and ∆ f s , d s is bounded.There exists a positive constant K that yields to: Through differentiating Equation (30) and combining Equation (32), we have: .
If defining the following bounds: where .
According to the Young's inequality, Equation (36) can be rewritten as follows: .
, then the following inequality is yielded to: Integrating Equation (36) directly, then we have: According to the proof in [22], t 0 .
χ[ξ N(χ) − 1]e Cτ dτ is bounded, and V(t) is bounded, which implies that all the error is bounded.

Numerical Simulations
In this section, numerical simulations under five scenarios are provided to illustrate the control schemes proposed in the previous sections.First of all, the initial flight condition in terminal phase is set as follows: The IGC model is built in the following first three scenarios with aerodynamics coefficients based on table in Appendix A. The aerodynamics coefficients in Scenario 4 (S4) and Scenario 5 (S5) are the same as Scenario 3 (S3).The control surface deflection is saturated within δ z,max = 10 deg, and the third-order actuator dynamics is given by the following transfer function: Scenario 1: Energies 2018, 11, 2605 9 of 17 Scenario 2: Scenario 4: Furthermore, the un-modeled measurement uncertainty of .
q is considered to be a first order transfer function with gain K sensor = 0.75, and time constant τ sensor = 0.05.
Scenario 5: Dead-zone is also part of control input nonlinearity, which can be defined by the following piecewise function of δ z,max : The time-varying disturbances in different channels of IGC model are given as follows: Secondly, the three NESOs in Equation ( 7) are set.According to the form of nonlinear function to be estimated, the first neural network input is chosen as [x 1 , x 2 , v].The elements of the center vector are set based on the bound of the states.The second neural network input is chosen as [x 2 , v], according to the bound of the states, the following center vector is employed: The third neural network input is chosen as [x 2 , x 3 , v], according to the bound of the states, the following center vector is employed: The gains on estimate error in Equation ( 7) are set with β = [30, 30, 30].Thirdly, the proportional reaching law of each step in Section 4 is set with constants k 1 = 1.2, k 2 = 1.45, k 3 = 1.45.The low-pass command filters are built with time constants τ 1 = τ 2 = τ 3 = 0.0198s.The third-order differentiator in Equation ( 23) is set with ε = 0.04.The Nussbaum function-based command filter is confirmed by Equation (22) with γ χ = 10 −6 .The terminal sliding surface in Equation ( 26) is determined by c 1 = 5, c 2 = 1, and the uncertainty and time-varying disturbance in Equation ( 28) are bounded by F = 10, D = 3.

Nominal Simulation
The nominal simulation results of Scenario 1 (S1) indicate the feasibility of the proposed method.As Figure 2a shows, LOS rate .q finally converges to zero under the saturated control surface deflection.In Figure 2d, the control surface deflection of the IGC scheme is well bounded within [−10, 10]deg.The state estimate errors of NESO converge well in very small range, the NESO can estimate the system states with good accuracy, and all the states converge under saturated control surface deflection.The results indicate that the proposed control scheme performs very well.

Nominal Simulation
The nominal simulation results of Scenario 1 (S1) indicate the feasibility of the proposed method.
As Figure 2a shows, LOS rate q  finally converges to zero under the saturated control surface deflection.In Figure 2d, the control surface deflection of the IGC scheme is well bounded within [ ]   Figure 3 shows the comparison simulation results between S2 and S1.The proposed control scheme is able to cancel out influence of the nonlinear uncertainty in aerodynamics coefficients of lift and moment.Figure 4 shows the comparison simulation results between S3 and S1.The nonlinear parts of control surface deflection in aerodynamics coefficients have a limited impact on the convergence with small overshoot in Figure 4a,b.It is indicated that the performance improves under complex nonlinearity with the help of the great approximation ability of neural network.

Comparison Simulations
Scenario 2 (S2) and Scenario 3 (S3) are compared to test the performance of NESO-based IGC under large aerodynamics coefficient uncertainties.In S2, L α contains nonlinear items of AOA, while M α , M ω are also nonlinear functions of AOA and pitch angle rate.In S3, L α , M α are influenced by control surface deflection.Besides, simulation in S4 provides comprehensive testing under complex model uncertainties.
Figure 3 shows the comparison simulation results between S2 and S1.The proposed control scheme is able to cancel out influence of the nonlinear uncertainty in aerodynamics coefficients of lift and moment.Figure 4 shows the comparison simulation results between S3 and S1.The nonlinear parts of control surface deflection in aerodynamics coefficients have a limited impact on the convergence with small overshoot in Figure 4a,b.It is indicated that the performance improves under complex nonlinearity with the help of the great approximation ability of neural network.
under complex model uncertainties.
Figure 3 shows the comparison simulation results between S2 and S1.The proposed control scheme is able to cancel out influence of the nonlinear uncertainty in aerodynamics coefficients of lift and moment.Figure 4 shows the comparison simulation results between S3 and S1.The nonlinear parts of control surface deflection in aerodynamics coefficients have a limited impact on the convergence with small overshoot in Figure 4a,b.It is indicated that the performance improves under complex nonlinearity with the help of the great approximation ability of neural network.The comparison results of S3 and S4 are shown in Figure 5, which indicate that the proposed method can cancel out the complex uncertainties and disturbances with unknown measurement dynamics.The comparison results of S3 and S4 are shown in Figure 5, which indicate that the proposed method can cancel out the complex uncertainties and disturbances with unknown measurement dynamics.
(c) ( The comparison results of S3 and S4 are shown in Figure 5, which indicate that the proposed method can cancel out the complex uncertainties and disturbances with unknown measurement dynamics.In Figure 6, simulation results in S1 and S5 are compared.Figure 5d shows that the control surface deflection of S5 has 1 deg dead-zone.It is seen that the NESO approximation-based IGC control scheme guarantees that the closed-loop system is stable under saturated control surface deflection with dead-zone nonlinearity.In Figure 6, simulation results in S1 and S5 are compared.Figure 5d shows that the control surface deflection of S5 has 1 deg dead-zone.It is seen that the NESO approximation-based IGC control scheme guarantees that the closed-loop system is stable under saturated control surface deflection with dead-zone nonlinearity.

Simulations with Velocity, Aerodynamics Coefficients Bias and Measurement Noises
With the help of Monte Carlo theory, a large number of numerical simulations under large uncertainties of velocity, aerodynamics coefficients and measurement noises is made to test robustness of the proposed design.In Figure 7, the simulation results under velocity dispersion between 1000 m/s and 1400 m/s are presented.The mean value is 1200 m/s, and the red curves shows the corresponding results.It is seen that the parameter settings of the proposed controller can satisfy large variations of initial flight velocity.In Figure 6, simulation results in S1 and S5 are compared.Figure 5d shows that the control surface deflection of S5 has 1 deg dead-zone.It is seen that the NESO approximation-based IGC control scheme guarantees that the closed-loop system is stable under saturated control surface deflection with dead-zone nonlinearity.

Simulations with Velocity, Aerodynamics Coefficients Bias and Measurement Noises
With the help of Monte Carlo theory, a large number of numerical simulations under large uncertainties of velocity, aerodynamics coefficients and measurement noises is made to test robustness of the proposed design.In Figure 7, the simulation results under velocity dispersion between 1000 m/s and 1400 m/s are presented.The mean value is 1200 m/s, and the red curves shows

Conclusions
A novel composite IGC scheme in the presence of high-order nonlinear actuator dynamics is developed to address nonlinear hypersonic flight control with control input constraint and multiple uncertainties.The complex aerodynamics uncertainties and time-varying disturbances in different channels are well estimated by extended state observers improved by neural networks.Hybrid In Figure 9, the Monte Carlo simulation results under measurement noises with ±20% uncertainties in gain and time constant in Equation (40d) are presented.It can be concluded that the proposed controller has good robustness in the presence of large uncertainties of measurement dynamics.

Conclusions
A novel composite IGC scheme in the presence of high-order nonlinear actuator dynamics is developed to address nonlinear hypersonic flight control with control input constraint and multiple uncertainties.The complex aerodynamics uncertainties and time-varying disturbances in different channels are well estimated by extended state observers improved by neural networks.Hybrid

Conclusions
A novel composite IGC scheme in the presence of high-order nonlinear actuator dynamics is developed to address nonlinear hypersonic flight control with control input constraint and multiple uncertainties.The complex aerodynamics uncertainties and time-varying disturbances in different channels are well estimated by extended state observers improved by neural networks.Hybrid differentiators are employed to avoid directly differentiating the virtual control inputs in the back-stepping process.Thus, the peaking phenomenon and chasing of back-stepping sliding mode controller is greatly depressed.Nussbaum function is introduced to deal with the unknown direction of saturated control surface deflection.Series of numerical simulations indicate that the NESO-based IGC scheme can cancel out large uncertainties and disturbances.In conclusion, the proposed approach extends IGC application in hypersonic flight more practically and efficiently.
, which can be seen in Figure A1 and Table A1.

Figure 1 .
Figure 1.The planar geometry in terminal phase.

.
Then the state space form of integrated guidance and control model can be built as follows:

Figure 1 .
Figure 1.The planar geometry in terminal phase.

Figure 2 .
Figure 2. states and estimate errors of S1: (a) comparison between q  and NESOI, (b) comparison between α and NESOII, (c) comparison between z ω and NESOIII, and (d) saturated z δ curve.

Scenario 2 (
S2) and Scenario 3 (S3) are compared to test the performance of NESO-based IGC under large aerodynamics coefficient uncertainties.In S2, L α contains nonlinear items of AOA, while , M M α ω are also nonlinear functions of AOA and pitch angle rate.In S3, surface deflection.Besides, simulation in S4 provides comprehensive testing under complex model uncertainties.
δ z of S5 sat δ z of S1 zero offset

Figure 6 .Figure 7 .
Figure 6.Comparison between S1 and S5: (a) curves of .q, (b) curves of α, (c) curves of ω z , and (d) curves of saturated δ z .Energies 2018, 11, x FOR PEER REVIEW 14 of 18 the corresponding results.It is seen that the parameter settings of the proposed controller can satisfy large variations of initial flight velocity.

Figure 8
Figure 8 presents the Monte Carlo simulation results under aerodynamics coefficients bias with 20% ± uncertainties in , o L L α simultion of aerodynamics coefficients α (deg)

Figure 8
Figure 8 presents the Monte Carlo simulation results under aerodynamics coefficients bias with ±20% uncertainties in L α , L o of S3.The results indicate that the proposed IGC scheme is able to cancel out the great aerodynamics uncertainties in hypersonic flight.

Figure 7 .
Figure 7. Monte Carlo Simulation under V dispersion: (a) curves of q  , (b) curves of α , (c) curves of z ω , and (d) curves of saturated z δ .

Figure 8 Figure 8 .Figure 9 .
Figure 8 presents the Monte Carlo simulation results under aerodynamics coefficients bias with 20% ± uncertainties in , o L L α

Figure 8 .Figure 9 .
Figure 8. Monte Carlo simulation under aerodynamics coefficients bias: (a) curves of q  , (b) curves of α , (c) curves of z ω , and (d) curves of saturated z δ .In Figure9, the Monte Carlo simulation results under measurement noises with 20% ± uncertainties in gain and time constant in Equation (40d) are presented.It can be concluded that the proposed controller has good robustness in the presence of large uncertainties of measurement dynamics.
vector, and P = p, − k 2 e 2 2 − k 3 e 2 3 − k 4 e 2 4 + e 1 e 2 + e 2 e 3 + e 1 ε 1 +e 2 . The state estimate errors of NESO converge well in very small range, the NESO can estimate the system states with good accuracy, and all the states converge under saturated control surface deflection.The results indicate that the proposed control scheme performs very well.