Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2013, Article ID 282168, 10 pages
http://dx.doi.org/10.1155/2013/282168
Research Article

Intermittent Fault Detection for Uncertain Networked Systems

1Department of Automation, TNList, Tsinghua University, Beijing 100084, China
2School of Automation, Beijing Institute of Technology, Beijing 100081, China
3Beijing Engineering Research Center of Industrial Spectrum Imaging, School of Automation and Electrical Engineering, University of Science and Technology of Beijing, Beijing 100083, China
4Key Laboratory for Advanced Control of Iron and Steel Process, School of Automation and Electrical Engineering, University of Science and Technology of Beijing, Beijing 100083, China

Received 12 September 2013; Accepted 25 September 2013

Academic Editor: Hongli Dong

Copyright © 2013 Xiao He 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

This paper investigates intermittent fault detection problem for a class of networked systems with multiple state delays and unknown input. Polytopic-type parameter uncertainty in the state-space model matrices is considered. A novel measurement model is employed to account for both the random measurement delays and the stochastic data missing (package dropout) phenomenon, which are typically resulted from the limited capacity of the communication networks. We aim to design an uncertainty-dependent fault detection filter such that, for all unknown input, all possible parameter uncertainties, and all incomplete measurements, the error between residual and weighted fault is made as small as possible. By converting the addressed robust fault detection problem into an alternative robust filtering problem of a certain Markovian jumping system (MJS), a sufficient condition for the existence of the desired robust fault detection filter is derived. A residual evaluation within an incremental form is brought forward to make the whole method suitable for intermittent fault detection. A numerical example is utilized to demonstrate the effectiveness of the proposed approach.

1. Introduction

For traditional control systems with point-to-point data transmission, a variety of methods have been proposed to deal with the modeling, identification, estimation, and control problems [13]. During the past decades, the rapid developments in network technologies have led to more and more feedback control systems with control loops closed via digital communication channels. Compared with the traditional point-to-point wiring, in networked systems, serial communication networks are used to exchange information (reference input, plant output, control input, etc.) among control system components (sensors, controller, actuators, etc.) [4]. The use of the communication channels can reduce the costs of cables and power, simplify the installation and maintenance of the whole system, and increase the reliability, so network-based analysis and designs have many industrial applications such as in automobiles, manufacturing plants, aircrafts, and HVAC systems. However, the insertion of the communication channels raises new interesting and challenging problems such as network-induced delays or packets dropout, see [46] for some representative works.

With the increasing demand for higher performance, higher safety, and reliability standards, fault detection and isolation (FDI) has been an active field of research over the past decades [7, 8]. The main purpose of fault detection is to construct a residual signal which can then be compared with a predefined threshold. When the residual exceeds the threshold, the fault is detected and an alarm is generated. Among different approaches for residual generation, the model-based approaches to FDI problems for dynamic systems have received more attention. For example, in [9], thenorm of transfer function matrix from unknown input to residual has been designed to be small, while thenorm (or the smallest nonzero singular value) of transfer function matrix from fault to residual has been guaranteed to be large. In [10], the error between residual and weighted fault has been made as small as possible, and then the FDI problem can be solved by using thefiltering approach.

Due to the popularization of the using of network cables, it is necessary and interesting to consider the FDI problem for networked systems with network-induced delays or data missing, see [11, 12] and the references therein. Since network-induced delays and data missing (dropout) phenomenon are inherently random and time-varying [13], they have been modeled in various probabilistic ways [14]. One of the attractive approaches is to use binary switching sequence viewed as a Bernoulli distributed white sequence taking on values of 0 and 1, since such a representation is very effective to describe network-induced delays [15] or data missing [16]. Very recently, in [17], the network-induced delay and data dropout problems have been investigated in an integrated way within a unified framework and the robust filtering problem with polytopic uncertainties has been thoroughly studied. Note that in all the aforementioned results, it has been assumed that the delay or missing characteristics are statistically mutually independent from transfer to transfer. In [18], the fault detection problem for systems with missing measurements has been discussed by characterizing the residual dynamics by a discrete-time MJS. In [19], the diagnosis of intermittent faults in dynamic systems modeled as discrete event systems has been considered. So far, to the best of the authors’ knowledge, the robust intermittent fault detection problems in the presence of parameter uncertainty for networked systems with simultaneous measurement delays and data missing have not been fully investigated, which constitutes the main focus of this paper.

