Research Article | Open Access

Zhanshan Wang, Huanxin Guan, Chengde Zheng, "Fault Diagnosis Observer Design for Discrete-Time Delayed Complex Interconnected Networks with Linear Coupling", *Mathematical Problems in Engineering*, vol. 2012, Article ID 860489, 22 pages, 2012. https://doi.org/10.1155/2012/860489

# Fault Diagnosis Observer Design for Discrete-Time Delayed Complex Interconnected Networks with Linear Coupling

**Academic Editor:**Bin Jiang

#### Abstract

Fault diagnosis for a class of discrete-time delayed complex interconnected networks with linear coupling in the case of actuator fault is studied. For the case of unavailability of network state, a state observer is first designed. Then a fault diagnosis observer is designed to detect the actuator fault on the basis of online adaptive approximator, which can approximate the unmodeled dynamics of the complex networks. Lastly, by choosing a suitable threshold, the actuator fault can be detected. A numerical simulation is used to show the effectiveness of the proposed method.

#### 1. Introduction

Fault diagnosis for dynamic systems has received more and more attention in the last decades due to the increasing demand for higher performance, higher safety and reliability standards. The main task of the fault diagnosis can be roughly described as the early detection and diagnosis of faulty elements in a system, as well as the time of the detection. Many fault diagnosis approaches have been reported in the existing literature [1–3]. All these proposed fault diagnosis approaches have been applied to practical processes successfully [4–11], to name a few. Among these techniques, one important type is the observer-based approach, where fault diagnosis observers are used to generate residual signals to detect the fault.

Model-based schemes have emerged as prominent approaches to fault diagnosis of continuous and discrete-time systems [1, 2, 10–15]. This approach is built on a mathematical model of the process that must be monitored, so that a number of residuals can be computed by taking the difference between the estimated value of the system output variables and their measured value. The residuals are then compared to suitable thresholds by detection and isolation logics in order to provide fault decisions regarding the health of the system [16–18]. Model-based approaches are well suited to monitor the centralized and small systems, but suffer from scalability and robustness issues when distributed and large-scale systems are concerned, for which there have been significant research activities recently [19]. In [20], a distributed fault detection scheme was introduced for a class of continuous-time systems based on an overlapping decomposition technique, in which the subsystems were monitored by a network of interconnected local fault diagnosis mechanism.

Note that all above concerned systems are mainly focused on isolated node systems, which are the counterparts of the complex interconnected network systems [21]. Over the past few years, complex networks have been gaining increasing research attention because of their potential applications in many real-world systems from a variety of fields such as biology, social systems, linguistic networks, and technological systems [22–27]. In particular, rapidly growing research results have been reported in the literature that have focused on the structural properties between the coupled nodes in the complex networks such as stability and robustness. As one of the mostly investigated dynamical behaviors, synchronization in complex networks with or without time delays has drawn significant research interest in recent years [22, 23, 26, 28, 29]. It is well known that the faults may happen in the complex network systems, which is composed of the same or different node system. Therefore, the fault diagnosis problem for the complex interconnected networks plays an important role in the fault diagnosis theory.

In the recent papers [30–32], the authors first studied the fault diagnosis problems in the complex interconnected networks, and the adaptive control methods were adopted to realize the fault-tolerant synchronization in the case of network deterioration. The main object in [30–32] is to adjust the coupling strength using the adaptive method. Obviously, this kind of deteriorated networks can have robust synchronization stability in the case of uncertainties. In [33], a state feedback-based decentralized control scheme was proposed for a class of networked large-scale systems with output sensor failures. However, the designed controller required all the state combination of the whole networks, which made the method in [33] restrictive in practice. In contrast to the fault diagnosis problem in the conventional fault diagnosis theory [16, 17], there are many problems to be resolved in the complex networks, for example, how to detect the fault in the isolated node system such as actuator fault and sensor fault. Such kind of problems is not studied deeply at present.

