Abstract

In order to investigate the attitude containment control problem for a microsatellite cluster, an event-triggered adaptive sliding mode attitude containment control algorithm is proposed for the satellite cluster flight system under directed topology, so that attitude of followers asymptotically converges to the convex hull formed by the leaders’ orientations. At first, the event-triggered control strategy is introduced into the attitude containment control problem for the microsatellite cluster. The triggering condition consisting of state-dependent and time-dependent function is designed to adjust control period and avoid the Zeno behaviour. When the function value meets the triggering condition, the event is triggered, state information is sampled, control law is computed, and actuators are updated, while the control action performed in nontriggering time is the same as the previous triggering instant. Then, in the presence of model uncertainties and external disturbances, an event-triggered adaptive sliding mode attitude containment control algorithm is presented under directed topology, and sufficient and necessary conditions for the followers to enter into the target area formed by the leaders are given. Furthermore, cell partitions from graph theory are employed to investigate the influence of information topology on steady states of followers, which provides theoretical basis for orientation design of cluster satellites. Finally, simulation results show that the proposed control strategy could reduce control execution frequency, as well as ensure the similar control performance with the time-triggered one, and followers belonging to the same cell have the same steady states.

1. Introduction

Satellite cluster has received growing attention in the recent years [1, 2] for its advantages of greater flexibility, faster response, higher reliability, lower cost, and better reconfigurability [3]. Contrary to satellite formation flight, cluster flight does not impose strict limits on the geometry of the cluster [4] and is hence more suitable for implementation by multiple microsatellite systems. Satellite cluster has been deployed by many institutes, such as ANTS [5], Breakthrough Starshot project [6], and so on.

Several technical barriers need to be broken down to pave the way for cluster satellites to come into being. Coordinated attitude control of cluster satellites has been identified as one of the enabling technologies. Although there has been lots of results on coordinated control problem for multiple satellites systems, we note that most of the existing research studies consider leaderless [7, 8] or one leader case [9, 10], where there exists no group objective or only a single group objective.

However, the presence of multiple leaders is more attractive for the satellite cluster system, owing to the fact that such strategies provide attractive solutions to cluster problems, both in terms of complexity and computational load. As a kind of extended consensus problem, the case with multiple leaders is what we call containment control [11]. The objective of containment control is to guarantee that all the followers asymptotically converge into the convex hull formed by the leaders through information interaction and coordinated control. The containment control problem received significant research interest due to its various applications, such as mobile robots [12], unmanned aerial vehicles (UAVs) [13], underwater vehicles [14], and satellite formation systems [15, 16].

By now, there has been lots of research studies on distributed containment control, and research studies could be divided into several types from different perspectives. According to the dynamics, the research studies include first-order systems [12, 17, 18], second-order systems [1921], linear systems [12, 17, 18], nonlinear systems [15, 2224], homogeneous and heterogeneous multiagent systems [25], etc. In view of information topology, fixed topology [15, 19, 22, 23], and switching topology [18], undirected graph [20] and directed graph [17, 23] are, respectively, considered. In many research studies, containment control is combined with other novel control strategies, such as finite-time control [14, 15, 22], adaptive control [14, 16], neural network [16], event-triggered control [26], and so on. For the design of leaders, there exists the case of stationary leaders [20, 24], dynamic leaders (constant velocity or time-varying velocity) [17], leaders formation [25, 27], and so on. In addition, other problems such as time delay [11, 21], model uncertainties [15, 16, 23], external disturbances [15, 16], collision avoidance [27], and unmeasured relative velocity [28] are also discussed.

Recently, distributed attitude containment control strategies have gained increased attention in satellite coordinated control community. In [22], the distributed finite-time attitude containment control problem for multiple rigid bodies was addressed. For multiple stationary leaders, a model-independent control law was proposed to guarantee the attitudes of followers converge to the stationary convex hull formed by leaders in finite time by using both the one-hop and two-hop neighbours’ information. Then, for dynamic leaders, a distributed sliding mode estimator and a nonsingular sliding surface were given to guarantee the attitudes and angular velocities of followers converge, respectively, to the dynamic convex hull formed by those of the leaders in finite time. Under undirected fixed connected graph, Weng et al. [29] investigated distributed robust finite-time attitude containment control for multiple rigid bodies with uncertainties including parametric uncertainties, external disturbances, and actuator failures. In [30], a distributed control strategy combined with finite-time command filtered backstepping (FTCFB) and an adaptive technique was established to solve the attitude containment control problem of spacecraft formation flying (SFF) with unknown external disturbances. The proposed novel distributed adaptive FTCFB approach could guarantee the containment errors of attitudes between leader spacecraft and follower spacecraft reach the desired neighbourhood in finite time under undirected topology. However, the information topology of cluster satellites may be directional in actual space missions. Because only a fraction of satellites was equipped with necessary sensors or communication equipment to measure relative state in cluster system, obviously, the directional information topology is more general. Ma et al. [31] studied the distributed finite-time attitude containment control problem for multiple rigid body systems with multiple stationary and dynamic leaders under directed graph. Based on sliding mode observer, adaptive attitude control algorithms were given, and the necessary and sufficient conditions were achieved which rendered all the followers converge to the convex hull spanned by the static and dynamic leaders in finite time. In [24], an attitude containment control algorithm was proposed in the case of undirected angle information topology and directed angular velocity information topology, and the case of unavailable relative angle velocity was also investigated.

Continuous or periodic sampled control scheme is usually applied to the aforementioned attitude containment control problem of satellite formation, whose results belong to time-triggered control. The state information of cluster satellites is usually sampled with a fixed and high sampling frequency, and the actuators are updated at each sampling instant, which increase the pressure of the whole network communication and lead to the wear of actuators and unnecessary energy consumption, thus seriously shortening the in-orbit operation life of cluster satellites. Moreover, in time-triggered control schemes, control action updates periodically even when the system has reached the desired state with satisfactory accuracy. Computation and communication pressure will be greatly increased, while resource and network bandwidth of the microsatellite cluster are extremely limited.

