## Servicios Personalizados

## Revista

## Articulo

## Indicadores

- Citado por SciELO

## Links relacionados

- Similares en SciELO

## Compartir

## Latin American applied research

##
*versión impresa* ISSN 0327-0793

### Lat. Am. appl. res. v.35 n.2 Bahía Blanca abr./jun. 2005

**Discrete event control of time-varying plants**

**E. Kofman ^{1}**

^{1} *Laboratorio de Sistemas Dinámicos. FCEIA - UNR - CONICET. Riobamba 245 bis - (2000) Rosario - Argentina Email: kofman@fceia.unr.edu.ar*

*Abstract* ¾ This paper studies the use of Quantized State Control (QSC) in Time-Varying (TV) plants. Making use of a Lyapunov analysis, the stability properties of Time Invariant QSC are extended to the non- stationary case. Then, based on the resulting stability theorem, a design algorithm is developed. Finally, the use of this algorithm -which allows the design of QSC controllers according to stability and convergence speed features- is shown with the design and the simulation of an illustrative example.

*Keywords* ¾ Digital Control. Nonlinear Control. Discrete Event Systems.

**I. INTRODUCTION**

Quantized State Control (Kofman, 2003) is a methodology which allows the digital implementation of previously designed continuous controllers based on the quantization of their state variables. The resulting digital controller -which can be described by a discrete event system- is completed with an asynchronous sampling scheme (Sayiner *et al.*, 1993) so that, ignoring the temporal errors introduced by the clock of the digital device, no time discretization is performed.

The use of QSC instead of classic discrete time digital controllers permits to conserve the region of attraction in nonlinear systems improving also the dynamic response and reducing the computational costs as well as the information transmitted between controllers and plants. In the LTI case, digital QSC controllers also ensure that the resulting trayectories do not differ from the ones obtained with the ideal continuous controllers in more than a bound which can be calculated with a closed formula (Kofman, 2002).

Since QSC takes into account the presence of quantization, it does not guarantee asymptotic stability but *ultimately boundedness* of the solutions (Khalil, 1996). This fact, which seems to be a disadvantage, is common to all the quantization based approaches (Delchamps, 1990; Brockett and Liberzon, 2000; Elia and Mitter, 2001). Moreover, in any sampled data control system where the effects of A/D and D/A converters are taken into account the theory can only talk about ultimate bounds (Miller *et al.*, 1989; Hou *et al.*, 1997).

The original definitions of QSC and the mentioned properties were established under the assumption of a Time-Invariant (TI) plant with a TI controller. Although the restriction on the controller cannot be avoided^{1}, there is no reason to consider only TI plants.

Taking into account the last remark, this work attempts to extend the properties and design algorithms of TI QSC to the Time Varying case.

The paper is organized as follows:

After recalling the principles and properties of QSC (Section II), the main results of QSC with TV-plants -a stability Theorem and a design algorithm- are presented in Section III. Finally, these results are illustrated with the design and simulation of a simple example.

**II. QUANTIZED STATE CONTROL**

QSC is a digital control scheme in which a continuous plant is controlled by a discrete event system obtained with the quantization of a continuous controller.

**A. Quantized State Systems**

Consider the State Equation System (SES) given by:

(1) |

Related to this system, an associated Quantized State System (QSS) (Kofman and Junco, 2001) is given by:

(2) |

In this system, *q*(*t*) (quantized variables) and *x*(*t*) (state variables) are related (componentwise) by hysteretic quantization functions, which are defined as follows:

**Definition 1.** *Hysteretic Quantization Function.*

*Let Q *= {*Q*_{0}*,Q*_{1}, ..., *Q _{r}*}

*be a set of real numbers where Q*

_{k}_{-1}

*< Q*1 £

_{k}with*k*£

*r. Let*Ω

*be the set of piecewise continuous real valued trajectories and let x*Î Ω

_{i}*be a continuous trajectory. Let b*:Ω ® Ω

*be a mapping and let q*(

_{i}= b*x*)

_{i}*where the trajectory q*

_{i}satisfies:(3) |

*and*

*Then, the map b is a* hysteretic quantization function.

The discrete values *Q _{k}* are called

*quantization levels*and the distance

*Q*

_{k}_{+1}-

*Q*is defined as the

_{k}*quantum*, which is usually constant. The width of the hysteresis window is ε and, as it was shown in Kofman

