Abstract

A robust attitude motion synchronization problem is investigated for multiple 3-degrees-of-freedom (3-DOF) helicopters with input disturbances. The communication topology among the helicopters is modeled by a directed graph, and each helicopter can only access the angular position measurements of itself and its neighbors. The desired trajectories are generated online and not accessible to all helicopters. The problem is solved by embedding in each helicopter some finite-time convergent (FTC) estimators and a distributed controller with integral action. The FTC estimators generate the estimates of desired angular acceleration and the derivative of the local neighborhood synchronization errors. The distributed controller stabilizes the tracking errors and attenuates the effects of input disturbances. The conditions under which the tracking error of each helicopter converges asymptotically to zero are identified, and, for the cases with nonzero tracking errors, some inequalities are derived to show the relationship between the ultimate bounds of tracking errors and the design parameters. Simulation and experimental results are presented to demonstrate the performance of the controllers.

1. Introduction

In the field of multivehicle cooperative control, robust consensus tracking under model uncertainties and exogenous disturbances has received increasing attention in recent years, where the output (or state) of each vehicle is required to robustly track a common, desired trajectory. For instance, switching controllers or sliding-mode controllers were designed in works [13] to reject input disturbances; a sliding-mode disturbance observer is combined with a consensus tracking algorithm in recent work [4] to improve the robustness and control accuracy of a multimotor system. Alternative robust control approaches include the ones based on uncertainty and disturbance estimators as in works [57], adaptive control approach [8], and the output-regulation approach [9].

In practice, integral control (IC) is widely used to attenuate disturbances in various (single) vehicle systems. This mainly owes to its structural simplicity and the well-known performance property that IC can asymptotically reject constant input disturbances. Noting these facts, many researchers begin to study IC-based robust control schemes for multivehicle systems (MVSs) as in [1012]. In particular, the recent work [10] shows that PI controllers successfully attenuate constant disturbances in the network of multiple single-integrator dynamics or the network of multiple double-integrator dynamics.

Another practical issue encountered in many control systems is the lack of sensors. As a result, state observers are often needed to generate the estimates of some necessary states. State observers can be roughly classified into two types: model-dependent ones and model-independent ones. As two representative model-dependent observers, Luenberger observer and Kalman filter suffer from the limitation that the estimation accuracy cannot be guaranteed when the system model suffers from severe uncertainties. To deal with this problem, many model-independent observers are proposed. For instance, some higher-order sliding modes (HOSM) observers (differentiators) were designed in work [13, 14] to ensure finite-time convergence even in the presence of input disturbances.

The main objective of this paper is to use integral control to improve the robustness of a distributed control algorithm for consensus tracking without velocity measurements. The effectiveness of the approach is proved by showing that (a) the resulting tracking errors are ultimately bounded for any input disturbances satisfying a simple Lipschitz-constant condition, and (b) zero-error asymptotic tracking is achieved for a constant input disturbance. The key technical differences between this paper and work [10] are summarized as follows:(1)The paper [10] considers the consensus problem without a common reference. In contrast, this paper considers a consensus tracking problem, where a common desired trajectory exists (and is supposed to be accessible only to partial vehicles in the group). In [10], velocity signals were used in the control design for second-order systems. We here assume that neither the velocity of leader nor the velocity of neighboring vehicles is accessible for control design.(2)Concerning the robustness improvement owing to integral control, the discussion in work [10] is restricted to the rather special case with constant input disturbances. This is not the case in this paper. Actually, we will use the concepts of input-to-state stability to study the general cases where the disturbances are nonconstants and are not completely rejected.(3)In work [10], the controller performance is verified by numerical simulation. In this paper, both numerical simulation results and experimental results on three 3-DOF helicopters are presented, to demonstrate the performance improvement owing to the use of integral control.

The experimental platform of “three 3-DOF helicopters” used in this paper is shown in Figure 1. The single laboratory 3-DOF helicopter with a so-called active disturbance system (ADS) is the same as that in [5] and is shown in Figure 2. This experimental apparatus was developed by Quanser Consulting Inc. for the purpose of control education and research [1521].

The rest of this paper is organized as follows. In Section 2, some preliminary knowledge is described and the problem is formulated. The robust distributed consensus tracking controllers are designed in Section 3. Numerical simulation and experimental results are presented in Section 4. Some concluding remarks are drawn in Section 5.

2. Preliminaries

2.1. Notation and Graph Theory

For matrix denotes its inversion and denotes its rank. refers to the identity matrix. For vector and matrix , , , , and , where with being eigenvalues. refers to -dimensional column vector with all elements being . denotes absolute value (modulus) of real number .

