Fuzzy Filter-Based FDD Design for Non-Gaussian Stochastic Distribution Processes Using TS Fuzzy Modeling

This paper studies the fuzzy modeling problem and the fault detection and diagnosis (FDD) algorithm for non-Gaussian stochastic distribution systems based on the nonlinear fuzzy filter design. Following spline function approximation for output probability density functions (PDFs), the T-S fuzzy model is built as a nonlinear identifier to describe the dynamic relationship between the control input and the weight vector. By combining the designed filter and the threshold value, the fault in T-S weight model can be detected and the stability of error system can also be guaranteed. Moreover, the novel adaptive fuzzy filter based on stochastic distribution function is designed to estimate the size of system fault. Finally, the simulation results can well verify the effectiveness of the proposed algorithm for the constant fault and the time-varying fault, respectively.


Introduction
In order to improve the stability and the security of system, the fault detection and diagnosis (FDD) algorithm for the complex systems has been an important part in the field of control engineering.Many significant approaches have also been presented and applied to practical processes successfully (see [1][2][3][4] and references therein).On the other hand, a series of modeling and control strategies which control the shape of output probability density functions (PDFs) for non-Gaussian stochastic processes have received considerable attention (see [5][6][7][8]).However, most of the existing FDD results for stochastic systems were only concerned with Gaussian variables, where mean or variance was the objective for optimization and control [9].Since 2000, the FDD algorithm for non-Gaussian stochastic distribution systems has begun to be discussed based on filter theory [10].In these results, the modeling problem is often ignored and only linear model or common nonlinear model is considered in FDD (see [10][11][12]).
In recent years, the well-known T-S fuzzy model was viewed as a popular and powerful modeling tool since it is a powerful solution that bridges the gap between linear and nonlinear control systems (see [13][14][15]).By introducing a family of fuzzy IF-THEN rules that represent local linear input-output relations of the system, many complex nonlinear models, such as descriptor systems [16], stochastic systems [17], and time-delay model [14], can be described or approximated by T-S fuzzy modeling.Meanwhile, some typical control problems, including dynamic tracking control [14], sliding-mode control [18], and filter design [19,20], have also been considered through T-S fuzzy modeling.
In this paper, we provide a novel FDD approach for non-Gaussian stochastic distribution systems.Based on Bspline approximation and T-S fuzzy modeling, the concerned FDD problem of the dynamic non-Gaussian systems can be transferred into a special nonlinear FDD problem for deterministic T-S weight dynamics.Instead of common nonlinear observer or filter in [10][11][12], the nonlinear fuzzy filter and the adaptive filter are conducted by involving the measured output PDFs with the T-S fuzzy weight dynamics.Moreover, by optimizing a series of linear matrix inequalities (LMIs), the fault existing in the stochastic processes can be detected with a defined threshold and the satisfactory estimation value for the size of fault can also be guaranteed.It is noted that the T-S fuzzy model is first applied into the FDD for stochastic distribution systems which solves the nonlinear modeling difficulty in previous results [10][11][12].This represents a significant extension of the previous results that only common linear/nonlinear weight dynamic models are considered and has also an independent significance in the field of fuzzy FDD.
In this paper, for a square matrix , we denote that sym() =  +   .The identity and zero matrices are expressed by  and 0, respectively.For a symmetric matrix , the notation  > (≥)0 is used to denote that  is positive definite (positive semidefinite).The case for  < (≤)0 follows similarly.Moreover, for a vector V(), define Euclidean norm by ‖V()‖ = V  ()V().

Problem Formulation and T-S Fuzzy Modeling
For a complex non-Gaussian stochastic process, we denote that () ∈   is the control input, () ∈ [, ] is the stochastic output of the concerned system, and  represents the fault to be detected and is supposed as a constant vector.
The conditional probability of output () lying inside [, ] can be defined as follows: where (, (), ) stands for the output PDF with fault term under the effect of the control input ().
Due to the PDF constraint condition ∫   (, (), ) = 1, only  − 1 weights are independent.Thus, the output PDF can be further expressed as where Similarly, with [7], the nonlinear term ℎ(()) satisfies the following equality: where It is obvious that the nonlinear term ℎ(()) satisfies the Lipschitz condition within its operation region; that is, for any  1 () and  2 (), there exists a known matrix  1 such that In the following, we will find the dynamic relationship between the control input () and the weight vectors ((), ), which corresponds to a further modeling procedure.It is well known that the T-S fuzzy model is a powerful solution for identifying complex nonlinear dynamics by a blending of some local linear system models.Compared with those results that only consider common linear/nonlinear weight dynamic model (see [10][11][12]), we will use the T-S model to describe the nonlinear weight dynamics and the th rule of the T-S fuzzy model is of the following form.
where () := ((), ) ∈  −1 is the independent weight vector.() ∈   is the unmeasured state and ( − ()) represents the time-delay state.The time-varying delay () is assumed to satisfy 0 < τ () <  < 1, where  is a known constant.  ,   ,   ,   , and   are known constant real matrices of appropriate dimension.  and   ( = 1, . . ., ,  = 1, . . ., ) are, respectively, the premise variables and the fuzzy sets,  is the number of the If-Then rules, and  is the number of the premise variables.By fuzzy blending, the overall fuzzy model can be defined as follows: where Moreover, we have for any .Hence, ℎ  () satisfies

