Hostname: page-component-78c5997874-m6dg7 Total loading time: 0 Render date: 2024-11-05T01:05:04.930Z Has data issue: false hasContentIssue false

Distributed fixed-time control for six-degree-of-freedom spacecraft formation with event-triggered observer

Published online by Cambridge University Press:  13 June 2023

J. Zhang
Affiliation:
Space Control and Inertial Technology Research Center, School of Astronautics, Harbin Institute of Technology, Harbin, China
H. Xia*
Affiliation:
Space Control and Inertial Technology Research Center, School of Astronautics, Harbin Institute of Technology, Harbin, China
Z. Li
Affiliation:
Shanghai Aerospace Control Technology Institute, Shanghai Academy of Spaceflight Technology, Shanghai, China
*
*Corresponding author. Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The paper addresses the six-degree-of-freedom coupled control problem for spacecraft formation flying subject to actuator saturation and input quantisation whilst considering limited communication resources. Firstly, a novel event-triggered distributed observer without continuous communications is presented to recover the information of the virtual leader. Remarkably, by embedding a hyperbolic tangent function-based nonlinear term into the triggering condition, the event-based observer realises a more reasonable trigger threshold. Subsequently, an adding-a-power-integrator-based fixed-time control algorithm is proposed for the follower spacecraft. Further, the control scheme ingeniously compensates for the actuator saturation and the input quantisation problems without embedding auxiliary systems. Finally, numerical simulations are carried out to highlight the advantages of the theoretical results.

Type
Research Article
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of Royal Aeronautical Society

Nomenclature

${\boldsymbol{{x}}_i}$

ith follower spacecraft’s position states vector

${\boldsymbol{{v}}_i}$

ith follower spacecraft’s velocity states vector

${\boldsymbol{{u}}_i}$

ith follower spacecraft’s control inputs vector

$\delta ,{u_{min}}$

quantiser parameters

${\boldsymbol\chi _s}$

decomposition matrix of the actuator saturation

${\boldsymbol\chi _q}$

decomposition matrix of the input quantisation

${\chi _{min}}$

lower bound for elements in the decomposition matrix ${\boldsymbol\chi _{\!s}}$

${\boldsymbol{{x}}_0}$

the virtual leader’s position states vector

${\boldsymbol{{v}}_0}$

the virtual leader’s velocity states vector

${\hat{\boldsymbol{{v}}}_{i,0}}$

the ith follower spacecraft’s estimation of ${\boldsymbol{{v}}_0}$

${\widetilde{\boldsymbol{{v}}}_{i,0}}$

estimation error for the ith follower spacecraft, ${\widetilde{\boldsymbol{{v}}}_{i,0}} = {\hat{\boldsymbol{{v}}}_{i,0}} - {\boldsymbol{{v}}_0}$

$\boldsymbol{{H}}$

the communication topology matrix of the formation

${\boldsymbol\zeta _i}$

the auxiliary variable related to the observational errors

$\alpha ,\beta $

fixed-time parameters

${o_1},{o_2},{o_3},{o_4}$

observer gains

${k_{o1}},{k_{o2}},{k_{o3}},\iota $

parameters of the event-trigger mechanism

${B_3}$

the upper bound of $\left\| {{{\ddot{\boldsymbol{{x}}}}_0}} \right\|$

${\boldsymbol\Delta _{i,0}}$

the ith follower spacecraft’s expected deviation formation vector

${\boldsymbol{{e}}_{xi}}$

tracking errors for ith follower spacecraft, ${\boldsymbol{{e}}_{xi}} = {\boldsymbol{{x}}_i} - {\boldsymbol\Delta _{i,0}} - {\boldsymbol{{x}}_0}$

${\boldsymbol{{e}}_{vi}}$

ith follower spacecraft’s tracking errors derivative, ${\boldsymbol{{e}}_{vi}} = {\boldsymbol{{v}}_i} - {\dot{\boldsymbol\Delta}_{i,0}} - {\boldsymbol{{v}}_0}$

${\boldsymbol\sigma _i},{\boldsymbol\gamma _i},{\boldsymbol\gamma _{di}},{\boldsymbol\kappa _i}$

the auxiliary variable related to the tracking errors

${k_1},{k_2},{k_3},{k_4},{k_5}$

controller gains

Abbreviations

SFF

spacecraft formation flying

ETM

event-trigger mechanism

ASKAE

attitude station keeping absolute error

PSKAE

position station keeping absolute error

1.0 Introduction

The spacecraft formation flying (SFF) system has drawn extensive attention owing to its costs, robustness and flexibility [Reference Mauro, Spiller, Bevilacqua and D’Amico1] in recent years, which can overcome the limitations of using a single spacecraft for mission goal accomplishment. Numerous space missions by deploying SFF system such as monitoring [Reference Shaw, Miller and Hastings2], remote sensing [Reference Ren and Beard3], and on-orbit services [Reference Ahn and Kim4] have been successfully launched.

For the SFF tracking control problem, the convergence rate is an important performance indicator reflecting the proposed algorithm’s effectiveness. Most of the existing studies [Reference Jiakang, Guangfu and Qinglei5, Reference Wu, Wang and Poh6] achieve asymptotic stability results only, of which convergence time tends to infinity. In contrast to the asymptotic control algorithm, the finite-time control algorithm can provide faster convergence rate, higher precision, and better robustness, with the adding a power integrator technique [Reference Zhou, Hu and Friswell7] the terminal sliding mode theorem [Reference Zou and Kumar8], and the homogeneity theorem [Reference Zou, de Ruiter and Kumar9] as the primary design methods. However, the prior estimation of the finite-time control algorithm’s convergence time is decided by the initial states. As an extension of finite-time method, the fixed-time control algorithm was introduced in [Reference Polyakov10], where the settling time is uniformly bounded with regard to the initial states. One of the basic synchronisation algorithms for the formation keeping and reconstruction problem is leader-follower formation, where the leader runs in a predetermined trajectory while the followers need coordinate. Following the above frameworks, Ren first proposed an SFF control algorithm [Reference Ren11] considering the two cases of whether all follower spacecraft can access the reference trajectory. In light of this, Gao proposed a fixed-time coordinated algorithm for the SFF system [Reference Han, Yuanqing, Zhang and Zhang12]. However, most scholars concentrated their research on the attitude coordinated control or modeled the translational and rotational motions separately, ignoring the mutual coupling between the orbit and attitude.

In practice, continuous interactive communication of formation spacecraft is challenging due to the limitation of onboard resources, especially communication channel bandwidth. To this end, the seminal work [Reference Tabuada13, Reference Dimarogonas, Frazzoli and Johansson14] reported an event-triggered method premised on asynchronous, aperiodic communication. Thereafter, the event-based control has drawn considerable attentions [Reference Nowzari, Garcia and Cortés15Reference Di, Li, Guo, Xie and Wang20] (just to name a few). In particular, Nowzari provides a comprehensive account of the motivation behind applying event-triggered strategies in multi-agent systems [Reference Nowzari, Garcia and Cortés15]. Further, continuous communication can be avoided both in the fixed-time controller algorithm update and in the triggering condition monitoring in [Reference Liu, Zhang, Yu and Sun17]. Despite this, it is worth noting that the lower limit of the trigger threshold for most existing works has a stable value, whether in the transient or steady state of the system. In other words, the lower limit of the trigger threshold does not increase when the error is large; on the contrary, the lower limit does not decrease when the error is small, which results in a waste of communication resources. Bearing this in mind, by utilising the characteristics of the hyperbolic tangent function, the change of the lower limit of the trigger threshold with the error is discussed in [Reference Di, Li, Guo, Xie and Wang20]. Nevertheless, it is noticed that the event-trigger mechanism (ETM) in [Reference Di, Li, Guo, Xie and Wang20] depends on continuous communications. As such, how to design an ETM for multiple spacecraft systems, which has a more flexible trigger threshold lower limit and eliminates continuous communication, still remains open and awaits a breakthrough.

Another practical problem for spacecraft is that the onboard actuators are susceptible to suffering from quantisation errors and magnitude constraints. Data transmission between the attitude control module and actuator module would introduce input quantisation errors, which may lead to the degradation of the control performance. To alleviate this concern, control algorithms that take account of input quantisation were investigated for spacecraft systems. In [Reference Zhou, Wen and Yang21] and [Reference Liu, Wang, Zhang and Chen22], the functional relationship between the control input and the quantised input is established by the non-linear decomposition of the quantised signal. Suffering from the piece-wise quantised input, most relevant literatures have used the above or similar methods to overcome the technique difficulty. Additionally, actuators have saturation constraints due to physical structure and energy limitations, which will severely limit the performance of the closed-loop system if the actuator is always saturated. Various methods have been designed for the input saturation problem, such as adding a neural network-based compensator [Reference Hu and Xiao23], introducing an auxiliary system [Reference Han, Yuanqing, Zhang and Zhang12], constructing a command prefilter [Reference An, Liu, Wang and Wu24], and so on. However, the above methods of dealing with actuator saturation dramatically increased the complexity of the control algorithm. Although scholars have extensively studied the input quantisation and actuator saturation problem, nearly no research has addressed the fixed-time anti-saturation and anti-quantisation control problem for 6-DOF SFF.

Statistically, most SFF missions are carried out by small spacecraft operating in low Earth orbit (LEO) with a mass of less than 500kg [Reference Di Mauro, Lawn and Bevilacqua25]. The light mass of the spacecraft means limited onboard resources, such as communication bandwidth, actuator capability, and energy. Given these onboard resource constraints, it is worth investigating the design of a distributed controller to meet the demands of increasingly complicated SFF missions for the spacecraft’s rapid and precise manoeuver capabilities. In this paper, we consider the onboard resource constraints for the problem of limited communication resources, actuator saturation, and input quantisation. Here, limited communication resources implies that only a portion of the spacecraft can access the virtual leader. Besides, for SFF operating in LEO, simultaneous attitude and constellation control are essential for timely identification and counteracting time-critical orbit spacing violations. In other words, designing a six-degree-of-freedom controller for the SFF mission that considers the attitude-orbit coupling is necessary. Bearing the discussion above in mind, we address the attitude and position coupled tracking fixed-time control problem for SFF systems considering limited communication resources, actuator saturation and input quantisation. A distributed event-triggered observer is first proposed to recover the virtual leader’s velocity states. More importantly, the event-based observer realises a more reasonable trigger threshold and less communication resource usage by embedding a hyperbolic tangent function-based nonlinear term into the triggering condition. Subsequently, a fixed-time distributed control algorithm is developed for the SFF system, handling the actuator saturation and the input quantisation problems ingeniously without embedding auxiliary systems. The main contributions of this paper are summarised as follows.

  1. (1) A novel distributed event-based fixed-time observer is developed to estimate the states of the virtual leader for SFF system. Compared with the traditional observer presented in [Reference Han, Yuanqing, Zhang and Zhang12], the observer proposed in this paper alleviate the chattering caused by symbolic functions. Besides, an additional nonlinear term is adopted by the observer to guarantee the asymptotic stability of the whole closed-loop system.

  2. (2) Compared with the existing ETM in [Reference Liu, Zhang, Yu and Sun17], the ETM of the observer proposed in this paper realises a more reasonable trigger threshold and a lower bound of trigger threshold with characteristics of the exponential function and the hyperbolic tangent function in the transient and steady state.

  3. (3) A distributed fixed-time attitude control scheme is established in the presence of the actuator saturation and the input quantisation problems by utilising the adding a power integrator technique. In contrast to the existing methods in [Reference Hu, Shi and Shao26], the proposed control scheme compensates the actuator saturation and the input quantisation problems ingeniously without embedding auxiliary systems.

The remainder of this paper is as follows. Section II introduces some preliminaries. The observer-based distributed control scheme is provided in Section III. Section IV offers numerical simulations. Finally, Section V gives a conclusion.

2.0 Preliminaries

2.1 Notations and lemmas

Define $\left\| \cdot \right\|$ as the induced norm of a matrix or the Euclidean norm of a vector. ${\boldsymbol{{I}}_n}$ denotes the $n \times n$ identity matrix. $ \otimes $ is defined as the Kronecker product. ${\lambda _{\min }}\!\left( \cdot \right)$ and ${\lambda _{\max }}\!\left( \cdot \right)$ represent the minimum and maximum eigenvalues of a symmetric matrix. The basic operations of the dual number, quaternion and dual quaternion are defined in the appendix. Given a vector $\boldsymbol{{x}} = {\left[ {{x_1},{x_2}, \ldots ,{x_n}} \right]^{\rm{T}}} \in {\mathbb{R}^n}$ and $\alpha \gt 0$ , we denote ${\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol{{x}}) = {\left[ {{\rm{sign}}\!\left( {{x_1}} \right){{\left| {{x_1}} \right|}^\alpha },{\rm{sign}}\!\left( {{x_2}} \right){{\left| {{x_2}} \right|}^\alpha }, \ldots ,{\rm{sign}}\!\left( {{x_n}} \right){{\left| {{x_n}} \right|}^\alpha }} \right]^{\rm{T}}}$ and $|\boldsymbol{{x}}| = {\left[ {\left| {{x_1}} \right|,\left| {{x_2}} \right|, \ldots ,\left| {{x_n}} \right|} \right]^{\rm{T}}}$ , where ${\rm{sign}}( \cdot )$ is the signum function. ${\rm{diag}}\!\left( {{x_1},{x_2}, \ldots ,{x_n}} \right)$ denotes a block-diagonal matrix. ${\rm{max}}\!\left( {{x_1},{x_2}, \ldots ,{x_n}} \right)$ and ${\rm{min}}\!\left( {{x_1},{x_2}, \ldots ,{x_n}} \right)$ denote the maximum and minimum value in $\left( {{x_1},{x_2}, \ldots ,{x_n}} \right)$ , respectively.

In order to get the main results of this paper, the following lemmas are introduced.

Lemma 1. [Reference Polyakov27] If the time derivative of a Lyapunov function $V$ satisfies $\dot V \le - {\kappa _1}{V^{{\rho _1}}} - {\kappa _2}{V^{{\rho _2}}}$ where ${\kappa _1} \gt 0,{\kappa _2} \gt 0,0 \lt {\rho _1} \lt 1$ , and ${\rho _2} \gt 1$ are some constants, the value of $V$ converges to zero in fixed-time. The settling time satisfies $T \le 1/\left[ {{\kappa _1}\!\left( {1 - {\rho _1}} \right)} \right] + 1/\left[ {{\kappa _2}\!\left( {{\rho _2} - 1} \right)} \right]$ .

Lemma 2. [Reference Zou and Kumar28] For any $x \in \mathbb{R},y \in \mathbb{R}$ , the inequality $\left| {{\rm{si}}{{\rm{g}}^\alpha }(x) - {\rm{si}}{{\rm{g}}^\alpha }(y)} \right| \le {2^{1 - \alpha }}|x - y{|^\alpha }$ is tenable if $\alpha \in (0,1]$ holds.

Lemma 3. [Reference Zou and Kumar28] For any $x \in \mathbb{R},y \in \mathbb{R}$ , the following inequality is tenable when $\alpha \geqslant 1$ . $|x - y{|^\alpha } \le {2^{\alpha - 1}}\left| {{\rm{si}}{{\rm{g}}^\alpha }(x) - {\rm{si}}{{\rm{g}}^\alpha }(y)} \right|$

Lemma 4. [Reference Yu, Yu, Shirinzadeh and Man29] For any positive constant ${a_1},{a_2}, \cdots ,{a_n}$ , the following inequality is tenable if $p \in \left( {0,2} \right)$ holds. ${\left( {a_1^2 + a_2^2 + \cdots + a_n^2} \right)^p} \le {\left( {a_1^p + a_2^p + \cdots + a_n^p} \right)^2}$

