Robust Controller Designing for an Air-Breathing Hypersonic Vehicle with an HOSVD-Based LPV Model

This paper focuses on the linear parameter varying (LPV) modeling and controller design for a flexible air-breathing hypersonic vehicle (AHV). Firstly, by selecting the measurable altitude and velocity as gain-scheduled variables, the original longitudinal nonlinear model for AHV is transformed into the LPV model via average gridding division, vertex trimming, Jacobian linearization, and multiple linear regression within the entire flight envelope. Secondly, using the tensor product model transformation method, the obtained LPV model is converted into the polytopic LPV model via high-order singular value decomposition (HOSVD). Third, the validity and applicability of the HOSVD-based LPV model are further demonstrated by designing a robust controller for command tracking control during maneuvering flight over a large envelope.


Introduction
Since the 50s and 60s of the 20th century, with the series of events such as the advent of long-range ballistic missiles, the successful return of manned spacecraft, and the X-15 test aircraft velocity faster than Mach 6, it marked that mankind formally entered the era of hypersonic velocity. The hypersonic vehicle has taken advantages of both spacecraft and aircraft and has become the new technological commanding point in the 21st century aerospace field. It has been hailed as the third epoch-making milestone in the history of aviation in the world after the invention of aircraft and the breakthrough of the sound barrier [1][2][3].
Compared with a conventional aircraft, the complex aerodynamic characteristics of hypersonic vehicles and the application of technologies such as scramjet engines and integrated airframe/engine configuration make the coupling more powerful between elastic bodies, propulsion systems, aerodynamics, and structural dynamics of hypersonic vehicles; besides, the physical model is more nonlinear. At pres-ent, the main means of dynamic modeling for hypersonic vehicles are computational fluid dynamics (CFD), wind tunnel tests, aircraft tests, and other test methods. Because of the imperfection of hypersonic aerodynamic database and CFD software, the low degree of wind tunnel test approximation, the high cost and risk of flight tests, and the limitation of hypersonic vehicles modeling methods, the system modeling has great structural and parameter uncertainties. Second, the large-span flight envelope of hypersonic vehicles not only requires the flight control system to adapt to different flight environments but also requires to satisfy the strict constraints such as more complex dynamic pressure, overload, flight envelopes, and actual physical characteristics (e.g., actuator saturation). In summary, a hypersonic vehicle is a highly dynamic, strongly coupled, fast time-varying, and multivariable nonlinear system. Its control system design not only faces the influence of model uncertainty and external environment but also faces the constraints of actual existence and the time-delay dynamic limitation, which makes the control, stability, safety, rapidity, and accuracy of hypersonic vehicle systems face great challenges. It has become one of the key issues to be solved and is also a hot issue in the field of flight control, such as adaptive control strategies based on the multiple Lyapunov function method or using an auxiliary reaction control system for airbreathing hypersonic vehicles [4][5][6][7][8].
Recently, the LPV controller has emerged as a reliable alternative to classical gain scheduling approaches and provides guaranteed stability, robustness, or performance [9,10]. Gain-scheduled control strategies for LPV systems have been developed intensively and bring promising solutions in many fields, such as missiles, aircrafts, spacecrafts, winding systems, automated vehicles, and robotic systems [11][12][13]14]. Zhang handled the reference tracking for a class of flexible hypersonic vehicles using a novel switched LPV framework [15]. Yang proposed a fault-tolerant controller based on robust model-predictive control and A polytopic LPV model for A hypersonic vehicle with external disturbances and actuator loss of effectiveness faults [16]. Hu proposes a novel passive fault-tolerant control method using weighted tube-based model predictive control via polytopic LPV for air-breathing hypersonic vehicles [17]. However, the formulation of the polytopic LPV model for nonlinear systems and controller design for LPV systems with uncertainty and external disturbance remain difficult, which motivates our study.
This paper takes the flight control system design of the hypersonic vehicle as the background. In the process of modeling, it is necessary to deal with various factors that affect the performance of the control system. Aircraft modeling errors, flight parameter changes of large package flight, elastic vibration deformation of aircraft, and other factors are regarded as internal factors leading to system uncertainty. Factors such as the complex flight environment, dynamic pressure effect, aerodynamic thermal effect, and the frictional resistance effect are regarded as external disturbances. Then, the attitude stability and command tracking issues during the hypersonic vehicle maneuvering flight in large envelope are transformed into the robust controller synthesis and guaranteed performance control of the uncertain disturbance system. Firstly, for a nonlinear longitudinal elasticity model of air-breathing hypersonic vehicles, the altitude and velocity are selected as scheduling parameters, uniform meshing is performed within the range of parameter variation, and then, the LPV model of the aircraft is obtained by Jacobi linearization and multivariate linear fitting. Secondly, the LPV model of the aircraft is transformed into a polytopic LPV model by the tensor product model change method, which lays a foundation for the subsequent aircraft controller design. Finally, a robust controller is designed based on the obtained aircraft polytopic LPV model, which verifies the effectiveness of the aircraft polytopic LPV model in the command tracking control of aircraft maneuvering flight in large envelopes.  [18][19][20], which takes into account the influence of the elastic modality of the aircraft. The simplified longitudinal elastic nonlinear model of the aircraft is as follows:

Problem Formulation
In equation (1), h represents the aircraft altitude, V represents the aircraft velocity, α indicates the angle of attack, θ denotes the pitch angle, Q represents the pitch rate, g indicates the gravitational constant, M is the pitching moment, I yy represents the moment of inertia around the body axis, T is the engine thrust, D represents resistance, L represents lift, η i (i = 1, 2) represents the elastic mode, ω i (i = 1, 2) represents the natural frequency of the elastic mode, ς i (i = 1, 2) represents the damping coefficient of the elastic mode, and N i denotes the generalized force.
The approximate expressions of the above aerodynamics, moments, thrusts, and generalized forces fitted near the equilibrium point are shown in equation (2). In equation (2), ρindicates air density, S is the reference area, C L ðα, δ e Þ represents the lift coefficient related to the angle of attack and elevator deflection, C D ðα, δ e Þ indicates the drag coefficient related to the angle of attack and elevator deflection, C M,α ðαÞ indicates the torque coefficient related to the angle of attack, C M,δ e ðδ e Þ indicates the torque coefficient related to the elevator deflection, z T represents the coupling coefficient between thrust and pitching moment, C α j T represents the thrust coefficient related to α j (i = 1, 2, 3), N α j i represents the generalized force coefficients related to α j , and N δ e 2 denotes the influence coefficient of elevator declination to generalized force N 2 , and C 0 T and N 0 i represent constant items in thrust and generalized forces, respectively.
International Journal of Aerospace Engineering The relevant force and moment coefficients are shown in equation (3), where 1/h s is the air density decay rate, C 0 L and C 0 D are the constant term in C L ðα, δ e Þand C D ðα, δ e Þ, respectively, C is the constant term in C M,α ðαÞ, q is the dynamic pressure, and Ф is the stoichiometrically normalized fuel-to-air ratio, β i ðh, qÞ is the ith thrust fit parameter. The parameter reference values of the elastic nonlinear model are detailed in Hu et al. [21] and listed as Table 1.
Select xðtÞ = ½h, V, α, θ, Q, η 1 , _ η 1 , η 2 , _ η 2 T as state variables, select uðtÞ = ½Φ, δ e T as system inputs, and select zðtÞ = ½V, h T as system outputs. Substitute equations (2) and (3) into (1) to obtain the affine nonlinear model of the system as follows: In equation (4), f i ðx, tÞ (i = 1, 2, ⋯, 9) is shown on the right side of equation (1). The equilibrium state of a hypersonic vehicle is a specific flight condition under which the line velocity, and angular velocity are constant or zero and the acceleration is zero under the given constraints. Therefore, the system equilibrium point sequence can be solved by equation (5) as follows: Within the permissible range of parameters, the velocity and altitude space of the selected scheduling variables are gridded and each point on the grid is calculated and balanced according to equation (5) to obtain the equilibrium working point of the aircraft. A set of LTI systems is obtained by using the Jacobian linearization method. Finally, these linearized local systems are multivariate linearly fitted to obtain the LPV model. The equidistant gridding method is chosen in this paper. Five points are selected at equal intervals in the velocity interval ½7200 9600 ft/s, and four points are selected at equal −44:201 × 10 2 lb · ft −1 · slug −0:5 3 International Journal of Aerospace Engineering intervals in the altitude interval ½80000 110000 ft, totaling 20 points, as shown in Figure 1.
By means of parameter meshing, Jacobian linearization, and multivariate linear fitting, the LPV model of the hypersonic vehicle with respect to the scheduling parameters is shown in equation (6). The fitting results of partial LPV model parameters are shown in Figure 2, and its corresponding fitting results are listed as follows.  International Journal of Aerospace Engineering the scheduling parameters means that an infinite number of LMIs are to be solved over the entire parameter trajectory, which is obviously difficult to achieve, the scheduling parameters can be spatially meshed to solve the problem of LMI at a finite number of mesh vertices but the performance of this method is more dependent on the mesh density; besides, the determination of the mesh density not only lacks theoretical proof but also is time-consuming. The polytopic set is a convex set, which could describe the entire system. Through transformation of the polytopic model, the problem of solving infinitely many LMI problems can be attributed to solving the LMI problem at a limited number of vertices, which can greatly reduce the amount of computation. The controller obtained by this method has continuous gain throughout the workspace and has global action features ([11] [12,22]). The system equation (6) can be expressed as follows.
where SðρðtÞÞ is a parameter-varying object and ρðtÞ ∈ Ω is a time-varying N dimensional parameter vector within the closed hypercube Parameter ρðtÞ can also include some elements of xðtÞ.
The ðN + 2Þ-dimensional tensor S is constructed from LTI vertex systems S i 1 i 2 ,⋯,i N .
where the row vector A n ðρ n ðtÞÞ ∈ ℝ I n (n = 1, ⋯, N) contains one bounded variable and continuous weighting functions α n,i n ðρ n Þ (i n = 1 ⋯ I n ). The weighting function α n,i n ðρ n ðtÞÞ ∈ ½0, 1 is the i n weighting function defined on the nth dimension of Ω, and ρ n ðtÞ is the nth element of vector ρðtÞ and ∑ I n i n =1 α n,i n ðρ n ðtÞÞ = 1. I n (n = 1, ⋯, N) denotes the number of the weighting functions used in the nth dimension of Ω. The symbol ⊗ denotes the Kronecker product. The tensor product model transformation method based on high-order singular value decomposition is an effective LPV system polytopic model transformation method. The brief process is as follows. Firstly, mesh the variable parameters in the LPV model, and then, discretize the system to form tensors. Finally, we perform high-order singular value decomposition. The related theory of high-order singular value decomposition can refer to the literature [23][24][25].
On the basis of the LPV model of the hypersonic vehicle in equation (6), the tensor product model transformation method is used to obtain the polyhedral description of the LPV model of the hypersonic vehicle. By uniformly dividing the altitude variable h and velocity variable V in the LPV model described in equation (6), the original system is approached in the specified area (velocity interval ½7200 9600 ft/s, altitude interval ½80000 110000 ft), where the number of sampling points is set to 200 × 200. Using the high-order singular value method decomposition, the singular value decomposition of the expansion matrix is arranged in a descending order, where N = 2, I 1 = I 2 = 3. Table 2 shows the top 5 largest singular values. All 6 (the remaining singular values are relatively small or equal to 0) of the singular values are preserved. The resulting numerical model is a convincing approximation of the original model, which has 9 polytopic vertex systems. The resulting approximation model is as follows:   International Journal of Aerospace Engineering In equation (10), w k ðV, hÞ = w i ðVÞw j ðhÞ, k ∈ f1, 2,⋯,9g, i, j ∈ f1, 2, 3g. Based on the high-order singular value decomposition results, the maximum modeling error of the approximation model described in equation (8) can be calculated as 0.003679 and the average modeling error is 0.002499.
In order to verify the approximate equivalence between the hypersonic elastic nonlinear model described in equation (1) and the polytopic LPV model obtained by the transformation of the tensor product model described in equation (8), an equilibrium state in the parameter variation interval in Table 3 is used as an example. The system approximation matrix at the equilibrium point is obtained by the Jacobian linearization method which can locally approximate the elastic nonlinear model and the tensor product model transformation method in equation (8). The zero-pole distribution map of the system is shown in Figure 3.
Comparing the zero-pole positions of the two methods at the same equilibrium point, we can find that the zeropole positions of the two methods are very close, which indicates that the polytopic LPV model based on the transformation of the tensor product model accurately reflects the dynamic characteristics of longitudinal elasticity nonlinear models of the hypersonic vehicle.
The general distribution of the nine time-varying weight functions shows that when any time-varying weight is 1, the rest of the weights are 0, which indicates that the LPV model weight function obtained by the tensor product model transformation satisfies the requirement of convex coordinates.

