Next Article in Journal
Is the Increased Transmissibility of SARS-CoV-2 Variants Driven by within or Outside-Host Processes?
Next Article in Special Issue
Stability Analysis on the Moon’s Rotation in a Perturbed Binary Asteroid
Previous Article in Journal
Sizing and Sitting of Battery Energy Storage Systems in Distribution Networks with Transient Stability Consideration
Previous Article in Special Issue
Asteroids and Their Mathematical Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simultaneous-Fault Diagnosis of Satellite Power System Based on Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set

1
Institute of Reliability Engineering, Beihang University, Beijing 100191, China
2
Science & Technology on Reliability & Environmental Engineering Laboratory, Beijing 100191, China
3
School of Reliability and Systems Engineering, Beihang University, Beijing 100191, China
4
Beijing Institute of Radio Metrology and Measurement, China Aerospace Science and Industry Corporation Limited, Beijing 100039, China
5
Shanghai Institute of Satellite Engineering, China Aerospace Science and Technology Corporation, Shanghai 201109, China
6
Marine Design and Research Institute of China, China State Shipbuilding Corporation Limited, Shanghai 200011, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(19), 3414; https://doi.org/10.3390/math10193414
Submission received: 12 August 2022 / Revised: 3 September 2022 / Accepted: 7 September 2022 / Published: 20 September 2022
(This article belongs to the Special Issue Mathematical Problems in Aerospace)

Abstract

:
Due to the increasing complexity of the entire satellite system and the deteriorating orbital environment, multiple independent single faults may occur simultaneously in the satellite power system. However, two stumbling blocks hinder the effective diagnosis of simultaneous-fault, namely, the difficulty of obtaining the simultaneous-fault data and the extremely complicated mapping of the simultaneous-fault modes to the sensor data. To tackle the challenges, a fault diagnosis strategy based on a novel rough set model is proposed. Specifically, a novel rough set model named FNζDTRS by introducing a concise loss function matrix and fuzzy neighborhood relationship is proposed to accurately mine and characterize the relationship between fault and data. Furthermore, an attribute rule-based fault matching strategy is designed without using simultaneous-fault data as training samples. The numerical experiments demonstrate the effectiveness of the FNζDTRS model, and the diagnosis experiments performed on a satellite power system illustrate the superiority of the proposed approach.

1. Introduction

The power system is regarded as the heart of a satellite, whose health management is critical to the on-orbit operation of the entire satellite. The satellite power system is mainly composed of a solar array and battery pack. The solar array exposes in the outer space environment for a long time, and is very vulnerable to external environment intrusion. The battery pack is in a frequent and long-term working state with the periodic operation of the satellite. Therefore, with the increasing probability of space junk collisions, intense radiation of space particles, and striking temperature differences in space, the satellite power system may have multiple independent single faults occurring at the same time, which is called simultaneous-fault [1]. Accurate fault diagnosis is the basis for the health management of a satellite. At present, the research on the diagnosis of single-fault has achieved great success [2,3]. However, as satellites become more complex in their functional composition and longer in their mission time, the mode of simultaneous-fault has become the key factor affecting the normal on-orbit operation of satellites, and the risk and influence caused by such a fault mode cannot be ignored, because the development speed and destructive power of such a fault mode are far more than that of a single-fault mode. Therefore, it is necessary to diagnose the simultaneous-fault precisely to make sound decisions to enable satellites to perform their missions smoothly and safely. This is the core motivation of our work, namely, we try to solve the diagnosis problem of simultaneous-fault, which is more complex and more harmful than that of single-fault.
With respect to the diagnosis of simultaneous-fault, there are two major challenges that can be listed as follows: (1) The historical simultaneous-fault data are scarce, which greatly limits the effectiveness of data-driven models; (2) The simultaneous-faults would involve multiple sensors, and the mapping between sensor data and fault modes is complicated, which leads to the uncertainty in the diagnosis process. Therefore, new cognitive methods and further research are needed for simultaneous-fault cognition and diagnosis. These can be considered as the technical motivation of our research.
Regarding the first challenge issue, the absence of historical simultaneous-fault data is a thorny problem that needs to be solved urgently. Unlike the traditional fault diagnosis studies that require all kinds of samples during the training phase [4,5,6], some literature has shown that multi-label classification is expected to achieve simultaneous-fault diagnosis without historical simultaneous-fault data [1,7,8,9,10]. The multi-label classification task focuses on the problem where each training simple is represented by a single instance with a single label, and the task is to yield a model that can predict the proper label sets for unseen instances [11]. Multi-label classification methods can be divided into two categories, one of which is the problem transformation methods, including Binary Relevance [12], Classifier Chains [13], Calibrated Label Ranking [14], and other classical methods; the other includes the algorithm adaptation methods, including multi-label K-nearest neighbor (ML-KNN) [15], multi-label decision tree (ML-DT) [16], etc. However, in the face of complex problems, the above methods cannot effectively deal with the problem of insufficient data, and there is still a need for long-term and in-depth research.
For the second challenge issue, some data mining methods are good solutions. In terms of the cognition of things, rough set theory provides a perspective of knowledge and data fusion. This is the main reason why this paper chooses the rough set model as the basic model. The setting of condition attribute and decision attribute can provide multiple information for the characterization of fault, which is conducive to extracting the mapping information between sensor data and fault modes. Rough set theory initiated by Pawlak [17] provides an authoritative mathematical framework for analyzing and handling ambiguous and uncertain data, which can be used to attribute reduction [18,19,20,21,22], rule extraction [23,24,25,26], and uncertainty reasoning [22,27,28,29]. Among kinds of rough set models, the decision-theoretic rough set (DTRS) model has been proved to be a generalized model of many other rough set models [30,31]. At present, there have been related studies on various decision-theoretic rough set models for fault diagnosis, which have proved that the models can effectively select the fault attributes when the pair of the threshold parameters is set appropriately [30,31]. Nevertheless, how to determine the appropriate threshold parameters is the biggest difficulty in the research and application of DTRS. In our previous work [32], we have presented a single-parameter decision-theoretic rough set (SPDTRS) model by setting only one parameter named compensation coefficient rather than two or six, which facilitates the convenient application of the DTRS model. However, the setting of the compensation coefficient in this model is still not clear enough, and the setting of the loss function matrix is defective. In addition, this model lacks the consideration of uncertain information in data description, which makes it unable to deal with continuous data directly. Therefore, in order to make the rough set model (i.e., SPDTRS) more effective in dealing with the simultaneous-fault problem, we need to carry out more targeted improvement work. The details are as follows.
Motivated by the analyses mentioned above, in this work, we propose a fault matching strategy for simultaneous-fault diagnosis based on a revised DTRS named fuzzy neighborhood ζ-decision-theoretic rough set model (FNζDTRS). Since there is a coupling relationship of fault characteristics between a single-fault and its associated simultaneous-fault, this paper proposes the fault matching strategy based on this principle. The main idea of the proposed strategy is that when an unknown simultaneous-fault occurs, its fault attributes are first selected by the FNζDTRS and then classified according to the correlation between the obtained fault attributes and the fault attributes of each single-fault selected by the FNζDTRS model beforehand. Therefore, the main novelties and contributions of this study can be listed as follows.
(1)
A novel and concise data-driven loss function matrix is designed for DTRS.
(2)
A fuzzy neighborhood ζ-decision-theoretic rough set model is proposed with the help of the fuzzy neighborhood relationship and the proposed loss function matrix, which can deal with hybrid data common in engineering.
(3)
The proposed FNζDTRS model, used for attribute reduction, has a significant advantage in classification accuracy compared with other existing rough sets. This proves that it is more suitable for real fault diagnosis.
(4)
A diagnosis strategy of simultaneous-fault is put forward based on a coupling mapping relationship between single-fault and its associated simultaneous-fault. This ensures that our strategy can handle both single-fault and simultaneous-fault.
(5)
The proposed strategy is successfully applied to the simultaneous-fault diagnosis of the satellite power system and only requires single-fault samples in the training phase, which is highly feasible for practical applications.
The remainder of this paper starts with some preliminaries and related work, then puts forward the presentation of the FNζDTRS model in Section 2 and presents the basic framework of simultaneous-fault diagnosis in Section 3. The effectiveness and superiority of the FNζDTRS model is verified through some numerical experiments in Section 4, and further demonstrated by a comparative analysis with several baseline algorithms for simultaneous-fault diagnosis in Section 5. The paper closes with main conclusions in Section 6.

