Comparison of rational approximations of (49) measured by H_{∞} and *H*_{2} norms

## 1. Introduction

Linear time-invariant time delay systems (LTI TDS) have usually been assumed to contain delay elements in input-output relations only. All the system dynamics has been hence modeled by point accumulations in the form of a set of ordinary differential equations. The Laplace transform then results in a transfer function expressed by a serial combination of a delayless term and a delay element. However, this conception is somewhat restrictive in effort to fit the real plant dynamics because inner feedbacks are often of the time-distributed or delayed nature.

*Anisochronic* (or *hereditary*) TDS models, on the other hand, offer a more universal dynamics description applying both integrators and delay elements either in lumped or distributed form so that delays appear on the left side of a differential equation which is no longer *ordinary* (ODE) but rather *functional* (FDE) - this brings the concept of *internal* (or *state*) delays. In contrast to undelayed systems, the main difference in dynamics is that their spectra are infinite in general. In the further text, an abbreviation TDS means LTI TDS containing state delays with or without input-output delays.

Already in (Volterra, 1928) differential equations incorporating the past states when studying predator-pray models were formulated. The theory of these models has been then developed by many outstanding authors, see e.g.(Bellman & Cooke, 1963), (Krasovskii, 1963), (Kolmanovskii & Nosov, 1986), (Zítek, 1983), (Górecki et al., 1989), and especially (Hale & Verduyn Lunel, 1993) and (Nicolescu, 2001), to name a few. Aftereffect phenomenon is included in many processes, e.g. in chemical processes (Zítek & Hlava, 2001), heat exchange networks (Zítek, 1997), in models of mass flow in sugar factory (Findeisen et al., 1970), in metallurgic processes (Morávka & Michálek, 2008), etc. Plenty of references to examples of processes with internal delays, covering a wide range of human activities (e.g. biology, chemistry, economics, communication and information technologies, etc.) are introduced in (Kolmanovskii& Myshkis, 1999; Niculescu 2001; Richard, 2003). Capabilities and advantages of this class of models and controllers for modeling and process control were broadly discussed in (Manitius & Olbrot, 1979). TDS models can be used not only for description of those systems embodying internal delays but they are successfully capable to fit the dynamics of high-order systems and processes even without apparent delays (Zítek & Víteček, 1999; Vyhlídal & Zítek, 2001; Simuenovic, 2011), which simplifies the processes description.

Using the Laplace transform applied to FDEs of TDS, input-output models in the single-input single-output (SISO) case can be expressed by the transfer function in the form of a fraction of so-called quasipolynomials (El’sgol’ts&Norkin, 1973) which can be viewed as polynomials in complex variable *s* over the ring of exponentials polynomials (or over a ring of linear combinations of real numbers and exponential functions in *s*). The concept of pseudopolynomials (Brethé& Loiseau, 1998), meromorphic functions (Zítek & Kučera, 2003)or a special ring (Gluesing-Lueerssen, 1997; Conte & Perdon, 2000), instead of quasipolynomials, can be used as alternatives; however, from the description point of view all these models are equivalent – their significance emerges while controller design.

Delay in the feedback control significantly deteriorates control performance, namely stability and periodicity, and makes controller design more complex – mainly if delays are not approximated or ignored. Therefore, design a suitable control law for such systems is a challenging task solved by various techniques and approaches; a plentiful enumeration of them can be found e.g. in (Richard, 2003). Due to the complexity of many direct methods, one can approximate an infinite-dimensional model by a finite-dimensional one to obtain a rational transfer function which can be handled by any standard control approach for undelayed plants. Similarly, in case of delayed (anisochronic) controllers obtained by controller design without delay approximation, one way how to implement thesecontrol laws on discrete-time-working machines such as PLC or PC is to find a delayless (i.e. finite-dimensional) approximating model followed by a discretization (e.g. by the z-transform with a holder).

This chapter aims two problems. First, output controller design for TDS models in a special ring is introduced. Second, resulting anisochronic controllers are approximated using various techniques for rational approximation of transfer function of TDS, which is the crucial part of the chapter. Hence, let us to make a brief insight into the two tasks.