Robust Controller Design
The LPV model of a hypersonic vehicle is as follows.
In equation (11), θðtÞ = ½V h T denotes the scheduling variables, AðθðtÞÞ and BðθðtÞÞ both are the corresponding parameter-dependent system matrix, and the numerical expression is described in equation (6). C is the observation matrix, D is the interference matrix, and wðtÞ denotes time-varying exogenous disturbance. The reference velocity and altitude of the aircraft are recorded as rðtÞ = ½V com h com T and the control output is recorded as yðtÞ = ½V t h t T . The controller design goal of the hypersonic vehicle not only needs to ensure the stability of the aircraft but also requires that the output of the aircraft can track a given command signal in real time. Therefore, the design goal of the hybrid target robust controller is to design an effective controller within the whole flight envelope ensuring that the output tracking error of the aircraft is 0.
In order to evaluate the system's steady-state error more effectively, the following tracking error integral term is introduced to eliminate the steady-state error in the feedback control strategy.
According to the tracking error integral term and combining polytopic description and the LPV model of the hypersonic vehicle, the extended system can be obtained as follows.
where  According to the parallel distributed compensation (PDC), the designed state feedback controller is as follows.
The output tracking problem of the hypersonic vehicle is transformed into the robust stability problem of the above extended system. When the closed-loop system described in equation (17) is robustly stable, the command tracking target of the original system described in equation (9) can be achieved.   International Journal of Aerospace Engineering