2. Preliminaries and Related Work

This subsection will review some notions about rough sets that are relevant to the development of our theory.
Definition 1.
(Decision system) A binary group: D S = U , C D can describe a decision system. Among them, U = x 1 , x 2 , , x m is called the universe, which is a finite and nonempty set. D is the set of decision attributes which is a nonempty set. C is the collection of conditional attributes, C D = , D 31. Therefore, the relationship between each element in a decision system can be represented as shown in Figure 1. To better understand the above definition, we describe the above-mentioned elements in combination with the fault diagnosis problem. C represents the parameters output by the sensor or the extracted feature attributes, D represents the category of the failure mode, and U denotes the collected data.
The decision-theoretic rough set (DTRS) presented by Yao et al. [33]. provides a concise semantic interpretation through a loss function matrix. The loss function matrix is described in Table 1.
Consider that X is the subset of samples with the same label d k . The state Q suggests that a related sample defined as x is in X , and the state Q c suggests that x is not in X . The set of actions a P , a B , and a N indicate the classification of x into three regions, which are x P O S X , x B N D X , x N E G X . P O S X denotes the acceptance of the event x X . B N D X denotes the deferment of the event x X , also considers B N D X denotes the non-commitment of the event x X . N E G X denotes the rejection of x X . Furthermore, λ P denotes the loss caused by taking actions ( a P , a B , a N ) while x X . λ N is the loss caused by taking actions ( a P , a B , a N ) while x X .
Consider this scenario: the risk of delaying the execution of the correct action is increased compared to that of the correct action, and both are less than the loss of taking the wrong action, the DTRS model therefore made a reasonable assumption: 0 λ P P λ B P < λ N P and 0 λ N N λ B N < λ P N , which is the basis for generating this rough set model.
Thanks to the above assumption, a pair of threshold parameters is used to define the positive region   P O S X , the boundary region B N D X and the negative region N E G X to construct the DTRS model, which is guided by the Bayesian risk minimization principle and the three-way decision theory. Thus, we have the form of a DTRS model as follows 31:
P O S α , β X = x U | P X | x α ,
B N D α , β X = { x U | β < P X x < α } ,
N E G α , β X = x U | P X | x β .
The following equations represent the relationship between the two threshold parameters and the six loss functions:
α = λ P N λ B N λ P N λ B N + λ B P λ P P ,
β = λ B N λ N N λ B N λ N N + λ N P λ B P .
The key parts of the DTRS model are the loss function or threshold parameter ( α , β ) . To study and employ the DTRS model, an important issue is how to determine these parameters. Inspired by the idea of being data-driven, our previous work proposed a single-parameter decision-theoretic rough set (SPDTRS) model [32] that simplifies the traditional DTRS model. Specifically, the model requires only one parameter to be preset rather than the pair of ( α , β ) or the six parameters in the loss function matrix. However, the solution of employing two truncation functions utilized in the model calculation makes the model relatively complex. Moreover, the interpretability of the loss function matrix in the SPDTRS model is slightly insufficient. The above two disadvantages are the focus of this paper in proposing a new rough set model.

3. Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set

3.1. Granular Computing Based on Fuzzy Neighborhood Relationship

In order to be more applicable to practical problems, rough set models need to be able to handle a hybrid dataset, including continuous and discrete data. Fuzzy relationship and neighborhood relationship are two effective means to deal with the spatial relationship of samples. Their combined form is used by a variety of models [34]. The fuzzy neighborhood relationship can analyze the relationship between the entities in the decision system more precisely. Therefore, to overcome the inability of the SPDTRS model to handle the hybrid dataset, we introduce this fuzzy neighborhood relationship.
Definition 2.
(Fuzzy neighborhood relationship) Given a decision system D S = U , C D , for an arbitrary sample x U , the fuzzy neighborhood subset of x is defined as:
[ x ] δ = y U | r x , y δ ,
where δ is fuzzy neighborhood radius. The range of δ is 0 δ 1 . If x and y are continuous data, we have r x , y = 1 1 n i = 1 n x i y i 2 . While the two elements x and y are discrete data, we have
r x , y = 1 ,   if   x i = y i 0 ,   if   x i y i ,
Thus, the fuzzy neighborhood subset is also called equivalence class. On this basis, the fuzzy conditional probability of x could be described as:
P ˜ X | [ x ] δ = r x , y | y X [ x ] δ r x , z | z [ x ] δ ,
where X is the subset of samples with the same label d k . Under the assumption of the fuzzy neighborhood subset [ x ] δ X , we have 0 < P ˜ X | [ x ] δ 1 , while P ˜ X | [ x ] δ = 1 if and only if [ x ] δ X . { } represents the sum of all elements in its set.

3.2. Determination of the Two Threshold Parameters

Considering the disadvantage of the SPDTRS model, a new loss function matrix is proposed, which is under fuzzy neighborhood relationship by a concise loss function relationship to avoid introducing the truncation functions. The novel SPDTRS model avoids the discussion of multiple situations and reduces the computational complexity of the SPDTRS model.
In the new loss function matrix, the data-driven loss functions under fuzzy neighborhood relationship is shown in Table 2. Besides, P ˜ X | [ x ] δ is the fuzzy neighborhood conditional probability, which can be calculated by Equation (8). The compensation coefficient is ζ with 0 ζ < 1 . S ˜ X | [ x ] δ and S ˜ c X | [ x ] δ are the significance coefficients, which can be described as follows:
S ˜ X | [ x ] δ = P ˜ X | [ y ] δ | y X [ x ] δ P ˜ X | [ z ] δ | z X ,
S ˜ c X | [ x ] δ = P ˜ X c | [ y ] δ | y X c [ x ] δ P ˜ X | [ z ] δ | z X .
where { } represents the sum of all elements in its set.
Under the assumption of the equivalence class [ x ] δ X , the relationships S ˜ X | [ x ] δ > 0 and S ˜ c X | [ x ] δ 0 hold, and S ˜ c X | [ x ] δ = 0 if and only if [ x ] δ X .
Subsequently, we can conclude the pair of threshold parameters according to the fuzzy neighborhood data-driven loss function matrix, which can be represented as follows. In addition, we rewrite S ˜ X | [ x ] δ = S , S ˜ c X | [ x ] δ = S c , P ˜ X | [ x ] δ = P for the sake of convenience.
α f n = λ P N λ B N λ P N λ B N + λ B P λ P P = S C S C 1 P ζ S C S C 1 P ζ + S P ζ 0 = S C 1 ζ + P ζ S C 1 ζ + P ζ + S P ζ ,
β f n = λ B N λ N N λ B N λ N N + λ N P λ B P = S C 1 P ζ 0 S C 1 P ζ 0 + S S P ζ = S C 1 P ζ S C 1 P ζ + S 1 P ζ .
Subsequently, we can set up three-way decision rules as follows:
Rule (P): Decide x P O S X while P ˜ X | [ x ] δ α f n ;
Rule (B): Decide   x B N D X while β f n < P ˜ X | [ x ] δ < α f n ;
Rule (N): Decide x N E G X while P ˜ X | [ x ] δ β f n .
According to the decision rules, the following roots about α f n and β f n can be described in the following two cases.
Case 1:
0 < ζ < 1
From the rule (P), we can obtain
P 2 ζ S C S C S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C ,
P 2 ζ S C S C + S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C .
From the rule (N), we can obtain
P 2 ζ S C + S S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C ,
P 2 ζ S C + S + S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C .
From these results, we can only accept two roots because of the relationship 0 < P 1 . Thus, we can rewrite these two roots:
α 1 f n = 2 ζ S C S C + S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C ,
β 1 f n = 2 ζ S C + S S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C .
Case 2:
ζ = 0
The values of these two parameters are as follows:
α 2 f n = 1 ,
β 2 f n = 0 .
The model arising from this case corresponds to the Pawlak model. Both boundary loss functions are equal to 0. Therefore the model clearly exhibits a two-way decision-making characteristic. Thus, the Pawlak model is one of the specific examples of our model, and ζ = 0 is a necessary non-sufficient condition for it.
Theorem 1.
For the fuzzy neighborhood data-driven loss function matrix, assuming the equivalence class [ x ] δ X , when [ x ] δ X holds, namely, the concerned equivalence class is not a consistent class, then we have:
(a1)  λ P P λ B P < λ N P ,
(a2)  λ N N λ B N < λ P N .
When [ x ] δ X holds, i.e., the concerned equivalence class is a consistent class, then we have:
(b1)  λ P P λ B P < λ N P
(b2)  λ N N = λ B N = λ P N = 0
Proof. 
 
