New Stability Tests for Discretized Fractional-Order Systems Using the Al-Alaoui and Tustin Operators

Stanisławski, Rafał; Rydel, Marek; Latawiec, Krzysztof J.

doi:https://doi.org/10.1155/2018/2036809

Complexity

On this page

Abstract Introduction Preliminaries Conclusion Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 2036809 | https://doi.org/10.1155/2018/2036809

New Stability Tests for Discretized Fractional-Order Systems Using the Al-Alaoui and Tustin Operators

Rafał Stanisławski,¹Marek Rydel,¹and Krzysztof J. Latawiec¹

Academic Editor: Rongqing Zhang

Received17 May 2018

Accepted18 Sept 2018

Published07 Nov 2018

Abstract

This paper provides new results on a stable discretization of commensurate fractional-order continuous-time LTI systems using the Al-Alaoui and Tustin discretization methods. New, graphical, and analytical stability/instability conditions are given for discrete-time systems obtained by means of the Al-Alaoui discretization scheme. On this basis, an analytically driven stability condition for discrete-time systems using the Tustin-based approach is presented. Finally, the stability of discrete-time systems obtained by finite-length approximation of the Al-Alaoui and Tustin operators are discussed. Simulation experiments confirm the effectiveness of the introduced stability tests.

1. Introduction

Stable discretization of continuous-time fractional-order systems is an important issue in various areas of science and technology, including system science, signal processing, and control theory. In this field, we have three main discretization operators which can generate discrete-time counterparts for continuous-time fractional-order systems, in terms of the Euler, Tustin, and Al-Alaoui methods. There are two main problems to be solved during the discretization process for fractional-order systems. Firstly, the three discretization schemes lead to infinite complexity of rational, discrete-time counterparts of fractional-order derivative. Therefore, in practical applications various finite-length approximations of the discretization operators have been used, involving the most popular finite fractional difference (FFD) approximation in the Euler approach [1, 2] and finite-length implementations of the continuous fraction expansion (CFE) method in the Tustin and Al-Alaoui approaches [2–6]. Also, a number of papers have presented some other approximation/discretization methods for the fractional-order derivative [7–10].

Secondly, it is well known that the discretization process affects stability conditions for the discrete-time counterparts of continuous-time fractional-order systems [11–14]. Moreover, the stability conditions of fractional-order discrete-time systems depend on a type of operator used in the discretization process. This can be easily seen when we compare stability results for discrete-time systems obtained for the “forward-shifted” Euler operator [15, 16] with those for the classical backward Euler operator [13]. The first results in the area of stability analysis for discrete-time fractional-order systems have been developed in [11], and they concern sufficient conditions only. More complete results have been obtained in a special case of discrete-time fractional-order positive systems [17–20]. Simple, analytical, necessary, and sufficient stability results for discrete-time systems have been obtained for both “forward-shifted” and backward Euler operators [13, 15, 16]. Specific, numerical stability results for discrete-time fractional-order systems based on the Euler expansion have been presented in [21–23]. On the other hand, it is well known that finite-length implementations of the Euler and Tustin operators affect the stability conditions for the underlying discrete-time systems [2, 11, 13, 15].

In this paper, we introduce simple, either analytically driven or purely analytical stability tests for discrete-time systems obtained by the use of the Al-Alaoui operator. These results are then used to propose a stability test for systems based on the discretized Tustin operator, which can be considered as a special case of the Al-Alaoui approach. Also, practically oriented results for discretization using a finite-length implementation of the Tustin approach of [2] are extended to finite-length approximation using the Al-Alaoui method [24].

This paper is organized as follows. Having introduced in Section 1 the stability problem for discretized commensurate fractional-order continuous-time systems, the system representations based on the Euler, Tustin, and Al-Alaoui discretizers are given in Section 2. Section 3 presents new stability results involving both graphical and analytical criteria, which are the main results of the paper. Moreover, Section 3 adopts the analytical stability criterion for the Al-Alaoui operator to the Tustin-based one. Discussion on the stability of discretization based on finite-length approximations of the Al-Alaoui operator is presented in Section 4. Simulation examples of Section 5 confirm the effectiveness of the proposed stability results. Section 5 summarizes the achievements of the paper.

2. Preliminaries

Consider a linear continuous-time state space system of commensurate fractional order described by where denotes a fractional-order derivative of order ; , , are the state, control, and output vectors, respectively; and , , and are the system matrices, with and being the number of inputs and outputs, respectively. Here, the fractional-order derivative will be described by three various representations, involving the Caputo, Riemann-Liouville, or Grünwald-Letnikov definitions. But regardless of the specific definition of the fractional-order derivative, assuming the zero initial conditions in (1), that is , for any , the Laplace transform of system (1) is as follows:

In the specific SISO case, when , the system of (2) can be described by the transfer function where and are the coprime polynomials in the variable and , , and , , are the poles and zeros of , called -poles and -zeros, respectively (compare [16]). Note that the -poles , , are the eigenvalues of the state matrix .

In the discretization process, we seek for a discrete-time equivalent of the fractional-order system (2), in form of the -transform where will be used as a discrete-time model of . Alternatively, for the SISO case, we can obtain a discrete-time fractional-order transfer function in form of with , , and , , being the -poles and -zeros, respectively, of the system defined in (3). Three main discretizers are used here as discretization functions , namely, the Euler, Tustin, and Al-Alaoui operators. The Euler operator is used in two versions, that is, the backward Euler operator or the “forward-shifted” Euler one with being the sampling period. The Tustin operator has the form

The Al-Alaoui operator [25, 26] is obtained through an assumption that the integration rule is realized as a weighted sum of the Tustin and backward Euler rules. This leads to the Al-Alaoui operator [25, 26] where is the weighting coefficient.

It has been presented in several papers that selection of the discretization operator in the discretization process affects the stability conditions for discrete-time systems [11, 13, 16]. For instance, the two Euler methods represented by (6) and (7) lead to quite different stability results for the discrete-time systems. A detailed stability analysis for discrete-time systems obtained by the use of the two Euler operators has been presented in [13, 15, 16].

In this paper, a new stability analysis for systems obtained by the discretization process using the Al-Alaoui operator will be presented.

3. Main Results

Firstly, we introduce a graphical stability approach, based on the Al-Alaoui operator.

Theorem 1. The discrete-time fractional-order system (4) or (5), with as in (9), is asymptotically stable if and only if all -poles of the system, , are outside the instability area where is the fractional-order of system (4) or (5), is the weighting coefficient, and is the argument of the function .

Proof. Consider the discrete-time fractional-order system (4) or (5) with as in (9). Taking into account that the Al-Alaoui operator of (9) is meromorphic, the stability area with respect to -poles depends on the contour of function , [13, 15]. Therefore, the bound contour of the stability/instability regions is as follows: where , , and . Now, substituting and of (12) into the circle equation with the center in the point , we obtain So, for , the inner contour fulfills the circle equation with the center in and radius . Therefore, we can write the same inner contour as a function of another parameter , that is, and finally we can describe the stability/instability contour for the Al-Alaoui operator of (12) as a function of as The next steps of the proof follow the lines of Theorem 1 of [13]. Note that the characteristic equation of system (4) or (5) is asymptotically stable if and only if all elements , , do not generate unstable poles in the -domain. Accounting that is a meromorphic function ( is homomorphic in ) and the contour presented in (16) is a simply closed contour in the complex plane, we can prove the stability of the elements , , on the basis on Cauchy’s argument principle and Rouché’s theorem [27].
Now, consider the simple closed contour as in Figure 1 whose interior domain is such that and . Then, the winding number of about , , that counts the number of times the curve winds around the point is where and and denote the numbers of zeros and poles of inside the contour , respectively. Since the transformation has a single pole in 0 (see (9)), the coefficient is 0. Therefore, the winding number will describe a number of unstable elements of . The contour is based on (16) and is presented in Figure 2. Now:
() Assume that the system is asymptotically stable. Then, all poles for all elements , , are outside the interior domain of the contour . Therefore, and all , , are outside the area bounded by , which can be described by (10).
() Assume that all eigenvalues , , are outside the instability area (10). Then, all eigenvalues are outside the contour and the winding number . Therefore, the elements , , do not generate poles in interior domain of the contour . So the system is asymptotically stable. This completes the proof.

A graphical presentation of the stability/instability areas is shown in Figure 3.

It can be seen from Figure 3 that the stability area depends on the weighting parameter of the Al-Alaoui equation. For , we obtain the instability area as for the backward Euler operator (compare [13]) and increasing leads to enlargement of the instability area.

On the basis of Theorem 1, we can present a new analytical result as follows.

Theorem 2. The discrete-time fractional-order system (4) or (5) with described by the Al-Alaoui operator (9), with , is not asymptotically stable if and only if there exists any where and are the modulus and argument, respectively, of the -pole of the system and being the Al-Alaoui weighting coefficient.