Fault Detection for T-S Fuzzy Weight Model
To detect the fault existing in the output stochastic distribution, we construct the following fuzzy filter: where x() is the estimated state,  is the gain to be determined later, and () represents the residual signal and is defined as follows: where () can be regarded as the prespecified vector defined on [, ].By defining () = () − x(),  = ∑  =1 ℎ  ()  , the estimated error system can be expressed as where It can be seen that ‖Δ()‖ ≤ ξ, where ξ = ‖ ∫   ()‖ is a known constant.Meanwhile, the residual signal can be expressed as In the fault detection phase, our objective is to find the gain  such that the estimated error system ( 14) is stable if  = 0, which can be formulated in the following theorem.Theorem 1.For the known parameters  > 0,  > 0, suppose that there exist matrices  > 0,  > 0, , such that the following LMIs are solvable, where then the error system ( 14) is stable when  = 0, and, for any  ∈ [−(), +∞), the error variable () satisfies where the gain  can be computed by  =  −1 .
It is noted that Theorem 1 gives a necessary condition for fault detection.Based on the conclusion of Theorem 1, when the fault  = 0, the residual error signal satisfies the following inequality: By defining the threshold  =  0 (‖Γ 1 ‖ + ‖Γ 2 ‖‖‖) + δ, we can conclude that the system has no fault, when the residual error signal is less than or equal to the threshold.If the residual signal is greater than the threshold, the system causes the existence of fault.

Fault Diagnosis for T-S Fuzzy Weight Model
Based on the results regarding fault detection, this part will estimate the size of the fault if the system has fault.We construct the following adaptive fuzzy filter: where F is the estimate value of the fault ,   > 0, ( = 1, 2) are the designed learning rate.Denoting () = () − x(), F() =  − F(), the estimated error system can be expressed as follows: Theorem 2. Suppose that ‖‖ ≤ /2, ‖ F‖ ≤ /2; then one has ‖ F‖ ≤ .For the known parameter  > 0 and matrices   , ( = 1, 2), there exist matrices  > 0,  > 0,  and constants where (29) then the error system (27) is stable in the presence of .The estimation error satisfies for all  ∈ [−, ∞).The gain  of diagnosis filter (26) can be computed by  =  −1 .
In the simulation, it is supposed that δ = ‖Λ 1 ‖ −1 = 5.The fault  is defined as  For the T-S fuzzy weight dynamics, the rule  = 2 and the model parameters are given as follows: Furthermore, we choose the following Gaussian type functions as our member functions: where  = 0.8.By defining  =  = 1, The responses of the residual signal are shown in Figure 1.It can be clearly seen that when  > 5, the fault will occur.Figure 2 is the fault and its estimated value that shows the estimation error can be converged in a small field.Figure 3 shows the 3D mesh plot of the output PDFs and we can find that the fault occurs in front of two peaks.Moreover, when the fault  is redefined as a time-varying function

Conclusions
This paper presents a novel FDD algorithm for non-Gaussian stochastic distribution systems based on T-S fuzzy modeling  and nonlinear fuzzy filter design.A series of LMIs based solution is presented such that the estimation error system is stable and the fault can be detected through a known threshold.Moreover, the adaptive filter based on the T-S fuzzy model is designed to estimate the size of system fault by optimizing the solutions for the concerned LMIs.Simulations are given to demonstrate the efficiency of the proposed approach.

Figure 1 :
Figure 1: The responses of residual signal.

Figure 2 :
Figure 2: Fault and its estimation value.
time-varying fault and its estimated value are shown in Figure4.

Figure 4 :
Figure 4: Time-varying fault and its estimation value.