Open access

Predictive Control for Active Model and Its Applications on Unmanned Helicopters

Written By

Dalei Song, Juntong Qi, Jianda Han and Guangjun Liu

Submitted: October 22nd, 2010 Published: July 5th, 2011

DOI: 10.5772/17716

Chapter metrics overview

2,233 Chapter Downloads

View Full Metrics

1. Introduction

Unmanned helicopters are increasingly popular platforms for unmanned aerial vehicles (UAVs). With the abilities such as hovering, taking off and landing vertically, unmanned helicopters extend the potential applications of UAVs. However, due to the complex mechanism and complicated aero-flow during flight, it is almost impossible to accurately model the dynamics of an unmanned helicopter in full flight envelope, and the significant model uncertainties associated with a nominal model may degrade the performance and even stability of an onboard controller.

Due to the difficulty in obtaining a high fidelity full envelope model, the multi-mode modeling technique has been proposed for rotor aircrafts, such as tilt-rotor aircraft XV-15 [1], helicopter BO-105 [2], UH-60 [3], R-50 [4] and X-Cell [5]. The mode-dependent model, which is identified and simplified according to a specific flight mode, such as hovering, cruising, taking off and landing, can be used for control design for the corresponding flight mode. However, the mode-dependent control suffers from at least two problems: one is the difficulty in accommodating the mode transition dynamics, and the other is the compensation of the ‘model drift’ due to flight dynamics change within one particular mode. Up to now, for the purpose of practical implementation, the mode transition problem can be partially dealt with by limiting the mode switching conditions [6], e.g., mode change is made through hovering mode.

Robust and adaptive control techniques [7-8], on the other hand, have been used to deal with the ‘model-shift’ within a flight mode. However, such control schemes normally need to know the boundary of internal and external uncertainties and relative noise distribution, which are difficult to identify accurately for a helicopter in full flight envelope. Although online identification technology can be used to obtain the real-time dynamics and disturbance, it is a large burden for the flight computer to reconstruct the robust controllers and reach the requested control period (>50Hz) for sampling and actuating due to the complex calculation of the robust/adaptive optimization process [9-10] and the strict weight limits of micro flight computers.

Besides the model uncertainties, another critical problem that limits the control performance of a helicopter is the time delay between the actuator command and the generation of relative aerodynamic force/torque [11], which will be called aerodynamics-delay/time-delay in the following sections. Normally, this time delay may cause reduced feedback gain of a model-based controller and result in poor robustness [12-13], i.e., sensitive to disturbances.

In recent years, the encouraging achievement in sequential estimation makes it an important direction for online modeling and model-reference control [14]. Among stochastic estimations, the most popular one is the Kalman-type filters (KFs) [15, 16, and 17]. Although widely used, the KFs suffer from sensitivity to bias and divergence in the estimates, relying on assumptions on statistic distribution such as white noise and known mean or covariance for optimal estimation. In many cases, it is more practical to assume that the noises or uncertainties are unknown but bounded (UBB). In view of this, the set-membership filter (SMF), which computes a compact feasible set in which the true state or parameter lies only under the UBB noise assumption, provides an attractive alternative [18-19].

On the control issue, model predictive control (MPC) can compensate for the aerodynamics delay and does not require a high accuracy reference nonlinear model [20]. Among these methods, linear generalized predictive control (GPC) has become one of the most popular MPC methods in industry and academia. However, the normal GPC is sensitive to process noise and model errors [21], which are unknown but bounded for helicopters when sudden ‘mode change’ happen and model-drift in full flight envelope. This makes the prediction biased, and results in the non-optimal process of controller solving.

In this paper, for realizing the coupling control of unmanned helicopters in full flight envelope, an active modeling based controller is developed based on a modified generalized predictive control and adaptive set-membership filter estimation (ASMF). The time varying model error and its boundary are estimated by the adaptive set-member filter, which is first proposed in [19]. Incremental prediction process and dimension reduction method is embedded into traditional GPC, which can decreases the computation burden and maintain prediction unbiased when ‘mode change’ happens. Based on this active estimation and the modified GPC controller, a novel optimal strategy for on-line compensation of model error is developed. Thus, aggressive flight can be achieved only based on the hovering model with time-delay terms. Using the identified hovering dynamics model as nominal model for controller, flight experiments have been conducted to test the performance of the proposed controller in full flight envelope on our UAV platform, and experimental results have demonstrated the effectiveness of the proposed method.


