Nonlinear Dynamics 29: 281–296, 2002.

© 2002 Kluwer Academic Publishers. Printed in the Netherlands.

Analogue Realizations of Fractional-Order Controllers

I. PODLUBNY, and I. PETRÁŠ,

Department of Informatics and Process Control, BERG Faculty, Technical University of Košice,
04200 Košice, Slovak Republic

B. M. VINAGRE
Department of Electronics and Electromechanical Engineering, Industrial Engineering School,
University of Extramadura, E-06071 Badajoz, Spain

P. O’LEARY
Institute of Automation, Montanuniversität Leoben, A-8700 Leoben, Austria

L’. DORČÁK
Department of Informatics and Process Control, BERG Faculty, Technical University of Košice,
04200 Košice, Slovak Republic

(Received: 16 July 2001; accepted: 7 December 2001)

Abstract. An approach to the design of analogue circuits, implementing fractional-order controllers, is presented.
The suggested approach is based on the use of continued fraction expansions; in the case of negative coefficients
in a continued fraction expansion, the use of negative impedance converters is proposed. Several possible methods
for obtaining suitable rational appromixations and continued fraction expansions are discussed. An example of
realization of a fractional-order Iλ controller is presented and illustrated by obtained measurements. The suggested
approach can be used for the control of very fast processes, where the use of digital controllers is difficult or
impossible.

Keywords: Fractional calculus, fractional differentiation, fractional integration, fractional-order controller, real-
ization.

1. Introduction

Although digital controllers are used more and more frequently for controlling many types of
complex processes, the role of analogue controllers should not be undervalued. Indeed, digital
controllers have some natural limitations, coming from their discrete nature, such as the length
of the sampling period and the time of computation, which should be significantly less than
the length of the sampling period. This sometimes makes the use of digital controllers prac-
tically impossible, especially in case of fast processes, such as vibrations, and the alternative
approach to controlling fast processes is represented by analogue controllers.
In this paper we describe an approach to the design of analogue fractional-order controllers.
The paper is organized as follows. First, we recall some basic relationships for describing
fractional-order systems and fractional-order controllers. Then we discuss some uses of con-
tinued fraction expansions, including their applications in the control theory. Finally, we show
how continued fraction expansions can be used for designing analogue circuits, implementing
Partially supported by grant No. VEGA 1/7098/20 of the Slovak Grant Agency for Science.
Partially supported by FEDER Research Grant No. IFD97-0755-C02-01.
Partially supported by a research grant of the Austrian Institute for Central and Eastern Europe.
282 I. Podlubny et al.

fractional-order systems and controllers. We also give an example of implementation of an Iλ

controller.

2. Fractional-Order Systems and Controllers

General information about various approaches to fractional-order differentiation and integra-

tion can be found in the available monographs on this subject [1–4] and in some other articles
in this special issue. Because of this, we do not discuss general definitions here. Instead,
we recall only the expressions for describing fractional-order systems and PIλ Dµ controllers
[3, 5], which are subjects of our interest in this paper.

2.1. F RACTIONAL D IFFERENTIAL E QUATIONS AND T RANSFER F UNCTIONS

A fractional-order control system can be described by a fractional differential equation of the

form
an D αn y(t) + an−1 D αn−1 y(t) + · · · + a0 D α0 y(t)
= bm D βm u(t) + bm−1 D βm−1 u(t) + · · · + b0 D β0 u(t), (1)
or by a continuous transfer function of the form:
bm s βm + bm−1 s βm−1 + · · · + b0 s β0
G(s) = , (2)
an s αn + an−1 s αn−1 + · · · + a0 s α0
γ
where D γ ≡ 0 Dt denotes the Riemann-Liouville or Caputo fractional derivative [3]; ak
(k = 0, . . . , n), bk (k = 0, . . . , m) are constant; and αk (k = 0, . . . , n), βk (k = 0, . . . , m)
are arbitrary real numbers.
Without loss of generality we can assume that αn > αn−1 > . . . > α0 , and βm > βm−1 >
. . . > β0 .

2.2. PIλ Dµ C ONTROLLERS

The fractional-order PIλ Dµ controller was proposed in [3, 5, 6] as a generalization of the

