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 systemswithmultiorder 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.


Introduction
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 [1] 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][7][8].We use in the paper sums and operators with ℎ-differences.Since ℎ > 0 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 Caputotype, Riemann-Liouville-type, and Grünwald-Letnikov-type difference operators, with step ℎ = 1 and commensurate order, were studied in [9].Relations between these three types of operators were studied in [10].
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 ℎ = 1.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 () = ( 1 , . . .,   ), where   ∈ (0, 1] and multistep (ℎ) = (ℎ 1 , . . ., ℎ  ), ℎ  > 0 for each  = 1, . . ., .In our paper we give proofs for solutions and some properties of systems defined by Riemann-Liouvilletype fractional difference operator and using the formula of transformations between Riemann-Liouville and Grünwald-Letnikov types of operators presented in [10], 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.
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  () : Z → R by the following values: Note that the following relations hold: (1) (3) for all () ∈ R  :  () () ≥ 0 and  * () () ≥ 0. For the family of functions  () we have similar behaviour for fractional summations as for  * () so we can state the following proposition.
Similarly as in [14], taking ℎ = 1 we can write that for   ∈ N 0 and   > 0,  = 1, . . ., : where 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 ℎ > 0 can be translated into the according operator with step ℎ = 1.We state results about the inversions of operators in two cases for any positive ℎ and ℎ = 1.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 [11]).Let  ∈ (0, 1].The Caputo-type ℎ-difference operator  Δ  ℎ, * of order  for a function  : (ℎN)  → R is defined by where  ∈ (ℎN) +(1−)ℎ .
Note that  Δ  ℎ, *  : (ℎN) +(1−)ℎ → R.Moreover, for  = 1 the Caputo-type ℎ-difference operator takes the form , where  ∈ (ℎN)  .For the Caputo-type fractional difference operator there exists the inverse operator that is the tool in recurrence and direct solving fractional difference equations.
The next presented operator is called fractional ℎdifference Riemann-Liouville-type operator.The definition of the operator can be found, for example, in [2] (for ℎ = 1) or in [6] (for any ℎ > 0).Definition 10.Let  ∈ (0, 1].The Riemann-Liouville-type fractional ℎ-difference operator  Δ  ℎ  of order  for a function  : (ℎN)  → R is defined by For  = 1 we have: The next propositions give useful formula for transforming fractional difference equations into fractional summations.
The third type of the operators that we take into our consideration is the fractional ℎ-difference Grünwald-Letnikovtype operator; see, for example, [15][16][17] for case ℎ = 1 and also for case ℎ > 0.
Definition 13.Let  ∈ R. The Grünwald-Letnikov-type ℎdifference operator  Δ ℎ of order  for a function  : (ℎN)  → R is defined by where We can easily see the following comparison.
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.

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  = 1, . . .,  and  ∈ N 1 , with the initial condition Therein functions   : (ℎN (ℎ  + ℎ  ).We stress that for the Grünwald-Letnikov-type the left side of (27) is shifted; compare [10].By Lemma 9 we can write system (27) in the form where ]  =   /ℎ  =   − 1.

Theorem 15. The solution of system (27) with initial condition (28) is of the form
. . .

𝑥 𝑎 𝑛
] and For  = 0 we see that Φ(0) =  (0)  (0) (0) =   ; hence, the formula agrees with the initial condition.The only one doubt that we should check is the value of  (0) (0) for the case of Riemann-Liouville-type operator.In fact we have  (0) (0) = diag (  1 (0), . . .,    (0)) =   .By Proposition 8 we obtain the recursive formulas for the Caputo-type operator,  = 1, . . .,  and  ∈ N 1 , ]  =   − 1: Since for the Riemann-Liouville-type operator we start with ), from Proposition 11 we obtain Using matrices  (0) ,    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 ∑  =1   >  we have  () () = 0; then the values of the matrix Φ in the formula (31) can be written as an infinite summation: Let  = ∑  =1   .Then the right-hand side of the recursive formula (36) takes the form Hence the thesis holds.

Controllability
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  = 1, . . .,  and  ∈ N 1 and with initial condition Therein functions   : (ℎN Moreover,   ∈ R 1× are constant-row matrices.By Lemmas 9 and 12 we can write system (40) in the form (45) Proof.The steps in the proof are similar to the case with one value for the step ℎ = 1 like it is presented in [14].
Note that controllability means that the final state can be reached in -steps, where  ∈ N 1 .
Let us define the controllability matrix for the system (42) in the following way, as in [13]: 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 () = ( 1 ,  2 ) some conditions were proved in [13].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 () > 0 do not depend on the type of the operator.
Proposition 18 (see [13]).If the matrix   is nonsingular, then the following control function

Observability
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  ∈ N 0 .Note that for  ∈ N 0 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 (−1) = 0.
where  ∈ N 0 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  ∈ R × and so Φ() ∈ R × .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. . . .

𝐶Φ (𝑞)
] ] ] is called the (fractional) observability matrix.Directly from the fact that matrix O  is nonsingular and from the proof of Theorem 21 it follows that this matrix is positively defined.

Conclusions
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.

Definition 17 .
System (40) is (completely) controllable in steps from the initial state (41) to the final state

Theorem 22 .
System (40) with the output (52) is observable in  steps if and only if one of the following conditions is satisfied:(i) columns of the matrix Φ() are linearly independent;