Efforts to overcome these problems have led to the proposal of event-triggered control strategy. Information interaction and controller updates are not determined by time, but by the triggering condition (event). Event-triggered mechanism consists of two types: state-dependent events [32] and time-dependent events [33]. In the event-triggered control strategy, the control tasks, consisting of sampling state information of satellites, computing control law, and updating actuators, are executed when the well-designed triggering condition is satisfied. Thus, communication and computation resources are utilized only when “needed” to preserve desired control performance [34]. It makes event-triggered control favourable, especially for satellite cluster missions with limited bandwidth and resources.

So far, event-triggered control has been investigated in the multi-rigid body system model with nonlinear characteristics. In [35], an event-triggered distributed adaptive controller was proposed to study the leader-follower consensus problem for a directed network of Euler–Lagrange agents. For the attitude control problem of spacecraft, Wu et al. [36] investigated the problem of spacecraft attitude stabilization control system with limited communication and external disturbances based on an event-triggered control scheme. Sun et al. [37] introduced an event-triggered control (ETC) strategy for the spacecraft attitude stabilization problem from the view of cyber-physical systems (CPSs); a new quaternion-based nonlinear control algorithm was proposed to ensure attitude dynamics systems’ exponential stability, and parameter selection of event function and controllers was discussed in this paper. There are also research studies combining event-triggered schemes with containment control. In [38], distributed event-triggered cooperative attitude control of multiple rigid bodies with leader-follower architecture was investigated; under the designed controllers with the event-triggered strategies, the orientations of followers converge to the convex hull formed by the desired leaders’ orientations with zero angular velocities. Xu et al. [26] studied the distributed event-triggered adaptive containment control problem for multiple Euler–Lagrange systems with stationary/dynamic leaders over directed communication networks.

Although various novel control strategies have been investigated for the attitude containment control problem of a satellite cluster, which enables cluster members to converge to the convex hull formed by leaders with a faster convergence rate, little attention has been paid to the relationship between system performance and information topology design of the satellite cluster system. Note that the interaction between satellites need not be bidirectional in practice due to communication bandwidth or sensor capability. Constrained by intersatellite distance and performance of sensors, only parts of followers can receive information from leaders directly. To the best of the authors’ knowledge, the event-triggered attitude containment control for the microsatellite cluster system under directed topology is worth studying and awaits a breakthrough. However, the existence of system uncertainties and unavoidable external disturbances of cluster system results in an unsatisfactory performance [39, 40]. Thus, the sliding mode control (SMC) strategy, which is robust to external disturbances and model uncertainties, is employed. The adaptive control method is also combined to realize the online estimation of uncertain parameters in real time, and it would not destroy the robustness properties of SMC [41].

In this paper, the attitude containment control problem and information topology structure design for the microsatellite cluster are investigated. First, the triggering condition consisting of the relative state error is given to adjust controller update period. If and only if the triggering condition is satisfied, state information is sampled, control law is computed, and actuators are updated. Then, an event-triggered adaptive sliding mode attitude containment control algorithm is proposed in the pressure of inertia uncertainties and external disturbances, which makes attitude of cluster members to enter asymptotically into the convex hull formed by leaders’ orientations. Furthermore, cell partitions in the view of graph theory are employed to investigate the influence of information topology on orientation of followers, which provides theoretical basis for information topology design of satellite cluster missions. Finally, simulation results show that the proposed event-triggered adaptive sliding mode attitude containment controller could drive the followers to enter into the convex hull formed by the leaders’ orientations in the presence of inertia uncertainties and external disturbances, and followers belonging to the same cell have the same orientation. Compared with the existing results, the proposed results in this paper have the following advantages.(1)In the framework of the Euler–Lagrange system, this paper presents an event-triggered adaptive sliding mode attitude containment control algorithm for the microsatellite cluster system under directed topology, so that the followers asymptotically converge to the target area formed by the leaders’ orientation in the presence of inertia uncertainties and external disturbances. The controller is updated only when the triggering condition is satisfied, and the state information of the triggering instant is utilized at the nontriggering time. Therefore, the control input is a piecewise function and the controller update for each satellite is asynchronous. Compared with the time-triggered method, the event-triggered adaptive sliding mode attitude containment control algorithm not only ensures the similar control performance of cluster satellites but also effectively reduces information transmission and actuator update frequency of the satellite cluster system. Event-triggered attitude containment control is superior in microsatellite cluster missions with limited resource and lower precision.(2)The triggering condition of time-varying threshold is given in this paper. Most researches only study state-dependent or time-dependent triggering condition. However, the triggering condition in this paper is the combination of state-dependent and time-dependent. The time-dependent function is introduced to avoid the Zeno behaviour, i.e., the controller does not update infinitely in finite time, nor does it update in a periodical manner. When the triggering condition is satisfied, state information is sampled, control law is computed, and actuator is updated, which can effectively reduce the computation and actuator update frequency while ensuring the control performance of cluster system.(3)From the perspective of graph theory, the control algorithm is fully distributed in the sense that each satellite can select their control gains according to only local information. Then, the influence of information topology design on orientation of the microsatellite cluster is analysed. It is shown that in an ideal environment, the stable state of each follower is a convex combination of all leaders’ states it can access. Two kinds of cell partition of cluster information topology are given, and it is proved that satellites belonging to the same cell partition have the same stable state. It provides a theoretical basis for information topology design of microsatellite cluster missions with performance requirements such as full coverage of the target area.