Some authors, e.g. in (Brethé & Loiseau, 1998), pointed out that the use quasipolynomials does not permit to effectively handle some stabilization and control tasks, such as internal stability, controller properness etc. Hence, the quasipolynomial (meromorphic) description of TDS can be extended to a fractional one where a transfer function can be viewed as a field of fractions over a ring. The ring of stable and proper quasipolynomial (RQ) meromorphic functions (**R**_{MS}) (Zítek & Kučera, 2003; Pekař & Prokop, 2010; Pekař & Prokop, 2011a, 2011b; Pekař et al., 2011) is a suitable candidate. Although the ring can be used for a description of even neutral systems (Hale & Verduyn Lunel, 1993) after some definition adjustment, only systems with so-called retarded structure are considered as the admissible class of systems in this paper. Contrary to some other algebraic approaches, the ring enables to handle systems with non-commensurate delays, i.e. it is not necessary that all system delays can be expressed as integer multiples of the smallest one.Control philosophy in this ring then utilizes the Bézout identity, to obtain stable and proper controllers, along with the Youla-Kučera parameterization for reference tracking and load disturbance rejection.

Final controllers, in most cases, evince internal delays. Hence, a special effort to apply the control law in discrete time when controllers’ realization on computersmust be made. To name just a few TDS discretization methods, state space approaches can be found e.g. in (Engelborghs & Roose, 2002) based on the so-called solution operator, or in (Breda et al., 2005) via the so-called infinitesimal generator, orusing Taylor series expansion of the state vector in a neighborhood of the working point (Hofreiter, 2003). Input-output approaches include e.g. utilization of delta transform operators (Middleton & Goodwin, 1990; Zítek & Petrová, 2002). Nevertheless, all these methods consider a “very small” sampling period. An equivalent discrete-time model via the z-transform with a holder ought to be found when the sampling period is “higher”. To utilize the z-transform, one has to find a finite-dimensional continuous model. This task is the primary aim of this chapter, i.e. to approximate an infinitesimal model of an anisochronic controller such that the final finite-dimensional model can be subjected to the z-transform.

In recent decades a huge number of papers and works have been focused on model reduction or rational approximation of TDS, see e.g. (Makilla & Partington, 1999a, 1999b; Battle & Miralles, 2000). A fair overview of some methods and approaches has been published in (Partington, 2004). An overwhelming majority of these methods, however, deals with input-output delays only ignoring internal or state delays on the left-hand side of differential equations, i.e. those transfer functions with exponential terms in the denominator. In the contrary, this chapter focuses transfer function rational approximations for TDS with internal (state) delays since there are no theoretical analytic results about approximations convergence and accuracy for such models up today (distinct from single input-output delay which has been deeply analysed by rigorous analytic means many times) and one can thus expect interesting numerical comparative results. The basic question is whether the higher order approximation yields the better (i.e. more accurate) finite dimensional approximating model. The appropriate methods are chosen so that they are easy to deal with and anyone can use them effortlessly. Namely, Padé approximation, shift operator approximations – Laguerre and Kautz shift - and Fourier analysis based method are introduced and benchmarked. The common principle of all approaches consists in substitution of exponential terms in the transfer function of TDS with a rational fraction in the Laplace complex variable *s*.

We made no attempts to deal with analytic mathematical proofs in this chapter. From the practical point of view, the efficiency and accuracy of the selected approaches is measured and compared via norms in the Hardy space, namely, *H*_{∞} and *H*_{2} norms, the definitions of which and some calculation tricks are presented here as well (Štecha & Havlena, 2000).

A Matlab-Simulink user-interface application has been programmed to make a benchmark of approximations easier which enables to enter an approximated (nominal) plant transfer function, to select methods to be compared and the choice of norms by which the accuracy of approximation is calculated. It is also possible to specify an order of the approximation for each selected approximation approach. As the outputs, a table of calculated norms and the gain Bode plot (that is relevant to the norms) are displayed.Some study cases on approximations of anisochronic controllers’ derived by this algebraic method are introduced and compared in this contribution. For instance, control of an unstable system describing roller skater on a swaying bow (Zítek et al., 2008), or a laboratory heating system (Pekař et al., 2009), are taken as quasi-practical examples.

