Journal of Applied Mathematics

Volume 2014 (2014), Article ID 257049, 8 pages

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

Research Article

A Coiflets-Based Wavelet Laplace Method for Solving the Riccati Differential Equations

Xiaomin Wang

School of Engineering, Huazhong Agricultural University, Wuhan, Hubei 430070, China

Received 28 January 2014; Accepted 1 July 2014; Published 15 July 2014

Academic Editor: Changbum Chun

Copyright © 2014 Xiaomin Wang. 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.

A wavelet iterative method based on a numerical integration by using the Coiflets orthogonal wavelets for a nonlinear fractional differential equation is proposed. With the help of Laplace transform, the fractional differential equation was converted into equivalent integral equation of convolution type. By using the wavelet approximate scheme of a function, the undesired jump or wiggle phenomenon near the boundary points was avoided and the expansion constants in the approximation of arbitrary nonlinear term of the unknown function can be explicitly expressed in finite terms of the expansion ones of the approximation of the unknown function. Then a numerical integration method for the convolution is presented. As an example, an iterative method which can solve the singular nonlinear fractional Riccati equations is proposed. Numerical results are performed to show the efficiency of the method proposed.

1. Introduction

In the recent years, fractional differential equations have been found to be effective to describe some physical phenomena such as rheology, damping laws, fractional random walk, and fluid flow [1–6]; the analytical asymptotic techniques of solutions to various fractional differential equations have been studied and many new numerical techniques have been widely applied to the nonlinear problems. In this paper, we present numerical solutions for the fractional Riccati differential equation

with initial conditions

where is a known function, , , ,  , are arbitrary constants, and is a parameter describing the order of the fractional derivative. The general response expression contains a parameter describing the order of the fractional derivative that can be varied to obtain various responses. For the case of , the fractional equation reduces the classical Riccati differential equation.

In this paper, we will introduce a modified fractional differential operator proposed by Caputo and Mainardi [3]. We mention the important example: the Riccati equations, which play a significant role in many fields of engineering and applied science [1, 5], and the analytical solution for Riccati equation in an explicit form seems to be unlikely except for certain special situations. Therefore, the numerical solutions for the Riccati differential equations are important for engineers and scientists. For example, Mohammadi and Hosseini [5] and Momani and Shawagfeh [7] solved nonlinear partial differential equations of fractional order by the Adomian decomposition method (ADM), which is based on applying the inverse of the fractional operator and has complicated algorithms in calculating much Adomian polynomials for nonlinear problems, and the He’s variational iteration method (VIM) [8], the piecewise variational iteration method (VIM) [9], which does not need the so-called Adomian polynomials but depends on variational theory with the proper selection of the initial approximation , and they both need a large number of iterative calculations. The homotopy perturbation (HPM) [10, 11] provides an effective procedure for explicit and numerical solutions of differential equations; however, the rate of convergence of the series solution depends on the auxiliary parameter and an initial guess .

Despite the progresses outlined above, the literatures on high accuracy and easy to implement numerical techniques are suitable for solving nonlinear Riccati differential equations. The wavelet method, which is a general analytical method and has been widely used by many mathematicians and engineers to solve various functional equations, was applied to signal decompositions and reconstructions, Laplace inversions [12], differential equation solutions [13–15], and active vibration control of piezoelectric smart structures [16, 17]. However, it is somewhat surprising that only very few studies have focused on the solution of fractional differential equations by using wavelet methods [18, 19].

In this paper, we introduce a Coiflets-based wavelet Laplace method (CWLM) that can efficiently solve Riccati differential equation. This method depends on an explicit wavelet approximation scheme for the nonlinear terms of unknown function in the equation, in which series coefficients are just the function samplings at corresponding nodal points, and also by using Laplace transform, the equation with singular integral kernel was converted into equivalent nonsingular integral equation. At last, numerical simulations are performed to show the efficiency of the method proposed.

2. Preliminaries

2.1. Fractional Derivative Operator

Here we give some necessary definitions and mathematical preliminaries of the fractional calculus which are used in this paper. The two most commonly used definitions are the Riemann-Liouville and Caputo. The difference between the two definitions is in the order of evaluation. The Riemann-Liouville fractional derivative operator of order is defined as [2]

Since the Riemann-Liouville fractional derivative has some disadvantages when trying to model real-world phenomena with initial and boundary conditions of fractional differential equations, in this paper, we will introduce a modified fractional differential operator which is proposed by Caputo and Mainardi [3].