The remainder of this paper is organized as follows. Mission scenarios, relative attitude dynamics of satellite cluster, basic knowledge of graph theory, and attitude containment control problem are briefly given in Section 2. In Section 3, the event-triggered adaptive sliding mode attitude containment control algorithm is proposed for the satellite cluster in the presence of inertia uncertainties and external disturbances. Then, sufficient and necessary conditions for asymptotical convergence of the cluster system and the absence of the Zeno behaviour are derived. Orientation of cluster satellites based on information topology design is given in Section 4, providing theoretical basis for information topology design of the microsatellite cluster. Numerical simulations to verify the effectiveness of the proposed control algorithm and information topology structure design are completed in Section 5, and concluding comments are given Section 6.

2. Preliminaries and Problem Formulation

In this section, some problem descriptions about mission scenarios are introduced, and then preliminary knowledge about containment control strategies is given.

2.1. Mission Scenarios

Traditional single leader consensus is a group of satellites aiming to achieve an agreement through information interaction and coordinated control, where leader’s trajectory will be tracked by followers of the cluster system. In contrast, containment control can drive the followers to enter into the target area formed by the multiple leaders. Containment control is more robust in the case of leader failure and more practical in microsatellite cluster missions with limited resource and lower precision. Two typical microsatellite cluster flying scenarios which are closely related to attitude containment control are firstly given and analysed as follows.

2.1.1. Earth Observation or Deep Space Exploration of Microsatellite Cluster

Microsatellite cluster has been employed in many observation missions, such as the Earth observation or deep space exploration. It is necessary for satellites to obtain and maintain a certain relative attitude [32]. In the observation missions, such as expanding observation view or searching observation target, members are not required to converge to the same orientation, but to enter into a target area formed by the leaders’ orientations, which is called attitude containment control. One of the basic objectives is that a subset of the satellite set (leaders) stabilizes to a specific relative attitude, and the orientation of the rest members (followers) enters into and remains within the specific attitude, determined by a convex hull formed by the leaders’ orientation. Meanwhile, each satellite is only allowed to communicate (attitude or angular velocity) with a specific member of the set, and these constraints limit the information interaction between satellites.

2.1.2. Coordinated Attitude Control for Fractionated Spacecraft

Fractionated spacecraft distributes the functional capabilities of a monolithic spacecraft into multiple free-flying, wirelessly linking modules (service modules and different payloads) [42]. One of the main challenges of this architecture is cluster flight, keeping the various modules within bounded configurations. The fractionated spacecraft generally does not require precise relative orbit and attitude control, as long as the relative distance is within the range of communication and the relative attitude control enables power transmission links [43]. In this case, some (virtual) module spacecraft form an orientation area for the cluster members’ attitude, which makes modules to maintain communication and power transmission.

There is no precise requirement for the final attitude of cluster members in the aforementioned attitude control missions of the microsatellite cluster. The attitude of followers only needs to reach an area instead of the consensus state. It is necessary to form the target area using (virtual) leaders’ orientation and then control the followers to enter into the target area through intersatellite information interaction and properly designed coordinated attitude control protocol.

Constrained by unidirectional measurement characteristics of sensors or GPS-like radio frequency communication, as well as relative distance limitation between satellites, information topology of microsatellite cluster is unidirectional, asymmetric, and sparse. And the information topology structure generally remains fixed in a short time if there does not exist large disturbance. In the following, we will study the attitude containment control problem of the microsatellite cluster under fixed directed topology.

2.2. Relative Dynamics Model of Satellite Cluster

In this paper, the modified Rodriguez parameters (MRPs) are used to describe attitude motion of the satellite cluster. MRP is a kind of attitude description method without redundancy and singularity. Attitude vector of MRPs is expressed by Euler axis and Euler angle as follows:

The kinematic equation of follower iswhere and , respectively, are MRPs and attitude angular velocity of the fixed-body coordinate system of satellite i relative to the inertia coordinate system, , is the 3×3 identity matrix, and is the skew-symmetric matrix of and is expressed aswhere represents the th component of vector .

Attitude dynamics equation of satellite i is expressed aswhere is the positive inertia matrix, represents the control torques, and represents the external disturbance torques.

Attitude kinematic equation using MRPs as the attitude parameter can be converted into an unified Euler–Lagrange form by appropriate transformation [44]. The Euler–Lagrange equation can effectively introduce the linear relative attitude expression between different satellites, which plays a beneficial role in coordinated attitude control design of satellite cluster. Taking the derivative of (2) and substituting (3) into (2), the Euler–Lagrange form of attitude dynamics equation for satellite i could be obtained as follows:where , is symmetric positive-definite inertia matrix of satellite i, is Coriolis and centrifugal torques matrix, , is nominal inertia of follower i, is uncertainties part, is uncertainties causing by model uncertainties, is control torques, and is disturbance torques. According to Reference [45], the Euler–Lagrange equation has the following properties: (a) is a symmetric positive-definite matrix; (b) here we assume that there exist positive satisfying , , ; (c) for any , is a skew-symmetric matrix, and there exist , ; (d) we assume that the left-hand side of (5) can be parameterized in linear expression as , for any , where is the regression function matrix and is the unknown constant parameter vector.

Assumption 1. All leaders’ states and state derivatives are bounded, that is, there exist constants and such that ,. In this paper, we add an additional assumption that .

Assumption 2. The external disturbance torque is bounded, that is, there exists a positive constant such that , where is bounded and unknown.

2.3. Graph Theory

Graph theory is introduced as an useful mathematical tool to describe intersatellite information interaction, providing theoretical basis for control performance analysis of the cluster system.

Suppose there are N + M satellites in the cluster, consisting of N followers (denoted by 1, 2, …, N) and M leaders (denoted by N + 1, …, N + M). We assume that the (virtual) satellites belong to either one of the two subsets, namely, the subset of followers or the subset of leaders .

