Journal of Control Science and Engineering
Volume 2016 (2016), Article ID 7256039, 11 pages
http://dx.doi.org/10.1155/2016/7256039

Research Article

Optimal Inversion Feedforward and Robust Feedback Based 2DOF Control Approach for High Speed-Precision Positioning Systems

Chao Peng,1,2 Chongwei Han,² Jianxiao Zou,¹ and Guanghui Zhang²

¹School of Mechanical Engineering, University of Electronic Science and Technology of China, Chengdu, Sichuan 611731, China
²Northwest Institute of Mechanical & Electrical Engineering, Xianyang, Shanxi 712099, China

Received 24 July 2015; Revised 30 October 2015; Accepted 10 November 2015

Academic Editor: Petko Petkov

Copyright © 2016 Chao Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper proposed a novel optimal inversion feedforward and robust feedback based two-freedom-of-freedom (2DOF) control approach to address the positioning error caused by system uncertainties in high speed-precision positioning system. To minimize the norm of the positioning error in the presence of model uncertainty, a linear matrix inequality (LMI) synthesis approach for optimal inversion feedforward controller design is presented. The specification of position resolution, control width, robustness, and output signal magnitude imposed on the entire 2DOF control system are taken as optimization objectives of feedback controller design. The robust feedback controller design approach integrates with feedforward controller systematically and is obtained via LMI optimization. The proposed approach was illustrated through a simulation example of nanopositioning control in atomic force microscope (AFM); the experiment results demonstrated that the proposed 2DOF control approach not only achieves the performance specification but also could improve the positioning control performance compared with mixed sensitivity feedback control and inversion-based 2DOF control.

1. Introduction

The performance of high speed-precision positioning system such as piezoelectric and precision motor actuated system is required to meet the specification in scanning probe microscopy (SPM) [1, 2], microelectromechanical system [3, 4], optical-hard disk drive [5, 6], and so on. The positioning control performance could be characterized by positioning resolution, control bandwidth, and system robustness [1–6]. Various feedback control approaches have been studied and demonstrated that it could render requirements on positioning resolution and robustness satisfied at low frequencies [7–9]. However, the positioning resolution and control bandwidth of feedback control are limited for the Bode Sensitivity Integral theorem [10, 11]. It has been demonstrated that the feedforward control could increase the feedback control bandwidth [12, 13]. To improve the system positioning control performance, considerable 2DOF control approaches which combine the feedforward controller and feedback controller have been studied [14–18].

The inversion control technology could achieve exact trajectory tracking and precision positioning for linear and nonlinear system without system uncertainty [19, 20]. The stable-inversion [21], preview-based optimal inversion [22], and inversion-based iterative control [23] approach have improved the practicality of inversion technology. However, the positioning control performance of inversion feedforward control is limited by modeling uncertainties and disturbances [24]. Thus, recently, the inversion feedforward based 2DOF control approaches which utilize the robust feedback control to compensated the limitation of inversion feedforward controller have been presented and demonstrated that they could achieve good positioning performance [25–27]. However, the feedforward controller and feedback controller in the above 2DOF control approaches are designed separately; the feedback controller is designed by conventional mixed sensitivity feedback control approach without taking the effect of feedforward controller into account. There exist challenges in inversion feedforward based 2DOF control: (1) inversion feedforward controller design optimization with considering the system uncertainty and (2) how to integrate the feedforward controller and feedback controller design systematically to obtain the desired requirement positioning resolution, control bandwidth, and robustness of the entire 2DOF control system.

To minimize the adverse effect caused by system uncertainties and make the requirements of performance of the entire 2DOF control system satisfied, a novel optimal inversion feedforward and robust feedback based 2DOF control approach is proposed in this paper. The main contribution of this paper is as follows. (1) A new LMI representation for optimal inversion feedforward controller design which takes minimizing the norm of the positioning error caused by feedforward control in presence of the system uncertainty as optimization objective. (2) A 2DOF control positioning performance optimization problem based on mixed sensitivity is formulated and a LMI optimization based robust feedback controller design approach integrating with feedforward controller systematically. The proposed 2DOF control design approach is evaluated through the experiment.

