Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 697345, 9 pages

http://dx.doi.org/10.1155/2013/697345

## Fault Detection for Linear Discrete Time-Varying Systems with Measurement Packet Dropping

^{1}School of Electrical Engineering, University of Jinan, Jinan 250022, China^{2}School of Control Science and Engineering, Shandong University, Jinan 250061, China^{3}School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798

Received 24 January 2013; Accepted 8 April 2013

Academic Editor: Wuquan Li

Copyright © 2013 Yueyang Li 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

The fault detection (FD) problem for linear discrete time-varying (LDTV) systems with measurement packet dropouts is considered. The objective is to design a new observer-based fault detection filter (FDF) as a residual generator through employing packet dropout information on the measurement sequence. Based on some new defined input-to-output operators, the FD problem is formulated in a framework of maximizing stochastic or performance index. By introducing an adjoint-operator-based optimization method, the analytical optimal solution can be derived in terms of solving a modified Riccati equation. A numerical example is provided to demonstrate the effectiveness of the proposed approach.

#### 1. Introduction

The problem of model-based fault detection and isolation has attracted much attention during the last two decades and bred numerous results. One of the most useful strategies is the observer-based fault detection (FD) approach which involves constructing an observer as a fault detection filter (FDF) for generating a residual signal that is sensitive to fault while insensitive to unknown input [1–7]. Among the developments of FD techniques, the filtering scheme which formulates the FD problem into the framework of filtering is widely used for systems with -norm bounded unknown inputs and faults [8–11]. Another efficient way is the so-called optimization scheme which employs norm as a robustness measure while the or norm is introduced as a fault sensitivity measure, and, then, the FD problem can be formulated as an or optimization problem. In [12–14], the global optimal solution to and optimization FD problem is given for linear time-invariant (LTI) systems by coprime factorization and linear matrix inequality (LMI) techniques, respectively. Recent years, some new results are dedicated to this problem for linear time-varying (LTV) systems, especially for linear discrete time-varying (LDTV) systems, under the background that real systems are intrinsically time varying and most industrial processes are operated by digital devices. In [15, 16], the unified approach in [12] is extended to optimization FD problem for linear discrete-time periodic systems. Finite horizon and FD formulation for LDTV systems is proposed in [17–19] and the optimal solution is derived by solving a standard Kalman-like Riccati equation.

On another frontline, research on systems with intermittent measurements has gained growing interests due to that the communication networks are widely utilized in control systems [20–25]. In the literature of control and filtering, stochastic parameters such as the Bernoulli distributed variables and Markov chain are introduced to model the phenomenon of missing measurements, which can be divided into two fundamentally different categories. One is the UDP-like protocol-based model where there is no acknowledgement mechanism of successful delivery of data packets. The other model is based on the TCP-like protocol where successful transmissions of packets are acknowledged at the receiver [20]. With the increasing demands for system safety and reliability, FD for systems with measurement packet dropping has become one of the most critical issues and many related results have emerged with the help of LMI technique. In [26], a robust FDF is designed for LTI systems with missing measurements. In [27–29], the same idea is extended to Markovian jump singular systems, T-S fuzzy systems, and discrete-time switched systems under intermittent measurements condition, respectively. The filtering FD scheme is also applied to networked control systems with both induced delays and incomplete measurements through augmenting retarded vectors into delay-free variables and numerical solutions are derived, we refer to [30–34] and references therein. However, it should be noted that most of the existing results are based on filtering scheme under the assumption that the measurements are transmitted using UDP-like protocol-based model and only infinite-horizon cases are considered. There remain some problems to be figured out in the open area; for example, the existing LMI-based approach does not efficiently work for time-varying systems owing to large computational burden brought by optimization iterative algorithm at each time instant, and, consequently, how to detect a fault for LDTV systems with missing measurements in the TCP-like protocols scenario.

