A Novel Strong Tracking Fault Prognosis Algorithm

Zhang, Qi; Jiang, Wei; Li, Tian-Mei; Zheng, Jian-Fei

doi:https://doi.org/10.1155/2015/676289

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusions Acknowledgments References Copyright Related Articles

Special Issue

Modeling for Prognostics and Health Management: Methods and Applications

View this Special Issue

Research Article | Open Access

Volume 2015 | Article ID 676289 | https://doi.org/10.1155/2015/676289

A Novel Strong Tracking Fault Prognosis Algorithm

Qi Zhang,¹Wei Jiang,¹Tian-Mei Li,¹and Jian-Fei Zheng¹

Academic Editor: Gang Li

Received27 Nov 2014

Revised21 Dec 2014

Accepted21 Dec 2014

Published08 Jun 2015

Abstract

Improving the ability to track abruptly changing states and resolving the degeneracy are two difficult problems to particle filter applied to fault prognosis. In this paper, a novel strong tracking fault prognosis algorithm is proposed to settle the above problems. In the proposed algorithm, the artificial immunity algorithm is first introduced to resolve the degeneracy problem, and then the strong tracking filter is introduced to enhance the ability to track abruptly changing states. The particles are updated by strong tracking filter, and better particles are selected by utilizing the artificial immune algorithm to estimate states. As a result, the degeneracy problem is resolved and the accuracy of the proposed fault prognosis algorithm is improved accordingly. The feasibility and validity of the proposed algorithm are demonstrated by the simulation results of the standard validation model and the DTS200 system.

1. Introduction

Particle filter (PF) is a leading and powerful algorithm for estimating the states of nonlinear or non-Gaussian systems. The past decades have witnessed a wide range of applications, including target tracking [1–4], data detection [5], modeling [6, 7], price forecasting, and fault detection [8–11]. On the other hand, a great number of investigators are interested in enriching particle filtering framework, and many new particle filters are proposed in recent years [12–15]. In these studies, it is found that resolving the degeneracy problem and improving the ability to track abruptly changing states are two difficult problems to particle filter applied to fault prognosis [16–20]. The degeneracy problem means that most particles are assigned to zero weights. As a result, the performance of the particle filter deteriorates because most computational resource is wasted. It is noted that, however, degeneracy can be reduced by resampling or choosing good importance sampling functions. Along this line of research, many resampling algorithms have been proposed for reducing the degeneracy. In these resampling algorithms, sequential importance resampling (SIR) is the representation which largely copies the particles with larger weights to replace the particles with smaller weights [21, 22]. The degeneracy problem is partially addressed, but the sample impoverishment is not fully concerned. Sample impoverishment means that most particles are the same in the set of particles since the particles with larger weights are largely copied. In this circumstance, choosing good importance sampling functions deserves further studies, and many investigators are interested in this question. For example, extended Kalman filter (EKF) was introduced to propose extended particle filter (EPF) by De Freitas et al. [6], and unscented Kalman filter (UKF) was introduced to propose unscented particle filter (UPF) by van der Merwe et al. [23]. Both EPF and UPF can resolve the degeneracy problem, but they cannot track abruptly changing states due to the disadvantages of EKF and UKF.

Strong tracking filter (STF) has good performance for tracking abruptly changing states, and thus it can be used to update particles. In addition, it is known that artificial immunity (AI) can search for the best one from all the range, and thus it can be used to clone and vary particles. Therefore, in this paper STF and AI algorithms are utilized jointly to improve particle filter algorithm. As a result, a novel fault prognosis algorithm based on strong tracking artificial immunity particle filter (STAIPF) is proposed to settle the above discussed problems. In the proposed algorithm, the artificial immunity algorithm is first introduced to resolve the degeneracy problem, and then the strong tracking filter is introduced to enhance the ability to track abruptly changing states. More specifically, the particles are updated by strong tracking filter, and better particles for states estimation are selected by utilizing the artificial immune algorithm to enhance the diversity of samples. Therefore, the degeneracy problem and sample impoverishment are resolved simultaneously, and the accuracy of the proposed fault prognosis algorithm is improved as well. Finally, the feasibility and validity of the proposed algorithm are demonstrated by the simulation results of the standard validation model and the DTS200 system.