Proof. Taking into account that the stability/instability bounding contour is given by (15), which can be presented in form of and introducing the resultant argument of , that is, , we obtain Note that (21) presents the modulus of as a function of its argument . Also note that for , the function is redefined in . For the -poles of the system with , , the -poles are outside the instability area (see the stability/instability areas in Theorem 1) and the system is asymptotically stable.
() Assume that the system is not asymptotically stable. Then, at least one , , lies inside or in the bound of the instability area described in Theorem 1. Since the bound function is redefined in , the argument of unstable eigenvalue has to be . For , , the -poles are inside the instability area when the modulus of the th -pole is lower (or equal) than the modulus of the stability/instability bounding contour presented in (21) for the same angle . Therefore, there is at least one -pole passing condition (18).
() Assume that an -pole , , of the system passes condition (18). Then, taking into account the stability/instability bound equation (21), the -pole is inside (or in the bound) of the instability area. Taking into account Theorem 1, the system is not asymptotically stable. This completes the proof.

On the basis of Theorem 2, we can immediately present the following.

Theorem 3. The discrete-time fractional-order system (4) or (5) with described by the Al-Alaoui operator (9), with , is asymptotically stable if and only if where and are the modulus and argument, respectively, of the -pole of the system and being the Al-Alaoui weighting coefficient.

Proof. Immediate from Theorem 2.

The results of Theorems 2 and 3 for the zero weighting coefficient of the Al-Alaoui operator () specialize to those for the backward Euler-based discretization scheme presented in [13]. Also, on the basis of the results, we can easily present the stability criterion for the Tustin-based discretization scheme.

Theorem 4. The discrete-time fractional-order system (4) or (5) with described by the Tustin operator (8), with , is asymptotically stable if and only if where is the argument of the -pole of the system.

Proof. Consider the stability Theorem 3. Accounting that the Tustin-based approach is a special case of the Al-Alaoui approach with , we arrive at the following modulus condition: Taking into account that then we have , , and we arrive at condition (23).

Remark 1. Note that the results of Theorem 4 are the same as those for continuous-time systems [28, 29]. Thus, using the Tustin-based discretization scheme does not affect the stability conditions for discrete-time systems. However, the Tustin-based discretization is inferior with respect to approximation accuracy [10].

Remark 2. The results of Theorems 2 and 3 show that discretization of a stable continuous-time system using the Al-Alaoui-based approach guarantees the asymptotic stability of a discrete-time system. Moreover, for the unstable continuous-time system, we can select such a sampling period and weighting coefficient that lead to a stable discrete-time system.

4. Finite-Length Implementation of Discretization Operators

It is important that the “ideal” Al-Alaoui and Tustin operators are not applicable in practice due to infinite-length implementation of the discretization equations. Therefore, in practical applications, we approximate the two above operators by use of rational finite-length discrete-time transfer functions. Usually, the approximations are obtained through the continuous fraction expansion (CFE) method, but in the Tustin-approach, the Muir recursion can also be applied [10]. Regardless of the approximation method, the rational discrete-time transfer function approximator to is obtained as follows:

In this case, the contour of the approximation in the complex plane is (see [2, 24])

Now, if curve (27) constitutes a simple closed curve in the complex plane, then the stability/instability areas with respect to , , are separated from each other by the contour. Exemplary stability/instability areas for , , , and various implementation lengths are presented in Figure 4.

5. Simulation Examples

Example 1. Consider a commensurate continuous-time fractional-order state space system as in (1) with and . The system has three -poles, i.e., and . Note that since , the system is unstable. The system is discretized by the use of the Al-Alaoui approach with the sampling period and four different values of the weighting coefficient , 0.6, 0.7, and 0.8. The stability properties of the discretized system designed by use of Theorem 3 are presented in Table 1. Moreover, Table 1 shows the stability results for the Tustin method as a special case of the Al-Alaoui approach with . The -poles with the stability/instability areas are depicted in Figure 5.
It can be seen from Table 1 that even though the continuous-time system is unstable, the Al-Alaoui method with the weighting coefficients and can lead to the stable discrete-time system, but for the weighting coefficients and , the discrete-time system still remains unstable. It is worth mentioning that the stability/instability of the discrete-time system is related to the unstable -poles and . Since the -pole is stable, the discrete-time counterpart will still be stable with respect to . The outcomes of Table 1 for the Tustin-based discretization scheme confirm the result of Theorem 4. Since the -poles and are outside the instability areas for and and inside the instability areas for and , the results of Figure 5 and Theorem 1 fully confirm the stability results specified in Table 1.

Example 2. Consider a commensurate continuous-time fractional-order state space system as in (1) with and whose eigenvalues are . Note that the system is unstable due to . The system is discretized by the use of the Al-Alaoui operator with the sampling period and two different values of the weighting coefficient and 0.8. Stability results for the discretized system are presented in Table 2.
It can be seen from Table 2 that the discrete-time system is stable for but unstable for .

Example 3. Consider fractional-order continuous-time system with and three -poles, i.e., and . Note that and the system are stable. The system is discretized with the sampling period and the approximation of the Al-Alaoui operator as in (26), with and different implementation lengths and 5
The stability areas and -poles of the system are presented in Figure 6.
It can be seen from Figure 6 that although the original system is stable, the discretized one by use of finite implementation of the Al-Alaoui operator is unstable for and stable for . Therefore, in contrast to the infinite-length Al-Alaoui operator, in some specific cases, the finite-length implementation of the Al-Alaoui discretizer for a stable continuous time system can lead to an unstable discrete-time system.

All the results presented in the above examples have been confirmed by a variety of BIBO stability experiments.

Remark 3. Matlab-scripted files for the above examples are available from the web: doi:10.5281/zenodo.1414176.

6. Conclusion

This paper has offered new, simple, graphical, and analytical stability/instability conditions for continuous-time commensurate fractional-order systems discretized by the use of the Al-Alaoui and Tustin operators. Firstly, theoretical stability condition for the Al-Alaoui operator has been given in a graphical way, which is then used in simple, analytical stability tests for both Al-Alaoui and Tustin approaches. Finally, the stability of discrete-time systems by use of finite-length approximations of the Al-Alaoui and Tustin discretizers has been discussed. Simulation examples fully confirm the original stability results obtained.

It is important that the fractional-order stability analysis presented in the paper is based on the eigenvalues of the state matrix, exactly as for the integer-order systems. This method is very useful; however, in case of high dimensions of the state matrices, calculation of eigenvalues may lead to numerical problems.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

I. Petráš, Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation, Nonlinear Physical Science, Springer, Berlin, Heidelberg, 2011.
View at: Publisher Site
M. Siami, M. S. Tavazoei, and M. Haeri, “Stability preservation analysis in direct discretization of fractional order transfer functions,” Signal Processing, vol. 91, no. 3, pp. 508–512, 2011.
View at: Publisher Site | Google Scholar
Y. Q. Chen and K. L. Moore, “Discretization schemes for fractional-order differentiators and integrators,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 49, no. 3, pp. 363–367, 2002.
View at: Publisher Site | Google Scholar
Y. Q. Chen and B. M. Vinagre, “A new IIR-type digital fractional order differentiator,” Signal Processing, vol. 83, no. 11, pp. 2359–2365, 2003.
View at: Publisher Site | Google Scholar
R. Stanisławski, M. Rydel, and M. Gałek, “Implementation issues in discretization of fractional-order derivative using the Al-Alaoui operator,” in 2016 21st International Conference on Methods and Models in Automation and Robotics (MMAR), pp. 844–847, Miedzyzdroje, Poland, 2016.
View at: Publisher Site | Google Scholar
B. M. Vinagre, Y. Q. Chen, and I. Petráš, “Two direct Tustin discretization methods for fractional-order differentiator/integrator,” Journal of the Franklin Institute, vol. 340, no. 5, pp. 349–362, 2003.
View at: Publisher Site | Google Scholar
P. Bania and J. Baranowski, “Laguerre polynomial approximation of fractional order linear systems,” in Advances in the Theory and Applications of Non-integer Order Systems, vol. 257 of Lecture Notes in Electrical Engineering, pp. 171–182, Springer, Heidelberg, Germany, 2013.
View at: Publisher Site | Google Scholar
J. Baranowski, W. Bauer, M. Zagórowska, and P. Piatek, “On digital realizations of non-integer order filters,” Circuits, Systems, and Signal Processing, vol. 35, no. 6, pp. 2083–2107, 2016.
View at: Publisher Site | Google Scholar
M. Bishehniasar, S. Salahshour, A. Ahmadian, F. Ismail, and D. Baleanu, “An accurate approximate-analytical technique for solving time-fractional partial differential equations,” Complexity, vol. 2017, Article ID 8718209, 12 pages, 2017.
View at: Publisher Site | Google Scholar
R. Stanisławski, K. J. Latawiec, and M. Łukaniszyn, “A comparative analysis of Laguerre-based approximators to the Grünwald-Letnikov fractional-order difference,” Mathematical Problems in Engineering, vol. 2015, Article ID 512104, 10 pages, 2015.
View at: Publisher Site | Google Scholar
A. Dzieliński and D. Sierociuk, “Stability of discrete fractional order state-space systems,” Journal of Vibration and Control, vol. 14, no. 9-10, pp. 1543–1556, 2008.
View at: Publisher Site | Google Scholar
S. Guermah, S. Djennoune, and M. Bettayeb, “A new approach for stability analysis of linear discrete-time fractional-order systems,” in New Trends in Nanotechnology and Fractional Calculus Applications, D. Baleanu, Z. B. Guvenc, and J. A. T. Machado, Eds., pp. 151–162, Springer, Dordrecht, Netherlands, 2010.
View at: Publisher Site | Google Scholar
R. Stanisławski, “New results in stability analysis for LTI SISO systems modeled by GL-discretized fractional-order transfer functions,” Fractional Calculus and Applied Analysis, vol. 20, no. 1, 2017.
View at: Publisher Site | Google Scholar
S. B. Stojanović and D. L. Debeljković, “Simple stability conditions of linear discrete time systems with multiple delay,” Serbian Journal of Electrical Engineering, vol. 7, no. 1, pp. 69–79, 2010.
View at: Publisher Site | Google Scholar
R. Stanisławski and K. J. Latawiec, “Stability analysis for discrete-time fractional-order LTI state-space systems. Part I: new necessary and sufficient conditions for the asymptotic stability,” Bulletin of the Polish Academy of Sciences Technical Sciences, vol. 61, no. 2, pp. 353–361, 2013.
View at: Publisher Site | Google Scholar
R. Stanisławski and K. J. Latawiec, “Stability analysis for discrete-time fractional-order LTI state-space systems. Part II: new stability criterion for FD-based systems,” Bulletin of the Polish Academy of Sciences Technical Sciences, vol. 61, no. 2, pp. 363–370, 2013.
View at: Publisher Site | Google Scholar
M. Busłowicz, “Robust stability of positive discrete-time linear systems of fractional order,” Bulletin of the Polish Academy of Sciences Technical Sciences, vol. 58, no. 4, pp. 567–572, 2010.
View at: Publisher Site | Google Scholar
M. Busłowicz and T. Kaczorek, “Simple conditions for practical stability of positive fractional discrete-time linear systems,” International Journal of Applied Mathematics and Computer Science, vol. 19, no. 2, pp. 263–269, 2009.
View at: Publisher Site | Google Scholar
T. Kaczorek, “New stability tests of positive standard and fractional linear systems,” Circuits and Systems, vol. 2, no. 4, pp. 261–268, 2011.
View at: Publisher Site | Google Scholar
T. Kaczorek, Selected Problems of Fractional Systems Theory, Springer, Berlin, Heidelberg, 2011.
View at: Publisher Site
D. Mozyrska and M. Wyrwas, “The -transform method and delta type fractional difference operators,” Discrete Dynamics in Nature and Society, vol. 2015, Article ID 852734, 12 pages, 2015.
View at: Publisher Site | Google Scholar
P. Ostalczyk, “Equivalent descriptions of a discrete-time fractional-order linear system and its stability domains,” International Journal of Applied Mathematics and Computer Science, vol. 22, no. 3, pp. 533–538, 2012.
View at: Publisher Site | Google Scholar
P. Ostalczyk, Discrete Fractional Calculus: Applications in Control and Image Processing, World Scientific, 2016.
View at: Publisher Site
K. Oprzedkiewicz, R. Stanisławski, E. Gawin, and W. Mitkowski, “A new algorithm for a CFE-approximated solution of a discrete-time noninteger-order state equation,” Bulletin of the Polish Academy of Sciences Technical Sciences, vol. 65, no. 4, pp. 429–437, 2017.
View at: Publisher Site | Google Scholar
M. A. Al-Alaoui, “Novel digital integrator and differentiator,” Electronics Letters, vol. 29, no. 4, pp. 376–378, 1993.
View at: Publisher Site | Google Scholar
M. A. Al-Alaoui, “Al-Alaoui operator and the new transformation polynomials for discretization of analogue systems,” Electrical Engineering, vol. 90, no. 6, pp. 455–467, 2008.
View at: Publisher Site | Google Scholar
R. Lopez, Advanced Engineering Mathematics, Addison Wesley Publishing Company, Boston, MA, USA, Second edition, 2001.
D. Matignon, “Stability results for fractional differential equations with applications to control processing,” in Computational Engineering in Systems and Applications Multiconference, pp. 963–968, Lille, France, 1996.
View at: Google Scholar
D. Matignon, “Stability properties for generalized fractional differential systems,” ESAIM: Proceedings, vol. 5, pp. 145–158, 1998.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2018 Rafał Stanisławski et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

759

Downloads

1203

Citations