For to be the smallest integer that exceeds , the Caputo time-fractional derivative operator of order is defined as [3]

where is -order derivative of the function and is the Gamma function.

The Laplace transform of Caputo fractional derivative operator is [3]

For more details on fractional calculus see [2, 3].

2.2. Coiflets Orthogonality Wavelet

Function is called orthogonal scaling function, if has the following properties.

(1)There is a sequence set , and it satisfies

(2)Set . Satisfy the orthogonality conditions; that is,

 where

 is Kronecker-Delta function.(3)For all , note

 where

 Function space is

Satisfy the following relations:

(i);(ii);(iii);(iv), is dense in ;(v);(vi)setting forming a Riese base of ; that is, there exist constants and , satisfy

For the square sum of sequence space , arbitrary sequence . Where coefficient is low-pass filter coefficient, it usually only has a finite number of nonzero value.

For any function , we have the following approximation:

where the coefficient

For the Coiflets wavelet scaling function, Zhou and Wang [20] have proposed the generalized Gaussian integral method of wavelet, then we have

where is the first-order moment of scaling function, which can be obtained accurately by filter coefficients in the two-scale equations. The approximation accuracy of (13) depends on the corresponding wavelet function vanishing moments, when the corresponding scaling function has the order of vanishing moment , and , and we have [20]

It means that the error of (13) becomes smaller while the scale and the maximum order of vanishing moments of scaling function approximation increase. Substitute (15) into (13)

The form equation (17) of the function of single-point reconstruction formula has the following characteristics [20].

(1)For the composite function of the function or the functional operations , its reconstruction (or approximation) formula is

 where ; the establishment of (18) is obvious, because (18) is for the establishment of an arbitrary function form, of course, including its complex functional form. Therefore, simply square integrable complex functions can be treated as in (18). In this way, furthermore, the importance of (18) is, for an arbitrary function over the function operator (which may be nonlinear operator) transformation, the computing of its right-hand side is equivalent to a linear operator computing. That is, the coefficient of the scaling function can be obtained by put operator on the role of the coefficients .(2)For the integral operator or the derivative operator on the function , we have

 That is, we just need to put the role of operator on in (17).

As we know, wavelet series approximation is a square integrable function defined in the infinite interval, when the approximating function is defined only in a finite interval, and we need to truncate the wavelet series, which may introduce the boundary effect significantly, and the corresponding numerical method of calculation led to decreased accuracy. Traditionally, general treatment of boundary conditions is by using the zero-extension, symmetric, or periodic extension and so on. To some extent, these approaches can effectively inhibit the jitter of the border when it is a special form of approximation function, but not universal, and does not consider the wavelet expansion to meet the boundary conditions. Different from the past expansion of function, in this paper, based on Taylor series expansion of the boundary extension which is applied on the function defined on a finite interval [14, 15], then the function needs to satisfy the corresponding boundary conditions which is embedded in the corresponding Taylor series.

First, we assume that the function , using Taylor series expansion on the borders, forms as follows:

where and ,  , , are expressed by the numerical difference from the interior point, and ,   are the numerical difference coefficients. When , , according to four-point-Malkoff numerical difference formulas [14, 15], we know

where and . For the specific boundary conditions, once the boundary conditions are given in (21), certain elements of the matrix values will be adjusted. For example, for boundary conditions , , we just need to make , , and other take the values given in (21) without change.

Then, (18) can be expressed as

where , , and .

Consider Coiflets scaling function with the support set , , applying the scaling function series approximation equation (17), and is an integral; for the function , we have the following form:

Substituting (22) into (23), then (23) can be rewritten as

where , and we denote

Thus, when specific boundary conditions are given, the differential coefficients of extension can be determined in accordance with the above process, and the corresponding improved scaling function is given by (25).

3. Numerical Algorithm

In this section, we will consider the modified wavelet approximation scheme to solve nonlinear Riccati differential equation (1) with the initial condition expressed by (2).

First, applying the Laplace transform on the time variable in (1) leads to

that is,

in which . or , and the formula of the Laplace transform of the fractional order derivative is considered [3]. When the function obtained by the inverse Laplace transform of is singular, we take , else we take . For example, for , , we just need to take ; then , whose inverse Laplace transform is , and it is nonsingular. Taking the inverse Laplace transform of (27), we obtain

where