*et al.*(2001), it should be taken equal to the quantum to improve the accuracy. The values

*Q*

_{0}and

*Q*are the lower and upper saturation bounds.

_{r} A fundamental property of a Quantization Function with hysteresis when *t* ³ *t*_{0} and *Q*_{0} £ *x _{i}*(

*t*) £

*Q*is given by the following inequality

_{r}. (4) |

The role of the quantization functions in (2) is to convert the state trajectories (*x*(*t*)) into piecewise constant ones (*q*(*t*)). In that way, provided that *u*(*t*)is piecewise constant, it results that the derivatives (*t*) are also piecewise constant and then the state trajectories are piecewise linear.

As a consequence of these properties of the trajectory forms, QSS can be exactly represented by discrete event models within the DEVS formalism framework (Zeigler *et al.*, 2000). The DEVS model related to a generic QSS and the proof of the mentioned properties can be found in (Kofman and Junco, 2001).

The possibility of representing QSS by DEVS models and the fact that DEVS models can be simulated in real time by digital devices^{2} (Zeigler and Kim, 1993) suggested the use of QSS as digital controllers and the definition of QSC.

**B. QSC Definition**

Consider the Continuous Control System (CCS)

(5) | |

(6) | |

(7) |

consisting of plant (5), controller (6) and their (ideal) interconnection (7).

It is being considered here that the plant could be time varying. When it comes to the controller, it is assumed that it is stationary and it has an input reference *u _{r}*(

*t*).

**Definition 2.** *Quantized State Controller.*

*A QSS associated to a continuous controller* (6) *is called Quantized State Controller (QSC controller). *

**Definition 3.** *QSC System.*

*A QSC system is a control scheme composed by a continuous plant and a QSC controller connected through asynchronous A/D and D/A converters.*

Figure 1. Block Diagram of the QSC system

Figure 1 shows a block diagram representation of a QSC system. Despite its state-equation-like representation, the controller is in fact a DEVS model.

The QSC implementation of the controller transforms (6) into the new set of equations:

(8) |

The asynchronous sampling scheme (Sayiner *et al.*, 1993) implies that the A/D conversions are performed only when the analog input and the digital output of the converters differ in a quantity corresponding to one quantization interval. Then, they can be seen as quantization functions with hysteresis where the quantization intervals and the hysteresis windows have the same size. Similarly, D/A converters can be represented by quantization functions (but without hysteresis). Thus, the presence of the asynchronous converters transforms (7) into:

(9) |

where variables (*t*)and (*t*) are quantized versions of the plant and the controller output variables

**C. QSC Properties**

The CCS closed loop equations can be derived from Eqs.(5)-(7) arriving to

(10) |

Let us define

Δx(_{c}t) = q(_{c}t) - x(_{c}t), | (11a) |

Δy(_{p}t) = (t) - y(_{p}t), | (11b) |

Δy(_{c}t) = (t) - y(_{c}t). | (11c) |

Thus, from these definitions and Eqs.(5), (8) and (9), the QSC closed loop equations can be written as:

(12) |

Then, the QSC system (12) can be seen as a perturbed version of the original CCS (10).

From the property of the quantization functions - provided that the variables *x _{c}*,

*y*and

_{p}*y*do not reach the corresponding saturation bounds

_{c}^{3}- it also results that the perturbation terms are bounded by (4). Thus, the properties of the QSC implementation of a CCS can be studied by looking at the effects of bounded perturbations in the original closed loop system.

Based on this remark, the stability theorem of TI QSC (Kofman, 2003) and the stability and error bound estimation of LTI QSC (Kofman, 2002) were proven.

The stability theorem of TI QSC tells that, given a CCS in which the origin is an asymptotically stable equilibrium point, an appropriate quantization can be found so that the resulting QSC controller ensures ultimately boundedness of the solutions (for any given ultimate bound) conserving also the estimate region of attraction.

That appropriate quantization can be found according to desired ultimate bound following an algorithm which makes use of a Lyapunov function of the original CCS.

On the other hand, the theorem for LTI QSC tells that the QSC implementation of an asymptotically stable LTI CCS is always globally and ultimately bounded. The ultimate bound can be estimated with a closed formula. Moreover, the trajectories of the CCS and the resulting QSC never differ from each other in more than the mentioned bound.

