Variance-Constrained Robust Estimation for Discrete-Time Systems with Communication Constraints

Wang, Baofeng; Guo, Ge; Gao, Xiue

doi:https://doi.org/10.1155/2014/980753

Mathematical Problems in Engineering

On this page

Abstract Introduction Preliminaries Proofs Numerical Example Conclusion Acknowledgments References Copyright Related Articles

Special Issue

Sensor/Actuator Networks and Networked Control Systems

View this Special Issue

Research Article | Open Access

Volume 2014 | Article ID 980753 | https://doi.org/10.1155/2014/980753

Variance-Constrained Robust Estimation for Discrete-Time Systems with Communication Constraints

Baofeng Wang,^1,2Ge Guo,²and Xiue Gao¹

Academic Editor: Fuzhong Nian

Received06 Sept 2013

Accepted22 Dec 2013

Published14 Jan 2014

Abstract

This paper is concerned with a new filtering problem in networked control systems (NCSs) subject to limited communication capacity, which includes measurement quantization, random transmission delay, and packets loss. The measurements are first quantized via a logarithmic quantizer and then transmitted through a digital communication network with random delay and packet loss. The three communication constraints phenomena which can be seen as a class of uncertainties are formulated by a stochastic parameter uncertainty system. The purpose of the paper is to design a linear filter such that, for all the communication constraints, the error state of the filtering process is mean square bounded and the steady-state variance of the estimation error for each state is not more than the individual prescribed upper bound. It is shown that the desired filtering can effectively be solved if there are positive definite solutions to a couple of algebraic Riccati-like inequalities or linear matrix inequalities. Finally, an illustrative numerical example is presented to demonstrate the effectiveness and flexibility of the proposed design approach.

1. Introduction

State estimation has long been a significant interesting problem in the control and information fields. In recent years, with technological advances in MEMS, DSP capabilities, computing, and communication technology, networked control systems (NCSs) are built massively, in which sensors, actuators, and controllers are not physically collocated [1, 2]. Hence, there are various uncertainties caused by limited communication capacity when exchanging information via a digital communication network [3–5]. With NCSs becoming more and more popular in various applications, state estimation problem with limited communication capacity has attracted recurring interests. Generally, there are three important issues that should be addressed for communication-based systems, that is, signal transmission delay, packet loss, and quantization effect. For example, the packet loss has been described by the binary Bernoulli distribution approach [6, 7] as well as Markovian jump approach [8], packet delay has been investigated as constant delay [9] and bounded delay [5] along with unbounded packet delay [10], and the quantization has been dealt with via a sector bound approach where the quantization errors are regarded as a class of uncertainties [11]. However, the aforementioned pieces of literature still suffer from some limitations such as the fact that most results only investigate one or two aspects of the communication issues mentioned above, especially quantization effect has been paid little attention in network-based system; some papers are concentrated on controller design problem, while state estimation problem has been neglected, which is equally important in applications and so on.

For networked environment, much effort has been made towards proposing useful alternative methods in different contexts [6, 8, 10, 12–15], among which filtering and robust filtering approaches are prominent to improve the robustness. Generally, there are two popular approaches used to design a robust filter, that is, Riccati equations approach [16, 17] and linear matrix inequalities approach [18]. However, in some practical engineering such as the tracking of a maneuvering target, the filtering performance requirements are described as the upper bounds in the error variances of estimation, that is, the steady-state error variance is not required to be the minimum but should not be more than the specified upper bound constraint. Based on the error covariance assignment (ECA) theory [19], which provided a closed form solution for directly placing the specified steady-state estimation error covariance to a linear system, variance-constrained filtering method has been developed and applied to some more realistic system. For example, sampled-data system was studied in [20] and bilinear systems in [21]. In [7, 22], Wang et al. considered the filtering problem for linear uncertain discrete-time system with missing measurements and randomly varying sensor delay, respectively. However, quantization effect and transmission delay or packet dropouts were neglected.