Information interaction topology of the cluster system can be modelled by a digraph with node set denoting satellite with dynamics or kinematic characteristics and edge set denoting intersatellite information interaction. Each edge means satellite j can access state information from satellite i, and satellite i is called a neighbour of satellite j. The neighbour set of satellite i is a collection of all the satellites from which satellite i can access state information. If , bidirectional edge indicates that satellites i and j can access state information from each other. A directed path from node i to node j is a sequence of edges of the form (i, i1), (i1, i2), …, (in, j), in a directed graph. A directed tree is a directed graph, where every node has exactly one parent except for one node, called the root, and the root has directed paths to every other node. A directed spanning tree of a directed graph is a direct tree that contains all nodes of the directed graph. A directed graph has or contains a directed spanning tree, if there exists a directed spanning tree as a subset of the directed graph.

Two matrices are frequently used to represent interaction topology: one is adjacency matrix with , if . In this paper, we assume that self-edges are not allowed, i.e., . And the other is Laplacian matrix with , .

The system Laplacian matrix L can also be written as block matrix:where is the N × N submatrix composed by rows and columns in which the followers of L are located, and it represents information interaction relationship among followers. is the submatrix composed by rows and columns in which the followers and leaders of L are, respectively, located, and it represents the information flow from leaders to followers. denotes followers and information interaction between followers.

Several graph theory tools are given to provide theoretical basis for information topology design of the cluster system.

Definition 1 (see [46]). For directed graph , a k partition of is composed of k cells with , and . Cj is a neighbour of Ci if there exist and such that is a neighbour of .

Definition 2 (see [47]). Suppose that is a k partition of ; if any two distinct vertices in have the same number of neighbours in Cj for all , then is called a partition. In particular, denotes the cell partition containing as m cells of .

Definition 3 (see [46]). Suppose is a k partition of ; if for each vertex in , it has the same number of neighbours in all neighbour cells of , then is called a partition.

2.4. Containment Control Model

The following definitions, assumptions, and lemmas related to containment control strategy are needed to derive the main results of this paper.

2.4.1. Convex Hull

There may exist multiple leaders in microsatellite cluster missions, and all followers are required to enter into the target area formed by the leaders’ state instead of reaching a consensus state. Several vertices on the boundary of the target area are selected so that the target area can be approximately replaced by convex polygons which are formed by these vertices [48]. Generally, the more the vertices are selected, the higher the approximate accuracy is.

Suppose that the target area (moving or stationary) can be approximated by a convex hull formed by M vertices.

The definition of convex hull is given as follows.

Definition 4 (see [23]). Let be a set in a real vector space . The set is convex if, for any x and y in , the point for any . The convex hull for a set of points in is the minimal convex set containing all points in . We use to denote the convex hull of X, that is, .
Vertex information of the target area could be provided by the Earth station or could be autonomously generated by cluster members which have strong sense, communication, and information processing capability.
Once the convex hull which approximates the target area is selected, vertex information of the convex hull can be seen as (virtual) leaders of the cluster system, while cluster members are regarded as followers.

2.4.2. Distributed Attitude Containment Control Formulation

Followers need to generate control decisions based on absolute state of itself and relative state information of its neighbours, which makes attitude of followers to converge to the convex hull formed by the leaders’ orientation.

Containment control protocol for followers could be written in the following form:where , , , and , respectively, are the attitudes and angular velocities of followers and leaders.

We assume that there is no information interaction between leaders. The trajectories of leaders are not affected by other members, while followers need to generate control instructions with neighbours’ information [49]. To drive all the states of followers to enter into the convex hull formed by the leaders’ orientation, it must be ensured that each follower can receive information from leaders directly or indirectly. Otherwise, there will exist followers whose motion is not affected by any leader, nor will converge to the convex hull formed by the leaders’ orientation. Thus, the information topology of the cluster system needs to meet the following assumption.

Assumption 3 (see [20]). Suppose that for each follower i, there exists at least one leader that has a path to the follower i.

Lemma 1 (see [20]). If Assumption 3 holds, then is invertible; each entry of is nonnegative and each row sum of is equal to one.

Lemma 2 (see [23]). Let be the Laplacian matrix associated with a directed graph . Then, has a single zero eigenvalue and all other eigenvalues have positive real parts if and only if has a directed spanning tree.

Lemma 3 (see [50]). Consider the system , where is continuously differentiable and globally uniform Lipschitz in for all t. If the unforced system has a globally exponentially stable equilibrium point at the origin , then the system is input-to-state stable.

According to Definition 4 and Lemma 1, the desired state is convex weighted average of the leaders' attitude and angular velocity, which can be written aswhere and , respectively, are the desired attitude and desired angular velocity.

Our aim in this paper is to propose appropriate distributed attitude control algorithm for the followers (i.e., those indexed from 1 to N), so that in an asymptotic manner, they can travel into the convex hull formed by the leaders (i.e., those satellites indexed from N + 1 to N + M). We will also analyse under what conditions the containment behaviours can be guaranteed and perform rigorous convergence with less sampled information and control action, that is, , , . On this basis, the influence of information topology design on orientation of followers is analysed, which provides a theoretical reference for orientation design of cluster members.

3. Distributed Attitude Containment Control for Microsatellite Cluster

For attitude containment control problem of leader-follower satellite cluster, an event-triggered adaptive sliding mode control algorithm is proposed in the presence of inertia uncertainties and external disturbances. The triggering condition based on relative attitude error is set for each follower, that is, the entire system is triggered asynchronously and triggering time of each follower is different. State information is sampled, control law is computed, and actuators are updated if and only if triggering conditions are met. At nontriggering time, controller of followers uses the state information of triggering instant.

The event-triggered adaptive sliding mode controller and triggering condition are designed to make attitude of followers asymptotically enter into the convex hull formed by the leaders and attitude angular velocity converge to 0; meanwhile, the update frequency of control tasks is reduced.

3.1. Event-Triggered Formulation for Attitude Containment Control