2. Active model based control scheme and reference model of a helicopter

Fig. 1 illustrates the active model based control scheme. The error between the reference model and the actual dynamics of the controlled plant is estimated by an on-line modeling strategy. The control, which is designed according to the reference model, should be able to compensate the estimated model error and it in real time. In the followings of this paper, we use the ASMF as the active modeling algorithm and the modified GPC as the control.

For normal missions of an unmanned helicopter, the flight modes include hovering (velocity under 5m/s), cruising (velocity above 5m/s), taking off and landing (distance to the ground is below 3m while significant ground effect exists) and the transitions among these modes. A reference model is typically obtained by linearizing the nonlinear dynamics of a helicopter at one flying mode. The model errors from linearization, external disturbance, simplification, and un-modeled dynamics can be considered as additional process noise [22]. Thus, a linearized state-space model for helicopter dynamics in full flight envelope can be formulated as

Figure 1.

The scheme of active model based control


where XR13 is the state, including 3-axis velocity, pitch and roll angle, 3-axis angle rate, flapping angles of main rotor and stabilizer bar, and the feedback of yaw gyro. YtR8is the output, including 3-axis velocity, pitch and roll angle and 3-axis angle rate, A0and B0 contain parameters that can be identified in different flight modes, and we use them to describe the parameters in hovering mode. UR4is the control input vector. CR13×8is the output matrix, kRis the time-delay for the driving system. The detail of building the nominal model and physical meaning of parameters is explanted in Appendix A.

To describe the dynamics change, in equation (1), here, we introduce f(Xt,X˙t,Wt)R13 to represent the time varying model error in full flight envelope, and WtR13 is the process noise.

The following two sections, based on model (1) will describe the way to estimate f(Xt,X˙t,Wt) and to compensate for model errors from process noise, parameters change, control delay and flight mode change in real applications.


3. ASMF based active model error estimation

As illustrated in Fig.1, adopting the active modeling process to get the model error f and system state X is the basis for elimination of the model error. Controller can only work based on nominal model and feedback of state and model error from active modeling process. In this section, the active modeling process is built based on an adaptive set-membership filter (ASMF) [19] since the UBB process noise.

First, we must obtain the reference equation for estimation. Compared with the sampling frequency (often >50Hz for flight control) of the control system, the model error f(X,X˙,W) can be considered as a slow-varying vector, which means


where ftis the sampling value of f(X,X˙,W)at sampling time t, and htis the assumed unknown but bounded (UBB) process noise.

Let the extended sampling state


Then, we can obtain the discrete equation from Eq. (1) as


whereAda=(AdBf013×13I13×13), Bda=(Bd013×4), Cda=(Cd08×13), Wta=(WtThtT)T, Bf=I13×13and ft is a 13×1vector for model errors. Here, t is the sampling time, Im×mis the m×m unit matrix and 0m×n is the m×n zero matrix. {Ad,Bd,Cd}is the discrete expression of system{A0,B0,C}. Here, time-delay k is ignored during the estimate process, and the compensation method will be discussed in the next part on modified GPC.

The model error i.e., f in Eq. (1), comes from the linearization while neglecting the coupling dynamics and uncertainties, and also the A0 and B0 because they are identified with respect to a specific flight mode, here hovering mode is selected as nominal flight mode since easy identification. Therefore, both the model error and the process noise Wa are vehicle dynamics and flight states dependent, and do the following assumption

Assumption:Wadoes not necessarily have a normal distribution.

Thus, the Kalman type filter cannot be applied, and adaptive set-membership filter, which is developed for UUB process noise and can get the uncertain boundaries of the states, is considered to estimate the states and model errors here.

In this section we only present the result of ASMF and please refer to [19] for the details about ASMF. With respect to Eq. (2), we can build the adaptive set-membership filter as Eq. (3), where Qa and Ra are the initial elliptical boundary of process and measurement noise respectively, rmis the maximum eigenvalue of R, pmis the maximum eigenvalue ofCdaPt|t1CdaT, Tr(•) is the trace of a matrix, δtand βt are the adaptive parameters of the filter. We can also obtain the boundary of the ith element X^iaof extended state X^t|taas(X^iaPii,X^ia+Pii), where Pii is the i-th diagonal element of matrix P.


4. Modified GPC for unmanned helicopters