In this paper, we investigate the filtering problem for discrete-time system in networked control systems (NCSs) subjected to limited communication capacity including the three aspects previously discussed. Through transmitting the communication limitations into a stochastic uncertainty system representation, variance-constrained filter problems are considered. It is shown that the filtering can effectively be derived if there are positive definite solutions to a couple of algebraic Riccati-like inequalities or linear matrix inequalities. Relative to the traditional optimal robust or filtering, for variance-constrained filtering approach, there still remains much freedom to be used for achieving other desired performance requirements after assigning to the system a specified variance upper bound. An illustrative numerical example is used to demonstrate the usefulness and flexibility of the proposed design approach.

The rest of this paper is arranged as follows. In Section 2, the robust filter design problem for discrete-time systems with communication limitations is formulated. Some preliminary results are given in Section 3, and the solution of this problem is developed in Section 4. In Section 5, the effectiveness of the proposed theory is demonstrated by an illustrative example. Some concluding remarks and future research directions are draw in Section 6.

2. Problem Formulation and Preliminaries

Consider the following linear discrete-time stochastic system: where is the state vector, is the measured output vector, and and are the process noise and measurement noise, respectively. , , and are known real constant matrices with appropriate dimensions.

We make the following assumptions on the statistical properties of the initial state, process noise, and measurement noise, which are also standard in Kalman filtering.

Assumption 1. For all integers ,(1), (2), (3)We point out that if and/or are colored noises with known statistics, they can be prewhitened.

Furthermore, throughout the paper, the following assumption is made.

Assumption 2. The matrix is Schur stable (i.e., all eigenvalues of are located within the unit circle in the complex plane).
As aforementioned analysis, generally, three effects need to be taken into consideration: measurement quantization, packet transmission delays, and packet loss. We model them mathematically as follows.

(i) Measurement Quantization. It is assumed that the sampled measurements of are first quantized via a quantizer, encapsulated into a packet, and then transmitted through a digital communication network. In this paper, we consider the logarithmic quantizer. According to [3, 23, 24], the quantizer is denoted as , which is symmetric; that is, , . For each , the set of quantized levels is characterized by and the associated quantizer is defined as follows: where Then, we have where denotes the measurements after quantization, and Furthermore, define For convenience of later analysis, we make the following assumption: It is clear that satisfies

(ii) Randomly Varying Packet Delays and Dropouts. In the paper, we assume that the digital communication channel is subjected to the randomly bounded packet delays and dropouts, which is modeled as following where represents the system output subject to randomly varying delays and dropouts, and the stochastic variable is a Bernoulli distributed white sequence taking values on 0 and 1 with following statistical properties.

Assumption 3. For stochastic variable , we have(1), , where is a known positive scalar.(2) is independent of , and .(3) are mutually independent.
Directly, from (1), the following formula can be derived, .
By substituting (6) into (11), we obtain And a compact representation of systems (1) and (12) can be given as follows: where , and For compact systems (13)-(14), observe that the matrices and contain stochastic parameter ; in view of Assumptions 1 and 3, we have Lemma 4.

Lemma 4. For systems (13) and (14), the initial state is uncorrelated with both and and has the following properties:(1), ,(2), (3) and are mutually independent,(4).

Remark 5. Mode (11) properly represents the unreliable networked communication where the time delay and packet dropouts are unavoidable, and the mode is based on the rationale that the induced data latency from the sensor to the controller is restricted not to exceed the sampling period. We can see that if , that is, , it denotes randomly varying packet delay, which was first introduced in [25], and subsequently it has been studied in paper [22]. And if , that is, , then there is no packet delay and loss. While and there is no packet received by the estimator site; that is, the packet is lost. Also, mode (11) also can be regarded as a special case of paper [26].

In this paper, for systems (13) and (14), the full-order filter considered is described by where stands for the state estimate of the stochastic systems (13) and (14) and the constant matrices and are filter parameters to be determined.

Define the steady-state estimation error covariance as follows

The objective of this paper is to determine the filter parameters and such that, for communication limitations including the three aspects previously analyzed. (1) the state of compact system (13) is mean square bounded and (2) the steady-state error covariance of the sequence satisfies where means the steady-state variance of the th error state and denotes the prespecified steady-state error-estimation variance constraint on the th state.

Remark 6. In (18), individual upper bound constraint, which can be obtained according to the engineering requirement, is imposed on individual steady-state estimation error variance. Note that conventional minimum variance filtering method aims to minimize a weighted scalar sum of the estimation error variances and is therefore not able to ensure that the multiple variance requirements will be satisfied [27].

3. Some Preliminary Results

In this section, we will first establish the state-space model of augmented system followed from compact systems (13) and (14) and filter (16) and then give the upper bound of the variance of the state estimation error under the three aspects of communication limitations.

Define the estimation error by Then, from (14) and (16), we have and, subsequently,

Again, we define an augmented system, whose state is defined as Then, combining (13) and (21), the dynamics of can be expressed by where , , and , where

It is clear that augmented system (23) comprises not only uncertainties but also stochastic matrix sequences and . This makes the augmented system (23) a stochastic parameter system, which reflects the characteristic of limited communication capacity. Now, we denote the state covariance matrix of augmented system (23) by

Since augmented system (23) involves both uncertainties and stochastic matrix sequences and caused by limit communication capacity, the computation of state covariance matrix is complex because of the complexity of the statistics analysis. Hence, the Lyapunov equation that governs the evolution of the covariance matrix is given by the following.

Lemma 7. The evolution Lyapunov equation of the covariance matrix from (23) can be written as where , , , , , , and satisfy (25) and , where

Proof. For the sake of convenience of discussion, we introduce the following new stochastic sequences: then
From Assumption 3, we have the statistical properties relating to stochastic variables , , , and as follows:
Now, in order to find the evolution Lyapunov equation of from augmented system (23), we have to solve two problems: (1) and (2) . (1)For , we have (i)For the first expression on the right hand of the above equation, Noticing that , where , we have So, (ii)Similarly, for the other expressions on the right hand of (32), we obtain that Hence, we have (2)For , from the statistical properties of and , it is easy to obtain that Consequently, (27) holds. This completes the proof of this lemma.

Define the steady-state covariance of (27) as follows: Then, from [7, 22, 28], we have the following useful lemma.

Lemma 8. There exists a unique symmetric positive semidefinite solution to the following discrete-time equation: that is, the convergence of in (27) is guaranteed to a constant value , and the state of (23) is mean square bounded if and only if where is the spectral radius and is the Kronecker product.

4. Main Results and Proofs

The main results are presented in this section. To start with, we first recall some lemmas that will be needed in the proof of our main results.

Lemma 9 (see [29]). Given matrices , , , and with compatible dimensions such that , let be a positive definite matrix and an arbitrary constant such that . Then, we have

Lemma 10 (see [30]). Given matrices and , there exists a solution to the matrix equation if and only if , where denotes the Moore-Penrose inverse of . Furthermore, all solutions can be parameterized by where is an arbitrary matrix.

Lemma 11 (see [31]). For a given negative definite matrix , there always exists a matrix such that .

Lemma 12 (Schur complement). Given constant matrices , , and , where and , then if and only if

Here, in the following theorem, a two-step approach will be developed to solve the filer problem. Firstly, we will characterize an upper bound on the steady-state error covariance satisfying (40) in terms of some free parameters, let this upper bound meet prespecified variance constraints (18), and then we parameterize all desired filter gains with which the resulting steady-state error covariance is not more than the obtained upper bound. It is shown that the design of variance-constraint robust filtering for uncertain systems with limited communication capacity is related to two quadratic matrix inequalities. For presentation convenience, we first define

Theorem 13. Let , and let be given positive scalar sequences and inequalities and hold. If there exist positive scalars , such that the following two quadratic matrix inequalities respectively, have positive-definite solutions , , and where denotes the Moore-Penrose inverse of , and in (50) Moreover, let be an arbitrary matrix satisfying (see Lemma 11) and an arbitrary orthogonal matrix (i.e., ). Then, filter (16) with the parameters determined by (51) and will be such that, for all admissible uncertainties caused by the limited communication capacity, (1) the state of augmented system (23) is mean square bounded and (2) the steady-state error covariance meets .

Proof. Define . Then, it follows directly from Lemma 9 and definitions (45)–(48) that where
Making use of Matrix Inverse Theory, we have and therefore, from inequality (49), we can obtain that .
Next, notice that the matrix is singular. It then follows from Lemma 10 that, there exists a solution such that if and only if (51) holds. Furthermore, if (51) is true, (52) gives a solution. Hence, substituting the expression of in (52) into (56) leads to easily.
Furthermore, we now consider , by using definitions (45)–(48), and we can rearrange (57) as follows:
Noticing the expression of in (53) and the fact that , we have
Hence, it follows from (59) and in view of the definition of the matrix and inequality (50), we can obtain that and the .
To this end, we can conclude that . Thus, it follows from (47) that which leads to (41). From Lemma 8, we know that the state of augmented system (23) is mean square bounded and there exists a symmetric positive semidefinite solution to (40). Then, the first claim of this theorem is proved.
Furthermore, subtract (40) from (61) to give which indicates again from Lemma 8 that , and therefore This completes the proof of the theorem.

Remark 14. From all of above discussions, we know that, if the conditions of Theorem 13 are all met and the positive-definite solutions satisfy then the design objective of robust filter for uncertain systems with limited communication capacity will be accomplished. It is mentionable that the existence of a positive-definite solution to (49) implies the asymptotical Schur stability of system matrix , and this means that Assumption 2 should hold.

Remark 15. It can be seen, from theory 1, that there exits much explicit freedom, such as the choice of parameter , the orthogonal matrix for our present design approach. We could use the freedom feature to take the more expected performance constraints into account within the same framework (e.g., the transient requirement and reliability behavior on the filtering process), which provide us with one possible future research direction (see Section 5).
Eventually, by using the Schur Lemma (Lemma 12), we present the main results which simultaneously presented a solution to matrix inequalities (49) and (50) as follows.

Corollary 16. If there exist positive scalars and two positive-definite matrices , such that the following LMIs (65) and (66) and the matrix inequality (50) hold and satisfies , then filter (16) determined by (52) and (53) will achieve the desired robust filtering performance for uncertain systems with limited communication capacity which is discussed previously. Cosider

We now briefly discuss the solvability of quadratic matrix inequalities (49) and (50), which play a key role in designing the expected filters. Since the parameter of (50) is not included in (49), from the previous corollary, for and , algebraic matrix inequality (49) is equal to the linear matrix inequalities (LMI) (65)-(66) which can be effectively solved by using LMI toolbox in MATLAB. Then, after and are gained and positive scalar sequences and are given satisfying that and , matrix inequality (50) becomes a standard Riccati-like matrix inequality for which can be easily solved in terms of the existed approach.

5. Numerical Example

In this section, we demonstrate the theory developed in this paper by means of a simple example. Consider the linear discrete-time system described by (1) with parameters given by

Furthermore, we assume the following.(1)The parameter for the quantizer is ; then, ; that is, .(2)The stochastic variable is Bernoulli distributed white sequences taking values 0 and 1 with and , respectively.(3)The constant and are given by and .(4)The process noise and measurement noise are zero-mean Gaussian white noise sequences with unity covariance.

The purpose of this example is to design the filter parameters and such that, for all admissible perturbations and multiple packet dropouts, augmented system (16) is mean square bounded, and the steady-state error covariance satisfies

We employ the standard LMI techniques to check the solvability of original matrix inequality (49) for and and solve the standard Riccati-like matrix for . Therefore, we obtain

From (52), one of the filter parameters, , is calculated as follows:

To obtain another parameter , we choose meeting and select the orthogonal matrix as follows: Then, we can get the parameter from (53) as follows:

From the above simulated results, it is not difficult to get that and ; hence, the specified mean square bounded is achieved. Furthermore, in Figure 1, the above curve is the bound and the under curve is the steady-state error variance; it is obvious that the steady-state error variance is not more than the individual upper bound, which verifies that steady-state error variance constraint is also achieved. Therefore, it is shown that the theory developed in this paper is effective and practical.

6. Conclusion