Motivated by the above discussions, the FD problem for LDTV systems with measurement packet dropout in TCP-like protocol frame will be investigated in this paper. The main contributions of the paper are twofold.(1)A more generalized FD problem description in the framework of optimizing stochastic or performance index is formulated by defining input-to-output operator that maps from fault or unknown input to residual signal. Then, a new observer-based FDF is proposed for generating residuals based on the acknowledgement mechanism of data packet transmission.(2)An analytical solution to the aforementioned problem is explicitly expressed by solving a modified Riccati equation through a proposed adjoint-operator-based optimization method.

The remainder of this paper is organized as follows. In Section 2, the FD problem is formulated in the sense of maximizing a stochastic sensitivity/robustness ratio. In Section 3, a unified solution to optimizing stochastic or performance index problem is derived by using the proposed adjoint-operator-based optimization scheme. Finally, an illustrative example is given in Section 4 to demonstrate the effectiveness of our proposed approach.

*Notations. *Throughout this paper, for a matrix , and stand for the transpose and inverse of , respectively. means the set of -dimensional real vectors. and denote identity matrix and zero matrix with appropriate dimensions, respectively. means that is positive (negative) definite. E represents the mathematical expectation of . denotes the Euclidean norm of . stands for the deterministic -norm of with , while for the stochastic case with , where is a positive integer. gives the definition of the inner product on a Hilbert space for -norm bounded vector and with appropriate dimensions, where is a positive-definite initial weighting matrix. means the occurrence probability of the event “”.

#### 2. Problem Formulation

Consider the following LDTV system with measurement packet dropping: where , , , , and denote the state, measurement output, control input, unknown input, and fault, respectively. , , and are known time-varying matrices with appropriate dimensions. and are -norm bounded signals. depicts the measurement packet dropouts and is assumed to be a scalar Bernoulli distributed binary stochastic variable; that is, where is a known constant which can be obtained as prior knowledge through statistical test method [32].

Before describing the underlying problems, the following definition and assumptions are introduced.

*Definition 1. *System (1) is said exponentially stable in mean square sense (ESMS) with , , and , if there exist and such that
for any initial condition .

*Assumption 2. * is stabilizable.

*Assumption 3. * is exactly observable.

*Assumption 4. * is available at each time instant .

The primary object of FD is to construct an FDF for generating a residual signal which is robust to unknown input while sensitive to fault . Note that the indicator is known due to Assumption 4, and hence can be involved in the following observer-based FDF: where is an estimation of , is the generated residual, is a guess of initial state, is the observer gain matrices, and is the (regular) post-filters to be determined.

Let ; it follows from (1) and (4) that

Denote where and .

Similar to the perturbed operator defined by -norm of input and output signals in [35], a linear operator that maps and a linear operator that maps can be defined as follows by the linearity of (1): with where is a positive-definite initial weighting matrix.

Based on the definitions above, we are now in the position to introduce an auxiliary FDF design problem for residual generation: find a suitable observer gain matrix and a regular post-filter such that system (5) is ESMS and satisfies the following performance:

After designing the FDF, the remaining task is to evaluate the generated residual. In this paper, we adopt the residual evaluation function in the following forms: where denotes the length of the evaluation time window [2]. The corresponding threshold is chosen as and hence the occurrence of faults can then be recognized based on the following rule:

*Remark 5. *Assumptions 2 and 3 are given without loss of generality to guarantee the stability of the proposed FDF. For more details about the definitions on stabilizable and exactly observable, we refer to [24, 36] and references therein. With the aid of time-stamp technique, Assumption 4 is also reasonable due to the acknowledgements of successful transmissions of data packets for TCP-like protocols. Furthermore, under the condition that Assumption 4 holds, our proposed FDF is different from the ones given in [26–29] since the new FDF exploits additional information on the observation sequence.