To eliminate the negative influence of model errors and control delay in flight, besides the active estimation algorithm like ASMF that does not require a normal distribution assumption, an effective control algorithm has to be designed according to the reference model of Eq. (1) while adopting the on-line estimation of f as compensation.

We describe the normal GPC in Section 4.1, and then, the modified scheme is proposed in Section 4.2 & 4.3 to eliminate the negative influence of model errors in real applications.

4.1. Preliminary work for generalized predictive control

Generally, for a linear system with actuator time delay like,


where XtRn×1is the system state vector at sampling time t, ytRl×1is the output vector, utRm×1is the control input vector, k is the actuators’ time-delay and Wt is process noise; traditional Generalized Predictive Control (GPC) [23] can be designed as:

Step I: Make prediction

Firstly, for the case that predictive step i is less than time-delay k (i.e., the time instant that system behavior cannot be regulated through current and future control action), prediction can be denoted as following equation,


whereX^t+i|t is the prediction state at time t+i, the superscript 1 denotes that the part of predicted variable that is independent of the current and future’s control actions.

Secondly, for the case that prediction step i is larger than the time delay k,


where p is the prediction range; similarly, X^t+k+i1|t1denotes the sub-variable of X^t+k+i1|tthat is independent of the current and future’s control actions.

Step II: Receding horizon optimization

After making prediction, the control vector can be obtained by minimize the following cost function:


And the optimal control inputs can be denoted as,


where G0is the predictive matrix, Xtvis the predictive state vector,Xt1 is the known vector insideXtv, λis the weight of control input, and Rtxis the reference of system states. The detailed definition of these matrixes can be referenced in [23].

Step III: Control implementation

The first element of vector Ut*is used as the control to the real plant. After that, go back to step I at the next time instant.

However, with application to the unmanned helicopters, this kind of GPC algorithm has the following three disadvantages, which will be solved in the next two sections:

  1. It cannot reject the influence of working mode changes, i.e., if


where (x0,u0)is the current operation point, which cannot be ensured on-line, π(x0,u0)is the valid range for model linearization and xtis the absolute state at time t, utis the absolute control input at time t. The biased prediction, due to the changing operation point(x0,u0), will bring steady errors for velocity tracking.

  1. Normal GPC is sensitive to mismatch of the nominal model, which means slow change in parameters (Ad,Bd)may result in prediction error and unstable control.

  2. The transient model errors of the nominal model from external disturbance, estimated by ASMF, cannot be eliminated. And this will also result in the non-minimum variance and the instability of the closed control loop.

4.2. Stationary increment predictive control

To reject the influence of working mode change and sensitivity to nominal parameters change in real application, i.e. the problem 1) and 2) in Section 4.1, we assume that the process noiseWt’s increment in Eq. (4) is a stationary random process, which means


is normal distribution. Where Δ=1q1 is the difference operator; q1is one-step delay factor. Thus, Eq. (4) can be rewritten as follows,




if behavior prediction is made based on Eq. (11), only the absolute statextand control inputut, which can be measured or estimated directly from sensors, are used and the current operation point (x0,u0) disappears in prediction. Thus, the problem of biased prediction due to changing of working point, i.e., problem 1), can be solved.

Otherwise, according to the process of traditional GPC, the set-point Rtx must be obtained for every prediction step, and this is often set as current reference states. However, for helicopter system, only measurable outputs are cared, such as position, velocity and etc; and the internal states, such as rotor’s pitch angle and yaw gyro’s feedback and so on, are coupled with the measurable states/outputs, and cannot be set independently. Others, this reference input often comes from position track planning, which changes quickly for flight and often cause a step-like signal for tracking. To avoid the step signal reference tracking, which is dangerous for unmanned helicopter system, we use a low pass filter to calculate the set-point inputs of the output in the future i-th step, i=1, …, p.

Let SPtRl×1 be the set-point input at time t, then we have


where α is the cut-off frequency of the filter, the initial valuert+k=y^t+k|t, rt+k+iis the i-th set-point input, and y^t+k|tis the estimate of output at time t+k.

Thus, the set-point problem is solved and the output prediction can be implanted based on increment model (11) as follows:

When the prediction step i is less than time-delay k,


When the prediction step is larger than time-delay k, let




Hence, the above problem 1), which comes from working mode change, is solved because x0 disappears in predictive equation (14).