Lemma 5. [Reference Qian and Lin30] For any $x \in \mathbb{R},y \in \mathbb{R},c \gt 0,d \gt 0,\gamma \gt 0$ , the inequality $|x{|^c}|y{|^d} \le \frac{{c\gamma |x{|^{c + d}}}}{{c + d}} + \frac{{d|y{|^{c + d}}}}{{{\gamma ^{c/d}}(c + d)}}$ is tenable.

Lemma 6. [Reference Han, Yuanqing, Zhang and Zhang12] For any ${x_i} \in \mathbb{R},i = 1,2, \ldots ,n$ , the following inequality is tenable if $v \in (0,1]$ . ${\left( {\sum\nolimits_{i = 1}^n \left| {{x_i}} \right|} \right)^v} \le \sum\nolimits_{i = 1}^n {\left| {{x_i}} \right|^v} \le {n^{1 - v}}{\left( {\sum\nolimits_{i = 1}^n \left| {{x_i}} \right|} \right)^v}$ .

Lemma 7. [Reference Zou, de Ruiter and Kumar9] For any ${x_i} \in \mathbb{R},i = 1,2, \ldots ,n$ , the following inequality is tenable if $v \gt 1$ . $\sum\nolimits_{i = 1}^n {\left| {{x_i}} \right|^v} \le {\left( {\sum\nolimits_{i = 1}^n \left| {{x_i}} \right|} \right)^v} \le {n^{v - 1}}\sum\nolimits_{i = 1}^n {\left| {{x_i}} \right|^v}$ .

Lemma 8. [Reference Tan, Yu and Man31] For matrices $\boldsymbol{{X}},\boldsymbol{{Y}}$ with equal dimensions and positive constant $\eta $ , the following inequality satisfies: $2{\boldsymbol{{X}}^{{\text{T}}}}\boldsymbol{{Y}} \le \frac{1}{\eta }{\boldsymbol{{X}}^{{\text{T}}}}\boldsymbol{{X}} + \eta {\boldsymbol{{Y}}^{{\text{T}}}}\boldsymbol{{Y}}$ .

2.2 Mathematical model

The coordinate system used in this paper is shown in Fig. 1, which defines the Earth-centred inertial frame ${O_E}{X_I}{Y_I}{Z_I}$ and the jth spacecraft body frame ${O_{Bj}}{X_{Bj}}{Y_{Bj}}{Z_{Bj}}$ . ${O_E}$ is the centre of the Earth and ${O_{Bj}}$ is the centre of mass of the jth deputy.

Figure 1. Coordinate system.

Based on dual quaternion, the states of ith spacecraft can be described as ${\hat{\boldsymbol{{q}}}_i}$ and ${\hat{\boldsymbol\omega}_i}$ , where ${\hat{\boldsymbol{{q}}}_i} = {\boldsymbol{{q}}_i} + \varepsilon \frac{1}{2}{\boldsymbol{{q}}_i} \circ {\boldsymbol{{r}}_i}$ , ${\hat{\boldsymbol\omega}_i} = {\boldsymbol\omega _i} + \varepsilon\!\left( {{{\dot{\boldsymbol{{r}}}}_i} + {\boldsymbol\omega _i} \times {\boldsymbol{{r}}_i}} \right)$ , $\varepsilon $ is the dual unit that satisfies ${\varepsilon ^2} = 0,\varepsilon \ne 0$ . ${\boldsymbol{{q}}_i}$ is the quaternion describing the rotation of the ith spacecraft body frame relative to the inertial frame. ${\boldsymbol{{r}}_i} = {\left[ {0,{{\left( {{{\vec{\boldsymbol{{r}}}}_i}} \right)}^{\rm{T}}}} \right]^{\rm{T}}} \in {\mathbb{R}^4}$ , ${\vec{\boldsymbol{{r}}}_i} \in {\mathbb{R}^3}$ denotes the position of the ith spacecraft expressed in the spacecraft body frame. ${\boldsymbol\omega _i} = {\left[ {0,{{\left( {{{\vec{\boldsymbol\omega}}_i}} \right)}^{\rm{T}}}} \right]^{\rm{T}}} \in {\mathbb{R}^4}$ , ${\vec{\boldsymbol\omega}_i} \in {\mathbb{R}^3}$ is the angular velocity of the ith spacecraft body frame relative to the inertial frame expressed in the body frame. The dynamic equations of six-DOF motion are [Reference Wu, Liu and Han32]

(1) \begin{align}{\dot{\hat{\boldsymbol{{q}}}}_i} = \frac{1}{2}{\hat{\boldsymbol{{q}}}_i} \circ {\hat{\boldsymbol\omega}_i}\end{align}
(2) \begin{align}\dot{\hat{\boldsymbol\omega}}_{i} = M_{i}^{-1}{\ast}\!\left\{\hat{\boldsymbol{{F}}}_{ig} + \boldsymbol{{sat}}\!\left[ \boldsymbol{{Q}}\!\left( \hat{\boldsymbol{{F}}}_{ic} \right) \right] + \hat{\boldsymbol{{F}}}_{ip} - \hat{\boldsymbol\omega}_{i} \times \left( \boldsymbol{{M}}_{i}\,{\ast}\,\hat{\boldsymbol\omega}_{i} \right) \right\}\end{align}

where ${\hat{\boldsymbol{{F}}}_{ig}} = {\boldsymbol{{f}}_{ig}} + \varepsilon {\boldsymbol\tau _{ig}} = - \mu {m_i}{\boldsymbol{{r}}_i}/{\left\| {{\boldsymbol{{r}}_i}} \right\|^3} + \varepsilon 3\mu {\boldsymbol{{r}}_i} \times \left( {{{\overline{\boldsymbol{{J}}} }_i}{\boldsymbol{{r}}_i}} \right)/{\left\| {{\boldsymbol{{r}}_i}} \right\|^5}$ denotes the gravity dual force, ${\hat{\boldsymbol{{F}}}_{ic}} = {\boldsymbol{{f}}_{ic}} + \varepsilon {\boldsymbol\tau _{ic}}$ denotes the control dual force that is yet to be designed, $\boldsymbol{{sat}}\!\left[ {\boldsymbol{{Q}}\!\left( {{{\hat{\boldsymbol{{F}}}}_{ic}}} \right)} \right]$ denotes the control dual force with actuator saturation and input quantisation, ${\hat{\boldsymbol{{F}}}_{ip}} = {\boldsymbol{{f}}_{ip}} + \varepsilon {\boldsymbol\tau _{ip}}$ denotes the disturbing dual force. $\mu $ is the gravitational parameter of Earth. Define that ${m_i}$ and ${\boldsymbol{{J}}_i} \in {\mathbb{R}^{3 \times 3}}$ are the mass and inertia matrix of the ith deputy, respectively. The matrix ${\boldsymbol{{M}}_i}$ is

(3) \begin{align}{\boldsymbol{{M}}_i} = \left[ {\begin{array}{c@{\quad}c}{{{\bf{0}}_{4 \times 4}}} & {{{\overline{\boldsymbol{{m}}} }_i}}\\[5pt] {{{\overline{\boldsymbol{{J}}}}_i}} & {{{\bf{0}}_{4 \times 4}}}\end{array}} \right] \in {\mathbb{R}^{8 \times 8}}\end{align}

where

\begin{align*}{\bar{\boldsymbol{{J}}}_i} = \left[ {\begin{array}{c@{\quad}c}1 & {{{\bf{0}}_{1 \times 3}}}\\[5pt] {{{\bf{0}}_{3 \times 1}}} & {{\boldsymbol{{J}}_i}}\end{array}} \right] \in {\mathbb{R}^{4 \times 4}},{\overline{\boldsymbol{{m}}}_i} = \left[ {\begin{array}{c@{\quad}c}1 & {{{\bf{0}}_{1 \times 3}}}\\[5pt] {{{\bf{0}}_{3 \times 1}}} & {{m_i}{\boldsymbol{{I}}_{3 \times 3}}}\end{array}} \right] \in {\mathbb{R}^{4 \times 4}}\end{align*}

Remark 1. As an extension of quaternion, the dual quaternion is one of the most commonly used methods for modeling the six-degrees-of-freedom coupled spacecraft. In contrast with other six-degree-of-freedom modeling methods, such as the Vectrix approach [Reference Ploen, Hadaegh and Scharf33] and the Cayley description [Reference Sinclair, Hurtado and Junkins34], the dual quaternion approach is more explicit in physical meaning and compact in its mathematical description. In particular, the introduction of dual numbers allows us to describe the six-degrees-of-freedom coupled dynamics with a unified equation, implying that we do not have to design separate controllers and Lyapunov functions for the rotational and translational parts of spacecraft.

Equation (1) and Equation (2) can be rewritten as

(4) \begin{align}\left\{ {\begin{array}{l}{{{\dot{\boldsymbol{{x}}}}_i} = {\boldsymbol{{v}}_i}}\\[5pt] {{{\dot{\boldsymbol{{v}}}}_i} = {\boldsymbol{{f}}_i}\!\left( {{\boldsymbol{{x}}_i},{\boldsymbol{{v}}_i}} \right) + {\boldsymbol{{g}}_i}\!\left( {{\boldsymbol{{x}}_i}} \right)\boldsymbol{{sat}}\!\left[ {\boldsymbol{{Q}}\!\left( {{\boldsymbol{{u}}_i}} \right)} \right] + {\boldsymbol{{d}}_i}}\end{array}} \right.\end{align}

where ${\boldsymbol{{x}}_i} = {\hat{\boldsymbol{{q}}}_i} = {\left[ {\begin{array}{l}{{{\hat{\boldsymbol{{q}}}}_{ir}}^{\rm{T}}}\quad {{{\hat{\boldsymbol{{q}}}}_{id}}^{\rm{T}}}\end{array}} \right]^{\rm{T}}},{\boldsymbol{{v}}_i} = {\dot{\hat{\boldsymbol{{q}}}}_i} = {\left[ {\begin{array}{l}{{{\dot{\hat{\boldsymbol{{q}}}}}_{ir}}^{\rm{T}}}\quad {{{\dot{\hat{\boldsymbol{{q}}}}}_{id}}^{\rm{T}}}\end{array}} \right]^{\rm{T}}}$ , ${\boldsymbol{{f}}_i}\!\left( {{\boldsymbol{{x}}_i},{\boldsymbol{{v}}_i}} \right) = \dot{\boldsymbol\Gamma}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol\Gamma ^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol{{v}}_i} - \boldsymbol\Gamma\!\left( {{\boldsymbol{{x}}_i}} \right)\boldsymbol{{M}}_i^{ - 1}\boldsymbol{{Z}} ({{\boldsymbol{{x}}_i}, {\boldsymbol{{v}}_i}}) \boldsymbol\Gamma\!\left( {{\boldsymbol{{x}}_i}} \right)\boldsymbol{{M}}_i^{ - 1}\boldsymbol{{G}}$ , ${\boldsymbol{{g}}_i}\!\left( {{\boldsymbol{{x}}_i}} \right) = {\boldsymbol\Gamma}\!\left( {{\boldsymbol{{x}}_i}} \right)\boldsymbol{{M}}_i^{ - 1}$ , ${\boldsymbol{{u}}_i} = {\left[ {\begin{array}{l}{{\boldsymbol{{f}}_{ic}}^{\rm{T}}}\quad {{\boldsymbol\tau _{ic}}^{\rm{T}}}\end{array}} \right]^{\rm{T}}}$ , ${\boldsymbol{{d}}_i} = {\boldsymbol{{g}}_i}\!\left( {{\boldsymbol{{x}}_i}} \right){\left[ {\begin{array}{l}{{\boldsymbol{{f}}_{ip}}^{\rm{T}}}\quad {{\boldsymbol\tau _{ip}}^{\rm{T}}}\end{array}} \right]^{\rm{T}}}$ , $\boldsymbol{{Z}}\!\left( {{\boldsymbol{{x}}_i},{\boldsymbol{{v}}_i}} \right) = \boldsymbol\Gamma \left[ {{\boldsymbol\Gamma ^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol{{v}}_i}} \right]{\boldsymbol{{M}}_i}{\boldsymbol\Gamma ^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol{{v}}_i} - \boldsymbol\Gamma \left\{ {{{\left[ {{\boldsymbol{{M}}_i}{\Gamma ^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol{{v}}_i}} \right]}^*}} \right\}{\left[ {{\boldsymbol\Gamma ^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right){\boldsymbol{{v}}_i}} \right]^*}$ , $G = {\left[ {\begin{array}{l}{{\boldsymbol{{f}}_{ig}}^{\rm{T}}}\quad {{\boldsymbol\tau _{ig}}^{\rm{T}}}\end{array}} \right]^{\rm{T}}}$ , the transformation matrix $\boldsymbol\Gamma\!\left( {{\boldsymbol{{x}}_i}} \right)$ is defined in Equation (63). ${\boldsymbol{{g}}_i}\!\left( {{\boldsymbol{{x}}_i}} \right)$ is invertible because $\boldsymbol\Gamma ({\boldsymbol{{x}}_i})$ and ${\boldsymbol{{M}}_i}$ are invertible.

The quantiser in Equation (4) is given as $\boldsymbol{{Q}}\!\left( {{\boldsymbol{{u}}_i}} \right) = {\left[ {Q\!\left( {{u_{i1}}} \right),Q\!\left( {{u_{i2}}} \right), \ldots ,Q\!\left( {{u_{i8}}} \right)} \right]^{\rm{T}}}$ . In our work, $Q\!\left( u \right)$ denotes the hysteretic quantiser, as shown in Equation (5). The subscript i and j are omitted for brevity.

