Active fault-diagnosis method using adaptive allocator and fault-tolerant adaptive control system design

A new approach to the active fault diagnosis (AFD) for input redundant plants is presented in this paper. A test signal which helps the diagnosis is injected to the plant in addition to nominal control input in the AFD typically. One feature of the proposed AFD is adaptive allocation of the test signal. The adaptive allocator which distributes the test signal and injects it to the redundant actuators is introduced in this paper. A fault-diagnosis (FD) system that estimates fault location and its magnitude from adaptive parameter of the adaptive allocator is constructed with a simple neural network (NN) model. Furthermore, a A fault-tolerant adaptive control system which includes the adaptive test signal allocator is designed based on the model reference adaptive control (MRAC) technique. The adaptive laws of the adaptive allocator and controller are derived by using a suitable Lyapunov function. By performing experiments using a two-input redundant plant, we show the effectiveness and applicability of the proposed AFD and MRAC system.


Introduction
Many machine-based systems consist of actuators, sensors, computers, and mechanical elements. The actuators play an important role in manipulating the plant states [1]. To improve the actuation and/or fault tolerance performance, input-redundant systems have been developed which have one or more additional actuators [2][3][4]. Based on the fault tolerance of redundant plants, it is possible to eliminate the faulty influence on the system. For instance, a self-repairing control system converts faulty actuators into healthy ones, and fault accommodation is also performed with the remaining healthy actuators in adaptive control systems [5][6][7][8][9]. From the perspective of system safety and maintainability, there is an increasing need to obtain information about the system health during operation [10,11]. Therefore, a number of studies that address the fault-diagnosis (FD) problem have been developed [11][12][13][14][15][16]. Most of the fault-diagnosis methods that have been developed employ the passive approach. In the passive approach, by observing the system behavior using the available signals, it is possible to monitor the plant without affecting the system states.
On the contrary, the active fault diagnosis (AFD) method has been proposed, where an auxiliary signal (test signal) is used to improve the precision and speed of the diagnosis [12,17]. In contrast to the passive fault diagnosis method, relatively few studies have focused on the active method [17]. In the literature [18], Takahashi proposed an algorithm for fault detection using a weak test signal for a system that has a redundant actuator. In the FD system which was developed by Ducard et al., the supervision module uses a test signal to help the FD process [19]. However, when the test signal is used for diagnosis, it also decreases the performance of the control system [20,21].
On the other hand, it is known that it is possible to inject an arbitrary signal into an input-redundant plant without affecting the plant states [22][23][24]. In the literature [22], Elgarsma et al. proposed a method to identify the plant parameters by using the test signal without affecting the attitude of the aircraft, based on the input redundancy. This idea has been applied to the AFD technique, and there are several related studies of this application in the field of aircraft control [25][26][27][28]. Bateman et al. proposed an AFD method in which the test signal is inputted into the null-space of the input coefficient matrix [25]. In the FD system reported in literature [25], the test signal affects the plant states only when a fault occurs. A supervisor performs the fault diagnosis using the test signal response that becomes visible immediately after the fault occurrence. The AFD method with input redundancy is described schematically in "Outline of conventional AFD" section. One drawback of this method is that the null-space injection requires the precise plant parameters related to the input coefficient matrix. Moreover, it is not clear how the system stability is influenced by the test signal during the fault occurrence and its compensation. In the literature [28], Boskovic proposed a control system design that ensures the plant stability under the test signal injection, even if an actuator has a fault. However, the study in [28] did not refer to a specific AFD method.
We propose a new AFD method with input redundancy in this paper based on the, and design a control system proposed in the literature [29] to implement it in this paper. In the literature [29], we have designed an adaptive control system that guarantees reference model tracking while applying arbitrary bounded test signals for the input redundant system with unknown parameters. A feature of the proposed AFD method is adaptive allocation of the test signals in the control system. One adaptive test signal allocator is placed for a pair of redundant actuators. Therefore, the number p of adaptive allocators is p = m!/(2!(m − 2)!) when the controlled plant has m redundant actuators. The effectiveness ratio of a pair of redundant actuators can be estimated from the adaptive parameter in the corresponding adaptive allocator. We define the ratio as a feature value R for the proposed AFD process. The fault location can be estimated from an increase or decrease of the feature value R. And fault magnitude also can be estimated from the convergent value of R. A simple NN (neural network) model that describes the relation between the feature value R and both the estimated fault location and its magnitude is introduced in the proposed AFD system. The AFD system design is described in "Design of AFD system based on adaptive parameters" section. The novelty of the proposed AFD method is to place the adaptive allocator for each pair of the redundant actuators. The fault diagnosis for the individual actuator can be achieved with such the placement of adaptive allocators. The feature of the proposed AFD compared with the conventional one is described in "Proposed AFD using adaptive allocator" section. Furthermore, a fault-tolerant adaptive control system which includes the adaptive test signal allocator is designed for the faulty model in "Plant representation with actuator fault" section. In "Plant representation with actuator fault" section, we introduce the plant representation with actuator fault. The proposed control system design for the faulty model is based on the model reference adaptive control (MRAC) technique. In "Fault-tolerant adaptive control design with adaptive allocator" section, we show the derivation of the adaptive law for the controller and allocator using a suitable Lyapunov function [29]. The adaptive laws of the adaptive allocator and controller are derived by using a suitable Lyapunov function, in "Fault-tolerant adaptive control design with adaptive allocator" section.
We verify the effectiveness and applicability of the designed control system and proposed AFD system, through the experiments using a two-input redundant machine-based plant, in "Experimental example" section. The conclusions are described in "Conclusion" section.

