Abstract and Applied Analysis

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

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

Research Article

Reproducing Kernel Method for Fractional Riccati Differential Equations

X. Y. Li,¹ B. Y. Wu,² and R. T. Wang¹

¹Department of Mathematics, Changshu Institute of Technology, Suzhou, Jiangsu 215500, China

²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

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.

Table 1: Comparison of the numerical solutions with the other methods for .

Ours
[10]
[11]

0.20
0.469516
0.428892
0.584307

0.40
0.933596
0.891404
1.024974

0.50
1.14488
1.132763
1.198621

0.60
1.33098
1.370240
1.349150

0.80
1.62153
1.794879
1.599235

1.00
1.81865
2.087384
1.801763

Table 2: Comparison of the numerical solutions with the other methods for .

Ours
[13]
IABMM [12]
MHPM [10]

0.20
0.312985
0.314869
—
—

0.40
0.695357
0.697544
—
—

0.50
0.901484
0.903695
0.8621
0.9010

0.60
1.10576
1.107866
—
—

0.80
1.47606
1.477707
—
—

1.00
1.76417
1.764520
1.7356
1.8720

Table 3: Numerical results for on .

Ours ()
Ours ()
Exact solution ()

0.5
0.769552
0.756446
0.756014

1.0
1.69828
1.68978
1.689500

1.5
2.19159
2.19599
2.195630

2.0
2.35063
2.35838
2.357770

2.5
2.39373
2.40046
2.40028

3.0
2.40539
2.41058
2.41081

3.5
2.40899
2.41302
2.41339

4.0
2.41047
2.41374
2.41401

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.

Table 4: Comparison of the numerical solutions with the other methods for .

Ours
[13]
IABMM [12]
MHPM [10]

0.20
0.307359
0.3099755
0.3117
0.3138

0.40
0.480346
0.4816318
0.4855
0.4929

0.60
0.597542
0.5977827
0.6045
0.5974

0.80
0.679657
0.6788495
0.6880
0.6604

1.00
0.738213
0.7368368
0.7478
0.7183

Table 5: Comparison of the numerical solutions with the other methods for .

Ours
[13]
IABMM [12]
MHPM [10]

0.20
0.237652
0.2387891
0.2393
0.2391

0.40
0.421766
0.4225830
0.4234
0.4229

0.60
0.565673
0.5661715
0.5679
0.5653

0.80
0.674464
0.6746270
0.6774
0.6740

1.00
0.754632
0.7545890
0.7584
0.7569

Table 6: Comparison of the numerical solutions with the other methods for .

Exact
Ours
[13]
MHPM [10]

0.20
0.197375
0.19738
0.197375
0.197375

0.40
0.379949
0.379956
0.379948
0.379944

0.60
0.53705
0.537061
0.537049
0.536857

0.80
0.664037
0.664053
0.664036
0.661706

1.00
0.761594
0.761618
0.761594
0.746032

(a)

(b)

(c)

Figure 1: The behavior of approximate solution with different values of ((a) ; (b) ; (c) ).

(a)

(b)

Figure 2: Comparison of approximate solutions with the exact solutions for ((a) exact solution; (b) absolute errors).

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

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

Acknowledgment

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

References

W. T. Reid, Riccati Differential Equations, Academic Press, New York, NY, USA, 1972.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

M. G. Cui and Y. Z. Lin, Nonlinear Numerical Analysis in Reproducing Kernel Space, Nova Science Publishers, 2009.

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.

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.

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.

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.

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.

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.

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.

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.