Motivated by the above discussions, this paper will study the fault diagnosis problem for a class of discrete-time complex interconnected networks in the presence of actuator fault. Considering the unavailability of the system states, a state observer is first constructed to estimate the states via the linear matrix inequality (LMI) method. When the system states are available, a fault diagnosis observer is designed based on an adaptive technique, in which the uncertainty is approximated by an online approximator. The adopted method is a generalization of the isolated node system to the complex system, in which the difficulty lies in dealing with the coupling of the interconnected subsystems. By constructing a set of suitable observers via the spirit of distributed control concept, we have discussed a special class of actuator faults in the complex system. A numerical example is used to show the effectiveness of the obtained result.

#### 2. Problem Description and Preliminaries

In this paper we will discuss the following complex interconnected networks with identical nodes: where is a positive integer, , and is the dimensional number of node system state . is the linearly dominated matrix with appropriate dimension, and is the known nonlinear function. is the control input matrix; is the control input. Positive integer is the discrete time-varying delay satisfying , where , are known integers and is a positive diagonal matrix. is the coupled configuration matrix of the networks with but not all zeros, , , denotes the disturbances and uncertainties, . is the system output, and is the system output matrix with appropriate dimension.

*Assumption 2.1. *The nonlinear function is bounded and continuous, which satisfies , where is a positive constant,
for any , , , and and are constant matrices.

Note that condition (2.3) can be changed into the following form:
or
Combining inequalities (2.4) and (2.5), we have

It should be noticed that the above model (2.1)-(2.2) may represent the discrete-time equivalent of a continuous model obtained via the well-known Euler method. The class of failure considered in this work is that of actuator fault. This class of failure can be represented as an unknown additive disturbance on the nominal input to the system . Hence, an actuator fault occurring at the th time step results in a faulty input given by
where represents the time profile of the unknown fault. Therefore, the nominal dynamics (2.1) in the presence of faults becomes
where the fault vector is given by .

The uncertain term is assumed to depend on the nominal input and on the parameter vector . If is linear in the parameter vector, it can be expressed as , where the matrix is assumed to be known, while is usually unknown (or partially known). If is not linear in the parameter or its structure is not exactly known, an approximation scheme can be obtained by resorting to the so-called online interpolators [34–38] (e.g., neural networks, fuzzy logic, and splines). By choosing a linear-in-the-parameters interpolator structure, the uncertain term can be expressed as follows:
where represents the interpolation error.

The following assumption is required.

*Assumption 2.2. * is a column full rank matrix, .

*Assumption 2.3. *The norm of the matrix is uniformly bounded by a constant , and the norm of the interpolation error is uniformly bounded by a constant .

Usually, the network output is a linear combination of the information about the network nodes. We wish to design an observer to estimate the network states through the available network output. In order to estimate the states, the following state observer is designed, :

Therefore, in view of (2.8) and (2.10), the state estimation error is given by
where , , , , , .

Before we give the main results, we need the following preliminaries.

Lemma 2.4 (see [26, 28]). *Let , , , , , . If and each row sum of is zero, then
**
where symbol denotes the Kronecker product. *

Lemma 2.5 (Schur complement, see [22, 27]). *Given constant matrices , and , where and is a positive definite symmetric matrix, then
**
if and only if
*

Lemma 2.6 (see [22, 29]). *Let and , , , be matrices with appropriate dimensions. The following statements about Kronecker product are true:*(1)*;
*(2)*;*(3)*;*(4)*;*(5)*.*

#### 3. State Observer Design

Now we state our main results in this section.

Suppose that the online approximator has met the precision of modelling of the uncertainty, then approaches to zero. For the case of no fault and no approximation error in (2.11), that is, , we will study the stability of the following state observer, : or in a compact form where , , .

Theorem 3.1. *Suppose that Assumption 2.1 holds. The observer error system (3.1) or (3.2) is globally asymptotically stable if, for a given observer gain matrix , there exist positive constant , symmetric positive definite matrices and with appropriate dimensions, such that the following condition holds, :
**
where is an identity matrix with appropriate dimension,
** is the elements of matrix .*

*Proof. *Let us consider the Lyapunov functional , where
where

Calculating the difference of along the system (3.2), we have
where

Note that the following equalities hold according to Lemma 2.6:

Substituting (3.8) and (3.9) into (3.7), it yields

According to Lemma 2.4 and inequality (2.6), (3.10) can be changed into the following form:
where is a free parameter to be adjusted, , denotes the maximum eigenvalue of a matrix , , , , and . Considering the negativeness of matrix in (3.3), we have . Letting , we have , and then
Let be a positive integer; then from (3.12) one has
which implies that
By letting , we can deduce that the series is convergent for , and therefore, , namely, , which completes the proof of Theorem 3.1.

In order to solve the observer gain matrix , according to Lemma 2.5, we have the following result.

Theorem 3.2. *Suppose that Assumption 2.1 holds. The observer error system (3.1) or (3.2) is globally asymptotically stable if there exist positive constant , symmetric positive definite matrices and with appropriate dimensions, and appropriately dimensioned matrix , such that the following condition holds, :
**
where
**
and the observer gain matrix .*

*Remark 3.3. *Note that we are mainly concerned with the complex interconnected networks with identical node system in this paper. This kind of complex system can represent many practical systems such as secure communication, chaos generators design, and harmonic oscillation generation [23–26] (and references cited therein). Synchronization problems for such kind of complex networks have been paid much attention in recent years. However, for the fault diagnosis problem, there are few results to be reported. As a complex dynamical system, it is inevitable for a fault to occur in the normal operation. Therefore, combining the knowledge of fault diagnosis theory for isolated node system and the complex networks theory, the authors try to extend the diagnosis method in node system to the complex system. Correspondingly, we only consider a special kind of faults, that is, actuator fault, which has the same form in each node system. The key point of the proposed method lies in the choice of suitable state observer and diagnosis observer. For other kinds of faults, such as sensor faults, parameter faults, and hybrid faults, it still needs to be further studied in the future.

*Remark 3.4. *The LMI condition in Theorem 3.1 requires the information of the interconnected couplings , , the node system parameters , besides a given common observer gain matrix . For a prescribed matrix , if the LMI condition (3.3) holds, the designed observer exists and an estimated state can be achieved. In contrast, Theorem 3.2 aims to present an observer gain matrix directly, which has the same effect on the state estimation.

#### 4. Adaptive Observer-Based Fault Diagnosis

Since all the states of the concerned system can be measured, in this section, we will consider the following fault diagnosis observer, which is used to detect the fault, :

In view of (2.8) and (4.1), the diagnosis error dynamics is given by where represents the uncertainties estimation error, , and is chosen to make the eigenvalues of in the unit circle, .

The residuals vector can be chosen as which can be rewritten as where .

Note that the residual vector is affected by the fault vector and the estimation error of the uncertain term. If an accurate estimation of uncertain term is achieved, the fault signature on the residual (i.e., its effect on the residuals) becomes more evident. If a parametric model of the uncertainties is available, an adaptive estimation algorithm of the unknown parameters can be set up. It is worth remarking that such a paradigm has been keenly exploited for adaptive fault identification. However, in this work the same concept is exploited in order to adaptively compensate for the uncertainties, so as to obtain small values of the residuals in the absence of faults. In this case, the uncertain term can be indirectly evaluated through the estimation of . Thus, an adaptive update law for the parameters estimate can be chosen as where is a positive definite symmetric matrix.

In view of (2.9), the uncertainties estimation error can be written as where is the parameters estimation error.

Therefore, we have the following diagnosis error system:

In order to implement the fault detection, we first consider the stability problem of diagnosis error system (4.8) in the absence of fault and interpolation error, that is,

Theorem 4.1. *Under Assumption 2.3, the diagnosis error system (4.9) is globally uniformly stable, and the error converges asymptotically to zero, .*

*Proof. *Consider the following functional:
where is a positive constant and is the solution to the Ricatti equation for a given symmetric and positive definite matrix .

The difference of can be calculated as follows, with the consideration of (4.6):
Considering (4.7) with , from (4.11) we can deduce
where . holds if satisfies the inequality

Since is a decreasing and nonnegative function, it converges to a constant value as , hence . This implies that both and remain bounded for all , and approaches to zero.

Once the residual vector is computed at each step, a fault is declared if each component of exceeds a suitably selected threshold. A priori selection of each threshold should be based on the expressions of the residuals vector. Namely, proper setting of the thresholds requires an accurate knowledge of the uncertainties influence on the residuals. However, this approach often leads to extremely conservative results. Therefore, an empirical approach may be pursued to set the residuals thresholds in alternative to (or in combination with) the approach based on a priori knowledge of the uncertainties.

#### 5. Numerical Simulation

In this section, we will use an example to show the effectiveness of the proposed result.

*Example 5.1. *Let us consider a complex system composed by systems (2.1) and (2.2) with and . The system parameters in (2.1) and (2.2) are as follows, , , , if , and if . , . The nonlinear vector-value function is given by , . Obviously, , .

By using the Matlab LMI Toolbox, using Theorem 3.2 to solve inequalities (3.15), the feasible solutions are as follows: , , , and . Correspondingly, the observer gain matrix .

When , and the initial states are randomly chosen, using the state observer (2.10), the state curves are depicted in Figures 1–3, respectively. Obviously, the estimated states can track the system states accurately.

When , , and the gain matrix in (4.1) is chosen as (in which the eigenvalues of are and , resp.), the states of the diagnosis observe (4.1) and the system states (2.1) are depicted in Figures 4–6, respectively.

When the following faults occur in the first node system:
the state errors of , diagnosis observer states , and the system state are depicted in Figures 6–8, respectively. The residual curves are depicted in Figure 9. If we choose suitably the fault threshold, we can accurately detect the fault.

Note that we assume that the fault occurs in the first node system. Since three node systems constitute a complex interconnected network, the fault in the first node system will affect the other two systems. The state response curves of the other systems and are depicted in Figures 10–11, respectively. Obviously, the states in the other two systems will be influenced by the first faulty system . In this case, if we judge whether the system is healthy or not, it is not sufficient to detect the fault by the other system. However, if we observe the residual curves of the other two systems, see Figure 12, we may judge the fault location by choosing a suitable logic rule, which is out of the coverage in the paper.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

#### 6. Conclusions

In this paper, a state observer is constructed for a class of discrete-time delayed complex interconnected networks with linear coupling based on LMI method. Then on the availability of the system states, an adaptive fault diagnosis observer is designed to realize the actuator fault detection. The uncertainty of the system is modelled by an online approximation technique, and the stability of the residual system and approximation system are analyzed. Finally, a numerical example is used to validate the proposed result. Note that the proposed state observer and the fault diagnosis observer can only implement the diagnosis task. How to realize the fault identification will be further studied.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grants 61074073 and 61034005, Program for New Century Excellent Talents in University of China (NCET-10-0306), and the Fundamental Research Funds for the Central Universities under Grants N110504001 and N100104102.

#### References

- J. Chen and R. Patton,
*Robust Model-Based Fault Diagnosis for Dynamic Systems*, Kluwer, Boston, Mass, USA, 1999. - S. Narasimhan and G. Biswas, “Model-based diagnosis of hybrid systems,”
*IEEE Transactions on Systems, Man, and Cybernetics A*, vol. 37, no. 3, pp. 348–361, 2007. View at: Publisher Site | Google Scholar - J. Stoustrup and H. H. Niemann, “Fault estimation—a standard problem approach,”
*International Journal of Robust and Nonlinear Control*, vol. 12, no. 8, pp. 649–673, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH - F. Basile, P. Chiacchio, and G. De Tommasi, “An efficient approach for online diagnosis of discrete event systems,”
*IEEE Transactions on Automatic Control*, vol. 54, no. 4, pp. 748–759, 2009. View at: Publisher Site | Google Scholar - N. Meskin and K. Khorasani, “Robust fault detection and isolation of time-delay systems using a geometric approach,”
*Automatica*, vol. 45, no. 6, pp. 1567–1573, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - Q. Yang and Y. Sun, “Automated fault accommodation for discrete-time systems using online approximators,” in
*Proceedings of the 30th Chinese Control Conference*, Yantai, China, July 2011. View at: Google Scholar - R. M. G. Ferrari, T. Parisini, and M. M. Polycarpou, “Distributed fault diagnosis of large-scale discrete-time nonlinear systems: new results on the isolation problem,” in
*Proceedings of the 49th IEEE Conference on Decision and Control (CDC '10)*, pp. 1619–1626, Hilton Atlanta Hotel, Atlanta, Ga, USA, December 2010. View at: Publisher Site | Google Scholar - H. Yang, B. Jiang, and V. Cocquempot, “Supervisory fault-tolerant regulation for nonlinear systems,”
*Nonlinear Analysis*, vol. 12, no. 2, pp. 789–798, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - I. Samy, I. Postlethwaite, and D.-W. Gu, “Survey and application of sensor fault detection and isolation schemes,”
*Control Engineering Practice*, vol. 19, no. 7, pp. 658–674, 2011. View at: Publisher Site | Google Scholar - R. Patton, P. Frank, and R. Clark,
*Fault Diagnosis in Dynamic Systems: Theory and Application*, Prentice Hall, 1989. - M. Blanke, M. Kinnaert, J. Lunze, and M. Staroswiecki,
*Diagnosis and Fault Tolerant Control*, Berlin, Germany, 2006. - R. Isermann, “Process fault detection based on modeling and estimation methods—a survey,”
*Automatica*, vol. 20, no. 4, pp. 387–404, 1984. View at: Google Scholar - J. J. Gertler, “Survey of model-based failure detection and isolation in complex plants,”
*IEEE Control Systems Magazine*, vol. 8, no. 6, pp. 3–11, 1988. View at: Google Scholar - P. M. Frank, “Fault diagnosis in dynamic systems using analytical and knowledge-based redundancy. A survey and some new results,”
*Automatica*, vol. 26, no. 3, pp. 459–474, 1990. View at: Publisher Site | Google Scholar - V. Venkatasubramanian, R. Rengaswamy, K. Yin, and S. N. Kavuri, “A review of process fault detection and diagnosis part I: quantitative model-based methods,”
*Computers and Chemical Engineering*, vol. 27, no. 3, pp. 293–311, 2003. View at: Publisher Site | Google Scholar - B. Jiang, M. Staroswiecki, and V. Cocquempot, “Fault accommodation for nonlinear dynamic systems,”
*IEEE Transactions on Automatic Control*, vol. 51, no. 9, pp. 1578–1583, 2006. View at: Publisher Site | Google Scholar - B. Jiang and F. N. Chowdhury, “Fault estimation and accommodation for linear MIMO discrete-time systems,”
*IEEE Transactions on Control Systems Technology*, vol. 13, no. 3, pp. 493–499, 2005. View at: Publisher Site | Google Scholar - K. Zhang, B. Jiang, and P. Shi, “Fast fault estimation and accommodation for dynamical systems,”
*IET Control Theory and Applications*, vol. 3, no. 2, pp. 189–199, 2009. View at: Publisher Site | Google Scholar - C. Langbort, R. S. Chandra, and R. D'Andrea, “Distributed control design for systems interconnected over an arbitrary graph,”
*IEEE Transactions on Automatic Control*, vol. 49, no. 9, pp. 1502–1519, 2004. View at: Publisher Site | Google Scholar - R. M. G. Ferrari, T. Parisini, and M. M. Polycarpou, “Distributed fault diagnosis with overlapping decompositions: an adaptive approximation approach,”
*IEEE Transactions on Automatic Control*, vol. 54, no. 4, pp. 794–799, 2009. View at: Publisher Site | Google Scholar - H. Zhang, D. Yang, and T. Chai, “Guaranteed cost networked control for T-S fuzzy systems with time delays,”
*IEEE Transactions on Systems, Man and Cybernetics C*, vol. 37, no. 2, pp. 160–172, 2007. View at: Publisher Site | Google Scholar - Z. Wang, Y. Wang, and Y. Liu, “Global synchronization for discrete-time stochastic complex networks with randomly occurred nonlinearities and mixed time delays,”
*IEEE Transactions on Neural Networks*, vol. 21, no. 1, pp. 11–25, 2010. View at: Publisher Site | Google Scholar - Z. Li and G. Chen, “Global synchronization and asymptotic stability of complex dynamical networks,”
*IEEE Transactions on Circuits and Systems II*, vol. 53, no. 1, pp. 28–33, 2006. View at: Publisher Site | Google Scholar - D. Ding, Z. Wang, B. Shen, and H. Shu, “${H}_{\infty}$ state estimation for discrete-time complex networks with randomly occurring sensor saturations and randomly varying sensor delays,”
*IEEE Transactions Neural Networks and Learning Systems*, vol. 23, no. 5, pp. 725–726, 2012. View at: Google Scholar - X. F. Wang and G. Chen, “Complex networks: small-world, scale-free and beyond,”
*IEEE Circuits and Systems Magazine*, vol. 3, no. 1, pp. 6–20, 2003. View at: Publisher Site | Google Scholar - Y. Liu, Z. Wang, J. Liang, and X. Liu, “Synchronization and state estimation for discrete-time complex networks with distributed delays,”
*IEEE Transactions on Systems, Man, and Cybernetics B*, vol. 38, no. 5, pp. 1314–1325, 2008. View at: Publisher Site | Google Scholar - Z. Wang, H. Zhang, and B. Jiang, “LMI-based approach for global asymptotic stability analysis of recurrent neural networks with various delays and structures,”
*IEEE Transactions on Neural Networks*, vol. 22, no. 7, pp. 1032–1045, 2011. View at: Publisher Site | Google Scholar - D. Gong, H. Zhang, Z. Wang, and D. Ma, “Synchronization criteria for an array of neutral-type neural networks with hybrid coupling: a novel analysis approach,”
*Neural Processing Letters*, vol. 35, no. 1, pp. 29–45, 2012. View at: Google Scholar - D. Gong, H. Zhang, Z. Wang, and B. Huang, “Novel synchronization analysis for complex networks with hybrid coupling by handling multitude Kronecker product terms,”
*Neurocomputing*, vol. 82, pp. 14–20, 2012. View at: Google Scholar - X. Jin and G. Yang, “Adaptive synchronization of a class of uncertain complex networks against network deterioration,”
*IEEE Transactions on Circuits and Systems I*, vol. 58, no. 6, pp. 1396–1409, 2011. View at: Publisher Site | Google Scholar - X. Jin, G. Yang, and W. Che, “Adaptive pinning control of deteriorated nonlinear coupling networks with circuit realization,”
*IEEE Transactions on Neural Networks and Learning Systems*, vol. 23, no. 9, pp. 1345–1355, 2012. View at: Google Scholar - Y.-W. Wang, J.-W. Xiao, and H. O. Wang, “Global synchronization of complex dynamical networks with network failures,”
*International Journal of Robust and Nonlinear Control*, vol. 20, no. 15, pp. 1667–1677, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. Zhang and L. Li, “Robust adaptive decentralized control for a class of networked large-scale systems with sensor network failures,” in
*Proceeding of the Chinese Control and Decision Conference (CCDC '11)*, pp. 2366–2370, Mianyang, China, 2011. View at: Google Scholar - H. Zhang and Y. Quan, “Modeling, identification, and control of a class of nonlinear systems,”
*IEEE Transactions on Fuzzy Systems*, vol. 9, no. 2, pp. 349–354, 2001. View at: Publisher Site | Google Scholar - M. M. Polycarpou and A. J. Helmicki, “Automated fault detection and accommodation: a learning systems approach,”
*IEEE Transactions on Systems, Man and Cybernetics*, vol. 25, no. 11, pp. 1447–1458, 1995. View at: Publisher Site | Google Scholar - A. T. Vemuri and M. M. Polycarpou, “Robust nonlinear fault diagnosis in input-output systems,”
*International Journal of Control*, vol. 68, no. 2, pp. 343–360, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH - A. B. Trunov and M. M. Polycarpou, “Automated fault diagnosis in nonlinear multivariable systems using a learning methodology,”
*IEEE Transactions on Neural Networks*, vol. 11, no. 1, pp. 91–101, 2000. View at: Publisher Site | Google Scholar - A. T. Vemuri, “Sensor bias fault diagnosis in a class of nonlinear systems,”
*IEEE Transactions on Automatic Control*, vol. 46, no. 6, pp. 949–954, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH

#### Copyright

Copyright © 2012 Zhanshan Wang 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.