The notations related to th 3-DOF helicopter (see Figure 2), , are listed in Notations.

The communication networks of helicopters can be modeled by directed graph , where , , and nonnegative matrix denote the set of nodes, the set of edges, and the weighted adjacency matrix of , respectively. Node () represents th helicopter, and an edge denotes that th helicopter can obtain information from th helicopter; that is, is a neighbor of . Use to denote all the neighbors of , and . A directed path from to is a sequence of ordered edges of form , with distinct nodes , and is said to be reachable from . A node is called the root node and is said to have a directed spanning tree, if all the other nodes are reachable from this node.

The elements of the weighted adjacency matrix satisfy and ( and ) if and only if . The Laplacian matrix of is denoted by , where , if , and . Note that the desired trajectory information is not accessible to all helicopters in this paper. Use constant matrix , which is defined as if th helicopter can access the desired trajectory information and otherwise , to describe this fact by allowing for some .

2.2. Problem Formulation

The elevation and pitch motions of th helicopter can be modeled as follows (see [5] or [21]): whereThe pitch angle is limited to within mechanically.

We made the following assumption on and .

Assumption 1. For each , the first-order derivatives of and with respect to exist and are piecewise continuous and bounded for all .

As shown in recent work [5], by applying the standard normalization and feedback linearization technique, the above nonlinear helicopter model can be reduced to the following 2-dimensional, disturbed double integrators:where , , and the new disturbances .

Under Assumption 1, , also exist and are bounded. Definewhere and are positive scalars.

For each , we define and use to denote the desired attitude trajectory for , which may be time-varying but are second-order differentiable with respect to . Then, the objective of this paper is to design for (3) to achieve robust attitude synchronization; that is, for each as .

3. Design of Distributed Controllers without Velocity Measurements

3.1. Controller Design

Noting that and for (3) can be designed in the same way, we thus introduce new unified variable and defineIn [8, 22], error is called local neighborhood synchronization error (LNSE).

For the considered helicopter, only angular position sensors (encoders) are equipped. Besides, we assume that the desired velocity is not accessible to any helicopter. To deal with this problem, we consider the following distributed controllers:where , , and are positive control gains, , , and denote the estimates of and , respectively. We construct the following systems to generate and . Specifically, the third-order FTC estimatorsare for , and the following second-order FTC estimators are for :where , , and are the estimates of , , and , respectively; and are consisting of a locally bounded Lebesgue-measurable noise with unknown features and an unknown base signal and whose second and first derivative have an known Lipschitz constant and , respectively; and , are the design parameters of the FTC estimators.

This estimation scheme is illustrated in Figure 3.

We have the following result for the above two FTC estimators [14, Theorem  1].

Lemma 2. Let and in (7) and (8) be recursively chosen as in Theorem of [13, 14]. Then, the estimation errors achieve zero in the absence of input noises after finite time of transient process. The convergence of the estimation errors is uniform in the sense that the convergence time is uniformly bounded by finite time which is a locally bounded function of initial estimation errors.

Remark 3. According to Lemma 2, in the absence of input noises, after finite time (), the following equalities hold:

Remark 4. Note that the neighbors’ control signals () are used in (6) and then a possible implementation loop issue arises in practical applications. As discussed in work [5], if the sample frequency is high enough, this problem can be resolved with using the neighbors’ control signals obtained during the previous sampling period, that is, , where denotes the fixed sampling step and in the following numerical simulations and experiments.

3.2. Analysis of the Closed-Loop Stability

Before the closed-loop stability analysis, we need to establish and analyse the relationship between tracking errors and LNSEs . DefineFrom (5c), ; then from (5d), for each . Then, with and as defined in Section 2.1, the following relationship equation between the tracking errors and LNSEs is derived:where is determined by the communication topology which satisfies the following condition in this paper.

Assumption 5. The desired trajectory information , has directed paths to all helicopters.

According to [23, Lemma  1.6], all eigenvalues of have positive real parts under Assumption 5. Then, the following lemma is obtained directly.

Lemma 6 (see [23, Lemma  1.6]). Under Assumption 5, matrix is full rank, that is, .

Thus, the relationship equation (11) implies that, under Assumption 5, the objective and as , , can be achieved by driving and to zero as , respectively.

Applying (6) to (3) gives

With (5d) and (12), we can further getwhere and is defined as

DefineThen, the following result is clear based on Lemma 2.

Lemma 7. There exists such that, in the absence of input noises,

Therefore, after finite time , (13) is the same as follows:DefineThen, (20) can be written asFurthermore, based on Lemma 7, by applying Routh-Hurwitz stability criterion and Lemma A.1 in Appendix to forced systems (13) with as inputs, the following theorem is obtained.

Theorem 8. Consider systems (13) under Assumptions 1 and 5 and the following parameter condition:Then, for each , in the absence of noise,(1)the system trajectories , , , and are globally bounded;(2), , , and as if as ;(3)for the general cases with nonvanishing , there exists such that, under any initial condition,whereand , are the solutions of the Lyapunov equationswith defined by (21).

Proof. Using (14), (16), and (17), the relationships between and can be written in a compact form asThen, with (4), it follows thatMoreover, for , based on Lemma 6, under Assumption 5, as if and only if as . Thus, under Assumptions 1 and 5 and parameter conditions (24), it is straightforward to derive the points ()-() with these results and Lemmas 6 and 7.
The detailed proof of point () is the same as [5, Theorem ] and hence omitted.

Remark 9. As shown in Theorem 8, by FTC differentiators (7) and (8) in the absence of input noises, for the case without angular velocity measurements, if condition is satisfied, the distributed controllers (6) are capable of achieving zero-error attitude-trajectory tracking for each helicopter in the group of helicopters. Thus, the distributed consensus attitude-trajectory zero-error tracking for the group of helicopters is achieved under this condition (). If this condition is not satisfied, the ultimate bounds of attitude-trajectory tracking errors resulting from controllers (6) are globally bounded for and ultimately bounded with ultimate bounds depending on parameters which can make the term small enough. Future efforts will be devoted to the optimization problem of the three parameters such that the term is minimized.

Remark 10. If delay in neighbors’ control inputs is systematically considered, the obtained error equations (13) are disturbed neutral delay systems as in [24]. For these kinds of delayed systems, [25, Lemma  1] can be used to examine the stability of the associated nominal delayed systems, where three easily testable conditions are included. We also note that these conditions ensure that the considered systems are delay-independent stable. However, the delay-independent stability is not the general case. Despite this fact, since the case with is exponentially stable, we readily derive that the system with an enough small time-delay is also exponentially stable, which ensure that the correspondingly disturbed equations are input-to-state stable provided that the delay is sufficiently small. In our future work, we will investigate these problems and systematically assess the effect of the introduced delay on the system stability as well as the convergence speed of tracking errors (as in the work [26]).

4. Numerical Simulation and Experimental Results

In this section, to demonstrate the effectiveness of the proposed scheme for robust motion synchronization, the numerical simulation results and experimental results of the attitude-trajectory consensus tracking of three 3-DOF helicopters which are labeled as H1 to H3 are presented. Note that we focus on the analysis of the results of elevation-trajectory consensus tracking hereafter. The results of pitch channel are similar to that of elevation channel and hence omitted.

The nominal parameters of the four helicopters (including virtual helicopter H0) are the same, which are presented in Table 1. Their initial states are specified as follows:

All the results are obtained on the directed communication graph shown in Figure 4. In this graph, helicopters and have access to desired trajectories which are the responses of H0 to desired signal , and Assumption 5 holds with and all other entries of and being . Then,

Without loss of generality, the following nonstep desired trajectories provided to H0 are adopted for numerical simulations and experiments:where all the angles are given in degrees. For the parameters of control signal which is determined by (6) with , , , and , choosewhere H0 is a virtual helicopter without the acting of LUDs; hence the integral term in the control law (6) is not adopted by H0 (i.e., ) for simplicity.

The design parameters chosen for the FTC estimators (7) and (8) are shown in Table 2. The initial state values of the differentiators in (7) and (8) are set to zero; that is, and

In the following, for comparison purposes, six cases (Cases ) are considered. In each case, are determined by (6). In order to theoretically verify the results of Theorem 8, Cases are designed for the numerical simulations. Specifically, to demonstrate that both tracking and synchronization performance without IC are not acceptable, IC is not adopted in Case ; points ()-() of Theorem 8 are verified by Case , where the helicopters with IC are subject to different constant disturbances; the helicopters with IC are subject to different time-varying disturbances in Case to show the validity of points () and () of Theorem 8. The results of Theorem 8 are experimentally verified by Cases on the experiment platform shown in Figure 1. IC is adopted in Cases - while, in Case , it is not. The ADSs of helicopters 2 and 3 are activated in Case while in Cases - are not. Note that, since the disturbances are unknown in the experiments, the control parameters of in experiments are not the same as those in numerical simulations, where the disturbances are specified and used to choose the control parameters.

The detailed differences among Cases are as follows:

() For Case , the numerical example with constant disturbance (, ) is simulated. The parameters of are chosen as , which corresponds to the situation without IC.

() Case is the same as Case , except for . That is, IC is applied for the numerical example with constant disturbance.

() Case is the same as Case , except for time-varying that are given as follows:Note the three points: (a) as shown in (35), the three helicopters are subject to different disturbances; (b) initial disturbances ; (c) to verify the results of Theorem 8, the disturbances used in Case satisfy conditions for each , while those used in Case do not.

() For Case , without IC are applied to the experiment platform shown in Figure 1. The parameters of are chosen as .

() Case is the same as Case , except for . That is, IC is applied to the experiment platform.

() For Case , in the experiments, the ADSs of helicopters 2 and 3 are activated and the same control laws as for Case are applied.

It is worthwhile noting that ADS has a dramatic effect on the motion of helicopter body, whether it is static or moving. This point is easily seen from the mechanical structure of helicopter as shown in Figure 2 and is also verified by the following experiments.

4.1. Case 1: Numerical Simulation without Using IC

The numerical simulation results for this case are presented in Figure 5. As shown in subfigures (a) and (b) therein, because IC is not adopted in control design to attenuate constant disturbances, neither tracking error nor synchronization error converges to a small neighborhood of zero. More specifically,(1)for the elevation axis, the magnitude of tracking error (of each helicopter) is greater than when and is steady as increases. The synchronization error between any pair of the three helicopters is of a magnitude greater than when and is steady as increases. The reason is that the constant disturbance of the helicopter is different from each other and is not attenuated by IC.(2)as shown in subfigures (c) and (d) of Figure 5, the voltage applied to either front motor or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), has a magnitude smaller than  V for all (smaller than  V for all ).(3)as shown in subfigures (e) and (f), the attained accuracies of the FTC estimators are and that corresponds to Lemma 2 for FTC estimators (7) and (8) with finite time (in this case, and for (7) and (8), resp.).

4.2. Case 2: Numerical Simulation with Using IC for Constant Disturbances

In this case, IC is adopted in control design to attenuate constant disturbances. The corresponding simulation results are presented in Figure 6, demonstrating that both tracking and synchronization performance are dramatically improved, compared with that achieved in Case . A detailed analysis of the results is given as follows:(1)For the elevation axis, when , the magnitude of tracking error (of each helicopter) is smaller than and the synchronization error between any two helicopters has a magnitude smaller than . Points () and () of Theorem 8 are verified by these results.(2)The voltage applied to either front motor or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), is of a magnitude less than  V for all .(3)The result of FTC estimator (7) is the same as that in Case and hence omitted. The attained accuracies of FTC estimator (8) are and the finite convergence time satisfies .

4.3. Case 3: Numerical Simulation with Using IC for Time-Varying Disturbances

In this case, IC is adopted in control design to attenuate time-varying disturbances. The corresponding simulation results are presented in Figure 7. Although subject to the time-varying disturbances shown in (35), both tracking and synchronization performance are improved compared with that achieved in Case , where the disturbances are constants. However, compared with that achieved in Case , where the disturbances are constants and also attenuated by IC, both tracking and synchronization performance are degraded observably. A detailed analysis of the results is given as follows:

() For the elevation axis, when , the magnitude of tracking error (of each helicopter) is smaller than ; that is,and the synchronization error between any two helicopters has a magnitude smaller than .

With parameters and (21), the following matrix can be obtained by resolving the Lyapunov equations (28):With (4) and (35), we have ; then, along with (32) and (37), we get in (27) that satisfies and the ultimate bound of the tracking errors in (26) satisfies . Then, points () and () of Theorem 8 are verified by (36).

() The voltage applied to either front motor or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), is of a magnitude less than  V for all .

() The result of the FTC estimators is the same as that in Case and hence omitted here.

4.4. Case 4: Experiment without Using IC

The experimental results for this case are presented in Figure 8. It is seen that the tracking error (of each helicopter) does not converge to a small neighborhood of zero and the synchronization performance is better than tracking performance. A detailed analysis of the results is given as follows:(1)For the elevation axis, the magnitude of tracking error (of each helicopter) is greater than for all ; the synchronization error (between any two helicopters) is , , and for all ; the synchronization error between H2 and H3 has a smaller magnitude than that between H1 and H2 (or H1 and H3), which is an immediate consequence of the fact that H2 and H3 can obtain information from each other and both of them are equipped with ADSs, whereas H1 is not.(2)The voltage applied to either the front or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), has a magnitude smaller than  V for all .

4.5. Case 5: Experiment with Using IC

The experimental results for this case are presented in Figure 9. It is observed that, compared with Case , both tracking and synchronization performance are improved because the disturbances are attenuated by IC in this case. More specifically,(1)for the elevation axis, the magnitude of tracking error (of each helicopter) is smaller than for all (, , and ); the synchronization error (between any two helicopters) is , , and for all ;(2)the voltage applied to either the front or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), has a magnitude smaller than  V for all .

4.6. Case 6: With Activated ADS and the Same Controller as for Case

The ADSs equipped on helicopters 2 and 3 are activated in this case, whose dynamic positions along the arms of helicopters are shown in Figure 10, with the same initial position m (m is the farthest possible position from propellers). This is different from the previous Case and Case where they are both fixed at m. The experimental results for this case are presented in Figure 11. It is seen that motion of ADSs does not lead to dramatic performance degradation in the steady state of tracking and synchronization errors, compared with the performance achieved in Case . A detailed comparative analysis of the results is given as follows:(1)For the elevation axis, the ultimate bound of either tracking error or synchronization error is nearly the same with that achieved in Case ; the transient performance is degraded a little (this degradation is relatively obvious for helicopter 2 because the effect of the moving ADS on helicopter 4 is delivered to helicopter 2 which has no access to the desired trajectory as shown in Figure 4).(2)The voltage applied to either the front or back motor, with a value equal to the maximum acceptable voltage  V at the beginning (), has a magnitude smaller than  V for all .

To demonstrate the different steady-state tracking performance clearly, the maximum magnitudes of tracking errors that appeared in each case during time interval are shown in Table 3 (Cases ) and Table 4 (Cases ). As shown in these tables, the tracking accuracy is improved in numerical simulations and experiments by IC.

5. Concluding Remarks

The robust distributed consensus tracking controllers with integral action for multiple 3-DOF experimental helicopters without velocity measurements have been studied under the condition that only the desired angular position measurements are accessible to a small subset of the helicopters. Motivated by the effectiveness of the tracking controller with integral action in disturbances rejection for the single vehicle, the distributed controllers have been proposed by combining FTC estimators with distributed integral controllers. With using the FTC estimators, great accuracy and finite-time convergence have been achieved for the estimation of the lacking information. Meanwhile, the distributed controllers with integral action stabilized the tracking errors and rejected the input disturbances. Through analysing the closed-loop stability, the conditions ensuring zero-error tracking and the ultimate bound of errors for the general cases with nonzero error have been derived. It has been verified through the results of numerical simulations and experiments on platform of “three 3-DOF helicopters" that the tracking and synchronization accuracy have been improved by the proposed controllers with proper parameters.

Future work will focus on the design and experimental verification of bounded distributed controller for the platform under the directed communication graph with time-delay. The control parameters in (6) also need to be optimized to get the minimized ultimate bound of tracking errors.

Appendix

A useful lemma used to show boundedness and ultimate bound of the solutions of some disturbed state equations is given as follows.

Lemma A.1 (see [5, Lemma  2]). Consider state solution of the linear time-invariant equation:where is the state, is the continuously differentiable input, and matrices and . If is Hurwitz, then(1)system (A.1) is globally input-to-state stable (GISS); that is, if is bounded for all , then with any initial state is also bounded for all ;(2)there exist class function and time (dependent on and ) such that with any satisfieswhere , , and are the maximum and minimum eigenvalues of the symmetric positive-definite matrix which is the solution of the Lyapunov equation:

Notations

, The elevation angle and pitch angle, respectively (rad)
, The elevation angle rate and pitch angle rate, respectively (rad/s)
, :The moments of inertia about elevation axis and pitch axis, respectively (kg·m2)
:The force constant of motor-propeller combination (N/V)
:The distance from travel axis and elevation axis to the center of helicopter body (m)
:The distance from pitch axis to either motor (m)
:The effective mass of helicopter body (kg)
The gravitational acceleration constant, (9.81 m/s2)
, :The lumped uncertainties and disturbances (LUDs) acting on elevation and pitch channels, respectively (N·m)
, :The voltages applied to the front motor and back motor, respectively; the voltage limit for the motors is 24 V (i.e., ,) (V)
, :The sum and the difference of and , respectively (V).

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grants nos. 61304016 and 61305132, the Fundamental Research Funds for the Central Universities under Grant no. ZYGX2016KYQD102. The authors would like to thank Professor Hugh. T. Liu at UTIAS for generously providing the 3-DOF helicopter platform at FSC lab to verify the controllers developed in this paper.