Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 272198, 7 pages

http://dx.doi.org/10.1155/2014/272198

Research Article

The Oscillation of a Class of the Fractional-Order Delay Differential Equations

Qianli Lu and Feng Cen

Department of Control Science and Engineering, Tongji University, Shanghai 200092, China

Received 5 June 2014; Accepted 26 June 2014; Published 3 August 2014

Academic Editor: Sanling Yuan

Copyright © 2014 Qianli Lu and Feng Cen. 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.

Several oscillation results are proposed including necessary and sufficient conditions for the oscillation of fractional-order delay differential equations with constant coefficients, the sufficient or necessary and sufficient conditions for the oscillation of fractional-order delay differential equations by analysis method, and the sufficient or necessary and sufficient conditions for the oscillation of delay partial differential equation with three different boundary conditions. For this, α-exponential function which is a kind of functions that play the same role of the classical exponential functions of fractional-order derivatives is used.

1. Introduction

For the past three centuries, fractional calculus has been dealt with by mathematicians. In the recent decades, it is widely used in the fields of engineering, science, and economy. Recently, fractional differential systems have gained scholar’s attention [1–3]. Many researchers demonstrated applications of fractional calculus in the frequency dependent damping behavior of many viscoelastic materials, dynamics of interfaces between nanoparticles and substrates, the nonlinear oscillation of earthquakes, bioengineering, continuum and statistical mechanics, signal processing, filter design, circuit theory, and robotics. In some practical systems, such as economic systems, biological systems, and space-light industry systems, due to the transmission of the signal or the mechanical transmission, the research on fractional differential systems with delay is desired.

In [2], the authors described the application of the background and significance of the differential delay equations oscillation and fractional differential equations. The delay differential equation is different from ordinary differential equations. Moreover, compared with the periodic solution, oscillatory solution has more extensive application prospect and theory value in terms of the general solution of equation. In this paper, the necessary or sufficient conditions for the oscillation of several fractional-order delay differential equations are analyzed and discussed. Some oscillation results are proposed for the fractional-order linear delay differential equation as follows:

where . Some oscillation results are also given for the delay parabolic differential equation as follows:

For the purpose, known methods are expanded for our results.

2. Preliminaries

We first give some definitions and make some preliminaries. The fractional-order derivative can be defined in different ways, wherein Riemann-Liouville’s definition and Caputo’s definition are two main expressions as follows.

Definition 1. Riemann-Liouville’s fractional order derivative [3, 4] is defined by

where is the Gamma function.

Caputo’s fractional order derivative is defined by

Definition 2. -exponential function is as follows [2]:

where satisfies the differential equation , .

The derivative has the property of

Mittag-Lefflerfunction [3, 4] can be defined by

The first order delay differential equation including fractional-order derivative [2] is expressed by

where , . Since is a positive rational number, let satisfy . The expression for (1) with constant coefficients is given by

where

Equation (10) is the characteristic polynomial of (8). The contribution in [2] is considered.

Lemma 3. The following statements are equivalent.

(1)Every solution of (8) oscillates.(2)The characteristic equation (10) has no real roots.

In this paper, the main results are derived from the application of the basic theory in [2]. The oscillation or nonoscillation of solutions for several fractional-order delay differential equations can be obtained based on these main results.

3. Main Result

3.1. The Fractional-Order Delay Differential Equation with Constant Coefficients

The general form of the fractional-order linear delay differential equation is

Since , the characteristic polynomial of (11) is as follows:

The detailed derivation can be found in Appendices A, B, and C. Using Lemma 3, we have the necessary and sufficient conditions for oscillation of all solutions of (11) as follows.

Corollary 4. The following statements are equivalent.

(1)Every solution of (11) oscillates.(2)The characteristic equation (12) has no real roots.

In order to use Corollary 4 more effectively, Theorem 6 is proposed. Before presenting the theorem, the function is given (see Appendices A, B, and C for details), and the maximum of is as follows:

Lemma 5.

(1)If or , (12) has a real root.(2)If , (12) has not a real root.

Proof. (1) If , it is obvious that (12) has a real root.

As , from the characteristic equation (12),

So or , and (12) has a real root.

(2) In the same way, we get , and (12) has not a real root.

Theorem 6.

(1)If or , the solution of (11) is nonoscillatory.(2)If , the solution of (11) is oscillatory.

Proof. The equation has a solution , if and only if is a root of characteristic of equation . Using the conditions for (11), we have the following.

(1)If or , (12) has a real root, and the solution of (11) is nonoscillatory.(2)If , (12) does not have a real root, and the solution is oscillatory. So we can obtain , and the solutions of (11) are oscillatory. Therefore, the proof of the theorem can be completed.