*Remark 6. *It should be emphasized that since stochastic characteristic is introduced by , the existing FD problem formulation and results in [17–19] for deterministic LTV systems would not be reasonably established. Our proposed performance index or is a generalized stochastic version comparing to that of the existing contributions. In what follows, a novel adjoint-operator-based optimization approach will be proposed with the well-defined linear operator and , and we will show that detailed interpretation of and is not necessary owing to our proposed method.

#### 3. Main Results

Before deriving the main results of this paper, the following definitions and lemmas should be given.

*Definition 7 (see [37]). *Let denote an operator or a system mapping from -norm bounded space to -norm bounded space . An operator is called to be the adjoint operator of from space to if , for all and , where stands for the initial vector.

*Definition 8 (see [38]). *Let denote an operator or a system mapping from -norm bounded input space to -norm bounded output space ; then, is said to be coisometric if for all .

Lemma 9. *For two stochastic operators and , where , and are -norm bounded signals, one has
*

*Proof. *Let and ; we can conclude that
which completes the proof.

*Lemma 10. Consider the following residual generators:
where is the observer gain matrix such that system (15) is ESMS and is the post-filter. Then
where is an operator that maps .*

*Proof. *By applying Lemma 1 in [18] to (15), we have that for the following residual generators
where is the observer gain matrix that ensures the stability of (15); an operator that guarantees exists, which can be realized as
The regularity of indicates there exists such that
where denotes the left inverse of , and then we have
Thus, from (18), can be represented as
where . This completes the proof.

*Lemma 11. For system (5), consider the operator that maps which is realized by the following discrete-time system:
where , , , and . Let be the adjoint operator of ; then is coisometric if there exists a positive definite matrix satisfying the following equation:
*

*Proof. *Denote . From [39], the state space representation of can be obtained as

For (24), define

Moreover, based on (22) and (24), we have

Remembering that and , then, from (26) and Definition 8, if
that is, (23) holds, is coisometric, which completes the proof.

*In view of the proposed lemmas above, we now present solutions to the optimization problem (9).*

*Theorem 12. If Assumptions 2 and 3 hold, then the following matrix pair
with
give a solution to the auxiliary FDF design problem and is the solution to the following equation with :
*

*Proof. *We first give the proof of the optimal solution to (9). Let be the optimal generated residual in the sense of (9). Since system (1) is linear, by applying Lemma 10, we know that there exists an operator such that
where and , which shows that

On the other hand, consider the operator that maps in (5); we have
where , which concludes that

Similarly, we have
where . Hence,

Based on (24) and according to [40, 41] and Definition 8; that is, if is coisometric, we have that
and the following inequality is immediately established based on Lemma 9:
which gives the optimal value of maximizing the performance index .

Following the same idea, we can prove that

Furthermore, from Lemma 11, by solving the following equations in (23):
we can obtain (28) and consequently
turns to (30).

Next, we discuss the stability of the proposed FDF. In fact, by introducing satisfying
where and means the Kronecker delta function, an artificial stochastic backward control system associated with (24) can be defined as follows:
where , , , and . From [36], if Assumptions 2 and 3 hold, there exists a unique positive solution to (30) in the sense of finding a linear quadratic (LQ) control law to (43). Furthermore, following [42], the proposed FDF is ESMS. This completes the proof.

*Theorem 12 addresses a static solution to the two-objective optimization problem (9). A more generalized dynamic operator form of this theorem can be derived as follows.*

*Theorem 13. Let be any observer gain that guarantees the stability of the proposed FDF (4); then that maps in the following way presents a dynamic post-filter for the auxiliary FDF design problem:
where and are determined in (28) and (30).*

*Proof. *The proof is straightforward as a result of Lemma 10. Hence, it is omitted here.

*Theorem 12 can be readily extended to the stationary FDF design over an infinite horizon. In this case, the system matrices , and will be replaced by constant matrices , , and , respectively. The following result, which is the stationary version of Theorem 12 for LTI systems with measurement packet dropping, can be obtained by redefining in Lemma 11 with a constant and replacing the Riccati difference equation (30) by a corresponding algebraic Riccati equation. In addition, from [24, 42], the ESMS of the proposed FDF can be guaranteed when Assumptions 2 and 3 are satisfied.*

