1. Introduction
In the realm of space exploration and the development of space infrastructure, spacecraft are pivotal assets, with the attitude control system (ACS) playing an essential role [
1]. The challenge of spacecraft attitude tracking, influenced by unpredictable disturbances and a variety of environmental effects in space, has attracted significant research interest [
2].
Recent advancements have seen an increasing focus on fixed-time attitude control for spacecraft. This approach, distinct from finite-time control, provides a predetermined system convergence time [
3]. However, the complex relationship between system parameters and fixed-time convergence has remained intricate and ambiguous. To overcome these challenges, the concept of predefined time has been introduced as a conservative estimate for the upper limit of fixed time. This method effectively addresses the limitations of fixed-time control, enhancing the predictability and reliability of system behavior [
4].
In response to the operational challenges faced by spacecraft, a range of sophisticated control methods have been developed to meet the stringent mission performance requirements. Optimal control theory has been advanced by Hu et al. [
5], and composite adaptive attitude controllers have been innovatively designed by Liu et al. [
6]. Lyapunov-based design strategies have been further developed by Hu et al. and Wu et al. [
7,
8], and robust control along with its integrated applications has been thoroughly explored in several studies [
9,
10,
11]. Additionally, dual-loop controllers incorporating sliding mode control (SMC) and model predictive control (MPC) have proven effective in addressing the complexities of spacecraft attitude tracking [
12,
13,
14].
The exploration of robust control algorithms is critical for mitigating the effects of uncertainty and disturbances encountered by spacecraft in orbit. Among the various robust control methods, sliding mode control (SMC) and its derivatives have shown exceptional capability in managing unknown interferences and uncertainties. Based on this foundation, finite-time sliding mode control methodologies have been developed to ensure the stable convergence of attitude tracking within predefined timeframes, receiving widespread attention for their practical engineering applications [
15,
16,
17]. Adaptive backstepping techniques, combined with second-order sliding mode control laws, have been proposed by Pukdeboon et al. [
18] to address finite-time attitude tracking challenges in rigid spacecraft without considering the boundary information about uncertainties and interference. Furthermore, Cao et al. [
19] introduced a combination of nonsingular terminal sliding mode control (NTSMC) and the Lyapunov function method to solve the attitude tracking control issues of flexible spacecraft equipped with redundant reaction wheels. Additionally, Pukdeboon et al. [
20] proposed a new finite-time attitude tracking controller that uses a finite-time disturbance observer to estimate disturbances under the presence of external interference, verifying that the tracking error converges to zero within a finite time. Moreover, Chen et al. [
21] employed a new polynomial finite-time performance function to reduce computational requirements and adopted a simpler error transformation to enhance the speed and accuracy of tracking error convergence. From another standpoint, different types of time-varying sliding mode control were employed to solve the tracking problem of expected attitude trajectories that are known or measurable, and the genetic algorithm optimized the parameters to ensure the strong robustness of the controller [
22]. Despite the advances in SMC, there are still notable challenges that need to be addressed. SMC can suffer from chattering effects, which can degrade control performance, making it difficult to ensure fast convergence and high accuracy.
Furthermore, advanced control methods like model predictive control (MPC) have been extensively applied to spacecraft attitude control [
14,
23,
24]. Saki et al. [
25] proposed an adaptive soft switching law based on Lyapunov theory to dynamically adjust controller parameters. Multiple optimal sub-models have been utilized, and direct adaptive structure multi-model predictive control has been employed for multi-operation point and maneuvering range attitude tracking. Zhang et al. [
26] designed a double-layer nonlinear controller based on the high precision and constraint properties of MPC, ensuring robust system dynamics and resilience to external disturbances. Nonetheless, the inherent uncertainties of these systems pose challenges in obtaining high-precision nominal models, rendering MPC overly conservative regarding disturbance robustness, often based on worst-case scenarios. To address these issues, an ideal Generalized Predictive Control (GPC) framework has been employed, which accommodates unknown terms and ensures the optimization of closed-loop system performance while maintaining the efficiency and accuracy of the control method, making it particularly suitable for analyzing trajectory tracking control problems amid parameter uncertainty and unknown constraint interference [
27,
28,
29]. Predictive control methods often require precise system models and can be computationally intensive, making real-time implementation challenging.
To overcome these challenges, this paper introduces a novel predefined-time attitude control scheme for rigid spacecraft, designed to achieve accurate attitude tracking within a predefined timeframe using quaternion representations for spacecraft attitude. The scheme leverages dynamic prediction techniques and the fastest optimizing prediction period to enhance control robustness against parameter uncertainties and disturbances. The main contributions of this study include:
Adaptive Global Optimization for Period Updates: The introduction of new theorems allows for the use of an adaptive global optimization approach to achieve period updates, relaxing the original restrictive conditions in control schemes proposed by Sun G et al. [
30].
Ease of Parameter Adjustment in Predefined-Time Control: Compared to fixed-time control methods, the parameters related to predefined-time control are more easily adjustable, providing a nonconservative upper bound on system stability. This makes the determination of stabilization time for the closed-loop system more straightforward and practical.
Continuous and Robust Control via Terminal Sliding Mode Control and Predictive Methods: This work diverges from the approach of Chen et al. [
31], which employs terminal sliding mode control technology to design predefined-time attitude controllers. Instead, it utilizes a method of optimizing prediction periods and predictive control to derive the control law, ensuring that the proposed control scheme is both continuous and robust.
For the sake of clarity and to facilitate subsequent mathematical expressions, we introduce the following notations in advance:
The symbols R, , N, , and denote the set of real numbers, the set of positive real numbers, the set of natural numbers, the set of positive integers, and the n-dimensional real number space, respectively.
For any matrix , and represent the maximum and minimum eigenvalues of A, respectively.
For any vector , denotes the Euclidean norm of x.
For any vector
, the notation
denotes the corresponding skew-symmetric matrix:
For any vector
and any positive scalar
,
denotes the vector:
where
, for
, is defined as:
The symbol denotes the set of all differentiable functions whose i-th derivative is continuous. If a function has partial derivatives and is continuous up to the k-th order, , then the function f is called a function. Here, represents continuous functions, and represents smooth functions, i.e., functions that have continuous partial derivatives of any order.
2. Spacecraft Attitude Control Problem
As is well known, Euler angles are frequently used to parameterize spacecraft attitude due to their clear and straightforward physical concept. However, Euler angles may encounter singularity or gimbal lock issues when describing the parameterized attitude of a spacecraft. To address these challenges, this paper adopts the quaternion representation for spacecraft attitude. The spacecraft, modeled as a rigid body, operates in low Earth orbit as depicted in
Figure 1a. Its actuators are momentum wheels (MWs), which provide torque about three orthogonal axes defined with respect to the body coordinate system
B.
In the spacecraft attitude control system, when MWs are used as actuators, they contribute additional angular momentum, denoted as
, where
is the inertia matrix and
is the angular velocity of the MWs. Furthermore, the kinematics and dynamics of the spacecraft are governed by the following nonlinear equations [
14,
32]:
where
represents the quaternion describing the orientation of the spacecraft body frame
B relative to the inertial coordinate system
. The vector
and the scalar
ensure
. The angular velocity
is defined in the body frame
B. The inertia matrix of the spacecraft is denoted by
, and
is defined as
.
The torque applied to the system, due to the orbital gravity gradient and unknown, bounded disturbances, is represented by
and
, respectively. Here,
, where
denotes the orbital velocity of a spacecraft in a near-Earth orbit,
represents the gravitational constant of the Earth,
is the distance between the spacecraft and the Earth’s center of mass, and
denotes the third column of the direction cosine matrix of the spacecraft [
15].
Based on the nonlinear dynamics model of the rigid spacecraft’s attitude given by (
1) and (
2), and considering the actual control process of the spacecraft, the following assumptions can be made.
Assumption 1. We assume that the form of the inertia matrix J in (2) is , where and represent the known nominal matrix and uncertain part of J, and satisfies , where is a known positive constant. Assumption 2. We assume that the disturbance is differentiable and bounded, i.e., , where is a known positive constant.
Assumption 3. We assume that the system input u is bounded.
In this study, the desired spacecraft attitude relative to the inertial coordinate system is represented by the quaternion , and the desired attitude angular velocity by . The attitude tracking error between the desired and actual attitudes is denoted as , with components satisfying the norm constraint .
The error quaternion
is defined by the equation:
To aid in the controller design, it is assumed that both
and
are bounded. The angular velocity error
is defined as
, where the matrix
is determined from the tracking error and is expressed as
. In addition,
satisfies
. The kinematic equation governing the attitude tracking error
is given by:
Here, impacts the dynamics of , influencing the attitude control process.
Combining (
2) and (
4), the differential control equation for the attitude tracking angular velocity error
is as follows:
where
,
.
The upper bound of the uncertain/disturbance part of the system can be determined by the following expression:
in which
,
, and
are unknown positive constants, which satisfy that
is strictly greater than the constant term of
,
is strictly greater than the coefficient of the term
of
, and
is strictly greater than the coefficient of the term
of
.
The primary objective of this study is to design a robust control law
(as illustrated in
Figure 1b), ensuring that both the attitude tracking error
and the angular velocity error
converge asymptotically to zero under disturbance-free conditions, or within a small bounded region when subject to external disturbances. Critically, the predefined stabilization time is designed to be invariant with respect to the initial system states, affirming the control’s robustness across varying operational scenarios.
Before proceeding to the main results of this paper, some preliminaries need to be presented.
4. Main Results
In this section, we propose a novel predetermined time attitude tracking control scheme for rigid spacecraft. To facilitate the controller design, we introduce a sliding mode variable:
where the auxiliary variable
and relevant parameters are defined according to Lemma 1. The controller development process is based on dynamic prediction techniques, and can be explained as follows.
Step 1. Dynamic prediction model for the system input-output.
Design the controller using a general form of optimization method to achieve asymptotic stability of the closed-loop system. Ensure that the spacecraft output
of the lower-triangular system converges to the origin in an optimal way according to the following performance index:
where
is the prediction period.
The above problem will be solved by a variable prediction period
, where
, and the update formula for
L (as shown in
Figure 2a) is given by:
where
is the update learning rate. The steepest descent update rule is transformed into a hidden Lipschitz nonlinear hybrid optimization form, ensuring that the control process achieves global optimization rather than merely local improvements. This transformation is critical for enhancing the robustness and effectiveness of the control strategy across all operational regimes. The optimization approach is equivalent to the adaptive change rate (
17), which dynamically adjusts based on the evolving system conditions, further ensuring the comprehensive applicability and efficiency of the control method [
39,
40,
41].
Based on the configuration of system (
15), the future output
within the prediction horizon (i.e.,
) is predicted using the following Taylor series expansion:
where
,
,
,
denotes the control sequence, and
is the defined constants related to the kinematic-dynamic properties of the spacecraft.
By the form of (
18), the performance index for prediction (
15) can be obtained in the following form:
where
,
,
.
Taking the partial derivative of with respect to U yields . In this paper, is a positive definite matrix. By setting and ensuring that , the optimal control sequence is obtained as .
Assuming that
and
have such forms, the simplified form of the control law can be obtained as follows:
where
,
and
are constants that depend only on
r and
n. The above derivation still holds when the prediction horizon is variable or state dependent [
42].
For the convenience of the proof of the predictive control described in this paper, let
, and the system (
15)–(
20) can be constructed into a more compact expression as follows:
where
I denotes the identity matrix,
,
,
.
According to the relevant theory of optimal control [
40,
41], the stability of the nominal system (
19) depends only on the control order
r. By choosing the control order
r appropriately, there exists a positive definite and symmetric matrix
such that
. In this paper, the update of the prediction horizon is presented in the following adaptive form:
where
and
are adjustable parameters that satisfy
,
. The speediest update method can be shown to be consistent with adaptive characteristics through (
22).
Considering the computational complexity limitations of the dynamics predictive controller for rigid spacecraft, the optimization order
is set as proposed in (
20). Higher-order residuals and remaining parameters are adjusted using the adaptive parameter adjustment described in Equation (
22). Consequently, the following can be derived:
where
,
and
.
The
represent the unknown estimation laws defined in (
22), satisfying the condition that
.
Theorem 1. For the attitude control system (4) and (5) of a rigid spacecraft, the closed-loop system is globally asymptotically stable under dynamic predictive control in the form of (23). The dynamic prediction period is determined by (22), and thus, the spacecraft’s attitudes error and attitude angular velocity error exist. Proof. The closed-loop system discussed in this proof section satisfies the following conditions:
It is evident that all states in the closed-loop system are uniformly bounded and that the system state converges to the origin, i.e., .
According to Lemma 3, for each function , there exist smooth scalar functions and such that .
Considering that
v is an unknown variable that depends on
, which needs to be evaluated in the controller (as shown in (
23)), the following definition is provided for the convenience of the following proof.
Let
and we have
, where
and
are known constants. We define
,
, where
,
[
42].
The domain of the closed-loop system under consideration in this proof section is defined as , where is a positive constant for any initial state that satisfies the following domain condition: .
We construct a Lyapunov function as , where , is an estimate of v, and is a design parameter. Based on the above definitions, we obtain the following constraint: , , s.t..
We can differentiate
with respect to the form of system (
21):
First, we estimate the last two terms on the right side of (
24).
When
holds, by the mean value lemma, we have:
Under the condition that
, where in
. Thus, we can obtain:
where
is a known constant. Thus, we can obtain:
By constructing the proof through the selection of adjustable parameters that satisfy the Hurwitz condition, it can be shown that
P is a positive definite matrix, which yields the following result:
where
.
Combining (
26)–(
28), we can transform (
24) to:
Furthermore, it can be derived that:
Let
and
. Taking (
23) into (
30) yields that domain condition holds.
For the convenience of the proof, suppose that . When , there exists with such that and (if , must be ensured). To simplify the subsequent proof, we assume that the Lyapunov function has the form and satisfies .
Through the aforementioned simplification, the form of the domain is defined as
,
. There exists
such that the Lyapunov function satisfies
. Hence, the following relational equation necessarily exists as follows:
holds for
;
;
holds for
where
(as shown in
Figure 2b).
Therefore, we can conclude that:
As , we have . In fact, z, , and L are uniformly bounded, i.e., .
Let ; we can deduce that . Defining , is uniformly bounded with F. According to Lemma 5, we can obtain that and . This can prove that the proposed update method for dynamic predictive control is locally stable and converges. Hence, the proof of Theorem 1 is completed. □
Remark 1. The control laws designed for the unknown terms that the system is subjected to are denoted by and in (22) in this study. Step 2. Design of a control law for u.
Based on the sliding mode variable
s and the dynamic prediction form of the control torque defined in (
23), the control law for spacecraft control torque is defined as follows:
where
is the predefined-time constant for system stability, and
,
, and
are predefined positive constants that satisfy the definition in (
6).
Theorem 2. Consider the spacecraft described by (1) and (2) and assume that the initial relative attitude does not contain singularities. If the control torque is given by (23) and (32), then the origin of the attitude tracking error system described by (4) and (5) is predefined time stable. Proof. Taking the time derivative of the proposed sliding mode variable
s:
Substituting the dynamic model of the system (
5) into (
33), we obtain:
Substituting (
34) into (
32) gives:
Based on the definition in (
6), there exists
, and we can derive from (
35) that:
Consider the following Lyapunov function candidate
. Taking the time derivative of
yields:
According to Lemma 1, in the convergence phase, the system is predefined time stable, and the predefined time is .
Once the tracking trajectory of the system undergoes a phase of stabilization converging to the origin, i.e.,
when
, by defining the variables, we can obtain:
where
is also the predefined-time constant for system stability.
Consider the following Lyapunov function candidate
. Taking the time derivative of
yields:
Similarly, according to Lemma 1, the system converges stably to the origin in the predefined time
. Moreover, the system under the control of (
23) and (
32) converges stably to the origin within the predefined-time constant
. Hence, the proof of Theorem 2 is completed. □
Therefore, from Theorems 1 and 2, we conclude that and converges to the origin within the predefined time.
Remark 2. The controller designed in this paper relaxes the torque constraint required by general dynamic prediction techniques. Moreover, the upper bound of the predefined-time constant for the closed-loop system in the results is independent of the initial conditions, making the parameter selection of the controller more convenient.
Remark 3. Considering that the sign function will lead to system chattering, the in the proposed controller (32) is modified to the following form:where is a small positive constant. Remark 4. The derivatives of the Lyapunov function in this study are scaled exclusively in the presence of system uncertainty. In contrast, other control schemes undergo multiple scaling steps in their derivation process. As a result, the proposed upper bound on the convergence time presented in this paper is notably less conservative.
5. Simulation Results
In this section, we validate the efficacy of the proposed control scheme through numerical simulation. The nominal inertia matrix of the spacecraft is defined with the following parameters: kg· m2. In the simulation, 30% of the inertia matrix is given as a deviation, which is represented as the uncertain part of the inertia matrix . The spacecraft is equipped with MWs drive, which is a rotating body fixed inside the spacecraft, with kg· m2. The reference attitude is . The control parameters are chosen as , , , , , , , , , , and . The saturation problem of the actuator has been considered during the system simulation, and the torque of the actuator is limit to N· m.
Remark 5. The control parameters in the proposed predefined-time control strategy are critical for achieving optimal performance. Gain parameters (k, α, β) influence the convergence speed and robustness of the control system. Typically, higher values of α and β enhance the system’s ability to reject disturbances and handle uncertainties, but they may also increase control effort and cause chattering. The parameter k should be chosen to balance fast convergence and system stability. Time constants (, , and ) define the predefined stabilization time. Smaller values lead to faster convergence but require more precise control actions. A moderate starting value is recommended, which can then be fine-tuned based on the desired system performance and response. Disturbance parameters (, , and ) estimate the upper bounds of disturbances and uncertainties. Accurate estimation is crucial for maintaining robustness. Initial values can be based on prior system knowledge and adjusted according to observed disturbance levels during operation. Adaptive parameters (, ) control the adaptation rate of the dynamic predictive control scheme. Higher values result in quicker adaptation to changing conditions but may risk instability if set too high. By adjusting the corresponding control parameters, optimal performance can be achieved under different spacecraft configurations and mission requirements.
First, it is assumed that there are no disturbance and measurement noises (define this as Case A). In this case, the initial conditions of the system are
and
deg/s. The simulation results are shown in
Figure 3. It can be observed that the attitude tracking errors of the spacecraft in Euler angles with a 3−1−2 rotation sequence (i.e.,
,
and
) are less than 0.002 degrees when
, which indicates that the proposed control scheme in this paper can achieve perfect spacecraft attitude tracking.
Next, the case with disturbances is addressed (define this as Case B). In this case, the elements of the initial attitude
and initial angular velocity
deg/s are given, respectively. To investigate the robustness performance of the closed-loop system, an external disturbance
,
N· m is given. The results are shown in
Figure 4, indicating that the attitude tracking errors of in Euler angles are less than 0.008 degree at the steady-state stage.
Furthermore, we discuss the case where measurement noise is present (define this as Case C). In this case, the initial values of the attitude
and angular velocity
are randomly selected from the range
and
deg/s, respectively. The measurement noise signals defined in the controllers during the simulation calculations are given according to An-Min Zou et al.’s method [
37]. The results in
Figure 5 show that at the steady state, the attitude tracking errors in the Euler angles are less than 0.01 degrees. The excellent performance of the proposed control scheme is further evidenced by the root mean square error (RMES) of the attitude tracking presented in
Table 1.
Finally, a comprehensive comparison is conducted between the proposed predefined-time dynamic predictive controller, traditional PID control (as shown in (
41)), and the existing fixed-time control (as shown in (
42)) methods under the same conditions:
where
,
, and
are the proportional, integral, and derivative gains, respectively. The PID controller is widely used due to its simplicity and effectiveness, but it has significant limitations when dealing with uncertainties, bounded actuators, and external disturbances:
where
,
(
) are constants, and
is a diagonal matrix. This fixed-time method provides a guaranteed convergence time but faces challenges in parameter tuning and robustness under high levels of uncertainties and disturbances.
As shown in
Figure 6, even with simultaneous uncertain factors and measurement noise affecting the spacecraft system, the proposed predefined-time dynamic predictive method exhibits superior performance in terms of spacecraft attitude angle tracking, angular velocity response, and control torque jitter compared to both PID and fixed-time controllers. The PID controller shows significant performance degradation under uncertainties and disturbances. The fixed-time controller offers improvements but still lacks the flexibility and robustness of the proposed method.
In conclusion, the proposed control strategy demonstrates strong robustness to existing uncertainties and disturbances, achieving the control goal of attitude tracking more effectively than traditional PID and fixed-time control methods. This highlights the practical benefits of the proposed method for advanced spacecraft attitude control applications.