In this paper, intermittent fault detection problem for a class of uncertain networked systems with multiple state delays and incomplete measurement is investigated. A sequence varying in a Markov fashion is employed in the measurement model so that both the measurement delays and data missing can be simultaneously represented. Polytopic-type parameter uncertainty in state-space model matrices is considered. After augmenting the state, the addressed robust fault detection problem is converted to an equivalent robustfiltering problem for a certain Markovian jumping system (MJS), and a sufficient condition for the existence of the desired robust fault detection filter is brought forward. By introducing the new residual evaluation function within an incremental form, the proposed method can detect the possible intermittent fault.

Notation. The notations used throughout the paper are fairly standard.anddenote, respectively, the-dimensional Euclidean space and the set of allreal matrices.means thatis real symmetric and positive definite. The subscript “” denotes the matrix transpose.represents the occurrence probability of the event “”, and whenandare both stochastic variables,stands for the mathematical expectation of.is the space of all square-summable vector functions over, withbeing the standardnorm of, that is,.is the set of proper and stable rational functions with real coefficients. In symmetric block matrices, we use “” to represent a term that is induced by symmetry, andstands for a block-diagonal matrix. Matrices, if their dimensions are not explicitly stated, are assumed to be compatible for algebraic operations and, sometimes, the arguments of a function will be omitted in the analysis when no confusion can arise.

2. Problem Formulation and Preliminaries

Consider the following class of discrete-time linear networked systems with multiple delays in the state: wherestands for the state vector,is the unknown input belonging to, andis the fault to be detected.() are integer time delays.is the measured output vector, which may contain random communication delays and stochastic data missing induced by the limited capacity of the communication networks. All system matrices in (1) are assumed to have appropriate dimensions.is a given real initial sequence on, andstands for the Kronecker delta; that is, Furthermore,is a stochastic variable whose role is to determine, at time, the size of the occurred delay as well as the possibility of data missing. In this paper,is assumed to be a discrete-time homogeneous Markov chain taking values in the following finite state space: and stationary transition probability matrix, where

Remark 1. The assumption that the switching between difference modes abides by a Markovian chain seems realistic since, in network-based signal transmissions, time delays and data dropouts typically occur in a batch mode, and the status of the network varies slower than the sampling period and the characteristics of the network at a certain time is usually dependent on the superior time instant. The transition from one mode to another may obey certain probability distribution, and our “Markov jumping” assumption of describes this phenomenon properly.

Remark 2. In our measurement model, the event means that the measurements are ideally transmitted over the network without any delays or data missing, the event() corresponds to the case that the-step measurement delay occurs. Without loss of generality and for the convenience in denoting the transition probability matrix, we set the state in the state space of the Markov chain corresponding to the measurement missing situation to, someans that the measurements are missing andconsists of pure noise only.
In this paper, we are interested in the problem of robust fault detection for uncertain system described by (1) with incomplete measurements. The system matricesare assumed to be uncertain but belong to a known convex compact set of polytopic type,that is, whereis a given convex bounded polyhedral domain described byvertices as follows: where,, denotes theth vertex of the polytope.
Consider a full-order fault detection filter of the following form: whereis the filter state vector andis the so-called residual that is compatible with the fault vector. For, we denote matrices,, andas,, and. Our main aim is to make the error between residualand fault signalas small as possible.
For the purpose of fault detection, it is not necessary to estimate the fault. Sometimes one is more interested in the fault signal of a certain frequency interval, which can be formulated as the weighted fault as follows: whereis a prescribed weighting matrix.

