Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. Robot. AI, 13 June 2017 | https://doi.org/10.3389/frobt.2017.00023

Characteristic Model-Based Robust Model Predictive Control for Hypersonic Vehicles with Constraints

  • 1Department of Electrical and Information Engineering, Jiangsu University, Zhenjiang, China
  • 2Department of Biomedical Engineering, National University of Singapore, Singapore, Singapore

Designing robust control for hypersonic vehicles in reentry is difficult, due to the features of the vehicles including strong coupling, non-linearity, and multiple constraints. This paper proposed a characteristic model-based robust model predictive control (MPC) for hypersonic vehicles with reentry constraints. First, the hypersonic vehicle is modeled by a characteristic model composed of a linear time-varying system and a lumped disturbance. Then, the identification data are regenerated by the accumulative sum idea in the gray theory, which weakens effects of the random noises and strengthens regularity of the identification data. Based on the regenerated data, the time-varying parameters and the disturbance are online estimated according to the gray identification. At last, the mixed H2/H robust predictive control law is proposed based on linear matrix inequalities (LMIs) and receding horizon optimization techniques. Using active tackling system constraints of MPC, the input and state constraints are satisfied in the closed-loop control system. The validity of the proposed control is verified theoretically according to Lyapunov theory and illustrated by simulation results.

Introduction

Hypersonic vehicle, a type of aircrafts with a flying speed of over 5 Mach, is a hot research topic due to its importance in national defense and military affairs. The hypersonic reentry flight usually covers a large flight envelope, during which the environmental and aerodynamic characteristics undergo huge variations (Harpold, 1979; Lu, 2014). The reentry dynamic model has fast time variability and strong coupling, is highly non-linear and is with strong disturbances and uncertainties (Wu and Meng, 2009). Thus, it is difficult to adapt the large envelope and multitasking flight movement using the traditional PID control, as the selection of the operating points and the design of the scheduling algorithm remains a time-consuming procedure for engineers. Therefore, domestic and foreign researchers study the advanced controller design to make the vehicles have stable flight characteristics and strong robustness (Zhang and Hu, 2012). Since there exist thermal and mechanical constraints in safety boundaries of reentry flight corridors, the hypersonic reentry control system must satisfy strict input and state constraints (Lu, 1999; Shen and Lu, 2003), including the constraints on aerodynamic angles, angular velocities, and control moments. If these complex physical demanding constraints are violated, the performance and stability of the control system will be deteriorated, then the safety and reliability of reentry flight can be affected seriously.

In the past decades, non-linear control methods such as adaptive backstepping control, non-linear disturbance observer, and robust control were used to improve robustness of the reentry control systems (Xu et al., 2012; Li et al., 2016; Wu and Meng, 2016). But, the system input and output constraints were not considered in Xu et al. (2012), Li et al. (2016), and Wu and Meng (2016). Well-known methods for constrained control include anti-windup control and reference governors. In Alicia and Andrea (2009), the combined reference governor and anti-windup control for hypersonic vehicles can suppress input saturation effectively. However, the compensator will be executed only when the saturation phenomenon appears, which affects the control accuracy and strong maneuverability of reentry hypersonic vehicles.

Model predictive control (MPC) has been widely accepted as an effective approach for control of constrained systems and applied in many industries, for its ability to handle hard constraints actively in receding horizon optimization (Chiristopher et al., 2013; Zhang and Sun, 2016). But high online-computation complexity brings difficulties to its applications in real-time control of hypersonic vehicles with high non-linear systems. Thus, much progress has been made to reduce computation complexity of MPC for hypersonic vehicles. In Van (2006), feedback linearization and linearized MPC were combined to design the control law for hypersonic vehicles. By transformation of the non-linear model into a linear model, it simplified the complexity of the predictive model. But, the computation was still complicated, due to the convention of original linear constraints into non-linear constraints. In Qin et al. (2016), a fuzzy predictive control strategy was designed for hypersonic vehicles with actuator constraints and input delays. Though a linearized predictive model is adopted to decrease model complexity, computational complex of the MPC is still high due to complex fuzzy rules.

Based on the above analysis, more research deserves to be done on MPC for hypersonic vehicles with permitted complexity and high robustness. Characteristic modeling, which was proposed by Wu et al. (2007), is constructed based on dynamical characteristics, environmental characteristics, and control performance, no only based on accurate system dynamics in conventional modeling (Wu et al., 2009). The main advantages of the characteristic modeling are depicting high-order linear time-varying systems and non-linear systems by low-order linear time-varying systems with a simple structure. The features including non-linear, strong coupling, and uncertainty for hypersonic flight make it hard to establish an accurate dynamic model, but the characteristic model can be well applied to this kind of complex system. Thus, such modeling method highly decreases difficulties in control design. A golden section adaptive attitude controller is presented in Gong and Wu (2010) and Xu (2012) based on the characteristic model. But the coupling in states and the coupling in controllers were not considered. The control (Gong and Wu, 2010; Xu, 2012) is difficult to be applied in the hypersonic vehicles in reentry stage since the hard constraints for the reentry vehicle were not considered.

