Abstract and Applied Analysis

Volume 2013, Article ID 475350, 11 pages

http://dx.doi.org/10.1155/2013/475350

## On Solutions to Fractional Discrete Systems with Sequential -Differences

Faculty of Computer Science, Bialystok University of Technology, Wiejska 45A, 15-351 Białystok, Poland

Received 28 August 2013; Accepted 19 October 2013

Academic Editor: Ali H. Bhrawy

Copyright © 2013 Małgorzata Wyrwas 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.

#### Abstract

We study the subject of a behaviour of the solutions of systems with sequential fractional -differences. We give formulas for the unique solutions to initial value problems for systems in linear and semilinear cases. Moreover, the sufficient condition that guaranties the positivity of considered systems is presented.

#### 1. Introduction

The first definition of the fractional derivative was introduced by Liouville and Riemann at the end of the 19th century. Later on, in the late 1960s, this idea was used by engineers for modeling various processes. Thus the fractional calculus started to be exploited since that time. This calculus is a field of mathematics that grows out of the traditional definitions of integrals, derivatives, and difference operators and deals with fractional integrals, derivatives, and differences of any order. Many authors prove that fractional differential and difference equations are more adequate for modeling physical and chemical processes than integer-order equations. Fractional differential and difference equations describe many phenomena arising in engineering, physics, economics, and science. In fact, several applications can be found in viscoelasticity, electrochemistry, electromagnetic, and so forth. For example, Machado [1] gave a novel method for the design of fractional order digital controllers. Fractional difference calculus has been investigated by many authors, for example, [2–12] and others. In particular, different delta and nabla type fractional differences have been studied in [13–15], where the authors relate these differences by deriving some dual identities. The calculus of fractional -differences was given, for instance, in [5, 10, 16–19]. The properties of systems defined by the fractional difference equations where studied, for example, in [20–24]. There exist definitions of sequential operators in continuous case with different types of derivatives like Caputo type, Riemann-Liouville type, and Hadamard type, see, for example, [9, 12, 25–27]. In parallel with this paper we developed the theory of fractional -differences with sequential operators in the paper [28], where the approximation of continuous fractional sequential derivative is considered. We compute different formulas of solutions and then we try to check the system's behaviour, precisely the positivity of solutions. As far as we know the subject of positivity is well developed for fractional linear systems with continuous time; see [8, 29–32]. However, positivity of fractional discrete systems with sequential -differences is still a field to be examined. In the present paper we open our studies in this field. We give formulas for the unique solutions to initial problems for systems in linear and semilinear cases. Moreover, the positivity of systems with sequential -differences is considered. We consider systems with sequential -differences of Caputo type, while in [32] systems with Grünwald-Letnikov operator are studied. It is worth to add that in [32] (and references within) the considered systems are not of sequential type. Additionally, we have the exact formulas of the solutions of the systems both with the sequential differences and with the ordinary differences while in [32] the recurrence form of solution for discrete systems with Grünwald-Letnikov difference is given.

The paper is organized as follows. In Section 2 all preliminary definitions, facts, and notations are gathered. Section 3 presents systems with sequential fractional differences with results on uniqueness of solutions. We include semilinear systems in Section 4. The Section 5 concerns positivity of considered systems. Finally, the illustrative example is presented.

#### 2. Preliminaries

Let us denote by the set of all real valued functions defined on . Let , and put , where . Let Then the operator is defined by . The next two definitions of -difference operator were originally given in [16, 17].

*Definition 1. *For a function the *forward **-difference operator* is defined as
while the *-difference sum* is given by
where , and .

*Definition 2. *For arbitrary the *-factorial function* is defined by
where is the Euler gamma function, that is, for all such that , , and we use the convention that division at a pole yields zero.

Notice that if we use the general binomial coefficient , then (4) can be rewritten as In the sequel we need the following technical properties.

Proposition 3 (see [12]). *Let . *(1)*For ** one has *, *where **, ** are such that both sides are well defined; *(2)*For ** one has *;(3)*For ** one has *.

The next definition with another notation was stated in [17]. Here we use more suitable summations.

*Definition 4. *For a function the *fractional **-sum of order * is given by
where , . Moreover, we define .

It is important to notice that the operator changes the domains of functions.

