`International Journal of Differential EquationsVolume 2012, Article ID 346089, 17 pageshttp://dx.doi.org/10.1155/2012/346089`
Research Article

## On the Solutions Fractional Riccati Differential Equation with Modified Riemann-Liouville Derivative

Department of Mathematics Engineering, Gümüşhane University, 29100 Gümüşhane, Turkey

Received 10 December 2011; Accepted 6 March 2012

Academic Editor: Ebrahim Momoniat

Copyright © 2012 Mehmet Merdan. 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

Fractional variational iteration method (FVIM) is performed to give an approximate analytical solution of nonlinear fractional Riccati differential equation. Fractional derivatives are described in the Riemann-Liouville derivative. A new application of fractional variational iteration method (FVIM) was extended to derive analytical solutions in the form of a series for these equations. The behavior of the solutions and the effects of different values of fractional order are indicated graphically. The results obtained by the FVIM reveal that the method is very reliable, convenient, and effective method for nonlinear differential equations with modified Riemann-Liouville derivative

#### 1. Introduction

In recent years, fractional calculus used in many areas such as electrical networks, control theory of dynamical systems, probability and statistics, electrochemistry of corrosion, chemical physics, optics, engineering, acoustics, viscoelasticity, material science and signal processing can be successfully modelled by linear or nonlinear fractional order differential equations . As it is well known, Riccati differential equations concerned with applications in pattern formation in dynamic games, linear systems with Markovian jumps, river flows, econometric models, stochastic control, theory, diffusion problems, and invariant embedding . Many studies have been conducted on solutions of the Riccati differential equations. Some of them, the approximate solution of ordinary Riccati differential equation obtained from homotopy perturbation method (HPM) , homotopy analysis method (HAM) , and variational iteration method proposed by He . The He’s homotopy perturbation method proposed by He  the variational iteration method  and Adomian decomposition method (ADM)  to solve quadratic Riccati differential equation of fractional order.

The variational iteration method (VIM), which proposed by He [28, 29], was successfully applied to autonomous ordinary and partial differential equations and other fields. He  was the first to apply the variational iteration method to fractional differential equations. In recent years, a new modified Riemann-Liouville left derivative is suggested by Jumarie . Recently, the fractional Riccati differential equation is solved with help of new homotopy perturbation method (HPM) .

In this paper, we extend the application of the VIM in order to derive analytical approximate solutions to nonlinear fractional Riccati differential equation: subject to the initial conditions where is fractional derivative order, is an integer, , , and are known real functions, and is a constant.

The goal of this paper is to extend the application of the variational iteration method to solve fractional nonlinear Riccati differential equations with modified Riemann-Liouville derivative.

The paper is organized as follows: In Section 2, we give definitions related to the fractional calculus theory briefly. In Section 3, we define the solution procedure of the fractional variational iteration method to show inefficiency of this method, we present the application of the FVIM for the fractional nonlinear Riccati differential equations with modified Riemann-Liouville derivative and numerical results in Section 4. The conclusions are then given in the final Section 5.

#### 2. Basic Definitions

Here, some basic definitions and properties of the fractional calculus theory which can be found in .

Definition 2.1. Assume denote a continuous (but not necessarily differentiable) function, and let the partition in the interval . Jumarie’s derivative is defined through the fractional difference : where . Fractional derivative is defined as the following limit form [1, 7]: This definition is close to the standard definition of derivatives (calculus for beginners), and as a direct result, the th derivative of a constant, , is zero.

Definition 2.2. The left-sided Riemann-Liouville fractional integral operator of order , of a function is defined as The properties of the operator can be found in [1, 7, 36].

Definition 2.3. The modified Riemann-Liouville derivative [33, 34] is defined as where .
In addition, we want to give as in the following some properties of the fractional modified Riemann-Liouville derivative.
Fractional Leibniz product law:
Fractional Leibniz Formulation:
Fractional the integration of part:

Definition 2.4. Fractional derivative of compounded functions [33, 34] is defined as

Definition 2.5. The integral with respect to [33, 34] is defined as the solution of the fractional differential equation:

Lemma 2.6. Let denote a continuous function [33, 34] then the solution of (2.5) is defined as For example, in (2.10) one obtains

Definition 2.7. Assume that the continuous function has a fractional derivative of order , for any positive integer and any , ; then the following equality holds, which is On making the substitution and , we obtain the fractional Mc-Laurin series:

#### 3. Fractional Variational Iteration Method

To describe the solution procedure of the fractional variational iteration method , we consider the following fractional Riccati differential equation: According to the VIM, we can build a correct functional for (3.1) as follows: Using (2.3), we obtain a new correction functional: It is obvious that the sequential approximations , can be established by determining , a general Lagrange’s multiplier, which can be identified optimally with the variational theory. The function is a restricted variation which means . Therefore, we first designate the Lagrange multiplier that will be identified optimally via integration by parts. The successive approximations of the solution will be readily obtained upon using the obtained Lagrange multiplier and by using any selective function . The initial values are usually used for choosing the zeroth approximation . With determined, then several approximations , follow immediately . Consequently, the exact solution may be procured by using