PID controller with integrator of real order λ and differentiator of real order µ. The transfer
function of such type the controller in Laplace domain has form:
U (s)
Gc (s) = = K + Ti s −λ + Td s µ (λ, µ > 0), (3)
E(s)
where K is the proportional constant, Ti is the integration constant and Td is the differentiation
constant. As we can see (Figure 1), the internal structure of the fractional-order controller con-
sists of the parallel connection the proportional, integration, and derivative part [7]. Transfer
function (3) corresponds in time domain with fractional differential equation (4)

u(t) = Ke(t) + Ti 0 Dt−λ e(t) + Td 0 Dtµ. (4)

Taking λ = 1 and µ = 1, we obtain a classical PID controller. If λ = 0 and/or Ti = 0, we

obtain a PDµ controller, etc. All these types of controllers are particular cases of the fractional-
order controller, which is more flexible and gives an opportunity to better adjust the dynamical
properties of the fractional-order control system.
Analogue Realizations of Fractional-Order Controllers 283

Figure 1. General structure of a PIλ Dµ controller.

As we see from Figure 1, a PIλ Dµ controller can be easily implemented in analogue

form if we know how to build an analogue circuit corresponding to s α , α ∈ R. Below we
will demonstrate how this can be done using rational approximations and continued fraction
expansions.
It can also be mentioned that the other kind of fractional-order controller, which character-
ized by the band-limited lead effect, can be found in the available literature [8, 9]:

1 + τ sr
Gc (s) = C , r ∈ R, C ∈ R, τ < τ. (5)
1 + τ sr
This type of controller can be realized using a recursive distribution of poles and zeros [10].

3. Some Uses of Continued Fractions

In this section we discuss some applications of continued fractions. First we recall their use
for approximating functions and investigating stability of linear systems. Then we introduce a
new relationship between continued fractions and multiple nested-loop systems.

3.1. CFE S AND A PPROXIMATIONS OF F UNCTIONS

It is well known that the Continued Fraction Expansions (CFE) is a method for evaluation
of functions, that frequently converges much more rapidly than power series expansions, and
converges in a much larger domain in the complex plane [11]. The result of such approxima-
tion for an irrational function, G(s), can be expressed in the form:
b1 (s)
G(s) a0 (s) +
a1 (s) + b2 (s)
b3 (s)
a2 (s)+ a (s)+···
3

b1 (s) b2 (s) b3 (s)

= a0 (s) + ··· , (6)
a1 (s)+ a2 (s)+ a3 (s)+
where ai s and bi s are rational functions of the variable s, or are constant. The application
of the method yields a rational function, G (s), which is an approximation of the irrational
function G(s).
On the other hand, for interpolation purposes, rational functions are sometimes superior to
polynomials. This is, roughly speaking, due to their ability to model functions with poles. (As
it can be seen later, branch points can be considered as accumulations of interlaced poles and
zeros). These techniques are based on the approximations of an irrational function, G(s), by
284 I. Podlubny et al.

a rational function defined by the quotient of two polynomials in the variable s:

Pµ (s)
G(s) Ri(i+1)...(i+m) = ,
Qν (s)
p0 + p1 s + · · · + pµ s µ
= (7)
q0 + q1 s + · · · + qν s ν
m + 1 = µ + ν + 1,
passing through the points (si , G(si )), . . . , (si+m , G(si+m )).

3.2. CFE AND S TABILITY OF L INEAR S YSTEMS

It is also known that continuous fraction expansions can be used for investigating stability of
linear systems. For this, the characteristic polynomial Q(s) of the differential equation of the
system should be divided in two parts, the ‘even’ part (containing even powers of s) and the
‘odd’ part (containing odd powers of s):

Q(s) = m(s) + n(s).

Then these two parts of the characteristic polynomial are used for creating its test function in
the form of a fraction, in which the highest power of s is contained in the denominator:

m(s) n(s)
R(s) = or R(s) = .
n(s) m(s)
The rational function R(s) should be written in the form of a continuous fraction:
1
R(s) = . (8)
1
b1 s +
1
b2 s +
............
1
bn−1 s +
bn s
If bk > 0, k = 1, . . . , n, then the system is stable. If some bk is negative, then the system
is unstable.
Considering the continued fraction (8) as a tool for designing a corresponding LC circuit,
we can conclude that stability of a linear system is equivalent to realizability of its test function
R(s) with the help of only passive electric components.

3.3. CFE AND N ESTED M ULTIPLE -L OOP C ONTROL S YSTEMS

Let us now establish an interesting new relationship between continued fractions and nested
multiple-loop control systems.
We first recall the known fact that the transfer function R(s) of the control loop with a
negative feedback shown in Figure 2 is given by [7]
G(s)
R(s) = . (9)
1 + G(s)H (s)
Analogue Realizations of Fractional-Order Controllers 285

Figure 2. A control loop with a negative feedback.

Figure 3. Nested multiple-loop control system – level 1.

From (9) it immediately follows that the transfer function of the circuit shown in Figure 3
is
1 1
P2n (s) = ∗ = , (10)
1 + 1 · Y2n (s) Y2n(s)
∗
where Y2n (s) = Y2n (s) + 1.
Using Equations (9) and (10) we obtain the transfer function of the system shown in
Figure 4:
1
Q2n−1 (s) = Z2n−1 (s) + P2n (s) = Z2n−1 (s) + . (11)
Y2n(s)
Combining Equations (9) and (10) we find the transfer function of the nested multiple-loop
system shown in Figure 5:
Q2n−1 (s) 1
P2n−2 (s) = =
1 + Q2n−1 (s)Y2n−2 (s) 1
Y2n−2 (s) +
Q2n−1 (s)
1
= (12)
1
Y2n−2 (s) +
1
Z2n−1 (s) + .
Y2n (s)
The transfer function of the system shown in Figure 6 is then given by the relationship

Q2n−3 (s) = Z2n−3 (s) + P2n−2 (s)

1
= Z2n−3 (s) + (13)
1
Y2n−2 (s) +
1
Z2n−1 (s) + .
Y2n(s)
286 I. Podlubny et al.

Figure 4. Nested multiple-loop control system – level 2.

Figure 5. Nested multiple-loop control system – level 3.

Continuing this process, we obtain the transfer function of the nested multiple-loop control
system shown in Figure 7 in the form of a continued fraction expansion, which is identical
with the Equation (24):
1
Z(s) = Z1 (s) +
1
Y2 (s) +
1
Z3 (s) +
1
Y4 (s) +
.................................
1
1
Y2n−2 (s) +
1
Z2n−1 (s) + .
Y2n (s)
Similarly to the above considerations, we can obtain a continued fraction expansion of
the transfer function of the other interesting type of a nested multiple-loop control system,
depicted in Figure 8:
1
Z(s) = (14)
1
Z1 (s) +
1
Y2 (s) +
1
Z3 (s) +
.................................
1
1
Y2n−2 (s) +
1
Z2n−1 (s) + .
Y2n (s)
Analogue Realizations of Fractional-Order Controllers 287

Figure 6. Nested multiple-loop control system – level 4.

Figure 7. Nested multiple-loop control system of the first type.

Figure 8. Nested multiple-loop control system of the second type.

288 I. Podlubny et al.

Both types of nested multiple-loop systems, presented in this section, can be used for
simulations and realizations of arbitrary transcendental transfer functions. For this, the transfer
function should be developed in a continued fraction, which after truncation can be represen-
ted by a nested multiple-loop system shown in Figures 7 or 8.

4. CFE and Rational Approximations of s α

In general [12], a rational approximation of the function G(s) = s −α , 0 < α < 1 (the frac-
tional integral operator in the Laplace domain) can be obtained by performing the CFE of the
functions:
1
Gh (s) = , (15)
(1 + sT )α

1 α
Gl (s) = 1 + , (16)
s
where Gh (s) is the approximation for high frequencies (ωT
1), and Gl (s) the approxima-
tion for low frequencies (ω 1).

EXAMPLE 1. Performing the CFE of the function (15), with T = 1, α = 0.5, we obtain
0.3513s 4 + 1.405s 3 + 0.8433s 2 + 0.1574s + 0.008995
H1 (s) = .
s 4 + 1.333s 3 + 0.478s 2 + 0.064s + 0.002844
EXAMPLE 2. Performing the CFE of the function (16), with T = 1, α = 0.5, we obtain
s 4 + 4s 3 + 2.4s 2 + 0.448s + 0.0256
H2 (s) = .
9s 4 + 12s 3 + 4.32s 2 + 0.576s + 0.0256