We can obtain the following prediction matrix for the output, which is often cared in helicopter tracking problem, from Eq. (12) and (13):


where Yt1 is the known part of p steps’ prediction, which cannot be influenced by current control input, and matrix G has the following form:


Compared with the normal GPC, the prediction of SIPC has better characteristics that can be described by the following theorem, which solves the above problem 2) in Section IV.A.

Theorem: for nominal model (11), when the nominal model parameters(Ad,Bd)change into(Adr,Bdr).

  1. M,N>0R, let the matrix norms satisfy


  2. DefineRmax{}is the operator for the maximum of eigenvalue of matrix.

Thus, if


Then, the state prediction obtained by Eqs. (13-14) maintains unbiased, and the characteristic is also guaranteed in traditional GPC conditions, i.e. Eq. (4), where Wtis normal distribution.

Proof: See Appendix B.

In Eq. (14),ΔU , including p control inputs, need to be optimized, while only the first one is used for control. This will occupy a great deal of computation resource and result in very low computational efficiency, especially with respect to the fast applications.

In order to reduce the computational burden of Eq. (14), we propose here a ‘step plan’ technique,


where β is an m×m diagonal matrix presenting the length of one step, which will be a parameter to be selected. Then, we can simplify Eq. (14) by only calculating the unknown control, which has smaller dimensions.


where Im×m is an m×m unit matrix. Thus, the number of the unknown control input vector (from current time t to the future time t+p-1) is reduced from p to 1, and the dimension of predictive matrix is changed from pl×pmtopl×m. This reduction brings low computer memory consuming and simplifies the receding horizon optimization in the following calculation.

To complete the horizon optimization and obtain the control input, the cost function of the stationary increment predictive control is designed as:


whereRt=(rt+k+1Trt+k+2T...rt+k+pT)T, WRlp×lpis the weight matrix for tracking error, and λRm×mis the weight matrix of the control increment.

In order to minimize the cost function of Eq. (19), we can calculate the control vector as follows:


where Kf=(G2TWG2+λ)1G2TW can be completed offline.

Consequently, the proposed stationary increment predictive controller (SIPC) can be designed as followings.

Step I: Make increment prediction

Based on the current and history measure value, use Eqs. (13-15) to obtain the prediction for future output Y^tand initial plan point


Step II: Plan for the set-point input

Use Eq. (12) to plan the future set-points, and obtain


Step III: Receding horizon optimization

Calculate the control incrementΔut, based on Eq. (20).

Step IV: Control implementation

Current control inputut=ut1+Δut, which is used as the control to the real plant. After that, go back to step I at the next time instant.

Thus, for real implementation, only the prediction of Eq. (13-15), the intenerating of Eq. (12), and the control law (20) need to be calculated online, thus the real time computation load, and steady tracking error are both reduced greatly compared with GPC, and the real test in section V has shown its feasibility.

The model error, problem 3), will be compensated by an online optimal strategy, which will be described later.

4.3. Optimal strategy for model error compensation

In order to compensate the model error in Eq. (1), the control vector has to match the following equation, which can be directly obtained from Eq. (1):


where Ut0 is the control vector need to be calculated by the predictive controller in section 4.2, designed based on the original model (1) without the model error f.

The control input at sampling time t cannot be solved directly from Eq. (21), because:

  1. Eq. (21) is difficult to be implemented because the dimension ofUt is less than that offt. Thus, only the approximate solution can be obtained with respect to (21);

  2. ftis actually an uncertainty set, an static optimal problem must be considered.

Thus, we introduce the following cost function with quadratic form to solve the above problem 1).


where H is a weight matrix, which can be selected.

On the other hand, ftis obtained from the ASMF algorithm introduced in section III, thus its convergence is very important for the validity of the whole controller. Actually, the convergence of ASMF algorithm is also influenced by the control actionUt. This is because the stability of the ASMF can be represented by the filter parameterδt, while δt in Eq. (3) can be rewritten as follows,


In [19], it has been shown the stability of the ASMF can be represented by the filter parameterδt, i.e., the ASMF is stable whenδt>0.

Firstly, define


Thus, from Eq. (23), in order to maintainδt+1>0, the maximum value of Jtδ(Ut,Yt+1) with respect to X^t|ta should be less than or equal to 1, i.e.,