The chapter is organized as follows. In Section 2 a brief general input-output description of TDS is introduced together with the coprime factorization for the **R**_{MS} ring representation. Fundamentals of algebraic controller design in **R**_{MS} using a simple control feedback scheme are presented in Section 3. Section 4 introduces selected rational approximation methods for anisochronic controllers’ transfer functions. Definitions and calculation issues regarding to norms for approximations’ performance comparison are focused in Section 5. The application part of the chapter, i.e. Matlab-Simulink user interface, is described in Section 6. Finally, illustrative benchmark examples are presented in Section 7.

## 2. TDS models

Since the chapter applies finite-dimensional approximations of (controllers’) transfer functions, there is no point to introduce TDS systems description in the state space, regardless of the fact that state-space models appear as a result of system modelling and it would be natural to mention them first. For further details about state-space TDS models the reader is referred to (Richard, 2003).The second aim of this section is to present the concept of fields of fractions as an extension of meromorphic description in the form of quasipolynomial fractions, namely the ring of quasipolynomial meromorphic functions, **R**_{MS}, which is suitable to meet some elementary. So-called retarded TDS are primarily focused; nevertheless, neutral ones have to be mentioned as well, since it is necessary to make the ring definition complete.

### 2.1. Input-output quasipolynomial model

Since the authors’ interest lies in single-input single-output (SISO) TDS and their input-output models, they are concerned here. Namely, transfer functions in the form of quasipolynomial fractions giving rise to the meromorphic representation are taken as initial models to be approximated.

For both lumped (point-wise) and distributed input-output and internal delays in the system or model, the Laplace transform of a state space model (considering zero initial conditions) can be formulated as follows

where

where *n* is the order of a quasipolynomial which usually agrees with the system degree (of a state-space model), *neutral*; otherwise, the system is *retarded*.

### 2.2. TDS stability

A retarded system (1) is said to be *asymptotically* (exponentially) *stable* if all poles are located in the open left half plane, *s* satisfying

Condition *not* system poles.

A rather more complicated TDS stability conditions are given regarding to neutral systems, since there may exist vertical strips of system poles tending to the imaginary axis of infinite amplitude.Moreover, these strips can be sensitive to even infinitesimally small deviations in delays, i.e. the position of system poles in the real axis is not continuous with respect to delays. A neutral system (1) is exponentially stable if there is no *s* such that

for arbitrarily

A system is *strongly stable* if so-called associated difference equation (in state-space formulation) remains asymptotically stable when subjected to small variations in delays (Hale &Verduyn Lunel, 1993), i.e. a TDS remains *formally stable*, see details in (Byrnes et al., 1984; Loiseau et al., 2002). Formal stability, roughly speaking, means that the rightmost vertical strip of poles does not cross the imaginary axis. If this holds also under small delay changes, the system is strongly stable. A necessary and sufficient strong stability (and thus also formal stability) condition in the Laplace transform can be formulated as

according to e.g. (Zítek & Vyhlídal, 2008), where *s*-power in the denominator

### 2.3. TDS model over RMS ring

Algebraic approaches for analysis and control of TDS can be performed either in the state space or in the realm of input-output models where fields, modules and rings as principal algebraic notions and tools are utilized. Usually, commensurate delays, i.e. those which can be expressed as integer multiples of the smallest one, are assumed; however, delays are naturally real-valued and thus this assumption is rather restrictive for real applications.Non-commensurate or rationally unapproximated delays results in a fraction of quasipolynomials as presented above. However, these transfer function representations are not suitable in order to satisfy some basic control requirements, e.g. controller feasibility, closed-loop (Hurwitz) and formal (strong) stability.

Rather more general approaches (Vidyasagar, 1985; Kučera, 1993) utilize a field of fractions where a transfer function is expressed as a ratio of two coprime elements of a suitable ring. A ring is a set closed for addition and multiplication, with a unit element for addition and multiplication and an inverse element for addition, i.e. division is not generally allowed. A powerful algebraic tool ensuring some basic control requirements, such as internal closed-loop stability and controller properness, is a ring of stable and proper RQ-meromorphic functions (**R**_{MS}). Since the original definition of **R**_{MS} in (Zítek & Kučera, 2003) does not constitute a ring, some minor changes in the definition were made e.g. in (Pekař & Prokop, 2009). Namely, although the retarded structure of TDS is considered only, the minimal ring conditions require the use of neutral quasipolynomials at least in the numerator as well. In this chapter, the ring definition is reformulated once more to comprise models of neutral type, distributed delays and formal stability.