Example 7. We consider delay differential equation with fractional-order derivatives

All the conditions of Theorem 6 are satisfied. Fractional derivative of some special functions is as . Thus, all the solutions of (15) are oscillatory. One of such solutions is .

We consider the following fractional-order linear delay differential equation:

: , .

Theorem 8. Assume that holds, and are satisfied. Every solution of the delay differential equation (17) oscillates.

Proof. We may assume that (17) has eventually positive solution. Then, there exists , and if , , , and .

We consider the eventually positive solution of the following inequality:

Assume the inequality (18) has a solution :

If , the inequality does not have a real root. That means every solution of the inequality (18) is oscillatory. Thus, if , every solution of the delay differential equation (17) oscillates.

3.2. The Fractional-Order Delay Differential Equation with Oscillating Coefficients

We consider oscillation of solutions of the fractional-order delay differential equation with oscillating coefficients:

Let . If is sufficiently large,

We have the sufficient conditions for oscillation of all solutions of (20) as follows.

Theorem 9. If is sufficiently large, . Then, the solution of (20) is oscillatory.

Proof. Assume that (20) has eventually positive solution. If there exists such that , then , ,

Similar to the proof process of Theorem 8, every solution of the inequality (24) is oscillatory due to . Thus, the inequality (24) does not have eventually positive solution. Therefore, if is sufficiently large, , and every solution of (20) is oscillatory. The proof of the theorem is complete.

We consider the following fractional-order linear delay differential equations:

Their characteristic equations are expressed by (27) and (28), respectively:

where , , , , , are nonnegative constants.

Based on Lemma 3, the statements are equivalent for oscillation of all solutions of (25) and (26) as follows.

Corollary 10. The following statements are equivalent.

(1)Every solution of (25) and (26) oscillates.(2)The characteristic equations of (27) and (28) have no real roots.

4. The Fractional-Order Delay Partial Differential Equation

To research the oscillation properties of solutions for delay partial differential equations, we consider the delay parabolic differential equations [5]:

where . is the bounded region with piecewise smooth in .

The hypotheses are always true as follows.

 , and is a nonnegative constant.  are nonnegative constants.

Consider the boundary conditions as follows:

where is the unit exterior normal vector in .

If there are no special instructions in this section, we only discuss the oscillation of solutions of the fractional-order delay partial differential equation with respect to time variable .

Theorem 11. Assume that and hold, and is satisfied.

(1)As or , the solution of (29) and is nonoscillatory.(2)As , the solution of (29) and is oscillatory.

Proof. Let . Equation (29) on both sides of the integral of in is expressed by

By using Green’s formula and boundary conditions ,

Let . Substituting (30) with (31),

Based on Corollary 10, (1) and (2) in Theorem 14 are proved.

Example 12. Consider the parabolic equation:

Because of , all the conditions in Theorem 11 (1) are satisfied. By using Theorem 11, every solution of (33), with

oscillates in . One of such solutions is .

For the more general equation of (29) with coefficients,

There are the boundary conditions as follows:

where is the unit exterior normal vector in .

Assume that

    , , , are nonnegative constants.

Some sufficient conditions for oscillation of all solutions of (35) with and are established in Theorem 15. Before presenting the theorem, the lemma is given as follows.

Lemma 13. For the Dirichlet problem,

where is a constant. The smallest eigenvalue of problem (36) is positive, and the corresponding eigenfunction is also positive on [5]. With each solution of problem (35) with the boundary condition , the function is defined as

Theorem 14. Assume that and hold, and is satisfied. Every solution of the delay differential equation (26) with several coefficients oscillates. Then, every solution of (35) oscillates.

Proof. Let . If we multiply both sides of (35) by and integrate in , we obtain

where , which satisfies Lemma 13.

By using Green’s formula and boundary condition , we obtain

where is the surface element on . Substituting with (37), we obtain

Based on Corollary 10, the proof of the theorem can be immediately completed. Let be defined as , , where denotes a solution of problem (35) with the boundary condition or . Then, we obtain

Theorem 15. Assume that hold and is satisfied. Every solution of the delay differential equation with several coefficients (25) oscillates. Then, every solution of (35) with the boundary condition oscillates.

Proof. Let on the be the definite integral in . Assume that (35) with the boundary condition has a nonoscillatory solution and (35) has an eventually positive solution:

By using Green’s formula and boundary condition , we obtain

where is the surface element on .

Let () be substituted into (42). We obtain

The detail of the proof process is similar to the proof process of Theorem 15.

