Abstract

This paper presents new results on modeling and analysis of dynamics of fractional-order discrete-time linear time-invariant single-input single-output (LTI SISO) systems by means of new, two-layer, “fractional-order discrete-time Laguerre filters.” It is interesting that the fractionality of the filters at the upper system dynamics layer is directly projected from the lower Laguerre-based approximation layer for the Grünwald-Letnikov difference. A new stability criterion for discrete-time fractional-order Laguerre-based LTI SISO systems is introduced and supplemented with a stability preservation analysis. Both the stability criterion and the stability preservation analysis bring up rather surprising results, which is illustrated with simulation examples.

1. Introduction

It is well known that orthonormal basis functions (OBF) can be a good alternative to, for example, state space or ARX models when solving the problems in modeling and identification of some classes of “regular” or integer-order dynamical systems, both linear [13] and nonlinear ones [4, 5]. Recently, the OBF machinery has increasingly been used in modeling and identification tasks for noninteger or fractional-order systems [69], a broad research field attracting a huge interest both from the academia and industrial environments [1017]. In this paper we limit ourselves to discrete-time fractional-order systems. There are two main research directions for potential use of OBF, in particular Laguerre functions, in the fractional-order modeling context, namely, Laguerre-based approximation of fractional-order Grünwald-Letnikov (GL) derivatives [8, 9] and Laguerre-based modeling of fractional-order dynamical systems [18, 19]. We pursue both research directions, employing discrete-time Laguerre filters to effectively model both the GL fractional-order difference (FD) [8, 9] and dynamics of a discrete-time fractional-order plant to be covered here. Such a two-layer Laguerre modeling idea is a new contribution to modeling of linear discrete-time fractional-order systems. Our new Laguerre-based approximators to FD, called FLD and FFLD [8, 9], follow the lines of the introduction of the “classical” (truncated or) finite FD, that is, FFD.

The FD and FFD [20] are often used in the context of modeling of discrete-time fractional-order LTI state space systems [14, 2022]. It is the asymptotic stability that is considered the most important aspect of an analysis of fractional-order LTI systems [2328]. As for continuous-time systems, the celebrated Matignon’s stability criterion [29] has established a simple argument condition for eigenvalues of a state matrix. It is only recently that a thorough asymptotic stability analysis has been made for discrete-time fractional-order LTI state space systems [30, 31] and the long-awaited stability criterion has been introduced for such systems [31]. That criterion will be reformulated here for LTI systems described by a new fractional-order model based on “fractionalized” Laguerre filters.

In this paper, a new concept of “fractional-order discrete-time Laguerre filters” is introduced for modeling of fractional-order LTI SISO systems and a new analytical stability criterion is announced. Also, the stability preservation problem is solved.

The remainder of this paper is organized as follows. Fundamentals of the regular OBF system description, in particular the Laguerre filters, are recalled in Section 2, with an important factorization of the Laguerre transfer functions presented. The first significant result of this paper is provided in a unified framework in Section 3, where a new term of “fractional-order discrete-time Laguerre filters” is coined, supported with a series of preliminary outcomes preparing for the main result of Section 4. That section provides a new general criterion for the asymptotic stability of discrete-time LTI SISO fractional-order Laguerre systems, being a nice supplementation of the criterion of [31]. The criterion is followed by a stability preservation analysis in Section 5. The parameter estimation problem for fractional-order Laguerre systems is outlined in Section 6 and the reorthonormalization issue is addressed in Section 7. Conclusions of Section 8 summarize the achievements of the paper.

2. Regular OBF System Description

An output of a “regular” (or integer-order) discrete-time OBF-modeled dynamical system or shortly OBF system can be described aswhere is the backward shift operator, and are the system input and output, respectively, at discrete-time , and and , are orthonormal transfer functions and weighting parameters, respectively. For discrete-time Laguerre filters there iswith and being a dominant pole. We proceed with the practically justified case of . The unknown parameters , can be easily estimated via, for example, Recursive Least Squares (RLS) or Least Mean Squares (LMS) algorithms using the linear regression formalism. We discriminate between the parameters of the Laguerre model of a dynamical system and those for the Laguerre-based fractional-order difference as in [8, 9].

Remark 1. Note that we do not need to account for the sampling period in as it will be indirectly included in the estimated parameters .

Remark 2. Note that pursuing an optimal Laguerre pole has been well established [2, 3234].

Depending on the time-domain or -domain contexts, we will interchangeably use the notations or , respectively. The same manner will be applied to other functions of or .