(5) \begin{align}Q(u(t)) = \left\{ {\begin{array}{l@{\quad}l}{{u_p}{\rm{sign}}(u),} & {\dfrac{{{u_p}}}{{1 + \delta }} \lt |u| \le {u_p},\dot u \lt 0,{\rm{or}}}\\[10pt]{} & {{u_i} \lt |u| \le \dfrac{{{u_p}}}{{1 - \delta }},\dot u \gt 0}\\[10pt] {{u_p}(1 + \delta ){\rm{sign}}(u),} & {{u_i} \lt |u| \le \dfrac{{{u_p}}}{{1 - \delta }},\dot u \lt 0,{\rm{or}}}\\[10pt] {} & {\dfrac{{{u_i}}}{{1 - \delta }} \lt |u| \le \dfrac{{{u_i}(1 + \delta )}}{{(1 - \delta )}},\dot u \gt 0}\\[10pt] {0,} & {0 \le |u| \lt \dfrac{{{u_{\min }}}}{{1 + \delta }},\dot u \lt 0,{\rm{or}}}\\[10pt] {} & {\dfrac{{{u_{\min }}}}{{1 + \delta }} \le |u| \le {u_{\min }},\dot u \gt 0,}\\[14pt] {Q\!\left( {u\left( {{t^ - }} \right)} \right)} & {\dot u = 0}\end{array}} \right.\end{align}

where ${u_p} = {\rho ^{1 - p}}{u_{\min }}(p = 1,2, \ldots )$ , $\delta = (1 - \rho )/(1 + \rho )$ , ${u_{\min }} \gt 0$ , $0 \lt \rho \lt 1$ , and $Q(u(t)) \in \left\{ {0, \pm {u_p}, \pm {u_p}(1 + \delta ),p = 1,2, \ldots } \right\}$ . Besides, ${u_{\min }} \gt 0$ denotes the range of the dead-zone for $Q(u(t))$ , and $\rho \gt 0$ represents a measure of quantisation density. We decompose the quantiser Equation (5) into $Q\!\left( {{u_{ij}}} \right) = {\chi _{qij}}\!\left( {{u_{ij}}} \right){u_{ij}} + {d_{qij}},j = 1,2, \cdots ,8$ , where

\begin{align*} \begin{array}{l@{\quad}l}{{\chi _{qij}}\!\left( {{u_{ij}}} \right) = \left\{ {\begin{array}{l@{\quad}l}{\dfrac{{Q\!\left( {{u_{ij}}} \right)}}{{{u_{ij}}}},} & {Q\!\left( {{u_{ij}}} \right) \ne 0,}\\[10pt] {1,} & {Q\!\left( {{u_{ij}}} \right) = 0,}\end{array}} \right.}\\[22pt] {{d_{qij}}\!\left( {{u_{ij}}} \right) = \left\{ {\begin{array}{l@{\quad}l}{0,} & {q\!\left( {{u_{ij}}} \right) \ne 0,}\\[5pt] { - {u_{ij}},} & {Q\!\left( {{u_{ij}}} \right) = 0.}\end{array}} \right.}\end{array}\end{align*}

It can be easily proven that the control coefficient ${\chi _{qij}}$ and the disturbance-like term ${d_{qij}}$ satisfy $1 - \delta \le {\chi _{qij}} \le 1 + \delta ,\left| {{d_{qij}}} \right| \le {u_{\min }}$ . Thus, we can rewritten the quantiser as

(6) \begin{align}\boldsymbol{{Q}}\!\left( {{\boldsymbol{{u}}_i}} \right) = {\boldsymbol\chi _q}\!\left( {{\boldsymbol{{u}}_i}} \right){\boldsymbol{{u}}_i} + {\boldsymbol{{d}}_q}\end{align}

where ${\boldsymbol\chi _q}\!\left( {{\boldsymbol{{u}}_i}} \right) = {\rm{diag}}{\left[ {{\chi _{qi1}}\!\left( {{u_{i1}}} \right)\!,{\chi _{qi2}}\!\left( {{u_{i2}}} \right)\!, \cdots ,{\chi _{qi8}}\!\left( {{u_{i8}}} \right)} \right]^{\rm{T}}}$ , ${\boldsymbol{{d}}_q}\!\left( {{\boldsymbol{{u}}_i}} \right) = {\rm{diag}}\big[ {{d_{qi1}}\!\left( {{u_{i1}}} \right)}, {d_{qi2}}\!\left( {{u_{i2}}} \right), \cdots , {d_{qi8}} \left( {{u_{i8}}} \right)\big]^{\rm{T}}$ .

Define $[\underline M ,\overline M ]$ as the magnitude constraint of the actuator. The saturation nonlinearity in Equation (4) is given as

(7) \begin{align}\boldsymbol{{sat}}\!\left( \boldsymbol{{x}} \right) = {\boldsymbol\chi _s}\!\left( \boldsymbol{{x}} \right)\boldsymbol{{x}} = {\left[ {sat({x_1}),sat({x_2}), \cdots ,sat({x_k})} \right]^{\rm{T}}}\end{align}

where ${\boldsymbol\chi _s}\!\left( \boldsymbol{{x}} \right) = {\rm{diag}}{\left[ {{\chi _{s1}}\!\left( {{x_1}} \right)\!,{\chi _{s2}}\!\left( {{x_2}} \right)\!, \cdots ,{\chi _{sk}}\!\left( {{x_k}} \right)} \right]^{\rm{T}}}$ with

\begin{align*}{\chi _{sk}}\!\left( x \right) = \left\{ {\begin{array}{l@{\quad}l}{\overline{M} /x,} & {{\rm{if}}\;x \geq {{\overline{M}}_u}}\\[5pt] {1,} & {{\rm{if}}\;{{\underline{M}}_u} \le x \le {{\overline{M}}_u}}\\[5pt] {{{\underline{M}}_u}/x,} & {{\rm{if}}\;x \leq {{\underline{M}}_u}}\end{array}} \right.\end{align*}

such that $0 \lt {\chi _{sk}} \lt 1$ . We can rewritten the term $\boldsymbol{{sat}}\!\left[ {\boldsymbol{{Q}}\!\left( {{\boldsymbol{{u}}_i}} \right)} \right]$ in Equation (4) as

(8) \begin{align}\boldsymbol{{sat}}\!\left[ {\boldsymbol{{Q}}\!\left( {{\boldsymbol{{u}}_i}} \right)} \right] = \boldsymbol\chi\!\left( {{\boldsymbol{{u}}_i}} \right){\boldsymbol{{u}}_i} + {\boldsymbol{{d}}_{sq}}\end{align}

where $\boldsymbol\chi\!\left( {{\boldsymbol{{u}}_i}} \right) = {\boldsymbol\chi _s}\!\left( {{\boldsymbol{{u}}_i}} \right){\boldsymbol\chi _q}\!\left( {{\boldsymbol{{u}}_i}} \right),{\boldsymbol{{d}}_{sq}} = {\boldsymbol\chi _s}{\boldsymbol{{d}}_q}$ .

2.3 Graph theory

We use an undirected graph $\boldsymbol{{G}} = \left\{ {\upsilon ,\boldsymbol\varsigma ,\boldsymbol{{A}}} \right\}$ to describe the communication topology of a formation with $n$ spacecrafts, where the node set $\upsilon = \left\{ {{\upsilon _1},{\upsilon _2}, \cdots ,{\upsilon _n}} \right\}$ , the edge set $\varsigma \subseteq \upsilon \times \upsilon $ , and the adjacency matrix $\boldsymbol{{A}} = \left[ {{a_{ij}}} \right] \in {\mathbb{R}^{n \times n}}(i = 1,2, \cdots ,n;\;j = 1,2, \cdots ,n)$ . If node ${\upsilon _i}$ can directly obtain the information of node ${\upsilon _j}$ , there is an edge in the graph from ${\upsilon _j}$ points to ${\upsilon _i}$ , denoted as $({\upsilon _i},{\upsilon _j}) \in \boldsymbol\varsigma $ , and ${\upsilon _j}$ is the child node of ${\upsilon _i}$ . For an undirected graph, if an edge connects ${\upsilon _i}$ and ${\upsilon _j}$ , then the two are parent-child nodes, that is, if $({\upsilon _i},{\upsilon _j}) \in \boldsymbol\varsigma $ , then $({\upsilon _j},{\upsilon _i}) \in \boldsymbol\varsigma $ . In the adjacency matrix $\boldsymbol{{A}}$ , ${a_{ij}} \gt 0$ if $({\upsilon _i},{\upsilon _j}) \in \boldsymbol\varsigma $ ; otherwise, ${a_{ij}} = 0$ . Generally, it is assumed that ${a_{ii}} = 0$ . The Laplace matrix $\boldsymbol{{L}}$ of graph $\boldsymbol{{G}}$ is defined as $\boldsymbol{{L}} = {\rm{diag}}\!\left\{ {\sum\limits_{j = 1}^n {a_{1j}}, \cdots ,\sum\limits_{j = 1}^n {a_{nj}}} \right\} - \boldsymbol{{A}}$ . For undirected graphs, the Laplace matrix $\boldsymbol{{L}}$ is symmetric.

Additionally, denote the virtual leader as spacecraft 0, of which states are given by ${\boldsymbol{{x}}_0}$ and ${\boldsymbol{{v}}_0}$ . A leader-following graph $\overline{\boldsymbol{{G}}}$ contains the virtual leader and the original graph $\boldsymbol{{G}}$ of n follower spacecraft. ${b_i} = 1$ if the ith follower spacecraft can access the virtual leader spacecraft directly, and ${b_i} = 0$ , otherwise. The graph $\overline{\boldsymbol{{G}}}$ is connected if there exists a path in $\overline{\boldsymbol{{G}}}$ from the leader node ${\upsilon _0}$ to every node ${\upsilon _i}$ . If $\overline{\boldsymbol{{G}}} $ is connected, the matrix $\boldsymbol{{L}} + \boldsymbol{{B}}$ associated with $\overline{\boldsymbol{{G}}}$ is symmetric and positive definite [Reference Hong, Hu and Gao35], where $\boldsymbol{{B}} = {\rm{diag}}({b_1}, \cdots ,{b_2})$ .

2.4 Problem statement

This paper committed to developing a distributed control scheme such that all follower spacecraft can track the virtual leader with the desired formation shape considering limited communication, i.e., $\mathop {{\rm{lim}}}\limits_{t \to T} \left\| {{\boldsymbol{{x}}_i} - {\boldsymbol\Delta _{i,0}} - {\boldsymbol{{x}}_0}} \right\| \le {o_1},\mathop {{\rm{lim}}}\limits_{t \to T} \left\| {{\boldsymbol{{v}}_i} - {{\dot{\boldsymbol\Delta} }_{i,0}} - {\boldsymbol{{v}}_0}} \right\| \le {o_2}$ , where ${\boldsymbol\Delta _{i,0}}$ is the ith follower spacecraft’s expected deviation formation vector relative to the virtual leader, and ${o_1},{o_2}$ are compact sets around zero.

Remark 2. We consider the desired position deviation formation vector ${\boldsymbol\Delta _{i,0}}$ to be an internal parameter of the ith follower spacecraft, which can be obtained without external communication.

In order to get the main results of this paper, the following assumptions are adopted.

Assumption 1. The desired states ${\boldsymbol{{x}}_0}$ and its first two derivatives are uniformly bounded such that $\left\| {{\boldsymbol{{x}}_0}} \right\| \le {B_1}$ , $\left\| {{{\dot{\boldsymbol{{x}}}}_0}} \right\| \le {B_2}$ , $\left\| {{{\ddot{\boldsymbol{{x}}}}_0}} \right\| \le {B_3}$ , where ${B_1}$ , ${B_2}$ , ${B_3}$ are positive constants.

Assumption 2. The disturbance ${\boldsymbol{{d}}_i}$ is bounded and satisfies $\left\| {{\boldsymbol{{d}}_i}} \right\| \le {d_{max}}$ , where ${d_{max}}$ is a positive constant.

Assumption 3. [Reference Zhou, Zhou, Shi, Li and Kan36] For a constant ${\chi _{min}}$ , we consider the control input restricted in a region such that ${\chi _{min}} \le {\chi _{sk}} \lt 1$ .

Assumption 4. [Reference Kun, Mou, Qingxian and Ronggang37, Reference Shi, Zhou, Zhou, Li and Chen38] For the practical SFF system described by Equation (4) with input saturation, ${\boldsymbol{{g}}_i}({\boldsymbol{{x}}_i})$ should be bounded and satisfies $\parallel {\boldsymbol{{g}}_i}({\boldsymbol{{x}}_i})\parallel \le {g_{max}}$ .

3.0 Distributed control scheme design

The procedure for designing the distributed control scheme is presented in this section. Firstly, a distributed event-triggered observer is designed with continuous communications, based on which we proposed an event-based distributed observer without continuous communications. Further, a distributed fixed-time controller is proposed for the follower spacecraft in the presence of the actuator saturation and the input quantisation problems. The block diagram of the distributed control algorithm is given in Fig. 2.

Figure 2. Schematic of the distributed control algorithm.

Remark 3. For the following reasons, we only add ETM in the observer, even though both the observer and the controller occupy communication resources. Firstly, the communication resources saved by ETM come at the cost of system performance; secondly, since the observer converges much faster than the controller, embedding the event-triggered mechanism in the observer has minimal impact on the system performance. If ETM is also embedded in the controller, it will save communication resources, but the system performance will be significantly reduced.

3.1 Event-triggered fixed-time observer design

This subsection develops the event-triggered fixed-time observer with continuous and without continuous communications to reconstruct the virtual leader spacecraft’s information.

Define a auxiliary variable as ${\boldsymbol\zeta _i} = \sum\nolimits_{j = 1}^n {a_{ij}}\!\left( {{{\hat{\boldsymbol{{v}}}}_{i,0}} - {{\hat{\boldsymbol{{v}}}}_{j,0}}} \right) + {b_i}\!\left( {{{\hat{\boldsymbol{{v}}}}_{i,0}} - {\boldsymbol{{v}}_0}} \right)$ where ${\hat{\boldsymbol{{v}}}_{i,0}}$ denotes the estimation of ${\boldsymbol{{v}}_0}$ of the ith follower spacecraft, ${\widetilde{\boldsymbol{{v}}}_i} = {\hat{\boldsymbol{{v}}}_{i,0}} - {\boldsymbol{{v}}_0}$ is the estimation error, $\widetilde{\boldsymbol{{v}}} = {\left[ {{{\widetilde{\boldsymbol{{v}}}}_1}^{\rm{T}},{{\widetilde{\boldsymbol{{v}}}}_2}^{\rm{T}}, \ldots ,{{\widetilde{\boldsymbol{{v}}}}_n}^{\rm{T}}} \right]^{\rm{T}}}$ , $\boldsymbol\zeta = \boldsymbol{{H}}\widetilde{\boldsymbol{{v}}} = {\left[ {{\boldsymbol\zeta _1}^{\!\!\!\rm{T}},{\boldsymbol\zeta _2}^{\!\!\!\rm{T}}, \ldots ,{\boldsymbol\zeta _n}^{\!\!\!\rm{T}}} \right]^{\rm{T}}}$ .

3.1.1 Event-triggered fixed-time observer with continuous communications

The event-triggered fixed-time observer with continuous communications for the ith follower spacecraft is described as:

(9) \begin{align}{{{\dot{\hat{\boldsymbol{{v}}}}}_{i,0}} = {{\dot{\hat{\boldsymbol{{v}}}}}_{i,0}}({t_k})} &= { - {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right] - {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right]}\nonumber \\[5pt] & \quad { - {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right] - {o_4}{\rm{obsat}}\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right]}\end{align}
(10) \begin{align}{t_{k + 1}} = \inf\!\left\{ {t \gt {t_k}\;:\;{h_{ij}} \gt 0,\forall j,j = 1, \cdots ,8} \right\}\end{align}

where

(11) \begin{align}\begin{array}{l}{{h_{ij}}\!\left( t \right) = } {\left| {{E_{ij}}\!\left( t \right)} \right| - \iota \left| {{v_{uij}}} \right| - {k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \left( {{k_{o3}}\left| {{\zeta _{ij}}(t)} \right|} \right)}}} }\end{array}\end{align}

With $0 \lt \alpha \lt 1$ , $\beta \gt 1$ , $\iota \in (0,1),{o_1},{o_2},{o_3},{o_4},{k_{o1}},{k_{o2}},{k_{o3}}$ are positive observer parameters, which will be chosen later. ${\boldsymbol{{E}}_i}\!\left( t \right) = {\left[ {{E_{i1}}, \cdots ,{E_{ij}}, \cdots } \right]^{\rm{T}}} = {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_4}{\rm{obsat}}\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] - {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_4}{\rm{obsat}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] = - {\dot{\hat{\boldsymbol{{v}}}}_{i,0}}({t_k}) - {\boldsymbol{{v}}_{ui}}(t)$ where ${\boldsymbol{{v}}_{ui}}(t) = {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] + {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] + {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] + {o_4}{\rm{obsat}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right]$ . ${\rm{obsat}}(\boldsymbol{{x}}) = [{\rm{obsat}}({x_1}), \cdots , {\rm{obsat}}({x_n}{)]^{\rm{T}}}$ where ${\rm{obsat}}(x)$ is designed as:

(12) \begin{align}{\rm{obsat}}(x) = \left\{ {\begin{array}{l@{\quad}l}{x/{\delta _o},} & {|x| \le {\delta _o}}\\[5pt] {{\rm{sign}}(x),} & {|x| \gt {\delta _o}}\end{array}} \right.\end{align}

Theorem 1. Consider the system Equation (4) with the designed observer Equation (9) and the event-triggered updated rule Equation (10). Suppose Assumption 1 holds. ${\hat{\boldsymbol{{v}}}_i}$ will estimate ${\boldsymbol{{v}}_0}$ in fixed time ${T_1}$ if the parameters satisfy $(1 - \iota ){o_4} \gt {B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} $ . Meanwhile, there always exists a non-zero bound time for any two triggers, which means ${t_{k + 1}} - {t_k} \gt \tau \gt 0$ , and $\tau $ is a positive constant.

Proof. Consider the following Lyapunov function candidate:

(13) \begin{align}{V_\zeta } = \frac{1}{2}{\widetilde{\boldsymbol{{v}}}^{\rm{T}}}\boldsymbol{{H}}\widetilde{\boldsymbol{{v}}}\end{align}

where $\boldsymbol{{H}} = (\boldsymbol{{L}} + \boldsymbol{{B}}) \otimes {\boldsymbol{{I}}_8}$ describes the communication topology of the formation. The time derivative of Equation (13) along with Equation (9) is given by

(14) \begin{align}{{{\dot V}_{\boldsymbol\zeta} }} &= {{\boldsymbol\zeta ^{\rm{T}}}\dot{\tilde{\boldsymbol{{v}}}}(t) = {\boldsymbol\zeta ^{\rm{T}}}\!\left[ {\dot{\hat{\boldsymbol{{v}}}}(t) - {{\dot{\boldsymbol{{v}}}}_{n0}}} \right]}\nonumber \\[8pt] &= {- {\boldsymbol\zeta ^{\rm{T}}}\!\left[ {\boldsymbol{{E}} + {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\boldsymbol\zeta ) + {o_2}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\zeta ) + {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left( \boldsymbol\zeta \right) + {o_4}{\rm{obsat}}\boldsymbol\zeta ) + {{\dot{\boldsymbol{{v}}}}_{n0}}} \right]}\nonumber \\[8pt] & \quad \le { |\boldsymbol\zeta {|^{\rm{T}}}|\boldsymbol{{E}}| - {o_1}{\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\boldsymbol\zeta ) - {o_2}{\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\zeta )}\nonumber \\[8pt] & \quad { - {o_3}{\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }\!\left( \boldsymbol\zeta \right) - {o_4}{\boldsymbol\zeta ^{\rm{T}}}{\rm{obsat}}(\boldsymbol\zeta ) - {\boldsymbol\zeta ^{\rm{T}}}{{\dot{\boldsymbol{{v}}}}_{n0}}}\end{align}

where $E = {\left[ {{E_1},{E_2}, \cdots ,{E_n}} \right]^{\rm{T}}}$ , ${v_{n0}} = {\left[ {{v_0},{v_0}, \cdots ,{v_0}} \right]^{\rm{T}}} \in {\mathbb{R}^{8n \times 1}}$ .

Considering the event-triggered update rule Equation (11) and Assumption 1 , we can imply the following relationship from Equation (14):

(15) \begin{align}{{{\dot V}_{\boldsymbol\zeta} }} & \le { \boldsymbol\zeta {|^{\rm{T}}}|\boldsymbol{{E}}| - {o_1}{\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\boldsymbol\zeta ) - {o_2}{\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\zeta )}\nonumber \\[5pt]& { - {o_3}{\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\zeta ) - {o_4}{\boldsymbol\zeta ^{\rm{T}}}{\rm{obsat}}(\boldsymbol\zeta ) + {B_3}\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\boldsymbol\zeta _{ij}}} \right|}\nonumber\\[5pt]& \le { - {o_1}(1 - \iota ){\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\boldsymbol\zeta ) - {o_2}(1 - \iota ){\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\zeta ) - {o_3}(1 - \iota ){\boldsymbol\zeta ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\zeta )}\nonumber\\[5pt]& { - {o_4}(1 - \iota )\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 {\zeta _{ij}}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right) + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} } \right)\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\zeta _{ij}}} \right|}\nonumber\\[5pt]& \le { - {o_1}(1 - \iota )(8n{)^{\frac{{\alpha - 1}}{{2\alpha }}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \frac{1}{\alpha }}} - {o_2}(1 - \iota )(8n{)^{\frac{{1 - \beta }}{2}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \beta }} - {o_3}(1 - \iota )\parallel \boldsymbol\zeta {\parallel ^{1 + \alpha }}}\nonumber\\[5pt]& { - {o_4}(1 - \iota )\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 {\zeta _{ij}}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right) + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} } \right)\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\zeta _{ij}}} \right|}\end{align}

where $ - {\boldsymbol\zeta ^T}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\zeta ) = - \sum\nolimits_{i = 1}^n \sum\nolimits_{j = 1}^8 {\left( {\zeta _{ij}^2} \right)^{\frac{{1 + \alpha }}{2}}} \le - \parallel \boldsymbol\zeta {\parallel ^{1 + \alpha }}$ and $ - {\boldsymbol\zeta ^T}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\zeta ) = - \sum\nolimits_{i = 1}^n \sum\nolimits_{j = 1}^8 {\left( {\zeta _{ij}^2} \right)^{\frac{{1 + \beta }}{2}}} \le - {(8n)^{\frac{{1 - \beta }}{2}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \beta }}$ are derived by Lemma 6 and Lemma 7 .

Noting that $\sum\nolimits_{i = 1}^n \sum\nolimits_{j = 1}^8 {\zeta _{ij}}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right) = \frac{1}{{{\delta _o}}}\sum \left( {{{\bar \zeta }_{ij}}^2} \right) + \sum \left| {{{\bar \zeta }_{ij}}} \right|$ where ${\bar \zeta _{ij}} = \left\{ {{\zeta _{ij}}:\left| {{\zeta _{ij}}} \right| \le {\delta _o}} \right\},{\bar \zeta _{ij}} = \left\{ {{\zeta _{ij}}:\left| {{\zeta _{ij}}} \right| \gt {\delta _o}} \right\}$ , $k \in [0,8n]$ .

For Equation (15), one has

(16) \begin{align}&\quad { - {o_4}(1 - \iota )\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 {\zeta _{ij}}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right) + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{02}}}}} } \right)\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\zeta _{ij}}} \right|}\nonumber \\[5pt] &\le { - {o_4}(1 - \iota )\left[ {\sum \left| {{\zeta _{ij}}} \right| + \frac{1}{{{\delta _o}k}}{{\left( {\sum \left| {{\zeta _{ij}}} \right|} \right)}^2}} \right]}\nonumber \\[5pt] & \quad + { \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{02}}}}} } \right)\sum \left| {{\zeta _{ij}}} \right| + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} } \right)\sum \left| {{\zeta _{ij}}} \right|}\nonumber \\[5pt] & \le { - {o_4}(1 - \iota )\frac{1}{{{\delta _o}k}}{{\left( {\sum \left| {{\zeta _{ij}}} \right|} \right)}^2} + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} } \right)\sum \left| {{\zeta _{ij}}} \right|}\nonumber {}\\[5pt] & = { - {o_4}(1 - \iota )\sum \left| {{\zeta _{ij}}} \right| + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{02}}}}} } \right)\sum \left| {{\zeta _{ij}}} \right|}\nonumber \\[5pt] & \quad { - {o_4}(1 - \iota )\frac{1}{{{\delta _o}k}}{{\left( {\sum \left| {{\zeta _{ij}}} \right|} \right)}^2} + {o_4}(1 - \iota )\sum \left| {{\zeta _{ij}}} \right|}\nonumber \\[5pt] & \le {\frac{{{o_4}(1 - \iota ){\delta _o}k}}{4}}\end{align}

when $k \ne 0$ , where the inequalities $k\sum \left( {{{\bar \zeta }_{ij}}^2} \right) \geqslant {\left( {\sum \left| {{{\bar \zeta }_{ij}}} \right|} \right)^2}$ which are derived by Lemma 6 has been used. It should be noticed that $ - {o_4}(1 - \iota )\sum\nolimits_{i = 1}^n \sum\nolimits_{j = 1}^8 {\zeta _{ij}}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right) + \left( {{B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{02}}}}} } \right)\sum\nolimits_{i = 1}^n \sum\nolimits_{j = 1}^8 \left| {{\zeta _{ij}}} \right| \le 0$ when $k = 0$ . It seems that the discussion of whether $k = 0$ in Equation (16) is redundant from the results, but the application of Lemma 6 introduces a singular term $\frac{1}{k}$ . Therefore, the discussion above is necessary. We can rewrite Equation (15) as

(17) \begin{align}{{{\dot V}_{\boldsymbol\zeta}} } &\le { - {o_1}(1 - \iota )(8n{)^{\frac{{\alpha - 1}}{{2\alpha }}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \frac{1}{\alpha }}} - {o_2}(1 - \iota )(8n{)^{\frac{{1 - \beta }}{2}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \beta }}}\nonumber \\[5pt] & \quad { - {o_3}(1 - \iota )\parallel \boldsymbol\zeta {\parallel ^{1 + \alpha }} + \frac{{{o_4}(1 - \iota ){\delta _o}k}}{4}}\end{align}

Noting that ${V_\zeta } = \frac{1}{2}{\tilde{\boldsymbol{{v}}}^{\rm{T}}}H\tilde{\boldsymbol{{v}}} = \frac{1}{2}{\boldsymbol\zeta ^{\rm{T}}}{H^{ - 1}}\boldsymbol\zeta \le \frac{1}{{2{\lambda _{{\rm{min}}}}\!\left( H \right)}}\parallel \boldsymbol\zeta {\parallel ^2}$ , Equation (17) becomes

(18) \begin{align}{{{\dot V}_{\boldsymbol\zeta} } } &\le { - {o_2}(1 - \iota )(8n{)^{\frac{{1 - \beta }}{2}}}{{\left[ {2{\lambda _{\min }}(\boldsymbol{{H}})} \right]}^{\frac{{1 + \beta }}{2}}}V_{\boldsymbol\zeta}^{\frac{{1 + \beta }}{2}}}\nonumber \\[5pt] &\quad { - {o_3}(1 - \iota ){{\left[ {2{\lambda _{\min }}(\boldsymbol{{H}})} \right]}^{\frac{{1 + \alpha }}{2}}}V_{\boldsymbol\zeta}^{\frac{{1 + \alpha }}{2}} + \frac{{{o_4}(1 - \iota ){\boldsymbol\delta _o}k}}{4}}\end{align}

which implies that ${V_{\boldsymbol\zeta} }$ converges to a small neighbourhood of zero in fixed time ${T_1} \le 1/\left[ {{\kappa _1}\!\left( {1 - \frac{{1 + \alpha }}{2}} \right)} \right] + 1/\left[ {{\kappa _2}\!\left( {\frac{{1 + \beta }}{2} - 1} \right)} \right]$ where ${\kappa _1} = {o_3}(1 - \iota ){\left[ {2{\lambda _{\min }}(\boldsymbol{{H}})} \right]^{\frac{{1 + \alpha }}{2}}}$ , ${\kappa _2} = {o_2}(1 - \iota )(8n{)^{\frac{{1 - \beta }}{2}}}{\left[ {2{\lambda _{\min }}(\boldsymbol{{H}})} \right]^{\frac{{1 + \beta }}{2}}}$ . Thus, we conclude that the estimation error $\widetilde {{\boldsymbol{{v}}_i}} = {\hat{\boldsymbol{{v}}}_i} - {\boldsymbol{{v}}_0}$ converges to a small neighbourhood of zero in fixed time.

For ${\boldsymbol{{E}}_i}\!\left( t \right) = {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] + {o_4}{\rm{sign}}\!\left[ {{\boldsymbol\zeta _i}\!\left( {{t_k}} \right)} \right] - {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right] - {o_4}{\rm{sign}}\!\left[ {{\boldsymbol\zeta _i}\!\left( t \right)} \right]$ , we have

(19) \begin{align}{\left| {{{\dot E}_{ij}}(t)} \right|} &= \left| \left[ {{o_1}{\rm{si}}{\rm{g}}^{\frac{1}{\alpha }}\!\left(\zeta _{ij}\right) + {o_2}{\rm{si}}{\rm{g}^{\beta}}\!\left(\zeta _{ij} \right) + {o_3}{\rm{si}}{{\rm{g}}^{\alpha }}\!\left( \zeta _{ij} \right) + {o_4}{\rm{obsat}}\!\left( {{\zeta _{ij}}} \right)} \right]' \right|\nonumber \\[5pt] & \le { \left| {{o_1}\frac{1}{\alpha }{\rm{si}}{{\rm{g}}^{\frac{{1 - \alpha }}{\alpha }}}\!\left( {{\zeta _{ij}}} \right) + {o_2}\beta {\rm{si}}{{\rm{g}}^{\beta - 1}}\!\left( {{\zeta _{ij}}} \right) + {o_3}\alpha {\rm{si}}{{\rm{g}}^{\alpha - 1}}\!\left( {{\zeta _{ij}}} \right) + \frac{1}{{{\delta _o}}}} \right|\left| {{{\dot \zeta }_{ij}}} \right|}\end{align}

From Equation (18) we know that ${\zeta _{ij}}$ is bounded, then we define the upper bound of $\left| {{o_1}\frac{1}{\alpha }{\rm{si}}{{\rm{g}}^{\frac{{1 - \alpha }}{\alpha }}}\!\left( {{\zeta _{ij}}} \right) + {o_2}\beta {\rm{si}}{{\rm{g}}^{\beta - 1}}\!\left( {{\zeta _{ij}}} \right) + {o_3}\alpha {\rm{si}}{{\rm{g}}^{\alpha - 1}}\!\left( {{\zeta _{ij}}} \right) + \frac{1}{{{\delta _o}}}} \right|$ as ${\eta _1}$ . Equation (19) can be transformed into:

(20) \begin{align}{\left| {{{\dot E}_{ij}}(t)} \right|} &\le { {\eta _1}\left| {{{\dot \zeta }_{ij}}} \right| \le {\eta _1}{h_{\max }}\!\left[ {{{\dot{\hat{v}}}_{ij}}(t) - {{\dot v}_{0ij}}\!\left( {{t_k}} \right)} \right]}\nonumber \\[5pt] & \le { {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|}\end{align}

where ${h_{max}}$ is the largest element in matrix $\boldsymbol{{H}}$ , ${\eta _{2ij}}\!\left( {{t_k}} \right) = {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right] + {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left[ {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right] + {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left[ {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right] + {o_4}{\rm{sign}}\!\left[ {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right]$ . Since ${\dot E_{ij}}({t_k}) = 0$ , it follows when $t \gt {t_k}$ :

(21) \begin{align}\left| {{E_{ij}}\!\!\left( t \right)} \right| \le \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s\end{align}

According to the triggering condition Equation (11), one has that the next event of state j of agent i will not be triggered before ${h_{ij}}(t) = 0$ or equivalently $\left| {{E_{ij}}\!\!\left( t \right)} \right| = \iota {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left( {\left| {{\zeta _{ij}}\!\left( t \right)} \right|} \right) + \iota {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left( {\left| {{\zeta _{ij}}\!\left( t \right)} \right|} \right) + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {\left| {{\zeta _{ij}}(t)} \right|} \right)}}} $ . From Equation (21), one has:

(22) \begin{align}{\left| {{E_{ij}}\!\!\left( {{t_{k + 1}}} \right)} \right|} &= { \iota {o_2}{\rm{si}}{{\rm{g}}^\beta }\!\left( {\left| {{\zeta _{ij}}\!\left( {{t_{k + 1}}} \right)} \right|} \right) + \iota {o_3}{\rm{si}}{{\rm{g}}^\alpha }\!\left( {\left| {{\zeta _{ij}}\!\left( {{t_{k + 1}}} \right)} \right|} \right)}\nonumber \\[5pt] & \quad + { {k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {\!\left( {{\xi _{ij}}\!\left( {{t_{k + 1}}} \right)|} \right)} \right.}}} }\nonumber \\[5pt] & \le { \int_{{t_k}}^{{t_{k + 1}}} {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s}\nonumber \\[5pt] & \le { \int_{{t_k}}^{{t_{k + 1}}} {\eta _1}{h_{\max }}\left| {{\eta _{2ij\max }} + {B_3}} \right|{\rm{d}}s}\end{align}

where ${\eta _{2ij\max }}$ is the upper bound of ${\eta _{2ij}}\!\left( {{t_k}} \right)$ . Equation (22) yields that ${t_{k + 1}} - {t_k} \geqslant {k_{o1}}\sqrt e / \!\left( {{\eta _1}{h_{\max }}\left| {{\eta _{2ij\max }} + {B_3}} \right|} \right)$ , which is shown that Zeno-free is guaranteed.

3.1.2 Event-triggered fixed-time observer without continuous communications

This subsection proposes the event-triggered fixed-time observer without continuous communications based on subsection 3.1.1.

For the ith follower spacecraft, $\dot{\hat{v}}_{i}$ is the same as Equation (9). The triggering condition is described as:

(23) \begin{align}{t_{k + 1}} = \inf\!\left\{ {t \gt {t_k}\;:\;{g_{ij}} \gt 0,\forall j,j = 1, \cdots ,8} \right\}\end{align}

where

(24) \begin{align}{{g_{ij}}(t)} & = { \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s}\nonumber \\[5pt] & { - \frac{\iota }{{1 + \iota }}\left| { - {{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + \frac{{{k_{o1}}}}{\iota }\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {{k_{o3}}\left| {{\xi _{ij}}\!\left( {{t_k}} \right)} \right|} \right)}}} } \right|}\end{align}

Theorem 2. Consider the system Equation (4) with the designed observer Equation (9) and the event-triggered updated rule Equation (24). Suppose Assumption 1 holds. ${\hat{\boldsymbol{{v}}}_i}$ will estimate ${\boldsymbol{{v}}_0}$ in fixed time ${T_1}$ if the parameters satisfy $(1 - \iota ){o_4} \gt {B_3} + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} $ .. Meanwhile, there always exists a non-zero bound time for any two triggers, which means ${t_{k + 1}} - {t_k} \gt \tau \gt 0$ , and $\tau $ is a positive constant.

Proof. According to Equation (21), we can obtain $\left| {{E_{ij}}\!\!\left( t \right)} \right| \le \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|ds$ , $t \in \left[ {{t_k},{t_{k + 1}}} \right)$ . Further, the triggering function Equation (24) enforces

(25) \begin{align}\left| {{E_{ij}}\!\!\left( t \right)} \right| \le \frac{\iota }{{1 + \iota }}\left| { - {{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + \frac{{{k_{o1}}}}{\iota }\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {\left| {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right|} \right)}}} } \right|\end{align}

Noting that ${E_{ij}} = - {\dot{\hat{v}}_{ij}}\!\left( {{t_k}} \right) - {v_{uij}}$ , Equation (25) can be transformed into

(26) \begin{align}\left| {{E_{ij}}(t)} \right| + \iota \left| {{{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + {v_{ui}}} \right| \le \iota \left| {{{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right)} \right| + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {\left| {{\zeta _{ij}}\!\left( {{t_k}} \right)} \right|} \right)}}} \end{align}

which is is a sufficient condition for

(27) \begin{align}\begin{array}{l}{} {\left| {{E_{ij}}\!\!\left( t \right)} \right| \le \iota \left| {{v_{uij}}} \right| + {k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {\left| {{\zeta _{ij}}(t)} \right|} \right)}}} }\end{array}\end{align}

Equation (27) is the same as the event-triggered updated rule Equation (11). Hence, similar to the proof of Theorem 1 , the estimation error $\widetilde {{\boldsymbol{{v}}_i}} = {\hat{\boldsymbol{{v}}}_i} - {\boldsymbol{{v}}_0}$ converges to a small neighbourhood of zero in fixed time ${T_1} \le 1/\left[ {{\kappa _1}\!\left( {1 - \frac{{1 + \alpha }}{2}} \right)} \right] + 1/\left[ {{\kappa _2}\!\left( {\frac{{1 + \beta }}{2} - 1} \right)} \right]$ , and Zeno-free is guaranteed.

Remark 4. Compared with the distributed fixed-time observer in [Reference Han, Yuanqing, Zhang and Zhang12], an additional term ${\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right]$ is adopted in Equation (9) to handle certain nonlinear term generated by the fixed-time controller. As shown in Equation (43), we guarantee the fixed time stability of the entire closed-loop system without recourse to the separation principle. Besides, the observer introduces an ETM without continuous communications, reducing energy consumption significantly.

Remark 5. The sign function results in adverse chattering. Differing from the work of [Reference Han, Yuanqing, Zhang and Zhang12], the saturation function is used in place of the sign function for the alleviation of chattering. The constant term $\frac{{{o_4}(1 - \iota )\delta k}}{4}$ is thus introduced in Equation (18), implying that ${V_{\boldsymbol\zeta} }$ converges to a region. Nonetheless, the region can be as small as desired by selecting parameter $\delta $ . The region is correspondingly scaled down to a smaller neighbourhood of the origin and the error caused by the approximation is negligible.

Remark 6. The ETM without continuous communications in Equation (24) is inspired by [Reference Liu, Zhang, Yu and Sun17], where the event-triggered mechanism is used in a consensus control scheme. Differently, Equation (24) embeds a bounded nonlinear term $\frac{{{k_{o1}}}}{\iota }\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {{k_{o3}}\left| {{\xi _{ij}}\!\left( {{t_k}} \right)} \right|} \right)}}} $ , which realises a more reasonable trigger threshold with characteristics of the exponential function and the hyperbolic tangent function in the transient and steady state.

Remark 7. In view of Equation (24), the key rules in selecting the parameters are: The trigger interval can be increased by increasing the parameter ${k_{o1}}$ . Subsequently, tuning the parameter ${k_{o2}}$ can adjust the exponential convergence rate of the triggering condition during the adjustment period. Meanwhile, the selection of ${k_{o3}}$ should ensure that ${k_{o3}}\left| {{\xi _{ij}}\!\left( {{t_k}} \right)} \right|$ can cover the unsaturated area of the hyperbolic tangent function.

3.2 Fixed-time controller design

In this subsection, a fixed-time distributed control scheme is designed for each follower spacecraft, handling the actuator saturation and the input quantisation problems ingeniously without embedding auxiliary systems. Define a auxiliary vector for ith follower spacecraft as

(28) \begin{align}{{{\boldsymbol\sigma}_i}} &= { \sum\limits_{j = 1}^n {a_{ij}}\!\left[ {({\boldsymbol{{x}}_i} - {\boldsymbol\Delta _{i,0}}) - ({\boldsymbol{{x}}_j} - {\boldsymbol\Delta _{j,str}})} \right] + {b_i}\!\left( {{\boldsymbol{{x}}_i} - {\boldsymbol\Delta _{i,0}} - {\boldsymbol{{x}}_0}} \right)}\nonumber \\[5pt] & = { \sum\limits_{j = 1}^n {l_{ij}}{\boldsymbol{{e}}_{xj}} + {b_i}{\boldsymbol{{e}}_{xi}}}\end{align}

where ${\boldsymbol{{e}}_{xi}} = {\boldsymbol{{x}}_i} - {\boldsymbol\Delta _{i,0}} - {\boldsymbol{{x}}_0}$ . Denote $\boldsymbol\sigma = {\left[ {\boldsymbol\sigma _1^{\rm{T}},\boldsymbol\sigma _2^{\rm{T}}, \ldots ,\boldsymbol\sigma _n^{\rm{T}}} \right]^{\rm{T}}}$ , ${\boldsymbol{{e}}_x} = \left[ {\boldsymbol{{e}}_{x1}^{\rm{T}},\boldsymbol{{e}}_{x2}^{\rm{T}},} \right.$ ${\left. { \ldots ,\boldsymbol{{e}}_{xn}^{\rm{T}}} \right]^{\rm{T}}}$ , then we have $\boldsymbol\sigma = \boldsymbol{{H}}{\boldsymbol{{e}}_x}$ . Based on the backstepping technique, the controller is designed as follows:

Step 1. Design of a virtual control scheme for ${\boldsymbol\gamma _i} = {\boldsymbol{{v}}_i} - {\hat{\boldsymbol{{v}}}_{i,0}} - {\dot{\boldsymbol\Delta}_{i,0}} + {k_1}{\rm{si}}{{\rm{g}}^\beta }\!\left( {{\boldsymbol\sigma _i}} \right)$ , where ${k_1}$ is a positive parameter.

Consider ${\boldsymbol\gamma _i}$ as a virtual control input for the system $\dot{\boldsymbol\sigma} = \boldsymbol{{H}}{\dot{\boldsymbol{{e}}}_x} = \boldsymbol{{H}}{\boldsymbol{{e}}_v} = \boldsymbol{{H}}\!\left( {{\boldsymbol{{v}}_i} - {{\hat{\boldsymbol{{v}}}}_{i,0}} - {{\dot{\boldsymbol\Delta}}_{i,0}} + {{\widetilde{\boldsymbol{{v}}}}_i}} \right)$ , A virtual control scheme for ${\boldsymbol\gamma _i}$ is designed as ${\boldsymbol\gamma _{di}} = - {k_2}{\rm{si}}{{\rm{g}}^\alpha }\!\left( {{\boldsymbol\sigma _i}} \right)$ , where ${k_1}$ is a positive parameter, ${\boldsymbol{{e}}_v} = {\left[ {\boldsymbol{{e}}_{v1}^{\rm{T}},\boldsymbol{{e}}_{v2}^{\rm{T}}, \ldots ,\boldsymbol{{e}}_{vn}^{\rm{T}}} \right]^{\rm{T}}}$ , ${\boldsymbol{{e}}_{vi}} = {\boldsymbol{{v}}_i} - {\dot{\boldsymbol\Delta}_{i,0}} - {\boldsymbol{{v}}_0}$ .

Consider the following Lyapunov function candidate

(29) \begin{align}{V_1} = \frac{1}{2}{\boldsymbol\sigma ^{\rm{T}}}{\boldsymbol{{H}}^{ - 1}}\boldsymbol\sigma \end{align}

and the time derivative of ${V_1}$ is given by

(30) \begin{align}{{{\dot V}_1}} &= { {\boldsymbol\sigma ^{\rm{T}}}{\boldsymbol{{e}}_v} = {\boldsymbol\sigma ^{\rm{T}}}\!\left( {\boldsymbol{{v}} - \hat{\boldsymbol{{v}}} - {{\dot{\boldsymbol\Delta}}_{str}} + \tilde{\boldsymbol{{v}}}} \right) = {\boldsymbol\sigma ^T}(\boldsymbol\gamma + \tilde{\boldsymbol{{v}}})}\nonumber \\[5pt] & = { {\boldsymbol\sigma ^{\rm{T}}}\!\left( {\boldsymbol\gamma - {\boldsymbol\gamma _d}} \right) + {\boldsymbol\sigma ^{\rm{T}}}\tilde{\boldsymbol{{v}}} - {k_1}{\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\sigma ) - {k_2}{\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma )}\end{align}

where $\boldsymbol\gamma = {\left[ {{\boldsymbol\gamma _1}^{\rm{T}},{\boldsymbol\gamma _2}^{\rm{T}}, \ldots ,{\boldsymbol\gamma _n}^{\rm{T}}} \right]^{\rm{T}}}$ , ${\boldsymbol\gamma _d} = {\left[ {{\boldsymbol\gamma _{d1}}^{\rm{T}},{\boldsymbol\gamma _{d2}}^{\rm{T}}, \ldots ,{\boldsymbol\gamma _{dn}}^{\rm{T}}} \right]^{\rm{T}}}$ .

According to Lemma 2 and Lemma 5 , we can imply the following unequal relationship:

(31) \begin{align}{{\boldsymbol\sigma ^{\rm{T}}}\!\left( {\boldsymbol\gamma - {\boldsymbol\gamma _d}} \right)} & \le { \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\sigma _{ij}}} \right|\left| {{\gamma _{ij}} - {\gamma _{dij}}} \right| \le \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 {2^{1 - \alpha }}{{\left| {{\sigma _{ij}}\parallel {\kappa _{ij}}} \right|}^\alpha }}\nonumber \\[5pt] & \le { \frac{{{2^{1 - \alpha }}}}{{1 + \alpha }}{\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma ) + \frac{{{2^{1 - \alpha }}\alpha }}{{1 + \alpha }}{\boldsymbol\kappa ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\kappa )}\end{align}
(32) \begin{align}{\boldsymbol\sigma ^{\rm{T}}}\tilde{\boldsymbol{{v}}} \le \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \left| {{\sigma _{ij}}} \right|{\!\left( {{{\left| {{{\tilde v}_{ij}}} \right|}^{\frac{1}{\alpha }}}} \right)^\alpha } \le \frac{1}{{1 + \alpha }}{\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma ) + \frac{\alpha }{{1 + \alpha }}{\tilde{\boldsymbol{{v}}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\tilde{\boldsymbol{{v}}})\end{align}

where ${\kappa _{ij}} = {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{ij}}} \right) - {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{dij}}} \right)$ .

Substitution of Equation (31) and Equation (32) into Equation (30) yields that

(33) \begin{align}{{{\dot V}_1}} & \le { - {k_1}{\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\sigma ) - \!\left( {{k_2} - \frac{{1 + {2^{1 - \alpha }}}}{{1 + \alpha }}} \right){\boldsymbol\sigma ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma )}\nonumber \\[5pt] & + { \frac{{{2^{1 - \alpha }}\alpha }}{{1 + \alpha }}{\boldsymbol\kappa ^{\rm{T}}}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\kappa ) + \frac{\alpha }{{1 + \alpha }}{{\tilde{\boldsymbol{{v}}}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}(\tilde{\boldsymbol{{v}}})}\nonumber \\[5pt] & \le { - {k_1}\parallel \boldsymbol\sigma {\parallel ^{1 + \beta }} - \!\left( {{k_2} - \frac{{1 + {2^{1 - \alpha }}}}{{1 + \alpha }}} \right)\parallel \boldsymbol\sigma {\parallel ^{1 + \alpha }}}\nonumber\\[5pt] & + { \frac{{{2^{1 - \alpha }}\alpha }}{{1 + \alpha }}\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }} + \frac{\alpha }{{1 + \alpha }}\tilde{\boldsymbol{{v}}}{\parallel ^{1 + \frac{1}{\alpha }}}}\end{align}

Step 2. Design of the control scheme for ${\boldsymbol{{u}}_i}$ .

The distributed fixed-time control scheme is described as:

(34) \begin{align}{{\boldsymbol{{u}}_i}} &= { {\boldsymbol{{g}}^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right)\left[ { - \frac{{{{\left\| {{\boldsymbol{{F}}_i}} \right\|}^2}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}}{{2{\iota _f}}} - {k_3}{\rm{si}}{{\rm{g}}^{2\alpha - 1}}\!\left( {{\boldsymbol\kappa _i}} \right)} \right]}\nonumber \\[5pt] & \quad +{ {\boldsymbol{{g}}^{ - 1}}\!\left( {{\boldsymbol{{x}}_i}} \right)\left[ { - {k_4}{\rm{si}}{{\rm{g}}^{\beta - 1 + \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right) - {k_5}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right]}\end{align}

where ${\boldsymbol{{F}}_i} = {\boldsymbol{{f}}_i}\!\left( {{\boldsymbol{{x}}_i},{\boldsymbol{{v}}_i}} \right) - {\dot{\hat{\boldsymbol{{v}}}}_{i,0}} - {\ddot{\boldsymbol\Delta}_{i,0}} + {k_1}\beta {\rm{si}}{{\rm{g}}^{\beta - 1}}\!\left( {{\boldsymbol\sigma _i}} \right){\dot{\boldsymbol\sigma}_i}$ , ${k_3},{k_4},{k_5}$ and ${\iota _f}$ are positive parameters.

Theorem 3. Suppose Assumption 2, 3 hold. With the implementation of the fixed-time observer Equation (9) and the control scheme Equation (34). The control parameters satisfy $2{k_5}{\chi _{min}}\!\left( {1 - \delta } \right) \gt {\iota _d}$ , ${k_1} \gt \frac{{{c_1}{k_1}\beta }}{{1 + \beta }}$ , ${k_2} \gt \frac{{1 + {2^{1 - \alpha }}}}{{1 + \alpha }} + \frac{{{c_3}\alpha }}{{1 + \alpha }}$ , ${k_4}{\chi _{min}}(1 - \delta )(8n{)^{\frac{{1 + \beta }}{2}}} \gt \frac{{{c_1}{k_1}}}{{1 + \beta }}$ , ${k_3}{\chi _{min}}(1 - \delta ) \gt {c_1}{c_2} + \frac{{{c_1} + {c_3}}}{{1 + \alpha }} + \frac{{{2^{1 - \alpha }}\alpha }}{{1 + \alpha }}$ , ${o_1}(1 - \iota )(8n{)^{\frac{{\alpha - 1}}{{2\alpha }}}} \gt \frac{{({c_1} + 1)\alpha }}{{(1 + \alpha ){\lambda _{\min }}{{(H)}^{1 + \frac{1}{\alpha }}}}}$ . Then the tracking errors ${\boldsymbol{{e}}_{xi}}$ , ${\boldsymbol{{e}}_{vi}}$ for each follower spacecraft can converge to a region in fixed time ${T_2}$ .

Proof. Consider the following Lyapunov function candidate [Reference Qian and Lin30]:

(35) \begin{align}{V_2} = \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \int\limits_{{\gamma _{dij}}}^{{\gamma _{ij}}} {\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( s \right) - {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{dij}}} \right)} \right){\rm{d}}s\end{align}

then the time derivative of ${V_2}$ is given by

(36) \begin{align}{\dot V_2} = {\dot \gamma ^T}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( \boldsymbol\kappa \right) + \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \int\limits_{{\gamma _{dij}}}^{{\gamma _{ij}}} \frac{{\rm{d}}}{{{\rm{d}}t}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( s \right) - {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{dij}}} \right)} \right){\rm{d}}s\end{align}

For the second term in the Equation (36), we could rewrite it as

(37) \begin{align}\sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \int\limits_{{\gamma _{dij}}}^{{\gamma _{ij}}} \frac{{\rm{d}}}{{{\rm{d}}t}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( s \right) - {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{dij}}} \right)} \right){\rm{d}}s = \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 \!\left( {2 - \alpha } \right){k_2}^{\frac{1}{\alpha }}{\dot \sigma _{ij}}{\Upsilon _{ij}}\end{align}

where

(38) \begin{align}{\Upsilon _{ij}} = \int\limits_{{\gamma _{ij}}}^{{\gamma _{ij}}} {\rm{si}}{{\rm{g}}^{1 - \alpha }}\!\left( {{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( s \right) - {\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left( {{\gamma _{dij}}} \right)} \right){\rm{d}}s \le {2^{1 - \alpha }}\left| {{\kappa _{ij}}} \right|\end{align}

which is obtained by Lemma 2 and Lemma 5 .

Substituting Equation (37) and Equation (38) into Equation (36), it can yield

(39) \begin{align}{{{\dot V}_2}} &\le { {{\dot{\boldsymbol\gamma}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}(\boldsymbol\kappa ) + (2 - \alpha )k_2^{\frac{1}{\alpha }}{2^{1 - \alpha }}{{\dot{\boldsymbol\sigma}}^{\rm{T}}}|\kappa |}\nonumber \\[5pt] &\le { {{\dot{\boldsymbol\gamma}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}(\kappa )}\nonumber\\[5pt] & +{ (2 - \alpha )k_2^{\frac{1}{\alpha }}{\lambda _{\max }}(\boldsymbol{{H}})\parallel \boldsymbol\kappa \parallel \left\| {\boldsymbol\gamma - {\boldsymbol\gamma _d} + \tilde{\boldsymbol{{v}}} - {k_1}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\sigma ) - {k_2}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma )} \right\|}\end{align}

Noting that $\left\| {\boldsymbol\gamma - {\boldsymbol\gamma _d}} \right\| \le \sum\limits_{i = 1}^n \sum\limits_{j = 1}^8 {2^{1 - \alpha }}{\left| {{\boldsymbol\kappa _{ij}}} \right|^\alpha } \le {2^{1 - \alpha }}{\!\left( {8n} \right)^{\frac{{1 - \alpha }}{2}}}{\left\| \boldsymbol\kappa \right\|^\alpha }$ , $\left\| {{k_1}{\rm{si}}{{\rm{g}}^\beta }\!\left( \boldsymbol\sigma \right)} \right\| \le {k_1}{\left\| \boldsymbol\sigma \right\|^\beta }$ , $\left\| {{k_2}{\rm{si}}{{\rm{g}}^\alpha }\!\left( \boldsymbol\sigma \right)} \right\| \le {k_2}{(8n)^{\frac{{1 - \alpha }}{2}}}{\left\| \boldsymbol\sigma \right\|^\alpha }$ where Lemma 2 , Lemma 5 and Lemma 6 are used, we have

(40) \begin{align}\dot{V}_{2}&\le{ {{\dot{\boldsymbol\gamma}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}(\boldsymbol\kappa ) + {c_1}{c_2}\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }} + {c_1}\parallel \boldsymbol\kappa \parallel \parallel \tilde{\boldsymbol{{v}}}\parallel + {c_1}{k_1}\parallel \boldsymbol\kappa \parallel \parallel \boldsymbol\sigma {\parallel ^\beta } + {c_3}\parallel \boldsymbol\kappa \parallel \parallel \boldsymbol\sigma {\parallel ^\alpha }}\nonumber \\[5pt]&\le { {{\dot{\boldsymbol\gamma}}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}(\boldsymbol\kappa ) + {c_1}{c_2}\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }} + {c_1}\!\left( {\frac{1}{{1 + \alpha }}\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }} + \frac{\alpha }{{1 + \alpha }}\parallel \tilde{\boldsymbol{{v}}}{\parallel ^{1 + \frac{1}{\alpha }}}} \right)}\nonumber \\[5pt] &+{ {c_1}{k_1}\!\left( {\frac{1}{{1 + \beta }}\parallel \boldsymbol\kappa {\parallel ^{1 + \beta }} + \frac{\beta }{{1 + \beta }}\parallel \boldsymbol\sigma {\parallel ^{1 + \beta }}} \right) + {c_3}\!\left( {\frac{1}{{1 + \alpha }}\boldsymbol\kappa {\parallel ^{1 + \alpha }} + \frac{\alpha }{{1 + \alpha }}\parallel \boldsymbol\sigma {\parallel ^{1 + \alpha }}} \right)}\nonumber \\[5pt] &= { {{\dot{\boldsymbol\gamma }}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}(\boldsymbol\kappa ) + \!\left( {{c_1}{c_2} + \frac{{{c_1} + {c_3}}}{{1 + \alpha }}} \right)\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }} + \frac{{{c_1}{k_1}}}{{1 + \beta }}\parallel \boldsymbol\kappa {\parallel ^{1 + \beta }}}\nonumber \\[5pt] &+{ \frac{{{c_3}\alpha }}{{1 + \alpha }}\parallel \boldsymbol\sigma {\parallel ^{1 + \alpha }} + \frac{{{c_1}{k_1}\beta }}{{1 + \beta }}\parallel \boldsymbol\sigma {\parallel ^{1 + \beta }} + \frac{{{c_1}\alpha }}{{1 + \alpha }}\parallel \tilde{\boldsymbol{{v}}}{\parallel ^{1 + \frac{1}{\alpha }}}}\end{align}

where ${c_1} = \!\left( {2 - \alpha } \right){k_2}^{\frac{1}{\alpha }}{\lambda _{\max }}\!\left( \boldsymbol{{H}} \right)$ , ${c_2}= { 2^{1 - \alpha }}{\!\left( {8n} \right)^{\frac{{1 - \alpha }}{2}}}$ , ${c_3} = {c_1}{k_2}{\left( {8n} \right)^{\frac{{1 - \alpha }}{2}}}$ .

Noting that ${\boldsymbol\gamma _i} = {\boldsymbol{{v}}_i} - {\hat{\boldsymbol{{v}}}_i} - {\dot{\boldsymbol\Delta}_{i,0}} + {k_1}{\rm{si}}{{\rm{g}}^\beta }\!\left( {{\boldsymbol\sigma _i}} \right)$ , we substitute ${\boldsymbol\gamma _i}$ , Equation (4) and control scheme Equation (34) into the first term in the left hand of Equation (40) and we can get