*Remark 5. *Note that .

According to the definition of -factorial function the formula given in Definition 4 can be rewritten as
for , . Observe that and for we have again (3).

*Remark 6. *In [7] one can find the following form of the fractional -sum of order :
that can be useful in implementation.

The following definition can be found in [33] for or in [10] for an arbitrary .

*Definition 7. *Let . The *Caputo **-difference operator * of order for a function is defined by

*Remark 8. *Note that , where .

We need the power rule formulas in the sequel. Firstly, we easily notice that for the well-defined -factorial functions have the following property: More properties of -factorial functions can be found in [10]. In our consideration the crucial role plays the power rule formula presented in [16], that is, where , , . Note that using the general binomial coefficient one can write (11) as If , then we have for , and Let us define special functions that we use in the next section to write the formula for solutions.

*Definition 9. *For we define
where belongs to the set of integers , and .

*Remark 10. *It is worthy to notice that for (a);(b) and the values are neglected for ;(c) and as the division by pole gives zero, the formula works also for , ;(d), .

We also need the following proposition.

Proposition 11. *Let and . Then for , *

*Proof. *We show only equality (15), as (16) is a symmetric one.

Let . For we define the following -factorial function . Since
for and for , by (11) we get
where . Hence

From the application of the power rule follows the rule for composing two fractional -sums. The proof for the case can be found in [7]. For any positive the proof is presented in [10].

Proposition 12. *Let be a real valued function defined on , where , . For the following equalities hold:
**where .*

The next proposition gives a useful identity of transforming Caputo fractional difference equations into fractional summations for the case when an order is from the interval .

Proposition 13 (see [10]). *Let , , and let be a real valued function defined on . The following formula holds:
*

The operators presented above can be extended to vectors in a componentwise manner.

#### 3. Solutions of Systems with Sequential Fractional Differences

Let and . Moreover, let us take and . Then we define Note that . Then we apply the next difference operator of order on the new function and consider here an initial value problem stated by the system where , with initial values where are constant vectors from .

Solutions of the state equations (23) and (24) of the sequential fractional discrete-time system can be computed in the recursive way. From Definition 7 and by Proposition 3, (23) can be written as Since by Proposition 3 the following relation holds, one gets Repeating the same computation for (24) one gets Hence, Then using the Chu-Vandermonde identity, that is, , , , one gets Then by Proposition 3, we have Consequently, since , and for Using Proposition 3 we get for

Therefore the solution of (23) and (24) is given recursively by the following formula:

Another possibility of computing the solution of (23) and (24) is to use Proposition 13 twice and then, for , we get: Hence where . Nextly, where .

Firstly we prove the formula for the unique solution in linear case of (23) and (24): , where is a constant square matrix of degree .

Theorem 14. *The solution to the system
**with initial conditions* (25) and (26), *that is*, *and* , , *is given by the following formula:**for* .

*Proof. *Notice that for we get . Since and , we get .

For let us define the sequence in the following way:
where and with .

We calculate the first step. As , then . Going further,
which could be written as
and, using Proposition 11, we get

Taking tending to we get formula (41) as the solution of (39) and (40) with initial conditions (25) and (26).

##### 3.1. Semilinear Sequential Systems

Firstly we state a technical lemma and notations.

Lemma 15. *Let and . Let and for . Then for one gets
**
where .*

*Proof. *First let us consider the case . Then from Proposition 12 we can write
where , .

Let and . Then
Equation (46) for follows inductively.

Note that Similar to the procedure presented in the proof of Lemma 15 we can prove that for and , :

Taking and using formula (38) we can write (50) shortly in the following way:

Moreover, we can also write direct formula for values given in (38) for nonlinear problem. In fact using Definition 4 of fractional summation, formula (38) of functions and Proposition 12 we write (51) as follows: Using the power rule formula for and by Remark 10 we can write the recursive formula for the solution to nonlinear problem given by (23) and (24) and conditions (25) and (26): The given formula (53) also works for as . Then . We can check the next steps: For special semilinear case when we have . Then

Theorem 16. *The solution to the system
**
with initial conditions (25) and (26), that is, and , is given by
**
for .*

*Proof. *Note that , and for , so we get .