In general, larger δt often means more rapid convergence of ASMF algorithm. That is, we should select an Ut to make Jtδ*(Ut,Yt+1) small as far as possible, that is,


We introduce the following cost function Jt(Ut) with consideration of both (22) and (25) at the same time:


where α=1δtR are the positive definite weight matrix. To minimizeJt(Ut), consideringJt(Ut)>0, the control can be obtained atJt(Ut)Ut=0, i.e.,





Here H can be selected asH=δtCdTCd. Thus, we can obtain the optimal control that minimizes Jt(Ut) as:


For the unknown measurement at time t+1 in Eq. (24), we consider that the control system is stable, so,Yt+1Δ(Yt). Here, Δ(Yt)is the elliptical domain ofYt. Because Jtδ(Ut,Yt+1) in Eq. (25) is positive definite, its maximum value point must be on the boundary, which can be estimated by the ASMF. Thus, we first define array Stito include the estimate of the i-th element’s two boundary endpoints as


where Yti is the i-th element in the vectorYt, Y^t+1iis the corresponding outputYt+1’s endpoints estimation. For setSti, i{1,2,...,8}and h is 0 or 1 for every i, ||iis the operator for absolute value of the i-th element in vector, and the function Col{j} is defined as follows:


Then, we define a set Stto describe all possible endpoint vector of theYt+1 as


where Y^t+1EPis the possible endpoint (EP) for output Yt+1at next sampling time t+1.

Thus, the proposed active modeling based predictive controller can be implemented by using the following steps:

Step I: Make increment prediction

Based on the current estimated stateX^t|ta, use the stationary increment predictive controller, as in section 4.2, to obtain the nominal control inputUt0;

Step II: Model error estimation and elimination

Based onUt0, compute the optimal control input*Ut:

Estimate the values and boundaries of state Xt and model errorft, using ASMF in (3);

Calculate the corresponding Ut(Y^t+1EP) for every Y^t+1EP in setSt by Eq. (29);

For every Ut(Y^t+1EP)in step 1), use Eq. (24) to obtain the maximum of functionJtδ(Ut(Y^t+1EP),Y^t+1EP), and get the *Y^t+1EPto let


The corresponding Ut(*Y^t+1EP) is the optimal control *Ut at time t, i.e.*Ut=Ut(*Y^t+1EP).

Step III: Receding horizon strategy

Go back to step I at the next time instant.


5. Flight test

5.1. Flight test platform

All flight tests are conducted on the Servoheli-40 setup, which was developed in the State Key Laboratory, SIACAS. It is equipped with a 3-axis gyro, a 3-axis accelerometer, a compass and a GPS. The sensory data can be sampled and stored into an SD card through an onboard DSP. Tab.1 shows the physical characteristics of SERVOHELI-40 small-size helicopter. More details of this experimental platform can be found in [24].

Figure 2.

SERVOHELI-40 small-size helicopter platform

Main rotor diameter2.15m
Stabilizer bar diameter0.75m
Rotor speed1450rpm
Dry weight20kg
Engine2-stroke, air cooled
Flight time45 min

Table 1.

Physical characteristics of SERVOHELI-40 small-size helicopter

5.2. Experiment for the verification of model error estimate when mode-change

We use the identified hovering parameters, through frequency estimate [25], as the nominal model for hovering dynamics of the ServoHeli-40 platform. The model accuracy is verified in hovering mode (speed less than 3m/s) and cruising mode (speed more than 5m/s), the results for lateral velocity are shown in Fig.3a.

Fig. 3 further shows the model difference due to mode change, where the red lines are the results calculated by the identified model with the inputs of hovering and cruising actuations, respectively, and blue lines are the measurements of the onboard sensors. Comparison shows that the hovering model outputs match the hovering state closely, but clear differences occur while being compared to the cruising state, even though the cruising actuations are used as the model inputs. This is the model error when flight mode is changed.

Figure 3.

Model difference due to mode change: (a) hovering conditions; (b) cruising conditions

To verify the accuracy of the estimate of the model error, described in Fig.3, the following experiment is designed:

  1. Actuate the longitudinal control loop to keep the speed more than 5 meter per second;

  2. Get the lateral model error value and boundaries through ASMF, and add them to the hovering model we built above;

  3. Compare the model output before and after compensation for model error.

