Abstract and Applied Analysis

Volume 2014 (2014), Article ID 970967, 6 pages

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

## Reproducing Kernel Method for Fractional Riccati Differential Equations

^{1}Department of Mathematics, Changshu Institute of Technology, Suzhou, Jiangsu 215500, China^{2}Department of Mathematics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China

Received 29 December 2013; Accepted 7 April 2014; Published 27 April 2014

Academic Editor: Youyu Wang

Copyright © 2014 X. Y. Li 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

This paper is devoted to a new numerical method for fractional Riccati differential equations. The method combines the reproducing kernel method and the quasilinearization technique. Its main advantage is that it can produce good approximations in a larger interval, rather than a local vicinity of the initial position. Numerical results are compared with some existing methods to show the accuracy and effectiveness of the present method.

#### 1. Introduction

This paper deals with the numerical solution of the following fractional Riccati differential equation: where denotes the Caputo fractional derivative of order and Riccati differential equations arise in many fields [1]. The problem has attracted much attention and has been studied by many authors. However, deriving its analytical solution in an explicit form seems to be unlikely except for certain special situations. Recently, many numerical methods [2–9] have been proposed to solve integer order Riccati differential equations. However, the discussion on the numerical methods for fractional order Riccati differential equations is rare. Odibat and Momani [10] developed a modified homotopy perturbation method for fractional Riccati differential equations. Li [11] presented a numerical method for fractional differential equations based on Chebyshev wavelets. Hosseinnia et al. [12] applied an enhanced homotopy perturbation method for fractional Riccati differential equations. Yüzbaşı [13] introduced a numerical method for fractional Riccati differential equations using the Bernstein polynomial. Khader [14] developed the fractional Chebyshev finite difference method for fractional Riccati differential equations. Yang and Baleanu [15], Yang et al. [16], and Baleanu et al. [17] proposed local fractional variation iteration for fractional heat conduction and wave equations.

Recently, based on the reproducing kernel theory, Cui, Geng, and Lin presented the reproducing kernel method (RKM) for linear and nonlinear operator equations [18–21]. The method has been developed and applied to many problems [22–26].

The aim of this paper is to present a new method for fractional Riccati differential equations, based on the RKM and the quasilinearization technique.

The rest of the paper is organized as follows. In the next section, the quasilinearization technique is applied to fractional Riccati differential equation. The RKM for reduced linear fractional differential equations is introduced in Section 3. The numerical examples are presented in Section 4. Section 5 ends this paper with a brief conclusion.

#### 2. Quasilinearization of Riccati Differential Equation (1)

In this section, the quasilinearization technique is applied to reduce (1) to a series of linear fractional problems. Define . By choosing an appropriate initial approximation for the function in and expanding around , it follows that Generally, one can write for ( = iteration index) Therefore, the following iteration formula for (1) can be derived: where and , and is the initial approximation.

Clearly, to solve (1), it suffices for us to solve the series of linear problem (5).

#### 3. Method for Solving Linear Fractional Problem (5)

To illustrate how to solve (5) we consider the problem of solving where and are continuous.

By the definition of Caputo fractional derivative, (6) is equivalent to the following equation: If , then the improper integral exists and Define This gives us a continuous function on and then integral also exists and .

Therefore, (6) can be converted into

Applying Hermite's quadrature formula to , one obtains where , is the number of nodes, and , .

Then (6) can be further equivalently approximated to

To apply the RKM to (12), it is necessary to construct the following reproducing kernel Hilbert space .

*Definition 1. * is an absolutely continuous real value function, . The inner product and norm in are given, respectively, by

Theorem 2. * is a reproducing kernel space and its reproducing kernel is
**
where .*

*Definition 3. * is an absolutely continuous real value function, . The inner product and norm in are given, respectively, by

Theorem 4. * is a reproducing kernel space and its reproducing kernel is
*

*Put
Clearly, is a bounded linear operator. Put and , where is the adjoint operator of . The orthonormal system of can be derived from Gram-Schmidt orthogonalization process of ,
*

*Theorem 5. If is dense on , then is the complete system of .*

*Theorem 6. If is dense on and the solution of (12) is unique, then the solution of (12) is
*

*Now, an approximate solution of (6) can be obtained by the -term intercept of the exact solution and
*

*Similarly, the approximate solutions can be obtained:
where .*

*4. Numerical Examples*

*4. Numerical Examples*

*Example 1. *Consider the following fractional Riccati differential equation [10–13]:
The exact solution for can be easily determined to be
Applying the proposed method, taking , the numerical results compared with other methods are listed in Tables 1 and 2. Taking , the numerical results on [] are listed in Table 3. From Table 3, it is easily found that the present approximations are effective for a larger interval, rather than a local vicinity of the initial position.

*Example 2. *Consider the following fractional Riccati differential equation [10–14]:
The exact solution for can be easily determined to be
According to the present method, taking , the numerical results compared with other methods are given in Tables 4, 5, and 6. Taking , the numerical results on [] are shown in Figures 1 and 2. From these figures we can conclude that the obtained numerical solutions are in excellent agreement with the exact solution for a larger interval.

*5. Conclusion*

*5. Conclusion*