Remark 3. Similar to [10], the introduction of a suitable weighting matrixcan limit the frequency interval of interest, and the system performance can then be improved. In fact, the use of weighted faultis more general than using the original fault, because if we impose, we can obtain.
Suppose a minimal realization ofis whereis the weighted fault state,is the original fault, andis the weighted fault.,,, andare assumed to be known real constant matrices with appropriate dimensions.
By defining and again, denoting matrices,,andas,,, and, we have the overall fault detection dynamics governed by the following system: where
After the above manipulations, the admissible sensor delays and data missing can be reformulated as the jumping parameters of a Markovian jumping system (11) with the same transition probability matrix.
The matrices,,, and,, are uncertain, but they belong to prescribed matrix polytopes, whereare given convex bounded polyhedral domain described byvertices as follows: anddenotes theth vertices of the polytopes, where Matrices,,,,,,, andare the same as defined in (12). Note that from (14), matrices,,, andare affinely dependent on the matrices, then for any, the uncertainty polytopes (13) have the same number of vertices, as well as the same combination coefficientswith (6), but different vertices for different Markovian model.
Recall the following definition of mean square stability for MJSs.

Definition 4 (see [20]). System (11) withis said to be mean square stable if for any initial conditionand initial distribution.

We further introduce the following definition.

Definition 5. System (11) with uncertain,(resp.,) is robust mean square stable if (11) is mean square stable for every,(resp.,).

Assumption 6. System (1) withandis assumed to be robust mean square stable.
With Definition 5, we can transform the robust fault detection filter design problem of system (1) to a robustfiltering problem for MJS (11). What we need to do here is to find the filter parameters,and() such that the augmented fault detection dynamics (11) is robust mean square stable and the infimum ofis made small in the feasibility of the following:
We further adopt a residual evaluation stage including an incremental evaluation functionand a thresholdof the following form: whereanddenotes the length of time window for evaluation function.
Based on (17), the occurrence of faults can be detected by comparingwithaccording to the following rule:

Remark 7. By introducing the residual evaluation function with an incremental form (17), one can detect the possible intermittent fault for an uncertain networked system by analyzing the residual signal once it is generated by the fault detection filter (7). The reason is that the residual evaluation signal (11) is decreased to a small value over time once the fault disappears. This nature makes the proposed method in this paper be used for intermittent fault detection.

3. Fault Detection Filter Design

In this section, we shall discuss the robust fault detection filter design problem of system (1), under the existence of parameter uncertainty (13). We introduce the following lemma which is useful in deriving our main results in the sequel.

Lemma 8. Consider system (1) with system uncertainty (5), for a given fault detection filter of the form (7), the augmented dynamic (11) is robust mean square stable and satisfies the constraint (16), if there exist matrice ,,,, andsuch that the following LMIs hold for all,, where,,, andare defined in (14),is defined in (12),,,,are defined in (9) and,

Proof. Considering the structure of,,, and, and imposing for all, it can be concluded from the bounded real Lemma in [21] that a sufficient condition is that there exist matrices,,such that the following LMIs hold for any, where,,,, andare defined in (12),,,, andare defined in (9),is the same defined in (21) and Following the steps as the proof of Theorem  1 in [22], it can be shown that LMIs (23) are feasible if and only if there exist matrices,,, andsatisfying
We are now in the position to prove that for system (1) with uncertainty (13), (20) ensures the robust mean square stable as well as the constraint (16) of the augmented dynamic (11). For an arbitrary fixed uncertain system with system matrices, one can always find a set of coefficients, , such that both (6) and (13) hold. Note that LMIs (20) are affine in the matrices,,,, and, multiplying suitable inequalities of (20) by appropriate scalarsand summing up, it can be readily shown that (25) holds for everywith a matrix,. By using the Bounded Real Lemma in [21], it follows that system (1) is robust mean square stable and (16) is satisfied. This concludes the proof.

Next, we give the robust fault detection filter design result for system (1) with system uncertainty (13).

Theorem 9. Consider system (1) with uncertain matrices (13), letbe a given scalar, there exists an admissible full-order robust fault detection filter of the form (7) ensuring that the overall augmented dynamics (11) is robust mean square stable and the constraint (16) is satisfied, if there exist matrices,,,,,,,,,, and, such that the following LMIs hold forand, whereis defined in (14) andis the same in (21). Moreover, if (26) is feasible, the parameters of the desired robust fault detection filter can be given by whereis any invertible matrix (e.g.,could be set as).

