`Mathematical Problems in EngineeringVolume 2010, Article ID 842896, 20 pageshttp://dx.doi.org/10.1155/2010/842896`
Research Article

## Dynamic Tracking with Zero Variation and Disturbance Rejection Applied to Discrete-Time Systems

1Department of Electrical Engineering, Universidade Estadual Paulista (UNESP), Campus of Ilha Solteira, 15385-000 Ilha Solteira, Brazil
2Academic Department of Electronics, Federal Technological University of Paraná (UTFPR), 80230910 Curitiba, PR, Brazil

Received 1 December 2009; Revised 29 June 2010; Accepted 9 December 2010

Copyright © 2010 Renato de Aguiar Teixeira Mendes 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

The problem of signal tracking in discrete linear time invariant systems, in the presence of a disturbance signal in the plant, is solved using a new zero-variation methodology. A discrete-time dynamic output feedback controller is designed in order to minimize the norm between the exogen input and the output signal of the system, such that the effect of the disturbance is attenuated. Then, the zeros modification is used to minimize the norm from the reference input signal to the error signal. The error is taken as the difference between the reference and the output signal. The proposed design is formulated in linear matrix inequalities (LMIs) framework, such that the optimal solution of the stated problem is obtained. The method can be applied to plants with delay. The control of a delayed system illustrates the effectiveness of the proposed method.

#### 1. Introduction

In a control systems theory, the design of controller using pole placement of closed loop discrete-time systems can be easily done. In [1] a controller using pole placement is used to obtain an exact plot of complementary root locus, of biproper open-loop transfer functions, using only well-known root locus rules. However, the problem of zero placement is not very much studied by the control researchers. In [2] a discrete-time pole placement is obtained by a control design technique that uses simple and multirate sample. The methodology proposed in [3] preserves the state feedback controller optimality by pole placement in a plain region specified in design. In the field of discrete-time systems pole placement we find [4], where discrete adaptative controllers are designed considering arbitrary zero location. Also, in [5] a class of nonmodeled dynamics is controlled using a zero placement.

In [6] a methodology is proposed using zero and pole placement for discrete-time systems, to obtain the signal tracking and disturbance rejection, respectively. However, for the signal tracking problem, when a state feedback estimator is proposed, a modification occurs in -norm value obtained with the initial controller that provides the disturbance rejection. The methodology proposed in this paper has the advantage of maintaining the -norm value obtained with the initial controller for the signal tracking problem.

The problem of signal tracking, in the presence of disturbance signal for continuous-time plant, was solved in [7], using a zero variation methodology. A methodology with a simpler mathematic formulation is proposed in [8]. The signal tracking problem with disturbance rejection in discrete-time systems is solved by an analytic method in [9], however, the mathematic formulation is complex and a frequency selective tracking is not presented as proposed in this manuscript. In [10] the use of linear matrix inequalities (LMIs) is considered in design of controllers, filters and stability study. Also, in [11], LMIs are used for the design of a dynamic output feedback controller in order to guarantee the asymptotic stability of a continuous-time system and minimize the upper bound of a given quadratic cost function. Furthermore, the LMI formulation has been used in several engineering problems (see, e.g., [1220]).

This manuscript proposes a formulation of a signal tracking with disturbance rejection optimization problem for discrete-time systems in the linear matrix inequalities framework, such that the optimal solution of the stated control problem is obtained. The proposed method is simpler than the other tracking techniques, and the main result is that when the problem is feasible the optimal solution is obtained with small computation effort, as the LMIs can be solved using linear programming algorithms, with polynomial convergence. The software MATLAB [21] is used to find the LMI solutions, when the problem is feasible. The control of a delayed system illustrates the effectiveness of the proposed method.

#### 2. Statement of the Problems

Consider a controllable and observable linear time-invariant multi-input multi-output (MIMO) discrete-time system, where , , , , , is the state vector, is the output vector, is the control input and is the disturbance input (exogenous input).

Problem 1. The disturbance rejection problem for discrete-time systems, using the dynamic output feedback of the system described in (2.1), is the following: minimize the upper bound of norm from the exogenous input to the output . In this context, the objective is to design a controller , , that attenuates the effect of disturbance signal in the output of the system. And in the tracking process, it is needful to design the controllers and that minimize the norm between the reference input and the tracking error .

Remark 2.1. The problems of weighted disturbance reduction and weighted reference tracking are closely related to the above one and will be addressed in Section 5.