By [12], , can be obtained by numerical method. From [6], we know

where is the generalized Mittag-Leffler type functions defined in the power series

It can be seen from (30) that the integral kernel is a nonsingular smooth function with property of when ; applying (24), (25), we have

By using (19), we denote

Applying (24) to approximate the term with in (28), we have

where . Integrating both sides of (34), we obtain

where can be exactly obtained according to Wang [22]. Inserting (35) into (28) and then setting yield

where , () and . Then one can solve the algebraic equations to obtain the solution of the Riccati equation. According to (30), note that if there should be for , which implies that when , thus (36) can be further simplified into the form

It can be seen from (37) that the solution can be directly obtained step-by-step as the index increases. In this process, no matrix inversion is needed. The scale is bigger, and the precision of the results of (36) and (37) is better, which will be shown in examples. In fact, in most practice problems, (37) is possible, and even without the special nature of in some case, we can still solve the nonlinear algebraic equation (36), which is nonsingular and easy. Then we promote the method to the fractional Riccati equations.

4. Numerical Examples

In this section we will give two numerical experimentsto illustrate the efficiency and apply the approach proposed in this paper.

Example 1. Consider the following fractional Riccati equation given in [7, 19, 21]:

with initial conditions

and the exact solution of (38) is when .

Figure 1 is the comparison between the exact result and the numerical result for . We can see that the numerical solution is in very good agreement with the exact solution. Figure 2 shows the absolute errors of the exact solution and the results obtained by CWLM with different . One can see that approximate solution converges to the exact solution, and the absolute error goes down while is increased. Contrast of numerical solution of (38) for and is given in Figure 3.

Figure 1: Comparison between the numerical and exact solution of Example 1 for .

Figure 2: The absolute errors of the exact solution and the results obtained by CWLM with different for .

Figure 3: Solution of (38) obtained by the CWLM for and .

Table 1 is the solutions of Example 1 obtained by different methods when . Compared with results in [21] obtained by using HAM with homotopy-Pade approximations and HWOMM in [19], which need to solve a nonlinear system of nearly hundreds of algebraic equations.

Table 1: Solution of (38) obtained by different methods for .

Methods

HAM [21]
HWOMM [19]
CWLM

1.0
0.6982
0.6987
0.6667
0.6837
0.6883

2.0
0.7858
0.7857
0.7732
0.7796
0.7814

3.0
0.8268
0.8258
0.8190
0.8224
0.8234

4.0
0.8525
0.8499
0.8455
0.8477
0.8484

5.0
0.8714
0.8664
0.8633
0.8648
0.8652

6.0
0.8869
0.8785
0.8761
0.8773
0.8776

7.0
0.9004
0.8878
0.8860
0.8869
0.8871

8.0
0.9127
0.8953
0.8938
0.8946
0.8948

9.0
0.9243
0.9016
0.9003
0.9009
0.9011

10.0
0.9354
0.9068
0.9057
0.9062
0.9064

Example 2. Consider the following fractional Riccati equation [7, 19]:

With initial conditions

the exact solution is when .

Figure 4 is the comparison between the exact result and the numerical result for . One can see that CWLM is accurate and is able to solve this nonlinear Riccati differential equation in a very wider region. Figure 5 shows the absolute errors of the exact solution and the results obtained by CWLM with different for Example 2 when , in which it can be seen that the bigger the value of the scale is, the smaller absolute error is. Figure 6 is the solution of (40) obtained by the CWLM for and , compared with results obtained by other methods in [19, 23, 24], which are given in Table 2.

Table 2: Solution of (40) obtained by the EHPM, MHPM, HWOMM, and CWLM for .

Methods

EHPM [23]
MHPM [24]
HWOMM [19]
CWLM

1.0
1.9753
2.0874
1.8186
1.1571
1.1538

2.0
2.2062
0.7787
2.1796
1.8206
1.8198

3.0
2.2813
5.8102
2.2684
2.0566
2.0629

4.0
2.3153
−0.0807
2.3062
2.1539
2.1643

5.0
2.3340
−396.4145
2.3273
2.2038
2.2164

6.0
2.3460

2.3408
2.2337
2.2477

7.0
2.3544

2.3504
2.2537
2.2685

8.0
2.3608

2.3574
2.2679
2.2833

9.0
2.3657

2.3629
2.2787
2.2945

10.0
2.3679

2.3673
2.2870
2.3031