5. Other Rational Approximations for s α

Besides using continued fractions, there are also other methods [13] for obtaining rational
approximations of fractional-order systems. However, since a ratio of two polynomials can be
expressed in the form of a finite continued fraction, any rational approximation is equivalent
to a certain finite continued fraction.

5.1. C ARLSON ’ S M ETHOD

The method proposed by Carlson in [14], derived from a regular Newton process used for
iterative approximation of the α-th root, can be considered as belonging to this group. The
starting point of the method is the statement of the following relationships:
(H (s))1/α − (G(s)) = 0; H (s) = (G(s))α . (17)
Defining α = 1/q, m = q/2, in each iteration, starting from the initial value H0 (s) = 1,
an approximated rational function is obtained in the form:
(q − m) (Hi−1 (s))2 + (q + m)G(s)
Hi (s) = Hi−1 (s) . (18)
(q + m) (Hi−1 (s))2 + (q − m)G(s)
Analogue Realizations of Fractional-Order Controllers 289

EXAMPLE 3. Starting from H (s) = (1/s)1/2 , H0 (s) = 1, after two iterations, we obtain

s 4 + 36s 3 + 126s 2 + 84s + 9

H3 (s) = .
9s 4 + 84s 3 + 126s 2 + 36s + 1

5.2. M ATSUDA’ S M ETHOD

The method suggested in [15] is based on the approximation of an irrational function by a

rational one, obtained by CFE and fitting the original function in a set of logarithmically
spaced points. Assuming that the selected points are sk , k = 0, 1, 2, . . ., the approximation
takes on the form:
s − s0 s − s1 s − s2
H (s) = a0 + ··· , (19)
a1 + a2 + a3 +
where
s − si
ai = vi (si ), v0 (s) = H (s), vi+1 (s) = . (20)
vi (s) − ai

EXAMPLE 4. With G(s) = (1/s)1/2 , finitial = 1, ffinal = 100, fk = {1, 1.7783, 3.1623,
5.6234, 10, 17.783, 31.623, 56.234, 100}, we obtain

0.08549s 4 + 4.877s 3 + 20.84s 2 + 12.995s + 1

H4 (s) = .
s 4 + 13s 3 + 20.84s 2 + 4.876s + 0.08551

5.3. O USTALOUP ’ S M ETHOD

The method [8–10] is based on the approximation of a function of the form:

H (s) = s δ , δ ∈ R+ (21)

by a rational function

N
1 + s/ωk
(s) = C
H , (22)

k=−N
1 + s/ω k

using the following set of synthesis formulas:

ωk+1 ωk+1
ω0 = α −0.5 ωu , ω0 = α 0.5 ωu , = = αη > 1,
ωk ωk

ωk+1 ωk log (ωN /ω0 ) log α
= η > 0, = α > 0, N= , δ= , (23)
ωk ωk log (αη) log (αη)
with ωu being the unit gain frequency and the central frequency of a band of frequencies
√
geometrically distributed around it. That is, ωu = ωh ωb , ωh , ωb are the high and low
transitional frequencies.
290 I. Podlubny et al.

Figure 9. Finite ladder circuit.

EXAMPLE 5. Using the Oustaloup’s method with

ωh = 102 , ωb = 10−2 ,
from which we have α = η = 2.5119, the obtained approximation for s −1/2 is
s 5 + 74.97s 4 + 768.5s 3 + 1218s 2 + 298.5s + 10
H5 (s) = .
10s 5 + 298.5s 4 + 1218s 3 + 768.5s 2 + 74.97s + 1

6. Design of Fractances Based on Rational Approximations and CFEs

A circuit exhibiting fractional-order behaviour is called a fractance [3].

Design of fractances can be done easily using any of the aforementioned rational approxim-
ations or a truncated CFE, which also gives a rational approximation (see, for example, [16]).
Truncated CFE does not require any further transformation; a rational approximation based
on any other methods must be transformed to the form of a continued fraction. The values
of the electric elements, which are necessary for building a fractance, are then determined
from the obtained finite continued fraction. If all coefficients of the obtained finite continued
fraction are positive, then the fractance can be made of classical passive elements (resistors
and capacitors). If some of the coefficients are negative, then the fractance can be made with
the help of negative impedance converters (Section 6.2).