The aim of the paper is to design characteristic model-based robust MPC for hypersonic vehicles with reentry constraints. According to input–output equivalence, dynamics of the hypersonic vehicle is described by a characteristic model composed of a linear time-varying system and a lumped disturbance. The time-varying parameters and the disturbance in the characteristic model are estimated according to the gray identification, which highly reduces effects of the measurement noises. The proposed robust MPC is obtained by solving an optimization with the characteristic model as the predictive model and a mixed H2/H function as the cost function. In the control, the constraints are tackled by explicitly considering them in the optimization and the control robustness is improved by the use of the H2/Hfunction. The feasibility and robustness of the proposed control are verified theoretically according to Lyapunov theory and illustrated by simulation results. The main contributions of this paper include the following:

(1) design of the characteristic model with time-varying parameters as the predictive model, which highly reduces complexity of the predictive model and approximation errors in Taylor linearization.

(2) identification of the time-varying parameters according to the gray identification, which reduces the effects of stochastic noises to identification accuracy.

(3) design of the H2/H MPC, which satisfies reentry constraints and improves control robustness to uncertainties.

The design process of this paper is organized as follows. In Section “Model Description and Problem Formulation,” the objective of this paper is stated. In Section “Characteristic Model of the Reentry Hypersonic Vehicle,” a characteristic model composed of a linear time-varying system and a lumped disturbance are constructed for the hypersonic vehicle. In Section “Gray Identification for the Characteristic Model,” the time-varying parameters and the disturbance for characteristic model are online estimated according to the gray identification. In Section “Mixed H2/H Robust Predictive Control,” the mixed H2/H robust predictive control law is proposed based on linear matrix inequalities (LMIs) and receding horizon optimization techniques. Section “Simulation Analysis” presents simulation result. At last, conclusions are given in Section “Conclusion.”

Model Description and Problem Formulation

To study the attitude stability of the winged hypersonic vehicle, it is common to neglect the wind, the Earth rotation and the Earth curvature. Using linear aerodynamics hypothesis and neglecting the influences of control surface deflections on aerodynamic forces, the rigid body attitude dynamic equations are as follows (Zhang and Hu, 2012):

α˙=ωyωxcosαtanβωysinαtanββ˙=ωxsinαωzcosαφ˙=ωx+ωysinφtanθ+ωzcosφtanθ
ω˙x=I1ωyωz+I2ωxωy+kxxωx+kxzωz+kxαδa+krδr+M1ω˙y=I3ωxωz+I4(ωx2ωz2)+kyyωy+kyeδe+M2ω˙x=I5ωxωy+I6ωyωz+kzxωx+kzzωz+kzαδa+kzrδr+M3

and

I1=IyIzIz2Ixz2IxIzIxz2,I2=(IyIxIz)IxzIxIzIxz2,I3=IzIxIy,I4=IxzIy,I5=Ix2+Ixz2IxIyIxIzIxz2,I6=(Iz+IxIy)IxzIxIzIxz2

where the attack angle α is the angle between the velocity vector and its longitudinal symmetry plane; the side slip angle β is the angle between the projection of the velocity vector and the vehicle vertical axis; φ and θ are the roll angle and, respectively, the pitching angle; ωx, ωy, ωz are angular velocities of the flight vehicle projected to the body frame; δr, δe, δa are rudders deflection, elevators deflection, and left right elevator difference, respectively; kxx, kxz, kxα, kr, kyy, kye, kzx, kzz, kzα, kzr are the aerodynamic moment coefficient, the detailed definition can be seen in Liu and Chen (2010); M1, M2, M3 are external moment disturbances. The attack angle and the side slip angle are depicted in Figure 1. The matrix of inertia moments is given by I=Ix0Ixz0Iy0Ixz0Iz. The task of control system is as follows:

(1) Control goal: with the constraints satisfied, track the outer guidance commands α, φ and make β keep at 0 in the presence of parameter uncertainties and disturbance torques. The control errors of α and φ are within 0.5.

(2) The constraints of the reentry control include

① The saturation constraints of rudder deflection

δa2δamax,δe2δemax,δr2δrmax

② The state constraints

FIGURE 1
www.frontiersin.org

Figure 1. Relationship between the speed coordinates (Oxayaza) and the vehicle coordinates (Oxbybzb).

