Nonlinear Robust Control of a Hypersonic Flight Vehicle Using Fuzzy Disturbance Observer

This paper is concerned with a novel tracking controller design for a hypersonic flight vehicle in complex and volatile environment. The attitude control model is challengingly constructed withmultivariate uncertainties and external disturbances, such as structure dynamic and stochastic wind disturbance. In order to resist the influence of uncertainties and disturbances on the flight control system, nonlinear disturbance observer is introduced to estimate them. Moreover, for the sake of high accuracy and sensitivity, fuzzy theory is adopted to improve the performance of the nonlinear disturbance observer. After the total disturbance is eliminated by dynamic inversion method, a cascade system is obtained and then stabilized by a sliding-mode controller. Finally, simulation results show that the strong robust controller achieves excellent performance when the closed-loop control system is influenced by mass uncertainties and external disturbances.


Introduction
The study of the hypersonic flight vehicle (HFV) has been implemented for more than two decades all around the world.As early as 1980s, Chavez and Schmidt had made pivotal works of developing classic longitudinal model of the vehicle dynamics [1,2].Afterwards, Bolender and Doman developed the important implicit model: the first-principle model in [3].In recent years, flexible effect, which results in structure dynamic and may significantly influence the performance of HFV, attracts more and more attention.Some models introduce the flexible modes in the motion equations [4,5].The most significant characteristics of HFV are the highly coupled and nonlinear nature of its dynamic behavior.Additionally, complicated flight condition and variability of the vehicle characteristics bring about mass uncertainties for the modeling.Some researchers, such as McNamara and Friedmann, have stated the difficulties faced in the modeling and control in their works [6].At present, the models are unprecedented, complex, and much closer to the reality.
A large amount of different control methods has been successfully applied to stabilize the system of HFV and accomplish special tasks.Some typical approaches can be much weighted, such as linear parameter-varying control [7], adaptive control [8,9], sliding-mode control [10], gainscheduling, dynamic inversion [11], and H ∞ design [12].Some other control methods (e.g., linear matrix inequity (LMI)) are based on linearized model by using the small disturbance principle at given trim point [13].However, most of these methods have much conservation because the controller should be redesigned if the trim points of the system are altered and the calculation is amazingly tremendous.Thus, nonlinear control methods, such as dynamic inversion and adaptive control, seem more effective and present good relationship with physical process.
Despite of recent success of NASA's X-51A experiment, the design of robust guidance and control systems for HFV is still a challenging area [14].Due to the complex flight environment (e.g., high mach, low air density, and thermalelasticity), seeking for a strong robust controller is on edge.This kind of controller cannot only adapt to severe environment but also tolerate its own faults.As for HFV, one of the most important tasks is to track the reference trajectory.However, some unknown factors, such as gust wind, Mathematical Problems in Engineering structure dynamic, actuator failure, and sensor fault, bring about great uncertainties to the system performance.Some linearized models have been stabilized by robust controller based on LMI, in which Lyapunov stability theory is widely applied [15].In [16], energy norm indices have been used for the fault detection problem in discrete-time Markovian jump systems.Some other novel methods also successfully addressed the uncertainties and disturbance problems, such as in [17]; it used extended Kalman filter to solve the problems of stochastic nonlinearities and multiple missing measurements instead of using observing methods.Nevertheless, in [13,18], the parameters and states in linearized models are all the deviations from the trim point rather than the real ones.Thus, nonlinear analysis methods may be more meaningful in the practical engineering.The study of faulttolerant control, which is based on nonlinear observer, makes great contributions to the design of strong robust controller.There are lots of various state estimators at present.Such as in [19,20], state estimators were constructed with data from sensor networks for discrete-time systems to resist bad influence by sensor saturation, delay, and uncertainties.In [14], a fuzzy disturbance observer was proposed to estimate the uncertainties and disturbances.The nonlinear observer cannot only reconfigure the original states but also monitor the total disturbance.Then, the disturbance can be eliminated by using dynamic inverse method.As for this method, descriptions for disturbance and uncertainties are more accurate and rational for the reality.
In order to achieve excellent performance for the disturbance observer, in this paper, we introduce a novel fuzzy disturbance observer (FDO) to monitor the total disturbance, which can intelligently change its behavior according to the error of observation.What we contribute is that we redefine the disturbance as the sum of internal disturbance and external disturbance.This conception was first proposed in active disturbance rejection control (ADRC) by Han [21].This redefinition reduces the quantity of calculation compared to other online methods, such as adaptive control.Moreover, it becomes feasible to apply this technology to electronic programming [21].
In this paper, our control objective is a kind of HFV similar to High-altitude Test Vehicle-2 (HTV-2).The goal for control is to make the attitude states (angle of attack, angle of the sideslip, and angle of roll) track the reference trajectories as soon as possible, especially when the flight environment contains large scale of uncertainties and disturbances.The rest of the paper is arranged as follows.Section 2 demonstrates the attitude model of the flight vehicle with mass uncertainties in detail.And the disturbance observer based on fuzzy theory is introduced in Section 3. In Section 4, controller is designed based on the disturbance observer.Simulation results in Section 5 confirm the controller's excellent operation and how well the observer monitors the disturbance in real time.