AFD with input redundancy
Outline of conventional AFD Figure 1 illustrates the schematic of the AFD with input redundancy in a two-input redundant plant. In this figure, the test signal τ is injected through a pre-designed allocator α. This allocator α is chosen such that its direction is orthogonal to the plant input coefficient vector, b. When no actuators have any fault, only the part of control effort caused by the test signal is canceled out in the total control effort, E := b T u , and the states and output of the plant are not affected by the test signal. In the case where an actuator has a fault, it is seen that the coefficient vector b changes to a faulty one b ′ in the faulty model representation. Based on the assumption that the faulty vector b ′ satisfies b ′ T α = 0 , the test signal τ will be not canceled out, and its effect appears on the plant output. In other words, the test signal has a symptom that is caused by the actuator fault on the plant output. The conventional AFD is performed from the symptom, available signals in the control system, and prior knowledge about the plant. However, it also means that the plant output is affected by the test signal when a fault occurrs. There is a trade-off between the FD performance and the stability of the plant output, in the conventional AFD method.

Proposed AFD using adaptive allocator
The proposed AFD method is not based on the fault symptom on the plant output in contrast to the conventional one. In proposed AFD (Fig. 2), the test signals are injected consistently through the adaptive allocator to the redundant actuators. The proposed AFD is performed from the behavior of adaptive parameter of the allocator and available signals in the control system. The stability of the plant output is ensured by the suitable adjusting of the adaptive allocator and adaptive controller in the proposed method. Both to guarantee the stability of plant output and to perform the AFD can be achieved simultaneously, in the proposed method. Therefore, the proposed method can perform AFD for the input redundant plant without obtaining parameters of b.

Plant representation with actuator fault
We consider the following linear plant: where x(t) ∈ Rn and u(t) ∈ Rm are the plant output and input, respectively. Matrices A, B have appropriate dimensions, and can be represented as follows: where 0 denotes a zero matrix which has dimension (n − 1) × m , a 2 ∈ Rn and b ∈ Rm are unknown coefficient vectors, but that A 1 ∈ R (n−1)×n is a known submatrix, and the signs of each component in vector b are known. We assume that the known matrix A 1 makes the following square matrix A 0 stable, by suitable selection of the vector a m ∈ Rn.
For example, the (n − q)-th order systems that have q(< n) integrators, or n-th order systems with relative degree n can be considered as plants that satisfy the above conditions. In these systems, matrix A 1 becomes known from its specific structure, and the plant parameters that may be unknown are in the vectors a 2 and b.
The fault on the actuator is represented by u(t) as where u c ∈ Rm is the controller output, and ū ∈ Rm is the fault input signal, which depends on the type of fault described later. Λ ∈ R m×m is an unknown fault-magnitude matrix. We define each element i , (i ∈ {1, 2, · · · , m}) in diagonal matrix Λ as the effectiveness of the actuator.
In this study, we focus on the following two types of actuator fault.