Thus, the quantization choice which completes the QSC design can be made based on these properties, according to the desired ultimate bounds.

**III. TIME VARYING QSC**

The stability theorem and the corresponding design algorithm were based on a stationary Lyapunov analysis. This section introduces the main contributions of this work by extending those results to the time varying case.

**A. Stability of TV QSC**

When the reference trajectory *u _{r}*(

*t*) is zero (or constant), the QSC system (12) can be rewritten as

= *f*(*x* + Δ*x*, Δ*y*, *t*), (13)

where *x* [*x _{p}*,

*x*]

_{c}*, Δ*

^{T}*x*[0, Δ

*x*]

_{c}*, Δ*

^{T}*y*[Δ

*y*, Δ

_{p}*y*]

_{c}*, and*

^{T}*f*[

*f*,

_{p}*f*]

_{c}*. With these definitions, the CCS (10) becomes*

^{T}. (14)

Then, the relationship between the stability properties of (14) and (13) can be stated in following theorem.

**Theorem 1.** *Let the origin be an asymptotically stable equilibrium point of the closed loop CCS *(14)*.As-sume that function f is uniformly continuous and a uniformly continuously differentiable Lyapunov function V *(*x*, *t*) *is known with*

W_{1}(x) £ V(x,t) £ W_{2}(x) | (15) |

(16) |

*t* ³ 0, *x* Î *D being D a compact set which contains the origin and W _{i} are continuous positive definite functions in D. *

*Let = *{*x|W*_{2}(*x*) £ *a*}* with a being an arbitrary positive constant which is enough small so that* {*x|W*_{1}(*x*) £*a*} *is a closed region inside D. *

*Let* = {*x|W*_{1}(*x*) £*b*} *being b an arbitrary positive constant* (*b* < *a*) *enough small so that* Ì .

*Then, a quantization can be found so that the QSC system trajectories starting in finish inside , reaching this region in finite time.*

*Proof.* The derivative of *V*(*x*, *t*) along the solutions of the QSC system (13) is

Then, using (16) it results that

Consider the sets = {*x|W*_{2}(*x*) £ *b*} and = {*x|W*_{1}(*x*) £ *a*}. From the hypotesis made about *a *and *b*, it results that

Ì Ì Ì *D*. (17)

Since *W*_{3}(*x*) is positive definite in *D*, it is positive in Ω_{1,2} - . Moreover, it exists a positive constant *s* where

. (18) |

Let us define the following function

(19) |

*∂* The continuity in functions *W*_{3} and *f* and the fact that *V* is continuously differentiable implies that α is uniformly continuous. From the definition of α and Eq.(18) it results that

α(*x*, 0, 0, *t*) £ -*s* *x* ÎΩ_{1,2}, *t* ³ 0. (20)

Let α* _{M}* be the function defined by

. (21) |

It can be easily seen that *αM *is continuous and α* _{M}* (0, 0) £

*-s*. Then, for any positive number

*s*

_{1}<

*s*a positive constant

*r*can be found so that the condition

||(Δ*x*, Δ*y*)||£ *r* (22a)

implies that

α* _{M}* (Δ

*x*, Δ

*y*) £

*-s*

_{1}(22b)

and then it results that

(23)

in Ω_{1,2}.

From here to the end, the proof follows Theorem 5.3 of (Khalil, 1996).

Let Ω* _{a,t}* = {

*x|V*(

*x*,

*t*) £

*a*} and Ω

*= {*

_{b,t}*x|V*(

*x*,

*t*) £

*b*} be time variable sets. From (15) it results that

Ì Ω* _{b,t}* Ì Ì Ì Ω

*Ì Ì*

_{a,t}*D*(24)

The border of Ω* _{a,t}* is inside Ω

_{1,2}, where (

*x*,

*t*) is negative. It means that the trajectories of the QSC system (13) cannot abandon Ω

*.*

_{a,t}Then, any trajectory starting in cannot abandon .

The border of Ω* _{b,t}* is also inside Ω

_{1,2}. Then the trajectories cannot abandon this time variable set.

To complete the proof, we need to ensure that the trajectories initiated in Ì Ω* _{a,t}* reach Ω

*in a finite time.*

_{b,t} Let f(*t*) be a solution of (13) starting in Ω* _{a,t}* (i.e.

*V*(f(0), 0) £

*a*) and let us suppose that

