Research Article | Open Access
Multiparameter Fractional Difference Linear Control Systems
The Riemann-Liouville-, Caputo-, and Grünwald-Letnikov-type fractional order difference operators are discussed and used to state and solve the controllability and observability problems of linear fractional order discrete-time control systems with multiorder and multistep. It is shown that the obtained results do not depend on the type of fractional operators and steps. The comparison of systems is made under the number of steps needed, firstly to achieve a final point, and secondly to distinguish initial conditions for particular operator.
The main topic of the paper concerns systems with fractional difference operators. Fractional sums and operators are generalizations of th order differences. Since discrete dynamics are often used in modeling real phenomena and recently there is noticed great progress in using fractional dynamics for descriptions of some behaviours, the discussion of properties of possible discrete fractional control systems is strongly needed. Particularly, we consider systems with multiorder and multistep, which could help better in an attempt at simulation of processes and could be easily adapted by CAS programs. Properties of the fractional sums and operators were developed firstly by Miller and Ross  and continued by Atici and Eloe [2, 3] and Abdeljawad et al. [4, 5]. Another concept of the fractional sum/difference was introduced in [6–8]. We use in the paper sums and operators with -differences. Since can represent a sample step, the presence of in operators is interesting from both engineering and numerical points of view. We analyze the controllability and observability problem for multiorder and multistep linear systems while the problems of local controllability of nonlinear control systems defined by Caputo-type, Riemann-Liouville-type, and Grünwald-Letnikov-type difference operators, with step and commensurate order, were studied in . Relations between these three types of operators were studied in .
In this contribution, the Riemann-Liouville-, Caputo-, and Grünwald-Letnikov-type fractional order difference operators are discussed and used to state and solve the controllability and observability problems of linear fractional order discrete-time control systems with multiorder and multistep. Our goal is to state basic properties of fractional difference multiorder and multistep linear control systems and to compare using some particular examples if there are some varieties between operators and numbers of steps that we need to have a system being controllable or observable or achieve some exact target. It is shown that the obtained results do not depend on the type of fractional operators and steps. The comparison of systems is made under the number of steps needed, firstly to achieve a final point and secondly to distinguish initial conditions for particular operator. As our results are formulated for multistep systems, they remain also true for systems with the same step along coordinates of solutions and particularly for the common step . We can also conclude similar corollaries comparing systems according to different or the same orders.
In the paper definitions of operators for general step and orders are presented. In our systems we use different orders and steps, so we define the multiorder as , where and multistep , for each . In our paper we give proofs for solutions and some properties of systems defined by Riemann-Liouville-type fractional difference operator and using the formula of transformations between Riemann-Liouville and Grünwald-Letnikov types of operators presented in , we state results also for systems with Grünwald-Letnikov-type operator.
The paper is organized as follows. In Section 2, all preliminary definitions, facts, and notations are gathered. Section 3 presents initial value problems of fractional order systems of multiorder and multistep difference equations with three types of operators. There are considered systems with operators written in general form. In Section 4, the controllability problem for generalized form of systems is stated and solved, while in Section 5 the observability issue is discussed. At the end of Section 5 there is the example of flat fractional system. The problem of observability of this system is discussed.
Now we are listing the necessary definitions and technical propositions that we use in the sequel of the paper. Let , and put Then . For a function the forward -difference operator is defined as (see ) where . In many papers there are defined fractional summations using the special kind of function called factorial function; see, for example, [7, 8, 10]. In order to have here simpler notation we use in definition generalized binomial function , where and is the Euler gamma function. We use the convention that division at a pole yields zero.
Definition 1. For a function the fractional -sum of order is given by where , , and additionally .
For we write shortly instead of . Note that . Two fractional -sums can be composed as follows.
Proposition 2 (see ). Let be a real valued function defined on , where . For the following equalities hold: where .
In  the authors prove the following lemma that gives transition between fractional summation operators for any and .
Lemma 3. Let and . Then, where and .
In our consideration the crucial role is played by the power rule formula; see : where , , and . Then for we have . Next, taking , , we get
As the tool in the present consideration we can define the next family of functions. Let Moreover, let and . Then we define the family of functions by the following values: Note that the following relations hold:(1), , and , where is the unit vector in with 1 at the th position and is the zero vector in ;(2), where ;(3) and as the division by pole gives zero, the formula works also for ;(4).
The properties of two-indexed functions were given in [12, Proposition 2.5] and generalized for -indexed functions in [12, 13] with multiorder . In the paper therein we use functions with step which gives better tool for multistep systems.
Proposition 4. Let for : and , , . Then for holds
Family of functions is useful for solving systems with Caputo-type operators. We formulate also similar family of functions that are used in solutions of systems with Riemann-Liouville-type operators. We define the family of functions by the following values: Note that the following relations hold: (1);(2), , ;(3)for all : and .
For the family of functions we have similar behaviour for fractional summations as for so we can state the following proposition.
Proposition 5. Let for : and , , . Then for holds
Similarly as in , taking we can write that for and , : where , . Also for and we have
The next step is to present three main fractional difference operators: Caputo-type, Riemann-Liouville-type, and Grünwald-Letnikov-type -difference operators. We present all definitions for any positive and scalar function . In fact each of the operators with step can be translated into the according operator with step . We state results about the inversions of operators in two cases for any positive and . Only the case of the Grünwald-Letnikov-type -difference operators has no exact inverse operator, but we can transform an initial value problem stated for this operator into initial value problem for the Riemann-Liouville-type -difference operators, where we change the domain of the state.
Definition 7 (see ). Let . The Caputo-type -difference operator of order for a function is defined by where .
Note that . Moreover, for the Caputo-type -difference operator takes the form , where .
For the Caputo-type fractional difference operator there exists the inverse operator that is the tool in recurrence and direct solving fractional difference equations.
Proposition 8 (see ). Let , , , and be a real valued function defined on . The following formula holds:
We can also state, similarly as in Lemma 3, the transition formula for the Caputo-type operator between the cases for any and .
Lemma 9. Let and . Then, where and .
From this moment for the case we write and the same notation we use for two other operators and fractional summations. Then the result from Proposition 8 is stated as where , , and .
Definition 10. Let . The Riemann-Liouville-type fractional -difference operator of order for a function is defined by
For we have: .
The next propositions give useful formula for transforming fractional difference equations into fractional summations.
Proposition 11. Let , , , and be a real valued function defined on . The following formula holds:
The next Lemma gives the transition formula for the Riemann-Liouville-type operator between the cases for any and .
Lemma 12. Let and . Then, where and .
For the case we write . Then the result from Proposition 11 is stated as where , and with .
Definition 13. Let . The Grünwald-Letnikov-type -difference operator of order for a function is defined by where
We can easily see the following comparison.
Proposition 14 (see ). Let . Then where or for and .
From Proposition 14 we get in the next chapters the possibility of stating exact formula for solution of initial value problems with the Grünwald-Letnikov-type difference operators by the comparison with parallel problems with the Riemann-Liouville-type operators.
The operators presented in this section can be extended to operators acting on vector valued functions in a componentwise manner.
3. Multiorder and Multistep Fractional Difference Linear Systems
In this section we consider initial value problems of fractional order systems of multiorder and multistep difference equations with three types of operators. In the next section we state as the conclusion the solution of initial value problem for control systems with the Caputo-, the Riemann-Liouville-, and the Grünwald-Letnikov-type difference operators and formulate common results for the controllability as well for observability. We consider systems with operators written in the general form where and , with the initial condition Therein functions , , , if or and if . We stress that for the Grünwald-Letnikov-type the left side of (27) is shifted; compare . By Lemma 9 we can write system (27) in the form where .
Proof. The most important step in the proof is to use the properties of fractional summations with different orders of functions from families and . Let and . For we see that ; hence, the formula agrees with the initial condition. The only one doubt that we should check is the value of for the case of Riemann-Liouville-type operator. In fact we have . By Proposition 8 we obtain the recursive formulas for the Caputo-type operator, and , :
Since for the Riemann-Liouville-type operator we start with , from Proposition 11 we obtain
Using matrices and functions , recursive relations (33) and (34) could be rewritten in the form
or equivalently as
It is important to stress that the recursive formula (36) and the initial value problem given by the system (27) with initial condition (28) are equivalent. Moreover, the recursive formula (36) represents the unique solution to the initial value problem. In order to prove that formula (30) is the solution to the stated problem we can check that it satisfies the equality (36) and initial condition. The zero step for initial condition has been already checked. Then, the crucial role is played by the following property:
that follows from Propositions 4 and 5. Moreover, as for sequences such that we have ; then the values of the matrix in the formula (31) can be written as an infinite summation:
Let . Then the right-hand side of the recursive formula (36) takes the form Hence the thesis holds.
In this section we construct the solution of initial value problems and consider the controllability problem of multistep and multiorder difference linear control systems for different forms of operators.
We consider systems with operators written in the general form where and and with initial condition Therein functions , , with and if or and if . Moreover, are constant-row matrices. By Lemmas 9 and 12 we can write system (40) in the form
Theorem 16. Let be the th row which is equal to . Then, the solution of system (42) with the initial condition which is unique and given by where and
Proof. The steps in the proof are similar to the case with one value for the step like it is presented in .
We can consider also the value . Then for the solution of (40) can be rewritten in the form
Note that controllability means that the final state can be reached in -steps, where .
Let us define the controllability matrix for the system (42) in the following way, as in : where and denotes the transposition of the matrix . It can be noticed that is a symmetric -dimensional matrix.
The next propositions give the formula for the values of steering control and also rank condition for complete controllability. The results have the same meaning like in the classical theory. For the fractional difference systems with the Caputo-type operator with one step and two orders some conditions were proved in . The basic difference here is the form of the matrices and , where we have multiplication of rows by or only for , which can influence the rank of controllability matrix . See the example in the end of the next section. It is important to notice that the form of the controllability matrix and of the formula of values do not depend on the type of the operator.
Proposition 18 (see ). If the matrix is nonsingular, then the following control function transfers the initial state to the final state .
Proposition 19. Control system (40) is controllable in steps if and only if .
Let us consider multiorder and multistep difference linear control system of the form (40), written as a common form for each operator, where the output function is the same as initial condition given by (41), where . Note that for from the output relation (52) it follows that the output trajectory of this system, denoted by , is given by where is the initial state described by (41) and .
Definition 20. System (40) with the output function (52) is observable in steps, if having a control , and output sequence one can determinate the initial condition of this system.
Let us define a real matrix where and denotes the transposition of . Since it is defined in the similar manner as the classical observability Gramian, we call it the (fractional) observability Gramian. Observe that and so .
At the end we give the example in which we show that for a given system with the Riemann-Liouville-type operator we have faster recognition of the initial state than for the system with the Caputo-type operator.
The matrix is called the (fractional) observability matrix. Directly from the fact that matrix is nonsingular and from the proof of Theorem 21 it follows that this matrix is positively defined.
Example 23. Let us consider linear control system with two fractional orders and with two steps : with initial condition . Let , . (1)Firstly we consider the situation with the Caputo-type difference and . Note that and , . Then . This implies that so our system is not observable in step. For we have Now we are ready to state the value of fundamental matrix for : . Then . Hence for holds , so the system is observable in step.(2)Taking into account the Caputo-type operator and , we get . This implies that hence then system is observable in step.(3)Now we take the Riemann-Liouville-type operator; then we need values of matrix function . Observe that then . And for ; hence then system is observable in step.
From our consideration it follows that we could not distinguish which model/type of the system is better according to the controllability or the observability properties. In fact it depends on the matrices inside the systems. Moreover parameters that one can include in the description of some process involve the model and the type operator with ramifications of the used model. We cannot downplay the role of the used orders and steps. However this paper does not discuss any particular process. This is rather theoretical consideration of possible properties.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
The author is grateful to anonymous referees for valuable suggestions and comments, which improved the quality of the paper. The project was supported by the funds of National Science Centre granted on the bases of the decision number DEC-2011/03/B/ST7/03476. The work was supported by Bialystok University of Technology Grant G/WM/3/2012.
- K. S. Miller and B. Ross, “Fractional difference calculus,” in Proceedings of the International Symposium on Univalent Functions, Fractional Calculus and Their Applications, pp. 139–152, Nihon University, Kōriyama, Japan, 1989.
- F. M. Atici and P. W. Eloe, “A transform method in discrete fractional calculus,” International Journal of Difference Equations, vol. 2, no. 2, pp. 165–176, 2007.
- F. M. Atici and P. W. Eloe, “Initial value problems in discrete fractional calculus,” Proceedings of the American Mathematical Society, vol. 137, no. 3, pp. 981–989, 2009.
- T. Abdeljawad and D. Baleanu, “Fractional differences and integration by parts,” Journal of Computational Analysis and Applications, vol. 13, no. 3, pp. 574–582, 2011.
- T. Abdeljawad, F. Jarad, and D. Baleanu, “A semigroup-like property for discrete Mittag-Leffer functions,” Advances in Difference Equations, vol. 2012, article 72, 2012.
- N. R. O. Bastos, R. A. C. Ferreira, and D. F. M. Torres, “Necessary optimality conditions for fractional difference problems of the calculus of variations,” Discrete and Continuous Dynamical Systems, vol. 29, no. 2, pp. 417–437, 2011.
- R. A. C. Ferreira and D. F. M. Torres, “Fractional -difference equations arising from the calculus of variations,” Applicable Analysis and Discrete Mathematics, vol. 5, no. 1, pp. 110–121, 2011.
- N. R. O. Bastos, R. A. C. Ferreira, and D. F. M. Torres, “Discrete-time fractional variational problems,” Signal Processing, vol. 91, no. 3, pp. 513–524, 2011.
- D. Mozyrska and E. Pawluszewicz, “Local controllability of nonlinear discrete-time fractional order systems,” Bulletin of the Polish Academy of Sciences, vol. 61, no. 1, pp. 251–256, 2013.
- D. Mozyrska, E. Girejko, and M. Wyrwas, “Comparison of -difference fractional operators,” in Advances in the Theory and Applications of Non-integer Order Systems, W. Mitkowski, J. Kacprzyk, and J. Baranowski, Eds., vol. 257 of Lecture Notes in Electrical Engineering, pp. 191–197, Springer, 2013.
- D. Mozyrska and E. Girejko, “Overview of fractional -difference operators,” in Operator Theory: Advances and Applications, vol. 229, pp. 253–268, Birkhäuser, 2013.
- D. Mozyrska, “Controllability of -difference linear control systems with two fractional orders,” International Journal of Systems Science, 2013.
- E. Pawluszewicz and D. Mozyrska, “Constrained controllability of -difference linear systems with two fractional orders,” in Advances in the Thoery and Applications of Non-Integer Order Systems, W. Mitkowski, J. Kacprzyk, and J. Baranowski, Eds., vol. 257 of Lecture Notes in Electrical Engineering, pp. 67–76, Springer, 2013.
- D. Mozyrska and E. Pawluszewicz, “Controllability of -difference linear control systems with two fractional orders,” in Proceedings of the 13th International Carpathian Control Conference (ICCC '12), I. Petráš, I. Podlubny, K. Kostúr, J. Kačúr, and A. Mojžišová, Eds., pp. 501–506, May 2012.
- I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, San Diego, Calif, USA, 1999.
- T. Kaczorek, “Fractional positive linear systems,” Kybernetes, vol. 38, no. 7-8, pp. 1059–1078, 2009.
- T. Kaczorek, “Reachability of cone fractional continuous-time linear systems,” International Journal of Applied Mathematics and Computer Science, vol. 19, no. 1, pp. 89–93, 2009.
Copyright © 2014 Dorota Mozyrska. 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.