This process of experiment can be described by Fig.4, and the results are shown in Fig.5. Fig.5a shows that model output (red line) cannot describe the cruising dynamics due to the model error when ‘mode-change’, similar with Fig.3b; however, after compensation, shown in Fig.5b, the model output (red line) is very close with real cruising dynamics (blue line), and the uncertain boundaries can include the changing lateral speed, which mean that the proposed estimation method can obtain the model error and range accurately by ASMF when mode-change.

Figure 4.

The experiment process for model-error estimate

Figure 5.

Model output before/after compensation: (a) before compensation; (b) after compensation

5.3. Flight experiment for the comparison of GPC SIPC and AMSIPC when sudden mode-change

In Section 5.2, the model-error occurrence and the accuracy of the proposed method for estimation are verified. So, the next is the performance of the proposed controller in real flight. In this section, the performance of the modified GPC (Generalized Predictive Control, designed in Section 4.1), SIPC (Stationary Increment Predictive Control, designed in Section 4.2) and AMSIPC (Active Modeling Based Stationary Increment Predictive Control, designed in Section 4.3), are tested in sudden mode-change, and are compared with each other on the ServoHeli-40 test-bed. To complete this mission, the following experimental process is designed:

  1. Using large and step-like reference velocity, red line in Fig.6-8, input it to longitudinal loop, lateral loop and vertical loop;

  2. Based on the same inputted reference velocity, using the 3 types of control method, GPC, SIPC and AMSIPC to actuate the helicopter to change flight mode quickly;

  3. Record the data of position, velocity and reference speed for the 3 control loops, and obtain reference position by integrating the reference speed;

  4. Compare errors of velocity and position tracking of GPC, SIPC and AMSIPC, executively, in this sudden mode-change flight.

GPC, SIPC and AMSIPC are all tested in the same flight conditions, and the comparison results are shown in Figs. 6-8. We use the identified parameters in Section 5.2 to build the nominal model, based on the model structure in Appendix A, and parameters’ selection in Appendix C for controllers

It can be seen that, when the helicopter increases its longitudinal velocity and changes flight mode from hovering to cruising, GPC (brown line) has a steady velocity error and increasing position error because of the model errors. SIPC (blue line) has a smaller velocity error because it uses increment model to reject the influence of the changing operation point and dynamics’ slow change during the flight. The prediction is unbiased and obtains better tracking performance, which is verified by Theorem. However, the increment model may enlarge the model errors due to the uncertain parameters and sensor/process noises, resulting in the oscillations in the constant velocity period (clearly seen in Fig.6&7) because the error of its prediction is only unbiased, but not minimum variance. While for AMSIPC (green line), because the model error, which makes the predictive process non-minimum variance, has

Figure 6.

Longitudinal tracking results: (a) velocity; (b) position error (<50s hovering, >50s cruising)

Figure 7.

Lateral tracking results: (a) velocity; (b) position error (25s~80s cruising, others hovering)

Figure 8.

Vertical tracking results: (a) velocity; (b) position error (<5s hovering; >5s cruising)

been online estimated by the ASMF and compensated by the strategy in section 4.3, the proposed AMSIPC successfully reduces velocity oscillations and tracking errors together.


6. Conclusion

An active model based predictive control scheme was proposed in this paper to compensate model error due to flight mode change and model uncertainties, and realize full flight envelope control without multi-mode models and mode-dependent controls.

The ASMF was adopted as an active modeling technique to online estimate the error between reference model and real dynamics. Experimental results have demonstrated that the ASMF successfully estimated the model error even though it is both helicopter dynamics and flight-state dependent.In order to overcome the aerodynamics time-delay, also with the active estimation for optimal compensation, an active modeling based stationary increment predictive controller was designed and analyzed.

The proposed control scheme was implemented on our developed ServoHeli-40 unmanned helicopter. Experimental results have demonstrated clear improvements over the normal GPC without active modeling enhancement when sudden mode-change happens.

It should be noted that, at present, we have only tested the control scheme with respect to the flight mode change from hovering to cruising, and vice versa. Further mode change conditions will be flight-tested in near future.


7. Appendix

A. Helicopter dynamics

A helicopter in flight is free to simultaneously rotate and translate in six degrees of freedom. Fig. 9 shows the helicopter variables in a body-fixed frame with origin at the vehicle’s center of gravity.

Figure 9.

Helicopter with its body-fixed reference frame

Ref.[18] developed a semi-decoupled model for small-size helicopter, i.e.,


where δu, δv, δw are longitudinal, lateral and vertical velocity, δp, δq, δr are roll, pitch and yaw angle rates, δφ and δθ are the angles of roll and pitch, respectively, a and b are the first harmonic flapping angle of main rotor, c and d are the first harmonic flapping angle of stabilizer bar, δrfbis the feedback control value of the angular rate gyro, δlatis the lateral control input, δlonis the longitudinal control input, δpedis the yawing control input, and δcol is the vertical control input. All the symbols except gravity acceleration g inAlon, Alat, Ayawheave, Blon,Blat and Byawheave are unknown parameters to be identified. Thus, all of the states and control inputs in (33), (34) and (35) are physically meaningful and defined in body-axis.

B. Proof for the predictive theorem


Assume the real dynamics is described as:


which is different from the reference model of Eq. (11). In Eq. (36), Xtis system state, Adris the system matrix, Bdris the control matrix, Utis control input, Wtis process noise. The one-step prediction, according to Eq. (36), can be obtained by Eq. (13-14),




According to condition 1) and 2), prediction is bounded, then,


and, when the system of Eq. (36) works around a working point in steady state, the mean value of control inputs and states should be constant, so we can obtain:


Eq. (39) indicates that the one step prediction of Eq. (37) is unbiased.

Assuming that prediction at time i-1 is unbiased, i.e..


for the prediction at time i, there is


Therefore, the prediction at time i is also unbiased.

C. Parameters’ selection for estimate and control in flight experiment

  1. For Modeling

The identification results for hovering dynamics are listed in Tab.D-1.

Longitudinal LoopLateral LoopVertical Loop
Mu-1.258Lv15.84Yaw Loop

Table 2.

The parameters of hovering model

  1. For ASMF



where Im×m is the m×m unit matrix and 0m×n is the m×n zero matrix.

  1. For GPC

p=10, γ=2.32I40×40,

  1. For SIPC

p=10, γ=2.32I4×4,


W=I80×80, k=10,

  1. For AMSIPC

p=10, γ=2.32I4×4,