Loss of effectiveness
The loss of effectiveness (LOE) of the i-th actuator is represented with i as: where t f (i) ∈ R, (t f (i) > 0) is the unknown time at which the fault occurs. i ∈ {1, 2, · · · , m} is the unknown index which denotes the faulty actuator, d i ∈ R, (0 ≤ d i < 1) is also an unknown parameter that represents the effectiveness of the faulty actuator. For instance, t f (2) = 50[s] , d 2 = 0.6 means that the second actuator loses its effectiveness of 40% after 50[s].

Stuck actuator output
The stuck fault of the actuator may be caused by the freezing of the electronic motor driver and solidified valve for instance. The stuck fault is represented as follows.
From the above representation, we note that the effectiveness of the stuck actuator is always zero. Therefore, the case where the actuator loses its effectiveness completely (d i = 0,ū = 0) can be considered as a special case in the stuck fault. After this section, we apply the word perfect LOE fault to the fault represented by d i = 0 , and we also apply the word partial LOE fault to the fault represented by 0 < d i < 1.

Fault-tolerant adaptive control design with adaptive allocator
In this section, we describe the design of an fault-tolerant adaptive control system with an adaptive allocator for testsignal injection and actuator faults [29]. For the system design, we considered the following assumptions: 1. Each actuator may have a fault that is expressed by either (i) or (ii) in "Plant representation with actuator fault" section. 2. At least one actuator has not perfect LOE fault. 3. All of the plant outputs are available to generate the control input signal.
The control system design is shown in the following three subsections. First, a stable reference model and adaptive controller are introduced in "Output tracking error equation" section. The output tracking error equation which represents the difference of plant output and reference model output is derived using the reference model and the adaptive controller. In "Adaptive allocator design for test signals" section, the adaptive test signal allocator is introduced. Finally, the adaptive laws of the adaptive test signal allocator and adaptive controller are derived from a suitable Lyapunov function in subsection 4.3.

Output tracking error equation
Equation (1) can be represented using (2) and (4) as: where x 1 ∈ R n−1 and x 2 ∈ R are parts of the output vector x.
We introduce a stable reference model with known matrix A 1 , where x m ∈ R n is the output of the reference model, and x m(1) ∈ R n−1 and x m(2) ∈ R are the components of reference model output vector. r ∈ R is an arbitrary bounded input signal for the reference model, and b m ∈ R is an arbitrary coefficient constant. The vector a m ∈ R n is selected such that the matrix A 0 in (3) becomes stable.
Boskovic proposed the following controller to design a stable MRAC system for such the plant and the reference model which are represented in (7) and (9), in the literature [5].
where L ∈ R m is the constant allocation vector. θ 1 ∈ Rn , θ 2 ∈ R and ξ ∈ R are the adaptive parameter vector and adaptive parameter, respectively. The adaptive laws of the adaptive parameters θ 1 , θ 2 , and ξ have been presented in [5].
In this paper, we propose the following controller which is extended from (10).
where L ∈ Rm is a constant vector that satisfies b T ΛL > 0 . θ 1 ∈ Rn , θ 2 ∈ R and ξ ∈ R are also the adaptive parameter vector and adaptive parameters, respectively. Each element τ q (t),(q = 1, 2, · · · , p) in vector τ is an arbitrary bounded test signal, and they are independent of each other. α(Φ(t)) ∈ R m×p is the adaptive allocator matrix for the test signal, and Φ(t) ∈ Rp is the adaptive parameter vector.
Substituting (11) into (7), we obtain We define the tracking error e = x − x m , and derive its time derivative.
From (9) and (12), the first row of (13) can be written as: We calculate the second row of (13) as: We introduce the following variables where we set the variable κ as: From the assumption, it satisfies κ > 0. Using these variables, we rewrite (15) as: Hence (13) can be represented as:

Adaptive allocator design for test signals Basic idea for adaptive allocator
We now consider the following two-input redundant plant.
The part of the input vector u that is orthogonal to the plant coefficient vector b := [b 1 , b 2 ] T , becomes 0 in the step to generate the control effort E = b T u . In other words, the orthogonal part of the input vector for the plant-coefficient vector does not affect the plant output x . Therefore, aiming for null-space injection, we design the controller output as: where the function f (x, x m , r) generates a nominal control input, and where x , x m , and r are the feedback of the plant output, reference model output, and reference input, respectively. τ is the test signal that should be canceled out on the plant output, α(θ) is the adaptive (15) allocator, and where θ is adaptive parameter. Figure 3 shows an allocator vector α and coefficient vector b in a two-dimensional (2D) plane. As shown in this figure, the allocator vector α has a constant norm, and can change its direction. If the b⊥α(θ ) can be satisfied by adjusting the parameter θ , then the effect on the plant outputs caused by the test signal will be removed.
In this section, we design the test signal allocator for the m-input plant based on the 2D rotary allocator. In the allocator design, m redundant actuators are divided into all possible pairs, and the 2D allocators are assigned to these. Therefore, the number p of test-signal allocators and test signals is: For ease of explanation, we begin from the three-input plant in the next section.

Allocator design for three-input plant
In the case of ( m = 3 ), let the input coefficient vector b Corresponding to these three pairs, we introduce the 3rd order test signal vector as: Then, we choose an allocator α (3) as follows.
Using these definitions, the coefficient part of the test-signal term b T Λα in (19) can be expanded with Λ (3) = diag{ 1 , 2 , 3 } as: We now define the following 2D vectors, Substituting the elements in (26) by the inner products of these vectors, The Eq. (29) indicates that the coefficients from test signals to the control effort can be represented using 2D allocators ᾱ 1 , ᾱ 2 , and ᾱ 3 .

Allocator design for m-input plant
According to the result in "Allocator design for threeinput plant" section, an allocator matrix can be chosen as follows to assign the 2D allocators to the m redundant actuators.
The test signal vector is: From these difinitions, we expand the test signal term in a similar way in the case of (m = 3) as: In addition, each element can be represented by the inner product of 2D vectors: where the 2D vectors are defined as follows: Figure 4 shows vectors b q and ᾱ q on a 2D plane, where q ∈ {1, 2, . . . , p} . By using an angle φ q that is defined by the vectors shown in Fig. 4, (33) can be represented as: where the norm condition ᾱ q = 1, (∀q ∈ {1, 2, . . . , p}) is used. Furthermore, we define the φ q as follows.
b T Λα = b 1ᾱ1 cosφ 1b2ᾱ2 cosφ 2 · · · · · ·b p−1ᾱp−1 cosφ p−1bpᾱp cosφ p = b 1 cosφ 1b2 cosφ 2 · · · · · ·b p−1 cosφ p−1bp cosφ p Using this definition, we rewrite (36) as: Therefore, the test signal term can be represented as: where we used the denotations of: When the actuators have no fault, matrix N is positive definite from the assumption that all of the elements in vector b are non-zero. Even if (m − 1) actuators have partial/perfect LOE fault, then the matrix is positive or semipositive definite.

Adaptive law design
Equation (19) is rewritten using (39) as: where µ and b 0 are as follows.
From the assumption that matrix A 0 is stable, the positive definite symmetric matrix solution P ∈ R n×n exists in the Lyapunov equation: where Q ∈ R n×n is a positive definite symmetric matrix. We now define a scalar function V with matrix P as: where S ∈ R p is defined as follows.
� sin 1 2φ 1 cos 1 2φ 1 sin 1 2φ 2 cos 1 2φ 2 · · · · · · sin 1 2φ p cos 1 (52) Substituting (54) into (56), we obtain: It is clear that V ≤ 0 , because matrix P is positive definite, as given in (45). Therefore, the scalar function V is positive definite and non-increasing for time t , and the signals e, �θ 1 , �θ 2 , �ξ are bounded, as is the case with adaptive parameters θ 1 , θ 2 , ξ from (16). Based on the fact that V is positive definite and non-increasing, as well as (56), there exists which is a finite value, so e ∈ L 2 . Moreover, the plant output x is bounded, and from that, reference model output x m is bounded in e = x − x m . From (42) and the fact that signals �θ 1 , �θ 2 , �ξ , x, r and the test signal term S T N τ are bounded, ė is also bounded. By using Barbalat's lemma, we obtain e → 0 at time t → ∞ because e ∈ L 2 and its derivative ė are bounded. We now discuss the boundedness of the adaptive parameter vector Φ . From (54), Φ can be given as: where Φ(0) is an arbitrary bounded initial value vector. Hence, the test signal and the error signal e are bounded, and the first term on the right-hand side of the above equation is bounded. From these facts, it is understood that the vector Φ is bounded. We summarize the designed MRAC system, as shown in Fig. 5.

Design of AFD system based on adaptive parameters
In this section, we present a new approach for AFD with input redundancy, and the design process of the AFD system is discussed.

Basic idea
In (34), we introduced the vector b q , which corresponds to a pair of redundant actuators. Figure 6 shows b q on the 2D plane. We assume that the vector b q corresponds to a pair of the i th and j( = i) th actuators. Figure 6a shows the fault-free case, so each element of b q is equal to b i and b j .
We consider the case where either the i th or j th actuator has a fault. As shown in Fig. 6b, the nominal vector b q changes to faulty vector b′ q after the fault occurrence. This change in the vector differs depending on the faulty location and its magnitude. Therefore, it is possible to estimate the faulty location and its magnitude if the vector b′ q is known, and if its change can be monitored precisely.
However, the plant parameters b i , b j and faulty parameters i , j are generally unknown, so b q cannot be obtained directly. Therefore, our idea is to monitor only the direction of b q . Figure 7 again shows the changes of vector b q from the perspective of its direction. arg(·) refers to the argument of the vector, and N (·) is the null-space of the vector. From Fig. 7b, it is clear that the changes of direction also differ depending on the fault location and its magnitude.

Feature value for ratio of control effectiveness
To monitor the direction of vector b q indirectly, we used the adaptive parameter φ q , which refers to the argument of allocator vector ᾱ q . In this section, we introduce a feature value for fault diagnosis from the adaptive parameter φ q .
To discuss the feature value, we assume that the following are satisfied on the designed MRAC system in 1. When sufficient time has passed from the time the control is started, tracking error e converges to 0 , and its time derivative ė also converges to 0. 2. We assume that the reference input r is a step-like signal, so ṙ = 0 is satisfied when sufficient time has passed. 3. A fault that is represented in (5) or (6) may occur in one actuator. In other words, only one actuator may have a fault, and the remaining ones are fault-free at every moment.
Next, based on these assumptions, we show that the adaptive allocator ᾱ q indicates an orthogonal direction to vector b q .
When sufficient time has passed from the start of the control operation, the following equation is held from the assumption 1) and (42).
Then, based on assumption (1) and (54), each time derivative of the adaptive parameter is as follows.
The time derivative of (58) can be calculated from the above results and assumption (2) as follows: For the test signals in which its time derivatives are bounded and independent of each other, S = 0 is needed to satisfy the above equation. Therefore, the adaptive parameter satisfies, where n q ∈ Z is an integer corresponding to index q . From (61) and (37), we obtained the following relations.
We found that the allocator ᾱ q is orthogonal to the vector b q when the above equation is satisfied. Figure 8 shows the plot of ᾱ q which corresponds to vector b q . When ᾱ q is orthogonal to the vector b q , vector b q can be expressed as follows.
(62)  Using this expression, we can calculate the ratio R q , ∀q ∈ {1, 2, . . . , p} of the first element in the vector b q as: This result means that it is possible to obtain the ratio of i th control effectiveness in a pair of i th and j th redundant actuators using only the adaptive parameter φ q . Moreover, from the assumption (3), increasing or decreasing of a ratio R q can be considered as the fault in i th or j th actuator.