*Corollary 14. Consider system (1) with constant system matrices and zero initial state. If Assumptions 2 and 3 hold, then the following matrix pair
with
give a solution to (9) in infinite horizon and is the solution to the following equation:
*

*Remark 15. *Theorem 12 gives unified analytical solutions to the associated FDF design problem of optimizing the defined stochastic and performance indices. It can easily conclude the following.(i)If ; that is, there measurement packet dropping occurs, the solution to the FDF design problem is derived based on a modified Riccati equation (30).(ii)If ; that is, there is no measurement packet dropout, our results will coincide with the one given in [18, 19], where the modified Riccati equation (30) will reduce to a standard Kalman-like filtering Riccati equation.

*Remark 16. *Roughly speaking, the measurement packet dropping phenomenon depicted by the Bernoulli random variable can be divided into three categories, that is, measurement packet dropouts, stochastic finite step delay, and multiple packet dropouts [20]. It can be concluded that if the packet dropouts indicator (as in this paper) is on-line known, our proposed algorithm can be applied to these scenarios.

*Remark 17. *It should be pointed out that for system (1), the existing filtering-based fault estimation approach in [39] cannot achieve the fault detection goal since the fault distribution matrix in the measurement is zero. However, our solution is independent on the fault distribution matrices, which infers its more generalized availability. In addition, comparing to the LMI-based fault estimation approaches given in [26–29], the proposed results are in analytical form and load less computational burden.

*4. Numerical Example*

*4. Numerical Example**To illustrate the effectiveness of the proposed method, we consider (1) with the following parameter matrices:
*

*
The missing measurement rate is displayed in Figure 1 with mean and the unknown input is simulated in Figure 2. A stepwise fault signal is assumed to be as follows:
*

*Set . We design the FDF by directly applying Theorem 12. The generated residual is shown in Figure 3 and its residual evaluation function in (10) is displayed in Figure 4 with the corresponding threshold = 0.7457. It can be seen from the simulation results that which implies that the residual can deliver fault alarm in 1 second after it occurs.*

*5. Conclusions*

*5. Conclusions**In this paper, the problem of FD for LDTV systems with measurement packet dropping has been dealt with. Under the condition that the packet dropout indicator is on-line available, an observer-based FDF has been employed as a residual generator through exploiting the information brought by the packet dropout indicator. The design of FDF has been formulated in the framework of maximizing stochastic or performance index. An analytical optimal solution has been derived by solving a modified Riccati equation based on the proposed adjoint operator optimization method. The achieved result has been illustrated by a numerical example.*

*Acknowledgments*

*Acknowledgments**This work is partially supported by National Natural Science Foundation of China (no. 61203083 and no. 61074021), the Doctoral Foundation of University of Jinan (no. XBS1242).*

*References*

*References*