(a1) If [ x ] δ X , then S ˜ X | [ x ] δ > 0 , 0 < P ˜ X | [ x ] δ < 1 . Since 0 ζ < 1 , λ B P = S ˜ X | [ x ] δ P ˜ X | [ x ] δ ζ , then 0 λ B P < S ˜ X | [ x ] δ . Due to λ P P = 0 and λ N P = S ˜ X | [ x ] δ , hence λ P P λ B P < λ N P .
(a2) If [ x ] δ X , then S ˜ c X | [ x ] δ > 0 , 0 < P ˜ X | [ x ] δ < 1 . Since 0 ζ < 1 , λ B N = S ˜ c X | [ x ] δ 1 P ˜ X | [ x ] δ ζ , then 0 λ B N < S ˜ c X | [ x ] δ . Due to λ N N = 0 and λ P N = S ˜ c X | [ x ] δ , hence λ N N λ B N < λ P N .
(b1) If [ x ] δ X , then S ˜ X | [ x ] δ > 0 , P ˜ X | [ x ] δ = 1 . Since 0 ζ < 1 , λ B P = S ˜ X | [ x ] δ P ˜ X | [ x ] δ ζ , then 0 λ B P < S ˜ X | [ x ] δ . Due to λ P P = 0 and λ N P = S ˜ X | [ x ] δ , hence λ P P λ B P < λ N P .
(b2) If [ x ] δ X , then S ˜ c X | [ x ] δ = 0 , and P ˜ X | [ x ] δ = 1 . Since 0 ζ < 1 , λ B N = S ˜ c X | [ x ] δ 1 P ˜ X | [ x ] δ ζ , then λ B N = 0 . Due to λ N N = 0 and λ P N = S ˜ c X | [ x ] δ , hence λ N N = λ B N = λ P N = 0 . QED. □

3.3. Establishment of FNζDTRS

Reasoning by Section 3.2, the expressions of these two threshold functions lead to the following results: α f n = f S , S c , ζ and β f n = f S , S c , ζ under the above two conditions, where α f n = α 1 f n , α 2 f n and β f n = β 1 f n , β 2 f n . Based on Equations (9) and (10), it is easily to obtain S and S c by confirming the parameter δ and analyzing the distribution information of the original data. In summary, we can change these two threshold parameters to α f n = f δ , ζ , β f n = f δ , ζ . The rough set model below, which is including parameters δ and ζ is defined as fuzzy neighborhood ζ-decision-theoretic rough set (FNζDTRS):
P O S ˜ = x U | P ˜ X | [ x ] δ α f n ,
B N D ˜ = x U | β f n < P ˜ X | [ x ] δ < α f n ,
N E G ˜ = x U | P ˜ X | [ x ] δ β f n ,
where both threshold parameters have different descriptions under the following two conditions:
Case 1:
0 < ζ < 1 , then
α 1 f n = 2 ζ S C S C + S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C ,
β 1 f n = 2 ζ S C + S S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C .
Case 2:
ζ = 0 , then
α 2 f n = 1 ,
β 2 f n = 0 .
In the above FNζDTRS model, there are only two cases to discuss, in contrast to the SPDTRS model that requires four cases to discuss, which greatly reduces the computational complexity of the SPDTRS model due to the concise setting of loss functions.
Theorem 2.
In the FNζDTRS model, given two compensation coefficients ζ 1 and ζ 2 with 0 ζ 1 < 1 and 0 ζ 2 < 1 , and the parameter δ is fixed, if there exists ζ 1 ζ 2 , then α f n ζ 1 α f n ζ 2 and β f n ζ 1 β f n ζ 2 hold.
Proof. 
If the equivalence class [ x ] δ X , then S C = 0 , according to Equations (22)–(25), when 0 < ζ < 1 , α f n ζ 1 = α f n ζ 2 = 0 , β f n ζ 1 = β f n ζ 2 = 0 , when ζ = 0 , α f n ζ 1 = α f n ζ 2 = 1 , β f n ζ 1 = β f n ζ 2 = 0 . If the equivalence class [ x ] δ X , then its monotonicity relations will be proved by the following derivations.
Part I: For ζ 1 ζ 2 α f n ζ 1 α f n ζ 2 , two cases need to be considered.
Case 1:
0 < ζ < 1
Since α 1 f n = 2 ζ S C S C + S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C , we set η = 1 / ζ , due to 0 < ζ < 1 , then η > 1 and α 1 f n = 2 S C η S C + η 2 ( S C ) 2 + 4 η S S C 4 S S C 2 S + S C . Due to S > 0 , S C > 0 , we can set α ¯ 1 f n = 2 S C η S C + η 2 ( S C ) 2 + 4 η S S C 4 S S C for simplicity of derivation. We can find the partial derivative of it, denoted as f α 1 , which is
f α 1 = α ¯ 1 f n η = S C + η S C 2 + 2 S S C η 2 S C 2 + 4 η S S C 4 S S C .
Its second-order partial derivative is denoted as f α 2 :
f α 2 = 2 α ¯ 1 f n 2 η = 4 S S C 2 S + S C η 2 S C 2 + 4 η S S C 4 S S C 3 / 2 < 0 .
Because f α 2 is less than 0, f α 1 is monotonically decreasing. Since η > 1 , f α 1 2 S > 0 1 η , f α 1 0 η + , then f α 1 > 0 . Therefore, α 1 f n grows monotonically with respect to η . Hence, α 1 f n decreases monotonically with respect to ζ , that is ζ 1 ζ 2 α 1 f n ζ 1 α 1 f n ζ 2 .
Case 2: ζ = 0
In this case, we have α 2 f n = 1 . Thus, for ζ 1 ζ 2 , we have α 2 f n ζ 1 = α 2 f n ζ 2 .
Part II: For ζ 1 ζ 2 β f n ζ 1 β f n ζ 2 , two cases need to be considered as well.
Case 1:
0 < ζ < 1
Since β 1 f n = 2 ζ S C + S S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C , we also set η = 1 / ζ , due to 0 < ζ < 1 , then η > 1 and β 1 f n = 2 S C + η S η 2 S 2 + 4 η S S C 4 S S C 2 S + S C . Due to S > 0 , S C > 0 , we have the simple form β 1 f n = 2 S C + η S η 2 S 2 + 4 η S S C 4 S S C . We describe partial derivatives of β 1 f n like:
f β 1 = β ¯ 1 f n η = S η S 2 + 2 S S C η 2 S 2 + 4 η S S C 4 S S C , f β 2 = 2 β ¯ 1 f n 2 η = 4 S 2 S C S + S C η 2 S 2 + 4 η S S C 4 S S C 3 / 2 > 0 .
From f β 2 > 0 , we know that f β 1 increases monotonously along with η . Since η > 1 , f β 1 2 S C < 0 1 η , f β 1 0 η + , then f α 1 < 0 , and β 1 f n is monotonously decreasing with regard to η . Therefore, β 1 f n is monotonously increasing with ζ , that is ζ 1 ζ 2 β 1 f n ζ 1 β 1 f n ζ 2 . QED. □
Theorem 3.
In the FNζDTRS model, the relationship 0 β f n α f n 1 holds.
Proof. 
If the equivalence class [ x ] δ X , then S C = 0 , according to Equations (22)–(25), when 0 < ζ < 1 , α f n = β f n = 0 , when ζ = 0 , α f n = 0 , β f n = 1 , satisfying the inequality 0 β f n α f n 1 . If the equivalence class [ x ] δ X , then in following three parts we intend to prove the inequality.
Part I: Proof of 0 β f n under two cases.
Case 1:
0 < ζ < 1
According to Equation (23), we could get β 1 f n = 2 ζ S C + S ψ 2 ζ S + S C ,   ψ = ( S + 2 ζ S C ) 2 4 ζ 2 S C S + S C , and then we only need to prove 4 ζ 2 S C S + S C > 0 . Since S > 0 , S C > 0 , then 4 ζ 2 S C S + S C > 0 , so β 1 f n > 0 .
Case 2:
ζ = 0
In this case, we have β 2 f n = 0 .
Part II: Proof of the inequality β f n α f n in two cases.
Case 1:
0 < ζ < 1
According to Equations (22) and (23), we could get
α 1 f n β 1 f n = S C S + S C 4 ζ 2 S + 4 ζ S + S C + S 4 ζ 2 S C + 4 ζ S C + S 2 ζ S + S C .
Since S > 0 , S C 0 , then
S C 4 ζ 2 S + 4 ζ S + S C = 4 ζ S S C 1 ζ + S C 2 > S C 2 = S C , S 4 ζ 2 S C + 4 ζ S C + S = 4 ζ S S C 1 ζ + S 2 > S 2 = S ,
so α 1 f n β 1 f n > 0 , that is β 1 f n < α 1 f n .
Case 2:
ζ = 0
In this case, we have α 2 f n = 1 , β 2 f n = 0 , so β 2 f n < α 2 f n .
Part III: Proof of the inequality α f n 1 in two cases.
Case 1:
0 < ζ < 1
According to Equation (22), if we want to prove α 1 f n 1 , then we need to prove S C 4 ζ 2 S + 4 ζ S + S C 2 ζ S + S C < 0 , which means we need to prove S C 4 ζ 2 S + 4 ζ S + S C ( 2 ζ S + S C ) 2 < 0 . Since S > 0 , S C 0 , then
S C 4 ζ 2 S + 4 ζ S + S C ( 2 ζ S + S C ) 2 = 4 ζ 2 S S C 4 ζ 2 S 2 = 4 ζ 2 S S C + S < 0 .
Therefore, we have α 1 f n < 1 .
Case 2:
ζ = 0
In this case, we have α 2 f n = 1 . QED. □
Theorem 4.
For a decision system, which is described as D S = U , C D with a fixed parameter δ , and two parameters ζ 1 and ζ 2 with 0 ζ 1 ζ 2 < 1 , we have P O S ˜ 1 P O S ˜ 2 , N E G ˜ 1 N E G ˜ 2 , B N D ˜ 1 B N D ˜ 2 .
Proof. 
At the very beginning of the proof, we assume an arbitrary sample y to facilitate the proof of the theorem. While y P O S ˜ 1 , we have P ( X | y ) α f n ζ 1 . Since 0 ζ 1 ζ 2 < 1 , the relation α f n ζ 1 α f n ζ 2 holds according to Theorem 2. Thus, P ( X | y ) α f n ζ 2 and y P O S ˜ 2 hold. Hence, we conclude that P O S ˜ 1 P O S ˜ 2 .
Likewise, to conclude that N E G ˜ 1 N E G ˜ 2 and B N D ˜ 1 B N D ˜ 2 is easy via Theorem 2.
From the above, we can find that ζ is inversely correlated with the range of neutrality and positively correlated with the uncertainty of the decision. QED. □

