Journal of Applied Mathematics

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

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

Research Article

Generalized Kalman-Yakubovich-Popov Lemma Based I-PD Controller Design for Ball and Plate System

Shuichi Mochizuki¹ and Hiroyuki Ichihara²

¹Graduate School of Science and Technology, Meiji University, 1-1-1 Higashimita, Tama-ku, Kawasaki-shi, Kanagawa 214-8571, Japan

²Department of Science and Engineering, Meiji University, 1-1-1 Higashimita, Tama-ku, Kawasaki-shi, Kanagawa 214-8571, Japan

Received 23 June 2013; Accepted 9 September 2013

Academic Editor: Baocang Ding

Copyright © 2013 Shuichi Mochizuki and Hiroyuki Ichihara. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The ball-on-plate balancing system has a camera that captures the ball position and a plate whose inclination angles are limited. This paper proposes a PID controller design method for the ball and plate system based on the generalized Kalman-Yakubovich-Popov lemma. The design method has two features: first, the structure of the controller called I-PD prevents large input signals against major changes in the reference signal; second, a low-pass filter is introduced into the feedback loop to reduce the influence of the measurement noise produced by the camera. Both simulations and experiments are used to evaluate the effectiveness of the design method.

1. Introduction

The ball and plate [1] is an unstable underactuated nonlinear system that has double integrators at the origin and that has two control inputs against four degrees of freedom (DOF). A camera located above the plate captures the position of the ball, and two motors manipulate the inclination angles of the plate to keep the ball on the plate. The ball and plate system is an extension of the ball and beam system [2] from one to two dimensions. The system has demonstrated various controller design methods for positioning and trajectory tracking of the ball: proportional integral derivative (PID) control [3], fuzzy control [4], neural network control [3, 5], and model predictive control [6]. In particular, PID control has the benefits of simple implementation and fewer hardware requirements, and it has been applied in many successful designs [7]. Because PID control enables us a limited performance, optimizing the parameters in a PID controller satisfying design specifications is an important subject for study. In the controller design of the ball and plate, it requires to consider limitation of the inclination angles with good transient and steady-state responses. Although proportional and derivative controllers are required to improve transient responses, a jump in the reference signal generates a large input signal that reaches the limitation angle that degrades the transient responses [8]. In addition, ball position data from the camera include measurement noise that also degrades the steady-state responses.

To overcome the above issues, this paper proposes a PID controller design method for the ball and plate system by open-loop shaping based on the generalized Kalman-Yakubovich-Popov (GKYP) lemma [9]. The GKYP lemma is a generalization of the standard KYP lemma [10], which establishes the equivalence between a frequency domain inequality (FDI) for a transfer function and a linear matrix inequality (LMI) associated with its state-space realization. The standard KYP lemma is available for the infinite frequency range while the generalized one can limit the frequency range to be (semi) finite. By introducing the GKYP lemma to PID controller design, design specifications by FDIs in the finite frequency ranges for the open-loop transfer function result in LMIs [9]. The GKYP lemma gives a systematic open-loop shaping design method through optimization to realize desirable transient and steady-state responses. In this visual feedback system, we introduce a low-pass filter. Since the filter gives freedom in optimization, it allows better steady-state responses and reduces the influence of the measurement noise. To prevent large input signals from P and D controllers, we adopt the I-PD (integral-proportional derivative) structure, whose design is still in the framework of the GKYP lemma, because the open-loop transfer functions of PID and I-PD structures are fundamentally the same.

The paper is organized as follows. The description of the ball and plate system, including its modeling and the measurement noise, is presented in Section 2. The GKYP lemma based I-PD controller design method with a low-pass filter is provided in Section 3. The design of the I-PD controller is described in Section 4. Simulation and experimental results are presented in Sections 5 and 6, respectively. Finally, in Section 7, we present our conclusions.

The notation used is standard. For a matrix , the transpose and complex conjugate transpose are denoted by and , respectively. For a Hermitian matrix , and denote positive (semi) definiteness and negative (semi) definiteness, respectively. The symbol stands for the set of Hermitian matrices. The subscript is omitted if . The real and imaginary parts of are denoted by and . For matrices and , denotes the Kronecker product. represents the Laplace transform of a signal .

2. Ball and Plate System

The ball and plate, a QUANSER 2D Ball Balancer, is shown in Figure 1. The system consists of a plate, a ball, an overhead camera, and two servo units. The plate is allowed to swivel in both the - and -directions. The overhead CMOS digital camera, a Point Grey Research Inc. FFMV-03M2C-CS, measures the position of the ball. The two servo units located under the plate are QUANSER SRV02 devices, each of which has a peak time of approximately 200 ([ms]) and an overshoot of approximately 5%. Each of the devices is connected to a side of the plate through a two DOF gimbal. The sampling time of the control system and the frame rate provided by the camera are 1 ([ms]) and 60 ([fps]), respectively. Thus the image information is renewed approximately every 17 ([ms]). There is a constant time delay of less than 60 ([ms]) between the measurement of the ball position and the manipulation of the servo units in the visual feedback system.

