Abstract

The problem of unmeasured parameters estimation for distributed control systems is studied in this paper. The Takagi–Sugeno fuzzy model which can appropriate any nonlinear systems is employed, and based on the model, an observer-based fuzzy H∞ filter which has robustness against time-delay, external noise, and system uncertainties is designed. The sufficient condition for the existence of the desired filter is derived in terms of linear matrix inequalities (LMIs) solutions. Moreover, the underdetermined estimation problem in which the number of sensors available is typically less than the number of state variables to be estimated is specifically addressed. A systematic method is proposed to produce a model tuning parameter vector of appropriate dimension for the estimation of the filter, and the optimal transformation matrix is selected via iterative solution to minimize the estimated error. Finally, a simulation example for turbofan aeroengine is given to illustrate the effectiveness of the proposed method, and the estimated error is less than 2.5%.

1. Introduction

Along with the continuous development of computer and network technology, distributed control systems (DCSs) have gradually become a new trend and attract much attention [1, 2]. Compared with traditional point-to-point control systems, DCSs have advantages such as high reliability, reduced weight, low cost, and ease of maintenance [3, 4]. The DCSs provide the low-level processing function via intelligent unit and are more conductive to the implementation of complex control algorithms [5, 6]. However, the signal of DCSs is transmitted via bus network, which brings new challenges, such as networked-induced time-delay, packet dropout, and packet disordering.

In practice, many parameters are difficult to obtain due to the limitation of the system itself (for example, the sensors cannot be installed in poor working conditions) or cost constraints. Thus, the estimation of unmeasured system parameters is significantly important, since the estimated parameters can not only be used in the design of controller but also for online monitoring [7, 8]. For the estimation problem of unmeasured system parameters, there are two main methods: (1) high-precision mathematical model [9]; (2) filtering technique and its extension methods [10]. The first method is to establish the mathematical model based on the structure of system to solve the unmeasured parameters. However, the modeling process is often cumbersome and complicated. And for some systems, the structure is extremely complicated and cannot be accurately modeled. Moreover, since this method generally requires iterative calculation, the amount of calculation is large.

For the second method, the estimation of unmeasured parameters is based on the value of measurable parameters using the filtering technique. Because this method uses a recursive algorithm, it is easy to implement on computer and meet the requirements of accuracy and real-time. Among the filters, the H∞ filter stands out because it has robust stability against external noise without priori knowledge of noise and precise mathematical model [1113].

However, in some cases, the number of sensors available is typically less than the number of state variables to be estimated, which is an underdetermined estimation problem and cannot be solved directly. A common approach to address this shortcoming is to estimate a subset of the unknown parameters and assuming that others remain unchanged [1416]. Although this approach enables online filter-based estimation, it will introduce error in the accuracy of overall estimation application. If any of the parameters which are assumed unchanged moves away from their nominal values, the estimates can no longer represent the true parameters. Litt [17] presented a novel approach based on singular value decomposition (SVD) that selects a model tuning parameter vector of low enough dimension. The model tuning parameter vector is constructed as a linear combination of all unmeasured parameters and the transformation matrix is generated by selecting the k most significant terms of singular values, where k is the number of available sensors. Similarly, the SVD is also used in [18] to determine which parts of the systems are observable if the whole system is unobservable. However, the unaccounted terms of singular values may also have great importance on the accuracy of overall estimation application even if they are small under some circumstances.

Furthermore, in [19], the transformation matrix is selected to minimize the theoretical mean squared estimation error at a steady-state open-loop linear design point. But the estimation results are affected by the parameter perturbation, and there are many limitations in application. Moreover, all the approaches above do not consider the impact of time-delay and uncertainties of the system. Thus, it has great significance and practical value to propose a novel approach for the underdetermined estimation problem.

Besides, due to the strong nonlinear fitting characteristics of complex nonlinear systems, the Takagi–Sugeno (T-S) fuzzy model has been widely used in the study of nonlinear systems since it was proposed [2022]. Thus in this paper, an observer-based fuzzy H∞ filter is constructed, and the unmeasured parameters estimation problem for T-S fuzzy distributed control systems with time-delay and parameters perturbation is studied.

The contributions of this paper can be concluded as follows: (i) a model tuning parameter vector of appropriate dimension is produced for the estimation of the filter; (ii) a systematic method is proposed for the selection of optimal transformation matrix to minimize the estimated error via iterative solution; (iii) the parameter perturbation is considered, so that the designed filter has a certain range of margins of disturbance.

The remainder of this paper is organized as follows. Section 2 introduces the modeling method of DCSs and the observer-based fuzzy H∞ filter. The main results are presented in Section 3, including the analysis and synthesis of the filtering error system.Section 4 presents the approach of optimal transformation matrix selection. Numerical simulation results are shown in Section 5 and finally conclusion and a discussion of further application of the method are presented.