*In this paper, combining the RKM, the numerical integral, and quasilinearization techniques, a new numerical method is proposed for fractional Riccati differential equations. The main advantage of this method is that it can provide accurate numerical approximations on a larger interval. Numerical results compared with the existing methods show that the present method is a powerful method for solving fractional Riccati differential equations.*

*Conflict of Interests*

*Conflict of Interests*

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

*Acknowledgment*

*Acknowledgment*

*This work was supported by the National Natural Science Foundation of China (Grant nos. 11326237, 11271100, and 11126222).*

*References*

*References*

- W. T. Reid,
*Riccati Differential Equations*, Academic Press, New York, NY, USA, 1972. View at MathSciNet - M. A. El-Tawil, A. A. Bahnasawi, and A. Abdel-Naby, “Solving Riccati differential equation using Adomian's decomposition method,”
*Applied Mathematics and Computation*, vol. 157, no. 2, pp. 503–514, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Lakestani and M. Dehghan, “Numerical solution of Riccati equation using the cubic B-spline scaling functions and Chebyshev cardinal functions,”
*Computer Physics Communications*, vol. 181, no. 5, pp. 957–966, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Z. Geng, Y. Z. Lin, and M. G. Cui, “A piecewise variational iteration method for Riccati differential equations,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2518–2522, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B.-Q. Tang and X.-F. Li, “A new method for determining the solution of Riccati differential equations,”
*Applied Mathematics and Computation*, vol. 194, no. 2, pp. 431–440, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Ghorbani and S. Momani, “An effective variational iteration algorithm for solving Riccati differential equations,”
*Applied Mathematics Letters*, vol. 23, no. 8, pp. 922–927, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Mohammadi and M. M. Hosseini, “A comparative study of numerical methods for solving quadratic Riccati differential equations,”
*Journal of the Franklin Institute. Engineering and Applied Mathematics*, vol. 348, no. 2, pp. 156–164, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Odibat and S. Momani, “Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order,”
*Chaos, Solitons & Fractals*, vol. 36, no. 1, pp. 167–174, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. L. Li, “Solving a nonlinear fractional differential equation using Chebyshev wavelets,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 15, no. 9, pp. 2284–2292, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. H. Hosseinnia, A. Ranjbar, and S. Momani, “Using an enhanced homotopy perturbation method in fractional differential equations via deforming the linear part,”
*Computers & Mathematics with Applications*, vol. 56, no. 12, pp. 3138–3149, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yüzbaşı, “Numerical solutions of fractional Riccati type differential equations by means of the Bernstein polynomials,”
*Applied Mathematics and Computation*, vol. 219, no. 11, pp. 6328–6343, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. M. Khader, “Numerical treatment for solving fractional Riccati differential equation,”
*Journal of the Egyptian Mathematical Society*, vol. 21, no. 1, pp. 32–37, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. J. Yang and D. Baleanu, “Fractal heat conduction problem solved by local fractional variation iteration method,”
*Thermal Science*, vol. 17, no. 2, pp. 625–628, 2013. View at Google Scholar - X.-J. Yang, D. Baleanu, and W.-P. Zhong, “Approximate solutions for diffusion equations on Cantor space-time,”
*Proceedings of the Romanian Academy, A: Mathematics, Physics, Technical Sciences, Information Science*, vol. 14, no. 2, pp. 127–133, 2013. View at Google Scholar · View at MathSciNet - D. Baleanu, J. A. Tenreiro Machado, C. Cattani, M. C. Baleanu, and X.-J. Yang, “Local fractional variational iteration and decomposition methods for wave equation on Cantor sets within local fractional operators,”
*Abstract and Applied Analysis*, vol. 2014, Article ID 535048, 6 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - M. G. Cui and Y. Z. Lin,
*Nonlinear Numerical Analysis in Reproducing Kernel Space*, Nova Science Publishers, 2009. View at MathSciNet - F. Z. Geng and M. G. Cui, “Solving a nonlinear system of second order boundary value problems,”
*Journal of Mathematical Analysis and Applications*, vol. 327, no. 2, pp. 1167–1181, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Z. Geng, “A numerical algorithm for nonlinear multi-point boundary value problems,”
*Journal of Computational and Applied Mathematics*, vol. 236, no. 7, pp. 1789–1794, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Z. Geng and S. P. Qian, “Reproducing Kernel method for singularly perturbed turning point problems having twin boundary layers,”
*Applied Mathematics Letters*, vol. 26, no. 10, pp. 998–1004, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. Inc, A. Akgül, and A. Kiliçman, “A novel method for solving KdV equation based on reproducing Kernel Hilbert space method,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 578942, 11 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Mohammadi and R. Mokhtari, “Solving the generalized regularized long wave equation on the basis of a reproducing Kernel space,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 14, pp. 4003–4014, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Akram and H. U. Rehman, “Numerical solution of eighth order boundary value problems in reproducing Kernel space,”
*Numerical Algorithms*, vol. 62, no. 3, pp. 527–540, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Abu Arqub, M. Al-Smadi, and S. Momani, “Application of reproducing Kernel method for solving nonlinear Fredholm-Volterra integrodifferential equations,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 839836, 16 pages, 2012. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Y. Wang, B. Han, and M. Yamamoto, “Inverse heat problem of determining time-dependent source parameter in reproducing Kernel space,”
*Nonlinear Analysis. Real World Applications*, vol. 14, no. 1, pp. 875–887, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet

*
*