(41) \begin{align}& \quad {\dot{\boldsymbol\gamma}_i^{\rm{T}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}\nonumber\\[5pt] & = { - {{\left[ {{k_3}\boldsymbol\chi \!\left( {{\boldsymbol{{u}}_i}} \right){\rm{si}}{{\rm{g}}^{2\alpha - 1}}\!\left( {{\boldsymbol\kappa _i}} \right) + {k_4}\boldsymbol\chi \!\left( {{u_i}} \right){\rm{si}}{{\rm{g}}^{\beta - 1 + \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right]}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}\nonumber \\[5pt] & \quad + {} { {{\left[ {{k_5}\boldsymbol\chi \!\left( {{\boldsymbol{{u}}_i}} \right){\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right]}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}\nonumber \\[5pt] & \quad+ { {{\left[ {{\boldsymbol{{F}}_i} + {\boldsymbol{{d}}_i} + \boldsymbol{{g}}({x_i}){\boldsymbol{{d}}_{sqi}} - \boldsymbol\chi \!\left( {{u_i}} \right)\frac{{{{\left\| {{F_i}} \right\|}^2}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}}{{2\iota _f^2}}} \right]}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}\nonumber \\[5pt] & \le { - {k_3}{\chi _{min}}(1 - \delta )\boldsymbol\kappa _i^{\rm{T}}{\rm{si}}{{\rm{g}}^\alpha }\!\left( {{\boldsymbol\kappa _i}} \right) - {k_4}{\boldsymbol\chi _{min}}(1 - \delta )\boldsymbol\kappa _i^{\rm{T}}{\rm{si}}{{\rm{g}}^\beta }\!\left( {{\boldsymbol\kappa _i}} \right) - {k_5}{\boldsymbol\chi _{min}}(1 - \delta ){{\left\| {{\boldsymbol\kappa _i}} \right\|}^{4 - 2\alpha }}}\nonumber \\[5pt] & \quad - { {\chi _{min}}(1 - \delta )\frac{{{{\left\| {{\boldsymbol{{F}}_i}} \right\|}^2}{{\left\| {{\boldsymbol\kappa _i}} \right\|}^{4 - 2\alpha }}}}{{2{\iota _f}}} + {{\left| {{\boldsymbol{{F}}_i}} \right|}^{\rm{T}}}\left| {{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right| + {{\!\left( {{\boldsymbol{{d}}_i} + \boldsymbol{{g}}({\boldsymbol{{x}}_i}){\boldsymbol{{d}}_{sqi}}} \right)}^{\rm{T}}}\!\left[ {{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right]}\nonumber {}\\[5pt] &\le { - {k_3}{\chi _{min}}\!\left( {1 - \delta } \right){\boldsymbol\kappa _i}^{\rm{T}}{\rm{si}}{{\rm{g}}^\alpha }\!\left( {{\boldsymbol\kappa _i}} \right) - {k_4}{\chi _{min}}\!\left( {1 - \delta } \right){\boldsymbol\kappa _i}^{\rm{T}}{\rm{si}}{{\rm{g}}^\beta }\!\left( {{\boldsymbol\kappa _i}} \right)}\nonumber \\[5pt] & \quad + { \frac{{{\iota _f}}}{{2{\chi _{min}}(1 - \delta )}} + \frac{{d_{\max }^2 + 8({g_{max}}{u_{\min }}{)^2}}}{{2{\iota _d}}}}\end{align}

where the inequality ${\left| {{\boldsymbol{{F}}_i}} \right|^{\rm{T}}}\left| {{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)} \right| \le {\chi _{min}}\!\left( {1 - \delta } \right)\frac{{{{\left\| {{\boldsymbol{{F}}_i}} \right\|}^2}{{\left\| {{\boldsymbol\kappa _i}} \right\|}^{4 - 2\alpha }}}}{{2{\iota _f}}} + \frac{{{\iota _f}}}{{2{\chi _{min}}\!\left( {1 - \delta } \right)}}$ and Assumption 2 have been used.

Noting that ${\dot{\boldsymbol\gamma}^{\rm{T}}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( \boldsymbol\kappa \right) = \sum\limits_{i = 1}^n {\dot{\boldsymbol\gamma}_i}^{\rm{T}}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\kappa _i}} \right)$ , we can rewrite Equation (40) as

(42) \begin{align}{{{\dot V}_2}} &\le{ - \left[ {{k_3}{\chi _{min}}(1 - \delta ) - {c_1}{c_2} - \frac{{{c_1} + {c_3}}}{{1 + \alpha }}} \right]\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }}}\nonumber \\[5pt] & { - \left[ {{k_4}{\chi _{min}}(1 - \delta )(8n{)^{\frac{{1 + \beta }}{2}}} - \frac{{{c_1}{k_1}}}{{1 + \beta }}} \right]\parallel \boldsymbol\kappa {\parallel ^{1 + \beta }}}\nonumber \\[5pt] &+ { \frac{{{c_3}\alpha }}{{1 + \alpha }}\parallel \boldsymbol\sigma {\parallel ^{1 + \alpha }} + \frac{{{c_1}{k_1}\beta }}{{1 + \beta }}\parallel \boldsymbol\sigma {\parallel ^{1 + \beta }} + \frac{{{c_1}\alpha }}{{1 + \alpha }}\parallel \tilde{\boldsymbol{{v}}}{\parallel ^{1 + \frac{1}{\alpha }}}}\nonumber \\[5pt] &+{ \frac{{n{\iota _f}}}{{2{\chi _{min}}(1 - \delta )}} + \frac{{nd_{\max }^2 + 8n{{({g_{max}}{u_{\min }})}^2}}}{{2{\iota _d}}}}\end{align}

Consider the overall Lyapunov function $V = {V_{\zeta 1}} + {V_1} + {V_2}$ and the time derivative of $V$ is given by the sum of Equation (17), Equation (33) and Equation (42).

(43) \begin{align}{\dot V } &\le{ - {p_{\sigma \beta }}\parallel \boldsymbol\sigma {\parallel ^{1 + \beta }} - {p_{\sigma \alpha }}\parallel \boldsymbol\sigma {\parallel ^{1 + \alpha }} - {p_{\kappa \beta }}\parallel \boldsymbol\kappa {\parallel ^{1 + \beta }} - {p_{\kappa \alpha }}\parallel \boldsymbol\kappa {\parallel ^{1 + \alpha }}}\nonumber\\[5pt] & \quad - {{p_{\zeta \beta }}\parallel \boldsymbol\zeta {\parallel ^{1 + \beta }} - {p_{\zeta \alpha }}\parallel \boldsymbol\zeta {\parallel ^{1 + \alpha }} - {p_\zeta }\parallel \boldsymbol\zeta {\parallel ^{1 + \frac{1}{\alpha }}} + {p_b}}\nonumber \\[5pt] & \le { - {l_\alpha }{{\!\left( {\parallel \boldsymbol\zeta {\parallel ^2} + \parallel \boldsymbol\sigma {\parallel ^2} + \parallel \boldsymbol\kappa {\parallel ^2}} \right)}^{\frac{{1 + \alpha }}{2}}} - {l_\beta }{{\!\left( {\parallel \boldsymbol\zeta {\parallel ^2} + \parallel \boldsymbol\sigma {\parallel ^2} + \parallel \boldsymbol\kappa {\parallel ^2}} \right)}^{\frac{{1 + \beta }}{2}}} + {p_b}}\end{align}

where ${p_{\sigma \beta }} = {k_1} - \frac{{{c_1}{k_1}\beta }}{{1 + \beta }}$ , ${p_{\sigma \alpha }} = {k_2} - \frac{{1 + {2^{1 - \alpha }}}}{{1 + \alpha }} - \frac{{{c_3}\alpha }}{{1 + \alpha }}$ , ${p_{\kappa \beta }} = {k_4}{\chi _{min}}(1 - \delta )(8n{)^{\frac{{1 + \beta }}{2}}} - \frac{{{c_1}{k_1}}}{{1 + \beta }}$ , ${p_{\kappa \alpha }} = {k_3}{\chi _{min}}(1 - \delta ) - {c_1}{c_2} - \frac{{{c_1} + {c_3}}}{{1 + \alpha }} - \frac{{{2^{1 - \alpha }}\alpha }}{{1 + \alpha }}$ , ${p_{\zeta \beta }} = {o_2}(1 - \iota )(8n{)^{\frac{{1 - \beta }}{2}}}$ , ${p_{\zeta \alpha }} = {o_3}(1 - \iota )$ , ${p_\zeta } = {o_1}(1 - \iota )(8n{)^{\frac{{\alpha - 1}}{{2\alpha }}}} - \frac{{({c_1} + 1)\alpha }}{{(1 + \alpha ){\lambda _{\min }}{{(H)}^{1 + \frac{1}{\alpha }}}}}$ , ${p_b} = \frac{{n{\iota _f}}}{{2{\chi _{min}}(1 - \delta )}} + \frac{{nd_{\max }^2 + 8n{{({g_{max}}{u_{\min }})}^2}}}{{2{\iota _d}}} + \frac{{{o_4}(1 - \iota )\delta k}}{4}$ , ${l_\alpha } = {\rm{max}}\!\left( {{p_{\sigma \alpha }},{p_{\kappa \alpha }},{p_{\zeta \alpha }}} \right)$ , ${l_\beta } = {\rm{max}}\!\left( {{p_{\sigma \beta }},{p_{\kappa \beta }},{p_{\zeta \beta }}} \right)$ are positive constants.

Figure 3. Communication topology.

Similar to the method used in Equation (38), we can derive that ${V_2} \le {2^{1 - \alpha }}\parallel \boldsymbol\kappa {\parallel ^2}$ . Then, $V = {V_{**\zeta 1}} + {V_1} + {V_2}$ can be rewritten as

(44) \begin{align}V &\le{ \frac{1}{{2{\lambda _{{\rm{min}}}}(\boldsymbol{{H}})}}\parallel \boldsymbol\zeta {\parallel ^2} + \frac{1}{{2{\lambda _{{\rm{max}}}}(\boldsymbol{{H}})}}\parallel \boldsymbol\sigma {\parallel ^2} + {2^{1 - \alpha }}\parallel \boldsymbol\kappa {\parallel ^2}}\nonumber \\[5pt] &\le { l\!\left( {\parallel \boldsymbol\zeta {\parallel ^2} + \parallel \boldsymbol\sigma {\parallel ^2} + \parallel \boldsymbol\kappa {\parallel ^2}} \right)}\end{align}

where $l = {\rm{max}}\!\left( {\frac{1}{{2{\lambda _{{\rm{min}}}}(\boldsymbol{{H}})}}{{,2}^{1 - \alpha }}} \right)$ . Substituting Equation (44) into Equation (43), it can yield

(45) \begin{align}\dot V \le - \frac{{{l_\alpha }}}{{{l^{\frac{{1 + \alpha }}{2}}}}}{V^{\frac{{1 + \alpha }}{2}}} - \frac{{{l_\beta }}}{{{l^{\frac{{1 + \beta }}{2}}}}}{V^{\frac{{1 + \beta }}{2}}} + {p_b}\end{align}

Table 1. Initial orbital elements of the virtual leader

Table 2. Desired states of each follower spacecraft

Table 3. Initial states of each follower spacecraft

By using Lemma 1 , we can conclude that $\boldsymbol\sigma $ and $\boldsymbol\kappa $ converge to a region within fixed time ${T_2} \le 1/\left[ {\frac{{{l_\alpha }}}{{{l^{\frac{{1 + \alpha }}{2}}}}}\!\left( {1 - \frac{{1 + \alpha }}{2}} \right)} \right] + 1/\left[ {\frac{{{l_\beta }}}{{{l^{\frac{{1 + \alpha }}{2}}}}}\!\left( {\frac{{1 + \beta }}{2} - 1} \right)} \right]$ , which in turn implies that $\boldsymbol\gamma $ and ${\boldsymbol\gamma _d}$ converge to zero within fixed time as well. Since $\boldsymbol\sigma = \boldsymbol{{He}}_{x}$ and ${\boldsymbol{{e}}_v} = \boldsymbol\gamma - {\boldsymbol\gamma _d} + \tilde{\boldsymbol{{v}}} - {k_1}{\rm{si}}{{\rm{g}}^\beta }(\boldsymbol\sigma ) - {k_2}{\rm{si}}{{\rm{g}}^\alpha }(\boldsymbol\sigma )$ , we can obtain that the tracking errors ${\boldsymbol{{e}}_x}$ and ${\boldsymbol{{e}}_v}$ converge to the origin within fixed time.

Remark 8. It should be mentioned that Equation (8) is inspired by [Reference Li, Wang and Wang39], which only achieved an asymptotic convergence result. In contrast, we deal with the actuator saturation and input quantisation problems under the framework of the “adding a power integrator” technique in an ingenious way, by embedding $ - \frac{{{{\left\| {{\boldsymbol{{F}}_i}} \right\|}^2}{\rm{si}}{{\rm{g}}^{2 - \alpha }}\!\left( {{\boldsymbol\kappa _i}} \right)}}{{2{\iota _f}}}$ in Equation (34).

Remark 9. Recalling the time derivative of ${V_1}$ in Equation (33), it is noted that $\frac{\alpha }{{1 + \alpha }}\parallel \tilde{\boldsymbol{{v}}}{\parallel ^{1 + \frac{1}{\alpha }}}$ would be compensated by $ - {o_1}(1 - \iota )(8n{)^{\frac{{\alpha - 1}}{{2\alpha }}}}\parallel \boldsymbol\zeta {\parallel ^{1 + \frac{1}{\alpha }}}$ in Equation (17), which is generated by the nonlinear term $ - {o_1}{\rm{si}}{{\rm{g}}^{\frac{1}{\alpha }}}\!\left[ {{\boldsymbol\zeta _i}({t_k})} \right]$ in Equation (9). By doing so, we guarantee the fixed time stability of the entire closed-loop system without recourse to the separation principle.

4.0 Numerical simulations and results

The performance of the proposed observer Equation (9), Equation (24) and the control scheme Equation (34) are tested in Numerical simulations. Four spacecrafts together with one virtual leader are considered, as shown in Fig. 3. The adjacency matrix $A$ is set as $[0,0.1,0,0;\;0.1,0,0.1,0;\;0,0.1,0,0.1;\;0,0,0.1,0]$ .

Figure 4. Errors of observation.

Figure 5. Errors of attitude and position.

Figure 6. Control of the orbit.

Statistics show that approximately 73% of SFF missions operate in low Earth orbit (LEO), and the vast majority of SFF missions are executed by small spacecraft with a mass of less than 500 kg [Reference Di Mauro, Lawn and Bevilacqua25]. The simulation environment is therefore set up, so that spacecraft with a mass of about 100 kg follow the virtual leader in LEO. For follower spacecrafts, the mass are chosen as ${m_1} = 102{\rm{kg}}$ , ${m_2} = 104{\rm{kg}}$ , ${m_3} = 103{\rm{kg}}$ , ${m_4} = 99{\rm{kg}}$ and the inertia matrices are chosen as ${\boldsymbol{{J}}_1} = [15.2, - 0.04,0.05; - 0.04,17.3,0.02;0.05,0.02,19.5]$ , ${\boldsymbol{{J}}_2} = [14.7,0.01, - 0.06;0.01, 17.3,0.03; - 0.06,0.03,19.5]$ , ${\boldsymbol{{J}}_3} = \left[ {15.4, - 0.043,0; 0.03,17.1,0.01;0,0.01,20.7} \right]$ , ${\boldsymbol{{J}}_4} = [ 14.9,-0.04,0; -0.04,16.9,0.02;0,0.02,19.8 ]$ . The external disturbances are selected as ${\boldsymbol{{f}}_{ip}} = 0.01[0.9{\rm{sin}}(0.075t + 0.8i - 0.3),1{\rm{sin}}(0.09t + 0.8i + 0.8),1.1{\rm{sin}}(0.06t + 0.8i + {0.2)]^T}N$ and ${\boldsymbol\tau _{ip}} = 0.001[0.9{\rm{sin}}(0.11t + 0.8i + 0.2),1{\rm{sin}}(0.125t + 0.8i + 1.2),1.1{\rm{sin}}(0.1125t + 0.8i + {2.2)]^T}N \cdot m$ . Table 1 gives the virtual leader’s initial orbital elements.

The controller and observer parameters of the ith follower spacecraft are chosen as $\alpha = 0.9$ , $\beta = 2$ , ${o_1} = 0.5$ , ${o_2} = 0.5$ , ${o_3} = 0.005$ , ${o_4} = 100$ , ${k_{o1}} = 0.0012$ , ${k_{o2}} = 7.8240$ , ${k_{o3}} = 3$ , $\iota = 0.1$ , ${\delta _o} = 0.001$ , ${k_1} = 1$ , ${k_2} = 0.05$ , ${k_3} = 0.1$ , ${k_4} = 0.5$ , ${k_5} = 0.5$ . The dead zone of the control force and the control torque are $0.005N$ and $0.0004N \cdot m$ , quantified density $\delta = 0.2$ . As shown in Table 2, we expect each spacecraft to be 200m from the others and to have the same attitude. ${\vec{\boldsymbol{{r}}}_0}$ represents the virtual leader’s position vector. ${\vec{\boldsymbol{{r}}}_i} - {\vec{\boldsymbol{{r}}}_0}$ denotes the position deviation for ith follower spacecraft from the virtual leader. Besides, Table 3 shows the initial states of each spacecraft. In order to describe the attitude more intuitively, we have used Euler angles ${\boldsymbol\theta _i}(0)$ and ${\boldsymbol\theta_i}(t)$ in the initial settings and the subsequent simulation figures. It is worth pointing out that we still use the dual quaternion in the control algorithm, where the Euler angles are used only to express the attitude change intuitively. The observer’s estimated initial value for the ith follower spacecraft is taken as its initial state value.

4.1 Evaluation of control performance and convergence

To display the observation errors intuitively, we depict ${\widetilde{\boldsymbol{{q}}}_i},{\widetilde{\boldsymbol{{r}}}_i}$ in Fig. 4 instead of ${\widetilde{\boldsymbol{{v}}}_i}$ , where ${\widetilde{\boldsymbol{{q}}}_i}$ and ${\widetilde{\boldsymbol{{r}}}_i}$ are observation errors of the virtual leader’s vector quaternion and position. It can be observed that all follower spacecrafts can estimate the virtual leader’s velocity in a fixed time.

The time responses of the attitude quaternion and position tracking errors ${\boldsymbol\theta _e},{\boldsymbol{{r}}_e}$ of each follower spacecraft are given in Fig. 5. The trajectory of the spacecraft’s position in the uncontrolled and the controlled case during an orbital period of 5,578 seconds is shown in Fig. 6. It is important to note that Fig. 6 uses an inhomogeneous coordinate system for more visual comparison. In addition, the minimum value, maximum value, and standard deviation of the controller convergence errors from 250 to 300 seconds are explicitly described in Table 4, which is needed while implementing the engineering solution of the mathematical treatment.

Figure 7 illustrate the control force and the control torque bounded by $0.52N$ and $0.02N \cdot m$ . As shown in these figures, the control output in the initial stage is saturated. Obviously, even in the presence of actuator saturation and input quantisation, the SFF system can estimate the states of the virtual leader effectively and track the trajectory of the virtual leader while maintaining the desired formation configuration.

Table 4. Desired states of each follower spacecraft

Figure 7. Control torque and force.

4.2 Comparison

To highlight the advantages of the proposed triggering condition Equation (24), we use the trigger conditions in [Reference Liu, Zhang, Yu and Sun17] as a comparison, which can be expressed as

(46) \begin{align}{{g_{ij}}(t)} = { \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s - \frac{\iota }{{1 + \iota }}\left| {{{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + \frac{{{k_{o4}}}}{\iota }} \right|}\end{align}

where ${k_{o4}}$ is positive parameter. Meanwhile, to obtain a fair comparison, take the two cases of ${k_{o4}} = {k_{o1}}\sqrt e $ and ${k_{o4}} = {k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} $ for comparison, where ${k_{o1}}\sqrt e $ and ${k_{o1}}\sqrt {{e^{1 + {k_{o2}}}}} $ are the upper and lower bounds of the nonlinear term ${k_{o1}}\sqrt {{e^{1 + {k_{o2}}\tanh \!\left( {{k_{o3}}\left| {{\xi _{ij}}\!\left( {{t_k}} \right)} \right|} \right)}}} $ in Equation (24), respectively. The two cases are expressed as

(47) \begin{align}{{g_{ij}}(t)} = { \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s - \frac{\iota }{{1 + \iota }}\left| {{{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + \frac{{{k_{o1}}}}{\iota }\sqrt e } \right|}\end{align}
(48) \begin{align}{{g_{ij}}(t)} = { \int_{{t_k}}^t {\eta _1}{h_{\max }}\left| {{\eta _{2ij}}\!\left( {{t_k}} \right) + {B_3}} \right|{\rm{d}}s - \frac{\iota }{{1 + \iota }}\left| {{{\dot{\hat{v}}}_{ij}}\!\left( {{t_k}} \right) + \frac{{{k_{o1}}}}{\iota }\sqrt {{e^{1 + {k_{o2}}}}} } \right|}\end{align}

According to the statistical analysis of the trigger simulation data, Figs. (8)–(9) can be obtained. Figure (8) plots the comparison results of the trigger threshold. The triggering intervals of the SFF system is presented in Fig. (9).

Figure 8. Curves of triggering condition thresholds.

Figure 9. Curves of triggering condition thresholds.

We divide the convergence process of the observation errors into the transient-state phase and the steady-state phase for discussion, aiming to demonstrate the superiority of the proposed ETM in detail. (A). In the transient-state phase. Equation (24) has a larger trigger threshold than Equation (47), but almost the same as Equation (48). Hence, the trigger times of Equation (24) is in the middle of the trigger times of Equation (47) and Equation (48). (B). In the steady-state phase. Equation (24) has a smaller trigger threshold than Equation (48), but almost the same as Equation (47). Therefore, the trigger times of Equation (24) and Equation (47) are almost the same. Further, The norm distribution of the observation errors of Equation (24) and Equation (47) are almost the same.

As can be observed in these discussions and figures, Equation (24) achieves the same steady-state performance as Equation (47) with less communication penalty. It is because that the lower limit of the trigger threshold of Equation (24) fluctuates with the accuracy compared with the stable lower limit of the trigger threshold of Equation (46). In other words, Equation (24) realises the amplification of the lower limit of trigger threshold when the error is large and the reduction of the lower limit of trigger threshold when the error is small, that is the proposed ETM achieves the same steady-state performance with less communication penalty.

However, there is no win-win situation in control systems. In fact, although Equation (24) has reduced communication penalty compared to Equation (47), it has increased fuel consumption. Figure (10) illustrates this with an example ${\widetilde{\boldsymbol\theta}_x}$ and ${\dot{\hat{\boldsymbol\theta}}_x}$ of the follower spacecraft 1. For an observer, the fuel it consumes is its control input ${\dot{\hat{\boldsymbol\theta}}_x}$ . The quantitative analysis of the different ETMs is given in Table 5, where errors distribution denotes the range of $\sum\nolimits_{i = 1}^4 \sum\nolimits_{j = 1}^6 |{z_{ij}}|/24$ with observer error ${z_{ij}}$ , and total input denotes $\sum\nolimits_{i = 1}^4 \sum\nolimits_{j = 1}^6 \int o{u_{ij}}(t)$ with observer control input $o{u_{ij}}$ . It can be calculated that ETM Equation (24) achieves a 8.98% saving in communication resources at a cost of 2.62% fuel consumption compared to ETM Equation (47).

Figure 10. Comparison of observer performance under different ETMs.

Remark 10. It should be noted that the execution of the observer is fully digital and does not really require the consumption of fuel, which is why fuel is quoted in the text above. In other words, the application of the ETM in the observer designed in this paper circumvents the cost of its fuel consumption.

Please note that ETM is only embedded in the observer and is not mounted on the controller. Figure (11) illustrates the impact of different ETMs on the controller’s performance with an example ${\boldsymbol{\theta}_{ex}}$ of the follower spacecraft 1. It can be seen that the observer embedding different ETMs has no effect on the transient performance of the controller, which is due to the fact that the observer converges much faster than the controller converges. Equation (24) and (47) have the same effect on the steady-state performance of the controller, and Equation (48) causes the steady-state performance of the controller to deteriorate. This is because the observation error of Equation (48) is higher than that of Equation (24) and Equation (47) (as shown in Fig. 10), and the observation error affects the convergence error of the controller.

Table 5. Communication analysis of different ETMs for all spacecrafts in 10 seconds

Figure 11. Impact of different ETMs on the controller’s performance.

Besides, to emphasise the advantages of the proposed control scheme, we reproduced the asymptotic controller [Reference Fan, Chen, Liu and Cao40] and the finite-time controller [Reference Zhang and Li41] in numerical simulations for comparison. We consider two performance indices for an intuitive comparison: (A) The Attitude Station-Keeping Absolute Error (ASKAE) $\sum\limits_{i = 1}^n \left\| {{\boldsymbol{{q}}_{e,i}}} \right\|$ . (B) The Position Station-Keeping Absolute Error (PSKAE) $\sum\limits_{i = 1}^n \left\| {{\boldsymbol{{r}}_{e,i}}} \right\|$ . Figure 12 gives the two performance indicators’ comparison curves for each control scheme. Obviously, the proposed fixed time controllers demonstrate faster convergence rates in both performance indices.

Figure 12. The comparison curves of the two performance indicators.

5.0 Conclusion

This paper investigated the attitude and position coupled tracking control problem of SFF systems with actuator saturation and input quantisation under an undirected communication graph. Since the communication bandwidth is limited, a resource-saving distributed control scheme was developed based on an event-triggered observer and the adding a power integrator technique. It is worth remarking that the novel observer realised the amplification of the trigger threshold when the error is large, but neither increased the convergence time nor affected the final accuracy. Besides, the control scheme compensated the actuator saturation and input quantisation in an ingenious way. Finally, a simulation of a four-spacecraft on LEO is carried out to illustrate the successful application of the proposed distributed control algorithm. The control scheme subject to the limited communication under the directed graph will be investigated in our future work.

Acknowledgements

This work was supported in part by the National Natural Science Foundation of China (grant number: 61304108) and National Key R&D Program of China (grant number: 2020YFC2200600).

Data availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

Mauro, G., Spiller, D., Bevilacqua, R. and D’Amico, S. Spacecraft formation flying reconfiguration with extended and impulsive maneuvers, J. Franklin Inst., 2019, 356, (6), pp 34743507.CrossRefGoogle Scholar
Shaw, G.B., Miller, D.W. and Hastings, D.E. Generalized characteristics of communication, sensing, and navigation satellite systems, J. Spacecr. Rockets, 2000, 37, (6), pp 801811.CrossRefGoogle Scholar
Ren, W. and Beard, R.W. Decentralized scheme for spacecraft formation flying via the virtual structure approach, J. Guid. Control Dyn., 2004, 27, (1), pp 7382.CrossRefGoogle Scholar
Ahn, C. and Kim, Y. Point targeting of multisatellites via a virtual structure formation flight scheme, J. Guid. Control Dyn., 2009, 32, (4), pp 13301344.CrossRefGoogle Scholar
Jiakang, Z.H.O.U., Guangfu, M.A. and Qinglei, H.U. Delay depending decentralized adaptive attitude synchronization tracking control of spacecraft formation, Chinese J. Aeronaut., 2012, 25, (3), pp 406415.Google Scholar
Wu, B., Wang, D. and Poh, E.K. Decentralized sliding-mode control for spacecraft attitude synchronization under actuator failures, Acta Astronaut., 2014, 105, (1), pp 333343.CrossRefGoogle Scholar
Zhou, J., Hu, Q. and Friswell, M.I. Decentralized finite time attitude synchronization control of satellite formation flying, J. Guid. Control Dyn., 2013, 36, (1), pp 185195.CrossRefGoogle Scholar
Zou, A. and Kumar, K.D. Distributed attitude coordination control for spacecraft formation flying, IEEE Trans. Aerosp. Electron. Syst., 2012, 48, (2), pp 13291346.CrossRefGoogle Scholar
Zou, A., de Ruiter, A.H.J. and Kumar, K.D. Distributed finite-time velocity-free attitude coordination control for spacecraft formations, Automatica, 2016, 67, pp 4653.CrossRefGoogle Scholar
Polyakov, A. Nonlinear feedback design for fixed-time stabilization of linear control systems, IEEE Trans. Autom. Control, 2012, 57, (8), pp 21062110.CrossRefGoogle Scholar
Ren, W. Formation keeping and attitude alignment for multiple spacecraft through local interactions, J. Guid. Control Dyn., 2007, 30, (2), pp 633638.CrossRefGoogle Scholar
Han, G.A.O., Yuanqing, X.I.A., Zhang, X. and Zhang, G. Distributed fixed-time attitude coordinated control for multiple spacecraft with actuator saturation, Chinese J. Aeronaut., 2022, 35, (4), pp 292302.Google Scholar
Tabuada, P. Event-triggered real-time scheduling of stabilizing control tasks, IEEE Trans. Autom. Control, 2007, 52, (9), pp 16801685.CrossRefGoogle Scholar
Dimarogonas, D.V., Frazzoli, E. and Johansson, K.H. Distributed event-triggered control for multi-agent systems, IEEE Tran. Autom. Control, 2012, 57, (5), pp 12911297.CrossRefGoogle Scholar
Nowzari, C., Garcia, E. and Cortés, J. Event-triggered communication and control of networked systems for multi-agent consensus, Automatica, 2019, 105, pp 127.CrossRefGoogle Scholar
Garcia, E., Cao, Y. and Casbeer, D.W. Periodic event-triggered synchronization of linear multi-agent systems with communication delays, IEEE Trans. Autom. Control, 2017, 62, (1), pp 366371.CrossRefGoogle Scholar
Liu, J., Zhang, Y., Yu, Y. and Sun, C. Fixed-time event-triggered consensus for nonlinear multiagent systems without continuous communications, IEEE Trans. Syst. Man Cybern. Syst., 2019, 49, (11), pp 22212229.CrossRefGoogle Scholar
Proskurnikov, A.V. and Mazo, M. Lyapunov event-triggered stabilization with a known convergence rate, IEEE Trans. Autom. Control, 2020, 65, (2), pp 507521.CrossRefGoogle Scholar
Zhang, J., Biggs, J.D., Ye, D. and Sun, Z. Extended-state-observer-based event-triggered orbit-attitude tracking for low-thrust spacecraft, IEEE Trans. Aerosp. Electron. Syst., 2020, 56, (4), pp 28722883.CrossRefGoogle Scholar
Di, F., Li, A., Guo, Y., Xie, C. and Wang, C. Event-triggered sliding mode attitude coordinated control for spacecraft formation flying system with disturbances, Acta Astronaut., 2021, 188, pp 121129.CrossRefGoogle Scholar
Zhou, J., Wen, C. and Yang, G. Adaptive backstepping stabilization of nonlinear uncertain systems with quantized input signal, IEEE Trans. Autom. Control, 2014, 59, (2), pp 460464.CrossRefGoogle Scholar
Liu, Z., Wang, F., Zhang, Y. and Chen, C.L.P. Fuzzy adaptive quantized control for a class of stochastic nonlinear uncertain systems, IEEE Trans. Cybern., 2016, 46, (2), pp 524534.CrossRefGoogle ScholarPubMed
Hu, Q. and Xiao, B. Intelligent proportional-derivative control for flexible spacecraft attitude stabilization with unknown input saturation, Aerosp. Sci. Technol., 2012, 23, (1), pp 6374.CrossRefGoogle Scholar
An, H., Liu, J., Wang, C. and Wu, L. Approximate back-stepping fault-tolerant control of the flexible air-breathing hypersonic vehicle, IEEE/ASME Trans. Mechatron., 2016, 21, (3), pp 16801691.CrossRefGoogle Scholar
Di Mauro, G., Lawn, M. and Bevilacqua, R. Survey on guidance navigation and control requirements for spacecraft formation-flying missions, J. Guid. Control Dyn., 2018, 41, (3), pp 581602.CrossRefGoogle Scholar
Hu, Q., Shi, Y. and Shao, X. Adaptive fault-tolerant attitude control for satellite reorientation under input saturation, Aerosp. Sci. Technol., 2018, 78, pp 171182.CrossRefGoogle Scholar
Polyakov, A. Nonlinear feedback design for fixed-time stabilization of linear control systems, IEEE Trans. Autom. Control, 2011, 57, (8), pp 21062110.CrossRefGoogle Scholar
Zou, A. and Kumar, K.D. Finite-time attitude control for rigid spacecraft subject to actuator saturation, Nonlinear Dyn., 2019, 96, (2), pp 10171035.CrossRefGoogle Scholar
Yu, S., Yu, X., Shirinzadeh, B. and Man, Z. Continuous finite-time control for robotic manipulators with terminal sliding mode, Automatica, 2005, 41, (11), pp 19571964.CrossRefGoogle Scholar
Qian, C. and Lin, W. A continuous feedback approach to global strong stabilization of nonlinear systems, IEEE Trans. Autom. Control, 2001, 46, (7), pp 10611079.CrossRefGoogle Scholar
Tan, C., Yu, X. and Man, Z. Terminal sliding mode observers for a class of nonlinear systems, Automatica, 2010, 46, (8), pp 14011404.CrossRefGoogle Scholar
Wu, J., Liu, K. and Han, D. Adaptive sliding mode control for six-dof relative motion of spacecraft with input constraint, Acta Astronaut., 2013, 87, pp 6476.CrossRefGoogle Scholar
Ploen, S.R., Hadaegh, F.Y. and Scharf, D.P. Rigid body equations of motion for modeling and control of spacecraft formations. part 1: Absolute equations of motion. In Proceedings of the 2004 American Control Conference, volume 4, pages 3646–3653, 2004.CrossRefGoogle Scholar
Sinclair, A.J., Hurtado, J.E. and Junkins, J.L. Application of the cayley form to general spacecraft motion, J. Guid. Control Dyn., 2006, 29, (2), pp 368373.CrossRefGoogle Scholar
Hong, Y., Hu, J. and Gao, L. Tracking control for multi-agent consensus with an active leader and variable topology, Automatica, 2006, 42, (7), pp 11771182.CrossRefGoogle Scholar
Zhou, Z., Zhou, D., Shi, X., Li, R. and Kan, B. Prescribed performance fixed-time tracking control for a class of second-order nonlinear systems with disturbances and actuator saturation, Int. J. Control, 2021, 94, (1), pp 223234.CrossRefGoogle Scholar
Kun, Y.A.N., Mou, C.H.E.N., Qingxian, W.U. and Ronggang, Z.H.U. Robust adaptive compensation control for unmanned autonomous helicopter with input saturation and actuator faults, Chinese J. Aeronaut., 2019, 32, (10), pp 22992310.Google Scholar
Shi, X., Zhou, Z., Zhou, D., Li, R. and Chen, X. Observer-based event-triggered fixed-time control for nonlinear system with full-state constraints and input saturation, Int. J. Control, 2022, 95, (2), pp 432446.CrossRefGoogle Scholar
Li, Y., Wang, C. and Wang, W. Adaptive quantized attitude control for spacecraft with input saturation. In 2017 13th IEEE International Conference on Control & Automation (ICCA), pages 449–454, 2017.CrossRefGoogle Scholar
Fan, R., Chen, X., Liu, M. and Cao, X. Attitude-orbit coupled sliding mode tracking control for spacecraft formation with event-triggered transmission, ISA Trans., 2020, 6, 338–348.Google Scholar
Zhang, B. and Li, F. Adaptive finite-time control for six-degree-of-freedom leader-following spacecraft formation using twistors, Adv. Space Res., 2022, 70, (5), pp 12971311.CrossRefGoogle Scholar
Figure 0

Figure 1. Coordinate system.

Figure 1

Figure 2. Schematic of the distributed control algorithm.

Figure 2

Figure 3. Communication topology.

Figure 3

Table 1. Initial orbital elements of the virtual leader

Figure 4

Table 2. Desired states of each follower spacecraft

Figure 5

Table 3. Initial states of each follower spacecraft

Figure 6

Figure 4. Errors of observation.

Figure 7

Figure 5. Errors of attitude and position.

Figure 8

Figure 6. Control of the orbit.

Figure 9

Table 4. Desired states of each follower spacecraft

Figure 10

Figure 7. Control torque and force.

Figure 11

Figure 8. Curves of triggering condition thresholds.

Figure 12

Figure 9. Curves of triggering condition thresholds.

Figure 13

Figure 10. Comparison of observer performance under different ETMs.

Figure 14

Table 5. Communication analysis of different ETMs for all spacecrafts in 10 seconds

Figure 15

Figure 11. Impact of different ETMs on the controller’s performance.

Figure 16

Figure 12. The comparison curves of the two performance indicators.