Proof. Considering (20), Letandbe two nonsingular matrices, and set Introducing the followind new matrices: where the entries,,, andare uniquely determined from the following relation: we further have the following relation
By defining we obtain
Performing congruence transformations to (20) by, define then, it can be easily shown that LMIs (26) together with the additional constraints (29) and (30) are equivalent to LMIs in (20). Hence, if there exist matrices,,,,,,,,, andsuch that LMIs (20) are feasible, the overall fault detection dynamic (11) is robust mean square stable and the constraint (16) is satisfied.
Furthermore, from LMIs (26), we have for, This indicates thatandare nonsingular and which implies thatis nonsingular and also ensures the existence of parameter matrices,, andin (28). The proof is completed.

Remark 10. In Theorem 9, uncertainty-dependent robust fault detection filter design result is provided, which reduces the conservatism than the uncertainty-independent results. If we impose to (26), the uncertainty-independent result can be recovered.

Remark 11. In most cases, we can know the size of the measurement delay or whether the data is missing at a certain time by using the time-stamp at the system node [5], and therefore the jumping parameters of the transformed MJS are accessible. In this sense, Theorem 9 provides us with network-status-dependent fault detection filter design methods. On the other hand, if the network status is not accessible; that is, the jumping parameters of the transformed MJS are unavailable, a network-status-independent result can be easily obtained by imposing in Theorem 9.
Note that (26) are LMIs over both the matrix variables and the prescribed scalar. This implies that (i) the robust full-order fault detection filter can be obtained from the solution of convex optimization problems in terms of LMIs, which can be solved via efficient interior-point algorithms [23]; (ii) the scalarcan be included as one of the optimization variables for LMIs (20), which makes it possible to obtain the minimum noise attenuation level bound for the fault detection dynamics (11). Then, the uncertainty-dependent suboptimal robust fault detection filter can be readily found by solving the following convex optimization problem.

Problem 12. Consider the parameter uncertainty (5), the sub-optimal robust fault detection filter for networked systems (1) with multiple state-delays and unknown inputs based on the idea of uncertainty dependence and network status dependence can be brought forward as follows:
For the problems mentioned above, the parameters of the sub-optimal robust fault detection filter can be determined by (28), and the sub-optimal robustattenuation level for fault detection dynamics is given by, whereare the sub-optimal solution of the corresponding convex optimization problems.
Here is a summary of the whole fault detection method for system (1).

Step 1. Determine the vertex of uncertain parameters of (1).

Step 2. Calculate the parameters of fault detection filter using Theorem 9.

Step 3. Get a appropriate threshold from experiments for a specific noise type.

Step 4. Generate a real-time residual signal from the fault detection filter designed in Step 2.

Step 5. Compare the evaluation function with the threshold and use the logic (19) to alarm a fault.

Remark 13. In the present work, a model-based approach is considered since there is a mathematical model for the system plant. However, when there is no such a model and only input and output data can be obtained in many complex systems, data-driven methods may work better than the model based ones since there is a leakage of prior information of system dynamics. Please see [24, 25] for typical data-driven fault detection methods.

4. A Numerical Example

To illustrate the effectiveness of the proposed method, we provide a numerical examples in this section. Consider system (1) with the following uncertain system parameters: whereis an uncertain real parameter satisfying. The initial state valuesare set to be. Let, so the state-space of the Markov chainis. The transition probability matrix is given by and the initial mode is set to be. For, the unknown inputis supposed to be a random noise uniformly distributed over, and the fault signalis of the following intermittent form:

The weighting matrix is supposed to be, with the following state space realization: whereandare shown in Figure 1.

282168.fig.001
Figure 1: Faultand weighting fault.

With the predefined parameters, from Theorem 9, Problem 12 can be solved by using the Matlab LMI toolbox [23]. As a result, the minimum noise attenuation level bound of the fault detection dynamic is, and the parameters of the sub-optimal fault detection filter in different modes are given by

If we impose the uncertainty-independent and network-status-independent method indicated in Remarks 10 and 11, a fault detection filter withcan be obtained.