Figure 1: 2D Ball Balancer.

2.1. Modeling

The -direction of the ball and plate system is illustrated in Figure 2. We assume that the ball is completely symmetric and homogeneous and does not slip on the plate and that all frictions are neglected. The plate rotates in the -Cartesian coordinates with the origin at the center of the plate. The equations of motion are

where is the position of the ball on the plate, and are the inclination angles of the plate to the - and -axis, respectively, is the mass of the ball, is the radius of the ball, is the gravitational acceleration, and is the inertia of the ball. In Figure 2, represents the angle of the load gear. The relationship between and is as follows:

where is the length of the side of the plate and is the length between the joint and the center of the load gear. The relationship of and is the same as (2), since both gear systems have the same hardware and the plate is symmetrical. The numerical values of the constant parameters in the equations of motion and (2) are shown in Table 1. Since and are limited as

from (2), the working ranges of and are

Table 1: Parameters of the ball and plate system.

Parameters
Numerical values

0.0252 [kg]

0.0170 [m]

9.81 [m/s²]

[kg·m]

0.275 [m]

0.0254 [m]

Figure 2: The ball and plate system.

If the angular velocities and are relatively low, the approximations

are often used. Linearizing the equations of motion at and , we have

Since the axes are independent of each other, we can focus on one axis, for example, the -axis. For the input and the output , the transfer function is given by

where , , and

2.2. Measurement Noise

In this visual feedback system, there is inevitable noise from the camera. To examine the noise level and frequencies, we observed the error signal between a fixed ball position and a measurement signal. The results are shown in Figure 3, where the upper part represents a time history of the error signal including noise and the lower part represents the fast Fourier transform (FFT) analysis of the error signal. The noise level in the error signals is relatively high at frequencies over 20 ([rad/s]).

Figure 3: Time history and FFT analysis results of measument signal.

3. I-PD Control by GKYP Lemma

This section describes an I-PD controller design method based on the GKYP lemma. The feedback control system is shown in Figure 4, where a filter is introduced into the control system.

Figure 4: Feedback control system of the ball and plate.

3.1. Low-Pass Filter

In the previous section, we showed that the measurement noise degrades the control performance. To reduce the influence of the noise, a low-pass filter is available in the controller design. According to the noise properties that we observed, it is sufficient to introduce a first-order low-pass filter into the output of the measurement, such that

where  ([rad/s]) is the cut-off frequency.

3.2. I-PD Controller

In the standard PID control, major changes in the reference signals cause large input signals to be generated by the proportional and derivative actions in the controller that saturate the actuator. Indeed, it is difficult to tune the parameters in the PID controller (Figure 4) such that the actuator in this visual feedback system is not saturated. The control system with the I-PD controller (Figure 4) has a structure whose inner loop includes the proportional and derivative actions [7, 8]. In this structure, the integral action alone acts on the error signal and prevents large signals being input to the actuator. The control input can be written as

where is the parameter to approximate the differentiator by a proper transfer function. , , and represent the proportional, integral, and derivative gains, respectively. The open-loop transfer function is

where

and are the integral time and derivative time, respectively. It should be noted that the open-loop transfer function of the I-PD structure is the same as that of the standard PID structure. To tune the parameters in the I-PD controller, we employ an open-loop shaping that realizes a desirable frequency response of the closed-loop system.

3.3. Generalized KYP Lemma

It is known that design specifications for an open-loop transfer function can be reduced to LMIs through the GKYP lemma [9]. We briefly review this lemma in the case of continuous-time systems.

The design specification consists of a frequency range and a desired property in that range. The frequency range can be represented by

where ,

The equality constraint in (13) distinguishes between continuous-time and discrete-time specifications. Since we address continuous-time systems in this paper, we use such that

The inequality constraint in (13) sets a frequency range . For example, a low frequency range is written as

where

Table 2 presents a summary of the choice of versus a type of the frequency range , where , , , and are real positive numbers, and . On the other hand, the desired property in a specific frequency range can be represented by

where

and are the input and output numbers of , respectively. For SISO systems, consider the requirement that in a frequency range is on the half plane under a straight line. That is, is under the straight line, such that

that is equivalent to (18) with

This requirement is designed to reduce sensitivity in a low frequency range. Another requirement is that in a frequency range is in the interior of the circle of radius with the center at . That is, is in the circle such that

which is equivalent to (18) with