*V*(f(*t*), *t*) > *b* *t* (25)

then we have (f(*t*), *t*) £ *-s*_{1} and after

(26) |

it results that *V*(f(*t*_{1}), *t*_{1}) £ *b* which yields a contradiction. Then, the region Ω* _{b,t}* must be reached before the finite time

*t*

_{1}.

Since Ω* _{b,t}* Ì the trajectory also reaches the set before that time.

Equation (22a) gives the maximum perturbation allowed to ensure the achievement of the proposed goal (i.e. region of attraction and ultimate bound in ). Since the maximum perturbation in each variable is given by the corresponding quantum, this equation should be used to choose the quantum at the different controller state variables and converters completing in that way the QSC design.

Observe that is also the estimation of the region of attraction of the CCS using the Lyapunov function*V*. Then, a QSC implementation can be found so that it conserves the estimated region of attraction.

**B. Design Algorithm for TV QSC **

The design of a QSC contoller can be divided in two steps. The first one is the design of the continuous controller, which can be done following any technique.

The second step is the choice of the quantization at each variable. The use of a very small quantum yields solutions which are very closed to the trajectories of the CCS. This is due to the property of convergence which tells that -under locally Lipchitz conditions on the functions- the solutions of the QSC system go to the solutions of the CCS when the quantization go to zero (Kofman, 2003). In that way, and also according to Theorem 1, the ultimate bound can be reduced to arbitrary small values.

However, the use of a small quantum increases the number of events at the controller and the digital device can fail in its attempt to give the correct output values at the required time.

Therefore, there is always a trade-off between accuracy and practical considerations related to the computational costs. Then the idea is to exploit Theorem 1 in order to choose the quantization according to some essential features (region of attraction and ultimate bound). In that way, the quantization adopted will be just as small as necessary to ensure those properties and -provided that the CCS is not too fast- the digital device will be able to correctly implement the resulting QSC controller.

The translation of these ideas into a design algorithm for QSC can be written as follows:

1. Design a continuous controller and calculate the Lyapunov function *V*(*x*,*t*) and the functions *W _{i}*(

*x*) according to (15)-(16) for the closed loop CCS.

2. Choose the QSC region of attraction and the ultimate region together with the constants

*a*and

*b*.

3. Obtain the perturbed closed loop function

*f*according to (13).

4. Obtain function α according to (19) and α

*according to (21).*

_{M}5. Calculate constant

*s*according to (18) and choose the positive constant

*s*

_{1}. If the goal is just to ensure ultimately boundedness

*s*

_{1}should be very small. Otherwise, if the speed of convergence is also important, it can be chosen taking into account (26).

6. Compute the value of

*r*according to (22).

7. Choose the quantization at the controller state variables and converters so that (22a) is satisfied.

It can be easily seen that this procedure leads to a QSC controller which ensures region of attraction and ultimate region .

**IV. EXAMPLES AND RESULTS **

The unstable time varying plant

can be stabilized by the controller

(27) |

The resulting closed loop system can be written as

Here, the Lyapunov candidate

verifies (15) with

(28) |

The orbital derivative is

which satisfies (16) with

Then, the closed loop CCS is asymptotically stable and the algorithm resulting from Theorem 1 can be used to design the QSC controller.

The first step for the QSC design consists in choosing the region of attraction and the ultimate bound.

Since inequalities (15)-(16) stand in Â^{2}(i.e. the CCS stability is global) it is not necessary to restrict the region of attraction except for chosing the saturation values. In this case, the choice of does not affect the calculations.

The ultimate bound will be chosen with *b *=0*.*5. Then, taking into account (28) it results that = {*x* Î Â^{2}| ||*x*|| £ 1}.

The perturbed equations (13) can be written as

and then, from (19) function α(*x*, Δ*x*, Δ*y*,*t*) results

Although the maximum α* _{M}* in (21) cannot be easily obtained, it can be bounded using the fact that ||

*x*|| > 1/3 in Ω

_{1,2}. Then, after some calculations it results that

Thus, taking the quantum equal to 0.018 in the controller and converters we can ensure that α* _{M}* £ -0.0093 in Ω

_{1,2}, which implies that the trajectories finish inside in finite time, with a minimum

*speed s*

_{1}= 0.0093.

Figures 2-5 show the simulation results for an initial condition *x _{p}* = 5,