Next, we consider the time-domain simulation using the obtained fault detection filter. we arbitrarily choose. Figure 2 shows the measurement mode with random delays and stochastic missing phenomenon.means that the measurement is missing, transmitted over the network ideally, and with one-step delay, respectively.

282168.fig.002
Figure 2: Measurement mode over network.

Figure 3 shows the generated residual signal, and the evolution ofis presented in Figure 4, where we choose the length of time window. After 400 times of simulations, we get a threshold. Figure 4 shows the real-time fault detection result:

282168.fig.003
Figure 3: Residual signal.
fig4
Figure 4: Evolution of.

From above inequations, we can observe the fault occurrence can be detected after 3, 5, 3 steps, while the disappearance of fault can be detected after 21, 25, 18 steps, respectively.

Remark 14. After introducing a novel residual evaluation function of the incremental type (17), one can detect not only the occurrence of a fault but also its disappearance. It can be observed from Figure 4 that the disappearance of the fault usually needs longer time to detect than the occurrence of faults. This is because the mathematical model of the system is based on the plant without the faults. This is of engineering significance and the proposed method in this paper can be applied to detect intermittent faults in many practical industrial processes.

Remark 15. The length of time window for evaluation functionin (17) is also a factor that can affect the fault detection performance. Longercan reduce the missing alarm rate and false alarm rate; however, it costs longer detection time. Shortcan provide a rapid fault detection but this results in larger missing alarm rate and false alarm rate.

5. Conclusions

In this paper, the robust intermittent fault detection problem has been investigated for a class of discrete-time uncertain networked systems with state delays and incomplete measurements. The random delay and stochastic missing phenomenon in the measurements have been simultaneously investigated. Polytopic-type parameter uncertainty in the state-space model matrices has been considered. By augmenting the states, the addressed robust fault detection problem has been converted to an auxiliary robustfiltering problem for a certain Markovian jumping system (MJS). A sufficient condition for the design of the desired robust fault detection filter has been established. A numerical example has been introduced to illustrate the effectiveness of the proposed methodology.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant 61074084 and 61074085 and an open Grant of the Key Laboratory of Hydro-Science and Engineering in Chongqing Jiaotong University under Grant SLK2010.