The rest of this paper is organized as follows. The design objective of optimal inversion feedforward and robust feedback based 2DOF control is formulated in Section 2. In Section 3, the LMI based optimal inversion feedforward controller design approach is presented. The design of robust feedback controller is given in Section 4. The proposed 2DOF control approach is illustrated through AFM simulation experiment in Section 5. Finally, the conclusion is discussed in Section 6.

2. Problem Formulation

2.1. System Description

Consider the SISO LTI system with parameter uncertainty:

where and are polynomials of Laplace variable s and is the uncertain parameter vector that parameterizes the transfer function .

Defining the nominal parameter vector of as , the nominal model of the system could be denoted as .

We assume that is known to lie in a box domain which is defined as

where are the variation of .

Assumption 1. System is invertible, and the relative degree of is at least two (i.e., has at least two more poles than zeros). System and its inverse are hyperbolic; that is, does not have poles or zeros on the imaginary axis of the complex plane.

Remark 2. The requirements that the system is invertible and hyperbolic and that the relative degree of is more than one are needed for computation [21] and robustness [22] of the exact inverse.

2.2. Optimal Inversion Feedforward and Robust Feedback Based 2DOF Control

Consider the 2DOF control system shown in Figure 1. In this figure, is the transfer function of a LTI plant as described in (1), and are feedforward controller and feedback controller, respectively. The signal represents the desired outputting and represents the input to the plant . The signal represents the actual output of the entire 2DOF control system.

Figure 1: Block diagram of 2DOF control system.

The transfer function from the desired trajectory to actual outputs and is given by

where is the feedback sensitivity function; that is,

The outputting error of entire 2DOF control system can be decoupled as multiplication of the feedforward path outputting error and the feedback sensitivity function :

Remark 3. By Bode Sensitivity Integral theorem, the small feedforward path output error could improve the bandwidth of feedback control and output error. If , the output error could be zero; however, it is impossible to find the exact inverse of the system for the modeling uncertainty.
The control block and output error of the 2DOF control system are now presented. The design goal of optimal feedforward and robust feedback based 2DOF control design is to make the following two optimization objects satisfied.
(1) Find an optimal robust inversion-based feedforward controller to minimize the norm of feedforward path output error in the presence of system uncertainties; that is,

(2) After determining the feedforward controller , find a robust feedback controller to minimize the norm of transfer function ; that is,

where

denotes the weighting function for output bandwidth and output error limitation, denotes the weighting function for robustness performance of system, denotes the weighting function for magnitude of output signal of feedback controller, and and are the transfer function from desired trajectory to output error and and to actual output .

3. Optimal Inversion Feedforward Controller Design

This section is devoted to optimal inversion feedforward controller design. At first feedforward controller based on optimal inverse of is designed by LMI optimization, and then the feedforward controller is implemented by preview-based inversion [21].

3.1. Optimal Inversion of System

A modulation function is introduced to inversion feedforward controller design as follows:

Consider Assumption 1; is improper; that is, the poles number is less than zeros number. Assume is proper; thus, inversion feedforward controller also could be decoupled as proper element :

where and are the denominator and numerator of , respectively, The order of is equivalent to the order of . Thus, the relative degree of is zero and the order of is equal to the relative degree of nominal model .

The transfer function of feedforward path output error could be rewritten as

Remark 4. Because the order of is equal to the relative degree of nominal model , the relative degree of is equal to zero. The modulation function could be obtained as .
Now, a proper transfer function whose relative degree is zero can be obtained by multiplying improper element with . Consider the parameter uncertainty domain of , the parameters of can also be known to lie in a box uncertainty domain, and the state space matrices of which depend on [28, 29]. The sate space representation of can be given as

The state space representation of is given as

Now, the state space representation of the feedforward path output error can be written as

where the state space matrices are given by

The objective of the inversion feedforward controller design in (6) could be transformed to find an optimal function; that is,

The bounded real lemma (BRL) [30] will be utilized for solving the optimal . The LMI representation lemma is given as follows.

Lemma 5. A proper element in inversion feedforward controller as in (10) is the solution to the minimization problem in (16) whose state space matrices , and and a symmetric matrix could minimize the scale in presence of uncertain parameters, subjected to following parameterized LMI, that is,

subject to

Proof. Consider the uncertainty domain and bounded real lemma; if there exist a positive-definite symmetric matrix and matrices , and which could render the parameterized LMIs in (18) satisfied, norm of feedforward path positioning error is less than in the presence of uncertain parameters ; that is, .
Thus, if matrices , , , and and a given and could satisfy the optimization objective in (17), a solution for minimizing problem in (16) exists. This completes the proof.

The solution to the optimization problem in (6) could be given by the following theorem.

Theorem 6. Assume the optimal inversion feedforward controller with form as (10) is .
(a) The controller that could solve the minimization problem in (6) exits, if there exists a solution where , for for the following optimization:

subject to

(b) The proper element in , , is given as follows:

Considering the partition of and , we introduce a partition of and its inverse :

Since ,

Define a full ranked and square matrix, which will be used in following matrix inequalities transformation:

To obtain a feasible solution which could render the LMIs in (18) satisfied and consider the partition of in (22), we introduce full ranked matrices , and and assume , and have following form:

Multiplying the inequalities in (18) by and on the left and on the right, respectively, we obtain the following LMI:

Now, the LMI item in (26) could be written as

Replacing the LMI items in (26) with (27), we obtain the LMI as in (20). This completes the proof.
(b) It is obvious that as (25) is the solution to (16) if the matrix variables , and are the solution for (19); and , could be deduced from (23). This completes the proof.

Now, the optimal inversion feedforward controller which satisfies the optimization problem in (6) can be obtained by Theorem 6.

3.2. Time Domain Implementation of Feedforward Controller

Note that the obtained optimal inversion feedforward controller as is improper; it will be realized by preview-based inversion [22]. The design steps are given as follows.

(1) The improper inversion feedforward controller can be transformed to proper form as by a transformed desired outputting which carries the preview knowledge of . The transformed desired outputting is obtained by multiplying the improper element and desired outputting . The output of feedforward controller is given as follows:

(2) By the proper and , the minimal state space realization of the optimal inversion feedforward controller can be given as

(3) Finally, the time domain representation of output of the optimal feedforward controller can be obtained as

By the above step, the optimal inversion is implemented in time domain.

4. Robust Feedback Controller Design

This section will discuss the solution to the robust feedback controller design problem in (7). The mixed sensitivity synthesis scheme of 2DOF control system is shown in Figure 2.

Figure 2: mixed sensitivity synthesis scheme.

By Figure 2, the sensitivity transfer functions of the entire 2DOF control system can be represented as follows:

Integrating with the feedforward which has been determined by design approach proposed in Section 3 and the weighting functions , and , which impose the requirements for positioning resolution performance, robustness of the entire 2DOF control system, and output signal magnitude of feedback controller, find a feedback controller which could make the optimization problem in (7) satisfied. The design problem of feedback controller could be transformed to seek a controller which could minimize the norm of transfer function from desired trajectory to weighting output as in Figure 2; that is,

The transfer function matrix of in Figure 2 is given as follows:

The state space representation of could be given as

The state space representation of feedback controller is given as follows:

Theorem 7. Consider the system as (34), if there exist symmetric matrices and and matrices with appropriate dimension , and that could satisfy the following optimization problem of positive scale , that is,

subject to

for ; the feedback controller which is the solution to the optimization problem

could be given as

where and are deduced from .

Proof. By BRL and the theorem in [30, 31], there exists a feedback controller such that for a certain parameter , if the LMIs in (37) hold, and the feedback controller is given in (39) for a system with certain parameter . It is obvious that if there exists a solution of the minimization problem for , the feedback controller which is given in (39) is the solution for the minimization problem in (38). This completes the proof.

Feedback controller which satisfies the optimization problem in (7) could be found by the following procedure.

Step 1. Discretize the uncertain parameters in box domain , set , , and initialize the parameter vector set ; ; , .

Step 2. By Theorem 7, obtain the feedback controller , which could satisfy the optimization problem and define .

Step 3. Define . By BRL, we can obtain . If , for , then and stop.

Step 4. Else , go to Step 2.

5. Implementation and Experiment

In this section, we will conduct a simulation example of nanopositioning in atomic force microscope (AFM) operation to illustrate the positioning control performance through the proposed 2DOF control design approach.

5.1. Experimental System Description

The AFM system utilizes piezo actuator to enable the - axis scanning of the AFM probe relative to the sample surface during AFM imaging as shown in Figure 3. The output error will cause AFM image distortion and sample or probe damaging during scanning. Thus, the positioning control performance of the piezo actuator is important in AFM which measures the surface properties.