Numerical Simulation
Under conditions of cruising flight of the hypersonic vehicle (altitude 85000 ft, velocity 7702 ft/s), a robust model predictive controller for state feedback of mixed performance indicators was designed to achieve tracking control of velocity commands and altitude commands [26]. In the actual control of the aircraft, the command must comprehensively consider the physical constraints of the aircraft such as structural stiffness and payload. In general, it is assumed that the given velocity and altitude reference command signals are generated by the second-order filters.
In equation (18), the velocity instruction signal filtering system natural frequency w V = 0:0195 and damping coefficient ζ V = 0:9. The altitude command signal filtering system natural frequency w h = 0:0219 and damping coefficient ζ h = 0:9. It is assumed that the initial equilibrium state of the longitudinal dynamics of the aircraft is as follows. The discretization sampling time T s of the continuous system is 0.01 s. The aerodynamic parameters in the aircraft model are mainly obtained by computational fluid dynamics simulations or wind tunnel experiments. Besides, the unknown external environment and complex dynamic influences such as dynamic pressure effects, viscous effects, lean gas effects, and atmospheric turbulence may be encountered in the flight process. Therefore, the abovementioned factors are assumed as uncertainties of model inertial parameters and aerodynamic parameters. Similar to the literature [27], the value of parameters lies within ±10% of the nominal value such as m = mð1 + △mÞ where j△mj ≤ 0:07 and ρ = ρð1 + △ρÞ where j△ρj ≤ 0:08. Compared with the ones in Qin et al. [27], the uncertainties are more aggressive and energy-bounded external disturbance is considered. When the desired velocity V e = 9200ft/s and desired altitude h e = 104000 ft, the simulation results using the state feedback controller design method are shown in Figures 5-9. The simulation results are compared with the results obtained through Qin's method [27] and Cai's method [28]. In the sequence, results obtained through the proposed method, Qin's method, and Cai's method are depicted by a red line marked by circles, blue dashed line marked by circles, and green dashed line marked by asterisks. Since the tracking curves are quite close, thus, only the tracking errors are given.     International Journal of Aerospace Engineering In Figures 5 and 6, it can be seen that the robust controller can ensure that a hypersonic vehicle achieves the fast and accurate velocity and altitude command tracking during the maneuvering flight of the large flight envelope. The above controller has a smaller tracking error and a relatively smooth state, which can effectively suppress the influence of serious aerodynamic parameter uncertainties and elastic modes in the flight process. Figures 7 and 8 show the variation curves of the angle of attack and pitch angle in the whole process, respectively. In order to achieve a high degree of climb while increasing the velocity, the hypersonic vehicle in the initial trim state adds engine fuel input to the controller output signal and the elevator is pulled up until the stable command tracking is achieved again to reach the equilibrium state. The angle of attack and the pitch angle will be stabilized in the new trim state. The curve of the control input shown in Figure 9 indicates that the engine fuel equivalence ratio and the elevator deflection angle all satisfy the system control input hard constraints.
In summary, in the case of external disturbances, parameter uncertainties, and elastic modal couplings during the maneuvering flight, the designed controller has strong robustness and ensures that the system can quickly track the command signal.

Conclusion
Aiming at a type of hypersonic vehicle longitudinal elasticity model disclosed in the existing literature, altitude and velocity are selected as the gain-scheduled variables, through a series of sequential steps such as meshing within the range of parameter variation, finding equilibrium points, Jacobian linearization, and multivariate linear fitting; the continuoustime LPV model of the aircraft is established. Then, using the tensor product model transformation method, the general LPV model is transformed into a polytopic LPV system model with 9 vertices and the rationality of the LPV model is verified. At the same time, regarding the maneuvering flight command tracking problem for hypersonic vehicles, a state feedback controller is used to translate the command tracking control problem into a robust stability problem for the corresponding extended closed-loop system combined with the obtained polytopic LPV system model. Simulation results demonstrate that the obtained polytopic LPV model based on high-order singular value decomposition in this paper is effective in solving the maneuvering flight command tracking and attitude stabilization control of general complex nonlinear systems, such as hypersonic vehicles. However, since the linearization process lose some beneficial nonlinearities, some gain perturbations due to finite word length or time delays of state or input might bring unpredictable complexities. LPV nonfragile control of air-breathing hypersonic vehicles with time delays needs further studies.

Data Availability
The longitudinal elastic nonlinear model of the hypersonic vehicle used in this paper is mainly based on the wingedcone model provided by NASA's Langley Research Center [18,19,20], which takes into account the influence of the elastic modality of the aircraft.

Ethical Approval
This article does not contain any studies with human participants or animals performed by any of the authors.

11
International Journal of Aerospace Engineering