6.1. D OMINO L ADDER C IRCUIT

Let us consider the circuit depicted in Figure 9, where Z2k−1 (s) and Y2k (s), k = 1, . . . , n, are
given impedances of the circuit elements. The resulting impedance Z(s) of the entire circuit
can be found easily, if we consider it in the right-to-left direction:
1
Z(s) = Z1 (s) + (24)
1
Y2 (s) +
1
Z3 (s) +
1
Y4 (s) +
.................................
1
1
Y2n−2 (s) +
1
Z2n−1 (s) + .
Y2n (s)
The relationship between the finite domino ladder network, shown in Figure 9, and the con-
tinued fraction (24) provides an easy method for designing a circuit with a given impedance
Analogue Realizations of Fractional-Order Controllers 291

Z(s). For this one has to obtain a continued fraction expansion for Z(s). Then the obtained
particular expressions for Z2k−1 (s) and Y2k (s), k = 1, . . . , n, will give the types of necessary
components of the circuit and their nominal values.

EXAMPLE 6. To design a circuit with the impedance

s 4 + 4s 2 + 1
Z(s) = , (25)
s3 + s
we have to develop Z(s) in continued fraction

s 4 + 4s 2 + 1 1
Z(s) = =s+ (26)
s +s
3 1
s+
1
3 9 1
s+
2 2
s.
3
From this expansion it follows that
9 1 2
Z1 (s) = s, Z3 (s) = s, Y2 (s) = s, Y4 (s) = s.
2 3 3
Therefore, for the analogue realization in the form of the first Cauer’s canonic LC circuit
[17] we have to choose the following values of coils and capacitors:
9 1 2
L1 = 1 [H ], L3 = [H ], C2 = [F ], C4 = [F ].
2 3 3
EXAMPLE 7. The function Z(s) given by Equation (25) can be written also in the form

s 4 + 4s 2 + 1 1 1
Z(s) = = + (27)
s +s
3 s 1
+
1
3s 9 1
+
2s 2
.
3s
From this expansion it follows that
1 9 1 2
Z1 (s) = , Z3 (s) = , Y2 (s) = , Y4 (s) = .
s 2s 3s 3s
Therefore, for the analogue realization in the form of the second Cauer’s canonic LC circuit
[17] we have to choose the following values of coils and capacitors:
2 3
C1 = 1 [F ], C3 = [F ], L2 = 3 [H ], L4 = [H ].
9 2
EXAMPLE 8. To design a circuit with the impedance

s 4 + 3s 2 + 8
Z(s) = , (28)
2s 3 + 4s
292 I. Podlubny et al.

one has to obtain a continuous fraction representation of the function Z(s),

s 4 + 3s 2 + 8 1 1
Z(s) = = s+ (29)
2s + 4s
3 2
2s +
1
1 1
− s+ .
12 3
− s
2
From this expansion it follows that
1 1 3
Z1 (s) = s, Z3 (s) = − s, Y2 (s) = 2s, Y4 (s) = − s.
2 12 2
Therefore, for the analogue realization in the form of the first Cauer’s canonic LC circuit
[17] we have to choose the following values of coils and capacitors:
1 1 3
L1 = [H ], L3 = − [H ], C2 = 2 [F ], C4 = − [F ].
2 12 2
Here we see negative inductances and capacitance. Such elements cannot be realized using
passive electric components. However, they can be realized with the help of active compon-
ents, namely operating amplifiers.

6.2. N EGATIVE -I MPEDANCE C ONVERTERS

The previous example shows that the use of CFE for analogue realization of arbitrary transfer
functions may lead to the appearance of negative impedances. This observation is not un-
known. For example, in the paper [12], Dutta Roy recalls Khovanskii’s continued fraction
expansion for x 1/2 found in [18] and makes a remark that
. . . if x is replaced by the complex frequency variable s, then the realization would
require a negative resistance. Thus, the [Khovanskii’s] CFEs do not seem to be useful
for realization of fractional inductor or capacitor.
Then he describes a method for circumventing this difficulty, which gives a continued fraction
expansion with positive coefficients.
However, the possibility of realization of negative impedances in electric circuits has been
pointed out by Bode [19, chapter IX]. Later, in 1970s, operational amplifiers appeared, which
significantly simplified creation of circuits exhibiting negative resistances, negative capa-
citances, and negative inductances. Such circuits are called negative-impedance converters
[20].
The simplest scheme of a negative-impedance converter (or current inverter) is shown in
Figure 10. The circuit consists of an operational amplifier, two resistors of equal resistance R,
and a component with the impedance Z. The entire circuit, considered as a single element,
has negative impedance −Z. This means that Iin = Vin /(−Z)).
For example, taking a resistor of resistance RZ instead of the element Z, we obtain a circuit,
which behaves like a negative resistance −RZ . The negative resistance means that if such an
element of negative resistance, for instance, −10 k: is connected in series with a classical
20 k: resistor, then the resistance of the resulting connection is 10 k:.
Let us now recall Example 8. Using negative-impedance converters, it is possible to design
a circuit with the required impedance Z(s), which will contain a negative capacitance C4 and
a negative inductance L3 .
Analogue Realizations of Fractional-Order Controllers 293

Figure 10. Negative-impedance converter.

Figure 11. Analogue fractional-order Iλ controller.

7. Example: Fractional-Order Iλ Controller

For experimental measurement we built a fractional-order Iλ controller which is a particular

case of the PIλ Dµ controller, (if K = 0 and Td = 0). The controller was realized in the
form of the finite domino ladder (n = 12), connected to feedback in operational amplifier
(Figure 11). It should be noted that the described methods work for arbitrary orders, but
the circuit elements with computed values are not usually available. Because of this, in our
experiment we proposed and realized the integrator with order λ = 0.5. It should be mentioned
that this simple case of the controller order can be realized also using the methods described
in [21–24], which do not involve explicit rational approximations.
In the case, if we will use identical resistors (R-series) and identical capacitors (C-shunt)
in the FDL, then the behaviour of the circuit will be as a half-order integrator/differentiator.
We used the resistor values R = 1k: (Rj = R, j = 1, . . . , n) and the capacitor values
C = 1µF (Cj = C, j = 1, . . . , n). For better measurement results we used two operational
amplifiers TL081CN in inverting connection.
A block diagram of the analogue fractional-order Iλ controller realization is shown in
Figure 11.
The resistors R1 and R2 are R1= R2 = 22 k:. The integration constant Ti can be
computed from relationship Ti = 1/ R/(Ri2 ∗ C), and for Ri = 22 k: we have Ti = 1.4374.
The transfer function of the realized analogue fractional-order Iλ controller is
Gc (s) = 1.4374 s−0.5 . (30)
Adjustment of the integration constant Ti of the fractional-order Iλ controller depicted
in Figure 11 was done by resistor Ri . If we change the resistor Ri , the integration constant
changes the value in the required interval.
In Figures 12 and 13 the measured characteristics of realized analogue fractional-order Iλ
controller are presented. In Figure 12 Bode plots a shown, and in Figure 13 is the time response
294 I. Podlubny et al.

Figure 12. Bode plots of the I1/2 controller (measurements).

Figure 13. Time response of the I1/2 controller to unit step input (measurements).

to the square input signal (unit step). We used frequency 100 Hz and amplitude ±10 V. It can
be seen from Figure 12 that the realized analogue of fractional-order Iλ controller provides a
good approximation in the frequency range [102 rad/sec, 5 · 102 rad/sec].
Measurements were done using IWATSU Digital Storagescope DS-8617 100 MHz, Hew-
lett Packard 35670A dynamic signal analyzer, Hewlett Packard 33120A 15 MHz function/
arbitrary waveform generator, power supply Thurlby-Thandar PL320QMD.
Analogue Realizations of Fractional-Order Controllers 295

8. Conclusion

In this paper we have demonstrated that the suggested use of continued fraction expansions
is a good general method for obtaining analogue devices (fractances) described by fractional
differential equations or by fractional-order transfer functions. Moreover, this approach can
be used for realization of other types of systems with transcendental transfer functions, which
can be developed in continued fractions. Furthermore, it has been shown that any rational
approximation of the transfer function can be used for designing the corresponding analogue
circuit, even if some of the coefficients of the resulting continued fraction are negative.
We have also introduced two types of nested multiple-loop systems, which can be easily
used for modelling, simulation, and realization of fractional-order systems and controllers,
and more generally for modelling, simulation and realization of systems, for which a rational
approximation of the transfer function can be obtained.
The exposition has been illustrated with several examples, including analogue realization
of an Iλ controller, for which experimental results were presented.

References

1. Miller, K. S. and Ross, B., An Introduction to the Fractional Calculus and Fractional Differential Equations,
Wiley, New York, 1993.
2. Oldham, K. B. and Spanier, J., The Fractional Calculus, Academic Press, New York, 1974.
3. Podlubny, I., Fractional Differential Equations, Academic Press, San Diego, CA, 1999.
4. Samko, S. G., Kilbas, A. A., and Maritchev, O. I., Integrals and Derivatives of the Fractional Order and
Some of Their Applications, Nauka i Tekhnika, Minsk, 1987 [in Russian].
5. Podlubny, I., ‘Fractional-order systems and PIλ Dµ -controllers’, IEEE Transactions on Automatic Control
44, 1999, 208–214.
6. Podlubny, I., ‘Fractional-order systems and fractional-order controllers’, UEF-03-94, Slovak Academy of
Sciences, Kosice, 1994.
7. Dorf, R. C. and Bishop, R. H., Modern Control Systems, Addison-Wesley, New York, 1990.
8. Oustaloup, A., Systèmes Asservis Linéaires d’Ordre Fractionnaire: Théorie et Pratique, Editions Masson,
Paris, 1983.
9. Oustaloup, A., La Dérivation non Entière, Hermès, Paris, 1995.
10. Oustaloup, A., Levron, F., Mathieu, B., and Nanot, F. M., ‘Frequency-band complex noninteger differenti-
ator: Characterization and synthesis’, IEEE Transactions on Circuit and Systems – I: Fundamental Theory
and Application 47, 2000, 25–39.
11. Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., Numerical Recipes in C. The Art of
Scientific Computing, 2nd edition, Cambridge University Press, Cambridge, 1992.
12. Dutta Roy, S. C., ‘On the realization of a constant-argument immitance of fractional operator’, IEEE
Transactions on Circuit Theory 14, 1967, 264–374.
13. Vinagre, B. M., Podlubny, I., Hernandez, A., and Feliu, V., ‘On realization of fractional-order controllers’,
in Proceedings of the Conference Internationale Francophone d’Automatique, Lille, Jule 5–8, P. Borne, J.-P.
Richard, and Ph. Vanheeghe (eds.), 2000, pp. 945–950.
14. Carlson, G. E. and Halijak, C. A., ‘Approximation of fractional capacitors (1/s)1/n by a regular Newton
process’, IEEE Transactions on Circuit Theory 11, 1964, 210–213.
15. Matsuda, K. and Fujii, H., ‘H∞ –optimized wave-absorbing control: analytical and experimental results’,
Journal of Guidance, Control, and Dynamics 16, 1993, 1146–1153.
16. Heymans, N. and Bauwens, J.-C., ‘Fractal rheological models and fractional differential equations for
viscoelastic behavior’, Rheologica Acta 33, 1994, 210–219.
17. Kvasil, J. and Čajka, J., An Introduction to Synthesis of Linear Circuits, SNTL/ALFA, Prague, 1981 [in
Czech].
18. Khovanskii, A. N., The Application of Continued Fractions and Their Generalizations to Problems in
Approximation Theory, Noordhoff, Groningen, 1963.
296 I. Podlubny et al.

19. Bode, H. W., Network Analysis and Feedback Amplifier Design, Tung Hwa Book Company, Shanghai, China,
1949.
20. Dostal, J., Operational Amplifiers, Butterworth-Heinemann, Boston, MA, 1993.
21. Nakagava, M. and Sorimachi, K., ‘Basic characteristics of a fractance device’, IEICE Transactions
Fundamentals E75-A, 1992, 1814–1818.
22. Oldham, K. B. and Zoski, C. G., ‘Analogue instrumentation for processing polarographic data’, Journal of
Electroanalytical Chemistry 157, 1983, 27–51.
23. Petráš, I. and Dorčák, L’., ‘Some possibilities for realization of fractional order controller’, Envirautom 4,
1999, 83–90 [in Slovak].
24. Wang, J. C., ‘Realizations of generalized Warburg impedance with RC ladder networks and transmission
lines’, Journal of Electrochemical Society 134, 1987, 1915–1920.