Figure 3: AFM operation diagram and scanning path.

To demonstrate the proposed 2DOF control design approach, we will take the -axis scanning motion control in the simulation experiment in MATLAB Simulink. The model of piezo actuator always is identified by several frequency response measurements. According to frequency response data, an approximate transfer function is obtained, which has parameter uncertainties compared with actual system. Considering the uncertainties, a second-order model of piezo actuator is given as follows [32]:

The variation bound is given as follows:

where , ; , .

The frequency responses of piezo actuator are shown in Figure 4.

(a)

(b)

Figure 4: The frequency response of piezo actuator.

5.2. Design of the Feedforward Controller

The inversion feedforward controller is selected as follows:

The improper element could be written as

The transfer function in (11) and its state matrices in (12) could be given as follows:

By using Theorem 6 and LMI toolbox in MATLAB, the optimal inversion feedforward controller is obtained as follows which could render

The output of optimal inversion feedforward controller could be obtained by the obtained and preview inversion design steps as in (28)–(30).

5.3. Design of the Feedback Controller

At first, we specify the weighting function. It is noted that the tracking error could be achieved to around 30% by optimal inversion feedforward controller. We chose the weighting function for control bandwidth and positioning error limitation as which specifies the positioning error to be less than 2.5% below frequency at 125 Hz. The frequency response of is shown in Figure 5.

Figure 5: The frequency response of .

Considering system robustness requirement and the parameter uncertainty domain, the weighting function for system robustness performance is chosen as . The frequency response of is shown in Figure 6.

Figure 6: The frequency response of .

Considering the limitation of the controller output signal, the weighting function for magnitude of output signal of feedback controller is chosen to be a constant, .

Through design procedure of robust feedback controller and selecting , we obtained the feedback controller by using LMI toolbox in MATLAB. It is given as follows:

where

By mixed sensitivity design approach without considering the feedforward controller and the above-selected weighting function, we obtained the controller by using MATLAB mixed sensitivity synthesis command “hinf”. It is given as follows:

5.4. Experiment and Result Discussion

The AFM -axis positioning control simulation experiment which is tracking triangular trajectory at three different frequencies (5 Hz, 50 Hz, and 125 Hz) was conducted. For comparison, the positioning performances of three control approaches were tested by simulation experiment. The approaches are as follows: the proposed optimal inversion feedforward and robust feedback 2DOF control approach, that is, , exact inversion feedforward-robust feedback based 2DOF control approach, that is, , and    mixed sensitivity feedback control only, that is, . The positioning control performance was evaluated by positioning maximum error and mean error, which are described as follows:

Figure 7 shows the positioning control simulation results and positioning errors by using different control approaches. The maximum positioning errors and mean errors for triangular trajectory scanning at different frequencies are shown in Table 1. The range of triangular trajectory was 1 μm. As shown, at scan rate 5 Hz, these three control approaches could achieve precision positioning; maximum positioning errors are 0.0019 μm, 0.0027 μm, and 0.0138 μm, respectively, and mean positioning errors are 0.0002 μm, 0.0003 μm, and 0.0073 μm, respectively. At scan rate 50 Hz, the proposed 2DOF control approach and the exact inversion feedforward and robust feedback based 2DOF control approach could obtain a better positioning control performance than mixed sensitivity feedback control only. The maximum positioning errors obtained by approaches (1) and (2) are 0.0034 μm and 0.0061 μm which are as small as 13.2% and 23.6% of that obtained by mixed sensitivity feedback control. The positioning mean errors obtained by approaches (1) and (2) are 0.0003 μm and 0.0034 μm which are as small as 2% and 22% of that obtained by approach (3). At scan rate 125 Hz, only the proposed 2DOF control could obtain precision positioning. The maximum positioning error obtained by approach (1) is 0.0057 μm, which is as small as 55% of that obtain by approach (2) and is as small as 12.7% of that obtain by approach (3). The positioning mean error obtained by approach (1) is 0.0008 μm, which is as small as 14.8% of that obtained by approach (2) and is as small as 2.6% of that obtained by approach (3). Through the experiment result analysis, it is obvious that the proposed 2DOF control could achieve better positioning precision than the other two control approaches both in high speed scanning.