- J. Chen and R. Patton,
*Robust Model-Based Fault Diagnosis for Dynamic Systems*, Kluwer Academic Publishers, Boston, Mass, USA, 1999. - S. X. Ding,
*Model-Based Fault Diagnosis Techniques*, Springer, Berlin, Germany, 2008. - B. Kulcsár and M. Verhaegen, “Robust inversion based fault estimation for discrete-time LPV systems,”
*IEEE Transactions on Automatic Control*, vol. 57, no. 6, pp. 1581–1586, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - Z. Wang, H. Guan, and C. 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. View at Publisher · View at Google Scholar · View at MathSciNet - J. Yu, G. Sun, and H. R. Karimi, “Fault-reconstruction-based cascaded sliding mode observers for descriptor linear systems,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 623426, 20 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Yin, P. Shi, F. Liu, and J.-S. Pan, “Gain-scheduled fault detection on stochastic nonlinear systems with partially known transition jump rates,”
*Nonlinear Analysis. Real World Applications*, vol. 13, no. 1, pp. 359–369, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - L. Wu, X. Yao, and W. X. Zheng, “Generalized ${H}_{2}$ fault detection for two-dimensional Markovian jump systems,”
*Automatica*, vol. 48, no. 8, pp. 1741–1750, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - J. Chen and R. Patton, “Standard ${H}_{\infty}$ filtering formulation of robust fault detection,” in
*Proceedings of SAFEPROCESS*, Budapest, Hungary, 2000. - M.-Y. Zhong, S. Liu, and H.-H. Zhao, “Krein space-based ${H}_{\infty}$ fault estimation for linear discrete time-varying systems,”
*Acta Automatica Sinica*, vol. 34, no. 12, pp. 1529–1533, 2008. View at Google Scholar · View at MathSciNet - Y. Li and M. Zhong, “Fault detection filter design for linear discrete time-varying systems with
multiplicative noise,”
*Journal of Systems Engineering and Electronics*, vol. 22, no. 6, pp. 982–990, 2011. View at Google Scholar - H. Dong, Z. Wang, and H. Gao, “Fault detection for Markovian jump systems with sensor saturations and randomly varying nonlinearities,”
*IEEE Transactions on Circuits and Systems I*, vol. 59, no. 10, pp. 2354–2362, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - 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 Google Scholar - J. Wang, G. Yang, and J. Liu, “An LMI approach to ${H}_{-}$ index and mixed ${H}_{-}/{H}_{\infty}$ fault detection
observer design,”
*Automatica*, vol. 43, no. 9, pp. 1656–1665, 2008. View at Google Scholar - A. Casavola, D. Famularo, and G. Franzè, “Robust fault detection of uncertain linear systems via quasi-LMIs,”
*Automatica*, vol. 44, no. 1, pp. 289–295, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - P. Zhang, S. X. Ding, G. Z. Wang, and D. H. Zhou, “Fault detection of linear discrete-time periodic systems,”
*IEEE Transactions on Automatic Control*, vol. 50, no. 2, pp. 239–244, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - P. Zhang, S. X. Ding, and P. Liu, “A lifting based approach to observer based fault detection of linear periodic systems,”
*IEEE Transactions on Automatic Control*, vol. 57, no. 2, pp. 457–462, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - X. Li,
*Fault Detection Filter Design for Linear Systems [Ph.D. thesis]*, Louisiana State University, 2009. - M. Zhong, S. X. Ding, and E. L. Ding, “Optimal fault detection for linear discrete time-varying systems,”
*Automatica*, vol. 46, no. 8, pp. 1395–1400, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - M. Zhong and S. X. Ding, “A projection-based method of fault detection for linear discrete time-varying system,”
*International Journal of Systems Science*, 2011. View at Google Scholar - J. Hespanha, P. Naghshtabrizi, and Y. Xu, “A survey of recent results in networked control systems,”
*Proceedings of the IEEE*, vol. 95, no. 1, pp. 138–162, 2007. View at Google Scholar - L. Bakule and M. de la Sen, “Decentralized resilient ${H}_{\infty}$ observer-based control for a class of uncertain interconnected networked systems,” in
*Proceedings of the American Control Conference*, Baltimore, Md, USA, 2010. - K. You, M. Fu, and L. Xie, “Mean square stability for Kalman filtering with Markovian packet losses,”
*Automatica*, vol. 47, no. 12, pp. 2647–2657, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Mo and B. Sinopoli, “Kalman filtering with intermittent observations: tail distribution and critical value,”
*IEEE Transactions on Automatic Control*, vol. 57, no. 3, pp. 677–689, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - H. Zhang, X. Song, and L. Shi, “Convergence and mean square stability of suboptimal estimator for systems with measurement packet dropping,”
*IEEE Transactions on Automatic Control*, vol. 57, no. 5, pp. 1248–1253, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Zhang, C. Liu, and X. Mu, “On stochastic finite-time control of discrete-time fuzzy systems with packet dropout,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 752950, 18 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - H. Gao, T. Chen, and L. Wang, “Robust fault detection with missing measurements,”
*International Journal of Control*, vol. 81, no. 5, pp. 804–819, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Zhao, J. Lam, and H. Gao, “Fault detection for fuzzy systems with intermittent measurements,”
*IEEE Transaction on Fuzzy Systems*, vol. 17, no. 2, pp. 398–410, 2009. View at Google Scholar - X. Yao, L. Wu, and W. X. Zheng, “Fault detection filter design for Markovian jump singular systems with intermittent measurements,”
*IEEE Transactions on Signal Processing*, vol. 59, no. 7, pp. 3099–3109, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - D. Du, B. Jiang, and P. Shi, “Fault detection for discrete-time switched systems with intermittent measurements,”
*International Journal of Control*, vol. 85, no. 1, pp. 78–87, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - X. He, Z. Wang, Y. Ji, and D. Zhou, “Robust fault detection for networked systems with distributed sensors,”
*IEEE Transaction on Aerospace and Electronic Systems*, vol. 47, no. 1, pp. 166–177, 2008. View at Google Scholar - X. He, Z. Wang, and D. H. Zhou, “Robust fault detection for networked systems with communication delay and data missing,”
*Automatica*, vol. 45, no. 11, pp. 2634–2639, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - D. Zhou, X. He, Z. Wang et al., “Leakage fault diagnosis for an internet-based three-tank system: an experimental study,”
*IEEE Transaction on Control Systems Technology*, vol. 20, no. 4, pp. 857–870, 2011. View at Google Scholar - X. Wan and H. Fang, “Fault detection for networked nonlinear systems with time delays and packet dropouts,”
*Circuits, Systems, and Signal Processing*, vol. 31, no. 1, pp. 329–345, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - X. Wan, H. Fang, and F. Yang, “Fault detection for a class of networked nonlinear systems subject to
imperfect measurements,”
*International Journal of Control, Automation, and Systems*, vol. 10, no. 2, pp. 265–274, 2012. View at Google Scholar - W. Zhang, Y. Huang, and L. Xie, “Infinite horizon stochastic ${H}_{2}/{H}_{\infty}$ control for discrete-time systems with state and disturbance dependent noise,”
*Automatica*, vol. 44, no. 9, pp. 2306–2316, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Huang, W. Zhang, and H. Zhang, “Infinite horizon LQ optimal control for discrete-time stochastic systems,” in
*Proceedings of the 6th World Congress on Intelligent Control and Automation*, Dalian, China, 2006. - M. Green and D. J. N. Limebeer,
*Linear Robust Control*, Prentice Hall, Upper Saddle River, NJ, USA, 1995. - X. Li and K. Zhou, “A time domain approach to robust fault detection of linear time-varying systems,”
*Automatica*, vol. 45, no. 1, pp. 94–102, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Li and M. Zhong, “On designing robust ${H}_{\infty}$ fault detection filter for linear discrete time-varying
systems with multiple packet dropouts,”
*Acta Automatica Sinica*, vol. 36, no. 12, pp. 1788–1796, 2010. View at Google Scholar - E. Gershon, U. Shaked, and I. Yaesh, “${H}_{\infty}$ control and filtering of discrete-time stochastic systems with multiplicative noise,”
*Automatica*, vol. 37, no. 3, pp. 409–417, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - A. Bobrowski,
*Functional Analysis for Probability and Stochastic Processes. An Introduction*, Cambridge University Press, Cambridge, UK, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - V. Dragan, T. Morozan, and A.-M. Stoica,
*Mathematical Methods in Robust Control of Discrete-Time Linear Stochastic Systems*, Springer, New York, NY, USA, 2010. View at Publisher · View at Google Scholar · View at MathSciNet

*
*