W=I80×80, k=10, β=0.8I4×4,



  1. 1. Tischler M.B., “Frequency-domain Identification 15 XV-15 Tilt-rotor Aircraft Dynamics in Hovering Flight,” Journal of the American Helicopter Society, 30 2), 3848 , 1985.
  2. 2. TischlerM. B.CauffmanM. G.Frequency-Response“.Methodfor.RotorcraftSystem.IdentificationFlight.Applicationto. B. O. 05 Coupled Rotor/Fuselage Dynamics,” Journal of the American Helicopter Society, 37 3), 317 , 1992.
  3. 3. FletcherJ. W.Identification“. 60 UH-60 Stability Derivative Models in Hover from Flight Test Data,” Journal of the American Helicopter Society, 40 1), 820 , 1995.
  4. 4. MettlerB.TischlerM. B.KanadeT.System“.Identificationof.Small-SizeUnmanned.HelicopterDynamics,”.AmericanHelicopter.Society 55 Annual Forum Proceedings, 2 17061717 , Montreal, Quebec, Canada, May 25-27, 1999.
  5. 5. GavriletsV.MetlterB.FeronE.Nonlinear“.Modelfor. a.Small-scaleAcrobatic.Helicopter,”Proceedings.ofthe.AmericanInstitute.ofAeronautics.GuidanceNavigation.ControlConference. 8 8 Montreal, Quebec, Canada, August 6-9, 2001.
  6. 6. MassimilianoM.ValerioS.Full“. A.EnvelopeSmall.CommercialAircraft.FlightControl.DesignUsing.Multivariable-IntegralProportional.Control,”I. E. E. E.Transactionson.ControlSystems.TechnologyVol. 161 (1), 169176 , January, 2008.
  7. 7. VoorsluijsM.MulderA.Parameter-dependent.robustcontrol.fora.rotorcraftU. A. V.GuidanceA. I. A. A.NavigationControlConference.Exhibit 111 111 San Francisco, California, USA, August 15-18, 2005.
  8. 8. Bijnens B., Chu Q.P. and Voorsluijs M., "Adaptive feedback linearization flight control for a helicopter UAV," AIAA Guidance, Navigation, and Control Conference and Exhibit, 110 110 San Francisco, California, USA, August 15-18, 2005.
  9. 9. Kahveci N.E., Ioannou P.A., Mirmirani M.D., “Adaptive LQ Control With Anti-Windup Augmentation to Optimize UAV Performance in Autonomous Soaring Applications,” IEEE Transactions on Control Systems Technology, Vol. 164 691707 , 2008
  10. 10. MacKunis W., Wilcox Z.D., Kaiser M.K., Dixon W.E., “Global Adaptive Output Feedback Tracking Control of an Unmanned Aerial Vehicle,” IEEE Transactions on Control Systems Technology, Vol. 186 13901397 , 2010.
  11. 11. CummingsM. L.MitchellP. J.Predicting“.ControllerCapacity.inSupervisory.Controlof.MultipleU. A.VsSystems,”I. E. E. E.Transactionson.ManCyberneticsPart. A.SystemsHumansVol. 382 451460 , 2008.
  12. 12. Jiang X., Han Q.L., “On guaranteed cost fuzzy control for nonlinear systems with interval time-varying delay,” Control Theory & Applications, IET, Vol. 16 17001710 , 2007.
  13. 13. Natori K., Oboe R., Ohnishi, K., “Stability Analysis and Practical Design Procedure of Time Delayed Control Systems With Communication Disturbance Observer,” IEEE Transactions on Industrial Informatics, Vol. 43 185197 , 2008.
  14. 14. Haykin and De Freitas N., “Special Issue on Sequential State Estimation,” Proceedings of the IEEE, 92 923 423574 , 2004.
  15. 15. LerroD.Bar-ShalomY. K.Tracking”.withDebiased.ConsistentConverted.vsMeasurementsE. K. EKF,” IEEE Transactions on Aerosp. Electron.System, 29 10151022 , 1993
  16. 16. JulierS.UhlmannJ.Unscented“.filteringnonlinearestimation,”.Proceedingsof.theI. E. E. E.Vol 923 401422 , 2004.
  17. 17. Song Q., Jiang Z., and Han J. D., “UKF-Based Active Model and Adaptive Inverse Dynamics Control for Mobile Robot,” IEEE International Conference on Robotics and Automation, 2007.
  18. 18. ShammaJ. S.TuK. Y.Approximate“,”I. E. E. E.Transactionson.AutomaticControl.Vol 425 648658 , 1997.
  19. 19. Zhou B., Han J.D. and Liu G., “A UD factorization-based nonlinear adaptive set-membership filter for ellipsoidal estimation,” International Journal of Robust and Nonlinear Control, 18 16), 15131531 , November 10, 2007.
  20. 20. Scholte E., Campbell M.E., “Robust Nonlinear Model Predictive Control With Partial State Information,” Control Systems Technology, IEEE Transactions on, Vol. 164 636651 , 2008.
  21. 21. DingB. C.XiY. G.Synthesis“. A.Approachof.On-lineConstrained.RobustModel.PredictiveControl.”.Automatica 40 401 163167 , 2004.
  22. 22. CrassidisJ. L.Robust“.Controlof.NonlinearSystems.Using-ErrorModel.ControlSynthesis,”.Journalof.guidancecontrol.dynamicsVol. 224 (4), 595601 , 1999.
  23. 23. GregorK.IgorS.Tracking-error“.Model-basedPredictive.Controlfor.MobileRobots.inreal.time.”Robotics.AutonomousSystems. 55 55 7 460469 , 2007.
  24. 24. Qi J.T., Song D.L., Dai. L., Han J.D., “The ServoHeli-20 Rotorcraft UAV Project,” International Conference on Mechatronics and Machine Vision in Practice, Auckland, New Zealand, pp.92-96, 2008.
  25. 25. Song D.L., Qi J.T., Dai. L., Han J.D. and Liu G., “Modeling a Small-size Unmanned Helicopter Using Optimal Estimation in The Frequency Domain,” International Conference on Mechatronics and Machine Vision in Practice, Auckland, New Zealand, December 24 97102 , 2008.
  26. 26. Song D.L., Qi J.T. and Han J.D., “Model Identification and Active Modeling Control for Small-Size Unmanned Helicopters: Theory and Experiment,” AIAA Guidance Navigation and Control, Toronto, Canada, AIAA-2010-7858, 2010.

Written By

Dalei Song, Juntong Qi, Jianda Han and Guangjun Liu

Submitted: October 22nd, 2010 Published: July 5th, 2011