Table 1: Comparison of positioning control error.

Frequency
(um)
(um)

(1)
(2)
(3)
(1)
(2)
(3)

5 HZ
0.0019
0.0027
0.0138
0.0002
0.0003
0.0073

50 HZ
0.0034
0.0061
0.0258
0.0003
0.0034
0.0148

125 HZ
0.0057
0.0103
0.0447
0.0008
0.0054
0.0303

Figure 7: AFM -axis positioning control simulation experiment result.

6. Conclusion

In this paper, a novel optimal inversion feedforward and robust feedback based 2DOF control approach for high speed-precision positioning systems is proposed. In this approach, an optimal inversion feedforward controller is designed to minimize the norm of the tracking error in the presence of model uncertainty via linear matrix inequality (LMI) synthesis. Integrating with the feedforward controller systematically, a robust feedback controller is designed to render the requirements of positioning resolution, control width, robustness, and output signal magnitude imposed on the entire 2DOF control system satisfied via LMI optimization. The proposed control design approach is implemented in an AFM system -axis positioning simulation experiment, and the experiment results demonstrated the effectiveness of the proposed 2DOF control approach.

Conflict of Interests

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

Acknowledgment

This work is supported by China Postdoctoral Science Foundation under Grant 2014M560710.

References

J.-C. Shen, W.-Y. Jywe, H.-K. Chiang, and Y.-L. Shu, “Precision tracking control of a piezoelectric-actuated system,” Precision Engineering, vol. 32, no. 2, pp. 71–78, 2008.

G. Schitter, F. Allgöwer, and A. Stemmer, “A new control strategy for high-speed atomic force microscopy,” Nanotechnology, vol. 15, no. 1, pp. 108–114, 2004.

M.-Y. Chen and J.-S. Lu, “High-precision motion control for a linear permanent magnet iron core synchronous motor drive in position platform,” IEEE Transactions on Industrial Informatics, vol. 10, no. 1, pp. 99–108, 2014.

B. E. Helfrich, C. Lee, D. A. Bristow et al., “Combined -feedback control and iterative learning control design with application to nanopositioning systems,” IEEE Transactions on Control Systems Technology, vol. 18, no. 2, pp. 336–351, 2010.

K. Yang, Y. Choi, and W. K. Chung, “On the tracking performance improvement of optical disk drive servo systems using error-based disturbance observer,” IEEE Transactions on Industrial Electronics, vol. 52, no. 1, pp. 270–279, 2005.

C. Du, L. Xie, J. N. Teoh, and G. Guo, “An improved mixed H ₂/H ∞ control design for hard disk drives,” IEEE Transactions on Control Systems Technology, vol. 13, no. 5, pp. 832–839, 2005.

Z. Z. Liu, F. L. Luo, and M. A. Rahman, “Robust and precision motion control system of linear-motor direct drive for high-speed X-Y table positioning mechanism,” IEEE Transactions on Industrial Electronics, vol. 52, no. 5, pp. 1357–1363, 2005.

B. Yao and L. Xu, “Adaptive robust motion control of linear motors for precision manufacturing,” Mechatronics, vol. 12, no. 4, pp. 595–616, 2002.

A. C. Shegaonkar and S. M. Salapaka, “Making high resolution positioning independent of scan rates: a feedback approach,” Applied Physics Letters, vol. 91, no. 20, Article ID 203513, 2007.

J. S. Freudenberg, C. V. Hollot, R. H. Middleton, and V. Toochinda, “Fundamental design limitations of the general control configuration,” IEEE Transactions on Automatic Control, vol. 48, no. 8, pp. 1355–1370, 2003.

A. Karimi, M. Kunze, and R. Longchamp, “Robust controller design by linear programming with application to a double-axis positioning system,” Control Engineering Practice, vol. 15, no. 2, pp. 197–208, 2007.

T. Uchihashi, N. Kodera, H. Itoh, H. Yamashita, and T. Ando, “Feed-forward compensation for high-speed atomic force microscopy imaging of biomolecules,” Japanese Journal of Applied Physics, vol. 45, no. 3, pp. 1904–1908, 2006.

D. Croft and S. Devasia, “Vibration compensation for high speed scanning tunneling microscopy,” Review of Scientific Instruments, vol. 70, no. 12, pp. 4600–4605, 1999.

C. Lee, G. Mohan, and S. Salapaka, “2DOF control design for nanopositioning,” in Control Technologies for Emerging Micro and Nanoscale Systems, vol. 413 of Lecture Notes in Control and Information Sciences, pp. 67–82, Springer, Berlin, Germany, 2011.

G. Schitter, A. Stemmer, and F. Allgower, “Robust two-degree-of-freedom control of an atomic force microscope,” Asian Journal of Control, vol. 6, no. 2, pp. 156–163, 2004.

C. Lee and S. M. Salapaka, “Robust broadband nanopositioning: fundamental trade-offs, analysis, and design in a two-degree-of-freedom control framework,” Nanotechnology, vol. 20, no. 3, Article ID 035501, 2009.

K. Takanori, M. Yoshihiro, and I. Makoto, “LMI-based 2-degree-of-freedom controller design for robust vibration suppression positioning,” IEEJ Transactions on Industry Application, vol. 131, no. 1, pp. 93–101, 2011.

M. Araki and H. Taguchi, “Two-degree-of-freedom PID controllers,” International Journal of Control, Automation, and Systems, vol. 1, no. 4, pp. 401–411, 2003.

S. Devasia, D. Chen, and B. Paden, “Nonlinear inversion-based output tracking,” IEEE Transactions on Automatic Control, vol. 41, no. 7, pp. 930–942, 1996.

L. R. Hunt, G. Meyer, and R. Su, “Noncausal inverses for linear systems,” IEEE Transactions on Automatic Control, vol. 41, no. 4, pp. 608–611, 1996.

S. Devasia and B. Paden, “Stable inversion for nonlinear nonminimum phase time-varying system,” IEEE Transactions on Automatic Control, vol. 43, no. 2, pp. 283–288, 1998.

Q. Zou and S. Devasia, “Preview-based optimal inversion for output tracking: application to scanning tunneling microscopy,” IEEE Transactions on Control Systems Technology, vol. 12, no. 3, pp. 375–386, 2004.

Y. Yan, Y. Wu, Q. Zou, and C. Su, “An integrated approach to piezoactuator positioning in high-speed atomic force microscope imaging,” Review of Scientific Instruments, vol. 79, no. 7, Article ID 073704, 9 pages, 2008.

S. Devasia, “Should model-based inverse inputs be used as feedforward under plant uncertainty?” IEEE Transactions on Automatic Control, vol. 47, no. 11, pp. 1865–1871, 2002.

Q. Zou, K. K. Leang, and E. Sadoun, “Control issues in high-speed AFM for biological applications: collagen imaging example,” Asian Journal of Control, vol. 6, no. 2, pp. 164–178, 2004.

C. G. L. Bianco and A. Piazzi, “A servo control system design using dynamic inversion,” Control Engineering Practice, vol. 10, no. 8, pp. 847–855, 2002.

K. K. Leang and S. Devasia, “Feedback-linearized inverse feedforward for creep, hysteresis, and vibration compensation in AFM piezoactuators,” IEEE Transactions on Control Systems Technology, vol. 15, no. 5, pp. 927–935, 2007.

M. Tahk and J. L. Speyer, “Modeling of parameter variations and asymptotic LQG synthesis,” IEEE Transactions on Automatic Control, vol. 32, no. 9, pp. 793–801, 1987.

G. Chesi, “Estimating the domain of attraction for uncertain polynomial systems,” Automatica, vol. 40, no. 11, pp. 1981–1986, 2004.

Y. He, M. Wu, and J.-H. She, “Improved bounded-real-lemma representation and control of systems with polytopic uncertainties,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 52, no. 7, pp. 380–383, 2005.

M.-N. Lee, J.-H. Moon, K. B. Jin, and M. J. Chung, “Robust H ∞ control with multiple constraints for the track-following system of an optical disk drive,” IEEE Transactions on Industrial Electronics, vol. 45, no. 4, pp. 638–645, 1998.

Q. Xu and Y. Li, “Model predictive discrete-time sliding mode control of a nanopositioning piezostage without modeling hysteresis,” IEEE Transactions on Control Systems Technology, vol. 20, no. 4, pp. 983–994, 2012.