#### 4. Applications

In this section, we present the solution of two examples of the Riccati differential equations as the applicability of FVIM.

Example 4.1. Let usconsider the fractional Riccati differential equation, we get with initial conditions: Construction the following functional: we have Similarly, we can get the coefficients of to zero: The generalized Lagrange multiplier can be identified by the above equations: substituting (4.6) into (4.3) produces the iteration formulation as follows: Taking the initial value , we can derive Then, the approximate solutions in a series form are As is The exact solution of (4.1) is , when .
Figure 1 indicates the solution obtained using FVIM versus the exact solution when . Figure 2 is plotted for approximate solution of time-fractional Riccati differential equation for , 0.8, 0.9, and 1. Equation (4.1) is solved by using the homotopy perturbation method (HPM) . FVIM solutions indicate that the present algorithm performs extreme efficiency, simplicity, and reliability. The results obtained from FVIM are fully compatible with those of the HPM.
Table 1 shows the approximate solutions for (4.1) obtained for different values of using the variational iteration method and HPM . From the numerical results in Table 1, it is clear that the approximate solutions are in high agreement with the exact solutions, when , and the solution continuously depends on the time-fractional derivative. Example 4.1 has been solved using HAM , ADM , VIM , and HPM .

Figure 2: Plots of approx. solution for different values of .

Example 4.2. Let usconsider the fractional Riccati differential equation, we get with initial conditions Construction the following functional: we have Similarly, we can get the coefficients of to zero: The generalized Lagrange multiplier can be identified by the above equations: substituting (4.16) into (4.13) produces the iteration formulation as follows: Taking the initial value , we can derive Then, the approximate solutions in a series form are As is The exact solution of (4.11) is , when .
Figure 3 is plotted for approximate solution of time-fractional Riccati differential equation found in Example 4.2. In Figure 4, we have shown the graphic of approximate solution of (4.11) for . Figures 2 and 4 show that a decrease in the fractional order causes an increase in the function.
Table 2 indicates the approximate solutions for (4.11) obtained for different values of using the variational iteration method and HPM . From the numerical results in Table 2, it is clear that the approximate solutions are in high agreement with the exact solutions, when , and the solution continuously depends on the time-fractional derivative.

Figure 4: Plots of approx. solution for different values of .

Example 4.3. Let usconsider the fractional Riccati differential equation, we get with initial conditions: Construction the following functional: we have Similarly, we can get the coefficients of to zero: The generalized Lagrange multiplier can be identified by the above equations: substituting (4.26) into (4.23) produces the iteration formulation as follows: Taking the initial value , we can derive Then, the approximate solutions in a series form are As is The exact solution of (4.21) is where is the Bessel function of first kind, when .
Figure 5 is plotted for approximate solution of time-fractional Riccati differential equation found in Example 4.3. In Figure 6, we have shown the graphic of approximate solution of (4.21) for . Figures 2, 4, and 6 show that a decrease in the fractional order causes an increase in the function.
Table 3 indicates the approximate solutions for (4.21) obtained for different values of using the HPM . From the numerical results in Table 3, it is clear that the approximate solutions are in substantial agreement with the exact solutions, when , and the solution continuously depends on the time-fractional derivative.

Figure 6: Plots of approx. solution for different values of .

#### 5. Conclusions

In this paper, variational iteration method having integral w.r.t. has been successfully implemented to finding approximate analytical solution of fractional Riccati differential equations. Variational iteration method known as very powerful and an effective method for solving fractional Riccati differential equation. It is also a promising method to solve other nonlinear equations. In this paper, we have discussed modified variational iteration method having integral w.r.t. used for the first time by Jumarie. The obtained results indicate that this method is powerful and meaningful for solving the nonlinear fractional differential equations. Three examples indicate that the results of variational iteration method having integral w.r.t. are in excellent agreement with those obtained by HPM, ADM, and HAM, which is available in the literature.

#### References

1. K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
2. I. Podlubny, Fractional Differential Equations, Academic Press, New York, NY, USA, 1999.
3. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, The Netherlands, 2006.
4. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.
5. M. Caputo, “Linear models of dissipation whose Q is almost frequency independent, part II,” Geophysical Journal International, vol. 13, no. 5, pp. 529–539, 1967.
6. A. A. Kilbas, H. H. Srivastava, and J. J. Trujillo, Theoryand Applications of Fractional Differential Equations, Elsevier, Amsterdam, The Netherlands, 2006.
7. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, New York, NY, USA, 1993.
8. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach, Yverdon, Switzerland, 1993.
9. G. M. Zaslavsky, Hamiltonian Chaosand Fractional Dynamics, Oxford University Press, 2005.
10. M. Merdan, A. Yıldırım, and A. Gökdoğan, “Numerical solution of time-fraction Modified Equal Width Wave Equation,” Engineering Computations. In press.
11. G. A. Einicke, L. B. White, and R. R. Bitmead, “The use of fake algebraic Riccati equations for co-channel demodulation,” IEEE Transactions on Signal Processing, vol. 51, no. 9, pp. 2288–2293, 2003.
12. M. Gerber, B. Hasselblatt, and D. Keesing, “The riccati equation: pinching of forcing and solutions,” Experimental Mathematics, vol. 12, no. 2, pp. 129–134, 2003.
13. R. E. Kalman, Y. C. Ho, and K. S. Narendra, “Controllability of linear dynamical systems,” Contributions to Differential Equations, vol. 1, pp. 189–213, 1963.
14. S. Bittanti, P. Colaneri, and G. De Nicolao, “The periodic Riccati equation,” in The Riccati Equation, Communications and Control Engineering, pp. 127–162, Springer, Berlin, Germany, 1991.
15. S. Bittanti, P. Colaneri, and G. O. Guardabassi, “Periodic solutions of periodic Riccati equations,” IEEE Transactions on Automatic Control, vol. 29, no. 7, pp. 665–667, 1984.
16. B. D. O. Anderson and J. B. Moore, Optimal Filtering, Prentice-Hall, Englewood Cliffs, NJ, USA, 1979.
17. W. T. Reid, Riccati Differential Equations: Mathematics in Science and Engineering, vol. 86, Academic Press, New York, NY, USA, 1972.
18. H. Aminikhah and M. Hemmatnezhad, “An efficient method for quadratic Riccati differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 4, pp. 835–839, 2010.
19. S. Abbasbandy, “Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian's decomposition method,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 485–490, 2006.
20. S. Abbasbandy, “Iterated He's homotopy perturbation method for quadratic Riccati differential equation,” Applied Mathematics and Computation, vol. 175, no. 1, pp. 581–589, 2006.
21. Y. Tan and S. Abbasbandy, “Homotopy analysis method for quadratic Riccati differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 3, pp. 539–546, 2008.
22. S. Abbasbandy, “A new application of He's variational iteration method for quadratic Riccati differential equation by using Adomian's polynomials,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 59–63, 2007.
23. N. A. Khan, A. Ara, and M. Jamil, “An efficient approach for solving the Riccati equation with fractional orders,” Computers and Mathematics with Applications, vol. 61, no. 9, pp. 2683–2689, 2011.
24. Z. Odibat and S. Momani, “Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order,” Chaos, Solitons and Fractals, vol. 36, no. 1, pp. 167–174, 2008.
25. J. Cang, Y. Tan, H. Xu, and S. J. Liao, “Series solutions of non-linear Riccati differential equations with fractional order,” Chaos, Solitons and Fractals, vol. 40, no. 1, pp. 1–9, 2009.
26. H. Jafari and H. Tajadodi, “He’s variational iteration method for solving fractional Riccati differential equation,” International Journal of Differential Equations, vol. 2010, Article ID 764738, 8 pages, 2010.
27. S. Momani and N. Shawagfeh, “Decomposition method for solving fractional Riccati differential equations,” Applied Mathematics and Computation, vol. 182, no. 2, pp. 1083–1092, 2006.
28. J. H. He, “Variational iteration method—a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.
29. J. H. He and X. H. Wu, “Variational iteration method: new development and applications,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.
30. J. H. He, “Some applications of nonlinear fractional differential equations and their approximations,” Bulletin of Science, Technology & Society, vol. 15, no. 2, pp. 86–90, 1999.
31. G. Jumarie, “Stochastic differential equations with fractional Brownian motion input,” International Journal of Systems Science, vol. 24, no. 6, pp. 1113–1131, 1993.
32. G. Jumarie, “New stochastic fractional models for Malthusian growth, the Poissonian birth process and optimal management of populations,” Mathematical and Computer Modelling, vol. 44, no. 3-4, pp. 231–254, 2006.
33. G. Jumarie, “Laplace's transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville derivative,” Applied Mathematics Letters, vol. 22, no. 11, pp. 1659–1664, 2009.
34. G. Jumarie, “Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functions,” Applied Mathematics Letters, vol. 22, no. 3, pp. 378–385, 2009.
35. G. Jumarie, “On the solution of the stochastic differential equation of exponential growth driven by fractional Brownian motion,” Applied Mathematics Letters, vol. 18, no. 7, pp. 817–826, 2005.
36. M.-J. Jang, C.-L. Chen, and Y.-C. Liu, “Two-dimensional differential transform for partial differential equations,” Applied Mathematics and Computation, vol. 121, no. 2-3, pp. 261–270, 2001.
37. N. Faraz, Y. Khan, H. Jafari, A. Yildirim, and M. Madani, “Fractional variational iteration method via modified Riemann-Liouville derivative,” Journal of King Saud University, vol. 23, no. 4, pp. 413–417, 2010.