FD system design
We calculated the proposed feature value R in (64) from the adaptive parameter, and the discussion for R has been based on the steady state of the system. In actual case, the adaptive parameter based feature value R may be moved in response to the disturbance, and there exists the transient state during system operation. Therefore, we should utilize other available signals in addition to the feature value R , in the FD process.
We constructed a classifier of an actuator fault using the machine-learning technique. Therefore, the proposed AFD system requires a preliminary experiment for the learning.
The overview of the proposed MRAC system with the AFD module is illustrated in Figs. 9 and 10a shows the schematic of the proposed AFD system. This AFD system outputs the fault diagnostic information f by monitoring the signals e, r,ẋ m and φ . The diagnostic information f has two elements as follows: where î , î ∈ {0, 1, . . . , m} is the estimated fault location, and ˆ i is the estimate of the fault magnitude of the îth actuator. When it is seen that no actuator has any fault, then the î indicates î = 0.

NN-based fault estimator
The core of the fault-diagnosis module is the NN model shown in Fig. 10b. Because the nonlinear relationship between fault magnitude and the feature value, we use the NN model which is known to be effective in approximating nonlinear functions as a failure estimator in this paper. In this study, we used the extreme learning machine (ELM) technique [30,31] to learn the NN model. ELM is a machine-learning method that was proposed by Huang et al. One of the major advantages of ELM is that it does not perform the repeated calculation for the learning. Another advantage is the simple structure of the NN model, as shown in Fig. 10b. NN model output Y = [y 1 , y 2 , . . . , y 2m ] T , Y ∈ R 2m is given as: T is the output weighting vector from the hidden layer to the output layer. z, (z ∈ Z, z > 2m) is the number of nodes in the hidden layer, and it is one of the design parameters of ELM. Input signal X ∈ Rv of the NN model is represented as: where k ∈ Z refers to the step number of discrete time.
The order v of vector X becomes v = (n + p + 2) from (67). In this study, we used the following sigmoid function as the activating function h(X).
(66) where A s(i) ∈ R v and B s(i) ∈ R v are the parameters that characterize the activating function. These parameters are randomly determined once, and are not reconstructed. The weighting vector β i is calculated as follows: where C is a positive scalar that is given by the designer. H and Ȳ i are respectively the output matrix and correct data vector corresponding to the learning data X , and they are given as follows: where n t is the number of samples that is included in the learning data, so H ∈ R n t ×z ,Ŷ i ∈ R n t .
In the ELM learning, the weighting parameters A s(i) , B s(i) among the input layer and hidden layer are not reconfigured, and the weighting parameters β i can be determined uniquely using (69). Therefore, if the order of the learning data is not large, the learning step can be completed quickly, and the computation time is short FD algorithm Figure 11 shows the FD algorithm that outputs the FD result from the output of the NN model. First, the algorithm selects the maximum of y j ∈ {1, 2, . . . , m} in Y , and names it i . If the sign of y i is positive, then the i-th actuator is assumed to be faulty and i is assigned to f 1 . Furthermore, y m+i , which refers to the estimated fault magnitude of the i th actuator, is assigned to f 2 , (i.e., f 2 = y m+i ).

Experimental setup and conditions
In this section, we perform the verification experiments of the proposed AFD system and the fault-tolerant adaptive control system. The main purpose of this experiment is to confirm the applicability of the proposed method in an actual plant. The experimental equipment is shown in Fig. 12. Using the experimental equipment shown in Fig. 12, we confirm the effectiveness and applicability of the proposed fault-tolerant adaptive control system and AFD system. Table 1 lists the specifications of the experimental plant. In this plant, one rotary shaft is driven by two actuators with a timing belt. Each actuator consists of a DC motor and a motor driver, and it generates a torque in response to individual command voltages. The experimental mechanical plant can be modeled as the following two-input first-order system.
where b p1 , b p2 > 0 . We provide the stable reference model as: And the constant vector L in (11) is, We chose the following periodic signal as the test signal for AFD. Equation (54) implies the adaptation of the parameter Φ proceeds by the tracking error which includes a test signal response of the plant. Thus, if a prior information about approximate dynamic characteristics of the plant is available, a test signal with a frequency that is hard to be decayed in the plant should be selected. A relatively low frequency signal was chosen for this experimental test signal due to the low-pass characteristics of the experimental plant. In addition, the test signal consists of sine waves with different periods and has a certain magnitude, for the speed of adaptation. Each adaptive parameter has the following initial value, and the adaptive gain is as follows: (75) τ (t) = 0.05(sin (0.8π t) + sin (1.5π t))  In the experiment, we implemented the actuator fault using the following equation.
where the fault parameters 1 , 2 and signal ū are calculated using (5) or (6). The design parameters of ELM learning in the FD system are as follows. (76)