Theorem 16. Assume that and hold and is satisfied. Every solution of the delay differential equation with the coefficients (25) oscillates. Then, every solution of (35) with the boundary condition oscillates.

Proof. Let on the be the definite integral in . For the sake of contradiction, assume that (35) with the boundary condition has a nonoscillatory solution and (35) with the boundary condition has an eventually positive solution:

By using Green’s formula and boundary condition , we obtain

where is the surface element on . Let () be substituted into (45). We obtain

The detail of the proof process is similar to the proof process of Theorem 15.

The delay parabolic differential equations are expressed by

Consider the boundary conditions, as shown in where is the unit exterior normal vector in .

  are nonnegative constants. , ; ; where is a nonnegative constant. , and , .

Theorem 17. Assume that , , and hold and is satisfied. If , the solution of (48) with the boundary condition is oscillatory.

Proof. Let . Assume that (48) with the boundary condition has no oscillation solutions, and (48) has eventually positive solution. If there exists such that , then ,  ,  . If is satisfied, (48) can be transformed into the inequality as follows:

The inequality (49) on both sides of the integral of in is expressed by

By using Green’s formula and boundary conditions , we obtain

Substituting (50) with (51) and , we obtain

By using Theorem 9, the proof of the theorem can be completed.

Appendices

A. Maximum of

Figure 1 illustrates the curves of the function with different . From bottom to top, , , , , , and 1, respectively. From Figure 1, it can be observed that there exists such that   . In detail, the values of with different are reported in Table 1. is the maximum for a fixed in Table 1.

Table 1: Values of with different .

0.0881072
0.1176864
0.1774485
0.2383790
0.2694375
0.3659126

0.0887969
0.1185872
0.1787170
0.2399068
0.2710383
0.3674039

0.089205
0.1191138
0.1794311
0.2407091
0.2718332
0.3678794 (m)

0.0892548
0.1191774
0.1795119
0.2407885
0.2719021
0.3678611

0.0892949
0.1192275
0.1795730
0.240842
0.2719425
0.3678068

0.0893251
0.1192646
0.1796147
0.2408705
0.2719547 (m)
0.3677171

0.0893456
0.1192889
0.1796375
0.2408740 (m)
0.2719395
0.3675928

0.0893567
0.1193007 (m)
0.1796417 (m)
0.2408533
0.2718973
0.3674346

0.0893585 (m)
0.1193001
0.1796277
0.2408087
0.2718286
0.3672431

0.0893512
0.1192875
0.1795959
0.2407409
0.2717341
0.3670191

0.0893351
0.1192631
0.1795465
0.2406501
0.2716142
0.3667631

0.0893102
0.1192271
0.1794801
0.2405370
0.2714694
0.3664759

0.0892768
0.1191797
0.1793968
0.2404020
0.2713003
0.3661581

0.0889875
0.1187809
0.17874
0.2394135
0.2701071
0.364132285

0.0885098
0.1181318
0.1777141
0.2379416
0.2683777
0.361433054

Figure 1: Curves of function () with different .

B. Supplementary Comment of

The curves of function and () are showed in Figure 2. The function is a supplementary comment of the maximum value of . From top to bottom, , , , , , and , respectively.

Figure 2: Curves of function and .

C. The Characteristic Equation

The Laplace transforms for two fractional derivatives with [6] are expressed as follows:

The fractional Laplace transform [2] is . There exist and such that . By taking the Laplace transforms on both sides of (11),

where . denotes the abscissa of convergence of ; that is, .

In the most general form , we can obtain and . and are both entire functions. Assume for all real . can be expressed as , . Therefore, is the characteristic polynomial of (11).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors thank Mr. Weisheng Xu for useful discussions on the proof process in the paper.

References

J. Wei, “The constant variation formulae for singular fractional differential systems with delay,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1184–1190, 2010.

Y. Bolat, “On the oscillation of fractional-order delay differential equations with constant coefficients,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 11, pp. 3988–3993, 2014.

A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,” Chaos, Solitons and Fractals, vol. 34, no. 5, pp. 1473–1481, 2007.

J. J. Trujillo, M. Rivero, and B. Bonilla, “On a Riemann-Liouville generalized Taylor's formula,” Journal of Mathematical Analysis and Applications, vol. 231, no. 1, pp. 255–265, 1999.

I. Kubiaczyk and S. H. Saker, “Oscillation of delay parabolic differential equations with several coefficients,” Journal of Computational and Applied Mathematics, vol. 147, no. 2, pp. 263–275, 2002.

M. S. Tavazoei and M. Haeri, “Chaotic attractors in incommensurate fractional order systems,” Physica D, vol. 237, no. 20, pp. 2628–2637, 2008.