3.4. FNζDTRS-Based Attribute Reduction Algorithm

Jia et al. [35] presented a reduction principle in response to the problem of attribute reduction by using the DTRS model. The core idea of it is minimizing the risk of the reduction subset. On this principle, our previous work has designed an attribute reduction algorithm based on the SPDTRS model [32,34]. It is built on minimizing the risk of overall decisions, which can be utilized for the attribute reduction in our proposed FNζDTRS model. Therefore, the detailed attribute reduction algorithm is not repeated in this paper. For details, the readers could refer to our previous work [32,34].

4. Strategy of Simultaneous-Fault Diagnosis

Under the assumption that the data in all single-fault modes are fully available, when an unknown fault occurs, we use the FNζDTRS model to mine the fault attributes of the unknown fault. If the fault attributes of the unknown fault are different from the fault attributes of the existing single-fault, then the unknown fault can be considered as a simultaneous-fault. Furthermore, we can use the attribute reduction results obtained from the FNζDTRS model to analyze and identify the corresponding fault modes of the simultaneous-fault. Finally, a strategy of simultaneous-fault diagnosis called fault matching strategy is formed, as shown in Figure 2.
The proposed fault matching strategy consists of two main parts, prior knowledge acquisition and rule matching. In the first part, the single-fault data with abnormal labels and normal data with normal labels are sent into the FNζDTRS model as the training data set, and the optimal fault attribute subsets of each single-fault are obtained by attribute reduction, used as the prior knowledge for subsequent diagnosis. In the second part, the simultaneous-fault data with abnormal labels and normal data with normal labels are combined to form the data to be diagnosed, and then the data are fed into the FNζDTRS model to obtain the optimal fault attribute subset. The optimal fault attribute subset is then obtained here and the optimal fault attribute subsets obtained in the first part are measured by using the Jaccard similarity coeffective. It is worth noting that there may be some single faults with the same fault attributes, therefore we set some rules based on the differences between attribute data to subdivide the faults and complete the fault matching, which can be seen in the subsequent experiments based on the satellite power system.
Based on the above description, we can write the core pseudo code in the above process, as shown in Table 3.

5. Numerical Experiment of Attribute Reduction