This requirement is designed to guarantee robustness in a high frequency range. Under these preparations, the generalized KYP lemma [9] is expressed as follows.

Table 2: Relation between and for continuous time.

(low)

(middle)

(high)

Lemma 1. Let be . in (13) and in (18) are given. Assume that for all . Then the finite frequency condition

holds if and only if there exist Hermitian matrices and such that the matrix inequality condition

is satisfied where

Equation (25) is affine with respect to , , , , and . In the case where and have affine design parameters, (25) is an LMI.

4. Control System Design

4.1. Filter Design

Considering the control performance and the noise level, we set the cut-off frequency in (9) as  ([rad/s]). We examine the effectiveness of this filter in the same experimental setup as given in Section 2.2. Figure 5 shows the spectral analysis results of the measurement and filtered signals, which are represented by the dotted and solid curves, respectively. From these results, it can be seen that the noise at frequencies over 20 ([rad/s]) has been reduced.

Figure 5: Time history and FFT analysis result of measurement signal with low-pass filter.

4.2. State-Space Realization of Open-Loop Transfer Function

To obtain an LMI based on the GKYP lemma, a state-space realization of is required to be affine with respect to a set of the design parameters . If we fix in at , the design parameters appear affinely in the numerator of . Indeed, the controllable canonical form of is written as

Realizations of and are also written as

respectively. By combining these realizations (27)–(29), we obtain a realization of as

where

Consequently, the state-space realization of is affine with respect to .

4.3. Specifications

To shape the Nyquist plot of , we require the following FDI specifications:

Specification (32) with a large ensures sensitivity reduction in the low frequency range by making the gain of high. Specification (33) requires the Nyquist plot to be outside a circle with its center at the point so that a certain stability margin is guaranteed. Specification (34) with a small ensures robustness against the unmodeled dynamics that typically exists in the high frequency range.

In addition to the above basic specifications, we also require the following FDIs that improve the property of trajectory tracking

Since the integral action alone works on the error between the output and the reference signals, the property of trajectory tracking depends directly on the integrator. Here we focus on the corner angular frequency by the integral action in the I-PD controller which is given by , where . We have a strong integral action, and the error is corrected quickly when the corner angular frequency is high, while too high a corner angular frequency causes overshoot and hunting. Thus we impose restrictions for the phase of the I-PD controller. It should be noted that the phase at lower frequencies is about −90 ([deg]) while the phase at the corner angular frequency is about 0 [deg]. If we find the frequency at a specific phase from −90 ([deg]) to 0 ([deg]), the corner angular frequency is greater than that frequency. Specifications (35) and (36) restrict the phase of the I-PD controller as well as the open-loop transfer function so that the corner angular frequency is greater than the frequency at the lowest point in the frequency range 0.3 ([rad/s]).

4.4. I-PD Controller Design

We design an I-PD controller by maximizing subject to Specifications (32)–(36) where

That is, the optimization problem is

Each of the design Specifications (32)–(36) is reduced to an LMI condition through Lemma 1 with the realization (30). The Specification (32) is modified to , where because includes the origin poles that prevent us from taking . Then the LMI optimization problem is to maximize subject to all these LMI conditions where , , and are the common decision variables, while and appear in the LMIs as independent decision variables. It should be noted that appears alone in in the LMI condition (25) corresponding to (36). In this sense, is also an independent decision variable. It should also be noted that and in (25) appear in each of the LMI conditions as the independent decision variables.

To solve this LMI optimization problem, we use YALMIP R20120806 [11], an LMI parser, and SPDT3 version 4.0 [12], an LMI solver, on MATLAB R2011b. The resulting optimal parameters in the I-PD controller and are

The Nyquist plots are shown in Figures 6 and 7 where satisfies the design specifications given in Section 4.3. Since is maximized, sensitivity is reduced in the frequency range. The Bode plots of , , and are shown in Figure 8 where the corner angular frequency of the I-PD controller is larger than 0.3 ([rad/s]).

Figure 6: Nyquist plot.

Figure 7: Nyquist plot (magnification).

Figure 8: Bode plot.

5. Simulation Results

The I-PD controller whose design is described in the previous section was evaluated by a simulation of the step response. To compare the response with that of a standard PID controller, we used the PID controller whose gain parameters are the same as those of the I-PD controller. Since both feedback systems have the same open-loop transfer function, their feedback properties must be the same, provided that each input signal does not saturate. The simulation results of the step response are shown in Figure 9 where the upper and lower parts are the output and input signals, respectively. The solid curves represent the responses given by the I-PD controller while the dashed curves represent those by the PID controller. One can see that the input signal given by the PID controller is saturated, while that by the I-PD controller is not saturated and satisfies the limitation (3). The output signal given by the I-PD controller settles down to the desired value without any overshoot.

