Abstract

One of the most important issues in application of noninteger order systems concerns their implementation. One of the possible approaches is the approximation of convolution operation with the impulse response of noninteger system. In this paper, new results on the Laguerre Impulse Response Approximation method are presented. Among the others, a new proof of convergence of approximation is given, allowing less strict assumptions. Additionally, more general results are given including one regarding functions that are in the joint part of and spaces. The method was also illustrated with examples of use: analysis of “fractional order lag” system, application to noninteger order filters design, and parametric optimization of fractional controllers.

1. Introduction

Noninteger order (fractional) systems are becoming more prevalent in control and signal processing applications. The fact that they cannot be directly implemented on digital or analog platforms is, however, a known issue. The issues come from difficulties with realization of noninteger derivatives and lack of semigroup property. The search for efficient approximation that will keep the beneficial properties of fractional system while allowing implementation is a topic of ongoing research.

In this paper, we discuss new developments in the Laguerre Impulse Response Approximation (LIRA) method. This method allows efficient approximation of wide class of noninteger order systems, in particular noninteger order filters. Those new developments include formal proofs of convergence in both general and particular cases and proofs of important properties needed for implementation.

Besides LIRA, there are also different methods of noninteger order system approximation. They can be divided into three groups. The first group is based on approximating operator in the frequency domain. The most popular approaches are Oustaloup’s method [1, 2] and CFE (Continuous Fraction Expansion) method [35]. These methods are based on different premises but both allow obtaining relatively easy approximations. CFE method is generally considered as worse of the two in the aspect of frequency response representation [6]. Oustaloup’s method gives a very good representation of frequency response at selected band at the cost of high numerical sensitivity. This sensitivity can be lessened with use of time domain realizations [7, 8]. Different approach to approximation uses discrete realizations based on truncation of series representations. In particular, this class of methods includes truncation of Grünwald-Letnikov derivative definitions and power series expansions (PSE) in variable (discrete frequency) domain. These methods give good approximations only at high frequencies. Examples of the method can be found in [9]. Improvements on PSE methods with better band of correct approximation were investigated by Ferdi [10, 11]. These two approaches, while being fundamentally different, have in common the fact that they are potentially effective but there are no proofs of their formal convergence to the actual noninteger order system. Another group of approximation methods uses diffusive realization of pseudodifferential operators [12] for approximation. First works in the area used finite-dimensional approximation with trapezoidal integration [13]. Later works used simple diagonal matrix time domain realization with modified Oustaloup nodes from frequency domain method [1416]. This method is especially useful in analysis of infinite-dimensional systems using operator theory. It can be used, for example, to prove such properties as stability of closed-loop system [17].

The method of approximation analyzed in this paper is the Laguerre Impulse Response Approximation (LIRA) which is based on approximation of noninteger order system impulse response with Laguerre functions. Early works in this area were unsuccessful in approximation of order integrators [18, 19]. In [20], it was shown that under certain assumptions method is convergent in and norms when approximating asymptotically stable noninteger transfer functions of relative degree equal to or greater than one. Most applications of the method can be found in filtering and parametric optimization [2123]. In [24], it was shown that in certain cases and at low orders of approximation LIRA outperforms Oustaloup method.

The main contribution of this paper is Theorem 6 which gives the justification of approximation of impulse responses that are in . It was analyzed earlier in [20]; however, proof given in that paper was based on analysis of integrals and was not fully precise. Moreover, there was an assumption that impulse response has to be bounded which in the light of results of this paper is not necessary. In addition to that, Theorem 6 can be considered in spaces of multivariable functions. Additional contribution is Theorem 2, which describes properties of functions in . There is also series of corollaries and lemmas which allows explicit construction of approximation method. Methods’ operation is illustrated with examples. The first example is qualitative where we explicitly verify the assumptions and other two come from applications: approximation of a filter and controller parameter optimization.

The rest of this paper is organized as follows. The first section contains the preliminaries necessary for theorems and lemmas presented further in the paper. In the main section, we prove the theorems concerning some properties of chosen types of functions and the main result concerning approximation in space. Then explicit formulas for Laguerre Impulse Response Approximation method are given along with three examples: “fractional order lag” system, noninteger order filtering, and parametric optimization of noninteger order closed-loop system. The paper ends with conclusions and proposition of further extensions of the method.

2. Preliminaries

In this section, we present necessary definitions and notations that will be used throughout the paper.

We focus on stable noninteger order systems given by transfer functions likewith zero initial conditions. Our interest is focused especially on their time domain representation in the form of convolution operator; that is,If realization in the form of noninteger differential equations is needed, we will, without the loss of generality, use Caputo definition of noninteger order derivative: where denotes the gamma function

We will consider function spaces and of real functions , where is an open subset of . Space is a Banach space of absolute integrable functions with the norm Space is a Hilbert space of square integrable functions, with norm induced by the scalar product Additionally, we will consider the Banach space of essentially bounded functions and a sequence space with norm

We will denote by ,, a set of orthonormal basis functions in . Additionally, for space , we will consider orthonormal Laguerre functions denoted by and given bywhere is an arbitrary positive constant and is th Laguerre polynomial of the form

3. Main Results

In this section, we will prove the theorem regarding convergence of the function approximation and its consequences that can be used for approximation of convolution operator. In order to do so certain assumptions and lemmas have to be stated.

We will start with proving an important result regarding functions in .

Remark 1. In this paper, we will take the domain of a quotient of two functions and as a set, where exists.

Theorem 2. The following statements are equivalent: (1).(2)There exists a function such that .

Proof. .
Let be defined as where is an arbitrary constant. Then is in . As , we have that .
.
We have so . Using Hölder’s inequality yields also the norm of in space: Hence, .

Remark 3. As an example of we can take bounded function in .

In order to prove the main result, we will need the following two lemmas.

Lemma 4. If there exists a function , such that and , , thenwhere and and .

Proof. As , we can write function as a series expansion of orthonormal basis functions in : with coefficients . Therefore, we have The same operation is done for and , since both are in as proven in Theorem 2. We have where and where . On the other hand, we haveLet us take partial sum of series (20). As it is convergent in , also norm of its partial sum is finite. We have thenAs norm in of left side of (21) is finite, we have which implies thesis.

Lemma 5. Let ,, be a family of sequences parametrized by , such that (1): when ,(2): .Then

Proof. Let us consider . The sequence is square-summable with respect to for all . Therefore, we can find a function such that Using the first assumption, we have From assumption (1), we know that we can choose , such that , for all and that Therefore, for any , there is , such that for all .

Now we can formulate and prove the intended theorem.

Theorem 6. If there exists a function , such that and , where are functions of orthonormal basis in , then

Proof. From Theorem 2 and assumptions, we know that , , and are in and as such may be written in form of infinite series: with coefficients , , and , respectively. Let The following estimate holds:Let us denote as a sequence with elements . From (32) and orthonormality of the basis , we have We will investigate the behavior of the sequence . Using linearity of inner product in , we have hence, Using the Schwartz inequality, we have and (36) imply that when , for all . From Lemma 4 sequence is square-summable for all . From properties of series expansion in , we can conclude that is also square-summable. As a consequence, the sequence is square-summable. The assertion of the theorem follows from inequality (33) and Lemma 5.

Remark 7. Although Theorem 2 guarantees the existence of for every function in , it does not ensure that for a given basis. Therefore, the choice of appropriate basis functions determines the existence of approximation.

As a direct consequence of Theorem 6, we have the following results.

Corollary 8. If and there exists a basis , such that for , with given by (11), then

Proof. From Theorem 2, we know that there exists such that . Using the assumption that , we have . So the assumptions of Theorem 6 are fulfilled which completes the proof.

Corollary 9 (convergence of Laguerre series). If and are orthonormal Laguerre functions, then where

Proof. Because for from Corollary 8 we have for , the assumptions of Corollary 8 are fulfilled.

Corollary 10 (convolution approximation in Laguerre basis). Let and , then for where with given by (39), we have

Proof. Using Hölder’s inequality, we have Applying Corollary 9 ends the proof.

4. Laguerre Impulse Response Approximation Method

Corollaries 9 and 10 give theoretical basis for approximating convolution operator (2) with a Laguerre series. In the following we present three lemmas that allow practical construction of approximant for systems described with SISO (single-input single-output) noninteger transfer functions.

Lemma 11 (computation of scalar products). If has a Laplace transform , then where is the th orthonormal Laguerre function with parameter :

Proof. Proof is given in [25, page 65].

Lemma 12 (choice of optimal ). Let and where and is the th orthonormal Laguerre function with parameter . The norm of the difference for given is minimal if

Proof. We will compute the square of error norm (arguments are dropped): It can be easily seen that the above expression is minimal when thesis is fulfilled.

Lemma 13 (convolution state space realization). The convolution of with given by where where is the th orthonormal Laguerre function with parameter , is equal to the solution of the following system of differential equations: with .

Proof. It can be easily seen that Let us set which along with (53) gives (52). From boundedness of and definition (54), we get .
We will now use two additional facts. First, the derivative of a Laguerre function with respect to first argument can be expressed with Laguerre functions of the same and lower order; that is, Additionally, the derivative of the convolution operator with respect to is given by Differentiating (54) and using (55) and (56), we obtain Collecting (57) for leads to (51).

Lemma 11 enables efficient computation of approximation coefficients. Lemma 12 shows how to choose the basis parameter allowing minimal approximation error for given approximation order. Finally, Lemma 13 gives a tool how to construct a state space (or if needed transfer function) realization of the approximant allowing among the others simulation, filtering, use in real time applications, or optimization.

We will now present three examples of using the proposed method for approximation of noninteger order systems.

4.1. “Fractional Order Lag”

This class of systems is widely used in various applications, usually without serious analysis. These systems are classified as noninteger order; however, their realization is in the form of noninteger differential equations (at least not their finite number). Their transfer function takes the formwhere is a constant. Impulse response for (58) can be found using inverse Laplace transform:where denotes gamma function [26]. It can be easily verified that (63) for is bounded and its norm can be easily computed. The impulse response is positive for all ; therefore, the absolute value in norm can be omitted. ConsiderTherefore and it is sufficient to use initial and final value theorems: Because of that and Remark 3, we can see that .

However, it can be shown that while for the impulse response is unbounded the assumptions of Theorem 6 are fulfilled.

Corollary 14. Approximation of function with Laguerre series is convergent in the sense of norm.

Proof. We will prove that function (67) fulfills the assumptions of Theorem 6.
First we will prove that . The following calculations are valid: where is the “upper” incomplete gamma function and .
The limit of (68) for can be written asand is bounded for ; that is, . The same analysis was conducted for the second limit :From formulas (70) and (71), one obtains that the norm of is bounded for .
Next, it can be easily seen that setting for gives and , for and , which completes the proof.

We will now present examples of approximation for different “fractional order lags,” one for each case. In this example and in the next example, we will analyze three types of plots. We will compare impulse responses of original noninteger order system with one of its approximants (with chosen according to Lemma 12) and with an approximant with nonoptimal . Analogous comparison will be performed with bode plots, which are a good method for transfer function evaluation. Finally, we will present the function with marked optimal and nonoptimal values of . In case of fractional order lag approximations presented are of order 20.

In Figures 1(a), 2(a), and 3(a), we analyze fractional order lag with and . As it can be seen in Figure 1(a), while both approximants are bounded, they rather closely approximate the unbounded impulse response. Differences with optimal and nonoptimal are especially visible in Figure 2(a), where magnitude and phase start to diverge for higher frequencies. It should be noted that phase approximation is strongly oscillatory. As it can be seen in Figure 3(a), function maximized with respect to is unimodal; however, the more complicated system becomes the function changes.

In Figures 1(b), 2(b), and 3(b), we can see analogous analysis for fractional order lag with and . In Figure 1(b), the differences between optimal and nonoptimal are more visible. Frequency responses for optimal are also better and less oscillatory (see Figure 2(b)). Finally, the analysis of optimization shows much more irregular shape of maximized function.

More detailed analysis of Laguerre Impulse Response Approximation method for fractional order lag can be found in [23].

4.2. Noninteger Order Filter

The main area of application of LIRA method investigated by authors was the realization of noninteger order filters. Our focus was mostly on generalizations of second-order filters developed by Radwan et al. [27]. Detailed analysis of method operation was presented in [24]. Here we will present only a brief description of the idea, with an example of low-pass filter. We consider a filter in a form where and fulfill the conditions of stability (see [28, 29]).

Remark 15. We do not have at the moment condition when the noninteger transfer function has an impulse response in . At the moment working hypothesis is that sufficient conditions are stability and transfer function’s relative degree (difference between max. degree of denominator and max. degree of numerator) is at least 1. Despite multiple trials, no counter example was found.

Equivalent representation of (73) is the realization in the form of a system of differential equations of order . This system can take the form (see [30])with the following matrices: The impulse response of filter (73) is given by [30] where is the two-parameter Mittag-Leffler function.

In order to illustrate the operation of LIRA for filter we have chosen one, for which poles are a conjugate pair in right complex half plane. Such filters are oscillatory and always have resonance peak. We compare both impulse responses and frequency responses. All approximations are of 25th order with optimally and nonoptimally chosen . The case is presented in Figure 4. As it can be seen in Figure 4(a), the optimal approximant of impulse response is almost indistinguishable from the exact value. Nonoptimal one well approximates initial behavior but for longer time stops at zero. In the frequency response (Figure 4(b)) the optimal approximation is again almost indistinguishable from noninteger order filter. However, nonoptimal causes oscillatory behavior on magnitude and phase responses. Moreover, phase response exhibits bias. It can be seen in Figure 4(c) that the difference between optimal and nonoptimal is on the level of 10% of optimized function value.

4.3. Optimization Problem

The third example chosen to show the use of approximation method consists of optimization of noninteger order closed-loop system. In contrast to classical integer order systems, the solution of noninteger order system yields problems with either exact form of the solution or its numerical implementation as it requires the whole “history” beginning with initial conditions (lack of semigroup property). Moreover, even though there are exact solutions presented, for example, in [26], it is worth noticing that they require special functions, that is, Mittag-Leffler function. Furthermore, taking into account that is a multivalued function for complex and noninteger , we can notice that direct calculation of integral of impulse response is complicated. Therefore, an approximation method is indispensable for effective optimization.

The main purpose is to find optimal controller parameters, , , and , in order to efficiently bring the closed-loop system to zero. For this purpose, we chose the following objective function where denotes the impulse response of the system. It can be observed (see, e.g., [21, 24]) that Thus, the optimization consists in minimizing the performance index (77) with respect to three parameters of controller (80). It is important to notice that during the procedure both types of parameters are optimized, and controller parameters.

Other examples of performance indicators were analyzed, for example, in [31].

4.3.1. Optimization Algorithm

The algorithm consists of two main steps described below:(1)Set starting points for controller parameters in such a way that the noninteger order system is asymptotically stable.(2)Perform the optimization routine calculating the chosen performance indicator.(a)For given values of parameters of the controller perform the stability check of the noninteger order system. If the parameters lie outside the asymptotic stability region, set the value of performance indicator to infinity and go to the third step. If the system is asymptotically stable, then go to second step.(b)In order to approximate the stable system with LIRA method, find current optimal for given controller parameters. If is infinite in this iteration, set the performance indicator to infinity. Otherwise, calculate the value of performance indicator using the approximation method.(c)On the basis of value of performance indicator, find the controller parameters for next step according to chosen optimization routine.

It is worth noticing that although the optimization is performed on approximated system, the stability test is performed for noninteger order one. To verify the stability, we used an original algorithm based on D-partition method for noninteger order systems (see, e.g., [32]).

4.3.2. Closed-Loop System

As an example, we analyzed an integer order system from Figure 5 with transfer function and noninteger order controller The transfer function of a closed-loop system is The poles of system (79) are real and negative. The similar system was widely analyzed in literature (see, e.g., [33]) in form of RC ladder system. Among the others, it is very useful as it allows approximating some distributed systems, for example, transmission line [34].

4.3.3. Optimization with respect to , , and

The system was analyzed in, for example, [22]. The parameters of approximation are highly nonlinear with respect to and therefore the performance indicator is not necessarily strictly convex.

The optimal values are , , and for . In Figure 6, we depicted the results in the form of impulse response of original system and the system with the controller.

4.3.4. Analysis of the Method

The approximation method proved to be useful for tuning noninteger order controllers. It allows approximating noninteger order system with an asymptotically stable system of integer order differential equations. Therefore, it is possible to use classical theory, especially for calculation of performance indicators. Although the approximation is stable, it does not ensure stability of noninteger order system.

5. Conclusions

In this paper, we have shown theoretical basis and actual effectiveness of LIRA method for approximating noninteger order systems. As it can be seen (documented in earlier works, i.e., [2124, 35]), it is an approximation of input-output relationships, adequate for such purposes like filtering or control loop optimization. The main contribution of this paper over the earlier results was relaxation of the assumptions regarding the impulse response of approximated system, as it is no longer required for it to be bounded.

The main weakness of the method lies in computation of approximation coefficients . While formula for computation of scalar products in Lemma 11 is useful for any order of approximation, it still requires computation of transfer function derivatives. The best results with respect to precision come from symbolic computation (particularly in Maple, MATLAB Symbolic toolbox gives worse results when converting from symbolic to numeric). It is, however, very time consuming, which is especially troublesome during optimization. Different approach is to use Leibnitz formula for derivative of th order described in [20]; unfortunately, this formulation, for high , introduces multiple factorials of high numbers which leads to numerical errors caused by factorial approximation (it was shown in [23]). Also changes in transfer function for this approach lead to rather advanced code modifications. Difference approximations of derivatives are too sensitive to be considered at all. One could consider computing scalar products directly from definition but it requires the time domain form of impulse response, which is usually either unavailable or expressed with special functions.

The other less severe problem is that there are still no definitive criteria on when an impulse response of noninteger order system is in . The working hypothesis that is described in Remark 15 even if proved will be only a sufficient condition. It was observed in [24] that approximation can be convergent also for relative degrees less than one.

Further development of the method will include addressing the mentioned flaws. However, the main interest of the authors is in applications. At the moment, we work on EEG signal filtering, vibration signal filtering, optimization of noninteger order controllers, and development of circuit realizations of noninteger order systems.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The work is realized in the scope of project titled “Design and Application of Non-Integer Order Subsystems in Control Systems.” Project was financed by National Science Centre on the base of Decision no. DEC-2013/09/D/ST7/03960.