According to the requirements of preventing heat, the following constraints should be satisfied:

α(k)2αmax

According to the structural safety requirements of vehicle mechanics, the following constraints of angular velocities need to be satisfied:

ωx2ωxmax,ωy2ωymax,ωz2ωzmax.

Characteristic Model of the Reentry Hypersonic Vehicle

The reentry systems (1) and (2) are strongly non-linear, time-varying, and uncertain multivariable coupling, so it is difficult to design controller directly for it. Based on the principle of input and output equivalence, we differentiate system outputs α, β, φ to second order, apply identical substitution and Taylor series approximate expansion (Wu et al., 2007). The discrete model of input–output time-varying second-order of α, β, and φ is as follows:

α(k+1)=fα1(k)α(k)+fα2(k)α(k1)+g1(k)δα(k)+Δα(k)
β(k+1)=fβ1(k)β(k)+fβ2(k)β(k1)+g2(k)δβ(k)+Δβ(k)
φ(k+1)=fφ1(k)φ(k)+fφ2(k)α(k1)+g3(k)δφ(k)+Δφ

where characteristic model parameters fi1, fi2 (i = α, β, φ) and g1, g2, g3 are time-varying, δα, δβ, δφ are equivalent control input. The characteristic model for the plant (Eqs 3–5) as given in Eqs 1 and 2 possesses the following properties: the coefficients fi1, fi2 (i = α, β, φ) and g1, g2, g3 are slowly time-varying, the ranges of these parameters can be decided a priori, in response to the same input, the output of the original model (Eqs 1 and 2) is identical to the output of the characteristic model (Eqs 3–5) at each sampling time. Sum of all the coefficients of (Eqs 3–5) equals 1 (Zhang and Hu, 2012) as follows:

fα1+fα21δαδefβ1+fβ21δβδrfφ1+fφ21δφδa

Define the lumped disturbance Δα(k) as follows:

Δα(k)=m11β(k)+m12φ(k)+Δα1(k)

where m11, m12, and Δα1(k) are coupling coefficient, Δα(k) includes dynamic coupling information with other channels, external disturbances, and unmodeled errors. To improve robustness of the control system, Δα(k) is estimated online and compensated feedforwardly in control. Similarly, the lumped disturbances in the other channels are described as follows:

Δβ=m21α(k)+m22φ(k)+Δβ1(k)Δφ=m31α(k)+m32β(k)+Δφ1(k)

Let δα ≈ δe, then the control constraints can be transformed into δα(k)2δemax [the readers can see Zhang and Hu (2012) for detail], the output constraints are α(k)2αmax and α˙2ωy2ωymax.

Gray Identification for the Characteristic Model

To accurately identify the parameters in the characteristic model is significant for control design.

In this section, the random noise in the measured reentry dynamic state can be approximately treated as white noise. Accumulated generating operation in the gray identification is used to weaken the impact of the random noise on original data sequences efficiently (Liu and Chen, 2010) and regularize original data.

Let the original reentry state sequence be

x(0)={x(0)(1),x(0)(2),........,x(0)(N)}

where N is the total number of the identification samples. The identification data can be regenerated as follows:

x(1)(k)=m=1kx(0)(m)

where m has the number of sampling data. Operating gray accumulation generation once, we can get more regular sequence, which eliminates the impact of random noise to some extent. Then, new input data for identification are depicted as follows:

{x1(1),x1(2),........,x1(N)}

In the following of this section, taking the attack-angle model (Eq. 3) as an example to illustrate identification of parameters. The identification of parameters in other models (Eqs 4 and 5) is similar and omitted. From the definition of gray data, the gray model of dynamics of the attack angle (Eq. 3) is given as follows:

α1(k+1)=fα1(k)α1(k)+fα2(k)α1(k1)+g1(k)δα1(k)+m11β1(k)+m12φ1(k)+Δα11(k)

During the identification time, the unmodeled error and disturbances are approximated as follows:

Δα11(1)=Δα1,Δα11(2)=2Δα11(1)=2Δα1,Δα11(N1)=(N1)Δα11(1)=(N1)Δα1.

Let

V=[fα1(k),fα2(k),g1(k),m11,m12,Δα1]T
B=B2B3...BN1=α1(2)α1(1)δα1(2)β1(2)ϕ1(2)2α1(3)α1(2)δα1(3)β1(3)ϕ1(3)3......α1(N1)α1(N2)δα1(N1)β1(N1)ϕ1(N1)N1

From the gray model, we have

DN(1)=BV

where DN(1)=[α1(3)α1(4)α1(N)]. Applying the least square method, we can get

V=(BTB)1BTDN(1)

To decrease effects of the lumped disturbances and improve robustness of the closed-loop control system, the feedforward input compensation is designed as follows:

uc(k)=g11(k)(m11β(k)+m12φ(k)+Δα1(k))

To accurately identify the parameters in the characteristic model is significant for control design.

Mixed H2/H Robust Predictive Control

Let Δf α1(k), Δf α2(k), Δg1(k) be the identification deviations of the actual characteristic model coefficients. ΔT is the compensation error after the compensation of reentry channel coupling information, unmodeled dynamics and disturbances in uc(k). It is obvious that a Δf α1(k), Δf α2(k), Δg1(k), and ΔT are bounded. Besides. Since α(k), δα(k) have definite region constraints, they are also bounded. As a result, the generalized disturbances d(k) are bounded and can be described as follows:

d(k)=Δfα1(k)α(k)+Δfα2(k)α(k1)+Δg1(k)δα(k)+ΔT

Then characteristic model after gray identification compensation is given as follows:

α(k+1)=fα1α(k)+fα2α(k1)+g1δα1(k)+d(k)
δα=uc+δα1

Since k0dT(k)d(k)dmax is satisfied, d(k) is energy bounded. Since the unmodeled errors and couplings are partly compensated by estimators proposed in the last section, robust control is still needed to improve robustness of the closed-loop control system.

In this section, based on mixed H2/H approach, a robust predictive control is proposed to have the better control performance under the premise that the general disturbance d(k) has been restrained efficiently (Patience and Orukpe, 2007).

Assume the tracking expectation is αc = 0. To reduce the tracking errors, the integral of attack angle e=[(ααc)dt] is considered as an extended state. Then, the augmented system is constructed as follows:

x(k+1)=A(k)x(k)+B(k)δα1(k)+Bdd(k)
A(k)=0100010fα2fα1,B(k)=00g1,Bd=001,x(k)=e(k)α(k)α(k+1)

In the subsequent control design, we use the characteristic model (Eq. 14) as a prediction model to simplify the complexity of predictive control. Considering the reentry thermal protection and overload safety, the constraints of the attack angle, the velocity and the input for the hypersonic vehicle control can be transformed into the following forms:

(1) constraints of the attack angle

α(k)2=[0  1  0]x(k)2αmax

(2) constraints of the velocity

ωy2α˙201 1x(k)1T2ωymax

(3) the input saturation

u=δα1+uc2δemaxδα12min{δemaxuc,δemaxuc}=umax

Define the indicator of control performance: z(k)=Q1x(k)Rδα1(k), where Q1 and R are weighting matrices of the state and control. Consider the following indexes (Jiang et al., 2017):

H performance index: for any arbitrary γ > 0, the transfer function from d(k) to z(k) satisfies

Tzdγ

② H2 performance index

z2=i=0zT(k)z(k)a

This paper is to design controller based on robust MPC and LMI, which is efficient to avoid adjusting parameters by experience, handle various reentry constraints and guarantee the robustness and stability of the time-varying characteristic system.

Theorem 1: Let x(k) be the state of the system (Eq. 14) measured at the sampling time k. Assume that there are constraints on the control input and plant output (Eqs 15–17). Then, the state feedback K in the control law δα1(k) = Kx(k) that minimizes the upper bound H2 performance at sampling time k is given by the following equation:

K=YQ1

where Q > 0 and Y are obtained from the solution (if it exists) of following linear objective minimization problem:

minY,Q,aa

subject to

1γ2dmaxaγ2dmaxIx(k)0Q0

and

Q00aγ2IA(k)Q+B(k)YaBdQQ1Q00aIRY000aI>0
αmax2I(A(k)Q+B(k)Y)TC1TC1(A(k)Q+B(k)Y)Q0C1=[0 1 0]
ωymax2I(A(k)Q+B(k)Y)TC2TC2(A(k)Q+B(k)Y)Q0C2=1T[01 1]
umax2IYYTQ0

where * represents terms readily inferred from symmetry.

Proof:

Stability constraints

Define V(x(k)) = xT(k)P(k)x(k) and V(x(k + 1)) = x (k + 1)TP(k)x(k + 1), where x(k + i) indicates the predictive state at k + i based on the measurement x(k) at sampling time k. δα1(k + i) is the ith predictive control at the sampling time k, we get the following formula:

V(x(k+i+1))V(x(k+i))=xT(k+i)(A(k)+B(k)K)TP(k+1)(A(k)+B(k)K)x(k+i)+dT(k+i)BdTP(k)Bdd(k+i)+2dT(k+i)BdTP(k)(A(k)+B(k)K)x(k+i)xT(k+i)P(k)x(k+i)+zT(k+i)z(k+i)zT(k+i)z(k+i)+γ2dT(k+i)d(k+i)γ2dT(k+i)d(k+i)=[xT(k+i)dT(k+i)]M(k)[xT(k+i)dT(k+i)]T+γ2dT(k+i)d(k+i)zT(k+i)z(k+i)

where

M(k)=[(A(k)+B(k)K)TP(k)(A(k)+B(k)K)P(k)+Q1TQ1+(RK)T(RK)(BdP(k)(A(k)+B(k)K))TBdP(k)(A(k)+B(k)KBdTP(k)Bdγ2I ]

From stability of gray identification, we can know that compensation error is close to 0. Then, if z(k + i) is also close to 0, then the system stability is obtained. If not, the following inequality is satisfied:

γ2dT(k+i)d(k+i)zT(k+i)z(k+i)<0

If M(k) < 0, then V(x(k + i + 1)) − V(x(k + i)) < 0, and x(∞) = 0, V(x(∞)) = 0. Summing the both sides of Eq. (20) from i = 0 to i = ∞, we can get the following:

xT(k)Px(k)=i=0[γ2dT(k+i)d(k+i)zT(k+i)z(k+i)]+i=0[xT(k+i)dT(k+i)]M(k)[xT(k+i)dT(k+i)]Ti=0[zT(k+i)z(k+i)]=xT(k)P0x(k) +i=0[γ2dT(k+i)d(k+i)]+i=0[xT(k+i)dT(k+i)]M(k)[xT(k+i)dT(k+i)]T

(I) When x(0) = 0, if M(k) < 0, i=0[γ2dT(k+i)d(k+i)zT(k+i)z(k+i)]0, the H performance can be satisfied.

(II) When M(k) < 0, i=0[zT(k+i)z(k+i)]xT(k)Px(k)+γ2dmaxa, the H2 performance can be satisfied.

Therefore, the necessary and sufficient condition to realize the control target is M(k) < 0, i.e.

M(k)=P(k)+Q1TQ1+(RK)T(RK)00γ2I+(A(k)+B(k)K)TBdTP(k)(A(k)+B(k)K)Bd<0

Using Schur complement, the LMI (Eq. 22) is equivalent to

P(k)00γ2IA(k)+B(k)KBdP(k)1Q100IRK000I>0

Let Q = aP−1, K = YQ−1. Pre-multiply and post-multiply Eq. 23 by diag{a1/2P−1, a1/2I, a1/2I, a1/2I, a1/2I}, we can get

Q00aγ2IA(k)Q+B(k)YaBdQQ1Q00aIRY000aI>0

Performance constraints

Using Schur complement property, H2 performance (19) is equivalent to:

1γ2dmaxaγ2dmaxIx(k)0Q0

Output constraints

Since M(k) < 0, V(x(k + i + 1)) ≤ V(x(k + i)), and {z: zT(k)Q−1z(k) ≤ 1} is an invariant ellipsoid for the predicted states of the uncertain system (Eq. 13). Then

maxi1y(k+i|k)2=maxi1C(A(k+i) +B(k+i)K)x(k+i|k)2maxzTz1C(A(k+i)+B(k+i)K)Q12z2=σ[C(A(k+i)+B(k+i)K)Q12],i0

Since y(k+i|k)2ymax, then

σ[C(A(k+i|k)+B(k+i|k)K)Q12]ymax

which is equivalent to

Q12(A(k+i)+B(k+i)K)CTC(A(k+i)+B(k+i)K)Q12ymax2I

Then, the attack-angle constraint (Eq. 15) is equivalent to

αmax2I(A(k)Q+B(k)Y)TC1TC1(A(k)Q+B(k)Y)Q0C1=[ 0 1 0 ]

and the angular velocity constraint (Eq. 16) is equivalent to

ωymax2I(A(k)Q+B(k)Y)TC2TC2(A(k)Q+B(k)Y)Q0C2=1T[ 0 1 1 ]

Input constraint

Physical limitations inherent in process equipment invariably impose hard constraints on the manipulated variable δα1(k + i). The amplitude constraint can be represented as follows:

δα1(k+i)22umax2,i0

Setting K = YQ−1, then δα1 = Kx = YQ−1x. Then, Eq. (31) is equivalent to

maxi0δα1(k+i)22=maxi0YQ1x(k+i)22maxzTz1YQ1z22=λmax(Q12YTYQ12)

Using Schur complement, Eq. (32) is equivalent to

umax2IYYTQ0

Theorem 2 (Feasibility): Any feasible solution of the optimization in Theorem 1 at a time k is also feasible for all times t > k. Thus, if the optimization problem in Theorem 1 is feasible at a time k, then it is feasible for all times t > k.

The proof of Theorem 2 is similar to Korhare et al. (1996) and Jiang et al. (2017), which is omitted here.

Theorem 3 (Robust stability): The feasible receding horizon state feedback control law obtained from Theorem 1 robustly asymptotically stabilizes the closed-loop system.

Proof: Choose the Lyapunov function V(x(k)) = xT(k)P(k)x(k), then we have

V(x(k+i+1))V(x(k+i))=[xT(k+i)dT(k+i)]M(k)[xT(k+i)dT(k+i)]T+γ2dT(k+i)d(k+i)zT(k+i)z(k+i)

From stability of gray identification, we can know that compensation error is close to 0. Then, if z(k + i) is also close to 0, then the system stability is obtained. If not, the following inequality is satisfied:

γ2dT(k+i)d(k+i)zT(k+i)z(k+i)<0

Since the constraint (Eq. 24) guarantees M(k) < 0 and V(x(k + i + 1)) − V(x(k + i)) < 0, V(x(k)) is a strictly decreasing Lyapunov function for the closed loop.

Algorithm realization

The process of the algorithm proposed in this paper is

① At the sampling time k, based on the real-time input and output data, identify the characteristic model parameters of attack-angle channel f α1(k), f α2(k), g1(k), m11(k), m12(k), and Δα1(k) according to gray identification.

② Design compensator uc to compensate three-channel coupling motion information, unmodeled information, and disturbance Δα1(k).

③ Compute robust predictive control law δα1 = Kx based on Theorem 1. Then, apply total control u = δα1 + uc into the reentry vehicle.

④ Repeat process ①.

Control algorithms in other channels are similar to the algorithm for the attach channel and are omitted.

Simulation Analysis

According to a certain flight state of a plane-symmetry near space hypersonic vehicle, such as (α, β, φ) = (10, 0, 0), obtain the corresponding characteristic model parameters range (Zhang and Hu, 2012). Let the flight height is 40 km, speed is 10 Ma, and the range of aerodynamic coefficient (Liu and Chen, 2010) be kxx ∈ [−0.3585 −1.4341], kxz ∈ [−0.0101 −0.0152], kyy ∈ [−0.1060 −0.2945], kzx ∈ [−0.0127 −0.0506], kzz ∈ [−0.0253 −0.0380], kxa ∈ [40.5909 34.3395], kxr ∈ [3.1819 3.1729], kye ∈ [25.3225 22.7902], kza ∈ [2.7196 2.2845], and kzr ∈ [4.5363 4.1338].

Let the maximum of attitude angular velocity is ωxmax=ωymax=ωzmax=60s, the maximum of input control is 30, sampling time is 0.01 s.

Angle of attack control channel

α(k+1)=fα1α(k)+fα2α(k1)+g1δα(k)+Δαfα1[1.99682.0032],fα2[1.00320.9968],g1[2.27542.5363]×103.

Side slip angle control channel

β(k+1)=fβ1β(k)+fβ2β(k1)+g2δβ(k)+Δβfβ1[1.99652.0035],fβ2[1.00350.9965],g2[3.37604.0577]×104.

Roll angle control channel

φ(k+1)=fφ1φ(k)+fφ2α(k1)+g3δφ(k)+Δφft1[1.98232.0181],ft2[1.01810.9823],g33.30554.1957×103.

The requirement of the control system is to keep side slip angle around 0, the tracking accuracy of attack angle <0.5, simulation is taken based on the parameters above.

(1) Attack-angle tracking without coupling characteristic model disturbance

Assume the identified parameter is A=011.99681.0032, the initial value is 5, and the expected value of attack angle is 2, the coupling characteristic model disturbance Δα is 0. The weighting matrices in the robust predictive control are Q1 = diag(1,000, 1,000, 1,000) and R = 10. The simulation results are as follows: the tracking curve of attack angle is shown in Figure 2. The predictive control and gray identification guarantee the attack-angle tracking accuracy. The input and angular velocity constraints are shown in Figures 3 and 4, which meet the constraints, respectively. The results show that the prediction based on the characteristics model is feasible.

(2) Attack-angle tracking with coupling characteristic model disturbance

FIGURE 2
www.frontiersin.org

Figure 2. Tracking attack angle.

FIGURE 3
www.frontiersin.org

Figure 3. Control input.

FIGURE 4
www.frontiersin.org

Figure 4. Angular velocity.

Assume the identified parameter is A=011.99681.0032, the initial value is 5, and the expected value of attack angle is 2, the coupling characteristic model disturbance is Δα = m11β(k) + m12ϕ(k) + Δα1(k) is 1*sin(0.2t). The weighting matrices in the robust predictive control are Q1 = diag(1,000, 1,000, 1,000) and R = 10. The simulation results are as follows: the tracking curve of attack angle is shown in Figure 5. The predictive control and gray estimation compensation resist the disturbances of the characteristic model and model uncertainties. Figure 6 presents the input, from which we can see the constraint are satisfied.

(3) Attack-angle tracking with random noise and time-varying disturbance

FIGURE 5
www.frontiersin.org

Figure 5. Tracking attack angle.

FIGURE 6
www.frontiersin.org

Figure 6. Control input.

Assume the attack-angle initial value is 5, and the expected value is 2. The coupling characteristic model disturbance is Δα = m11β(k) + m21ϕ(k) + Δα1(k) is 1*sin(0.2t) in the sinusoidal change form 1*sin(0.2t). The attack angle is of 0 mean. In the case of the random state noise with the variance of 0.1, the characteristic model parameter of the attack angle A=011.890.976,g1=2.367×103 is obtained using the gray identification. Then, the robust MPC is adopted in the simulation. Figure 7 is the tracking curve of attack angle, which presents that the disturbance and model uncertainty in the characteristic model are resisted with the common effects of H2/H predictive control and gray estimation compensation. Figures 8 and 9 are the curves of rudder input and angular velocity, respectively, showing the satisfaction of input and output constraints with LMIs (Eqs 29–31).

(4) Robustness analysis

FIGURE 7
www.frontiersin.org

Figure 7. Tracking attack angle.

FIGURE 8
www.frontiersin.org

Figure 8. Control input.

FIGURE 9
www.frontiersin.org

Figure 9. Angular velocity.

If the identified parameters model parameter is A = [0 1; 2 −1], but the actual model parameter is A = [0 1; 1.9968 −1.0032], and the coupling characteristic model disturbance is Δα = m11β(k) + m12ϕ(k) + Δα1(k) is 1*sin(0.2t), the tracking curve of attack angle is shown in Figure 10, which can show the good tracking performance. The side slip angle keeps near the 0 according to the curve of side slip angle shown in Figure 11. The results show that the robust predictive control has the strong robust tracking ability.

FIGURE 10
www.frontiersin.org

Figure 10. Tracking attack angle.

FIGURE 11
www.frontiersin.org

Figure 11. Tacking side slip angle.

Conclusion

To design control for constrained, strong coupled, uncertain non-linear reentry hypersonic vehicles, the characteristic model is constructed and used as the predictive model in MPC. The original complicated reentry kinetics and dynamics are simplified by identification of the time-varying parameters in the characteristic model-based gray identification. The robust predictive control law is proposed based on the receding moving horizon optimization and LMI, which ensures stability of the closed-loop system and guarantees constraints satisfaction. Theoretical proof and simulation results confirm the effectiveness of the proposed algorithm.

Author Contributions

YP proposed the main idea. JZ and TS derived the controller and carried out simulations. JZ, TS, and YP wrote the paper.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors gratefully acknowledge the reviewers for their constructive and insightful comments for further improve the quality of this work.

Funding

This work was supported by the National Natural Science Foundation of China (No. 51379044, No. 51375210, No. 61503158, and No. 51375307).

References

Alicia, Z., and Andrea, S. (2009). “Combined reference governor and anti-windup design for constrained hypersonic vehicles models,” in AIAA (New York, NY: American Institute of Aeronautics and Astronautics), 2009–6283.

Google Scholar

Chiristopher, P., Morgan, B., and Ilya, K. (2013). “Model predictive control guidance with extended command governor inner-loop flight control for hypersonic vehicles,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit (New York, NY: American Institute of Aeronautics and Astronautics), 5028–5048.

Google Scholar

Gong, Y. L., and Wu, H. X. (2010). Characteristic model-based adaptive attitude control for hypersonic vehicle. J. Astronaut. 31, 2122–2128. doi: 10.3873/j.issn.1000.1328

CrossRef Full Text | Google Scholar

Harpold, J. (1979). Shuttle entry guidance. J. Astronaut. Sci. 28, 239–268.

Google Scholar

Jiang, W., Wang, H. L., Lu, J. H., Qin, W. W., and Cai, G. B. (2017). Nonfragile robust model predictive control for uncertain constrained time-delayed system with compensations. J Franklin. Inst. 354, 2832–2855. doi:10.1016/j.jfranklin.2017.02.001

CrossRef Full Text | Google Scholar

Korhare, M. V., Balakrishnans, V., and Moraris, M. (1996). Robust constrained model predictive control using linear matrix inequalities. Auromotica 32, 1361–1379.

Google Scholar

Li, Z. Y., Zhou, W. J., and Liu, H. (2016). Nonlinear robust control of hypersonic aircrafts with interactions between flight dynamics and propulsion systems. ISA Trans. 64, 1–11. doi:10.1016/j.isatra.2016.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., and Chen, R. W. (2010). Modeling algorithm using accumulated generating operation and weighted product. Syst. Eng. Electron. 32, 976–979. doi:10.3969/j.issn.1001-506X

CrossRef Full Text | Google Scholar

Lu, P. (1999). Regulation about time-varying trajectories: precision entry guidance illustrated. J. Guid. Control Dyn. 22, 784–790. doi:10.2514/2.4479

CrossRef Full Text | Google Scholar

Lu, P. (2014). Entry Guidance: a unified method. J. Guid. Control Dyn. 37, 713–728. doi:10.2514/1.62605

CrossRef Full Text | Google Scholar

Patience, E., and Orukpe, I. M. (2007). “Model predictive control based on mixed H2/Hcontrol approach,” in Proc of the 2007 American Control Conf (New York, NY: IEEE Computer Society), 6147–6150.

Google Scholar

Qin, W. W., He, B., Liu, G., and Zhao, P. T. (2016). Robust model predictive tracking control of hypersonic vehicles in the presence of actuator constraints and input delays. J. Franklin Inst. 353, 4351–4367. doi:10.1016/j.jfranklin.2016.08.007

CrossRef Full Text | Google Scholar

Shen, Z. J., and Lu, P. (2003). On board generation of three-dimensional constrained entry trajectories. J. Guid. Control Dyn. 26, 111–121. doi:10.2514/2.5021

CrossRef Full Text | Google Scholar

Van, W. R. (2006). Combined feedback linearization and constrained model predictive control for entry flight. J. Guid. Control Dyn. 29, 427–434. doi:10.2514/1.14511

CrossRef Full Text | Google Scholar

Wu, G. H., and Meng, X. Y. (2016). Nonlinear disturbance observer based robust backstepping control for a flexible air-breathing hypersonic vehicle. Aerosp. Sci. Technol. 54, 174–182. doi:10.1016/j.ast.2016.04.018

CrossRef Full Text | Google Scholar

Wu, H., and Meng, B. (2009). Review of the control research of hypersonic vehicle. Adv. Mech. 39, 756–765. doi:10.3321/j.issn:1000-0992

CrossRef Full Text | Google Scholar

Wu, H. X., Hu, J., and Xie, Y. C. (2007). Characteristic model-based all-coefficient adaptive control method and its applications. IEEE Trans. Syst. Man Cybern. Syst. 37, 213–221. doi:10.1109/TSMCC.2006.887004

CrossRef Full Text | Google Scholar

Wu, H. X., Hu, J., and Xie, Y. C. (2009). Characteristic Model-Based Intelligent Adaptive Control. Beijing: China Science and Technology Publishing.

Google Scholar

Xu, B., Sun, F., Liu, H., and Ren, J. (2012). Adaptive Kriging controller design for hypersonic flight vehicle via back-stepping. IET Control Theory Appl. 4, 487–497. doi:10.1049/iet-cta.2011.0026

CrossRef Full Text | Google Scholar

Xu, L. J. (2012). Hypersonic vehicle control based on the all-channel-coupled characteristic model. Aerosp. Control Appl. 38, 12–18. doi:10.3969/j.issn.1674-1579

CrossRef Full Text | Google Scholar

Zhang, J., and Sun, T. (2016). Disturbance observer-based sliding manifold predictive control for reentry hypersonic vehicles with multi-constraint. Proc. Inst. Mech. Eng. 230, 485–495. doi:10.1177/0954410015593564

CrossRef Full Text | Google Scholar

Zhang, Z., and Hu, J. (2012). Stability analysis of a hypersonic vehicle controlled by the characteristic model based adaptive controller. Sci. China Inf. Sci. 55, 2243–2256. doi:10.1007/s11432-011-4455-9

CrossRef Full Text | Google Scholar

Keywords: hypersonic vehicle, reentry, characteristic model, system constraints, robust model predictive control, gray identification

Citation: Zhang J, Sun T and Pan Y (2017) Characteristic Model-Based Robust Model Predictive Control for Hypersonic Vehicles with Constraints. Front. Robot. AI 4:23. doi: 10.3389/frobt.2017.00023

Received: 09 October 2016; Accepted: 15 May 2017;
Published: 13 June 2017

Edited by:

Gianluca Antonelli, University of Cassino, Italy

Reviewed by:

Yanghai Nan, Glasgow Caledonian University, United Kingdom
Francesco Pierri, University of Basilicata, Italy

Copyright: © 2017 Zhang, Sun and Pan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yongping Pan, biepany@nus.edu.sg