The traditional time-triggered control scheme is usually updated periodically, and periodic sampling may be a better control scheme in the view of control scheme design and problem analysis. However, when the system is working normally, periodic control will cause unnecessary energy consumption and actuator update.

To reduce the pressure of computation, communication, and actuators and meet the constraints of mass, volume, and power on the microsatellites, the event-triggered mechanism is introduced to solve coordinated attitude control of the microsatellite cluster, while to ensure the control performance of the microsatellite cluster, the triggering condition is used to adjust update period of controller and reduce the amount of computation and communication pressure. The design of the event-triggered coordinated attitude controller could be divided into two main steps. First is the setting of triggering condition. The triggering function with relative attitude error is designed to adjust update period of controller, state information is sampled, control law is computed, and actuator is updated if and only if the event is triggered. Second, the event-triggered distributed attitude controller is designed in the presence of model uncertainties and external disturbances, which makes followers to converge to the convex hull formed by the leaders’ orientation. And the Zeno behaviour is excluded. Triggering frequency is reduced and resource is saved as well as control performance of cluster system is ensured. The triggering condition is given in this section, and distributed attitude containment control protocol will be proposed in the next section.

First, to rewrite the dynamics equation of the satellite cluster into the form which is convenient for stability analysis, auxiliary variable is introduced:

Design the sliding variable :where α is a positive constant and is the weight of adjacency matrix of directed graph .

According to the properties of equations (5) and (9), the dynamics system could be rewritten as

In coordinated attitude control of the satellite cluster, is the inertia matrix with unknown parameters, . For any , there exists . is the approximation error for the ith follower, . Linear operator is shown as follows:

Regression function matrix:

Let the triggering time sequence of follower i be , where and . denotes the kth event time of the ith satellite.

Two types of state errors are defined to facilitate the design of the triggering condition:where is the estimation of at time .

According to Reference [35], the triggering condition can be defined as follows:where is the symmetric positive-definite matrix, , , , and .

Remark 1. is a time-dependent function. The introduction of can completely avoid the Zeno behaviour. The initial attitude and angular velocity , , of followers and attitude and angular velocity , , of leaders are given at initial time . The system continuously monitors the auxiliary state and uncertain parameter of each follower. Once the triggering condition of satellite i is satisfied, that is, , the event of satellite i is triggered immediately, triggering time is , state information is sampled, control law is computed, and actuators are updated. The state errors and are reset to 0. In nontriggering time, the controller uses state information of the last triggering instant, i.e., for , . Control law is not computed until the next event is triggered, and thus the control input is piecewise constant. Under the action of the event-triggered attitude containment controller, the attitude of followers asymptotically converges to the convex hull formed by the leaders’ orientation.
In this paper, the event-triggered distributed attitude control algorithm could be written in the following form:It is noteworthy that the control law does not need be computed until the next event.
Attitude dynamics equation of followers could be rewritten in the following form:Correspondingly, triggering time is defined asOnce , the event is triggered, state information is sampled, control law is computed and actuators are updated for satellite i, and the corresponding triggering function is less than 0 again.
The event-triggered control scheme is shown in Figure 1. The scenarios we consider include sampling time and triggering time; the sampling time is determined by the fixed clock frequency of the sensor, and the latter is controlled by the triggering condition based on the state. It is noteworthy that the time interval between two consecutive events is usually variable and can be equal to the minimum interval (that is, the sampling period). When the interval between two consecutive events is a fixed sampling period, the system degenerates into a traditional periodic sampling control. For the triggering conditions in this paper, continuous detection is necessary, which may require additional equipment and is a waste of communication and computing resources.

3.2. Event-Triggered Adaptive Sliding Mode Attitude Containment Control

For each follower i, the event-triggered adaptive sliding mode control protocol is designed as

Adaptation law of and is defined as follows:where is a symmetric positive-definite matrix and is a positive constant.

According to (9) and (10), system (17) can be rewritten as

The following conclusion shows that event-triggered adaptive sliding mode attitude control protocol (19), adaption laws (20) and (21), and triggering condition (15) can realize attitude containment control of microsatellite cluster system (5) under directed information topology. When , the attitude and angular velocity of followers asymptotically enter into the convex hull formed by leaders, that is, and . Moreover, it can be further proved that the Zeno behaviour will not occur in the microsatellite cluster system under the action of the proposed event-triggered attitude control strategy.

Theorem 1. If Assumptions 1 and 2 hold, under the action of triggering condition (15), event-triggered adaptive sliding mode attitude control protocol (19), and adaption laws (20) and (21), the attitude of followers asymptotically enters into the convex hull formed by leaders’ orientation, that is, and , and then the sufficient and necessary condition is that Assumption 3 holds.

(1) Proof. Sufficiency: The proof includes the following two consecutive steps. (i) State trajectory asymptotically converges to the sliding surface, that is, . (ii) In the case of , the state of followers will reach , asymptotically.

First, consider the following Lyapunov candidate:where and the estimation error is

Taking the derivative of , we can obtain

According to property (c) of the Euler–Lagrange equation, is a skew-symmetric matrix; substituting (24) into (25), we can obtain that

Substituting definitions of state errors (14) into (26), it can be obtained that

According to Assumption 2, we can get . Thus,

In the time interval , we can know that the system trajectories start from a region where . And is a symmetric positive-definite matrix; then, the upper bound for can be expressed as

It is worth noting that triggering condition (15) guarantees holds throughout the evolution of cluster system (5). Thus, we can get

Because , holds. For , analysing the right-hand side of (24), the upper bound of can be expressed as

Integrating both sides of (25), for any , we can get

Because , , we obtain that , and thus is bounded. According to (23), and are bounded for every . When Assumption 3holds, we can know exists according to Lemma 1. Then, (10) can be rewritten as follows:

Let , and equation (33) can be written as