For based on the proof for linear case we can write the solution formula as follows:
where . Then taking into account formulas (50) and (51) we get the form (58) as the solution of (56) and (57) with initial conditions (25) and (26).

#### 4. Positivity

Let be the set of real matrices with the nonnegative entries and .

Based on [8, 31, 32] we consider the following definitions.

*Definition 17. *The fractional system (23) and (24) is called *positive* fractional system if and only if for any initial conditions .

Similarly as in [32] we will use the recursive formula (35) to show the positivity of the considered systems. Observe that the systems considered in this paper are of the sequential type while in [32] the sequential systems are not studied. Moreover, the Grünwald-Letnikov operator with the step equal to one is used in [32], whereas we study the systems with the -differences of Caputo type. So in our case the steps are equal to .

From [32] we have the following lemma.

Lemma 18. *If , then , .*

Moreover, for we have .

Using the properties of Euler gamma function one can show that if , then for , one has the following inequality:

Proposition 19. *Let . If and for all , **
then for all .*

*Proof. *The proof is by the induction principle. Assume that and both and . Since and , by (35) we have
Moreover,
where and . Let and denote the th coordinates of the vectors and , respectively. Then since for ,
Note that for we have , so . Then since and (61) holds for , we get
For we have
Since , for , and (61) holds for , we have
for , .

Now for we get Since for , , for , and (61) holds for , using (60) we have for , .

Assume that for . Using the properties of gamma function one can show that for . Then applying to (35) we get for By inequality (60) and the assumptions, that is, , , and by (61) one gets . Hence using the induction principle we get that for all .

Corollary 20. *Let and denote the identity matrix. If and for all , **
then for all that is, (39) and (40) is positive.*

*Remark 21. *In [32] the sufficient condition concerning the positivity of the linear discrete systems with Grünwald-Letnikov operator is as follows: . In our case since we have systems with the sequential fractional -difference, in our condition and both orders and appear. Note that taking the sufficient condition (71) has the form .

Let us now consider some examples that illustrate the solution of the considered systems.

*Example 22. *Let , , and . Then let us take and and consider the linear system with sequential difference in the following form:
with initial conditions (25) and (26), that is, and , . Moreover let and . Then the matrix has nonnegative values and the fact of staring from nonnegative initial conditions provides the remaining in the same positive cone. In Figure 1 we present the trajectory for steps on the -plane that start from and with the initial value . Let us notice that in the case one can see that in a few initial steps points are escaping from the positive area (see Figure 2). All calculations were done in the Maple program.

#### 5. Conclusions

The behaviour of the solutions of systems with sequential fractional difference is studied. We present recursive formulas for nonlinear systems and give the exact formulas for the unique solutions to systems in linear and semilinear cases. We prove the sufficient condition for the positivity of considered systems.

Our future goal is to study stability of systems with sequential fractional differences.

#### Acknowledgments

This project was supported by the foundayion of the National Science Centre granted on the bases of the decision. The DEC-2011/03/B/ST7/03476. This work was supported by Bialystok University of Technology Grant no. G/WM/3/2012.

#### References