*x*= 0.

_{c}

Figure 2. *x _{p}* vs.

*t*

Figure 5 also shows that the design was very conservative. The ultimate bound observed in the simulation is less than 0.03 (in norm 2) which is more than 30 times smaller than the estimated.

There are two reasons which can explain this. The first one is that Lyapunov analysis often leads to conservative estimations of the ultimate bounds. The second reason is that in this case the D/A converter does not introduce any perturbation since it is exactly matched with the quantizer of *x _{c}* and the A/D converter.

Figure 3. Final oscillations in *x _{p}* vs.

*t*

Figure 4. *x _{p}* vs.

*x*

_{c}

Figure 5. Final oscillations in *x _{p}* vs.

*x*

_{c} **V. CONCLUSIONS**

This work studied the properties of QSC in time varying systems showing that the main features and advantages of the methodology are also satisfied in these cases.

Theorem 1 proved that the QSC control of time varying plants allows to ensure ultimately bounded-ness conserving the region of attraction.

A design algorithm derived from this theorem was also introduced and its use was illustrated with an example.

When it comes to future work, it should be considered that the main theoretical properties were already proven for general nonlinear time varying systems and for LTI systems. There is an intermediate case which should be taken into account which corresponds to Linear Parameter-Varying (LPV) plants. The example given in fact corresponds to that category and although a Lyapunov analysis could be done for that case, the result was very conservative. If the geometrical analysis for LTI systems of (Kofman, 2002) were extended for LPV plants, less conservative results might be obtained.

^{1} In the general case it is not possible to obtain a discrete event description of a quantized time-varying system.

^{2} DEVS representation of QSS is exact. However, real time simulation of DEVS has errors related to temporal resolution and round-off in digital devices

^{3} It will be considered that the non-saturation region is enough large so that the trajectories do not leave it.

**REFERENCES**

1. Brockett, R. and D. Liberzon, "Quantized feedback stabilization of linear systems", *IEEE Trans. Automat. Contr.*, **45**, 1279-1289 (2000). [ Links ]

2. Delchamps, D., "Stabilizing a linear system with quantized state feedback", *IEEE Trans. Automat. Contr.*, **35**, 916-924 (1990). [ Links ]

3. Elia, N. and S. Mitter, "Stabilization of Linear Systems with Limited Information", *IEEE Trans. Automat. Contr.*, **46**(9), 1384-1400 (2001). [ Links ]

4. Hou, L., A. Michel, and H. Ye, "Some qualitative properties of sampled-data control systems", *IEEE Trans. Automat. Contr.*, **42**, 1721-1725 (1997). [ Links ]

5. Khalil, H., *Nonlinear Systems*, Prentice-Hall, New Jersey, 2nd edition (1996). [ Links ]

6. Kofman, E., "Quantized-State Control of Linear Systems",*Proceedings of AADECA 2002*, Buenos Aires, Argentina (2002). [ Links ]

7. Kofman, E., "Quantized-State Control. A Method for Discrete Event Control of Continuous Systems", *Latin American Applied Research*, **33**(4), 399-406 (2003). [ Links ]

8. Kofman, E. and S. Junco, "Quantized State Systems. A DEVS Approach for Continuous System Simulation", *Transactions of SCS*, **18**(3), 123-132 (2001). [ Links ]

9. Kofman, E., J. Lee, and B. Zeigler, "DEVS Representation of Differential Equation Systems. Review of Recent Advances", *Proceedings of ESS'01*, 591-595 (2001). [ Links ]

10. Miller, R., A. Michel, and J. Farrel, "Quantizer effects on steady-state error specifications of digital feedback control systems", *IEEE Trans. Automat. Contr.*, **34**, 651-654 (1989). [ Links ]

11. Sayiner, N., H. Sorensen, and T. Viswanathan, "A non-uniform sampling technique for A/D conversion", *Proceedings of ISCAS'93*, **2**, 1220-1223 (1993). [ Links ]

12. Zeigler, B. and J. Kim, "Extending the DEVS-Scheme knowledge-based simulation environment for real time event-based control", *IEEE Trans. on Robotics and Automation*, **3**(9), 351-356 (1993). [ Links ]

13. Zeigler, B., T. Kim, and H. Praehofer, *Theory of Modeling and Simulation. Second edition*, Academic Press, New York (2000). [ Links ]