Journal of Control Science and Engineering

Volume 2017, Article ID 6354208, 9 pages

https://doi.org/10.1155/2017/6354208

## Fault Diagnosis of Nonlinear Uncertain Systems with Triangular Form

Department of Automation, Xiamen University, Xiamen 361005, China

Correspondence should be addressed to Xiafu Peng; nc.anis@ufaixgnep

Received 31 March 2017; Revised 14 June 2017; Accepted 28 June 2017; Published 1 August 2017

Academic Editor: Chunhui Zhao

Copyright © 2017 Qi Ding et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A novel approach to fault diagnosis for a class of nonlinear uncertain systems with triangular form is proposed in this paper. It is based on the extended state observer (ESO) of the active disturbance rejection controller and linearization of dynamic compensation. Firstly, an ESO is designed to jointly estimate the states and the combination of uncertainty, faults, and nonlinear function of nonlinear uncertain systems. It can derive the estimation of nonlinear function via the state estimations and system model. Then, linearization of dynamic compensation is employed to linearize the system by offsetting nonlinear function mandatorily using its estimation. An observer-based residual generator is designed on the basis of the prior linearized model for fault diagnosis. Moreover, threshold treatment technique is adopted to improve the robustness of fault diagnosis. This method is utilizable and simple in construction and parameter tuning. And also we show the construction of ESO and give the corresponding convergence proof succinctly. Finally, a numerical example is presented to illustrate the validity of the proposed fault diagnosis scheme.

#### 1. Introduction

In the past few decades, fault diagnosis technology has been advanced by leaps and bounds, and there has been significant research activity in the design and analysis of fault diagnosis [1–4]. Furthermore, the method based on analytic models is one of the most pioneer and in-depth approaches in fault diagnosis area (see, e.g., [5–14]). Model-based diagnosis can make great use of the system information and has a high diagnosis precision for fault location, fault type judgement, and amplitude estimation of fault. It is mainly applicable for linear systems, whereas in many practical situations, nonlinear properties of the control systems cannot be ignored in dealing with fault diagnosis problem. Therefore, it has practical significance in research fault diagnosis method for nonlinear systems based on analytic models [15].

At present, model-based fault diagnosis for nonlinear systems is generally based upon two kinds of basic theory: One way is to linearize the original nonlinear system near its working point using state feedback linearization technique to obtain the linear model. And the residual is designed based on unknown input decoupled technique to minimize the influence of the modeling error which is treated as unknown input. The other is that nonlinear observer method or nonlinear parameter estimation method is applied to fault diagnosis for nonlinear systems directly, without resort to any linearization. In a word, model-based fault diagnosis can be defined as the determination of faults of a system from the comparison of available system measurements with a prior information represented by the system’s mathematical model through generation of residual quantities and their analysis [16]. A back-propagation neural network fault detection method was employed in [17] to detect the proton exchange membrane fuel cell (PEMFC) current fault type. According to the diagnosis results, the reconfiguration mechanism determined three nonlinear controllers based on feedback linearization approaches which are, respectively, built to adjust the voltage and pressure difference in the case of normal to be selected. The sliding mode observer for the nonlinear fault diagnosis was designed to observe the whole nonlinear system and actuator in [18]. Then the faults are diagnosed according to the residual value changes of two observers. In [15], the system model was decomposed by and the observer with high sensitivity for faults and strong robustness for unknown input disturbances was designed by using the decoupling technology. Reference [19] proposed a new approach based on adaptive sliding mode observer to detect, isolate, and identify faults to deal with the shortcomings that existing model-based fault diagnosis methods of aircraft engines require high precision engine model. Reference [20] established a time-varying failure model of aeroengine control and designed an adaptive observer to estimate the state, actuator, and sensor faults of the system simultaneously. In recent years, various diagnosis methods based on the knowledge have been proposed, including neural network method, rough set-based method, genetic algorithm-based method, expert system method, and wavelet transform method. But all these methods have well-defined control limits [21].

The novel method broke through some limitations in this paper, compared with the above-mentioned methods, as assumption of measurable system variables and amplitude and type limitation about faults. It can be concluded from theoretical analysis and simulation results that the proposed method is easily used to identify the fault and get its approximate function, which differs from traditional methods. Reference [22] investigated fault diagnosis method based on ESO which estimated several faults directly. References [23, 24] also estimated fault directly using ESO to get its approximate function. Without doubt, this method took uncertainty as part of the fault in uncertain systems, which led to large error of estimated fault and actual fault. Reference [25] used the same technique, but here too there was the same issue about misinformation and inaccurate fault estimation. In short, this method has poor robustness.

In this paper, system variables are estimated to derive the estimation of nonlinear function using ESO with high tracking accuracy, high estimation precision, and less sensitivity to model errors and external disturbances. Then, linearization of dynamic compensation is employed to linearize the system via estimated value. And the residual system is constructed based on full-order state observer for the prior linearized model to get approximate function of fault. Moreover, in order to strengthen the robustness of fault diagnosis, this paper adopts threshold value treatment technique. Eventually the complex fault diagnosis problem of nonlinear system is simplified as one of linear system to perform the task of fault diagnosis, and the method is robust.

The rest of this paper is organized as follows. In Section 2, the problem of nonlinear fault diagnosis is formulated. The proposed fault diagnosis scheme is presented in Section 3. In Section 4, the derivation of threshold for robust fault detection is given. Section 5 describes a simulation example to illustrate the effectiveness of the robust fault diagnosis method. Finally, Section 6 presents some concluding remarks.

#### 2. Problem Statement

Consider a class of uncertain nonlinear systems described by where , , and are the system state vector, input vector, and output vector, respectively. is the dynamic process fault of nonlinear system, and is the disturbance or uncertainty of the system. is a real nonlinear vector function and has the form as follows: is the input matrix. is the output matrix. In addition, is the fault coefficient matrix.

There are many practical applications of such systems, such as Vander Pol oscillator system and robot arm system. Research on fault diagnosis for such systems has both theoretical and practical importance. It is shown in [26] that the coordinate change transforms the system into the triangular formwhere represents nonlinear function of the transformed system (4).

##### 2.1. System Variables Estimation

In this paper, we estimate the system variables by designing ESO with high tracking accuracy and high estimation precision.

By defining a new system state in system (4) as , the augmented system can be obtained as

In order to track the system states efficiently, in this work it is proposed to extend uncertainty, faults, and nonlinear function as a new system state and design an ESO as follows:where and are ESO state and input. are observer parameters, and . Based on the experience, the observer parameter can be roughly selected and adjusted in the light of Fibonacci sequence. is sampling period. Also the function can be expressed aswhere and is range length of linear interval which is of a minor amount.

Let , . From (5) and (6), the observer estimation error can be shown as

Theorem 1. *Denoting the Jacobian matrix of system (8) as , let ; that is, . If there exist two symmetric positive definite matrices , , such that , the matrix is negative semidefinite in some neighborhood of the origin. Then system (8) is asymptotically stable.*

*Proof. *For the symmetric positive definite matrices , , let us consider a Lyapunov function candidateAssuming that is singular, then one can find a nonzero vector such that . Sincethe singularity of implies that ; we can get which contradicts the assumed negative definiteness of .

The invertibility and continuity of guarantee that the function can be uniquely inverted. This implies that for . Given the above result, the function is positive definite.

The derivative of can be writtenUsing the fact that , we can get Because is negative semidefinite and is positive definite, is negative definite. If as , system (8) is globally asymptotically stable based on Lyapunov theorem.

Note that as long as we select appropriately parameters , ESO can estimate both states of system (4) and the extended state very well; that is,Also ESO is evolved by the optimal control, and large quantities engineering experiences from practice show that ESO is convergent.

##### 2.2. Linearization Process

The estimation of nonlinear function can be derived and the extended state estimation can be got with system state estimations and system model. On this basis we substitute the estimation of nonlinear function for its value and then select appropriately control variable to achieve linearization of dynamic compensation. In this work, the control variable is given by then nonlinear system (4) is turned into linear integrator series control system (16); that is, where is the sum of not eliminated uncertainty and estimation errors.

It is easy to know that system state estimations contain disturbing information to make uncertainty weakened by linearization of dynamic compensation.

#### 3. Fault Diagnosis

As mentioned above, nonlinear system (4) is turned into linear integrator series control system (16) by linearization of dynamic compensation which simplifies fault diagnosis of nonlinear system into that of linear system, whereafter a full-order fault detection filter is designed for linearized system, achieving convergence of observation error as well as fault diagnosis based on residual error system. Linear integrator series control system (16) can be uniformly expressed in state equations of general forms as, distinctly, with and is an observable pair which is verifiable by .

Based on (17), an observer-based residual generator can be designed aswhere is the gain matrix which needs to be set.

Defining the state estimation and the residual, and , respectively, the dynamics of the residual generator are governed by Denote the sum of fault and error by . Note that is input of the residual generator, so the relation between residual and input can be definitely obtained to perform the task of fault diagnosis. In addition, the design and analysis in this paper are based on the assumptions as follows:(A1)The real time action is bounded.(A2) is a known constant.(A3) remains bounded before and after the occurrence of a fault; that is, .

With this method, the fault diagnosis problem can be transformed into the problem of pole placement. Now the object is to compute the gain to make residual generator asymptotically stable and fulfill the following performance indicators:(A1)Residual generator (20) is asymptotically stable when no fault occurs.(A2)For the required property , transfer function between the sum and residual signal satisfies inequality when a fault occurs.