Preliminary experiment
In the experimental plant, there is only one way of making a pair of redundant actuators. Therefore, the feature (79) Z = 46, C = 2 36 .
value that is introduced in (64) corresponds to a scalar R 1 . When the actuators have no fault,R 1 has an unknown value in 0 < R 1 < 1 , from the discussion in "Feature value for ratio of control effectiveness" section. As the fault magnitude of the first actuator becomes large, R 1 decreases toward 0 . On the other hand, as the fault magnitude of the second actuator becomes large, R 1 increases toward 1. Therefore, it can be seen that adjusting the R 1 such that R 1 = 0.5 is satisfied in the fault-free case has an advantage for fault diagnosis. For such an advantage, we give a weighting matrix W to the test signal allocator.
From (72) and (78), the control effort E τ that the test signal generates in the plant can be represented as: where we used Λ = I because we now consider the faultfree case. Equation (81) shows that the ratio wb p2 b p1 of elements in the vector b p can be adjusted by manipulating the weighting constant w . Applying the weighting matrix W to the test signal allocator corresponds to transforming the vector b p with positive definite diagonal matrix W . As a result, we note that the discussion about the MRAC system in Sect. 4 does not lose its generality.
We performed the preliminary experiments to determine the weight parameter w . We began the preliminary experiment with a weight parameter w = 1 , and (80) other experimental conditions are as shown in "Experimental setup and conditions" section. By performing the preliminary experiment, we obtained a desirable control result. We also determined that φ s = −7.39[rad] is the experimental convergence value of the adaptive parameter φ from the fault-free preliminary experiment. Figure 13 shows the adaptive allocator α(φ) with φ = φ s . The plant vector b p can be expressed as follows based on the assumption that the vector α(φ s ) is orthogonal to the vector b p .
where we assume that each element of the vector b p is positive. From this result, we determine the weight parameter as w = − cos φ s sin φ s = 1.998 for the main experiment.

Main experimental result
Because of space limitations, we present only a part of the experimental results here. Figure 14 shows an experimental result. In this case, the second actuator has a partial LOE fault with magnitude 2 = 0.4 at 150 (s) from the start of the experiment. In the Fig. 14, the upper plot shows the plant output and reference model output, while the middle and lower plots show the first input u c1 and second input u c2 of the plant, respectively. With respect to the period before 150(s), we see that the plant output converges to the model output. Then, the plant inputs u c1 , u c2 are consistently vibrated even after the plant output converges to the model output. At 150 (s), the second actuator has a partial LOE fault. We confirmed that the reference model tracking is maintained after the fault occurs.    Figure 15 shows the AFD result in the experiment of Fig. 14. Figure 15a shows a discrete-time series of elements in the ELM input vector X , and (b) shows a discrete-time series of the FD result. In plot (b), we can see that f 1 is almost 0 when the actuator has no fault. About 200 steps after the fault occurrence, f 1 is 2, and then f 2 gradually tends to 0.4.
An accuracy rate of f 1 , which indicates the fault detection and isolation result, is summarized in Table 2. The content of the table was obtained from six experiments in total, and it includes the result of Fig. 15. The number of samples of each experimental dataset which was used to plot Table 2 and to learn ELM is presented in Table 3. Table 3 shows the number of datasets used for training and verification for each fault condition. From this table, it can be seen the data corresponding to the condition in Fig. 15 ( = 0.4 ) does not exist in the learning data. Nevertheless, the fault diagnosis has been achieved as shown in Fig. 15b. This also shows the high generalization performance of the proposed FD system using ELM. From Table 2 and Fig. 15, we confirmed that the proposed FD system can estimate the fault occurrence and its location with high accuracy.
We performed an experiment on a simple two-input machine-based plant to verify the proposed method in this section. In recent years, systems with redundant actuators, such as those seen in drones and electric/hybrid-electric vehicles [32], have attracted much attention. The proposed FD method and the fault-tolerant control design also have the potential to be applied to such the input redundant systems.