Nonlinear Trajectory Controller with Improved Performances for Waveriders

1 e Research Institute of Pilotless Aircra , Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China 2Key Laboratory of Unmanned Aerial Vehicle Technology of Ministry of Industry and Information Technology, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China 3 e College of Astronautics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China


Introduction
Air-breathing waveriders are considered a possible means to access space, and these vehicles have commercial and military implications.Unlike rocket-based (combined cycle) systems, air-breathing waveriders can reach orbital speeds without carrying oxygen by applying the scramjet propulsion system [1].However, waveriders are characterized by extreme aerothermo-elastic-propulsion interactions and uncertainty.As a result, the design of satisfactory control systems for air-breathing waveriders becomes a critical issue to address, with strong couplings between aerodynamic and propulsive effects, while addressing significant uncertainties associated with the complicated flight conditions and special geometries required for these vehicles [2].
For the design of guidance and control systems of airbreathing waveriders, an important task is to construct dynamical models.To this end, a comprehensive analytical model of waverider was developed in [3] using Newtonian impact theory to estimate the pressures that act on the vehicle.Based on the integration of computational fluid dynamics (CFD) and analytical techniques, the stability and control derivatives were determined in [4] to build a model of the longitudinal dynamics of the waverider and investigate the propulsion/airframe integration.A first principlebased dynamic model of the waverider was provided in [5] to incorporate aero-thermo-elastic-propulsion interactions: this model applied inviscid compressible oblique shock-expansion theory to estimate aerodynamic forces and moments, as well as a 1D Rayleigh flow model and an Euler-Bernoulli beam model determined scramjet propulsion and structural flexibility, respectively.Given that the nonlinear physics-based dynamical model in [5] was highly complicated to be employed for control analysis and design, a simplified model of this version was given in [6] to remove a set of weak couplings of nonlinear dynamics.For a comprehensive aero-thermo-elastic-propulsion model of the six-degree-offreedom dynamics of waverider in [7], a control-oriented model was developed in [8] to allow stability, controllability, and robustness analyses and support the adaptive and Complexity nonlinear control law design.Based on these control-related models of the waverider, some nonlinear controllers have been developed, including the adaptive sliding mode control [9] and linear output feedback tracking control [10].
The control system design of waveriders should be integrated with the guidance law because any small attitude variation leads to a significant change in the trajectory on the hypersonic flight condition.Therefore, feasible trajectory generation is critical for the design of guidance and control of waveriders to improve the flight performance and minimize the control power [11].The periodic cruise can typically save more fuel than the steady-state cruise in the hypersonic flight phase [12]; thus, applying the periodic flight trajectory helps improve the hypersonic flight performance.Nevertheless, tracking such a trajectory is difficult for the guidance and control system while guaranteeing the respected control qualities, such as fast setting time, small overshot, and robustness during the hypersonic flight.To this end, Parker et al. [13] designed a control system of an air-breathing waverider based on approximate feedback linearization to achieve excellent tracking performance and robustness.Alternatively, a nonlinear adaptive dynamic inversion controller was developed in [14] for a waverider to achieve the desired reference track while being robust to the system uncertainties.The dynamic inversion-based flight control, however, has difficulty in ensuring robustness and good transient performance simultaneously for waveriders [15].
A so-called composite nonlinear feedback (CNF) control technique, presented by Lin et al. [16], can ameliorate transient performances in the control process.Turner et al. [17] applied CNF to higher-order and multiple-input systems.Chen et al. [18] developed a CNF control to a general class of systems with input saturation.This controller consisted of a linear feedback control law and a nonlinear feedback control law without any switching element; as a result, the CNF design can capture the time-optimal maneuver in asymptotically tracking situations [19].The CNF technique has also been used for nonlinear systems [20], uncertain chaotic systems [21], discrete-time linear systems [22], master/slave synchronization of nonlinear systems [23], and switched systems with input saturation [24].In particular, a design method of composite nonlinear feedback control technique was presented in [23] for the synchronization of master/slave nonlinear systems, and a new condition was derived for the master and slave systems.Also, a robust tracking problem was considered in [24] for the uncertain switched systems with input saturations, and a new form of the nonlinear function was generated based on a CNF control law to adapt the changes in the tracking targets.Furthermore, a variant factor technique was proposed in [25] to address input saturation in a class of nonlinear systems.If the flight control law is designed with the dynamic inversion and CNF technique, some important control qualities can be improved for waveriders, such as the transient performance and decoupling control capability.Accordingly, this study develops a novel flight control law for a general model of waverider using the CNF technique to guarantee integrated control performances, such as fast setting time, small overshot, and robustness on the hypersonic flight condition.Not only that, the periodic flight trajectory can be followed using this designed controller under uncertain flight conditions, and optimal energy consumption is guaranteed in the course of the waverider flight accordingly.
The rest of this paper is divided into several sections.Section 2 deals with the model properties and periodic cruise trajectory generation of the waverider.Section 3 considers the guidance and control law design for the periodic trajectory using feedback linearization theory and CNF technique.The proof of system stability and robustness is provided accordingly.Section 4 demonstrates an illustrative example, which tests the feasibility of the proposed controller for a general longitudinal model of the waverider.The response performances and fuel consumption savings are compared in relation to the dynamic inversion control and CNF control.Section 5 concludes.

Nonlinear Controller with Improved Performances for MIMO Model
A multi-input multioutput (MIMO) affine nonlinear model of waveriders is written with input saturation as follows [26]: where  ∈ R 7 ,  ∈ R 2 , and   ∈ R,  = 1, 2 are the state, output, and control variables, respectively. is a function defined as follows: for  = 1, 2, where   is the maximum of   .Following that, we give the following assumptions [27].
(2) The relative degree of the nonlinear model is The zero dynamics of the nonlinear model is stable on  if any.
Based on the differential geometry theories, for the nonlinear model there is a nonsingular transformation Φ() from  to × 0 which is a compact and connected set of R 7 [28], where for  = 1, ⋅ ⋅ ⋅ , 7 and  = 1, 2, for  = 1, 2, for  = +1, ⋅ ⋅ ⋅ , 7 and  = 1, 2, and   () is a smooth function on  for  =  + 1, ⋅ ⋅ ⋅ , 7. According to the nonlinear state transformation Φ(), the nonlinear model can be changed as follows: where  0 ∈ R 7− is the state variable of the zero dynamics and  ∈ R  .In addition,   (,  0 ) ∈ R (,  0 ) = (  () where Φ −1 represents the inverse of Φ, which is defined as a nonsingular function from  to  ×  0 .Accordingly, the system matrices are for  = 1, 2, where   ∈ R   ×  ,   ∈ R   ×1 ,   ∈ R 1×  , and  denotes the identify matrix with approximate dimensions.In addition, a reference generator needs to be designed for the waverider to produce reference in (69) to be tracked, and, based on the equivalent model in (8), it is constructed as follows: where  ∈ R  is the state of the reference generator,   denotes the feedback gain matrix,   represents the virtual signal source, and   indicates the initial value.This reference generator can produce an arbitrary type of output signal, including the sinusoidal signal and the ramp signal by choosing   ,   , and   .
Furthermore, we will present a CNF control law associated with the reference generator given in the above section in order that the proposed CNF control law can track the designed reference in (15).To this end, we first define   =  − , and then, based on ( 8) and ( 15), we get This error equation is applied in the design of the CNF control law, and the detailed design procedure will be introduced in the following steps.
Steps .Design a linear feedback control law where  is the chosen state feedback gain matrix R 2× such that (1)  +  is an asymptotically stable matrix and (2) the closed-loop system ( −  − ) −1  has the expected properties.The selection of  is to make the closed-loop poles of  +  have a dominating pair with a small damping ratio, leading to a fast rise time in the response process.
Steps .Given a positive definite symmetric matrix  ∈ R × , we solve the following Lyapunov equation: for  > 0. Accordingly, the nonlinear feedback portion of the CNF control law   is provided by where  = diag( 1 ,  2 ), and   ,  = 1, 2, represent the nonpositive functions regarding   , which are used to gradually adjust the damping ratio of the closed-loop system to obtain better tracking performance.
Steps .The linear feedback control law and nonlinear feedback part are constituted to get a generalized CNF control law: Steps .The control law for the equivalent linear model is designed as follows: Furthermore, for a set of scalars   ∈ (0, 1),  = 1, 2, let  > 0 be the largest positive scalar such that, for all  ∈   , where the following condition holds: where   is not zero as  is invertible, and    ∈ R 2 is a vector in which only the i-th element of   is 1 and the other is zero.
Beyond this, for getting further results, first we propose a lemma below.
Afterwards, we provide a Lagrange function as follows: where 2  is the Lagrange multiplier.Following that Actually, if   ̸ =   , there is no stationary point.As   =   , there is one stationary point,    = 1, and      −1    =     so that Thus, This completes proof of Lemma 1, and the result of Lemma 1 will be applied in the proof of the following Theorem.
(3) Zero dynamics of the nonlinear model in ( ) is stable on is nonsingular on .en, with the CNF control law comprising ( ) and ( ) will drive system output  to track arbitrary reference  from an initial state asymptotically without steady-state error, provided that the following properties are satisfied.
(3) Let   be the initial value of .en, the initial condition   of   satisfies Proof.Consider   =  − ; then When   ∈   , it shows that Thus, Afterwards, with ( 32) and ( 37), we have Following that, the input channels with the saturation can be expressed by [21]: where  =    −1     ,   = diag( 1 ,  2 ), and   ∈ [0, 1] for  = 1, 2. Substituting (39) into ( 16), we have As a result, the closed-loop system involving ( 16) and (40) can be expressed by The above equation is combined with two parts and the second part is assumed to be stable for the zero dynamics, so only stability of the first part needs to be proved accordingly.
Remark. the equivalent form of ( 41) is expressed as where   =  +  and   = .Based on the results in [28], the stability of the closed-loop system matrix   in the region of linear matrix inequality (LMI) is satisfied with where the symmetric matrix ,  and matrix  make Ω = { ∈  : L + Q + Q  < 0}.Accordingly, Ω indicates a LMI region.Therefore, the adjusted function  can ensure the system asymptotic stability to meet the stability condition of Linear Matrix Inequalities (LMIs) accordingly.To be specific, the main task of adding the nonlinear part is to alter the damping ratios of the closed-loop system as the outputs approach the periodic cruise trajectory commands, and the nonlinear gains can help to improve stability robustness and to suppress the effect of uncertainties in accordance with the adaptive regulation ability.

Nonlinear Trajectory Controller for Waveriders
This study investigates the design methods of the track controller using CNF technique as the waverider flies along with the periodic cruise trajectory.The vehicle shape chosen for this study is shown in Figure 1, where the vehicle geometry is assumed to have unity depth [5]. Figure 1 illustrates the basic geometry of the waverider.A combination of oblique shock and Prandtl-Meyer flow theory is used to determine the pressures on the vehicle surfaces, whereas a 1D Rayleigh flow scramjet propulsion model is applied with a variable geometry inlet.Specifically, an oblique shock occurs on the lower forebody as  ≥ − 1, .The respective airflow properties across the bow shock are determined by [29]: where   is the specific heat ratio for air,   denotes the shock wave angle,  indicates the angle of attack,  ∞ and  1 represent the flight Mach on the front and rear of the shock wave, respectively, and  ∞ and  1 indicate the pressures on the front and rear of the shock wave, respectively.Similarly, a Prandtl-Meyer expansion occurs on the upper surface as  ≤ − 1, , and the relations between the geometric parameters and aerodynamic characteristics are decided by [30]: Once the pressures acting on all vehicle surfaces are computed based on either oblique shock theory or Prandtl-Meyer expansion flow theory, the aerodynamic forces such as the lift , drag , and pitch moment   can be determined.Aside from these aerodynamic forces, the propulsive force  is calculated according to the momentum theorem from fluid mechanics and expressed by [31] where   and   indicate the exit area and inlet area of the propulsive system, respectively,   and   denote the exhaust exit pressure and velocity, respectively, and ṁ is the mass flow of air, which is determined by [32]: where ℎ 2 and ℎ 3 denote the total enthalpy at the combustor inlet and exit, respectively,  and ṁ are the equivalence ratio and mass flow of fuel, respectively,   represents the combustor efficiency, and   and   indicate the stoichiometric fuel-air ratio and fuel lower heating value for the given fuel, respectively.
When the lift , drag , pitch moment   , and thrust  are acquired, the longitudinal model of the waverider along with the velocity coordinate system is provided by [33]: where  denotes the flight velocity,  is the flight path angle,  represents the pitch angle velocity,  indicates the angle of attack, ℎ is the flight height, and  denotes the flight range, which is determined by ṙ = cos.In addition, , , and   are the airplane mass, gravitational constant, and moment of inertia, respectively.For (60), the elevon deflection angle   and equivalence ratio Φ are chosen as the control inputs, and then, based on the aforementioned force computation theories, the lift , drag , pitch moment   , and thrust  can be expressed as functions regarding the flight states and inputs by where  and  indicate the reference area and mean aerodynamic chord, respectively,  represents air density, and   ,   , and   denote the coefficients with regard to the lift, drag, and pitch moment, respectively.We consider the classical short-period approximation by applying the Jacobian linearization for (60) and (61).The transfer function between the flight path angle and elevator deflection is approximately described as [33]: where where   = (1/  )(  /),   = (1/  )(  /),    = (1/  )(  /  ),   = −(1/)(/), and    = −(1/)(/  ).In addition,   represents the transfer function gain, and the subscript  denotes the trim value.
For the basic geometry of waverider in Figure 1,   > 0 because of the large force acted on the lower surface, which presents an unstable pole according to (63), thereby leading to unstable short-period dynamics for the longitudinal model of the waverider [5].In addition, the lift corresponding to the elevator deflection results in the existence of a nonminimum phase zero because of    < 0. Fortunately,    is very small because of the flight in the very low density environment, such that this nonminimum phase zero locates at the far right.In sum, the model dynamics of waverider are unstable, nonminimum, and strong coupling and present large uncertain disturbances and unknown model dynamics.As a result, designing a satisfactory guidance and control system for waverider is a challenging task.
In this section, a guidance and control law using the dynamic inversion and CNF technique is derived for the waverider to track the optimized periodic cruise trajectory.This proposed design process integrates the guidance and control loops by merging the attitude, velocity, and altitude dynamics of the waverider.
The established nonlinear model of waverider in (60) and ( 61) is difficult to use for some nonlinear control techniques because the obtained lift, drag, and moment are nonanalytical.As a result, the control law design must conduct the model simplification associated with the design requirements.The critical dynamic characteristics concerning the initial model should be retained, although it is simplified to be analytically tractable.Thus, the core work in the model simplification course aims to identify the analytical expressions of the lift, drag, thrust, and pitching moment that need to not only conform to the initial model features but also be in accordance with the control-oriented design goal.
In this study, the forms for the force and moment functions in (61) are selected as [6]:  denote the needed identification parameters for the lift terms.The item parameters for the drag, pitching moment, and thrust also have to be acquired.The scramjet of the waverider strongly depends on the angle of attack and dynamic pressure along with the input Φ.The flight attitude is also dramatically influenced by the propulsive action because of the airflow compression and expansion at the inlet and outlet, which leads to the significant change in the lift and drag.As a result, the mutual effects between the flight attitude change and propulsive action bring a significant challenge for the design of the guidance and control system.
In further implementation, the parameters of the polynomial functions in (65) and (66) are acquired using the leastsquares approach based on the model datum derived from the anterior estimation approaches concerning the aerodynamic and propulsive forces.Theoretically, these parameters that have considerable coupling items may result in fairly accurate fitting results, but the complexity of these functions is unfavorable for the guidance and control law design.Thus, from the perspective of the control design for the waverider, we need a compromise selection of the function forms in (65) and (66).

Complexity
To design the integrated guidance and control law in the preceding part, the outputs of the nonlinear model in (60), (61), and (65) are selected as  1 =  and  2 = ℎ.In addition, a second-order actuator model is added to the input , and the propulsive commanded value   is selected as the new control input; that is, where   and   represent the damping ratio and natural frequency of the actuator, respectively.Thus, the control inputs are changed as  1 =   and  2 =   .Furthermore, we choose the flight state vector as  = [, , ℎ, , ,   , β  ]  such that the nonlinear model in combination with (60) and ( 67) becomes a 7D model; that is, this model has a relative degree of 7 over the operating range of interest.
The application of periodic cruise trajectories for hypersonic flight can achieve better fuel consumption savings than that of the steady-state cruise trajectories because of the more effective kinetic-potential energy interchange of the former [14].However, designing and tracking such periodic cruise trajectories are more difficult and complicated than the process for the steady-state trajectories.Therefore, planning and optimizing the feasible periodic trajectories are critical to improve the cruise performances for waverider.
For the steady-state trajectories, once the model properties, environmental parameters, and aerodynamic expressions are determined, the corresponding trim values with regard to the given velocity and altitude can be solved by For the periodic cruise trajectories, how to acquire considerable performance benefit needs to be discussed in view of fuel consumption savings.Compared with the steadystate trajectories, where the engine is always switched on, the periodic cruise trajectories allow the waverider to switch off its engine at some points.For simplicity, we consider that the engine works as  > 0 when the vehicle continues to climb.In turn, the engine turns off in the descent phase.We further adopt a parameterized periodic trajectory in [34], which is given by where ℎ  is the amplitude of the periodic trajectory,  denotes the damping term, and   indicates the period parameter regarding the flight range.ℎ  , , and   are the adjusted parameters that should be optimized to minimize the following cost function [14]: where   is the total flight time.This cost function minimizes the ratio of the fuel consumption over the flight range.
According to the parameterized periodic trajectory in (69), the dynamic balance values can be computed by where Based on (74), (75), and (76), the dynamic trim states can be obtained in relation to each flight speed and altitude.If the periodic trajectory is given, then the reference command associated with the guidance and control action is decided.The optimal periodic cruise problem in (73) is too complex to be solved analytically; thus, it is solved numerically in this study using the genetic algorithm toolbox in MATLAB.
Once the optimal values are obtained, we design the guidance and control law for the waverider to follow this optimized reference trajectory.Furthermore, the nonlinear controller in Section 1 is employed to the model-based optimized reference trajectory of the waverider, and the structure diagram of the trajectory controller using CNF technique is provided for the waverider in Figure 2.
Figure 2 shows that the presented controller combines the CNF technique with the inversion control structure for  the waverider.To be specific, the nonlinear model is first converted into the resulting equivalent model based on the differential geometry theories.Then, a linear feedback control law is designed to become a closed-loop system with faster rise time.Afterwards, the nonlinear feedback control law is applied to reduce overshoot and improve transient performance in tracking control.Not only that, a combination of the linear and nonlinear feedback control will guarantee closed system stability and robustness in the presence of the unknown model dynamics, while repressing the uncertain disturbances in the complicated flight environment.
The selection of the design parameters  and  in ( 18) and ( 19) is important to improve transient performance in tracking control, and this is because the nonlinear part of CNF controller will assist speeding up the settling time and decreasing the overshoot, or, equivalently, contributing to a large part to the control input as the tracking errors are small.Basically, the poles of the closed-loop system can be adjusted by the functions   as a result that these poles approach the invariant zeros of the auxiliary system which is similarly defined in [20] as |  | become larger and large.In fact, the larger |  | will lead to a larger damping ratio and yield a smaller overshoot, so one possible choice of   is provided as follows: where   and   are chosen positive scalars to yield a desired track performance.Evidently, when the track error   is large, the effect of the nonlinear part of the CNF controller is limited due to the small value of   , whereas the nonlinear part will be effective as   is small.Finally, we note that the choice of   is nonunique, and any function will work as if it has similar characteristics of that given in (77).
Besides the given guidance and control law design of waverider, the CNF control methods are also considered a general strategy for flight control.Correspondingly, the flow diagram with regard to the design steps using the CNF technique is provided in Figure 3.
According to Figure 3, the general design steps for the flight control law design are provided as follows.
Steps .The database of the aerodynamic forces and thrust is built based on the CFD and engineering estimation methods.The obtained aerodynamic forces and thrust are input into the nonlinear airplane model.A full simulation model derived from first principles is established to verify the presented control designs.
Steps .Based on the acquired databases, the polynomial expressions regarding aerodynamic forces and thrust are identified using fitting approaches.A curve-fitted model is accordingly obtained.The inherent dynamics of this model with the fitting expressions is compared with the full simulation model to demonstrate the effectiveness of the model simplification mean.
Steps .After determining the simplified model with the polynomial expressions, we obtain the resulting equivalent model using differential geometry theories.Unlike the Jacobian linear structure with the given flight point, the feedback linearization method can maintain the nonlinear model dynamics, broaden the resulting structure over a large Complexity flight range, and be combined with some advanced control methods.
Steps .The guidance and control law with linear feedback and a nonlinear control part is designed using the CNF technique and inversion control structure.This proposed guidance and control system can track general target reference trajectories with input saturation.The respected flight qualities, including fast setting time, small overshot, and robustness, are ensure.
In general, the developed guidance and control law using the CNF technique is generally suitable for the nonlinear model of any fixed-wing vehicle to track arbitrary flight trajectories in line with the aforementioned step-by-step design procedure.This process depends on the controloriented model with the curve-fitted mean instead of the linear model with the Jacobian linearization while driving the applicability of the CNF technique in combination with the dynamic inversion control.As a result, this proposed controller can ameliorate transient performance in tracking the anticipated reference trajectories and guarantee system robustness.

Illustrative Example
This study adopts the typical waverider geometry of waverider in [5] to validate the feasibility of the proposed control law.This configuration makes the vehicle forebody become a significant section of the inlet compression process, whereas its afterbody performs a large part of the nozzle, which generates considerable thrust [35].The propulsive system design encompasses the entire undersurface of the waverider.As a result, the engine and airframe become one [36].In a real application, the periodic flight trajectory can be designed to widen the cruise range and to save the consumed energy [12].For this reason, an illustrative example is provided for this study in a similar manner to the actual situation for waveriders [37].
Aerodynamic data are obtained for the geometry of the waverider in Figure 1 using a combination of oblique shock and Prandtl-Meyer flow theory in (51)-(56).The thrust is estimated based on the 1D Rayleigh flow scramjet propulsion model in (57)-(59).For simplicity, some shape parameters are suitably changeable in this study in comparison with those in [5] to ensure the rationality of the proposed model from the design viewpoint of flight control.When these model parameters are known, the polynomial expressions concerning these forces and moment can be identified.Additionally, we select the operation condition as  = 8, and ℎ = 27.The resulting aerodynamic and propulsive expressions are approximately acquired according to (65), and they are expressed below: Furthermore, the trim states in the vicinity of the operation condition based on (68) are plotted in Figure 4.
Figure 4 shows that the trim values dramatically change with the different flight Mach numbers and altitudes near the operation point, which demonstrates that the inherent dynamics of the waverider is sensitive to the change in the flight condition.Moreover, the flight dynamics is also unconventional due to the presence of the zero and pole in the right plane regarding (62), as displayed in Figure 5.
Figure 5 shows the zeros and poles on the right plane as the flight Mach number is changed from 7.8 to 8.2 and the altitude is changed from 26.5 km to 27.5 km.Obviously, these zeros and poles indicate that the nonlinear model of the waverider is unstable and a nonminimum phase as a result of the limitation for the control law design, such as the control bandwidth and effective control region.Furthermore, we select  = 0.5, ℎ  = 27000, and ℎ  = 300 in (69).The consumed fuel with the changes in the adjusted parameter   and flight range  is plotted in Figure 6.
According to (73) and (74), we can obtain the optimal adjusted parameter   = 62931 by using the genetic algorithm.Thus, the reference command of the optimal periodic cruise in (73) is acquired as The dynamic balance states can be solved by substituting (79) into (74).The results with the change in the flight range are provided in Figure 7.
Figure 7 demonstrates that the flight path angle, angle of attack, and control inputs change periodically to follow the expected flight trajectory in (79).Also, we consider the velocity change command with Δ  = 50/ from  = 8 in the control process.In this case,  1 and  2 are selected as the unit matrix, and  = diag( 1 ,  2 ) is given as The nonlinear gain functions  1 and  2 are selected as In the simulation, the response results of the structures with CNF and without CNF are considered for the waverider model, passing through 200s.The results in relation to the steady-state flight and periodic cruise are presented in Figure 8.          the nonlinear part of CNF makes the system outputs follow the periodic altitude command and step velocity command well while guaranteeing the control performances including small overshoot and fast settling time.In addition, the angle of attack is smaller using the CNF controller such that the propulsive force can keep stable and the drag can decrease accordingly.As a result, the consumed energy will reduce for the hypersonic vehicle due to the nonlinear part of CNF, and the transient performances will be improved for the periodic cruise trajectory of the hypersonic vehicle.
Figure 11 shows that the fuel consumed when flying with the periodic cruise trajectory is less than that flying with the steady-state cruise trajectory, whereas the consumed fuel using the CNF controller is minimal throughout the entire flight process.These results indicate that the flight control law using the CNF technique can not only improve the track performances but also achieve high fuel consumption savings for the waverider.

Conclusion
This study proposes a design method of the flight control law using the CNF technique to track the periodic cruise trajectory of the waverider.Compared with the well-known methods, the advantage of the proposed controller can guarantee the respected control qualities, including fast setting time, small overshot, and robustness over the hypersonic flight range.Furthermore, the simulation results show that the proposed control law improves the flight performance with input saturation.Fast setting time, small overshot, and robustness are ensured for the waverider while considering the tracking issue of the periodic cruise trajectory.However, we have assumed that all the flight states are available to feedback.Thus, future works should introduce the state observers for this presented control law so that these flight states can be estimated to satisfy the requirements in real applications.

Figure 3 :
Figure 3: General design steps for flight control law Using CNF technique.

Figure 5 :
Figure 5: Change surfaces of zero and pole on the right plane with regard to  = 7.8 − 8.2 and ℎ = 26.5− 27.5.

Figure 8
Figure8shows the contrast curves of the structures with CNF and without CNF.In consideration of the nonlinear part of the CNF control law in(21), the velocity and altitude output can rapidly follow the commands signals in contrast to that using only inversion control.

Figures 9
Figures 9 and 10 indicate the change curves of the angle of attack, flight path angle, and control inputs that gently return to the anticipated balance values due to the control action.Compared with the nonlinear control law provided in[38], the proposed controller which depends on

Figure 6 :
Figure 6: Consumed fuel with changes of adjusted parameter   and flight range .

Figure 8 :
Figure 8: Response curves in relation to level flight and periodic cruise between with CNF and without CNF.

Figure 9 :Figure 10 :
Figure 9: Angle of attack and flight path angle between with CNF and without CNF.
Structure diagram of guidance and control design using CNF technique for hypersonic vehicle.