The remainder of this paper is structured as follows. In Section 2, the particle filter is introduced. Section 3 provides the artificial immune algorithm. In Section 4, we present a strong tracking filter. A novel strong tracking fault prognosis algorithm is proposed in Section 5. Section 6 provides simulation results. This paper is concluded in Section 7.

2. Particle Filter

Actually, particle filter is a sequential Monte Carlo methodology. Its primary principle is to recursively compute relevant probability distributions by importance sampling and to approximate the probability distributions with discrete random variables. Detailed information of particle filter can be found in [21].

In general, the following state-space and observation equations are considered:where the subscript denotes time index, is an observations vector, is a state vector, is a system noise vector, is an observation noise vector, is a measurement function, and is a system transition function. The first equation is known as the state equation while the second one is known as the measurement equation. It is usually assumed that the analytical forms of the two functions and the distributions of the two noises in (1) are known. Then, the object is to recursively estimate based on the observations and the above assumptions.

Due to the general nature of model (1) and the impact of non-Gaussian noises, it is difficult to solve the above filtering problem in an analytical manner. In this case, particle filter is an effective alternative. In particle filter, the number of effective particles is commonly denoted by , which is used to weigh the degeneracy degree of the particles. In other words, the smaller implies the worse degeneracy degree. Here, is defined aswhere denotes the normalized weighting coefficient of particle at time , and denotes rounding to the nearest integer.

Based on the above descriptions, the SIR algorithm can be summarized as follows [21], which is one of the most popular particle filtering algorithms.

A Procedure Description of SIR

Step 1 (initialization). Draw particles according to initial importance function, and is assumed to denote the th particle, and then set .

Step 2. Update.where denotes the particle at time , which is updated according to formula (3).

Step 3. Assign the weighting coefficient of according to

Step 4. Normalize weighting coefficient by

Step 5. If , then run system resampling and assign the same weighting coefficient to all particles:End.

Step 6. Estimate the state at time according to

Step 7. Return to Step 2.

3. Artificial Immune Algorithm

Artificial immune algorithm (AIA), heredity algorithm, and evolutionary algorithm are all bionic algorithms which simulate the behaviors of natural organisms. Artificial immune algorithm is characterized by its diversity, self-regulation, clone, and mutation, which is a global searching algorithm based on natural immune systems [24, 25].

Antigen and antibody correspond to target function and possible value of the associated optimization problem, respectively. Generally, affinity has two forms. One is the appetency which denotes the matching degree between antigen and antibody. The other is the repellency which denotes the similar degree of antigen and antibody, and the repellency ensures the diversity of the antibodies. Immune algorithm is based on the memory cell and clone number calculated according to affinity, which can converge to the optimal solution ensured by the above steps. The basic procedure of immune algorithm is summarized as Figure 1.

At first, we should identify antigen and draw the initial antibody, then calculate the affinity between the antibodies and antigen, and update the memory cell. Next, clone and mutate the better antibodies to produce the new antibody set. Finally, return to calculate the affinity until the end.

4. Strong Tracking Filter

The usual EKF cannot track abruptly changing states when the system and the filter become steady, which may result in greater tracking error and even divergence. It is known that strong tracking filter has the advantage to track abruptly changing states. In this paper, the adopted strong tracking filter is a suboptimal fading extended Kalman filter (SFEKF) [26, 27].

Consider a kind of systems as follows:where the integer is a discrete time variable, is a state vector, is an input vector, and is an output vector. and are nonlinear functions which have first-order continuous partial derivatives with respect to the states. is a known matrix. is a q-dimension Gaussian white noise vector, and is an m-dimension Gaussian white noise vector. In this paper, we consider the case that the noises have the following statistical characteristics:where is a symmetric nonnegative definite matrix, is a symmetric positively definite matrix, and is the initial state. It is further assumed that is statistically independent of and .

For the considered system represented by (8), a procedure description of strong tracking filter is summarized as follows [26, 27].

A Procedure Description of Strong Tracking Filter

Step 1. Set , and choose the initial value of , and infirmness gene .

Step 2. Calculate according to

Step 3. Calculate and according to

Step 4. Calculate according to

Step 5. Calculate the fading factor aswhere the parameter is named forgetting factor, and . In general, we set .

Step 6. Estimate the value of by the following recursive formulas:

Step 7. Calculate by, andreturn to Step 2 until the end of circle.

5. Strong Tracking Fault Prognosis Algorithm

A novel strong tracking fault prognosis algorithm is proposed in this section, in which the artificial immunity algorithm is introduced to resolve the degeneracy problem, and the strong tracking filter is introduced to enhance the ability to track abruptly changing states. To be specific, the particles are updated by strong tracking filter, and better particles for estimation are selected by artificial immune algorithm to enhance the diversity of samples. As such, the proposed algorithm can predict fault more accurately.

The following is the detailed idea description of the proposed algorithm. First, draw particles as the initial antibodies according to the important sampling function. Then, we update the particles by STF and calculate the appetency and repellency of the antibodies to ascertain the clone number of every antibody. A new antibody set is made up of the initial antibodies and clone antibodies, in which some similar antibodies are discarded by respectively calculating affinities of all the antibodies. By this procedure, the good ones are chosen to estimate the state of next time. Thus, the degeneracy problem and the sample impoverishment are resolved, and the ability of tracking abruptly changing states is much improved because of the introduced STF in the proposed fault prognosis algorithm.

In this paper, we consider the nonlinear dynamic system as (1), and set 1 and as the antigen and antibody, respectively. Based on the STF algorithm and AI algorithm, we propose a novel strong tracking fault prognosis algorithm, which is illustrated step by step in the following.

A Procedure Description of the Proposed Fault Prognosis Algorithm

Step 1 (initialization). Draw particles according to the importance function at time , and is assumed to denote each particle; then set and choose the initial value of , , and infirmness gene .

Step 2 (update). Update the particles using the strong tracking filter, which calculates according to (10)–(15).

Step 3 (assign weights). Assign the weighting coefficients to according to (4) and then normalize weighting coefficients according to (5).

Step 4 (build initial antibodies). The updated particles are regarded as initial antibodies put in the memory cell.

Step 5 (calculate the affinity). Calculate and according to the following formulas which are the appetencies of every antibody and the repellencies between the antibodies.

Step 6 (clone the antibodies). Clone the antibodies according to the appetencies of every antibody. Here the clone number is ascertained by

Step 7 (mutation). Mutate the whole antibodies according to

Step 8 (choose the optimal particles). Calculate and of the mutated antibodies and discard one of the two antibodies if which is set according to the value range of particles. Sort the remainder antibodies in terms of and choose the anterior antibodies to replace the memory cell. Return to Step 3 until the ending criterion is satisfied.

Step 9 (estimate). Consider the antibodies in memory cell as new particles to estimate the states at according to (7).

Step 10. -step-ahead prediction.

Step 11. -step-ahead fault prognosis probability at time.where denotes the fault state.

Step 12. Calculate the weighted fault prognosis probability at time.

Step 13. If , resample the particles.

Step 14. Set , return to Step 2.

6. Simulation Experiments and Results Analysis

In this section, two simulation experiments are presented to demonstrate the effectiveness of the proposed algorithm. The first one is a standard validation model that has been widely used to validate the ability to track the states and the ability to resolve the degeneracy problem and sample impoverishment. The second experiment is used to validate whether or not the proposed algorithm can predict the system’s fault with abruptly changing states.

6.1. Experiment 1 of the Standard Validation Model

The model used in this section can be found in many studies [28], which is a standard validation model to validate the performance of various particle filters.

Specifically, the state equation and the measurement equation of the standard validation model are represented aswhereand is a Gaussian noise with mean 0 and variance 10, and is also a Gaussian noise with mean 0 and variance 1. The normal distribution is the initial importance function. For implementing the proposed algorithm, we set , , , and , where is the number of particles. In this simulation, we choose cycling 5 times as the terminal criterion. The estimated results of strong tracking artificial immunity particle filter for the standard validation model are shown in Figure 2.

(a)

(b)

(c)

(d)

Figure 2 shows that the estimate values of states are almost the same with the real values of states, so STAIPF can track states of the standard validation model and the proposed algorithm is feasible and effective.

The number of effective particles is used to weight the degeneracy degree of the particles. The smaller denotes more serious degeneracy degree. In this simulation, STAIPF and STPF run 10 times and the means of are respectively calculated to validate whether STAIPF can resolve the degeneracy problem or not to a certain extent. Table 1 is the simulation results of effective particles, in which and , respectively, denote the number of effective particles and the mean of STAIPF, and and , respectively, denote the number of effective particles and the mean of STPF.

Table 1 and the mean of of two algorithms show that is 2.4 times as much as . Therefore, STAIPF is better than STPF about the ability to resolve the degeneracy problem. As a result, STAIPF can resolve degeneracy problem because of introducing artificial immunity algorithm.

Further, the normalized weighting coefficients of particles at the final time are shown in Figure 3 to demonstrate that STAIPF can resolve the degeneracy problem better than STPF.

(a)

(b)

We can find that the normalized weighting coefficients of particles at the final time about STAIPF are all greater than zero, but most of the normalized weighting coefficients about STPF equals zero. As a result, STAIPF can resolve the degeneracy problem better than STPF.

6.2. Experiment 2 of the DTS200 System

6.2.1. System Description

A modified mathematical model of the DTS200 system [26] as shown in Figure 4 was produced by the Amira Automation Company in Germany. DTS200 system is composed of three cylinders tank 1 (), tank 2 (), and tank 3 () with the cross-section area . The cylinders are connected to each other by pipes with the cross-section area . The outflow valve with circular cross-section area is located at tank 2. The outflowing distilled water collected in a water tank supplies to the pump1 and pump2. A leak valve with the cross-section area is fixed on the bottom of each tank, which is usually closed.

The pump flows and denote the input signals, and the liquid heights , , and of , , and denote the output signals. The liquid heights and measured by 3 piezoresistive difference pressure sensors are controlled by and .

The system is modeled by

With the vectors defined by

The significance of the above mathematical symbols is summarized as follows: : outflow coefficients, : liquid heights [m], : flow rates [m³/s], , : pump flow rates [m³/s], : cross-section area of cylinder [m²], : cross-section area of connection pipe [m²],where , and . Considerwhere : earth’s acceleration of gravity [m/s²], : sign of the argument , : liquid height difference between two tanks connected to each other, : flow in the connecting pipe, and : outflow coefficient.

The results of the unknown quantities are calculated as

The technical parameters of the DTS200 system are given as follows:

For the above concerned system, it can be equivalently modeled with Euler discretization technique as the following nonlinear discrete model:where = 1 s is the sampling interval, the initial levels of , , and are m, m, m, respectively, and = = 4.5 × 10⁻⁵m³/s.

To achieve fault prognosis, in this experimental study, the fault state is defined as

6.2.2. Simulation Results and Analysis

The abrupt change fault is introduced to the system in order to verify the tracking ability of the proposed algorithm. We suppose that changes according to the following formula:

To implement the proposed algorithm, we set the number of particles , , , , and the process noise and measurement noiseare and , respectively. Further, we suppose that the simulation time is and the initial particle meets the condition as .

The simulation results of state estimation and fault prognosis probability of the normal system are shown as Figure 5, which are performed by the strong tracking fault prognosis algorithm.

(a)

(b)

(c)

(d)

We can find from Figure 5 that the estimated values of the system states are all around the real ones with time and the weighted fault prognosis probability is always zero, and the normalized weighting coefficients of most of particles are not zero. As a result, the proposed algorithm can accurately predict the normal system and resolve the degeneracy problem.

The simulation results of the state estimation and fault prognosis probability of the faulty system are shown as Figure 6, which are performed by the strong tracking fault prognosis algorithm.

(a)

(b)

(c)

(d)

In this case, it is found from Figure 6 that the system states deviate from the real value when changes with time, and the range ability of is the biggest. The system fault will appear when deviates from the real value to a certain point. The normalized weighting coefficients of most of particles are not zero. As a result, the proposed algorithm can accurately predict the system fault in time and resolve the degeneracy problem well.

Further, we extract the values of and fault prognosis probability from Figures 6(a) and 6(b). In this paper, the fault is declared when the fault prognosis probability prob is greater than 0.5. Specifically, the values of and the weighted fault prognosis probability prob during the period are shown in Table 2 to demonstrate that the proposed algorithm can accurately predict the system fault.

According to the definition of the fault state in (31), the system fault will happen when deviates from the real value to 0.2700 at time 73. It is observed from Table 2 that the proposed algorithm can predict the fault at time 71. As a result, the algorithm can accurately predict the system fault in time.

Together with the above simulation results, we can observe that the degeneracy problem and the sample impoverishment can be resolved simultaneously in the proposed fault prognosis algorithm, and thus the feasibility and validity of the proposed algorithm are demonstrated.

7. Conclusions

Considering that SIR and EPF cannot track abruptly changing states and STPF can track abruptly changing states, a novel strong tracking fault prognosis algorithm is proposed to settle the degeneracy problem and the sample impoverishment in particle filter, which is applied to predict the fault of systems with abruptly changing states. The simulation results of the standard validation model show that the strong tracking fault prognosis algorithm can resolve degeneracy problem, and the simulation results of the DTS200 system demonstrate that the proposed algorithm can track abruptly changing states and accurately predict the system fault in time. As a result, the proposed algorithm is feasible and effective for the applications in fault prognosis.

Conflict of Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant 61104223, 61174030, 61304103, 61374126, and 61473094 and the China’s Postdoctoral Science Fund 2014M552670.

References

W. Ng, J. Li, S. Godsill, and J. Vermaak, “A review of recent results in multiple target tracking,” in Proceedings of the 4th International Symposium on Image and Signal Processing and Analysis, pp. 40–45, September 2005.
View at: Google Scholar
J. A. Pocock, S. L. Dance, and A. S. Lawless, “State estimation using the particle filter with mode tracking,” Computers and Fluids, vol. 46, no. 1, pp. 392–397, 2011.
View at: Publisher Site | Google Scholar
A. Xydes, M. Moline, C. G. Lowe, T. J. Farrugia, and C. Clark, “Behavioral characterization and particle filter localization to improve temporal resolution and accuracy while tracking acoustically tagged fishes,” Ocean Engineering, vol. 61, no. 15, pp. 1–11, 2013.
View at: Publisher Site | Google Scholar
C. Chang and R. Ansari, “Kernel particle filter for visual tracking,” IEEE Signal Processing Letters, vol. 12, no. 3, pp. 242–245, 2005.
View at: Publisher Site | Google Scholar
T. Ghirmai, M. F. Bugallo, J. Míguez, and P. M. Djurić, “A sequential Monte Carlo method for adaptive blind timing estimation and data detection,” IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 2855–2865, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
J. F. G. De Freitas, M. Niranjan, A. H. Gee, and A. Doucet, “Sequential Monte Carlo methods to train neural network models,” Neural Computation, vol. 12, no. 4, pp. 955–993, 2000.
View at: Publisher Site | Google Scholar
D. Ryu, F. Liang, and B. K. Mallick, “Sea surface temperature modeling using radial basis function networks with a dynamically weighted particle filter,” The Journal of the American Statistical Association, vol. 108, no. 501, pp. 111–123, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
X. R. Li, C. W. Yu, S. Y. Ren, C. H. Chiu, and K. Meng, “Day-ahead electricity price forecasting based on panel cointegration and particle filter,” Electric Power Systems Research, vol. 95, no. 2, pp. 66–76, 2013.
View at: Publisher Site | Google Scholar
T. Wei, Y. F. Huang, and P. Chen, “Particle filtering for adaptive sensor fault detection and identification,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 3807–3812, Orlando, Fla, USA, May 2006.
View at: Publisher Site | Google Scholar
M. Z. Chen and D. H. Zhou, “Particle filter based fault prediction of nonlinear systems,” in IFAC Symposium Proceedings of Safe Process, pp. 2971–2977, Elsevier Science, Washington, DC, USA, 2001.
View at: Google Scholar
F. Cadini, E. Zio, and G. Peloni, “Particle filtering for the detection of fault onset time in hybrid dynamic systems with autonomous transitions,” IEEE Transactions on Reliability, vol. 61, no. 1, pp. 130–139, 2012.
View at: Publisher Site | Google Scholar
M. Ades and P. J. van Leeuwen, “An exploration of the equivalent weights particle filter,” Quarterly Journal of the Royal Meteorological Society, vol. 139, no. 672, pp. 820–840, 2013.
View at: Publisher Site | Google Scholar
P. M. Stano, Z. Lendek, and R. Babuska, “Saturated particle filter: almost sure convergence and improved resampling,” Automatica, vol. 49, no. 1, pp. 147–159, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
W. Hassan, N. Bangalore, P. Birch, R. Young, and C. Chatwin, “An adaptive sample count particle filter,” Computer Vision and Image Understanding, vol. 116, no. 12, pp. 1208–1222, 2012.
View at: Publisher Site | Google Scholar
S.-H. Park, Y.-J. Kim, and M.-T. Lim, “Novel adaptive particle filter using adjusted variance and its application,” International Journal of Control, Automation and Systems, vol. 8, no. 4, pp. 801–807, 2010.
View at: Publisher Site | Google Scholar
X. S. Si, C. H. Hu, X. Y. Kong, and D. H. Zhou, “A residual storage life prediction approach for systems with operation state switches,” IEEE Transactions on Industrial Electronics, vol. 61, no. 11, pp. 6304–6315, 2014.
View at: Publisher Site | Google Scholar
G. Li, S. J. Qin, Y. Ji, and D. Zhou, “Reconstruction based fault prognosis for continuous processes,” Control Engineering Practice, vol. 18, no. 10, pp. 1211–1219, 2010.
View at: Publisher Site | Google Scholar
X.-S. Si, M.-Y. Chen, W. Wang, C.-H. Hu, and D.-H. Zhou, “Specifying measurement errors for required lifetime estimation performance,” European Journal of Operational Research, vol. 231, no. 3, pp. 631–644, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
J. Zhang, C. Hu, X. Si et al., “Remaining useful life estimation for systems with time-varying mean and variance of degradation processes,” Quality and Reliability Engineering International, vol. 30, no. 6, pp. 829–841, 2014.
View at: Google Scholar
X.-S. Si, W. Wang, C.-H. Hu, M.-Y. Chen, and D.-H. Zhou, “A Wiener-process-based degradation model with a recursive filter algorithm for remaining useful life estimation,” Mechanical Systems and Signal Processing, vol. 35, no. 1-2, pp. 219–237, 2013.
View at: Publisher Site | Google Scholar
P. M. Djurić, J. H. Kotecha, J. Zhang et al., “Particle filtering,” IEEE Signal Processing Magazine, vol. 20, no. 5, pp. 19–38, 2003.
View at: Publisher Site | Google Scholar
M. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, “A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking,” IEEE Transactions on Signal Processing, vol. 50, no. 2, pp. 174–188, 2002.
View at: Publisher Site | Google Scholar
R. van der Merwe, A. Doucet, N. de Freitas, and E. Wan, “The unscented particle filter,” Tech. Rep. CUED/F-INFENG/TR 380, Cambridge University, London, UK, 2000.
View at: Google Scholar
D. Dasgupta, Artificial Immune Systems and Their Applications, Springer, Berlin, Germany, 1999.
View at: Publisher Site
X. W. Zhou, “Evolutionary algorithm and its application in artificial immune system,” in Proceedings of the 2nd International Symposium on Intelligent Information Technology Application (IITA '08), pp. 32–36, Shanghai, China, December 2008.
View at: Publisher Site | Google Scholar
D. H. Zhou and Y. Z. Ye, Modern Fault Diagnosis and Fault Tolerant Control, Tsinghua University press, Beijing, China, 2000.
W.-H. Wang, Z. Peng, and D.-H. Zhou, “State estimation of a class of hybrid systems in the presence of unknown mode transition,” Acta Automatica Sinica, vol. 31, no. 3, pp. 456–457, 2005.
View at: Google Scholar | MathSciNet
Y. W. Mo and D. Y. Xiao, “Evolutionary particle filter and its application,” Control Theory and Applications, vol. 22, no. 2, pp. 269–270, 2005.
View at: Google Scholar

Copyright

Copyright © 2015 Qi Zhang 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

687

Downloads

700

Citations