Research Article | Open Access
Germán A. Ramos, Josep M. Olm, Ramon Costa-Castelló, "Digital Repetitive Control under Nonuniform Sampling: An LMI Stability Analysis", Mathematical Problems in Engineering, vol. 2011, Article ID 120570, 16 pages, 2011. https://doi.org/10.1155/2011/120570
Digital Repetitive Control under Nonuniform Sampling: An LMI Stability Analysis
Digital repetitive control is a technique which allows tracking periodic references and/or rejecting periodic disturbances. Repetitive controllers are usually designed assuming a fixed fundamental frequency for the signals to be tracked/rejected and its main drawback being a dramatic performance decay when this frequency varies. A usual approach to overcome the problem consists of an adaptive change of the sampling period according to the reference/disturbance period variation. This paper presents a stability analysis of a digital repetitive controller working under time-varying sampling period by means of an LMI gridding approach. Theoretical developments are illustrated with experimental results, which are preceded by a detailed description of fundamental issues related to the implementation procedure.
Repetitive control [1, 2] is an internal model principle-based control technique  that allows both the tracking and rejection of periodic signals. Essentially, this is achieved including a generator of the reference/disturbance in the control loop. Its use has reported successful results in different control areas, such as CD and disk arm actuators , robotics , electro-hydraulics , electronic rectifiers , pulse-width modulated inverters [8, 9], and shunt active power filters .
It is usual to design repetitive controllers assuming a fixed period for the signals to be tracked/rejected. Then, a fixed sampling period is selected and, eventually, the value of the ratio is structurally embedded in the control algorithm. However, it is well known that even slight changes in the frequency of the tracked/rejected signals result in a dramatic decay of the controller performance .
A first set of proposals dealing with this problem maintain the initially selected sampling period, and robustness is achieved by means of large memory elements [11, 12] or introducing a fictitious sampler operating at a variable sampling rate and later using a fixed frequency internal model . Both ideas work well for small frequency variations at the cost of increasing the computational burden.
The second approach proposes to adapt the controller sampling rate according to the reference/disturbance period [14–16]. This allows preserving the steady-state performance while maintaining a low computational cost, but, on the other hand, the original linear time invariant (LTI) system becomes linear time varying (LTV). This structural change requires a new stability study, but no formal proofs regarding this issue are reported in the quoted references.
This paper analyzes the stability of a system containing a digital repetitive controller working under time-varying sampling period. The proposed methodology uses a linear matrix inequality (LMI) gridding approach [17, 18] that allows assessing the bounded input-bounded output (BIBO) stability of the closed-loop system in a known, bounded interval where the reference/disturbance period is assumed to vary. The theoretical results are experimentally validated through a mechatronic plant . A detailed description and discussion about fundamental issues related to the implementation procedure is also provided.
The structure of the paper is as follows. Section 2 contains a brief introduction to digital repetitive control and a study of stability issues in case of constant sampling period. Section 3 analyzes the stability of the system under a time-varying sampling period using LMI techniques. Experimental results are collected in Section 4, while conclusions and further research lines are presented in Section 5.
2. Digital Repetitive Control under Constant Sampling Period
Recall that the trigonometric Fourier series expansion of a -periodic signal reads as By the Internal Model Principle , the inclusion of the generator of signal (2.1) in the control loop results in the tracking/rejection of any -periodic reference/disturbance signal. Hence, following , the transfer function of a generator of (2.1) may be written as However, being a delay term with a gain , it is sufficient to include inside the control loop, which can be implemented as with a positive feedback, as depicted in Figure 1. Notice that the transfer function (2.2) has poles at , . Therefore, from a frequency point of view, exhibits infinite gain at frequencies , for all . This assures zero-tracking error at these frequencies in closed loop if the closed-loop system is stable.
It is also worth mentioning that some studies relate repetitive control to control learning techniques (see, e.g., ). This is due to the fact that the basic repetitive structure learns a signal of length and repeats it as a periodic signal of period if the input to the system is set to zero (see Figure 1).
The implementation of a time delay in continuous time is a complicated point. Fortunately, in discrete time it is an easier task; if the reference/disturbance signal period is a multiple of the sampling period , the digital implementation is reduced to a circular queue. Therefore, the discrete transfer function that should be included in the loop is where .
In addition to the constraint that represents the demand of a constant ratio between and , it is important to point out that should be selected taking into account that discrete-time implementations can only deal with those harmonics which are below the Nyquist frequency .
Repetitive controllers are composed of two main elements: the internal model, , and the stabilizing controller, . The internal model is the one in charge of guaranteeing null or small error in steady state, while the stabilizing controller assures closed-loop stability. Several types of internal models are used depending on the concrete periodic signal to deal with [19–22]. In this work the generic internal model is used, that is, is indeed taken as plays the role of a low-pass filter in charge of introducing robustness in the high frequency range . Although the internal model and the stabilizing controller can be arranged in different ways, most repetitive controllers are usually implemented in a “plug-in” fashion , as depicted in Figure 2: the repetitive compensator is used to augment an existing nominal controller, . This nominal compensator is designed to stabilize the plant, , and provides disturbance attenuation across a broad frequency spectrum. The signals denote, respectively, the -transforms of the input signal to be tracked and of the output disturbance to be rejected.
Assume that both and are constant, which makes also constant, and let stand for the corresponding -transform of . Sufficient stability criteria are given in the next Proposition.
Proposition 2.1. The closed-loop system of Figure 2 is stable if the following conditions are fulfilled [19, 20].(1)The closed-loop system without the repetitive controller, , is stable, where (2).(3), where is a design filter to be chosen.
Remark 2.2. These conditions hold for a proper design of , , and . Namely [19, 20],(i)it is advisable to design the controller with a high enough robustness margin; (ii) is designed to have gain close to 1 in the desired bandwidth and attenuate the gain out of it; (iii)a trivial structure which is often used for in case that is minimum phase is : Otherwise, alternative techniques should be applied in order to avoid closed-RHS plane zero-pole cancellations . Moreover, there is no problem with the improperness of because the internal model provides the repetitive controller with a high positive relative degree. Finally, as argued in , must be designed looking for a trade-off between robustness and transient response.
3. Digital Repetitive Control under Variable Sampling Period
The repetitive controller introduced in the previous section contains the ratio , which is embedded in the controller implementation. This is not a problem if the reference or disturbance periodic signal has a known constant period. However, the controller performance decays dramatically when a variation of appears . This paper propounds to adapt the controller sampling period following the reference/disturbance period , with the aim of maintaining a constant value for . Hence, on the one hand, , and are designed and implemented to provide closed-loop stability for a nominal sampling period , in accordance with Proposition 2.1 and Remark 2.2. Their structure remains always invariant; that is, it undergoes no further structural changes. On the other hand, the period of the sampler device that precedes the plant is accommodated to the variation of ; therefore, its discrete-time representation is that of a LTV system.
Remark 3.1. The above described operation may yield major changes in the dynamic response of the system. As an example, consider that the plant in Figure 2 is Let the repetitive controller be constructed assuming and a reference/disturbance nominal period s, this yielding a nominal sampling period of s. Let also , , and be selected as indicated in Section 4.2. Then, the sufficient stability conditions established in Proposition 2.1 for the nominal sampling period are fulfilled. Figure 3 depicts the maximum modulus eigenvalue of the closed-loop system matrix as a function of in a neighborhood of . This information allows stating that instability is ensured for ms.
The above discussion indicates that, although the proposed technique allows adapting the system to the specific signal frequency to be tracked/rejected without changing the digital controller, the sampling rate change may affect closed-loop stability. In what follows, an LMI-based approach to study closed-loop stability under varying sampling rate condition is developed.
Regarding the time-varying nature of the sampling period of the plant, the stability analysis is carried out in the state-space formalism. Let represent the continuous-time plant state-space representation, that is, Assume that is sampled at , with and , the sampling periods being . Let also and , where Therefore, the plant evolution at the sampling instants is given by the discrete-time LTV system: In case that remains constant, i.e., , for all , system (3.4) corresponds to a discrete-time LTI system with -transform transfer function . In an aperiodic sampling framework and vary with , and the -transform representation is no longer valid.
Similarly, under varying sampling period, the closed-loop system depicted in Figure 2 can be described in state space by a quadruple , which may be constructed combining and the state-space representations of , , and . Thus, and are constant matrices, while and depend continuously on . A methodology for studying the closed-loop system under time-varying sampling conditions will be introduced in the rest of this section.
Let the sampling period, , take values in a known compact subset .
Proposition 3.2. The uniform exponential stability of the zero state of implies the uniform BIBO stability of the system .
Proof. According to , the result follows if , , and are uniformly bounded matrices, for all , and this is indeed true: depends continuously on , which belongs to a compact set , while are constant matrices.
The preceding result allows reducing the stability analysis of to that of the zero state of .
Let us define
Proposition 3.3 (see ). If there exists a matrix such that where , then the zero state of is uniformly exponentially stable.
At this point, it is immediate to realize that relation (3.6) in Proposition 3.3 yields an infinite set of LMIs. The gridding approach introduced in [17, 18] allows a simplified stability analysis that may be performed in two stages, if necessary.
In a first stage, advantage is taken from the fact that , corresponding to evaluated in , is stable by construction, because recall that , and are designed with such purpose.
Proposition 3.4. Assume that the stability conditions of Proposition 2.1 are satisfied for a nominal sampling period . Then, (1)the zero state of the LTI system is uniformly exponentially stable,(2)the LMI problem with , is feasible,(3)let be a solution of the LMI problem (3.7) for a fixed . Then, there exists an open neighborhood of , say , such that is BIBO stable in .
Proof. It follows from the stability hypothesis that all the eigenvalues of have modulus less than 1, which yields immediately item 1 . Item 2 stems from the fact that the sufficient condition for uniform exponential stability established in Proposition 3.3 is also necessary for a discrete-time LTI system. Finally, item 3 follows immediately from Propositions 3.2 and 3.3 once the continuity of the matrix elements of with respect to is taken into account.
Assume that we are interested in analyzing the stability of for all sampling periods . Let then be a feasible solution of the LMI problem (3.7) for a fixed ; its existence is guaranteed by Proposition 3.4. Let also , with , be a sufficiently fine grid of . If , for all , stability for all may be probably inferred.
Otherwise, in case that there exists at least a single such that , the gridding procedure proposed in  may be carried out as follows. Let be a sorted set of candidate sampling periods suitably distributed in . Then, one may solve the following finite set of LMIs: for a fixed . In case that the problem is feasible and a solution, , is encountered, the negative semidefinite character of is to be checked for intermediate values of in each open subinterval . If this fails to be accomplished, (3.8) has to be solved again for a finer grid of .
Remark 3.5. The application of the procedures to the plant introduced in Remark 3.1 yields the following results. The evolution of the maximum modulus eigenvalue of being the solution of (3.7) for , is depicted in Figure 4. Hence, BIBO stability is ensured in ms. Alternatively, the approach that stems from (3.8) guarantees stability in ms; this interval is computed by successive extensions of an initial containing until (3.8), with , becomes definitely unfeasible.
Remark 3.6. Notice that(i)comparatively, in the first approach a shorter stability radius may be reasonably expected;(ii)the main drawback of both analysis paths is associated with the fact that stability is rigorously guaranteed just in open neighborhoods of each sampling period where either or are satisfied;(iii)the applicability of the above developed method to large scale or very high order systems is subject to the limitations of the numerical algorithms used to solve the LMI problems (3.7) or (3.8);(iv)a faulty estimation of and/or an error in the implementation of may yield important performance degradation . However, closed-loop BIBO stability is not threatened unless the plant is actually sampled with values lying outside the region where stability is guaranteed.
4. Experimental Setup and Results
4.1. Plant Description
Systems with rotary elements are usually affected by periodic disturbances due to the movement of these parts (e.g., electrical machines, CD players…). This kind of system is supposed to be moving, in some cases, at a fixed angular speed. Under these working conditions any friction, unbalance, or asymmetry appearing on the system generates a periodic disturbance that affects its dynamical behavior. Reference  contains a description of a mechatronic plant designed to reproduce this working conditions. This device is composed of a bar holding a permanent magnet in each end, each magnet being magnetically oriented in the opposite way, and attached to a DC motor and two fixed permanent magnets (see a sketch in Figure 5). The rotation of the DC motor causes a pulsating load torque () that depends on the mechanical angle of the motor axis. When the motor axis angular speed is constant (), the pulsating torque is a periodic signal with a fundamental period directly related to the axis speed: , with expressed in . The control goal for this plant is maintaining the motor axis angular speed constant at a desired value. However, this type of disturbances may not be completely rejected by a regular controller, so a repetitive controller is designed in next subsection.
Figure 6 also shows the open-loop time response of the plant containing the fixed magnets. It is important to note that the speed describes an almost periodic signal. This type of disturbances may not be completely rejected by a regular controller, so a repetitive controller is designed in next subsection.
Finally, it is worth mentioning that in this experimental setup the disturbance could be written as a nonlinear part of the system, and therefore nonlinear control techniques should be used instead of repetitive control. An interesting application of this procedure to speed control of mechatronic systems is reported in . However, for the present case this would require the use of complex models and, consequently, of complex nonlinear controllers which in general might not provide the same performance. Hence, we proceed assuming unmodeled or uncertainty components as disturbances, which is quite common in control theory and, at the same time, allows the use of repetitive control in this problem.
4.2. Control Design
A first-order model for the plant is experimentally derived from its open-loop time response without the fixed magnets, namely, The controller is constructed from (4.1), for a speed of rev/s and obtaining 25 samples per period, that is, . These conditions imply a nominal sampling period of ms. Under these assumptions the nominal discrete time plant is According to Remark 2.2, the following design issues have been taken into account: (i) provides a very robust inner loop, (ii)the first order null-phase FIR filter provides good performance in this case,(iii)the fact that is minimum phase allows , with .
These settings yield the control law with , where is the system output (speed) and is the reference.
4.3. Stability Analysis
Although the controller is designed to regulate the speed at 8 rev/s, in practice it will be necessary to move from this design point. Let us assume that we are interested in varying the speed reference in the interval rev/s; this entails a sampling period variation in the interval s.
The stability analysis that stems from Proposition 3.4 includes the solution of the LMI (3.7), which is known to be feasible, and the checking of the negative definite character of . Figure 7 shows the evolution of the maximum modulus eigenvalue of when solving for and also for 50000 uniformly distributed values of . Therefore, it can be presumed that the closed-loop system may operate in a speed range of rev/s with dynamically preserved stability. This speed interval is obviously very narrow and operation conditions are limited to a sampling period interval such that . However, it is important to recall that this test comes not from a necessary condition but from a sufficient condition, so moving out of this interval does not necessarily imply instability (although this step can be avoidable in most cases, its results are very useful to tune the LMI solver regarding next step).
In order to guarantee a broader stability interval, the second method described in Section 3 may be applied. Therefore, 50 uniformly distributed points are selected in s. These points are used to construct the set of LMIs (3.8), and a feasible solution with is obtained. Figure 7 depicts the maximum modulus eigenvalue of , designating with a star the 50 points leading to the LMI formulation. The maximum modulus eigenvalue of corresponding to a finer grid consisting of 50000 uniformly distributed point are also drawn in Figure 7. These points are used to check the sign of in the intervals between the points defining the LMI set. It can be seen that for every point in this finer grid of the interval ; hence, stability is dynamically preserved therein. This method extends the previously obtained stability interval rev/s, thus providing less conservative results. Further extensions of the new interval could also be feasible.
4.4. Implementation Issues
The first issue to deal with in the implementation of the control action is the measurement of the uncertain or time-varying period of the signal to be tracked/rejected. In the general case an adaptive scheme similar to those in [13, 15, 16] is used for this task; the complete controller architecture is depicted in Figure 8. The frequency values, , are obtained by means of a frequency observer using information from different sources, namely, reference profile, output signal, and control action. Then, the sampling rate is calculated as .
Nevertheless, as it has been previously stated, in mechanical turning systems the frequency to be tracked/rejected is directly related to the turning speed. Hence, for the analyzed mechatronic plant the disturbance frequency is straightforwardly computed from the turning speed reference. This allows us to study the effect of the adaptive approach on the closed-loop system decoupled from the frequency observer dynamics.
For the general case, the study of the effect of the frequency observer dynamics on the global system stability is out of the scope of this paper. However, it is also worth recalling that any fault in the estimation and/or implementation of causes a difference between the experimental ratio and the implemented value for which, in turn, may yield a performance degradation. As an example, Figure 9 shows the first harmonic gain factor evolution of the internal model against a relative deviation of the real sampling period with respect to the nominal ms, namely, where is defined in (2.5) and Notice that even small deviations of entail important gain reductions. However, as stated in Remark 3.6(iv), closed-loop stability is not threatened unless .
Finally, and also in accordance with Remark 3.6(iv), the frequency observer output has to be saturated so as to guarantee that .
The controllers have been implemented using RTLinuxFree-3.1, a real-time operating system (RTOS). This RTOS allows programing using standard POSIX system calls  with the user interface of a GNU/Linux Operating System (kernel 2.6.9). The controllers are encoded in C language, while periodic task handling has been implemented using RTLinuxFree capabilities.
The actuator is composed of a Johnson Electric HC615L (20 W, 5 A) DC motor, a full-bridge inverter (LMD18200, 3 A full-bridge) and pulse width modulation (PWM) hardware, included in a NuDQ PCI8133 board. The speed sensor is composed of a 500 pulse-per-revolution optical encoder Hewlett-Packard H9730A, a quadrature decoder and a counter, all of them are included in a NuDQ PCI8133 board. The interface between the controller real-time task and the NuDQ PCI8133 board has been implemented using low level commands.
4.5. Experimental Results
The reference profile used in the experimentation is the following. During the time interval s, the reference is maintained constant at the nominal value of rev/s. At s a ramp reference change, from rev/s to rev/s, is introduced in the system; then, the speed is kept constant for 10 s and finally at s the speed is gradually augmented at a constant acceleration until it reaches the value rev/s at s.
The experiment has been carried out for three different settings of the sampling time. In Figure 10(a) the sampling period is kept constant at the nominal value ms, for all ; it is important to realize that, in comparison to the uncompensated speed profile depicted in Figure 6, now disturbances are almost rejected in , that is, when both reference and sampling are at the nominal values. However, for s disturbances cannot be properly compensated, and performance is strongly degraded. Figure 10(b) depicts the output behavior when the sampling period is varied adaptively assuming that the estimation of uses a 2nd order frequency observer; the profile of the actually used is in Figure 11. Notice that the performance gets worse in the regions where there is estimation error. Finally, Figure 10(c) portrays the response under adaptive variation of from an exact estimation of . According to Sections 4.3 and 4.4, stability is preserved in the three situations because always belongs to . Performance in the steady-state depends on the accuracy of the estimation of , while performance during transients, although not guaranteed by repetitive control theory, follows here the same pattern as that observed for the steady state; that is, it improves with better estimations of .
This paper analyzed the BIBO stability of a closed-loop system containing a digital repetitive controller working under time-varying sampling period. The analysis was carried out using an LMI gridding approach that allows a stability assessment in a known, bounded interval where the reference/disturbance period is assumed to vary. The theoretically predicted results have been experimentally validated through a mechatronic plant in which rejection of periodic disturbances with time-varying period is successfully achieved. A detailed description of fundamental issues related to the implementation procedure has been provided.
Current research is being devoted to develop redesign methodologies to address the situations in which the sampling period falls out of the stability interval .
This work is partially supported by the spanish Ministerio de Ciencia e Innovación under Project DPI2010-15110.
- Y. Yamamoto, “Learning control and related problems in infinite-dimensional systems,” in Proceedings of the European Control Conference, pp. 191–222, June 1993.
- G. Hillerström and K. Walgama, “Repetitive control theory and applications—a survey,” in Proceedings of the 13th IFAC World Congress, vol. D, pp. 1–6, San Francisco, Calif, USA, July 1996.
- B. A. Francis and W. M. Wonham, “Internal model principle of control theory,” Automatica, vol. 12, no. 5, pp. 457–465, 1976.
- K. Chew and M. Tomizuka, “Digital control of repetitive errors in disk drive systems,” IEEE Control Systems Magazine, vol. 10, no. 1, pp. 16–20, 1990.
- M. Yamada, Z. Riadh, and Y. Funahashi, “Design of robust repetitive control system for multiple periods,” in Proceedings of the 39th IEEE Conference on Decision and Control, vol. 4, pp. 3739–3744, December 2000.
- D. H. Kim and T. C. Tsao, “Robust performance control of electrohydraulic actuators for electronic cam motion generation,” IEEE Transactions on Control Systems Technology, vol. 8, no. 2, pp. 220–227, 2000.
- K. Zhou, D. Wang, and G. Xu, “Repetitive controlled three-phase reversible PWM rectifier,” in Proceedings of the American Control Conference, pp. 125–129, June 2000.
- K. Zhou, D. Wang, and K. Low, “Periodic errors elimination in CVCF PWM DC/AC converter systems: repetitive control approach,” IEEE Proceedings, vol. 147, no. 6, pp. 694–700, 2000.
- K. Zhou and D. Wang, “Digital repetitive learning controller for three-phase CVCF PWM inverter,” IEEE Transactions on Industrial Electronics, vol. 48, no. 4, pp. 820–830, 2001.
- R. Costa-Castelló, R. Griñó, and E. Fossas, “odd-harmonic digital repetitive control of a single-phase current active filter,” IEEE Transactions on Power Electronics, vol. 19, no. 4, pp. 1060–1068, 2004.
- M. Steinbuch, “Repetitive control for systems with uncertain period-time,” Automatica, vol. 38, no. 12, pp. 2103–2109, 2002.
- M. Steinbuch, S. Weiland, and T. Singh, “Design of noise and period-time robust high-order repetitive control, with application to optical storage,” Automatica, vol. 43, no. 12, pp. 2086–2095, 2007.
- Z. Cao and G. F. Ledwich, “Adaptive repetitive control to track variable periodic signals with fixed sampling rate,” IEEE/ASME Transactions on Mechatronics, vol. 7, no. 3, pp. 378–384, 2002.
- G. Hillerström, On Repetitive Control, Luleå University of Technology, 1994.
- T. C. Tsao, Y. X. Qian, and M. Nemani, “Repetitive control for asymptotic tracking of periodic signals with an unknown period,” Journal of Dynamic Systems, Measurement and Control, vol. 122, no. 2, pp. 364–369, 2000.
- R. D. Hanson and T. C. Tsao, “Periodic sampling interval repetitive control and its application to variable spindle speed noncircular turning process,” Journal of Dynamic Systems, Measurement and Control, vol. 122, no. 3, pp. 560–566, 2000.
- P. Apkarian and R. J. Adams, “Advanced gain-scheduling techniques for uncertain systems,” IEEE Transactions on Control Systems Technology, vol. 6, no. 1, pp. 21–32, 1998.
- A. Sala, “Computer control under time-varying sampling period: an LMI gridding approach,” Automatica, vol. 41, no. 12, pp. 2077–2082, 2005.
- R. Costa-Castelló, J. Nebot, and R. Griñó, “Demonstration of the internal model principle by digital repetitive control of an educational laboratory plant,” IEEE Transactions on Education, vol. 48, no. 1, pp. 73–80, 2005.
- T. Inoue, M. Nakano, T. Kubo, S. Matsumoto, and H. Baba, “High accuracy control of a proton synchroton magnet power supply,” in Proceedings of the 8th IFAC World Congress, pp. 216–220, 1981.
- R. Griñó and R. Costa-Castelló, “Digital repetitive plug-in controller for odd-harmonic periodic references and disturbances,” Automatica, vol. 41, no. 1, pp. 153–157, 2005.
- G. Escobar, P. G. Hernandez-Briones, R. E. Torres-Olguin, and A. A. Valdez, “A repetitive-based controller for the compensation of harmonic components,” in Proceedings of the IEEE International Symposium on Industrial Electronics, pp. 3397–3402, June 2007.
- K. K. Chew and M. Tomizuka, “Steady-state and stochastic performance of a modified discrete-time prototype repetitive controller,” Journal of Dynamic Systems, Measurement and Control, vol. 112, no. 1, pp. 35–41, 1990.
- M. Tomizuka, T. C. Tsao, and K. K. Chew, “Analysis and synthesis of discrete-time repetitive controllers,” Journal of Dynamic Systems, Measurement and Control, vol. 111, no. 3, pp. 353–358, 1989.
- G. Hillerström and R. C. H. Lee, “Trade-offs in repetitive control,” Tech. Rep. CUED/F-INFENG/TR 294, University of Cambridge, 1997.
- W. J. Rugh, Linear System Theory, Prentice-Hall, Upper Saddle River, NJ, USA, 2nd edition, 1996.
- B. Mahawan and Z. H. Luo, “Repetitive control of tracking systems with time-varying periodic references,” International Journal of Control, vol. 73, no. 3, pp. 1–10, 2000.
- B. Gallmeister, POSIX.4 Programming for the Real World, O'Reilly & Associates, 1st edition, 1995.
Copyright © 2011 Germán A. Ramos 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.