Nonlinear Attitude Control Model
And the angle-rate equations are given as where   ,   ,   ,   ,   ,   ,   ,   , and   are related to one or four parameters (their expressions are listed below) and can be gotten by interpolation method: (, ,   ,   ) ( = , , , , , ) ,   () ( = , , ) . (3) Figure 1 shows the moment and force curves and their fitted functions versus  and , respectively.As the coefficients are determined by four variables: , ,   , and   .The inputs   and   are both frozen to 0.1 degree to draw the threedimensional pictures below.
Therefore, ( 1) and ( 2) can be rewritten to (4) by separation of variables: where Equation ( 4) is a typical system that can be stabilized by back-stepping method, in which Ω can be seen as a pseudoinput.In this system,   ,   , and   are all determined by the reference trajectory and known for controller design.

Problem Formulation.
As for describing the complex and volatile flight condition, a great amount of parameter uncertainties and stochastic wind disturbance is introduced into this model.The parameter uncertainties and external disturbance (bandlimited and varying slowly) are expressed in the equations as follows: In system (6), f 1 and f 2 contain lots of complicated structure information of this model, which are too complicated if they are written into formula expression.For simplicity, some items are rewritten as the following expression: thus, system (6) can be rewritten as where In particular, total disturbance is redefined as unknown and bandlimited, to represent the sum of partial structure information, model parameter uncertainties, and wind disturbance.From the above consideration, a cascade system with total disturbance is presented.As for the practical hypothesis, the total disturbance, hence, contains wide range of parameter uncertainties, such as force and moment coefficients deviation, inertia deviation, gravity center deviation, and atmosphere density deviation.As the goal for flight control, the goal is to converge the flight attitude and angle-rate (Θ(), Ω()) into the given reference trajectory (Θ * (), Ω * ()) in finite time.And the transient response performance should be as better as the system can achieve.Consider the following: Dynamic inversion method has been widely used to address such problem; however, the most critical and difficult work is the disturbance estimation.Disturbance estimation relies on identification via inputs and outputs data.As long as the total disturbance is estimated (notated as Π ( = 1, 2)), the problem surely becomes easy to disentangle by dynamic inversion method as the following expression; Ω behaves as a pseudoinput: Ω  .Consider the following: In (11), the inversion matrix B −1 cannot be obtained because B is not a square matrix.Thus, pseudoinversion B † is used to replace it, and the error which results from inaccurate dynamic inversion can be reduced by error-integral feedback.If Π can be perfectly monitored, the plant system (4) can be rewritten as, by replacing with (11), Assumption: the observation error ( Π −Π  ) is minor, and its effect on the stability of the closed loop system is negligible.This assumption is reasonable because it has been proved  that the observation error is uniformly ultimately bounded (UUB, ∈ L ∞ ) within a region of arbitrarily small size [14].Thus, we will neglect the observation error in the following parts.
The estimation of Π ( = 1, 2) plays an extremely important role in the whole control system.The total disturbance should be estimated accurately online, and fuzzy theory is quite suitable for this condition.After choosing the membership functions, the fuzzy disturbance observer will have different response time and convergence speed in accordance with which membership area the observation error is located in.This kind of intelligent observer will make the online identification process achieve excellent performance.

Disturbance Observer Based on
Fuzzy Theory  the system (8), the following NESO is designed to estimate the total disturbance: where In system (13), the nonlinear function G(Γ, u) is a fuzzy inference engine, which is used to estimate the total disturbance.Π is the extended state, which is defined in the subtitle.The G(Γ, u) can be also rewritten as G(Γ) because the input u(Γ) is determined by full-state feedback.The constant vector c, if chosen properly, will control the weight the reconfigured error has in the observer equation and it will enhance the stability of the whole system by negative feedback.The goal for this observer is to make the following equation hold: Figure 2 depicts the whole controlled system based on fuzzy disturbance observer (FDO).

Disturbance Observer Based on Fuzzy Theory. In [14],
Euntai gave detailed explanations on fuzzy disturbance observer.A typical fuzzy system follows the IF-THEN rules which are characterized by a set of condition → action rules, or in IF-THEN form as follows [22]: or IF  1 is   1 and, . . ., and   is    THEN   is .
In this paper, the principle of product inference, centeraverage and singleton fuzzifier is adopted.Then, the expression of the estimated total disturbance is obtained as follows: where θ ∈ R 27 × 1 , (Γ) ∈ R 27 × 1 .θ is the adjustable vector, which is composed of   , and (Γ) is the fuzzy basis function.Equation ( 18) is a general solution formula for every channel of the estimated disturbance vector Π.The adaptive parameter is shown as follows (take the channel of , e.g.): and the fuzzy basis functions are as follows: where . .,  27 ]  , and For the sake of reducing calculation amount, there are only three fuzzy membership functions for each channel of , , and , which are shown in Figure 3.
The disturbance observation error ΔΓ will be uniformly ultimately bounded (ΔΓ ∈  ∞ ), which has been proved in [14].And this method will accelerate the convergence speed of the FDO when the observation error is approaching zero.As for (19), the fuzzy basis functions (, , ) will be the same when dealing with the channels of , , and .But when estimating the total disturbance of the channels of   ,   , and   , the functions will change to (  ,   ,   ), and the method of obtaining it is similar to the former.

Controller Design Based on Sliding-Mode Method
After observing and removing the total disturbance, the cascade system (12) is obtained.In this paper, the slidingmode method is used to deal with the cascade system.Angle tracking errors are defined as follows: where Θ * is the given reference trajectory.Differentiating (22), where slaw  is the sliding-mode approaching law.In order to achieve good tracking performance, in this paper, we choose the exponential approaching law instead of the traditional proportional approaching law: In ( 24), the sliding-mode surface function s  satisfies threes conditions: (1) s  ( =  0 ) = 0, (2) s  is differentiable, and (3) s   ṡ  < 0. These conditions guarantee not only the stability of tracking process but also the steady tracking error  asymptotically approaching to zero.Substituting (25) into (24), we get the controlled input as Remark.According to the repeated trials, it is more likely to achieve good control performance when the parameters are chosen as  1 = 5,  1 = 10,  11 = 3, and  12 = 1.All of these parameters can significantly influence the stability, the rising time and the tracking error of this closed-loop system but no single parameter can determine any of them independently.Thus, so far the adjusting method is one of the best methods to determine these parameters.Substituting (26) into (11), we obtain the final controller expression with the negative feedback of the integral of the tracking error: where

Simulations
In this part, a wide range of parameter uncertainties and external disturbance is set for the plant system, and parameter uncertainties are realized by the distractions of relevant parameters in the model of plant system.The goal is to test the ability of the controlled system to resist disturbance and its own uncertainties.It is also expected that the disturbance observer will achieve good performance.

Parameters Distractions and External Wind Disturbance.
In Table 1, all of the distracted parameters are listed, in which numbers expressed in percentage mean relative distractions, while others indicate absolute distractions.All of these changes occur in the model of plant system, which represents the vehicle in the virtual reality environment for flight.
Then, the model for wind disturbance is given as follows.In general, both stationary wind and gust wind will not be taken into consideration until the simulation runs to 60 seconds.The speed of stationary wind will change according to the flight altitude, and the accurate value for this speed can be gotten by linear interpolation.Table 2 demonstrates the relationship between the wind speed and the flight altitude.
As for the gust wind model, a regular model with the period of 50 km is considered, which is illustrated in Figure 4.For each time, the gust wind starts at  0 and ends at  3 .After Mathematical Problems in Engineering 9 50 km flight, which is at  0 + 50 km, the same wind reappears again.
The directions of both stationary wind and gust wind are stochastic, and   denotes the azimuth angle of wind.Thus, the wind vector is as follows: where The wind vector mainly influences the system outputs ( and ) by changing the speed and direction of the incoming flow.So far, all the modeling work for total disturbance has been finished.In the next subsection, several pictures of simulation results will be shown.

Simulation Results.
The main goal for this paper is to make the HFV track reference trajectory with high performance in severe situation.At first, it only tests the performance of the controller based on sliding-mode method; thus, total disturbance is not added into the plant system.Simulation results are shown in Figure 5(a).It can be clearly learned that the controller has performed an excellent work for the tracking task.Although the tracking performance is not good at about 50 seconds because the angle of attack suddenly drops which results in the rapid declination of lift force, the whole system soon goes back to its stable states under the action of controller and the error is limited within two degrees.When the total disturbance is added into the system, the performance of the controlled system is shown in Figure 5(b).From these figures, it demonstrates that even though the disturbance exerts negative impact (marked in the figures) on the whole system, it has successfully finished the tracking task and left minor steady error.All these accomplishments benefit from the excellent performance of FDO. Figure 6 demonstrates that the total disturbance is superbly monitored by the FDO.In Figure 6(a), it clearly shows that the periodic gust wind disturbance is totally observed, which is in accordance with the real situation.

Conclusions
This study indicates that, although the flight condition for the hypersonic flight vehicle contains mass uncertainties and gust wind disturbance, the fuzzy disturbance observer monitors the total disturbance very well.The definition of total disturbance provides a new way of disentangling complex structure information.This method contains both the analytical dynamic inversion and the numerical online identification, which contributes to its excellent performance when the tracking trajectory varies drastically.In the future, this work is expected to focus on the study of the estimated total disturbance separation and analysis.

Figure 2 :
Figure 2: Controlled system based on FDO.
, , and  performance in tracking references without total disturbance , , and  performance in tracking references with total disturbance

Figure 5 :
Figure 5: Tracking results without and with total disturbance.
Estimated disturbance of angle-rate channels

Table 1 :
Distractions of parameters in the plant system model.