- J. A. T. Machado, “Analysis and design of fractional-order digital control systems,”
*Systems Analysis Modelling Simulation*, vol. 27, no. 2-3, pp. 107–122, 1997. View at Google Scholar · View at Scopus - 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. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. M. Atıcı and P. W. Eloe, “A transform method in discrete fractional calculus,”
*International Journal of Difference Equations*, vol. 2, no. 2, pp. 165–176, 2007. View at Google Scholar · View at MathSciNet - F. M. Atıcı 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Scopus - F. Chen, X. Luo, and Y. Zhou, “Existence results for nonlinear fractional difference equation,”
*Advances in Difference Equations*, vol. 2011, Article ID 713201, 12 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. T. Holm,
*The theory of discrete fractional calculus: development and application [Ph.D. thesis]*, University of Nebraska, Lincoln, Neb, USA, 2011. - T. Kaczorek, “Fractional positive continuous-time linear systems and their reachability,”
*International Journal of Applied Mathematics and Computer Science*, vol. 18, no. 2, pp. 223–228, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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, 1988. - D. Mozyrska and E. Girejko, “Overview of the fractional $h$-difference operators,” in
*Advances in Harmonic Analysis and Operator Theory*, vol. 229 of*Operator Theory: Advances and Applications*, pp. 253–268, Springer, New York, NY, USA, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. D. Ortigueira, “Fractional discrete-time linear systems,” in
*Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '97)*, vol. 3, pp. 2241–2244, IEEE Computer Society, Munich, Germany, 1997. View at Publisher · View at Google Scholar - I. Podlubny,
*Fractional Differential Equations*, Academic Press, San Diego, Calif, USA, 1999. View at MathSciNet - T. Abdeljawad, “Dual identities in fractional difference calculus within Riemann,”
*Advances in Difference Equations*, vol. 2013, article 36, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - T. Abdeljawad, “On delta and nabla Caputo fractional differences and dual identities,”
*Discrete Dynamics in Nature and Society*, vol. 2013, Article ID 406910, 12 pages, 2013. View at Publisher · View at Google Scholar - T. Abdeljawad, D. Baleanu, F. Jarad, and R. Agarwal, “Fractional sums and difference with binomial coefficients,”
*Discrete Dynamic Systems in Nature and Society*, vol. 2013, Article ID 104173, 6 pages, 2013. View at Publisher · View at Google Scholar - R. A. C. Ferreira and D. F. M. Torres, “Fractional $h$-difference equations arising from the calculus of variations,”
*Applicable Analysis and Discrete Mathematics*, vol. 5, no. 1, pp. 110–121, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Mozyrska and E. Pawłuszewicz, “Fractional discrete-time linear control systems with initialisation,”
*International Journal of Control*, vol. 85, no. 2, pp. 213–219, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - D. Mozyrska, E. Girejko, and M. Wyrwas, “Comparison of
*h*-difference fractional operators,” 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. 191–197, Springer, New York, NY, USA, 2013. View at Google Scholar - F. Jarad, T. Abdeljawad, D. Baleanu, and K. Biçen, “On the stability of some discrete fractional nonautonomous systems,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 476581, 9 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Mozyrska and E. Pawluszewicz, “Local controllability of nonlinear discrete-time fractional order systems,”
*Bulletin of the Polish Academy of Sciences: Technical Sciences*, vol. 61, no. 1, pp. 251–256, 2013. View at Publisher · View at Google Scholar - D. Mozyrska and E. Pawluszewicz, “Controllability of
*h*-difference linear control systems with two fractional orders,”*International Journal of Systems Science*, 2013. View at Publisher · View at Google Scholar - E. Pawluszewicz and D. Mozyrska, “Constrained controllability of
*h*-difference linear systems with two fractional orders,” 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. 67–75, Springer, New York, NY, USA, 2013. View at Publisher · View at Google Scholar - M. Wyrwas, E. Girejko, D. Mozyrska, and E. Pawłuszewicz, “Stability of fractional difference systems with two 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. 41–52, Springer, New York, NY, USA, 2013. View at Publisher · View at Google Scholar - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, John Wiley & Sons, New York, NY, USA, 1st edition, 1993. View at MathSciNet - M. Klimek, “Sequential fractional differential equations with hadamard derivative,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 12, pp. 4689–4697, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. M. Furati, “A Cauchy-type problem with a sequential fractional derivative in the space of continuous functions,”
*Boundary Value Problems*, vol. 2012, article 58, 2012. View at Publisher · View at Google Scholar - D. Mozyrska, E. Girejko, and M. Wyrwas, “Fractional nonlinear systems with sequential operators,”
*Central European Journal of Physics*, 2013. View at Publisher · View at Google Scholar - T. Kaczorek, “Reachability and controllability to zero tests for standard and positive fractional discrete-time systems,”
*Journal Européan des Systemès Automatisés*, vol. 42, no. 6–8, pp. 769–787, 2008. View at Publisher · View at Google Scholar · View at Scopus - T. Kaczorek, “Fractional positive linear systems,”
*Kybernetes*, vol. 38, no. 7-8, pp. 1059–1078, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Kaczorek,
*Selected Problems of Fractional Systems Theory*, vol. 411, Springer, Berlin, Germany, 2011. - T. Abdeljawad, “On Riemann and Caputo fractional differences,”
*Computers & Mathematics with Applications*, vol. 62, no. 3, pp. 1602–1611, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet