Remarks on Fractional Discrete Cone Control Systems with n-Orders and Their Stability

In the paper, fractional discrete cone control systems with n-orders are considered. Some relations between invariance and (asymptotic) stability properties of the presented systems are discussed. Operators employed to the considered systems are Caputo-, Riemman-Louville-, and Grünwald-Letnikov type ones. Cone systems with control, which are particular invariant systems with control, together with their stability and asymptotic stability properties are examined.


Introduction
Positive systems are the systems whose state and input variables are never negative, with a given nonnegative initial state. These systems appear frequently in practical applications as well as in real phenomena, among other in biology, medicine, economics, electrotechnics, control system design, etc. (see [6,7,13,25,26] and the references therein). The natural generalization of positive systems are cone systems, i.e., systems whose trajectories always remain in the given cone if they are initialized in this cone. Moreover, the special attention can be put on invariant cone systems, since they can be employed to design stabilizing controllers.
In recent years, there has been observed a growing interest in the theory and applications of fractional differential and difference equations. Many authors proved that such types of equations are more adequate for modeling physical and chemical processes than equations with integer order. Fractional differential and difference equations describe many phenomena arising in engineering, physics, or economics. In fact, one can find several applications in viscoelasticity, electrochemistry, electromagnetic, etc. For example, Machado in [19] gave a novel method for the design of fractional order digital controllers.
In both theory and applications, one can meet several definitions of the fractional derivatives, among which the most popular are Caputo-, Riemman-Louville-, and Grünwald-Letnikov operators, so there appears the problem how to deal with differences resulting from the application of these operators. The first steps in this topic were made in [20]. Properties of the fractional sum, Caputo-and Riemman-Louville-type difference operators, were developed in [1][2][3][4]21]. Basic information on fractional calculus concept, ideas, and applications of these operators can be found for example in [15,18,24]. In [8], there was adopted a more general fractional h-difference Riemman-Liouville operator, where on the one hand h represents sample step, on the other hand, for h tending to zero, the solutions of the fractional difference equation may be seen as approximations of the solutions of corresponding Riemann-Liouville fractional equations.
The goal of the paper is to examine under which conditions control systems are cone systems. To this aim, in Section 2, we present the needed notation and properties of the h-difference operators of fractional order (with arbitrary h > 0). Operators, which we consider are the three most important among all fractional operators: Caputo-, Riemman-Louville-, and Grünwald-Letnikov type difference operators. In [11], it is shown that these three types of h-difference fractional operators are related to each other. Moreover, the Grünwald-Letnikov-type fractional h-difference operator can be expressed by the Riemann-Liouville-type fractional h-difference operator. So, systems with these operators can be studied simultaneously. Taking into account this fact, in Section 3, there is introduced a discrete-time cone control system with fractional order and properties of its trajectories are discussed. Since a cone is a special case of a polyhedron, in Sections 4 and 5, basing on some properties of polyhedron contractiveness, the problems of stability and asymptotic controllability of class of consider systems are tackled.

Some Preliminaries
For α > 0, h > 0 and a ∈ R let For a function x : (hN) a → R, then the forward h-difference operator is denoted by while the h-difference sum is given by where t = a + (k + 1)h and k ∈ N 0 . Additionally, we define a where t h ∈ Z − := {−1, −2, −3, . . .}, and we use the convention that division at a pole yields zero. Notice that if we use the general binomial coefficient a b := (a+1) (b+1) (a−b+1) , then function (1) can be rewritten as Lemma 1 [10] If 0 < α ≤ 1, then For a function x ∈ (hN) a the fractional h-sum of order α > 0 is defined by h for any r ∈ (hN) a , then (see [8]) where t ∈ (hN) a+αh . Using the general binomial coefficient, one can rewrite Eq. (2) in the form Note that if ψ ≡ 1, then for μ = 0, a = (1 − α)h, and t = nh + a + αh, it holds Definition 2 [5] Let α ∈ (0, 1] and a ∈ R. The Riemann-Liouville-type fractional hdifference operator RL a α h of order α for a function x : (hN) a → R is defined by Definition 3 [21] Let α ∈ (0, 1] and a ∈ R. The Caputo-type fractional h-difference operator C a α h of order α for a function x : (hN) a → R is defined by where t ∈ (hN) a+(1−α)h .
Note that the operator C a α h for any α ∈ (0, 1] changes the domain of the function x, i.e., it maps real valued functions defined on the set (hN) a into real valued functions defined on the set (hN) a+ ( holds also for the Riemann-Liouville-type h-difference operator. Moreover, for α ∈ (0, 1], it holds (see [11]): where z ∈ C is a complex variable for which the series ∞ k=0 y(k)z −k converges absolutely.
The last operator that we take under our consideration is the fractional h-difference Grünwald-Letnikov-type operator.

Cone Systems
In order to define cone systems, let us discuss first the problem of existing of solutions of a nonlinear control system of fractional order. The reasoning is similar to the one given for finding solutions of the nonlinear autonomous system of the fractional order α i ∈ (0, 1] given in [23], so we present only main steps of it.
Let i = 1, . . . , n and 0 < α i ≤ 1. Let us consider the following fractional Caputo h-difference system with n orders α 1 , . . . , α n as follows with initial values where The set U is called the control space and satisfies the following property: U is such that U ⊆ int U and any two points in the same connected component of U can be jointed by a smooth curve lying in int U , except for end points. Let J 0 (m) denotes the set of all sequences u = (u 0 , u 1 , . . .), where u n := u(nh) ∈ U , n ∈ N 0 . We assume that function f depends on finite number of elements u i .  x 1 (a 1 + t), x 2 (a 2 + t), . . . , x n (a n + t), u(t)) . For simplicity, we state all definitions and theorems for the case when the equilibrium point is (0, 0) ∈ R n × U , i.e. x e i = 0, i = 1, . . . , n and u e = 0. There is no loss of generality in doing so because any equilibrium point can be shifted to the origin via certain change of variables.
Then we can write systems (4) or (6) respectively, in the following forms where F : N 0 × R n × J 0 (m) → R n . Therefore, systems (7a) and (7b) can be expressed in one compact form with the initial condition x n (a n + n 0 h) Then the solution of IVP given by Eqs. 8 and 9 can be obtain using the same reasoning as in [23] and in [10]. It is a uniquely defined map γ : R n × (hN) a × J 0 (m) → R n by initial state X 0 and control sequence u ∈ J 0 (m) and described by γ (X 0 , a, u(0)) = X 0 γ (X 0 , a + kh, u(kh)) = I k · X + So, γ (X 0 , ·, u) is defined by its values γ (X 0 , kh, u(kh)) = X(k), k ∈ N 0 , and denotes the state forward trajectory of system (8).
Let us take a nonsingular matrix with ith, i = 1, . . . , n, row given as p i = (p i1 , . . . , p in ). Then the set K P := x ∈ R n : ∀i = 1, . . . , n : is called a linear cone of state generated by the matrix P in R n (see [12,14]). If X n := {X : N 0 → R n }, then the set is called a linear cone of states with the vertex at 0 generated by the matrix P in the space X n where x i : N a i → R.

Definition 9
Let P ∈ R n×n be a given matrix. The nonlinear fractional difference system (8) together with initial condition (9) is called a P cone fractional system if X(·) ∈ P for any X 0 ∈ K P .

Theorem 10
Let K and P be given as in Eqs. 12 and 13 together with X 0 ∈ K P . If there exists a control u ∈ J 0 (m) such that for every x ∈ K P it holds then for every k ∈ N 0 system (4) or (6) is a P cone system.
Proof The proof uses the analogous reasoning as the one given in [10] for the similar result but without a control. For the proof, we use mathematical induction. For k = 1, the formula holds, since if X 0 ∈ K P , then we get X(1) = F (0, X 0 , u(0)) + 1,1 X 0 ∈ K P . Now, we assume that the hypothesis is true for some k, i.e., F (k, x, u) + k,k · x ∈ K P for every x ∈ K P and u ∈ J 0 (m). This means that X(k) ∈ K P . Then, by assumption and Lemma 1, it holds Hence from mathematical induction, the thesis holds for any natural k.
Note that for system (8) together with initial condition (9) and with the right-hand side autonomous also the implication "only if" in Theorem 10 is true.

Corollary 11
Let P be a nonsingular n × n matrix. Then system (8) together with initial condition (9) with the right-hand side F (k, X(k), u(k)) = AX(k) + Bu(k), where A ∈ R n×n and B ∈ R n×m , is a (P , Q) cone system 1 if and only if P · A + k,k P −1 ∈ R n×n + and P BQ −1 ∈ R n×m + for nonsingular matrix Q of the form with ith, i = 1, . . . , m, rows given as q i = (q i1 , . . . , q im ).
Let us define a feedback control law by where κ(x) is a vector function with values in R m . We say that a function κ : R n → J 0 (m) is an admissible feedback law for system (8) if for every k ∈ N 0 there exists a map γ (X 0 , ·, u) given by Eq. 11 such that u(k) = κ(γ (X 0 , a + kh, u(kh))). If κ is admissible feedback law for Eq. 8 (or respectively for Eqs. 7a or 7b)), then the closed loop system is of the form n 0 (α) X (k) =F (k, X(k), κ(X(k))) , where X(k) is given by Eq. 10. As an immediate consequence of Theorem 10, we have the following.

Corollary 12
Let K P and P be given as in Eqs. 12 and 13 together with X 0 ∈ K P . If there exists a feedback law (14) such that for every x ∈ K P , it holds then for every k ∈ N 0 system (15) is a P cone system.
In particular, if system (8) is a linear one, i.e., F (k, X(k), u(k)) = AX(k) + Bu(k), application of linear feedback law κ = F X(k) with a matrix F ∈ R m×n gives the closed loop system of the form
Recall that a set S is an invariant set for the system (15) if and only if every trajectory of this system starting within S remains inside it. Then λ-contractiveness with respect to system (15) implies invariance property for this system.
The one step admissible set to (H, w) with respect to system (15)  for any ε ∈ (0, 1] and natural q. Proof The result is an immediate consequence of Definition 13 and of polyhedron (18).
Assumption 1 Let us assume that functionF given in Eq. 15 is continuous in all variables and (classically) continuously differentiable at the equilibrium point of the given system. Under Assumption 1, let us define matrix and consider a linear fractional order discrete-time system System (20) is called a linear approximation of the nonlinear one given by Eq. 15. Note that for a given initial condition and for an arbitrary sequence of controls u ∈ J 0 (m), there exists the unique solution of linear approximation (20). Recall that a constant vector X e = (0, . . . , 0) is an equilibrium point of fractional difference system (20) if and only if ( n 0 (α) X e )(k) = AX e for all k ∈ N 0 . Let us notice that the trivial solution X ≡ 0 is an equilibrium point of system (20). The equilibrium point X e = 0 of Eq. 20 is said to be (a) stable if, for each > 0, there exists δ = δ ( ) > 0 such that X 0 < δ implies X(k) < , for all k ∈ N 0 ; (b) attractive if there exists δ > 0 such that X 0 < δ implies lim k→∞ X(k) = 0 ; (iv) asymptotically stable if it is stable and attractive.
System (20) is called stable/asymptotically stable if its equilibrium point X e = 0 is stable/asymptotically stable.
In the linear case, the effective characterization of asymptotic stability was given in [22].

Remarks on Asymptotic Controllability of Cone Systems
Let us draw our attention to the problem of controllability of system (8). Classically, controllability of the given system means that it is possible to transfer the considered system from a given initial state to a final state using controls from a certain set, see for example [16,17]. Suppose that a set V is a subset of the state space of system (8).
i. X is asymptotically controlled to a final state X f without leaving V if there exists a control u ∈ J 0 (m) such that lim k→∞ γ (X 0 , a + kh, u(kh)) = X f and γ (X 0 , a + kh, u(kh)) ∈ V for all k ∈ N 0 . ii. If X e is an equilibrium, then system (8) is asymptotically controlled to X e if for each neighborhood V of X e there is some neighborhood W of X e such that each X ∈ W can be asymptotically controlled to X e without leaving V .

Proposition 20
Let K P and P be given as in Eqs. 12 and 13 together with X 0 ∈ K P . Suppose that there exists some feedback law u(k) = κ(X(k)) so that X e is a local asymptotically stable state for P cone system (15). Then system (8) is asymptotically controlled to X e .
Proof Since asymptotical stability of P cone system (15) at X e means that lim k→∞ γ (X 0 , a + kh, κ(k)) = X e and γ (X 0 , a + kh, κ(X(k)) is in a neighborhood of X e for all k ∈ N 0 , hence the thesis follows from Definition 19.
Let us assume function F given in Eq. and consider a linear fractional order discrete-time system System (22) is called a linear approximation of the nonlinear one given by Eq. 8.
Proof If there is law u(k) = κ(X(k)) such that (H, ωελ) is λ-contractive set with respect to system (20) and (H, ωελ) is a compact polyhedron, then by Proposition 17 the closed loop system (15) is asymptotically stable. Hence, by Definition 19, system (22) is asymptotically controllable at X e = X 0 .

Conclusions
We examined fractional discrete cone control systems with n-orders. Some relations between invariance and (asymptotic) stability properties of the presented systems where discussed. Since there are several definitions and notations of the fractional derivatives, among which the most popular are Caputo-, Riemman-Louville-, and Grünwald-Letnikov operators, we employ right them as fractional discrete Caputo-, Riemman-Louville-, and Grünwald-Letnikov type operators to the systems. In the paper, there were considered cone systems with control, which are particular invariant systems with control, together with their stability and asymptotic stability properties.