Figure 4: Comparison between the numerical and exact solutions of Example 2 for .

Figure 5: The absolute errors of the exact solution and the results obtained by CWLM with different for .

Figure 6: Solution of (40) obtained by the CWLM for and .

5. Conclusions

In this paper, a numerical method based on the Coiflets wavelet operational method is applied to solve the fractional differential equations. In this method, the equation with fractional differential order is transferred to an integral equation of convolution type by the Laplace transform and then the solution is approximated by the modified wavelet approximate scheme. This simple method was established by Zhou et al. [14] (including the author) and has been applied to solve nonlinear equations of vibration, diffusion, and wave equations. Compared with the results for solving fractional Riccati differential equation by the other numerical methods [7, 19, 21, 23, 24], the results for numerical examples demonstrate that the present method can give a high accurate approximation in a larger region. This is also the advantage of the present method.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by National Natural Science Foundation of China (11302081) and Huazhong Agricultural University Scientific and Technological Self-innovation Foundation under Grant (52902-0900206074).

References

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

I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.

M. Caputo and F. Mainardi, “Linear models of dissipation in anelastic solids,” La Rivista del Nuovo Cimento, vol. 1, no. 2, pp. 161–198, 1971.

R. Metzler, E. Barkai, and J. Klafter, “Anomalous diffusion and relaxation close to thermal equilibrium: a fractional Fokker-Planck equation approach,” Physical Review Letters, vol. 82, no. 18, pp. 3563–3567, 1999.

F. Mohammadi and M. M. Hosseini, “A comparative study of numerical methods for solving quadratic Riccati differential equations,” Journal of the Franklin Institute, vol. 348, no. 2, pp. 156–164, 2011.

D. Hong, Y. Kim, and J. Wang, “A new approach for the analysis solution of dynamic systems containing fractional derivative,” Journal of Mechanical Science and Technology, vol. 20, no. 5, pp. 658–667, 2006.

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. Geng, Y. Lin, and M. Cui, “A piecewise variational iteration method for Riccati differential equations,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2518–2522, 2009.

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.

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.

J. Wang, Y. Zhou, and H. Gao, “Computation of the Laplace inverse transform by application of the wavelet theory,” Communications in Numerical Methods in Engineering, vol. 19, no. 12, pp. 959–975, 2003.

Y. 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.

Y. H. Zhou, X. M. Wang, J. Z. Wang, and X. J. Liu, “A wavelet numerical method for solving nonlinear fractional vibration, diffusion and wave equations,” CMES: Computer Modeling in Engineering and Sciences, vol. 77, no. 2, pp. 137–160, 2011.

X. J. Liu, J. Z. Wang, X. M. Wang, and Y. H. Zhou, “Exact solutions of multi-term fractional diffusion-wave equations with Robin type boundary conditions,” Applied Mathematics and Mechanics, vol. 35, no. 1, pp. 49–62, 2014.

Y. Zhou, J. Wang, X. J. Zheng, and Q. Jiang, “Vibration control of variable thickness plates with piezoelectric sensors and actuators based on wavelet theory,” Journal of Sound and Vibration, vol. 237, no. 3, pp. 395–410, 2000.

J. Wang, X. Wang, and Y. Zhou, “A wavelet approach for active-passive vibration control of laminated plates,” Acta Mechanica Sinica, vol. 28, no. 2, pp. 520–531, 2012.

Ü. Lepik, “Solving fractional integral equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 214, no. 2, pp. 468–478, 2009.

Y. L. Li, N. Sun, B. C. Zheng, Q. Wang, and Y. C. Zhang, “Wavelet operational matrix method for solving the Riccati differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 3, pp. 483–493, 2014.

Y. H. Zhou and J. Z. Wang, “A generalized Gaussian integral method for the calculation of scaling function transforms of wavelets and its applications,” Acta Mathematica Scientia A: Shuxue Wuli Xuebao, vol. 19, no. 3, pp. 293–300, 1999.

J. Cang, Y. Tan, H. Xu, and S. Liao, “Series solutions of non-linear Riccati differential equations with fractional order,” Chaos, Solitons and Fractals, vol. 40, no. 1, pp. 1–9, 2009.

J. Z. Wang, Generalized theory and arithmetic of orthogonal wavelets and applications to researches of mechanics including piezoelectric smart structures [Ph.D. dissertation], Lanzhou University, 2001.

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.

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.