2. Problem Formulation

Consider the DCSs with time-delay, which can be described by a class of T-S fuzzy model as follows:

Plant rule i: If is and is , thenwhere denotes the fuzzy set, r denotes the number of IF-THEN rules, and denotes the premise variable. is the vector of state variables; is the vector of control inputs; is the vector of measured outputs; is the vector of unmeasured outputs and is noise signal which belongs to . are matrices with appropriate dimension. are unknown matrices which represent the time-varying uncertainties of the system. is the health parameters of system, which represents the physical characteristics of each component.

Remark 1. The health parameters are considered in this paper, because the performance degradation of each component in the system is inevitable during the working process. And the system’s performance is affected by the level of degradation, which is generally described in terms of unmeasured health parameters such as efficiencies and capacities related to each major module in most cases. If the health parameters move away from their nominal values, the shift in other performance variables will be induced. They may be treated as a set of biases and can be augmented to the system states.
It is assumed that both sensors and actuators are time-driven. The data has timestamp and is transmitted in a single-packet, and incorrect order of the data packet does not exist. Since time-delay is related to the bus load at a certain moment, it is assumed that the time-delay has an upper bound, which is (under the above assumption, the upper bound of time-delay can be estimated; see [23] in detail). As a consequence, can be modeled as a finite state Markov stochastic process on a finite set . The transition probability from at time to at time iswhere and . is the transition probability rates from at time t to at time , and there is .
The state feedback controller is applied in this paper, which is . Then, by fuzzy blending, the global fuzzy augmented model can be obtained as follows:whereFor the underdetermined estimation problem, a low-dimensional model tuning parameter vector is produced to represent the information of high-dimensional health parameter vector . The model tuning parameter vector is constructed as a linear combination of all health parameters, given bywhere , and is a transformation matrix, which is applied to construct the tuning parameter vector. The selection of optimal is obtained in Section 4.
Thus, the estimation of the health parameters can be obtained aswhere is the pseudoinverse of .
The final augmented dynamic fuzzy model can be rewritten aswhere

Remark 2. If the model tuning parameter vector is of appropriate dimension, the augmented system (7) is observable and the observer-based filter can be used to estimate the unknown health and performance parameters. And the key of the estimation is to find the optimal transformation matrix so that the low dimension tuning vector can represent as much of the information as possible. Moreover, different from the method in [1719], the solution of the optimal transformation matrix is transformed into an optimization problem.
It is assumed that the uncertainties of the system can be described in the following form:where are known matrices and is a time-varying unknown matrix which satisfiesThen the following observer-based fuzzy filter is constructed:where is the filter state, and and are the filter outputs. is the filter parameter to be determined.
The estimation error of state isDefine and , and the filtering error system is given bywhereDefine the following parameters:The goal is to design a filter in form of (11) so that the filtering error system can meet the following requirements simultaneously:(1)The filtering error system is asymptotically stable when (2)Under the zero initial condition, the filtering error system satisfiesfor any nonzero .

3. Main Results

In this section, the sufficient condition for the existence of the desired filter is derived in terms of LMIs solutions. Before proceeding with the study, the following Lemma is needed.

Lemma 1. (see [24]). are real matrices with appropriate dimensions, and is a time-varying unknown matrix which satisfies . Then for a scalar , the following inequalityalways holds.

Theorem 1. If there exist positive matrices , and matrix such that the following inequality holds:wherethen the filtering error system is asymptotically stable and the prescribed H∞ performance can be guaranteed.

Proof. Select the Lyapunov function aswhereWhen , the derivation of isIt is obvious that there isThus, can be further rewritten aswhereBy the Schur complement, if the inequality (24) holds, it can be obtained that there is . Thus, there is and the filtering error system is asymptotically stable.
Secondly, define a new functionwhere is a scalar which satisfies .
Thus, under the zero initial condition, there iswhereAccording to the Schur complement, it follows thatThus, if the inequality holds, there is and equation (16) holds for any nonzero . The proof is completed.

Theorem 2. For a given scalar , the filtering problem is solvable under the conditions above with a H∞ performance level , if there exist positive matrices , matrices and scalars such that the following inequalities hold:whereMoreover, the filter parameter can be solved by

Proof. LetThen can be rewritten aswhereThen, for further analysis, based on equation (9) and Lemma 1, equation (35) can be decomposed into

Moreover, define

By applying the Schur complement, equation (37) is equivalent to

Thus, the inequalities (30) and (31) hold. The proof is completed.

4. Optimal Transformation Matrix Selection

It can readily be obtained that the estimation precision is directly affected by the estimation of the unknown states (health parameters). Thus, the key of the parameters estimation is to find an optimal transformation matrix , which can ensure that the low-dimensional tuning vector can represent as much of the information of the health parameters as possible.

In this paper, the optimization objective is to find an optimal transformation matrix , which can ensure that the estimation error of the measured outputs is minimum.

Thus, the optimization problem can be described as

And the method based on iterative solution is used in this paper, and the process of optimal transformation matrix selection is shown in Figure 1:

(1)To begin with, the initial value of is given randomly. And during the estimation process, the optimal of the last moment is used as the initial value to reduce the quantity of calculation. In order to avoid the calculation converging to a poorly scaled result, the Frobenius norm of must satisfy .(2)With reference to equation (7), construct the reduced-order state-space model.(3)Solve the parameters of the H∞ filter.(4)Calculate estimated error using equation (40).(5)Determine whether the estimated error achieves convergence within a tolerance (user-specified):(1)If converged, skip step 6 and proceed directly to step 7.(2)If not converged, proceed to step 6.(3)If the number of iterations exceeds 100, skip step 6 and proceed directly to step 7 to meet the real-time requirements.(6)Use the MATLAB lsqnonlin function to update , and the new value still needs to satisfy .(7)Return the optimal value of , and ends.

5. Simulation Example

Before proceeding with the simulation, the high-precision aero-thermodynamic mathematical model of aeroengine in full envelope is the basis of the simulation. It can not only build the small deviation dynamic state space models, but can also replace the real engine in simulation. The detailed modeling process is carried out by reference to [25, 26]. This paper takes a type of double-rotor turbofan engine as research object and uses modeling method based on component characteristics. The afterburner is not considered in this engine model. For greater adaptability, faster calculation speed, and stronger convergence, the self-tuning Broyden quasi-Newton method [27] is used to solve equilibrium equations in the component-level mathematical model. And the small deviation dynamic model is solved by fitting method [2830].

Assuming that the aeroengine works in the maximum state and then based on the methods of flight envelop division in [3133], a T-S fuzzy model is constructed and the initial membership function is shown in Figure 2. In this paper, the T-S fuzzy model has two state variables, four health parameters, two control inputs, five measured outputs, and three unmeasured outputs, all shown in Table 1.

Then the four T-S fuzzy rules with the selection of and can be obtained as follows:Rule 1. If is about 0 km and is about 0, thenRule 2. If is about 4.8 km and is about 1.24, thenRule 3. If is about 12.8 km and is about 0.7, thenRule 4. If is about 18.8 km and is about 1.72, then

Moreover, the uncertainties of the system are

The system sampling period is set as , and the external disturbance is

The state space of time-delay is , and the state transition matrix of is

The working condition is chosen as H = 5 km and Ma = 0.3. The given value of is .

Two types of performance degradation are selected to simulate, which are slow varying and sudden change of health parameters. For the sudden change type (Type I), there is a step change for health parameters, and this type is used to represent the sudden damage, such as foreign objects damage. And the slow varying type (Type II) is defined as the linear variation of health parameters to represent the gradual performance degradation during using process. These two types can cover almost all the health parameters changes in practice, and they are simulated respectively to illustrate the effectiveness of the method proposed in this paper.

The results are shown in Figures 36. It can be seen that the method proposed in this paper can effectively simulate the health and performance parameters of the aeroengine. The estimated values can accurately track the changes as the health parameters shift.

By comparing the results of the two types of health parameters change, it can be concluded that the estimated values for slow varying type are more accurate. For the type of sudden change, it takes a while for the estimated values to stabilize and the overshoot is slightly larger.

Table 2 shows the average estimated error of the health and performance parameters. It can be seen that the estimated error is within 2.5% for all types. At the same time, the results are more accurate except for the beginning of simulation.

Furthermore, the results show that the method proposed in this paper can accurately estimate most system parameters of interest. On this basis, the method can also be used in the investigations of direct control, fault diagnosis, health management, and online monitoring. And for aeroengine, there are already some researches which indicate that the changes of health parameters can characterize specific faults, allowing for online monitoring and alerting [3436].

6. Conclusion

In this paper, a method is proposed for the underdetermined estimation problem of the distributed control systems. First, the T-S fuzzy model for DCSs is constructed and a model tuning parameter vector of appropriate dimension is produced. Then an observer-based fuzzy filter is designed and the sufficient condition for the existence of the designed filter is derived in terms of LMIs solutions. Besides, the method based on iterative solution is used to select the optimal transformation matrix to minimize the estimated error. Finally, the results of the simulation show that the proposed method can effectively estimate unmeasured parameters of the aeroengine and the estimated error is less than 2.5%.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (51476187, 51606219, and 51506221).