A crucial factorization is recalled; that is, (2) is rewritten as [19, 35]withand the filter outputs obtained as and , withThe filters and of (4) can be described in terms of the backward differenceswhere and .

3. FD/FFD/LD/FLD/CFLD/FFLD-Based Laguerre Systems: Unified Framework

We refer to the Grünwald-Letnikov (GL) discrete-time fractional-order difference (FD) and its Laguerre-based equivalents LD and CFLD [8, 9] as well as to their approximators FFD [20], FLD, and FFLD [8, 9]. We strongly advocate the use of the FFLD approximator to FD, being an effective combination of FFD and FLD [9]. In FFLD, the FFD share is to precisely model the high-frequency properties of FD whereas FD’s long “tale” is excellently modeled by FLD [9]. We will shortly write “FD/FFD/LD/FLD/CFLD/FFLD-based Laguerre system” instead of “Laguerre model of a fractional-order system based on FD/FFD/LD/FLD/CFLD/FFLD.”

Here is the first new result of this paper, that is, the introduction of a concept of “fractional-order discrete-time Laguerre filters.”

Theorem 3. Consider the FD/FFD/LD/FLD/CFLD/FFLD fractional-order differences [8, 9] uniformly described aswith , and is the upper bound for when [9, 20], and , , being the Laguerre filters and the particular differences specified as (i)FD for and ,(ii)FFD for and ,(iii)LD for , , and calculated as in Theorem  1 of [8, 9],(iv)FLD for , , and calculated as in Theorem  1 of [8, 9],(v)CFLD for , , and calculated as in Theorem  2 of [8, 9],(vi)FFLD for and calculated as in Theorem  2 of [8, 9].Then FD/FFD/LD/FLD/CFLD/FFLD-based fractional-order LTI SISO Laguerre system can be described aswhere and are the system input and output, respectively, , , are the unknown parameters, and the “fractional-order discrete-time Laguerre filters” are as follows:with being the dominant Laguerre pole andwith andpertaining to specifications (i) to (vi) for the particular differences.

Proof. Accounting for (6) on the one hand and (8) on the other hand (compare [9]), the fraction-formalized filters and can now be described aswithand , is computed asFinally, the outputs from the two fractional filters arewith , calculated as in (17). Now, from (18) and (19) we can easily obtain transfer functions (11) and (12), respectively.
Conclusively, the output from the FD/FFD/LD/FLD/CFLD/FFLD-based fractional Laguerre system iswhere , and are given as in (10).

An outstanding value of Theorem 3 is at least threefold. Firstly, a new concept of “fractional-order discrete-time Laguerre filters” is introduced. Secondly, the factorization of the fractional-order Laguerre filters holds true in the same way as for the regular (nonfractional) ones, with (10) being still valid. This means that the block diagram for the fractional Laguerre system is identical with that for the regular one, with and substituted for and , respectively. This provides a nice technical tool for modeling in, for example, Matlab. Thirdly, the whole calculation process for the model output of the fractional Laguerre system is quite similar to that for the regular Laguerre system, with an additional fractional component involved in the former case (see (11)).

Remark 4. Note that FFLD as in (7) and (8) can be considered as one general unified model from which all the other considered models can be obtained according to specifications (i) through (v), with FFLD itself covering specification (vi).

Remark 5. Note that the whole “fractionality” of the otherwise (two-layer) fractional-order Laguerre system is contained in the (lower-layer) transfer function . Also note that for the FD-based fractional Laguerre system we have the familiar GL transfer function (compare [30]). Of course, for we arrive at the classical Laguerre filters and .

Remark 6. Let us emphasize that the transfer function is related to the function (or ) introduced in [30] in the asymptotic stability analysis of fractional-order discrete-time state space systems. In fact, we havewith being the imaginary unit and .

Remark 7. Possible accounting for the sampling period in FD/FFD/LD/FLD/CFLD/FFLD-based Laguerre systems when transferring from a continuous-time (Riemann-Liouville/Grünwald-Letnikov) derivative to the discrete-time difference results in the substitutions and . This finally leads to the substitution in (11).

Remark 8. In addition to computational simplicity, the Laguerre-based difference is attractive also in that it can model both well-damped () and oscillatory behaviors (), without referring to, for example, the Kautz filters in the latter case.

An exemplary block diagram of the output calculation process for the FLD-based fractional Laguerre system is presented in Figure 1, with the “internal” (or lower-layer) FLD being additionally zoomed.

4. Stability Results for Fractional-Order Laguerre Systems