A term **R**_{MS}ring is represented by a proper ratio of two quasipolynomials *n* and a numerator can be factorized as

where *l* and*s*-power. The element is analytic and bounded in *s*_{0} such that

It is said that

Notice, for instance, that *H*_{∞} unstable due to an unbounded gain at the imaginary axis.

Moreover,

In addition, the ratio is proper, i.e. *l**n*. More precisely, there exists a real number *R*> 0 for which holds that

### 2.4. Coprime factorization in RMS

Let the plant be initially described by the transfer function

where

where

## 3. Controller design in R_{MS}

The aim of this section is to outline controller design based on the algebraic approach in the **R**_{MS} ring satisfying the closed loop stability in that sense that all transfer functions in the feedback are from the ring (i.e. they lies in *H*_{∞} and are proper and formally stable) controller feasibility, reference tracking and load disturbance rejection. As a control system, the simple feedback loop is chosen for the simplicity, see Fig. 1.

For algebraic controller design in **R**_{MS} it is initially supposed that not only the plant is expressed by the transfer function over **R**_{MS} but a controller and all system signals are over the ring. Let

External inputs, reference and load disturbance signals, respectively, have forms

where **R**_{MS}.

The following important feedback transfer functions can be derived

where the controller transfer function is factorized as follows

and the common denominator meromorphic function reads

Meromorphic functions **R**_{MS} and the fraction (16) is (Bézout) coprime. The numerator of **R**_{MS} agrees to the characteristic quasipolynomial of the closed loop generally defined in (2).

A brief description of basic control design steps follows.

### 3.1. Closed loop stabilization

Given a Bézout coprime pair **R**_{MS} the closed-loop system is stable if and only if there exists a pair

A particular stabilizing solution of (17), say

where

The proof of the statement above can be done analogously as in (Zítek & Kučera, 2003) where a three-step proof for a similar ring was presented. Condition (12) ensures i.a. that there can exist the ring inversion of

### 3.2. Reference tracking

The task of this subsection is to find

where ^{+}) – this requirement also satisfies that

Alternatively,

Note that the controller approach fails for formally unstable controlled processes since then the feedback loop remains formally unstable (neutral term can not be affected by a controller).

### 3.3. Load disturbance rejection

The attenuation of the load disturbance signal entering a plant model can be done analogously as for reference tracking. Thus, *Z(s)* is chosen so that

where

One has to be careful when deciding about the form of

## 4. Selected transfer function rational approximations

Selected easy-handling transfer function rational approximation methodsbased on the substitution of exponential elements by a rational function are described in this section. The emphasis is put on the engineering usability of the methods.

### 4.1. Padé approximation

In the second half of the 19^{th} century, a French mathematician Henry Padé devised a simple and, nowadays, one of the most used and favorite rational approximations which is based on the comparison of derivatives of the approximating and approximated functions in zero. More precisely, let *n*-*n* Padé approximation is the function *n*th order with

Padé approximation of

where *n* is the order of the approximation. Obviously, one can approximate another function, e.g. the whole transfer function.

A method called *diagonal* Padé approximation,which is distinguished by some authors, see (Battle & Miralles, 2000;Richard, 2003) can be expressed as

However, it is easy to verify, that (21) and (22) represent the same approximations. In fact,

### 4.2. Shift operator approximations

Approaches based on operator shifting yields from the fact that a delay term *s* can be vied as a derivative operator.

Hence

A concise overview of some important shift operator approaches follows.

#### 4.2.1. Laguerre shift

The eventual relation for this shift operator approximant is given by the formula

where *n* is the order of the approximation again, see (Makilla & Partington, 1999a) for more details. The Laguerre shift is successfully used in robust control.

#### 4.2.2. Kautz shift

As presented in (Makilla & Partington, 1999a), the *n*th order Kautz shift reads

In the source referenced above, there has been analytically proved for input-output delays that the Kautz shift approximation is asymptotically twice more accurate than the Laguerre one.