Remark 2.2. The block diagram of the control process used in this manuscript to solve Problem 1 is given in Figure 1, where is a dynamic output feedback controller, the controllers and were used to solve the zero variation problem in order to obtain a tracking system, is the reference input signal.

Figure 1: Discrete-time optimal control system with pole placement and zero modification.

The state space equation of the control system shown in Figure 1 can be written as: Rewriting the system (2.2) in a compact form, it follows that: where Using the -transform in order to solve the system (2.3), consider the initial conditions equal to zero. One obtains the transfer function between input signals reference input and exogenous input and the measured output of the system as showed in the following equation: For the transfer function from to in (2.6), the minimization of norm is obtained with the initial design of controller, that implies in the minimization of the perturbation effect in to the system output.

Figure 1 shows the addition of the term in the structure of the controller. The purpose of the controller is only to change the zeros of the transfer function from to and it does not change the poles obtained in the initial design of . The transfer function from to is not changed by or , according to (2.4) and (2.5). In this way the performance of the norm controller is not affected.

For the optimal tracking design, the relation between error signal and reference signal described in (2.7) is considered, making the perturbation signal equal to zero in (2.6), In this case, using the zero modification one can design a tracking system that minimizes the norm between the reference input and the tracking error . In Section 4, motivated by the work in [22], we show that and modify the zeros from to . The process of the zeros modification does not interfere in the disturbance rejection. Therefore in agreement with (2.6), has no influence on the transfer function from to . In (2.6) one uses the zeros location, by the specifications of the and in , in the process of minimization of the norm of the transfer function between the reference signal and the tracking error.

#### 3. Dynamic Output Feedback Controller Design

The following theorem leads to a new method to design the in a LMI framework, and the goal is to attenuate the effects of exogenous signal in the output of discrete-time systems. By using [23], a pole placement constraint region with radius and center in is required and used in this work to provide the designer with an expedite way to keep the controller gains within appropriate bounds.

Theorem 3.1. Consider the system (2.1) with dynamic output feedback by the upper bound controller, . Then the optimal solution of the norm between the input and the output , with pole placement in a region of radius and center in can be obtained from the solution of the following LMI optimization problem, where, , , , , and in (3.1), (3.2), and (3.3) are the set of LMIs optimization variables. The radius and center in are the pole placement constraints illustred in Figure 2
For solution of (3.4): , where and can be obtained by L-U decomposition of [24].
The compensator dinamic matrices is obtained with the solution of equation (3.4).
In (3.1), (3.2) and (3.3), one has the following change of variables:

Figure 2: The pole placement region in the plane.

Proof. A realization of dynamic output feedback system is as followings: where The optimization problem below described in LMI framework [25], is used to design the compensator with pole placement constraints However, a high computational effort is needed to solve this problem, because the optimization problems (3.7) and (3.8) are described as a solution of BMIs. Then, using a linear transformation, the problem can be easily solved, based on LMI framework. First, the matrix and its inverse are considered as follows where , , and The condition (3.3) is obtained by considering the Lyapunov matrix and premultiplying and postmultiplying by and , respectively. After, the condition (3.1) is obtained premultiplying and postmultiplying the inequation (3.7) by (3.11) and (3.12), respectively. Then, pre and postmultiplying (3.8) by (3.13) and (3.14), respectively, results in condition (3.2) The zero modification is showed in the next section.

Remark 3.2. In this paper, the methodology adopted to solve the pole location problem affords the designer an expedite way to keep the controller gains within appropriate bounds, a key requisite for implementation purposes.

#### 4. Zeros Variation

Motivated by the work in [22] the present paper uses the zero variation in order to obtain the global optimum of norm to solve the tracking problem.

A controller that modifies the zeros from to is designed considering the system as shown in Figure 3. The zeros of closed-loop system are allocated at arbitrary places according to the and values, where and .

Figure 3: Zeros variation of the controlled system.

The plant is described by The transform of (4.1) with zero initial condition, is A zero of the system is a value of such that the system output is zero even with a nonzero state-and-input combination. Thus if we are able to find a nontrivial solution for and such that is null, then is a zero of the system [22]. Combining the two parts of (4.2) we must satisfy the following requirement: Also, the compensator can be described as follows where .

A more general method to introduce is to add a term to and also a term to the control equation , as shown in Figure 3. The controller, with these additions, becomes equal to: Now, considering the controller (4.5), if there exists a transmission zero from to , then necessarily there exists a transmission zero from to , unless a pole of the plant cancels the zero. The equation to obtain from to (we let because we are considering only the effects of , then ) in (4.5), is the following: Because the coefficient matrix in (4.6) is square, the condition for a nontrivial solution is that the determinant of this matrix must be zero. Thus we have Multiplying the second column of the matrix described in (4.7) at right by a nonzero matrix and then adding to the first column of (4.7) the product of by the last column, we have: And so, considering , where the modified zeros from to are the solutions . It is important to notice that the gain and the vector do not only modify the system zeros but also are used to obtain the optimal solution of the tracking problem.

#### 5. Tracking Design

The solution of the tracking problem is based on the design of the matrices of the controller and , that minimize the norm of . Weighted frequency is added in the tracking system in order to track signals in a frequency band specified in the design. In tracking design with wheigthed frequency, the goal is to find a global solution that optimize the problem described as follows: where is a dynamic system designed to specify wheighted frequency in the output. A stable, linear and time invariant system realization is considered as indicated in (2.7). Figure 4 illustrates the structure of inclusion of frequency wheighted in the design of tracking system.

Figure 4: System structure to the tracking with frequency wheighted.

The system (2.4) can be represented by state variables in function of and , as follows: In addition, a possible state space realization of is: A metodology for the tracking design problem solution with wheigthed band is proposed in Theorem 5.1 considering the elements of the compensator matrix already fixed.

Theorem 5.1. Considering the system with output filter (5.3), if there exist a solution for the LMI described in (5.4) and (5.5), then the gain and vector that minimize the norm from to can be obtained

Proof. For the tracking design with band weighted, substitute , , , and in (5.3) by (3.7). It results in the optimization problem described in (5.4) and (5.5). The gain and vector are obtained from this process. The matrix is partitioned in the form , .

The theorem proposes that one can obtain an optimal zero location that guarantee the global optimization of the tracking error norm. The proposed design is formulated in linear matrix inequalities (LMIs) framework, such that the optimal solution is obtained.

The controller and are the optimal solution of (5.4) and (5.5) and minimize the norm between the reference input signal and the tracking error signal . The filter has the goal of adjusting the controllers and for the selected frequency band. Then the zero variation described in LMI framework considers the filter dynamics to adjust the operation of the system tracking in this frequency band. However, in the implementation or simulation of the tracking process control, the filter is discarded.

The paper also proposes the statement of weighted disturbance minimization problem that is achieved based on Theorems 3.1 and 5.1 demonstrations and uses the following formulation: The following examples illustrate the effectiveness of the proposed method.

#### 6. Example 1

Consider a tank temperature control system described in [22], with a zero-order hold. The goal is to design a tracking system for flow control with disturbance attenuation. The model of the delayed system [22] is, A sampling period of 0.01 seconds is used in the design. The parameters and are adopted in the design.

We found , and therefore, and .

Using the process to find the transform of a delayed continuous-time function (6.1), the parameters are substituted and we obtain or, The state-space description of the system is where is the state vector, is the control signal and is a disturbance signal in the system.

The design of the tracking system must include operation for reference signals of low frequencies (smaller than 0.1 rad/s). In such a case, the following filter was considered Using Theorem 3.1 the controller is designed for the system described in (6.4) and showed in (6.6). This controller minimizes the -norm of to . In this design we obtain a disk of radius and center in is used as a pole placement constraint: The -norm of to for the closed-loop system is 0.0342, implying the attenuation of the effect of the disturbance signal in the system. Figure 5 illustrates the frequency response of .

Figure 5: Frequency response of .

To design a tracking system, the proposed zero-variation methodology given by (5.4) and (5.5) were used and the norm of to was minimized considering the signal with low frequencies (smaller than 0.1 rad/s). The obtained -norm for all frequency spectra was equal to 1.32, while for frequency band specified in the problem, the largest magnitude of frequency response was 3.13 × 10−3). This implies that the tracking system operated adequately in the frequency band specified in the problem.

Figure 6 illustrates the frequency response of and one can verify that the -norm in frequency band follows the characteristcs of a tracking system. The and optimum controller values were In the simulation, an unit step input is considered. The simulation result is illustrated in Figure 7.

Figure 6: Frequency response of .
Figure 7: Response of the output to a unit step input .

In the second simulation we use the following ramp signal: with . The simulation results are illustrated in Figure 8.

Figure 8: Output of the system for ramp input.

Finally, an input signal and a disturbance signal with random amplitudes were simulated. The maximum amplitude of the random signal was equal to 1. The simulation results are illustrated in Figure 9.

Figure 9: Output signal and input signal are almost overlapped.

For this example, the zeros of the system where −1; and . The poles of the system with a feedback controller where ; ; and . It is possible to see that the poles of the system with a feedback controller were allocated according to the circle constraint specification. The pole-zero map is shown in Figure 10.

Figure 10: Pole-zero map of the closed-loop system obtained with , and .

The example above shows the methodology effectiveness. The disturbance rejection and the minimization of the tracking error for specified frequency band were reached. It was showed that the methodology works properly for ramp, unit step and sinusoidal signals for any frequency in the specified frequency band.

#### 7. Example 2

Consider a discrete form of a continuous model plant that has a zero in the right-half plane. A sampling period of 0.01 seconds is used in design. The state-space description for the system is as follows:

Continuous Model:

Discrete Model:
where is the state vector, is the control signal and is a disturbance signal in the system.

The design of the tracking system must include operation for reference signals of low frequencies (down to 5 rad/s). In such a case, the filter was considered Using Theorem 3.1 the controller is designed for the system described in (7.2) and shown in (7.4). This controller minimizes the -norm of to . In this design a disk of radius 0.85 and center in −0.1 is used as a pole placement constraint The -norm of to for the closed-loop system was 1.577  10−3, implying the attenuation of the effect of the disturbance signal in the system. Figure 11 illustrates the frequency response of .

Figure 11: Frequency response of .

To design a tracking system, the proposed zero-variation methodology (5.4) was used in which the norm of to is minimized considering the signal with low frequencies (down to 5 rad/s). The obtained -norm for all frequency spectra was equal to 1.8, while for frequency band specified in the problem, the largest magnitude of frequency response was 0.031. This implies that the tracking system operated adequately in the frequency band specified in the problem.

Figure 12 illustrates the frequency response of and one can verify that the -norm in frequency band follows the characteristcs of a tracking system. The and optimum parameters values were Then, an input signal and a disturbance signal with random amplitudes were simulated, and it was found that the maximum amplitude of the random signal was equal to 1. The simulation results are illustrated in Figure 13.

Figure 12: Frequency response of .
Figure 13: Output signal and input signal are almost overlapped.

#### 8. Conclusion

In this manuscript, it is proposed a methodology to solve the tracking and disturbance rejection problem applied to discrete-time systems. Considering Figure 1 the disturbance signal acting in the plant can be attenuated by minimizing the -norm from to , by using a dynamic feedback compensation. In the tracking process, a zero-variation methodology is used in order to minimize the -norm between the reference signal and tracking error signal, where the tracking error is the diference between the reference signal and system output signal . In the tracking design with disturbance rejection, the pole placement is used to attenuate the disturbance signal effect, while the zero variation allows the tracking. The zero modification do not interfere in the design of the disturbance rejection. In the tracking process, the frequency band wheighted allows to choose the frequency band on the reference input signal. The tracking method and disturbance rejection are based on LMI framework. Then, when there exists a feasible solution the design can be obtained by convergence polynomial algorithms [23, 25] available in the literature.

#### Acknowledgments

The authors gratefully acknowledge the partial financial support by FAPESP, CAPES and CNPQ of Brazil.

#### References

1. M. C. M. Teixeira, E. Assunção, R. Cardim, N. A. P. da Silva, and E. R. M. D. Machado, “On complementary root locus of biproper transfer functions,” Mathematical Problems in Engineering, vol. 2009, Article ID 727908, 14 pages, 2009.
2. M. De la Sen, “Pole-placement in discrete systems by using single and multirate sampling,” Journal of the Franklin Institute B, vol. 333, no. 5, pp. 721–746, 1996.
3. A. Saberi, P. Sannuti, and A. A. Stoorvogel, “${H}_{2}$ optimal controllers with measurement feedback for discrete-time systems: flexibility in closed-loop pole placement,” Automatica, vol. 33, no. 3, pp. 289–304, 1997.
4. M. M'Saad, R. Ortega, and I. D. Landau, “Adaptive controllers for discrete-time systems with arbitrary zeros: an overview,” Automatica, vol. 21, no. 4, pp. 413–423, 1985.
5. W. C. Messner and C. J. Kempf, “Zero placement for designing discrete-time repetitive controllers,” Control Engineering Practice, vol. 4, no. 4, pp. 563–569, 1996.
6. R. A. T. Mendes, Controle Ótimo com Modificação de Zeros para o Problema do Rastreamento em Sistemas Discretos usando LMI, M.S. thesis, Unesp, São Paulo, Brazil, 2007.
7. C. Q. Andrea, Controle Ótimo 2 e com Alocação de Zeros para o Problema de Rastreamento usando LMI, M.S. thesis, UNESP, São Paulo, Brazil, 2002.
8. E. Assunção, C. Q. Andrea, and M. C. M. Teixeira, “2 and -optimal control for the tracking problem with zero variation,” IET Control Theory and Applications, vol. 1, no. 3, pp. 682–688, 2007.
9. B. M. Chen, Z. Lin, and K. Liu, “Robust and perfect tracking of discrete-time systems,” Automatica, vol. 38, no. 2, pp. 293–299, 2002.
10. M. C. Oliveira, Controle de Sistemas Lineares Baseado nas Desigualdades Matriciais Lineares, Ph.D. thesis, Unicamp, Campinas, Brazil, 1999.
11. J. H. Park, “On design of dynamic output feedback controller for GCS of large-scale systems with delays in interconnections: LMI optimization approach,” Applied Mathematics and Computation, vol. 161, no. 2, pp. 423–432, 2005.
12. E. Assunção and P. L. D. Peres, “A global optimization approach for the 2-norm model reduction problem,” in Proceedings of the 38th IEEE Conference on Decision and Control (CDC '99), pp. 1857–1862, Phoenix, Ariz, USA, December 1999.
13. E. Assunção, M. C. M. Teixeira, F. A. Faria, N. A. P. da Silva, and R. Cardim, “Robust state-derivative feedback LMI-based designs for multivariable linear systems,” International Journal of Control, vol. 80, no. 8, pp. 1260–1270, 2007.
14. M. C. M. Teixeira, E. Assunção, and R. G. Avellar, “On relaxed LMI-based designs for fuzzy regulators and fuzzy observers,” IEEE Transactions on Fuzzy Systems, vol. 11, no. 5, pp. 613–623, 2003.
15. M. C. M. Teixeira, E. Assunção, and R. M. Palhares, “Discussion on: “ output feedback control design for uncertain fuzzy systems with multiple time scales: an LMI approach”,” European Journal of Control, vol. 11, no. 2, pp. 167–169, 2005.
16. M. C. M. Teixeira, E. Assunção, and H. C. Pietrobom, “On relaxed LMI-based design for fuzzy regulators and fuzzy observers,” in Proceedings of the 6th European Control Conference, pp. 120–125, Porto, Portugal, 2001.
17. E. Teixeira, E. Assunção, and R. G. Avelar, “Design of SPR systems with dynamic compensators and output variable structure control,” in Proceedings of the International Workshop on Variable Structure Systems, vol. 1, pp. 328–333, Alghero, Italy, 2006.
18. F. A. Faria, E. Assunção, M. C. M. Teixeira, R. Cardim, and N. A. P. da Silva, “Robust state-derivative pole placement LMI-based designs for linear systems,” International Journal of Control, vol. 82, no. 1, pp. 1–12, 2009.
19. R. Cardim, M. C. M. Teixeira, E. Assunção, and M. R. Covacic, “Variable-structure control design of switched systems with an application to a DC-DC power converter,” IEEE Transactions on Industrial Electronics, vol. 56, no. 9, pp. 3505–3513, 2009.
20. F. A. Faria, E. Assunção, M. C. M. Teixeira, and R. Cardim, “Robust state-derivative feedback LMI-based designs for linear descriptor systems,” Mathematical Problems in Engineering, vol. 2010, Article ID 927362, 15 pages, 2010.
21. P. Gahinet, A. Nemirovsk, A. J. Laub, and M. Chiali, LMI Control Toolbox User's Guide, The Mathworks Inc., Natick, Mass, USA, 1995.
22. G. F. Franklin, J. D. Powell, and M. L. Workman, Digital Control of Dynamic Systems, Addison Wesley, New York, NY, USA, 2nd edition, 1990.
23. M. Chilali and P. Gahinet, “ design with pole placement constraints: an LMI approach,” IEEE Transactions on Automatic Control, vol. 41, no. 3, pp. 358–367, 1996.
24. E. L. Lima, Albebra Linear, Coleção Mateméatica Universitéaria, IMPA, Rio de Janeiro, Brazil, 4th edition, 2000.
25. R. M. Palhares, R. H. C. Takahashi, and P. L. D. Peres, “ and Ótimo 2 guaranteed costs computation for uncertain linear systems,” International Journal of Systems Science, vol. 28, no. 2, pp. 183–188, 1997.