4.1. Stability of FD/FFD/LD/FLD/CFLD/FFLD-Based Laguerre Systems: A Unified Framework

We adopt the stability results of [30, 31] to the Laguerre-based fractional systems to obtain an original outcome which is the main result of this paper.

Theorem 9. Consider the FD/FFD/LD/FLD/CFLD/FFLD-based Laguerre system described by (7) to (20) with and , computed from (15) and (16) under respective specifications (i) to (vi) using and , as in (18) and (19), respectively. The system is asymptotically stable if and only ifwhere and defined as in (21), with (i) for the FD/LD/CFLD-based systems,(ii) for the FFD-based system,(iii), with , as in Theorem  1 of [8, 9] for the FLD-based system,(iv) with , , as in Theorem  2 of [8, 9] for the FFLD-based system.

Proof. Recall the specific form of the fractional Laguerre filters (10) incorporating the cascade combination of the first-order filters and , whose poles are just , which is the (only) eigenvalue in the corresponding (scalar) state space model as in [30, 31]. Since the argument of that real-valued eigenvalue is zero, we consider the stability contour involving the function as in Theorem  7 of [30], but for equal to zero. Now, for the FD-based system we have under specification (i) of Theorem 3, so that and the stability condition is just . Since LD/CFLD are equivalent to FD [8, 9], we arrive at the same condition for the LD/CFLD-based systems. For the FFD-based system we have under specification (ii) of Theorem 3, so , thus ending up with condition (22). In a similar way, for the FLD and FFLD-based systems we have, respectively, , with , , as in Theorem  1 of [8, 9], and with , , as in Theorem  2 of [8, 9], leading again to condition (22). Finally, the cascade combination of the considered filters is asymptotically stable if and only if all the particular filters are asymptotically stable, which completes the proof.

It is rather surprising that the Laguerre-based FFD system can be asymptotically stable even for (and ). This stability result yielded from (22) is a nice theoretical confirmation of our earlier surprising simulation outcome of [35]. In fact, with under and finite , the pole can be admitted higher than unity. But on the other hand, for the term can be lower than so that even for the FFD system may be unstable, this being yet another nice stability result confirming our earlier surprising simulation observations.

Out of a plethora of simulations runs, we firstly present two selected examples.

Example 10. Consider the FFD-based fractional Laguerre system with and two different values of and . It can be concluded from Theorem 9 that for the system is asymptotically stable for . On the other hand, for the system is unstable for .

In a similar way, the FLD/FFLD-based fractional Laguerre systems can be stable even for , which is illustrated in the following.

Example 11. Consider FLD/FFLD-based Laguerre systems with . Let and in the FLD-based system and , , and in the FFLD-based one. For the FLD-based Laguerre system we have and for the FFLD-based one we have , so in both cases condition (22) is satisfied and the two systems are asymptotically stable.

Remark 12. Note that the “closedness” of to zero is an excellent indicator of the quality of the specific approximator FFD/FLD/FFLD to FD/LD/CFLD.

Remark 13. In order to account for the sampling period when transferring from the continuous-time derivative to FD/LD/CFLD/FFD/FLD/FFLD, in the stability analysis for the fractional Laguerre system we replace (22) withwith as in Theorem 9.

Remark 14. It is interesting that the stability condition for fractional-order Laguerre systems based on the FD difference (and equivalently, LD and CFLD ones [8, 9]), that is, , is identical with that for the regular integer-order Laguerre system. Also note that in this case () the stability condition is, rather surprisingly, independent of both the sampling period and the order . In contrast, criterion (23) for the FFD/FLD/FFLD-based fractional Laguerre systems is always dependent on both and , and this is because in that case.

5. Stability Preservation Analysis

The stability preservation issue for fractional-order systems under direct discretization methods has been considered in [36]. The problem reduces to the provision of a stable discretized system approximating its stable continuous-time original. Here we present our stability preservation results for the FFD/FLD/FFLD-based Laguerre systems.

Theorem 15. Consider the (discretized) FFD/FLD/FFLD-based Laguerre system as in Theorem 3, under respective specifications (ii), (iii), and (iv). Assume that (1) and ,(2) and ,(3) and , or(4) and .Then, respectively, (1)the system is asymptotically stable regardless of the sampling period,(2)the system is asymptotically stable if and only if ,(3)the system is asymptotically stable if and only if ,(4)the system is unstable regardless of the sampling period.

Proof. Immediate from condition (23).

Some remarks are due, now. Firstly, the results for variants 1 and 3 of Theorem 15 are quite surprising. Secondly, we have verified all the variants of Theorem 15 in a plethora of simulation runs, the selected examples of which are tabulated. Table 1 collects the results, with the particular variants of Theorem 15 superindexed in the Comment column. The results of Table 1 are self-explanatory.

6. Parameter Estimation for Fractional Laguerre Filter

When output (9) (or (20)) of the fractional Laguerre filter is corrupted with noise and the order is known, the model output can be described in a linear regression fashion , with and the unknown parameters estimated analytically, for example, (R)LS method. When is unknown the parameter vector can be estimated numerically using, for example, LS or, preferably, GA procedures.

Example 16. An FD-based discrete-time state space system , with its fractional (commensurate) order and and the output white noise of zero mean and standard deviation , is modeled by means of the FFD- (), FLD- (), and FFLD- (, ) based fractional Laguerre models with optimally selected (via a combined LS/GA algorithm) and . The parameter is selected (heuristically) to be equal to and , parameters are estimated by the RLS method. The input is a discrete-time zero mean unity-variance Gaussian random signal. The MSPEs for the outputs of the three Laguerre-based models with respect to that for the FD-based state space system are presented in Table 2.
Table 2 shows that the best performance is obtained for the FFLD-based Laguerre model. In this case, the FFLD-based fractional Laguerre model with 45 parameters gives essentially lower MSPE than the FFD-based one with 500 parameters. This shows that the FFLD is a very good approximation to FD. Also, a simple FLD-based fractional Laguerre model, with 35 parameters only, gives satisfactory results.
The frequency responses of the actual system and (hardly distinguishable) FFD/FLD/FFLD-based fractional Laguerre models are presented in Figure 2. Table 3 presents standard deviations of the magnitude and phase of the frequency responses for the FFD/FLD/FFLD-based fractional Laguerre models with respect to the FD-based system.

Remark 17. It is interesting that fractional Laguerre-based models may have their fractional orders () different from those for fractional state space originals (); see Example 16. The discrepancy is related to the spectrum of eigenvalues of . Specifically, for the “slightly” oscillatory eigenvalues as in Example 16 () we have closer to unity than the original . Also, “strong” oscillatory fractional behaviors may be Laguerre modeled with ; compare Remark 8.

7. The Reorthonormalization Issue

A side effect of factorizations (3) or (10) is that the orthonormality property is lost for fractional Laguerre systems; that is, the functions , are no longer orthonormal. However, the fractional Laguerre filters can be easily reorthonormalized using the standard procedures [18, 37], with the results illustrated with Example 18.

Example 18. Consider the system as in Example 16 and two FFLD-based fractional Laguerre models with parameters as in Example 16 with (a) nonorthonormalized and (b) reorthonormalized Laguerre filters. The model parameters are estimated using the RLS method.
The simulation experiments show that both models give quite the same MSPE errors equal to some . So, the reorthonormalization process may not improve the model accuracy. However, it is well known that the reorthonormalization process improves conditioning of the covariance matrix in the RLS estimation process. This is shown in Figure 3, which presents the condition numbers of the covariance matrix as a function of time for both orthonormalized and nonorthonormalized fractional Laguerre filters.

8. Conclusion

This paper has presented a series of new results in modeling and analysis of discrete-time fractional-order LTI SISO systems using Laguerre filters. In the “lower” modeling layer described in [8, 9], the Grünwald-Letnikov fractional-order difference has been effectively approximated by the Laguerre filters. This paper has produced new results related to modeling and analysis of the “upper” system dynamics layer. Firstly, new “fractional-order discrete-time Laguerre filters” have been shown to efficiently model discrete-time fractional-order LTI SISO systems. Secondly, a new analytical stability criterion has been announced for discrete-time fractional-order Laguerre systems. Thirdly, new results on stability preservation for discretized fractional-order Laguerre systems have been derived. New achievements of the paper have been illustrated with simulation examples. Our current and future research works are directed to predictive/adaptive control of Laguerre-based fractional-order discrete-time MIMO systems, both linear and nonlinear ones.

Nomenclature

CFLD:Combined fractional/Laguerre-based difference
FD:Grünwald-Letnikov fractional-order difference
FFD:Finite fractional difference
FFLD:Finite (combined) fractional/Laguerre-based difference
FLD:Finite Laguerre-based difference
GL:Grünwald-Letnikov
LD:Laguerre-based difference
OBF:Orthonormal basis functions.

Competing Interests

The authors declare that they have no competing interests.