In this paper, we have considered a variance-constraint robust filtering problem for discrete-time systems with limited communication capacity. It has been shown that the problem is solvable if a couple of LMIs or Riccati-like matrix inequalities have positive definite solutions. A numerical example is presented to demonstrate the effectiveness and flexibility of the proposed design approach.

Further extension of the present results to more complex systems such as nonlinear systems is a possible future research direction. Moreover, in NCSs, there usually exist simultaneously the parameter uncertainties and intermediate uncertainties; thus, simultaneously considering delay, quantization, observation losses, and parameter uncertaintes is another future avenue of research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by Liaoning Province Education Department Scientific Research Project (L2012445), the National Natural Science Foundation of China under Grant no. 60504017, Fok Ying Tung Education Foundation under Grant 111066, and Program for New Century Excellent Talents in University under Grant NCET-04-0982.

References

M. Moayedi, Y. K. Foo, and Y. C. Soh, “Adaptive Kalman filtering in networked systems with random sensor delays, multiple packet dropouts and missing measurements,” IEEE Transactions on Signal Processing, vol. 58, no. 3, pp. 1577–1588, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
W.-A. Zhang, L. Yu, and G. Feng, “Optimal linear estimation for networked systems with communication constraints,” Automatica, vol. 47, no. 9, pp. 1992–2000, 2011.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
H. J. Gao and T. W. Chen, “ $H_{\infty}$ estimation for uncertain systems with limited communication capacity,” IEEE Transactions on Automatic Control, vol. 52, no. 11, pp. 2070–2084, 2007.
View at: Publisher Site | Google Scholar
H. L. Dong, Z. D. Wang, and H. J. Gao, “Robust H_∞ filtering for a class of nonlinear networked systems with multiple stochastic communication delays and packet dropouts,” IEEE Transactions on Signal Processing, vol. 58, no. 4, pp. 1957–1966, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
S. Sun, “Linear minimum variance estimators for systems with bounded random measurement delays and packet dropouts,” Signal Processing, vol. 89, no. 7, pp. 1457–1466, 2009.
View at: Publisher Site | Google Scholar | Zentralblatt MATH
Z. D. Wang, F. W. Yang, D. W. C. Ho, and X. Liu, “Robust finite-horizon filtering for stochastic systems with missing measurements,” IEEE Signal Processing Letters, vol. 12, no. 6, pp. 437–440, 2005.
View at: Publisher Site | Google Scholar
Z. D. Wang, D. W. C. Ho, and X. H. 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 Site | Google Scholar | MathSciNet
S. C. Smith and P. Seiler, “Estimation with lossy measurements: jump estimators for jump systems,” IEEE Transactions on Automatic Control, vol. 48, no. 12, pp. 2163–2171, 2003.
View at: Publisher Site | Google Scholar | MathSciNet
K. Gu, V. Kharitonov, and J. Chen, Stability of Time-Delay Systems, Birkhauser, 1st edition, 2003.
View at: MathSciNet
L. Schenato, “Optimal estimation in networked control systems subject to random delay and packet loss,” in Proceedings of the 45th IEEE Conference on Decision & Control (CDC '06), pp. 5615–5620, 2006.
View at: Google Scholar
M. Fu and L. Xie, “The sector bound approach to quantized feedback control,” IEEE Transactions on Automatic Control, vol. 50, no. 11, pp. 1698–1711, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
V. Malyavej and A. V. Savkin, “The problem of optimal robust Kalman state estimation via limited capacity digital communication channels,” Systems & Control Letters, vol. 54, no. 3, pp. 283–292, 2005.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
J. Liang, Z. Wang, and X. Liu, “Distributed state estimation for discrete-time sensor networks with randomly varying nonlinearities and missing measurements,” IEEE Transactions on Neural Networks, vol. 22, no. 3, pp. 486–496, 2011.
View at: Publisher Site | Google Scholar
D. R. Ding, Z. D. Wang, H. Dong, and H. Shu, “Distributed H_∞ state estimation with stochastic parameters and nonlinearities through sensor networks: the finite-horizon case,” Automatica, vol. 48, no. 8, pp. 1575–1585, 2012.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
B. Wang and G. Guo, “Robust estimation for discrete time-varying systems with limited communication capacity,” Asian Journal of Control, vol. 14, no. 2, pp. 502–511, 2012.
View at: Publisher Site | Google Scholar | MathSciNet
F. W. Yang, Z. D. Wang, and Y. S. Hung, “Robust Kalman filtering for discrete time-varying uncertain systems with multiplicative noises,” IEEE Transactions on Automatic Control, vol. 47, no. 7, pp. 1179–1183, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
Z. Dong and Z. You, “Finite-horizon robust Kalman filtering for uncertain discrete time-varying systems with uncertain-covariance white noises,” IEEE Signal Processing Letters, vol. 13, no. 8, pp. 493–496, 2006.
View at: Publisher Site | Google Scholar
F. Wang and V. Balakrishnan, “Robust steady-state filtering for systems with deterministic and stochastic uncertainties,” IEEE Transactions on Signal Processing, vol. 51, no. 10, pp. 2550–2558, 2003.
View at: Publisher Site | Google Scholar
E. Yaz and R. E. Skelton, “Continuous and discrete state estimation with error covariance assignment,” in Proceedings of the IEEE Conference of Decision Control, pp. 3091–3092, Brighton, UK, 1991.
View at: Google Scholar
Z. D. Wang, B. Huang, and P. Huo, “Sampled-data filtering with error covariance assignment,” IEEE Transactions on Signal Processing, vol. 49, no. 3, pp. 666–670, 2001.
View at: Publisher Site | Google Scholar
Z. D. Wang and H. Qiao, “Robust filtering for bilinear uncertain stochastic discrete-time systems,” IEEE Transactions on Signal Processing, vol. 50, no. 3, pp. 560–567, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
Z. D. Wang, D. W. C. Ho, and X. H. Liu, “Robust filtering under randomly varying sensor delay with variance constraints,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 51, no. 6, pp. 320–326, 2004.
View at: Publisher Site | Google Scholar
N. Elia and S. K. Mitter, “Stabilization of linear systems with limited information,” IEEE Transactions on Automatic Control, vol. 46, no. 9, pp. 1384–1400, 2001.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Fu and L. Xie, “The sector bound approach to quantized feedback control,” IEEE Transactions on Automatic Control, vol. 50, no. 11, pp. 1698–1711, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
E. Yaz and A. Ray, “Linear unbiased state estimation for random models with sensor delay,” in Proceedings of the 35th IEEE Conference on Decision and Control, pp. 47–52, Kobe, Japan, 1996.
View at: Google Scholar
S. L. Sun, “Linear minimum variance estimators for systems with bounded random measurement delays and packet dropouts,” Signal Processing, vol. 89, no. 7, pp. 1457–1466, 2009.
View at: Publisher Site | Google Scholar | Zentralblatt MATH
R. F. Stengel, Stochastic Optimal Control: Theory and Application, John Wiley & Sons, New York, NY, USA, 1986.
W. L. de Koning, “Optimal estimation of linear discrete-time systems with stochastic parameters,” Automatica, vol. 20, no. 1, pp. 113–115, 1984.
View at: Publisher Site | Google Scholar | Zentralblatt MATH
L. Xie, Y. C. Soh, and C. E. de Souza, “Robust Kalman filtering for uncertain discrete-time systems,” IEEE Transactions on Automatic Control, vol. 39, no. 6, pp. 1310–1314, 1994.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
A. Ben-Israel and T. N. E. Greville, Generalized Inverses: Theory and Applications, John Wiley & Sons, New York, NY, USA, 1974.
Z. Wang and B. Huang, “Robust $H_{2} / H_{\infty}$ filtering for linear systems with error variance constraints,” IEEE Transactions on Signal Processing, vol. 48, no. 8, pp. 2463–2467, 2000.
View at: Publisher Site | Google Scholar | MathSciNet

Copyright

Copyright © 2014 Baofeng 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.

PDF Download Citation

Download other formats

Order printed copies

Views

808

Downloads

963

Citations