Figure 9: Step response simulation results.

It should be noted that the gain parameters in the designed controller are not tuned with the I-PD structure. In our experience, it is difficult not to saturate the input limitation for the PID structure using any design method.

6. Experimental Results

This section evaluates the I-PD controller whose design is given in Section 4 through an experiment of the step response. The PID controller used in Section 4 was also evaluated for comparison. The results of trajectory tracking control by the I-PD controller were also evaluated.

6.1. Step Response Experiment

The results of the step response experiment are shown in Figure 10 where the description of the figure is the same as that of Figure 9. In this experiment, the influence of the time delay appeared and the input signals were slightly larger than those obtained in the simulations. The rise and settling time results are, however, almost the same as those obtained in the simulation.

Figure 10: Step response experimental results.

6.2. Trajectory Tracking Experiments

We tested two kinds of trajectories for tracking control, a square and a circular trajectory. The side of the square trajectory was 0.1 ([m]), and the radius of the circular trajectory was 0.05 ([m]). The results are shown in Figure 11 where the left part shows a trajectory of the square trajectory tracking control, and the right part shows a trajectory of the circular trajectory tracking control. The time histories of the ball and input angles are shown in Figures 12 and 13. In the square trajectory tracking control experiment, the responses were similar to those in the step response experiment except for a slight vibration. Such vibration phenomena are noticeable in the responses of circular trajectory tracking control, in particular, the case when the input signal is relatively small. The reason for these phenomena could be the friction of the ball against the plate or a backlash of the gear system.

Figure 11: Experimental results of tracking control.

Figure 12: Experimental results of square response.

Figure 13: Experimental results of circle response.

7. Conclusions

This paper applied the GKYP lemma to an open-loop transfer function including an I-PD controller and a noise reduction filter for the ball and plate system. The multiple FDI specifications for the finite frequency ranges were satisfied by a solution of the LMI optimization problem. The solution includes the optimal parameters in the I-PD controller. The first-order low-pass filter reduced the noise in the high frequency range and improved the steady-state response. Both simulations and experiments evaluated the effectiveness of the design method by comparing the standard PID controller.

The PI-D (proportional integral-derivative) control system, which moved the derivative controller to the inner feedback loop, also has the same open-loop transfer function as the standard PID controller. Thus the approach in this paper can also be applied to the PI-D controller.

References

M. Moarref, M. Saadat, and G. Vossoughi, “Mechatronic design and position control of a novel ball and plate system,” in Proceedings of the Mediterranean Conference on Control and Automation (MED' 08), pp. 1071–1076, June 2008.

P. Kokotovic, “The joy of feedback: nonlinear and adaptive,” IEEE Control Systems Magazine, vol. 12, no. 3, pp. 7–17, 1992.

F. Zheng, X. Li, X. Qian, and S. Wang, “Modeling and PID neural network research for the ball and plate system,” in Proceedings of the International Conference on Electronics, Communications and Control (ICECC '11), pp. 331–334, September 2011.

X. Fan, N. Zhang, and S. Teng, “Trajectory planning and tracking of ball and plate system using hierarchical fuzzy control scheme,” Fuzzy Sets and Systems, vol. 144, no. 2, pp. 297–312, 2004.

M. Bai, Y. Tian, and Y. Wang, “Decoupled fuzzy sliding mode control to ball and plate system,” in Proceedings of the 2nd International Conference on Intelligent Control and Information Processing (ICICIP '11), pp. 685–690, July 2011.

F. Borrelli, Constrained Optimal Control of Linear and Hybrid Systems, vol. 290 of Lecture Notes in Control and Information Sciences, Springer, 2003.

K. J. Åstrőm and T. Hägglund, Advanced PID Control, Instrumentation Systems, 2005.

Q. Li and Z. Kemin, Introduction to Feedback Control, Prentice Hall, 2009.

T. Iwasaki and S. Hara, “Generalized KYP lemma: unified frequency domain inequalities with design applications,” IEEE Transactions on Automatic Control, vol. 50, no. 1, pp. 41–59, 2005.

A. Rantzer, “On the Kalman-Yakubovich-Popov lemma,” Systems & Control Letters, vol. 28, no. 1, pp. 7–10, 1996.

J. Lőfberg, “YALMIP: a toolbox for modeling and optimization in MATLAB,” in Proceedings of the IEEE International Symposium on Computer-Aided Control System Design (CACSD '11), Taipei, Taiwan, 2004.

K. C. Toh, M. J. Todd, and R. H. Tütüncü, “SDPT3—a MATLAB software package for semidefinite programming, version 1.3,” Optimization Methods and Software, vol. 11, no. 1, pp. 545–581, 1999.