References

  1. S. Elmadssia, K. Saadaoui, and M. Benrejeb, “New delay-dependent stability conditions for linear systems with delay,” Systems Science and Control Engineering, vol. 1, no. 1, pp. 2–11, 2013. View at Google Scholar
  2. Y. Chen and K. A. Hoo, “Stability analysis for closed-loop management of a reservoir based on identification of reduced-order nonlinear model,” Systems Science and Control Engineering, vol. 1, no. 1, pp. 12–19, 2013. View at Publisher · View at Google Scholar
  3. M. Darouach and M. Chadli, “Admissibility and control of switched discrete-time singular systems,” Systems Science and Control Engineering, vol. 1, no. 1, pp. 43–51, 2013. View at Google Scholar
  4. W. Zhang, M. S. Branicky, and S. M. Phillips, “Stability of networked control systems,” IEEE Control Systems Magazine, vol. 21, no. 1, pp. 84–99, 2001. View at Publisher · View at Google Scholar
  5. J. Nilsson, Real-time control systems with delays [Ph.D. dissertation], Department Automatic Control, Lund Institute of Technology, Lund, Sweden, 1998.
  6. D. Yue, E. Tian, and Q.-L. Han, “A delay system method for designing event-triggered controllers of networked control systems,” IEEE Transactions on Automatic Control, vol. 58, no. 2, pp. 475–481, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  7. P. M. Frank, S. X. Ding, and B. Koppen-Seliger, “Current developments in the theory of FDI,” in Proceedings of the 4th IFAC Symposium Fault Detection Supervision Safety Technical Processes, Proceedings of SAFEPRO-CESS, pp. 16–27, Budapest, Hungary, 2000.
  8. M. Kinnaert, “Fault diagnosis based on analytical models for linear and nonlinear systems,” in Proceedings of the 15th International Workshop on Principles of Diagnosis, Proceedings of SAFEPROCESS, pp. 37–50, Washington, DC, USA, 2003.
  9. S. X. Ding, T. Jeinsch, P. M. Frank, and E. L. Ding, “A unified approach to the optimization of fault detection systems,” International Journal of Adaptive Control and Signal Processing, vol. 14, no. 7, pp. 725–745, 2000. View at Publisher · View at Google Scholar
  10. M. Zhong, H. Ye, P. Shi, and G. Wang, “Fault detection for markovian jump systems,” IEE Proceedings Control Theory and Applications, vol. 152, no. 4, pp. 397–402, 2005. View at Publisher · View at Google Scholar
  11. H. Fang, H. Ye, and M. Zhong, “Fault diagnosis of networked control systems (MED ’09),” in Mediterranean Conference on Control and Automation, Proceedings of SAFEPROCESS, pp. 1–12, Beijing, China, June 2006. View at Publisher · View at Google Scholar
  12. M. Nader and K. Khashayar, “Actuator fault detection and isolation for a network of unmanned vehicles,” IEEE Transactions on Automatic Control, vol. 54, no. 4, pp. 835–840, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  13. R. Luck and A. Ray, “An observer-based compensator for distributed delays,” Automatica, vol. 26, no. 5, pp. 903–908, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  14. L. Zhang, Y. Shi, T. Chen, and B. Huang, “A new method for stabilization of networked control systems with random delays,” IEEE Transactions on Automatic Control, vol. 50, no. 8, pp. 1177–1181, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  15. F. Yang, Z. Wang, Y. S. Hung, and M. Gani, “H control for networked systems with random communication delays,” IEEE Transactions on Automatic Control, vol. 51, no. 3, pp. 511–518, 2006. View at Publisher · View at Google Scholar · View at MathSciNet
  16. Z. Wang, D. W. C. Ho, and X. Liu, “Variance-constrained filtering for uncertain stochastic systems with missing measurements,” IEEE Transactions on Automatic Control, vol. 48, no. 7, pp. 1254–1258, 2003. View at Publisher · View at Google Scholar · View at MathSciNet
  17. X. He, Z. Wang, and D. Zhou, “Robust H filtering for networked systems with multiple state-delays,” International Journal of Control, vol. 80, no. 8, pp. 1217–1232, 2007. View at Publisher · View at Google Scholar
  18. P. Zhang, S. X. Ding, P. M. Frank, and M. Sader, “Fault detection of networked control systems with missing measurements,” in Proceedings of Asian Control Conference, pp. 1258–1263, Melbourne, Australia, 2004.
  19. O. Contant, S. Lafortune, and D. Teneketzis, “Diagnosis of intermittent faults,” Discrete Event Dynamic Systems: Theory and Applications, vol. 14, no. 2, pp. 171–202, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  20. O. L. V. Costa and M. D. Fragoso, “Stability results for discrete-time linear systems with Markovian jumping parameters,” Journal of Mathematical Analysis and Applications, vol. 179, no. 1, pp. 154–178, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. O. L. V. Costa and R. P. Marques, “Mixed H2/H-control of discrete-time Markovian jump linear systems,” IEEE Transactions on Automatic Control, vol. 43, no. 1, pp. 95–100, 1998. View at Publisher · View at Google Scholar · View at MathSciNet
  22. M. C. de Oliveira, J. C. Geromel, and J. Bernussou, “Extended H2 and H norm characterizations and controller parametrizations for discrete-time systems,” International Journal of Control, vol. 75, no. 9, pp. 666–679, 2002. View at Publisher · View at Google Scholar · View at MathSciNet
  23. S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, vol. 15 of SIAM Studies in Applied Mathematics, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, USA, 1994. View at Publisher · View at Google Scholar · View at MathSciNet
  24. S. Yin, H. Luo, and S. X. Ding, “Real-time implementation of fault tolerant control system with performance optimization,” IEEE Transactions on Industrial Electronics, 2013. View at Publisher · View at Google Scholar
  25. S. Yin, S. X. Ding, A. Gaghani, H. Hao, and P. Zhang, “A comparison study of basic data-driven fault diagnosis and process monitoring methods on the benchmark Tennessee Eastman process,” Journal of Process Control, vol. 22, no. 9, pp. 1567–1581, 2012. View at Publisher · View at Google Scholar