According to Lemma 2, all eigenvalues of have positive real parts. According to Assumption 1, is bounded and . It thus follows that when , (34) is globally exponentially stable at the origin . We can conclude from Lemma 3 that cluster system (34) is input-to-state stable with respect to the input and the state . According to Lemma 3, because input is bounded, state is bounded. Thus, we can obtain that is bounded from (34). Equation (9) can be written as , and thus is bounded. It also follows from (10) that is bounded. Differentiating (9), we can get that is bounded.

According to property (b) of the Euler–Lagrange equation, , , , , and , are all bounded. Then, according to equation (11), we can obtain that is bounded. Because , we can obtain from (17) that is bounded.

Because (t) ≥ 0, and , then it can be derived from (32) thatwhich implies that is bounded, that is, .

, , and thus is uniformly continuous. According to Barbalat’s lemma [51], . Because (34) is input-to-state stable with respect to the input and state , it can be concluded that . As a result, it follows that , so similarly ; according to Assumption 1, , and thus . According to , we can obtain that . The attitude of followers asymptotically converges to the convex hull formed by leaders’ orientation.

According to Reference [52], we see that the system trajectories are attracted towards the sliding manifold as long as . This holds for all triggering instants whenever . As a result, the trajectories reach the neighbourhood of the sliding manifold due to and enter the region where the sign of changes before next triggering occurs. So, in the time interval, the decrement of cannot be guaranteed because zero crossing of occurs, and hence the trajectory increases. However, the sliding trajectories are ultimately bounded within this region. Therefore, the ultimate band is the region where . The size of this band can be calculated by finding the maximum deviation of sliding trajectory with zero crossing.

This ultimate region can now be derived as follows. We know that for any ,

This gives the maximum deviation of sliding variable from its immediate sampled value. Then, the maximum value of band can be obtained for the case and is given in (15).

According to (32), ; then, substituting (23) into it, we can obtain

Because , we can obtain that

Further,

That is,

Second, consider the following Lyapunov function candidate:

Taking the derivative of V2 along (34) gives

There exists a diagonal matrix with , , such that is symmetric positive definite. Because , holds. Then, (42) can be written as

Furthermore,

It can be obtained that

Let . Because , then according to (41), we can get

According to the aforementioned analysis, . According to Assumption 1, is bounded, and thus is bounded. Moreover, we assume that , and thus we can obtain that .

Substituting (41) into equation (45), we can get

From definitions of and , we can also obtain

Therefore, for any bounded initial states and , the states of (5) using equations (15) and (19)–(21) will always lie in the compact set .

(2) Proof. Necessity: Necessity is proven by contradiction. If Assumption 3 does not hold, then there must exist parts of followers which cannot receive information from any leader directly or indirectly (through other followers). According to Reference [53], the followers can be divided into two subsets, namely, one set with the followers that can receive the information from the leaders directly or indirectly, denoted by , and the follower set, which cannot receive information from leaders denoted by . The system Laplacian matrix L can also be in the following form:

Let the attitude of and be and . Auxiliary variables , are given with respect to and . Then, we can obtain

Because and are proved to converge to zero, the final states of and will depend on the equilibrium points. We can get from equation (50) that the trajectories of the followers in are independent of the leaders, and thus these followers cannot always converge to the convex hull formed by the leaders from any initial state. For the followers in , Assumption 3 holds. From the analysis of the sufficiency in this proof, the motion of followers in may also be affected by its neighbours including leaders and followers in , and thus there might exist some followers in that cannot converge to the convex hull formed by the leaders.

The proof is completed.

Remark 2. In order to avoid the chattering problem caused by the sign function, controller (19) can be substituted by saturation functionwhere is a small positive constant.

Remark 3. The attitude containment control problem with multiple leaders has been studied in the context. Obviously, if there exists one leader in control protocol (19), the multiple leader-follower coordinated attitude containment control problem becomes the coordinated attitude tracking problem with single leader.

3.3. The Absence of Zeno Behaviour

Zeno behaviour means the minimum time interval between two consecutive events is 0 and the event triggers infinite times in a finite time, which is forbidden in control tasks. Let denote interevent time interval of event-triggered attitude control protocol (19)–(21) and triggering condition (15). needs to be strictly positive to exclude the Zeno behaviour.

is satisfied in the context. But at , there may exist the case of . The attitude containment control is not reached at this moment. Therefore, the event is triggered in the following two cases:(1)If , the interevent time interval of this case is denoted as and triggering condition is .(2)If , the interevent time interval of this case is denoted as and triggering condition is .

Comparing the aforementioned two cases, we can know that for any , there exists , so it can be concluded that case (1) takes longer time, which implies . We just need to prove that there exists strictly positive to exclude the Zeno behaviour.

Theorem 2. If conditions of Theorem 1 hold, then the interevent time interval of triggering condition (15) has positive lower bound.

Proof. According to Reference [35], for any and , , the derivative of with respect to time satisfieswhere has been proved to be bounded. Let positive constant denote the maximum of , and we obtain .
Similar to equation (42), the derivative of satisfiesIt can be seen from the Section 3.1 that and are both bounded; according to Assumption 1 in Reference [35], we know that is bounded. Then, because and are bounded, is bounded according to adaption update law. Let positive constant denote the upper bound of , and can be obtained.
Let ; it follows thatAccording to case (2), the lower bound of interevent time interval can be derived:where B > 0. can get positive solution from equation (55). Thus, there exists time interval between two consecutive events, and the system avoids continuous control. It means that the interevent time interval of triggering condition (15) has positive lower bound.

4. The Influence of Information Topology Design on Followers’ Orientation

Compared with single leader/leaderless case, containment control has lower accuracy requirement for the final state of cluster members. However, in practical space missions, the orientation of satellites in the target area needs to meet certain constraints, such as multiple satellites observing the same orientation simultaneously, the observation field covering the entire target area, and so on.

It indicates in [46] that the steady state of each follower is a convex combination of all leaders’ states it can access, and the combination coefficient is a quantity related to the system Laplacian matrix. It can be concluded that the orientation of followers in the target area is determined by system information topology (including the weights that are assigned to edges).

In this section, approaches from graph theory to investigate influence of information topology on the distributions of followers are presented to provide a reference for orientation design of the microsatellite cluster.

4.1. The Constraints of Leader Reachable Set on Followers’ Orientation

Denote as the coefficient matrix of steady state of followers with respect to the state of leaders, where , . It reflects the relationship between system information topology and steady state of followers. According to the definition of reachable set [46], (1) if follower i is reachable from multiple leaders including , , the motion of i is affected by multiple leaders, and its stable state is correspondingly determined by the states of these leaders. (2) If follower i is only reachable from leader , , there exists a directed path from to i, that is, the motion of follower i will be only affected directly or indirectly by leader . (3) If follower i is not reachable from leader , , the motion of i will not be affected by leader , and thus the steady state of follower i is not related to .

is denoted as the reachable set of leader . is a set of all followers which is only reachable from leader . represents the unreachable set from leader . represents other follower set in without , and represents the subset without in .

Then, and can be written in the following block matrix forms:where .

Theorem 3. Assume all leaders have converged to the steady states . Then, under the action of triggering condition (15), adaptive sliding mode attitude containment control algorithm (19), and adaption laws (20) and (21), all followers of cluster system (5) will asymptotically converge to the steady state .

Proof. Sincewe haveand we can obtain thatFurther, each entry of is positive. Otherwise, , which means does not have neighbours in , which contradicts the definition of .
Now Theorem 3 has been proved.

4.2. The Constraints of Graph Differentiation on Followers’ Orientation

In microsatellite cluster flying missions, each member obtains the information from its neighbours through communication or relative state measurement. Due to the performance difference of sensors and communication equipment, as well as relative distance between members, neighbour satellite sets of each satellite are different. However, there may exist some commonalities among cluster members in information interaction, according to which cluster members can be divided into several subsets, and dynamic behaviour of cluster members belonging to the same subset may also have commonalities.

Although steady orientation of followers can be roughly estimated based on leader reachable sets, in some observation missions, there exist more constraints on followers’ orientation. The steady state of followers can be described using partition and partition of cluster information topology. These two graph theory tools can divide the cluster member into several subsets, and the number of neighbours in other subsets of cluster members belonging to the same subset has a certain commonality, which provides theoretical basis for information topology design of the microsatellite cluster.

At first, we prove that the satellites in the same cell partition belonging to partition have same steady state.

Theorem 4. For event-triggered adaptive sliding mode attitude containment control protocol (19), triggering condition (15), and adaption laws (20) and (21) of the microsatellite cluster system which satisfies , if system information topology has a cell partition , then all the followers that belong to the same cell have the same steady state.

Proof. According to Reference [46], suppose is a cell partition of with ; then, according to the cell partition of , by ordering the vertices appropriately, the graph Laplacian can be written in the following form:The block matrix at the lower left corner and the lower right corner is denoted by A, B, respectively. Since , is obtained. Then, by Lemma 2 in [46], . Assuming that each vertex in Ci has sij number of neighbours in , then, according to the definition of partition, each row sum of the submatrix equals sij. Combining with Corollary 1 in [46], it is given thatwhere is the cardinality of .
It follows thatBy using Lemma 3 in [46], it can be concluded that all the followers that belong to the same cell Ci(i = 1, …, k) have the same steady state.
Now Theorem 4 has been proved.
Then, another cell partition, partition, is provided to prove that satellites in the same cell partition belonging to partition have the same steady state.

Theorem 5. For event-triggered adaptive sliding mode containment control protocol (19)–(21) of microsatellite cluster which satisfies , if the system information topology has a cell partition with or , , then all the followers that belong to the same cell have the same steady state.

Proof. According to Reference [46], construct the deduced unweighted graph of as follows: use a follower to represent cell and draw an edge from to is a neighbour cell of . Then, has m reaches according to Assumption 3. Let be the Laplacian matrix of the deduced graph ; it can be derived that with if has unique solution. Finally, it can be verified that is the solution of .
Now Theorem 5 has been proved.

5. Simulation Results

5.1. Simulation Results and Analysis of Attitude Containment Control

In this section, simulations for multiple leader-follower satellite cluster are presented to illustrate the effectiveness of the proposed control protocol and information topology design. Suppose that in microsatellite cluster observation mission, six satellites, denoted by , are needed to obtain observation information from three different directions. Suppose the target area is defined by three attitude orientations, and the leaders are stationary relative to the reference frame. The attitude of leaders, respectively, is

The nominal inertia of six followers, respectively, is

The initial estimation parameter is .

The initial state of followers is, respectively, shown in Table 1.

System information topology is shown in Figure 2.

Control gain coefficients and adaptive parameters are

Triggering parameters are .

A microsatellite cluster in LEO is mainly affected by the gravity gradient torque, while the disturbances such as the solar radiation pressure torque will be dominant for a cluster in high-Earth orbits such as the geostationary orbit. All these torques are slowly varying and can be treated as signals composed of constants and periodic trigonometric functions. Taking into account these factors, the disturbances are chosen aswheredenote the constant disturbing torques andis used to simulate the residual disturbances including aerodynamic torques, solar radiation torques, and similar effects.

5.1.1. Event-Triggered Attitude Containment Control