#### 4.2.3. Padé shift

This type of shift approximation is based on the second order Padé approximation (Makilla & Partington, 1999b; Battle & Miralles, 2000) and it is also called Padé-2. It can be formulated e.g. as follows

#### 4.2.4. Fourier analysis based method

This approximation methodology results from the analysis of the delay-in-feedback step response or, generally, the response to the input of the form

To demonstrate the initial idea more precisely, consider a single delay term

Hence

and the output *p(t)*

Inserting (30) into (29) yields

From the Fourier series expansion of

The exponential term can be then approximated by the inverse of (28) as

The reader is referred to (Battle & Miralles, 2000) for more details about the derivation and methodology.

## 5. Approximation performance evaluation using norms

The accuracy of rational approximation methods and approaches has been usually measured and evaluated via conscientious and rigorous mathematic analytic proofs. The derivation of such results has been facilitated by the fact that delays in input-output relation only, ignoring the effect of internal delays, have been considered. This is the main gap which ought to be healed up; however, no analytic methods are used in this chapter anyway.

Transfer function norms in the Hardy space known from robust control for unstructured uncertainty measurement instead of any complex analytic method for the comparison of approximation accuracy are used in this contribution since this conception is sufficient for particular cases and for engineering practice. Moreover, analytic derivations of accuracies for internal delays would be much more difficult than those for input-output ones.

First of all, define the objective to be minimized in various meanings. The difference between the nominal

We chose two the most used norms in the Hardy space, namely, *H*_{2} and *H*_{} norms, which are effective to express the proximity of transfer functions in the frequency domain. The definitions and some basic properties are then followed by some calculation tricks when dealing with the norms.

### 5.1. H_{2} norm

The *H*_{2},sometimes called quadratic, norm of a stable strictly proper transfer function is defined as

The norm is finite for strictly proper stable systems having no pole on the imaginary axis, and the meaning of *H*_{2} is energy of

Note that for TDS with distributed delays, there can exist a denominator root (or roots) of

In computer practice, i.e. working with discrete samples, the integral in (36) is calculated as a sum within a finite range of nonnegative frequencies,

Residual expansion can be used when analytic (and continuous) calculation of

where

where *m* is the multiplicity of a pole.

### 5.2. H_{∞} norm