The effectiveness and advantage of the proposed FNζDTRS model is verified on several hybrid decision systems from the UCI (http://archive.ics.uci.edu/ml/index.php, accessed on 21 July 2021) and KEEL (https://sci2s.ugr.es/keel/datasets.php, accessed on 21 July 2021) datasets. As shown in Table 4, the test datasets include both discrete and continuous data. Specific comparative experiments regarding parameters test and attribute reduction are conducted on the same hard and soft platforms. Ten baseline classifiers, including NaiveBayes, REPTree, LogitBoost, SMO, Filtered, Bagging, PART, IBk, J48 and JRip, are employed with a 10-fold cross-validation in Weka (https://waikato.github.io/weka-wiki/downloading_weka/, accessed on 21 July 2021) software to demonstrate the accuracy of attribute selection. The input data are normalized into the range of [0, 1] during preprocessing.

5.1. Parameters Test for FNζDTRS

For the FNζDTRS model, two parameters ζ and δ need to be set in advance. As described in Section 2, the theoretic value field of ζ is 0 , 1 , and δ is 0 , 1 . Therefore, ζ is sampled with an interval of 0.05, and end at 0.99. δ is also sampled with an interval of 0.05, but end at 1. Figure 3 shows the experimental results.
The results show that the appropriate settings of ζ and δ range in [0.3, 0.99] and [0.85, 0.95], respectively. It could be explained by the fact that a smaller ζ will result in a larger boundary region, marking more samples as uncertain state. It means setting a smaller ζ for the FNζDTRS decision system will lead to greater uncertainty. When applied in real world, ζ should be adjusted appropriately according to the risk of wrong decision. When the danger of making a bad decision is low, ζ can be set larger, and vice versa. On the other hand, with the fuzzy neighborhood threshold δ closer to 1, the fuzzy neighborhood granules will be finer, allowing for the samples to be classified accurately into the appropriate regions. The above two parameters are the core parameters of the model proposed in this paper, and their setting values directly affect the test results. Therefore, when setting the above parameters, the values of the two parameters need to be adjusted according to the actual needs with the above principles.

5.2. Comparison Experiments on Attribute Reduction

In this part, seven related models, DTRS-EF [36], DTRS-SMDNS [37], SPDTRS-EF [32], SPDTRS-SMDNS [32], NDTRS [38], FDTRS [39] and FN3WD [34], are introduced into a contrastive analysis on the attribute reduction to demonstrate the superiority of the proposed FNζDTRS. The settings of the relevant parameters in these comparison models are the same as those of the corresponding models. The number of reduction attributes and the classification accuracy are the common evaluation indicators of the comparison experiment of attribute reduction [40,41]. The standard deviations of ten trials are also calculated, and the results are shown in Table 5 and Table 6.
According to the results in Table 4 and Table 5, the following analysis can be obtained:
(a)
The analysis based on the classification accuracy indicates that the FNζDTRS model is superior to other rough set models. The main reason may lie in the different methods to describe spatial granules. Discretization methods such as EF and SMDNS are commonly introduced to process continuous data in the traditional DTRS models, which results in the destruction of the spatial structure of granules. Using special measures (such as fuzzy relationship, neighborhood relationship, etc.) can avoid the distortion of the discretization method, but it also has some disadvantages, such as simple measurement, insufficient description ability, etc. The proposed FNζDTRS model utilizes fuzzy neighborhood relationships to overcome the above shortcomings. Compared with other DTRS models, the description of spatial granules is more precise in our model and results in the higher classification accuracy.
(b)
With respect to the number of reduction attributes, the FDTRS model has the least number of reduction attributes, but it fails to achieve a desired classification accuracy, whereas the FNζDTRS model can maintain high classification accuracy while keeping the number of reduction attributes small. The results show that the classification ability can be maintained or improved only when the reduction attributes are accurately selected. The above conclusion also conforms to the basic principle of attribute reduction, that is, in the operation of reduction, we want to get a relatively concise set, which can ensure that the original classification accuracy is not reduced, and the purpose is to improve the operation efficiency.
(c)
The standard deviation is used to measure the robustness of models. It is obvious that the standard deviation of the FNζDTRS model is the smallest regardless of the classification accuracy or the number of reduced attributes, which directly proves that the robustness of the FNζDTRS model is the highest compared to other models. The above robustness characteristics also show that we have a large selection range when setting our two parameters, which is conducive to the wide application of the model in practical projects.

6. Simultaneous-Fault Diagnosis of Satellite Power System

In-orbit faults of the power system should be avoided to the maximum extent for satellites. Therefore, simulation is the best platform to mine fault diagnosis knowledge. In this section, the effectiveness of the proposed simultaneous-fault diagnosis scheme is verified with the simulation model of a geosynchronous (GEO) satellite power system [3]. As shown in Figure 4, the power system works in a direct energy transfer mode during the simulation, and ten telemetry parameters can be measured in the marked position. The information of the telemetry parameters is shown in Table 7.
The raw data used for simultaneous-fault diagnosis is composed of the above-mentioned ten kinds of attributes, and all the data are selected in the stationary period. The dataset is stored in a time-series format, with each subset representing one of the scenarios shown in Table 8. There are a total of 12 scenarios, where scenario 0 represents the system without any fault. F1 represents open-circuit failure in solar array. F2 represents the short-circuit failure in the battery. F3 represents shunt regulator failure without shunt. F4 represents shunt regulator failure with constant shunt. The remaining 7 scenarios are concurrent failures composed of the above 4 single failures occurring at the same time, where F3 and F4 cannot occur simultaneously.
It can be seen that this means can effectively solve the problem of insufficient data in the simultaneous-fault diagnosis, which also responds to one the difficulties in the simultaneous-fault diagnosis introduced at the beginning of this paper.

6.1. Simultaneous-Fault Diagnosis Based on the Fault Matching Strategy

The two main parts of the simultaneous-fault diagnosis strategy, namely prior knowledge acquisition and rule matching, are equivalent to the training and testing process. We choose 4 kinds of single-fault data as the training set and the remaining 7 kinds of simultaneous-fault data as the testing set. The results obtained through the first step of prior knowledge acquisition are shown in Table 9. It can be found that the results of the output attribute subset of F3 and F4 are the same. Therefore, in order to distinguish F3 and F4, further information needs to be excavated. For attribute a3, its corresponding shunt current data can directly distinguish F3 from F4. The shunt current value of F3 fluctuates between 6.42–8.67 and that of F4 is between 12.45–14.77. Therefore, F3 and F4 can be distinguished by setting the threshold value of the shunt current average.
The results obtained through the second step of rule matching are shown in Table 10. For a simultaneous-fault, the Jaccard similarity coefficient between its output attribute subset and the output attribute subset of each single-fault obtained in the training process can be calculated in turn. If the Jaccard similarity coefficient is 0, the corresponding single fault can be eliminated preliminarily. Since F3 and F4 cannot be distinguished by the Jaccard similarity coefficient, it is necessary to further distinguish F3 and F4 through the set shunt current threshold and to finally obtain the matching result. Through the final matching result, it can be found that the diagnostic accuracy of the fault matching strategy is 100%. The above fault matching process comprehensively utilizes the similarity of attributes and expert knowledge, which can ensure that the obtained diagnosis results are more accurate.

6.2. Comparison Experiment on Simultaneous-Fault Diagnosis

6.2.1. Experimental Setup

The superiority of the proposed fault matching strategy is demonstrated through comparison experiments of simultaneous-fault diagnosis with several multi-label classification algorithms (Binary Relevance, Classifier Chain, Calibrated Label Ranking, ML-KNN, and ML-DT). In these comparison algorithms, the classifiers of the first three algorithms are all set to Random Forest, which is the best classifier after the pretest, and the value of k for ML-KNN is 12. Subset accuracy, hamming loss, precision, recall, and F1 are introduced as the metrics of the multi-label classification performance [11].
Subset   Accuracy = 1 n i = 1 n I y i = y ^ i . ,
Hamming   Loss = 1 n L i = 1 n j = 1 L I y i j y ^ i j ,
Precision = 1 n i = 1 n y i j = 1 y ^ i j = 1 y i j = 1 ,
Recall = 1 n i = 1 n y i j = 1 y ^ i j = 1 y ^ i j = 1 ,
F 1 = 1 n n i = 1 2 y i j = 1 y ^ i j = 1 y i j = 1 + y ^ i j = 1 ,
where y i represents the ground-truth label vector of the i sample, y ^ i represents the predicted label vector, y i j represents the ground-truth label of j position in i sample, y ^ i j represents the predicted label of j position in i sample. Intuitively, subset accuracy, precision, recall, F1 perform as the multi-label counterparts of traditional metrics. Hamming lose performs as a special metric of multi-label.
The training set of all algorithms uses single-fault data, while the test set uses simultaneous-fault data. Ten-fold cross validation is also employed in this part.

6.2.2. Experimental Results and Analysis

The results are shown in Table 11. It is obvious that the proposed method outperforms other methods in the simultaneous-faults diagnosis of a satellite power system, which benefits from both the FNζDTRS model and the fault matching strategy (FNζDTRS-FMS). Fundamentally, on the one hand, the attribute reduction results obtained from the FNζDTRS model can accurately represent the fault information, and on the other hand, the proposed matching rules are reasonable and reliable.
In addition to the results of our method, as for multi-label classification methods, the false alarm rate performs satisfactorily, while the missed diagnosis rate performs poorly. It can be seen that the recall rates have been lower than 50%, which is unacceptable in engineering applications. At the same time, the results of the F1 index generated by the multi-label classification methods are relatively low, only slightly higher than 60%. Similarly, the results obtained by the Calibrated Label Ranking method are similar to those of ML-KNN and ML-DT. The results obtained by the Classifier Chain method are in the middle level. In addition, Binary Relevance has the best performance among the multi-label classification methods, which means that there are no relevant dependencies between the single faults.

7. Conclusions

In this work, a novel DTRS model called FNζDTRS is proposed and the fault match strategy (FMS) based on the FNζDTRS model is designed to overcome three fundamental hurdles faced by simultaneous-fault diagnosis. The effectiveness and superiority of our methodology is demonstrated by both numerical experiments conducted on several standard datasets and comparison analysis of simultaneous-fault diagnosis performed on a simulation model of a satellite power system. Consequently, two main conclusions can be drawn, as follows.
(1)
The proposed FNζDTRS model performs attribute reduction more effectively compared with other models, and it has strong generalization ability. This benefits from the concise loss functions and the introduction of the fuzzy neighborhood relationships. The advantages of our model can greatly promote the smooth implementation of the model in simultaneous-fault diagnosis, which reflects the effectiveness and superiority of our selection of this model.
(2)
The proposed FNζDTRS–FMS does not require simultaneous-fault samples to accomplish training and performs excellently in simultaneous-fault diagnosis compared to classic multi-label classification algorithms. This is completely consistent with the real situation, that is, the existing data cannot completely cover all the imagined failure modes. Therefore, the diagnostic strategy proposed in this paper has stronger application value.
Although the model we proposed has the above advantages, it still has the problem of low computational efficiency compared with the classical rough set because of the use of fuzzy neighborhood computing. This is one of our future research directions. Furthermore, our future work will also focus on fusing rough set models and multi-label learning algorithms to make the simultaneous-fault diagnosis framework more general.

Author Contributions

Data curation, R.Z.; Formal analysis, Y.J. and T.Z.; Funding acquisition, C.L. and M.S.; Methodology, L.T. and C.W.; Writing—original draft, Y.C., Y.J. and M.S.; Validation, Y.J.; Investigation, T.Z.; Supervision, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (Grant Nos. 61903015 and 61973011), the Fundamental Research Funds for the Central Universities (Grant No. KG21003001), National key Laboratory of Science and Technology on Reliability and Environmental Engineering (Grant No. WDZC2019601A304), as well as the Capital Science & Technology Leading Talent Program (Grant No. Z191100006119029).

Data Availability Statement

The public data can be accessed through the following two websits: the UCI (http://archive.ics.uci.edu/ml/index.php, accessed on 21 July 2021) and KEEL (https://sci2s.ugr.es/keel/datasets.php, accessed on 21 July 2021). However, the data set of the satellite power system is not convenient for public display. If necessary, researchers can make further consultation via email.

Acknowledgments

Thank the reviewers for their comments on improving the quality of our papers, and thank the journal editors for their enthusiastic work.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

U = x 1 , x 2 , , x m the universe, which is a finite and nonempty set.
D the set of decision attributes that is a nonempty set.
C the collection of conditional attributes.
X the subset of samples with the same label d k .
a P , a B , and a N the classification of x into three regions, which are x P O S X , x B N D X , x N E G X .
P O S X the acceptance of the event x X .
B N D X the non-commitment of the event x X , denotes the deferment of the event x X .
N E G X the rejection of x X .
λ P the loss caused by taking actions ( a P , a B , a N ) while x X .
λ N the loss caused by taking actions ( a P , a B , a N ) while x X .
α , β the threshold parameters of the DTRS model.
δ fuzzy neighborhood radius.
ζ the compensation coefficient.

References

  1. Li, S.; Cao, H.; Yang, Y. Data-driven simultaneous fault diagnosis for solid oxide fuel cell system using multi-label pattern identification. J. Power Sources 2018, 378, 646–659. [Google Scholar] [CrossRef]
  2. Suo, M.; Tao, L.; Zhu, B.; Chen, Y.; Lu, C.; Ding, Y. Soft decision-making based on decision-theoretic rough set and Takagi-Sugeno fuzzy model with application to the autonomous fault diagnosis of satellite power system. Aerosp Sci. Technol. 2020, 106, 106108. [Google Scholar] [CrossRef]
  3. Suo, M.; Zhu, B.; An, R.; Sun, H.; Xu, S.; Yu, Z. Data-driven fault diagnosis of satellite power system using fuzzy Bayes risk and SVM. Aerosp Sci. Technol. 2019, 84, 1092–1105. [Google Scholar] [CrossRef]
  4. Asgari, S.; Gupta, R.; Puri, I.K.; Zheng, R. A data-driven approach to simultaneous fault detection and diagnosis in data centers. Appl. Soft Comput. 2021, 110, 107638. [Google Scholar] [CrossRef]
  5. Liang, P.; Deng, C.; Wu, J.; Yang, Z.; Zhu, J.; Zhang, Z. Single and simultaneous fault diagnosis of gearbox via a semi-supervised and high-accuracy adversarial learning framework. Knowl.-Based Syst. 2020, 198, 105895. [Google Scholar] [CrossRef]
  6. Zhang, Z.; Li, S.; Xiao, Y.; Yang, Y. Intelligent simultaneous fault diagnosis for solid oxide fuel cell system based on deep learning. Appl. Energ. 2019, 233, 930–942. [Google Scholar] [CrossRef]
  7. Pooyan, N.; Shahbazian, M.; Salahshoor, K.; Hadian, M. Simultaneous Fault Diagnosis using multi class support vector machine in a Dew Point process. J. Nat. Gas. Sci. Eng. 2015, 23, 373–379. [Google Scholar] [CrossRef]
  8. Vong, C.; Wong, P.; Ip, W. A New Framework of Simultaneous-Fault Diagnosis Using Pairwise Probabilistic Multi-Label Classification for Time-Dependent Patterns. Ieee T Ind. Electron. 2013, 60, 3372–3385. [Google Scholar] [CrossRef]
  9. Wong, P.K.; Zhong, J.; Yang, Z.; Vong, C.M. Sparse Bayesian extreme learning committee machine for engine simultaneous fault diagnosis. Neurocomputing 2016, 174, 331–343. [Google Scholar] [CrossRef]
  10. Wu, B.; Cai, W.; Chen, H.; Zhang, X. A hybrid data-driven simultaneous fault diagnosis model for air handling units. Energy Build. 2021, 245, 111069. [Google Scholar] [CrossRef]
  11. Zhang, M.; Zhou, Z. A Review on Multi-Label Learning Algorithms. IEEE Trans. Knowl. Data Eng. 2014, 26, 1819–1837. [Google Scholar] [CrossRef]
  12. Boutell, M.R.; Luo, J.B.; Shen, X.P.; Brown, C.M. Learning multi-label scene classification. Pattern Recognit. 2004, 37, 1757–1771. [Google Scholar] [CrossRef]
  13. Read, J.; Pfahringer, B.; Holmes, G.; Frank, E. Classifier chains for multi-label classification. Mach. Learn. 2011, 85, 333–359. [Google Scholar] [CrossRef]
  14. Fuernkranz, J.; Huellermeier, E.; Mencia, E.L.; Brinker, K. Multilabel classification via calibrated label ranking. Mach. Learn. 2008, 73, 133–153. [Google Scholar] [CrossRef]
  15. Zhang, M.; Zhou, Z. ML-KNN: A lazy learning approach to multi-label leaming. Pattern Recognit. 2007, 40, 2038–2048. [Google Scholar] [CrossRef]
  16. Clare, A.; King, R.D. Knowledge Discovery in Multi-label Phenotype Data. In European Conference on Principles of Data Mining and Knowledge Discovery; Springer: Berlin/Heidelberg, Germany, 2001; pp. 42–53. [Google Scholar]
  17. Pawlak, Z. Rough sets. Int. J. Comput. Inf. Sci. 1982, 11, 341–356. [Google Scholar] [CrossRef]
  18. Dong, L.; Chen, D.; Wang, N.; Lu, Z. Key energy-consumption feature selection of thermal power systems based on robust attribute reduction with rough sets. Inf. Sci. 2020, 532, 61–71. [Google Scholar] [CrossRef]
  19. Su, L.; Yu, F. Matrix approach to spanning matroids of rough sets and its application to attribute reduction. Theor. Comput. Sci. 2021, 893, 105–116. [Google Scholar] [CrossRef]
  20. Sahu, R.; Dash, S.R.; Das, S. Career selection of students using hybridized distance measure based on picture fuzzy set and rough set theory. Decis. Mak. Appl. Manag. Eng. 2021, 4, 104–126. [Google Scholar] [CrossRef]
  21. Dash, S.R.; Dehuri, S.; Sahoo, U.K. Interactions and Applications of Fuzzy, Rough, and Soft Set in Data Mining. Int. J. Fuzzy Syst. Appl. 2015, 3, 37–50. [Google Scholar] [CrossRef]
  22. Zhang, P.; Li, T.; Wang, G.; Luo, C.; Chen, H.; Zhang, J.; Wang, D.; Yu, Z. Multi-source information fusion based on rough set theory: A review. Inf. Fusion 2021, 68, 85–117. [Google Scholar] [CrossRef]
  23. Bai, H.; Ge, Y.; Wang, J.; Li, D.; Liao, Y.; Zheng, X. A method for extracting rules from spatial data based on rough fuzzy sets. Knowl.-Based Syst. 2014, 57, 28–40. [Google Scholar] [CrossRef]
  24. Landowski, M.; Landowska, A. Usage of the rough set theory for generating decision rules of number of traffic vehicles. Transp. Res. Procedia 2019, 39, 260–269. [Google Scholar] [CrossRef]
  25. Sharma, H.K.; Kumari, K.; Kar, S. A rough set theory application in forecasting models. Decis. Mak. Appl. Manag. Eng. 2020, 3, 1–21. [Google Scholar] [CrossRef]
  26. Guo, Y.; Tsang, E.C.C.; Xu, W.; Chen, D. Adaptive weighted generalized multi-granulation interval-valued decision-theoretic rough sets. Knowl.-Based Syst. 2020, 187, 104804. [Google Scholar] [CrossRef]
  27. Wang, T.; Liu, W.; Zhao, J.; Guo, X.; Terzija, V. A rough set-based bio-inspired fault diagnosis method for electrical substations. Int. J. Electr. Power Energy Syst. 2020, 119, 105961. [Google Scholar] [CrossRef]
  28. Sang, B.; Yang, L.; Chen, H.; Xu, W.; Guo, Y.; Yuan, Z. Generalized multi-granulation double-quantitative decision-theoretic rough set of multi-source information system. Int. J. Approx. Reason. 2019, 115, 157–179. [Google Scholar] [CrossRef]
  29. Zhang, P.F.; Li, T.R.; Yuan, Z.; Luo, C.; Liu, K.Y.; Yang, X.L. Heterogeneous Feature Selection Based on Neighborhood Combination Entropy. IEEE Trans. Neural Netw. Learn. Syst. 2022, 1–14. [Google Scholar] [CrossRef]
  30. Wang, L.; Shen, J.; Mei, X. Cost Sensitive Multi-Class Fuzzy Decision-theoretic Rough Set Based Fault Diagnosis. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017; pp. 6957–6961. [Google Scholar]
  31. Yu, J.; Ding, B.; He, Y. Rolling bearing fault diagnosis based on mean multigranulation decision-theoretic rough set and non-naive Bayesian classifier. J. Mech. Sci. Technol. 2018, 32, 5201–5211. [Google Scholar] [CrossRef]
  32. Suo, M.; Tao, L.; Zhu, B.; Miao, X.; Liang, Z.; Ding, Y.; Zhang, X.; Zhang, T. Single-parameter decision-theoretic rough set. Inf. Sci. 2020, 539, 49–80. [Google Scholar] [CrossRef]
  33. Yao, Y.Y.; Wong, S.K.M. A decision theoretic framework for approximating concepts. Int. J. Man Mach. Stud. 1992, 37, 793–809. [Google Scholar] [CrossRef]
  34. Suo, M.; Cheng, Y.; Zhuang, C.; Ding, Y.; Lu, C.; Tao, L. Extension of labeled multiple attribute decision making based on fuzzy neighborhood three-way decision. Neural Comput. Appl. 2020, 32, 17731–17758. [Google Scholar] [CrossRef]
  35. Jia, X.; Liao, W.; Tang, Z.; Shang, L. Minimum cost attribute reduction in decision-theoretic rough set models. Inf. Sci. 2013, 219, 151–167. [Google Scholar] [CrossRef]
  36. Yao, Y. Three-way decisions with probabilistic rough sets. Inf. Sci. 2010, 180, 341–353. [Google Scholar] [CrossRef]
  37. Jiang, F.; Sui, Y. A novel approach for discretization of continuous attributes in rough set theory. Knowl.-Based Syst. 2015, 73, 324–334. [Google Scholar] [CrossRef]
  38. Li, W.; Huang, Z.; Jia, X.; Cai, X. Neighborhood based decision-theoretic rough set models. Int. J. Approx. Reason. 2016, 69, 1–17. [Google Scholar] [CrossRef]
  39. Song, J.; Tsang, E.C.C.; Chen, D.; Yang, X. Minimal decision cost reduct in fuzzy decision-theoretic rough set model. Knowl.-Based Syst. 2017, 126, 104–112. [Google Scholar] [CrossRef]
  40. Wang, C.; Qi, Y.; Shao, M.; Hu, Q.; Chen, D.; Qian, Y.; Lin, Y. A Fitting Model for Feature Selection with Fuzzy Rough Sets. IEEE Trans. Fuzzy Syst. 2017, 25, 741–753. [Google Scholar] [CrossRef]
  41. Wang, C.; Shao, M.; He, Q.; Qian, Y.; Qi, Y. Feature subset selection based on fuzzy neighborhood rough sets. Knowl.-Based Syst. 2016, 111, 173–179. [Google Scholar] [CrossRef]
Figure 1. The illustration of a decision system.
Figure 1. The illustration of a decision system.
Mathematics 10 03414 g001
Figure 2. The procedure of fault diagnosis strategy for simultaneous-fault.
Figure 2. The procedure of fault diagnosis strategy for simultaneous-fault.
Mathematics 10 03414 g002
Figure 3. The accuracy results of ζ and δ with respect to FNζDTRS.
Figure 3. The accuracy results of ζ and δ with respect to FNζDTRS.
Mathematics 10 03414 g003aMathematics 10 03414 g003b
Figure 4. The schematic diagram of the power system.
Figure 4. The schematic diagram of the power system.
Mathematics 10 03414 g004
Table 1. The detailed information of a loss function matrix.
Table 1. The detailed information of a loss function matrix.
Q Q c
a P λ P P λ P N
a B λ B P λ B N
a N λ N P λ N N
Table 2. The fuzzy neighborhood data-driven loss function matrix.
Table 2. The fuzzy neighborhood data-driven loss function matrix.
Q Q c
a p λ P P = 0 λ P N = S ˜ c X | [ x ] δ
a B λ B P = S ˜ X | [ x ] δ P ˜ X | [ x ] δ ζ λ B N = S ˜ c X | [ x ] δ 1 P ˜ X | [ x ] δ ζ
a N λ N P = S ˜ X | [ x ] δ λ N N = 0
Table 3. The core pseudo code of the diagnosis strategy.
Table 3. The core pseudo code of the diagnosis strategy.
Input:Raw data of each single fault and normal state
Output:Fault mode
Part IPrior Knowledge Acquisition
For each DS regarding to single fault or state
   Initialized: red = ,   C l = C, Rred = H, //H is a large positive number.
   While Cl
        For c C l
             a = c   red, // a is a temporary set.
            Compute the risk generated by a .
        End For
        Find such a subset a = c   r e d with the minimum risk, i.e., R a .
        If R a < R r e d
          The subset a is the selected set
        End If
    End While
End For, return the reduction r e d set of each state.
Part IIRule Matching
Utilize the above code to obtain the reduction set r of the given fault data to be diagnosed.
For each r e d
   Compute the similarity between r e d and r .
End For
Find such a r e d with the maximum similarity, which could be considered as the similar fault mode f.
Return f
Table 4. The information of the employed datasets.
Table 4. The information of the employed datasets.
IDFull NameNameSamplesAttributeDiscreteContinuousClassSource
1Mutagenesis-AtomsAtoms161810822KEEL
2Australian Credit ApprovalAustralian69014862UCI
3Breast CancerBreast2779632UCI
4Heart Disease ClevelandCleve29613762UCI
5Statlog HeartHeart27013672UCI
6IrisIris1504043UCI
7Website PhishingPhishing1353101003UCI
8South African HearthSaheart4629182UCI
9Seismic-BumpsSeismic2584181262UCI
10Congressional Voting RecordsVote435161602UCI
Table 5. The classification accuracy of the reduction subset.
Table 5. The classification accuracy of the reduction subset.
IDNameDTRS-EFDTRS-SMDNSSPDTRS-EFSPDTRS-SMDNSNDTRSFDTRSFN3WDFNζDTRS
1Atoms69.65 ± 2.3871.08 ± 1.6670.94 ± 1.2071.27 ± 1.1470.42 ± 1.6170.87 ± 2.0771.71 ± 1.0072.08 ± 1.13
2Australian81.34 ± 5.0982.39 ± 5.6482.27 ± 2.8783.46 ± 0.8583.31 ± 2.7672.06 ± 12.3184.56 ± 0.3984.97 ± 0.42
3Breast72.31 ± 1.0672.53 ± 0.8272.67 ± 0.7072.65 ± 0.7272.72 ± 0.6770.38 ± 0.7573.21 ± 0.2974.01 ± 0.81
4Cleve79.13 ± 1.0978.34 ± 5.0278.99 ± 1.0179.37 ± 0.6177.69 ± 6.6866.64 ± 11.2980.13 ± 0.7881.34 ± 0.37
5Heart78.34 ± 3.3878.72 ± 5.2079.99 ± 1.8779.95 ± 0.8876.31 ± 7.1668.37 ± 10.4280.27 ± 2.2580.99 ± 1.23
6Iris94.85 ± 0.5094.85 ± 0.6294.93 ± 0.4794.87 ± 0.5594.95 ± 0.4162.96 ± 17.9894.82 ± 0.4095.27 ± 0.33
7Phishing84.23 ± 9.6286.05 ± 6.1187.08 ± 1.1087.08 ± 1.1086.40 ± 5.0869.93 ± 14.0687.15 ± 1.0687.14 ± 1.07
8Saheart69.23 ± 0.9969.01 ± 1.2969.49 ± 0.3869.51 ± 0.4369.35 ± 0.4067.93 ± 1.9969.63 ± 1.6170.35 ± 1.06
9Seismic92.64 ± 0.8092.33 ± 0.9092.53 ± 0.4691.94 ± 0.6991.91 ± 0.7092.01 ± 0.7692.56 ± 0.4893.21 ± 0.15
10Vote94.66 ± 0.3894.67 ± 0.3894.63 ± 0.3694.63 ± 0.3694.65 ± 0.3883.86 ± 15.3694.63 ± 0.3694.65 ± 0.31
Average81.64 ± 2.5382.00 ± 2.7682.35 ± 1.0482.47 ± 0.7381.77 ± 2.5972.50 ± 8.7082.87 ± 0.8683.40 ± 0.69
* Bolded indicates that the model achieves the best performance on this dataset.
Table 6. The number of reduction attributes.
Table 6. The number of reduction attributes.
IDNameDTRS-EFDTRS-SMDNSSPDTRS-EFSPDTRS-SMDNSNDTRSFDTRSFN3WDFNζDTRS
1Atoms5.5 ± 1.86.4 ± 1.86.8 ± 0.47.8 ± 0.45.2 ± 1.24.5 ± 1.31.2 ± 0.42.0 ± 0.0
2Australian11.2 ± 1.912.2 ± 3.011.3 ± 0.613.0 ± 0.112.7 ± 1.44.7 ± 2.710.8 ± 0.68.4 ± 0.5
3Breast7.7 ± 2.78.5 ± 2.09.0 ± 0.09.0 ± 0.19.0 ± 0.23.6 ± 1.08.0 ± 0.06.9 ± 0.6
4Cleve9.8 ± 0.512.5 ± 2.49.9 ± 0.513.0 ± 0.210.2 ± 2.64.1 ± 2.77.1 ± 1.63.0 ± 0.0
5Heart7.1 ± 2.012.1 ± 2.67.9 ± 0.612.7 ± 0.410.0 ± 3.04.8 ± 3.16.8 ± 1.33.0 ± 0.0
6Iris3.8 ± 0.73.6 ± 0.83.6 ± 0.54.0 ± 0.14.0 ± 0.01.7 ± 1.32.6 ± 0.51.0 ± 0.0
7Phishing8.3 ± 2.38.8 ± 1.49.0 ± 0.09.0 ± 0.08.8 ± 1.41.1 ± 0.29.0 ± 0.09.0 ± 0.0
8Saheart8.7 ± 1.68.5 ± 1.99.0 ± 0.09.0 ± 0.09.0 ± 0.06.4 ± 2.84.8 ± 0.42.6 ± 0.3
9Seismic6.6 ± 4.810.9 ± 6.34.0 ± 0.213.4 ± 0.913.8 ± 0.58.3 ± 0.72.0 ± 0.01.0 ± 0.0
10Vote8.5 ± 0.68.5 ± 0.68.5 ± 0.58.5 ± 0.58.5 ± 0.81.1 ± 0.38.5 ± 0.58.1 ± 0.3
Average7.7 ± 1.99.2 ± 2.37.9 ± 0.39.9 ± 0.39.1 ± 1.14.0 ± 1.66.1 ± 0.54.5 ± 0.2
* Bolded indicates that the model achieves the best performance on this dataset.
Table 7. Information of the telemetry parameters.
Table 7. Information of the telemetry parameters.
IDAttributeRate RangeData TypeUnit
a1Duty cycle0–1Continuous-
a2Bus current13.5–17.3ContinuousA
a3Shunt current5.3–12.4ContinuousA
a4Battery current3.6–19.4ContinuousA
a5Output power1070–1090ContinuousW
a6Battery pressure2.0–5.4ContinuousMPa
a7Battery quantity54.3–71.2ContinuousAh
a8Status word−1 0 1Discrete-
a9Bus voltage40.5–43.1ContinuousV
a10Battery voltage33.0–40.5ContinuousV
Table 8. Different scenarios for faults in satellite power system.
Table 8. Different scenarios for faults in satellite power system.
ScenarioFault NameScenarioFault Name
0----6F1-F3
1F17F1-F4
2F28F2-F3
3F39F2-F4
4F410F1-F2-F3
5F1-F211F1-F2-F4
Table 9. The results of the training process.
Table 9. The results of the training process.
Fault NameThe Output Attribute SubsetAverage Value of the Data for Attribute a3
F1a5-
F2a7-
F3a2, a3, a97.46
F4a2, a3, a913.47
Table 10. The results of the testing process.
Table 10. The results of the testing process.
Fault NameThe Output Attribute SubsetJaccard Similarity CoefficientAverage Value of the Data for Attribute a3Matching Result
F1F2F3F4
F1-F2a5, a70.500.5000-F1-F2
F1-F3a2, a3, a5, a90.2500.750.756.63F1-F3
F1-F4a2, a3, a5, a90.2500.750.7512.62F1-F4
F2-F3a2, a3, a7, a900.250.750.757.47F2-F3
F2-F4a2, a3, a7, a900.250.750.7513.46F2-F4
F1-F2-F3a2, a3, a5, a7, a90.200.200.600.606.63F1-F2-F3
F1-F2-F4a2, a3, a5, a6, a7, a90.170.170.500.5012.63F1-F2-F4
Table 11. The results of the simultaneous-fault diagnosis.
Table 11. The results of the simultaneous-fault diagnosis.
AlgorithmAccuracy/Subset AccuracyHamming LossPrecisionRecallF1
FNζDTRS–FMS100.0 ± 0.0-100.0 ± 0.0100.0 ± 0.0100.0 ± 0.0
Binary Relevance82.34 ± 5.394.43 ± 1.35100.0 ± 0.093.88 ± 1.8796.28 ± 1.14
Classifier Chain65.78 ± 8.939.30 ± 2.63100.0 ± 0.086.42 ± 4.4591.31 ± 3.07
Calibrated Label Ranking0.0 ± 0.032.14 ± 7.11100.0 ± 0.045.24 ± 0.061.90 ± 0.0
ML-KNN0.0 ± 0.032.15 ± 0.099.98 ± 0.045.23 ± 7.1161.89 ± 0.0
ML-DT0.0 ± 0.032.35 ± 0.3199.59 ± 0.6245.03 ± 0.3161.63 ± 0.42
* Bolded indicates that the model achieves the best performance on this metric.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tao, L.; Wang, C.; Jia, Y.; Zhou, R.; Zhang, T.; Chen, Y.; Lu, C.; Suo, M. Simultaneous-Fault Diagnosis of Satellite Power System Based on Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set. Mathematics 2022, 10, 3414. https://doi.org/10.3390/math10193414

AMA Style

Tao L, Wang C, Jia Y, Zhou R, Zhang T, Chen Y, Lu C, Suo M. Simultaneous-Fault Diagnosis of Satellite Power System Based on Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set. Mathematics. 2022; 10(19):3414. https://doi.org/10.3390/math10193414

Chicago/Turabian Style

Tao, Laifa, Chao Wang, Yuan Jia, Ruzhi Zhou, Tong Zhang, Yiling Chen, Chen Lu, and Mingliang Suo. 2022. "Simultaneous-Fault Diagnosis of Satellite Power System Based on Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set" Mathematics 10, no. 19: 3414. https://doi.org/10.3390/math10193414

APA Style

Tao, L., Wang, C., Jia, Y., Zhou, R., Zhang, T., Chen, Y., Lu, C., & Suo, M. (2022). Simultaneous-Fault Diagnosis of Satellite Power System Based on Fuzzy Neighborhood ζ-Decision-Theoretic Rough Set. Mathematics, 10(19), 3414. https://doi.org/10.3390/math10193414

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop