Mathematical Problems in Engineering

Volume 2019, Article ID 5839576, 8 pages

https://doi.org/10.1155/2019/5839576

## Sensor Fault Diagnosis and Fault-Tolerant Control for Non-Gaussian Stochastic Distribution Systems

School of Electrical Engineering, Zhengzhou University, Zhengzhou 450001, China

Correspondence should be addressed to Lina Yao; moc.621@qxnl_ellehcim

Received 3 December 2018; Accepted 23 January 2019; Published 11 February 2019

Academic Editor: Xiangyu Meng

Copyright © 2019 Hao Wang and Lina Yao. 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 sensor fault diagnosis method based on learning observer is proposed for non-Gaussian stochastic distribution control (SDC) systems. First, the system is modeled, and the linear B-spline is used to approximate the probability density function (PDF) of the system output. Then a new state variable is introduced, and the original system is transformed to an augmentation system. The observer is designed for the augmented system to estimate the fault. The observer gain and unknown parameters can be obtained by solving the linear matrix inequality (LMI). The fault influence can be compensated by the fault estimation information to achieve fault-tolerant control. Sliding mode control is used to make the PDF of the system output to track the desired distribution. MATLAB is used to verify the fault diagnosis and fault-tolerant control results.

#### 1. Introduction

With the rapid development of modern science and technology, control systems have become more complex, large-scale, and intelligent, which puts more demands on the control of engineering control systems [1, 2]. The safety and reliability of the system must be ensured during the operation of the system. Otherwise, huge personal and property losses will be caused. Luckily, fault diagnosis and fault-tolerant control play an important role in detecting and avoiding such accidents. It is essential to carry out the research of controlling stochastic distribution systems in the field of modern control theory, which is widely used in controlling of the ore particle distribution in the grinding process, controlling of pulp uniformity and controlling of particle uniformity in the papermaking process, and the polymer in the process of the chemical reaction as well as the flame distribution control during the combustion process of the boiler.

In reality, there are many production sites that are out of reach and are accompanied by high temperature, high pressure, and toxic environments. However, the fault is inevitable. If it fails to be dealt with in time, it will waste precious time and may lead to extremely serious consequences. Therefore, it is of great significance to carry out research on fault diagnosis and fault-tolerant control of stochastic distribution systems to improve its reliability and safety and avoid loss of personnel and property. Stochastic distribution control has been regarded an important topic in the control field in recent decades. In most control projects, practical systems are subjected to stochastic input. These inputs may be derived from noise, stochastic disturbances, or stochastic parameter changes. Professor Hong Wang put forward to the idea of stochastic distribution control [3]. Different from the model in the existing control systems, the overall shape of the probability density function of the output of the stochastic system is considered. The goal of the controller design is to select a good rigid control input so that the probability density function shape of the system output can track the given distribution.

At present, there are many research results for fault diagnosis and fault-tolerant control of actuator fault in non-Gaussian stochastic distribution systems. For fault diagnosis, fault diagnosis observers or filters based methods are usually used. In literatures [4–6], the adaptive fault diagnosis observer is used to diagnose the actuators fault, and the fault amplitude is accurately estimated. In literatures [7, 8], fault reconstruction based on learning observer is recorded. In literatures [9–11], sensor fault diagnosis and output feedback control are described. In literature [12], the filter based method has been shown to be effective method for fault diagnosis. For fault-tolerant control of non-Gaussian stochastic distribution systems, two situations are considered: (1) the desired PDF is known; (2) the desired PDF is not known in advance. An adaptive PI tracking fault-tolerant controller is designed to track the desired PDF in the literature [13]. In literatures [14, 15], the minimum entropy fault-tolerant control algorithm is constructed based on the unknown PDF, and the performance index function is designed to find the control input to make the performance index function be minimized. In literature [16], collaborative system fault diagnosis and model prediction fault-tolerant control for the stochastic distribution system are described. In literature [17], actuator fault diagnosis and PID tracking fault-tolerant control for n subsystems collaboration are introduced.

The research of sensor fault diagnosis and fault-tolerant control of non-Gaussian stochastic distribution system is rarely documented; however, the sensor fault is inevitable. Thus, it is very meaningful for the work to be carried out in this paper. In this paper, the learning observer is used to diagnose the fault of the non-Gaussian stochastic distribution system, and the fault is compensated using the fault estimation information. The sliding mode control algorithm is used to make the system output PDF to track the expected PDF.

#### 2. Model Description

The system output probability density function (PDF) is approximated by linear B-spline function. are basis functions defined in advance on the interval . are corresponding weights associated with the number of basis functions. can be expressed as

Since the integral of on is equal to 1, the following equation holds: where . Therefore, only weights are independent of each other, and the linear B-spline model is specifically given as follows: where , , . is the PDF approximation error that can be ignored. The non-Gaussian stochastic distribution system model can be expressed as where is the state vector, is the control input vector, , is the fault vector and the weight vector, respectively, is a full rank matrix. is observable; are known matrices with appropriate dimension. A new state variable is introduced for fault diagnosis [18]. where is a Hurwitz matrix, . Combined with (4) and (5), the augmented system model can be expressed as follows:where

#### 3. Fault Diagnosis

In order to estimate the size of the fault, the fault diagnosis observer is designed as follows: where is the estimation of state, is the estimation of weight vector, is a state variable, is the estimation of , and is the residual can be expressed as , is observable, and it is easy to know that is observable. The parameter is defined as the learning interval, which can be taken as an integer multiple of the sampling period or sampling period [19]. represents the state observation error. , , , and are gain matrices with appropriate dimension that need to be determined.

Lemma 1. *For two matrices and with approximate dimension, the following inequality holds [20]: *

*Assumption 2. *, where is a given positive constant. The observation error dynamic system obtained by (6), (7) and (8) is formulated as follows:

Theorem 3. *It is supposed that Assumption 2 holds. If there exist positive-definite symmetric matrices , , , and matrices , , , , the following inequalities and equation hold: where is a given positive constant. The observer gain can be obtained from . Equation (12) can be converted into the following LMI: where is a small positive number. Lemma 1 can ensure that the following inequalities hold: *

In order to prove the stability of the system (11), the following Lyapunov function is selected as The derivative can be obtained as follows: From (14) and (15), the following inequality can be obtained: For (20), when inequalities (12), (13), and (14) are satisfied, the following inequality can be further obtained: where . The proof is completed.

After diagnosing the system sensor fault, the original system is rebuilt to observe the original system PDF after the fault occurs and the subsequent fault-tolerant control is prepared. The original system state, weights, and output PDF and their observations can be obtained by the following coordinate transformation. Inequality (16) can be proved as follows: It can be formulated that It can be further obtained that

#### 4. Fault-Tolerant Control

The basic idea of fault-tolerant control (FTC) is to compensate the influence of fault on the system performance. Sensor fault may exist in many forms, such as constant, time-varying, or even unbounded. When the system sensor fault occurs, it usually does not work properly. A simple fault compensation scheme is designed in this paper. is the weight vector of the system output, which needs to be compensated in order to ensure the performance of the system after the fault occurs. The fault estimation can be obtained by the observer. The compensated weight can be expressed as follows:

Using a compensated sensor output , control algorithms are given to ensure that whether fault occurs or not, sensor fault compensation can be performed using fault-tolerant operation. To make sure that the probability density function (PDF) of the system can still track a given probability density function after a fault occurs, a fault-tolerant controller needs to be designed. A new weight error vector is defined as , where is the weight of the expected output. The following weight error dynamic system can be obtained.where and .

With the sliding mode control method, two necessary conditions should be satisfied: the accessibility of the system state and the asymptotic stability of the sliding mode dynamic process. The sliding mode control law is designed so that the state trajectory at any moment can reach the sliding surface during a limited time. The switching function is designed as follows:where the matrix is selected such that is a Hurwitz matrix. The matrix is chosen to make be a nonsingular matrix. After the state trajectory of the weight error dynamic system reaches the sliding surface, and should be satisfied simultaneously. The equivalent control law is shown as follows:Equation (29) is substituted into (27), and (27) can be further obtained as

Theorem 4. *If both positive-definite symmetric matrices and exist, the following inequality holds: Then the sliding mode dynamic system is stable. The following Lyapunov function is selected as The first-order derivative of is formulated as where and . Using (33), the following inequality can be obtained as when , , and the sliding mode dynamic system is stable. In order to ensure that the state trajectory starting at any position can reach the sliding surface during a limited time, the following sliding mode control law is designed: where . The sliding mode control law can ensure that the state trajectory reaches the sliding surface during a limited time, .**The following Lyapunov function is chosen as where is selected as , . The state trajectory reaches the sliding surface during a limited time. The available fault-tolerant control law is shown as follows: *

#### 5. A Simulation Example

In order to verify the effectiveness of the algorithm, the proposed method is applied to the process of molecular weight distribution (MWD) dynamic modeling and control, and a continuous stirring reactor (CSTR) is considered as an example. The closed-loop control diagram of the polymerization process is shown in Figure 1. The specific mathematical model is shown as follows: where is the average residence time of reactants (), is the initial concentration of initiator (); is the initiator concentration ( ); is the initial concentration of monomer ( ); is the monomer concentration ( ); are the reaction rate constants; are constants associated with the control input; are the free radical. The above specific physical meaning can be referred to the literature [21].