The H-infinity (

i.e. it expresses the supreme of the amplitude (gain) frequency characteristics of *H*_{∞}norm, it is said that it is *H*_{∞} stable and lies in the space

The norm is also called *L*_{2} gain, That is, the *H*_{∞} stable system has finite

The frequency characteristics supreme can be easily found by standard analytic means, or by mapping the values of

Note, for instance, that a transfer function having no pole on the imaginary axis but a sequence of poles with real part converging to zero can have an infinite *H*_{∞} norm due to an unbounded gain, see (Partington & Bonnet, 2004).

## 6. Matlab-Simulink user interface

The main application part of this chapter consists in development of a simulation program that can evaluate the quality of approximations. A user-friendly interface in Matlab-Simulink environment has been hence developed by the authors in order to provide testing and comparison of approximation approaches introduced above. The source code can be found in (Pekař & Kurečková, 2011).

Fig. 3 displays the working environment of the application. In the upper part of the user interface, one is allowed to enter a nominal plant transfer function as a numerator and denominator quasipolynomials. The user must be careful since the syntax scheme which is indicated by grey font (or by the help hint) in the edit box must be kept.

Then, a user can select one or more approximation methods to be compared and check the choice of norms by which the accuracy of approximation is calculated. It is also possible to specify an order of the approximation for each selected method.

By clicking the button *Get approximation accuracy*, the programme returns a table of calculated normsin the text form.A graphical comparison of amplitude frequency responses for all chosen approximations is obtained by clicking the button *Plot amplitude characteristics*.

Final approximation controller transfer functions are returned to the Matlab workspace to the command line.

## 7. Benchmark examples

The aim of this section is to verify introduced approximation methods using the user interface. The approximation accuracy with respect to the order of an approximation is a very attractive question as well. As will be seen from the following examples, some rather surprising resultscan be obtained.

As mentioned above, results of the benchmark ought to be primarily used by the authors for the rational approximation of controllers with internal delays (also called anisochronic) as a first step of controllers’ discretization via the z-transform, for the computer implementation.

Three examples of algebraic controller design with rational approximations using the user interface follow. The first one gives results for anisochronic controllers design in control of a hypothetic simple stable LTI with stepwise reference, to demonstrate mainly the design procedure. The second one presents the approximation of an anisochronic controller for an attractive unstable system of a skater on a swaying bow. The third example deals with control of a laboratory circuit heating system, as a typical real-life representative of retarded TDS.

### 7.1. Stable system with stepwise reference

Let the plant be described by the transfer function of a stable first order TDS model with both internal and input-output delays as

Condition

Coprime (Bézout) factorization is the first step of controller design in **R**_{MS}ring as follows

where

Consider the simplest practical case that both external inputs are from the class of step functions, hence

where

Find a stabilizing particular solution by (17). Set e.g.

Now parameterize the solution according to (18) to obtain controllers asymptotically rejecting the load disturbance and tracking the stepwise reference

The numerator of *P*(*s*) has to have at least one zero root. Moreover, it is appropriate to have *P*(*s*) in a simple form, which is fulfilled e.g. when

providing

Thus, final controller’s structure is the following

The obtained control structure can be easily compared with the well-known Smith predictor structure and note that the controller is of the anisochronic type because of delay in the transfer function denominator. It is naturally possible to take

i.e.

Model (41) can to fit the dynamics of a high order undelayed system; for instance, a tenth order system governed by the transfer function

can be estimated by model (41) with

*n*is the order of the approximation method, and test the accuracies for all methods introduced above. The best results for each of them measured by

*H*

_{∞}and

*H*

_{2}norms with the corresponding order are displayed in Table 1.

Obviously, the best result for *H*_{∞} is given alongside by the Padé approximation and Laguerre shift of the first order, whereas, amazingly, higher orders make results worse. The Fourier analysis based methods yields almost the same score for all studied orders.The benchmark results for the *H*_{2} norm with *H*_{∞}, i.e. the Padé approximation and Laguerre shift of the first order are the best and the Fourier analysis based methods gives almost the same results for all orders.The corresponding gain frequency responses for the approximations of orders as in the last column in Table 1 are displayed in Fig. 4.

Method | H_{(} | n | H_{2} | n |

Padé approx. | 0.243 | 1 | 0.237 | 1 |

Laguerre shift | 0.243 | 1 | 0.237 | 1 |

Kautz shift | 0.408 | 5 | 0.299 | 3 |

Padé shift | 0.277 | 1 | 0.29 | 1 |

Fourier analysis | 0.38 | 1 | 0.354 | 5 |

The approximating transfer function by Padé approximation and Laguerre shift with *n* = 1 is given by (51), which agrees with a conventional PID controller

### 7.2. Control of the roller skater on a swaying bow

Consider an unstable system describing roller skater on a swaying bow, (Zítek et al., 2008), governed by the transfer function

see Fig. 5, where *P*, delays *b*, *a* are real parameters. Skater controls the servo driving by remote signals into servo electronics.

Let the model parameters be *b* = 0.2, *a* = 1,

where *p*_{2}, *p*_{1}, *p*_{0}, *q*_{3}, *q*_{2}, *q*_{1}, *q*_{0}

and e.g.

The comparison of the best controller rational approximations can be found in Table 2. Again, the method based on the Fourier series expansion very slowly approaches the limit value of the *H*_{2}norm (*n*. The only method evincing the better asymptotical results with the higher order approximation is the Kautz shift. Again, the Padé and Laguerre approximations of the first order give very good results with the approximating controller transfer function

Method | H_{(} | n | H_{2} | n |

Padé approx. | 5.6674 | 1 | 8.577 | 1 |

Laguerre shift | 5.6674 | 1 | 8.577 | 1 |

Kautz shift | 6.0847 | 5 | 9.7446 | 5 |

Padé shift | 6.0668 | 1 | 9.674 | 1 |

Fourier analysis | 5.6674 | 1 | 7.8114 | 5 |

Fig. 6 displays Bode magnitude plots for the best orders for *H*_{2} (the last column in Table 2), which verifies a very good performance of all the approaches.

### 7.3. Control of a circuit heating plant

The laboratory heating plant, a photo and a sketch of which, respectively, are displayed in Fig. 7, was assembled at the Faculty of Applied Informatics of Tomas Bata University in Zlín in order to test control algorithms for systems with dead time. The original description of the apparatus and its electronic circuits can be found in (Dostálek et al., 2008).

The heat transferring fluid (namely distilled water) is transported using a continuously controllable DC pump {6} into a flow heater {1} with maximum power

Originally, it was intended to control input delays only; however, it was shown that the plant contains internal delays as well. A detailed mathematical model was presented in (Pekař et al., 2009). A linearized model of the relation between the power to the heater

where all real parameters in the model are complex algebraic functions of physical quantities in the circuit and input and output steady states, see details in the literature. It was determined that for a certain working point, the parameters are

The controller structure obtained by controller design in

Table 3 displays the best values for approximations at their orders*H*_{2} norm is measured within the frequency range

Method | H_{(} | n | H_{2} | n |

Padé approx. | 27994.7 | 4 | 1197.3 | 5 |

Laguerre shift | 28528.6 | 3 | 1194.3 | 2 |

Kautz shift | 28365.9 | 2 | 1196.4 | 2 |

Padé shift | 28650.2 | 3 | 1198.7 | 1 |

Fourier analysis | 28760.9 | 1 | 1205.1 | 1 |

Immense values of *H*_{∞} norm are caused by a pair controller poles which is very close to the imaginary axis (^{th} order provides the best value for *H*_{∞}, whereas Laguerre shift of the 2^{nd} order gives the best approximation measured by *H*_{2} norm. Again, this example indicates that the higher order of the approximation does not mean the more accurate result automatically.

The comparison of Bode magnitude plots for the best orders (which are placed in the last column in Table 3) for *H*_{2} can be seen in Fig. 8.

To conclude study cases above, it is startling that the best approximations measured by *H*_{2} and *H*_{∞} norms are mostly given by the well known and widely used Padé approximation of the first (or a low) order. By simulations, the higher order of an approximation does not generally yields the more accuracy finite dimensional model, which is in the contradiction with a general expectation and analytic results for rational approximations for TDS with input-outputs delays. Comparative Bode plots above indicate the usability and a very good efficiency of selected methods.

## 8. Conclusion

This chapter intended to propound the reader a methodology for algebraic controller design for systems with internal delays, followed by a comparison of several easy-handling techniques for rational (i.e. finite-dimensional) approximation of anisochronic (i.e. infinite-dimensional) controllers – or their transfer functions, more precisely.

The first, controller design, part was based on the solution of the Bézout identity, to obtain stable and proper controllers, along with the Youla-Kučera parameterization for reference tracking and load disturbance rejection, in **R**_{MS} ring.

The second, controller approximation, part of this chapter intended to use five quite simple methods, namely Padé approximation, Laguerre, Kautz shift and Fourier analysis based method, how to approximate a delay exponential term in the transfer function. Two norms in the Hardy space, namely, *H*_{2} and *H*_{∞} norms were used as a measure for the accuracy and efficiency of the approaches since there are no rigorous analytic results about the approximation efficiency and accuracy in the literature and these norms are usually sufficient in technical practice. Approximation methods were selected so that they are serviceable also for people with basic mathematical background knowledge. The authors plan to use approximated controllers in their digital implementation on PC or PLC.

Matlab with Simulink was a very useful assistant here. The authors programmed a simple user interface which enables the user to enter a nominal transfer function and select approximation methods to be used and their orders. As a result, the program returns the accuracies in both text and graphical forms.

Simulation experiments with the program were made. Control of a simple stable TDS, control of unstable TDS of a skater on the swaying bow and control of a laboratory circuits heating plant were benchmark examples. The results were very interesting and startling because the habitual Padé approximation proved to be very good and, moreover, the higher order approximation did not automatically mean the better result for systems with internal delays.

## Acknowledgement

The authors kindly appreciate the financial support which was provided by the Ministry of Education, Youth and Sports of the Czech Republic, under the project CEBIA-Tech No. CZ.1.05/2.1.00/03.0089.