1. Introduction
In recent years, the quest to stabilize flexible spacecraft has intensified, sparked by their advantages for space missions and evidenced by substantial research (e.g., [
1,
2]). These spacecraft exhibit complex dynamics, with a mix of flexible and rigid behaviors, and the critical system parameters that influence these dynamics are often hard to measure precisely. This poses significant challenges in control strategy design, affecting system performance. Fuzzy system models have emerged as a robust solution to handle such nonlinear complexities, effectively approximating nonlinear functions with minimal error (e.g., [
3,
4]). Leveraging the T-S fuzzy model framework, recent studies have developed advanced control laws that address the dynamic intricacies of nonlinear systems. For example, the problem of memory-event-triggered fault detection of networked IT2 T-S fuzzy systems was considered in [
5], and the authors in [
6] achieve decentralized adaptive event-triggered
filtering approaches for a class of networked nonlinear interconnected systems. Moreover, research detailed in [
7] applied an event-triggered
control to a truck-trailer model using T-S fuzzy models, while [
8] used advanced Lyapunov methods to design a controller for a type-2 T-S fuzzy system under external disturbances. These efforts represent significant strides in the control of complex, flexible spacecraft systems.
Furthermore, it is crucial to highlight the pivotal role that actuators play in control systems, as their failures can lead to suboptimal performance or even system instability [
9,
10]. This is particularly significant in the context of space vehicle applications, where reliability is of paramount importance. The research delineated in [
11] addressed the intricate issue of adaptive fuzzy control within a finite-time scope for dynamical systems that are not explicitly modeled, with a particular focus on accommodating actuator faults. Similarly, the study in [
12] presented an advanced actuator fault-tolerant control method employing adaptive barrier fast terminal sliding mode techniques for space vehicles. These works underscore the profound impact of actuator faults on the efficacy of control systems. Actuator malfunctions can markedly reduce the operational efficiency of space vehicles, with the potential for such deficits to escalate to catastrophic outcomes, emphasizing the criticality of robust fault-tolerant mechanisms in space systems.
Consequently, there has been a growing emphasis on the design of reliable control systems to mitigate actuator faults, ensuring high reliability and robustness, not only in the realm of space vehicles but also in industrial systems. Nonetheless, it is essential to highlight that numerous established dependable control methodologies operate with the presumption that faults are deterministic in nature [
13,
14]. In practical scenarios, faults can be both deterministic and stochastic, arising from factors such as aging or damage to the control components. To address this complexity, recent research efforts have explored various strategies to handle actuator faults with stochastic characteristics. As an illustration, in the paper referenced as [
15], scholars investigated the challenge of output feedback
control within the context of nonlinear spatially distributed systems that were impacted by actuator faults following a Markovian jumping pattern. It is imperative to acknowledge the pivotal role that the Markov jump process plays in the accurate representation of actuator faults. This stochastic framework is instrumental in characterizing the unpredictable nature of such faults, thereby providing a robust foundation for subsequent analyses. Consequently, an extensive body of scholarly work has been dedicated to studying various aspects and implications of this phenomenon within the context of control systems [
16,
17,
18,
19]. To be concrete, a robust state feedback control strategy was formulated for T-S fuzzy systems affected by sensor multiplicative faults, utilizing the dynamic parallel distributed compensation technique [
20]. In another study [
21], actuator failures were modeled as random variables subject to a Markov process. The primary innovation here was the proposal of a novel T-S fuzzy controller tailored for spacecraft dealing with stochastic actuator faults. The authors in [
22] attempted to address the challenges posed by actuator faults by introducing a constraint condition, namely achieving
performance. In a different approach [
23], a novel stabilization approach was utilized to address the challenges posed by actuator faults. This method featured a novel iterative linearization algorithm for designing controller gains. Moreover, in the paper by [
24], the problem of finite-time control in a nonlinear flexible spacecraft system affected by stochastic faults was addressed. However, it is important to note that these random faults were assumed to be homogeneous. In operational environments, the failure dynamics of actuators are inherently variable, manifesting time-dependent characteristics influenced by an array of factors including material degradation, ambient environmental conditions, and the intensity of operational demands. These variables engender fluctuating failure rates, rendering the supposition of time-invariant failure mechanisms unsuitable for authentic representation. Consequently, the premise of homogeneous actuator failures is a simplification that does not accurately encapsulate the complexities of real-world operational scenarios (see [
25,
26]). Consequently, there is a pressing need for research that investigates nonhomogeneous actuator faults in conjunction with stochastic variables, as this aligns more closely with the practical challenges encountered in control systems.
Conversely, the conventional
control approach is commonly utilized to examine the desired behavior of controlled dynamics and explore the asymptotic characteristics of system trajectories, typically under the assumption of an infinite time interval, as governed by Lyapunov stability theory. Nonetheless, there exist practical scenarios where the central focus pertains to the performance of dynamical systems within a specified finite time duration. As an example, circumstances may emerge where the influence of external disturbances necessitates the assurance that system states remain within permissible limits throughout this finite time span [
27,
28]. In such scenarios, the notion of finite-time boundedness assumes significant importance. Over recent decades, there has been a noticeable surge in interest regarding this concept, substantiated by a substantial body of research and the corresponding references [
29,
30,
31]. In these studies, the free-weighting matrix method has been commonly employed, often necessitating the introduction of slack matrices. However, it is important to recognize that the adjustability of these introduced matrices is inherently limited, and there is room for improvement in this regard.
With the above analysis, we delve into the subject of finite-time fuzzy fault-tolerant control (FTC) for spacecraft system with actuator faults in this paper. What distinguishes our research from the prior literature, notably [
22], is our consideration of actuator faults subject to a nonhomogeneous Markov distribution, a choice that aligns more closely with real-world scenarios. Recognizing that closed-loop control systems involve nonhomogeneous random data, we develop a nonhomogeneous Lyapunov functional that comprehensively incorporates all available information from the closed-loop system. It is worth noting that the nonhomogeneous Markov distribution can be simplified to the traditional distribution by setting the derivative of nonfeasible matrices to zero, thereby adding an extra layer of versatility to our approach. To augment the adaptability of our control strategy, we introduce a set of flexible parameters. This innovative controller design approach permits the independent selection of these parameters, offering a higher degree of freedom and customization when compared to existing methods. Importantly, our flexible parameters are predefined, distinct from the unknown slack matrices encountered in prior works [
22,
24]. When compared to existing studies such as [
22,
24], our control approach exhibits superior practicality and performance.
Notations: The symbols and represent, respectively, the vector space of dimension n and the space of matrices. The norm denoted by corresponds to the Euclidean norm for vectors, and denotes the linear space comprising square integrable functions defined over the time interval , where . For a matrix A, its transpose is denoted as , and its inverse as . The maximum and minimum eigenvalues of matrix A are symbolized as and , respectively. The sum of A and its transpose is succinctly represented as or , and denotes an n-dimensional diagonal matrix. For a vector , represents its skew-symmetric matrix, and refers to the mathematical expectation of a random variable. The symbol ⨂ is employed to denote the Kronecker product for matrices. The symbol ∗ signifies the symmetric term within a matrix. Unless explicitly mentioned, matrices are assumed to possess compatible dimensions.
2. System Description and Preliminaries
Consider a dynamic model for a flexible spacecraft, taking into account the following variables and parameters in
Table 1.
To represent the spacecraft’s rotational state, we introduce quaternion notation. Set
, and define the quaternion vector as
, with the constraint that
. The dynamic model that governs the behavior of this spacecraft system can be represented as follows:
where for
,
and
as the diagonal matrix
. Here,
represents the damping ratio for the
ith mode,
represents the modal frequency for the
ith mode, and
represents the total number of modal frequencies. The control input to the system is denoted as
, and
represents an external disturbance variable.
Remark 1. In Equations (1)–(4), the incorporation of system inertia, encompassing the dynamics of both the rigid body and its elastic attachments, alongside the incorporation of coupling matrices and the variables representing elastic coordinates, fosters a holistic and all-encompassing modeling approach. This comprehensive formulation of the spacecraft’s dynamics lays a robust groundwork for the rigorous evaluation and validation of the proposed control mechanisms within our study. As
,
, and
are interconnected, a more convenient approach is to define a consolidated disturbance term as follows:
Combining (
3) and (
4), it yields
assuming that
is satisfied for
.
Consider the state vector,
both of which are the controlled outputs. By combining Equations (
1)–(
5), we can deduce the subsequent dynamic system by employing an “IF-THEN” fuzzy rule.
Rule
i: IF
,
,...,
, THEN
where
, we have
, and
are T-S fuzzy sets. Within this framework, the system matrices
,
, and
can be expressed as follows:
Hence, the T-S fuzzy system can be delineated as
where
The sensor network acquires its signal from the system’s output, with an output feedback approach employed for controller design. In practical systems, the output signal, denoted as , is routinely sampled at predetermined time points such as utilizing a sampler. The discretely sampled output signal is subsequently directed to the sensor network’s input interface. When considering a specific sampling instant , it is imperative to acknowledge that the value of remains constant throughout the time span . This signifies that the signal received by the sensors within this interval (i.e., for ) is consistently represented by .
Traditionally, the conventional method for detecting output signals involved deploying a single sensor mode within the model. In this setup, each measurement output was individually sampled by dedicated sensor units and then transmitted collectively, aggregated into a unified packet for storage. However, it is crucial to recognize that this assumption may not hold in various real-world scenarios. A sensor network typically comprises a multitude of distributed sensing devices operating in coordination. The measurement outputs available to a particular sensor node, denoted as node
i, do not originate solely from that node itself. Instead, these outputs result from a combination of measurements from node
i and those from its neighboring nodes within the network of sensors. The network of sensors’ topology is defined as follows: The directed graph denoted as
characterizes the spatial configuration of sensors. Here, the set
functions as an indexing set encompassing
m sensor nodes,
denotes the collection of edges that establish connections between pairs of sensor nodes, and each edge within
is represented by an ordered pair
. Furthermore,
corresponds to the matrix containing weights associated with the connections in the network. Notably, the elements of the adjacency matrix are constrained to be positive for edges that exist in the graph (i.e.,
if
). Additionally, for nodes
, the matrix
has its diagonal elements
, which can be calculated as
The ideal control input is
the
sensor output is
and
, where
represents the output matrix of the
-th sensor, and
fulfills the following condition:
In this paper, we consider a scenario where the sampling interval is aperiodic but still bounded. This implies that for all consecutive sampling instances, we have . Let us define a function to represent the time elapsed since the last sampling instant . This function satisfies the conditions and for all time points except at the specific sampling instance .
Considering nonhomogeneous actuator failure, which can be denoted as
where the vector
sees (
8),
represents a bounded actuator fault. For constant
,
can be represented as
. Furthermore, we define
and
as the fault matrices, where
for
. These matrices are associated with the fault introduced into the system.
Remark 2. It is noted that the control gain K in (8) is a parameter designed to maintain system performance in the event of actuator faults and needs to be adjusted for potential fault impacts. The fault matrices in (9) describe how faults alter system dynamics, and the control gain K must be sufficiently robust to accommodate these changes as described by , ensuring system stability and performance criteria under various fault conditions. Hence, the design of K must take into account the characteristics of to achieve effective control. Considering the timeliness,
presents the nonhomogenous Markov process satisfying
and transition rates are given as follows:
where
and
.
In this paper, the nonhomogeneous transition rate matrix
is confined within a polytope defined by several vertices. This confinement is characterized by the following expression:
where
are given matrices that represent the vertices of the polytope, and
where we have
as a parameter vector. It is important to note that the rate of change of this vector, denoted as
, is presumed to be both bounded and known. Specifically,
is assumed to adhere to the following condition:
and
is known. It is worth noting that the condition
leads to the following implication:
Consequently, the upper bound on
can be expressed as:
Remark 3. The characterization of the nonhomogeneous transition rate matrix is considered to adhere to a polytopic representation. As an illustration, let us assume a nonhomogeneous transition rate matrix with two dimensions and two vertices, which can be expressed as follows:
Hence, the conditions satisfied by
and
are as follows:
Let us examine the following transition rates:
In this scenario, we define
and
as random functions generated using MATLAB, and they satisfy the constraint
with the additional conditions
and
. Under (
12) and (
13),
and
are depicted in
Figure 1. It is evident that the transition rates exhibit time variation. Furthermore, the vertices of the transition rate matrix are determined as follows:
Thus, we have
To simplify the notation, we define
,
, and the same for
with
being either 0 or 1. Additionally, in Equation (
9), we establish the bounds
, where
and
are the lower and upper limits of
, respectively. By introducing the matrices
,
,
, and
, we can express
as a combination of
and
scaled by
, where
with
.
By incorporating (
6) and (
9), we introduce the following fuzzy FTC strategy. Under the sampled-data input, this controller is designed as follows:
Rule
i: IF
,
,...,
, THEN
where
is the desired controller.
Combining (
14) and (
6), we can derive the following closed-loop system:
Furthermore, by recognizing that
we can rewrite (
15) as
where
and
In this paper, our primary aim is to develop a fuzzy fault-tolerant controller with the goal of satisfying the following performance criteria:
(i) Provided with scalars
, and a positive definite matrix
, we can conclude that system (
16) exhibits finite-time stability and data boundedness (FTSDB) concerning the parameters
if the following conditions hold:
(ii) For
and
, if the following condition is satisfied, then (
16) meets the
performance:
By defining w(s) and imposing the constraint , we devise a methodology that capitalizes on the recorded output from an individual sensor and its adjacent sensors. This approach yields a novel condition for maintaining system stability with the added constraint of performance. Additionally, we derive a method for designing controller gains to meet these requirements.
3. Main Result
To date, the Lyapunov stability analysis method has been extensively applied in the design of controllers for systems with delays, owing to its broad range of potential applications. In line with this approach, the primary findings of this paper have also been derived. To derive the FTC method, we will introduce two helpful preconditions.
Lemma 1 ([32]). Given constants , and variable x: , one has Lemma 2 ([22]). Given P, Q, and R, if there exists , then for , the following condition holds: Theorem 1. Given scalar values , , , and , system (16) is exponential stability in the mean-square sense and can be characterized as FTSDB concerning the parameters and satisfies the performance criterion if there exist scalar values , , , , , , , and matrices , , such thatwhere and Proof. Let us formulate the nonhomogeneous Lyapunov functional as follows:
where
Let
be the weak infinitesimal generator. Along the nonhomogenous Markov distribution, we have
When calculating
along the closed-loop system given in (
16), we obtain the following result:
Building upon Lemma 1, we derive a lower bound for
as follows:
For
,
and any free parameters
,
, one has
where
Together with
and
, it yields
Let
. One has
where
.
Note that
We have
Consider that
Thus, we obtain
It is important to note that (which belongs to the interval ) can be considered either as an upper or lower bound. Thus, can lead to the conclusion that .
Due to the presence of the nonhomogeneous terms
and
, it should be noted that the inequality in Equation (
18) is infinite and nonlinear in nature. Define
. Since
, if the matrices
and
satisfy the conditions (
18) and (
19), one has
Considering the inequality
, we have
Considering the interval
, one has
Moving forward, we will establish an upper bound for
, taking into account the continuity property
Considering that the initial sampling instant
may not coincide with 0, let us address
by accounting for the following: When
t lies within the interval
, we set
. Consequently, based on the following condition, we can derive the nonlinear spacecraft dynamics
which means
Based the norm space, it yields
It is clear that there are constants
,
, and
such that
For
, one has
which can be further deduced as
Set
in (
44), and let
be known. It is not difficult to see that
It is obtained from (
42) that
which can further deduce
Let
,
. It is easy to see that
Due to
, one has
Since
, it yields
Under the initial condition, it yields
From
, it yields
This completes the proof. □
Remark 4. In this research, the inclusion of the term within holds pivotal importance as it facilitates the integration of information pertaining to the nonhomogeneous Markov process. This augmentation significantly enhances the realism of our findings when compared to the existing literature. It is noteworthy that the transformation presented in Equation (40) plays a critical role in effectively managing the derivative of . Additionally, the introduction of boundary values denoted as serves to transform the non-feasible derivative into more tractable finite values. It is noted that Theorem 1 is not easy to settle due to the presence of nonlinearity. In the following, a feasible condition to reach the finite-time controller will be given.
Theorem 2. For any , and positive constants α, γ, , , , , , T, any matrix , we can establish the existence of positive scalars , , , , , , , , , and any matrices , , for such thatwhere , , symmetric matrix and We can conclude that the closed-loop fuzzy system described by (
16) exhibits exponential stability in the mean-square sense and stochastic FTSDB behavior with respect to the parameters
and simultaneously meets the criteria for
performance.
Proof. To begin, we multiply both sides of inequalities (
18) and (
19) by
and its transpose, respectively. Then, we have
,
,
,
,
,
,
. Using these, we can establish inequality Equation (
45) based on the inequalities Equations (
18) and (
19). Now, we will demonstrate that inequalities (
46)–(
50) imply inequality (
20). Let
and
. It is evident that
Then, it follows that
based on Schur’s complement theorem and inequalities (
46) and (
47). Consequently, we can deduce that
. Similarly, we can establish that
and
, which implies that
,
,
, and
. In other words, the formula (
50) holds. This completes the proof. □
In accordance with the conventions and terminology commonly employed in the academic literature, we shall elucidate the process of determining the parameters for the controller through the presentation of the Algorithm 1.
Remark 5. In this work, we build upon the traditional use of static slack matrices and , integrating them with dynamic variables and to establish a pioneering control design methodology. This dual approach is aimed at achieving two main objectives:
(i) It establishes an adaptive framework that enables precision calibration to match the specific dynamic traits of the system. This level of adaptability is crucial for advancing the effectiveness of the controller beyond what is reported in [22,24]. Through meticulous adjustment of and , we anticipate significant improvements in system robustness and performance. We propose employing sophisticated optimization techniques, such as genetic algorithms and the Nelder–CMead simplex method, to determine the most advantageous parameter configurations, as suggested by [33,34]. (ii) Our methodological proposition is mindful of the computational burden that large-scale matrix calculations entail. Strategic manipulation of the free matrices and effectively reduces the matrix dimensionality, thus alleviating computational loads. This aspect is especially beneficial for real-time applications where computational efficiency is paramount.
Algorithm 1 Fuzzy fault-tolerant controller solving algorithm [24]. |
Require: T-S fuzzy system matrix parameters , , ; |
nonhomogeneous fault matrix , ; |
Sensor network parameters , ; |
FTSDB parameters ; |
Ensure: ; |
1: Define matrix variable (, Q, R, W, U, D, , , , , ; |
2: for do |
3: Given stochastic fault matrix ; |
4: Let ; |
5: loop; |
6: for do |
7: for do |
8: Construction matrix , , and ; |
9: Establish the constraint Equation (45) of Theorem 2 |
10: end for |
11: end for |
12: while do |
13: Let ; |
14: goto loop. |
15: end while |
16: end for |
17: Establish the constraints Equations (46)–(50) of Theorem 2; |
18: Solve the constraints; |
19: Return . |
Remark 6. In this study, it is imperative that remains finite to establish a set of linear conditions. In cases where any of these conditions diverge to infinity, the only viable recourse is to set , resulting in the constructed Lyapunov function becoming homogeneous. It is noteworthy that a homogeneous Markov process represents a distinct case within the framework of the proposed nonhomogeneous model introduced in this paper. Specifically, when , the nonhomogeneous model effectively reduces to a homogeneous one in existing works [22,24]. In such scenarios, the findings and results outlined in this paper remain applicable and general. Remark 7. In this paper, through the T-S fuzzy framework, the nonlinear flexible spacecraft system is transformed into the T-S fuzzy model represented by Equation (6). A Lyapunov function tailored for the T-S fuzzy model (6) is then constructed, which addresses the fault-tolerant control issues inherent in nonlinear flexible spacecraft systems. Consequently, the functional approach detailed herein can be extended and applied to real-world systems that can be accurately described by model (6), exemplifying the versatility and broad applicability of the methods developed in this study to fields such as robotics or autonomous vehicles [35,36]. Remark 8. It should be acknowledged that in addition to the T-S fuzzy model framework, an array of alternative modeling methodologies can be judiciously utilized to mitigate actuator faults within nonlinear systems influenced by nonhomogeneous Markov processes. Techniques such as the backstepping method, neural-network-based control paradigms, and reinforcement learning algorithms are documented for their efficacy in this domain [37,38]. The intrinsic benefits and limitations of these methods need to be carefully assessed against the backdrop of the system’s unique attributes. The choice of a particular modeling technique is contingent upon system-specific demands, encompassing the extent of nonlinearity, computational resource allocation, and the requisite robustness and fault tolerance levels within the control architecture. As for empirically validating the proposed methodology, methods such as hardware-in-the-loop simulations or physical experiments with spacecraft hardware deserve further work. Remark 9. This paper primarily tackles the theoretical challenges associated with fault-tolerant control mechanisms within aerospace systems through a polytopic set theory approach. While the theoretical framework laid down provides a robust foundation for addressing system faults, it stops short of delving into practical applications. Future research will pivot towards translating these theoretical insights into tangible applications within the aerospace sector. This transition will necessitate a thorough assessment of the actual capabilities and limitations of spacecraft hardware and sensors to ensure the developed control strategies are not only theoretically sound but also practically viable. Subsequent studies will aim to refine the control paradigm to align precisely with the operational demands and performance benchmarks set by current spacecraft hardware configurations and sensor technology integrations.
Remark 10. The methodology proposed in this research relies heavily on the precision and timely response of actuators. In practice, however, this control approach may present compatibility issues with existing spacecraft systems. Additionally, the harsh space environment, characterized by radiation, extreme temperature fluctuations, or space debris, could lead to equipment failure or degraded performance. To mitigate these risks, it is prudent to refer to existing literature [39,40] for continuous testing and validation, along with the implementation of redundant sensor and actuator designs, to bolster system reliability and resilience. Such strategies can effectively ensure that theoretical research is translated into safe and reliable applications in actual space missions. 4. Numerical Example
Consider the flexible spacecraft model (
1), sourced from [
22]. The primary energy associated with elastic vibration predominantly resides in low-frequency modes. Consequently, it becomes necessary to account for only the first four elastic modes. The respective natural frequency vector for these modes is denoted as
rad/s, with corresponding damping ratios
. Then,
and
are specified as
and
is given as
The other condition are the same as Ref. [
24], which are omitted here.
Selecting three operating points, defining and setting , , and , the flexible spacecraft can be constructed as
Rule 1: IF
,
, THEN
Rule 2: IF
,
, THEN
Rule 3: IF
,
, THEN
The selection of lower and upper membership functions is as follows: when
j takes on values of
if
where
and
. The corresponding curves are depicted in
Figure 2 and
Figure 3.
Considering the previously mentioned three operating points, we can determine the values of the matrices
,
, and
for
in the presence of the external disturbance
. These matrices are illustrated in
Figure 4. In order to stabilize the flexible spacecraft system, matrix
is provided as follows:
The sensor output matrices are given as follows: , , and . For and , consider the following three cases.
Case II—flawed operational state with partial loss of functionality:
Case III—bias faulty mode:
The transition rate matrix for the nonhomogeneous Markov process is as follows:
where
Setting
,
,
, and
, and focusing solely on Case I, where
and
, we utilize Theorem 2 to determine that the maximum allowable value for the upper bound of
is
. When considering all three cases under the same theorem, this upper bound is refined to
. Notably, this finding surpasses the upper bound of
previously established in [
24]. From the above analysis, we draw three key conclusions:
C1: The comparison between
and the value reported in [
24] (
) clearly demonstrates the superiority in the current study, signifying a notable enhancement in the system’s performance.
C2: Observing that the upper bound remains at and across different faults, it is evident that the system’s performance is optimized as fault tolerance is reduced. Thus, to achieve better system performance, one must tolerate smaller faults, and vice versa.
C3: It is important to highlight that the homogeneous case represents a special instance of the nonhomogeneous case as given in (
52), with the latter being a more encompassing and general representation.
The three pivotal conclusions drawn highlight that the proposed method not only exhibits a more generalized framework but also outperforms existing models in terms of system performance.
Further, based on the condition specified in Equation (
53), Theorem 2 provides the feedback gains as follows:
When initializing the system with
, we can observe the state responses of the closed-loop fuzzy system in
Figure 5 and
Figure 6. These figures serve to highlight the effectiveness of the control design method, particularly when dealing with a nonhomogeneous Markovian process.