Theorem 2. *Consider the dynamics of residual generator (20) rewritten by under the assumption given above, suppose the expected eigenvalues of filter (19) are . If there exist two symmetric positive definite matrices and a state feedback matrix verifying Lyapunov equation and , respectively, then residual generator (19) is asymptotically stable, and also .*

*Proof. *For the symmetric positive definite matrices , let us consider a quadratic Lyapunov function candidateIf there are no fault and uncertainty, that is, and , we can get based on (20)Then we get . Since as , the residual generator (19) is globally asymptotically stable based on Lyapunov theorem.

As previously stated, is an observable pair. Thus, integrator series control system (16) is completely observable. And all the eigenvalues of residual generator (19) can be configured arbitrarily based on pole placement theorem; that is, there exists a state feedback matrix such that .

Based on duality principle, the observability of is equivalent to the controllability of . Computing dual coefficient matrices , , pole assignment method is employed to calculate state feedback matrix with and the expected eigenvalues . Therefore, the gain matrix is determined according to the conditions of convergence of the monitor and pole assignment method; that is, .

*Remark 3. *Let us denote the eigenvalues of the matrix by . The equality means that all the eigenvalues of filter (19) can be configured arbitrarily.

From the dynamics of the residual generator (20), we can transfer function between the sum and residual signal as follows:whereThen we haveBy applying the inverse Laplace transformation to (26), we have

The above equation shows the definite relationship between the sum and residual signal . Also Theorem 2 indicates that if there is no fault, filter (19) can accurately track the state and output of system (17). That is to say, the residual signal is asymptotically stable at zero. When a fault occurs, there is a sudden change in the state and output of system (17), but filter (19) cannot track them timely. Thus, the residual signal is large at first which can be used to detect fault. From (27), we can approximately estimate the fault based residual signal to identify fault. Now the comprehensive analysis proves that the proposed scheme is effective for fault diagnosis of system (1).

In addition, linearization of dynamic compensation ensures that the uncertainties would not interfere with fault diagnosis, because the uncertainties are offset by its estimation. References [23, 24] estimated fault directly using ESO to get its approximate function. Without doubt, this method took uncertainty as part of the fault in uncertain systems, which led to large error of estimated fault and actual fault such that the result of fault diagnosis is inaccurate. Compared with [23, 24], the method in this paper is complex to some extent. But the proposed fault diagnosis method based on linear system is mature and applicable. The most important thing is that the proposed method is more robust.

#### 4. Threshold Analysis

Robust technique is employed in this paper to enhance the reliability of fault detection and reduce the false alarm rate. There are two kinds of robust techniques: the unknown input decoupled observer technique and the limited threshold technique. For uncertainties of the system, threshold treatment technique is adopted to realize robust fault evaluation in the algorithm. When the residual signal exceeds the defined threshold, we can identify that there exists fault based on the mechanism. Otherwise, there is no fault.

Theorem 4. *For system (4), ESO (6) is designed to estimate the states, and also the observer-based residual generator (19) is designed to diagnose fault. If there is no fault, the residual signal satisfies the following inequality for the time domain :*

*Proof. *By defining the state transfer matrix , then we can get For , we have At the initial time , state estimations are the same as true values; that is, ; that is, .

Based on norm theory, we obtainBy substituting the above inequality into (20), then we havewhere is the length of the time domain; that is, . AndTherefore the threshold can be set as ; that is, . Then the fault detection can be performed using the following mechanism:

#### 5. Example and Simulation

Consider the following nonlinear system with a fault: From (35), we can know that nonlinear function is

The combination of uncertainty, fault, and nonlinear function is tracked by ESO which can be designed based on system (6). And so are the states , . According to the estimations of states and system structure, the estimation of nonlinear function can be obtained.

In the simulation example, is white noise with an energy of 0.01, permissible estimating error , time-domain length , and the fault detection threshold . The parameters of ESO are adopted as , and step length . And the faulty signal is nominated as having two types:

The simulation considers the effect of the two fault signals separately. When sine function fault occurs, the state estimation errors are given in Figures 1 and 2 showing the combination of uncertainty, fault, and nonlinear function tracking trajectory and it can be seen from this figure that a good tracking performance of the proposed ESO is achieved. Figure 3 shows the capability of the proposed method to track sine function fault and diagnose the fault. When ramp function fault occurs, the state estimation errors given by ESO are illustrated in Figure 4, and the tracking effects of ESO on the combination can be seen in Figure 5. Figure 6 shows the ramp function fault signal and its estimation. From Figures 1 and 4 and Figures 2 and 5, it is seen that all the true values are approached by the estimated ones fleetly on the system states and the combination which is estimated by the ESO. The existence of band-limited white noise influences the accuracy of estimation, but the fault estimations can also track them timely which we can know from Figures 3 and 6.