Simulation results of satellite cluster are shown in Figure 3. Figures 3(a) and 3(b) are the curves of relative attitude and relative angular velocity over time, respectively. It can be seen that the relative attitude converges to the convex hull formed by the leaders at about 600 s, and relative angular velocity converges to 0 within 700 s. The followers can converge to the convex hull even though there exists large disturbance torque. The interevent time of each follower is shown in Figures 3(c) and 3(d). At the initial stage, the state of cluster members is far from the desired state, then the event is triggered frequently, and the update of control input is relatively frequent, but when the system asymptotically converges to the desired state, fewer events are triggered, and interevent time increases and tends to be stable finally. We can find that if satellites are influenced by periodic disturbance, the interevent time also changes periodically. The control torques acting on each satellite are shown in Figure 3(e). It can be seen that control torques acting on the followers are piecewise function, and the control input is only updated at the next triggering instant. The control torques are limited within the range of . Attitude trajectories of followers are shown in Figure 3(f), from which we can see that followers asymptotically converge to the convex hull formed by the leaders. Figures 3(g) and 3(h) show the event error and triggering threshold for each follower in 200 s. It can be seen that if event error exceeds threshold, the event is triggered and the state is sampled.

In a word, under the action of the event-triggered adaptive sliding mode attitude controller, the attitude of each satellite asymptotically converges to the convex hull formed by the leaders, and angular velocity asymptotically converges to 0.

5.1.2. Traditional Time-Triggered Attitude Containment Control

According to Reference [53], the time-triggered adaptive sliding mode attitude containment control algorithm is

The simulation results of time-triggered distributed attitude adaptive sliding mode controller are shown in Figure 4. The attitude of the cluster system converges to the convex hull formed by the leaders at about 600 s, and attitude angular velocity converges to 0 within 700 s. The controller is continuously updated while the convergence rate and control accuracy are not better than the event-triggered one.

It can be clearly seen from Figures 3 and 4 that both time-triggered and event-triggered control strategies can realize attitude containment. In addition, it is noteworthy that event-triggered containment control is updated in an aperiodic manner, while time-triggered control is updated at a fixed interval of 0.01 s. Within 1200 s, the sampling and control input update times of the time-triggered control method are 120,000, while the event of each follower, respectively, is 391, 445, 405, 354, 402, and 480, from which the update frequency of control action is greatly reduced by the event-triggered control strategy. Compared with time-triggered attitude containment control protocol, event-triggered one can effectively reduce the control input update frequency while ensuring the similar control performance. Through the reasonable selection of control parameters and triggering function, event-triggered control can guarantee the convergence rate and control accuracy and reduce the amount of computation and communication, as well as the update frequency of actuators.

5.2. The Influence of Information Topology on Follower’s Orientation

Suppose in the Earth observation mission of the microsatellite cluster, twelve satellites are used to obtain the observation information of three different orientations. In order to meet the accuracy requirement, two satellites can be used to observe the same azimuth at the same time. Suppose the target area is formed by three attitude orientations. The initial state of leaders, respectively, is

Nominal inertia of followers, respectively, is

The initial estimation parameter of followers is .

The constant disturbance acting on the followers 1–6 and periodic disturbances are the same as Section 5.1, while constant disturbance acting on followers 7–9 is

The initial state of followers is, respectively, shown in Table 2.

Information topology of the cluster system is shown in Figure 5.

The information topology structure of the microsatellite cluster system is shown in Figure 5. It can be seen that Figure 5 belongs to partition. Satellites 1, 2, 7, and 8 belong to the same cell, 3 and 4 belong to the same cell, and 5 and 6 belong to the same cell. According to Theorem 4, satellites belonging to the same cell converge to the same azimuth under the action of control torques in the presence of model uncertainties and external disturbances. This property can be used to observe the same azimuth. Simulation results are shown in Figure 6. The attitude of follower converges to the convex hull formed by the leaders within 600 s, and attitude angular velocity converges to 0 within 700 s. The control torques are piecewise constant. Satellites 1 and 2 converge to the same azimuth, and 7 and 8 converge to same azimuth, respectively. Because satellites in C5 only receive information from C3, 3, 4, 5, and 6 converge to the same azimuth. Within 1200 s, the sampling and control input update times of each follower, respectively, are 232, 247, 522, 420, 162, 142, 196, 288, and 315.

6. Conclusion

In this paper, an event-triggered adaptive sliding mode attitude containment control protocol is proposed in the framework of the Euler–Lagrange model for the attitude containment control problem of the microsatellite cluster system. Considering the constraints of resource and power on the microsatellite cluster system, the event-triggered control strategy is introduced into the attitude containment control problem of the microsatellite cluster. The triggering condition consisting of state-dependent function and the time-dependent function is given to adjust controller update period and exclude the Zeno behaviour. If and only if the triggering condition is satisfied, state information is sampled, control law is computed, and actuators are updated. Then, under directed topology, the event-triggered adaptive sliding mode attitude containment control algorithm is proposed, which makes attitude of followers asymptotically enter into the convex hull formed by leaders’ orientation in the presence of inertia uncertainties and external disturbances. Numerical simulation is carried out to verify the effectiveness of the proposed event-triggered distributed attitude containment control algorithm. Then, compared with the time-triggered one, it can be seen that while ensuring the control performance of the cluster system, the designed event-triggered attitude containment controller only updates control law in the triggering instant, which effectively reduces the amount of computation and communication and update frequency of actuators and saves resources on board.

Furthermore, the influence of cluster information topology structure design on the stable state of containment control algorithm is also studied. An appropriate information topology is designed to meet the attitude orientation requirements of containment control. It is shown that the steady state of each follower is a convex combination of all leaders’ states it can access. The cell partition of cluster information topology is given based on the number of satellite’s neighbours, and it is proved that the cluster members belonging to the same cell have the same stable state. However, there is no requirement for the information interaction between satellites inside the cell, and the information links inside the cell will not affect the stable state of cluster members. It provides a theoretical basis for the design of information topology of the microsatellite cluster system. Numerical simulation is conducted to verify the influence of information topology on steady state of the microsatellite cluster system.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.