Iterative Methods for Nonlinear Equations or Systems and Their Applications 2014
View this Special IssueResearch Article  Open Access
Xiaomin Wang, "A CoifletsBased Wavelet Laplace Method for Solving the Riccati Differential Equations", Journal of Applied Mathematics, vol. 2014, Article ID 257049, 8 pages, 2014. https://doi.org/10.1155/2014/257049
A CoifletsBased Wavelet Laplace Method for Solving the Riccati Differential Equations
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 socalled 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 Coifletsbased 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 RiemannLiouville and Caputo. The difference between the two definitions is in the order of evaluation. The RiemannLiouville fractional derivative operator of order is defined as [2] Since the RiemannLiouville fractional derivative has some disadvantages when trying to model realworld 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 timefractional 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 KroneckerDelta 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 lowpass 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 firstorder moment of scaling function, which can be obtained accurately by filter coefficients in the twoscale 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 singlepoint 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 righthand 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 zeroextension, 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 fourpointMalkoff 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 MittagLeffler 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 stepbystep 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.
Table 1 is the solutions of Example 1 obtained by different methods when . Compared with results in [21] obtained by using HAM with homotopyPade approximations and HWOMM in [19], which need to solve a nonlinear system of nearly hundreds of algebraic equations.

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.

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 Selfinnovation Foundation under Grant (529020900206074).
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. View at: MathSciNet
 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. View at: Publisher Site  Google Scholar
 R. Metzler, E. Barkai, and J. Klafter, “Anomalous diffusion and relaxation close to thermal equilibrium: a fractional FokkerPlanck equation approach,” Physical Review Letters, vol. 82, no. 18, pp. 3563–3567, 1999. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar
 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 Site  Google Scholar  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 Site  Google Scholar  MathSciNet
 F. Geng, Y. Lin, and M. Cui, “A piecewise variational iteration method for Riccati differential equations,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2518–2522, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Google Scholar
 X. J. Liu, J. Z. Wang, X. M. Wang, and Y. H. Zhou, “Exact solutions of multiterm fractional diffusionwave equations with Robin type boundary conditions,” Applied Mathematics and Mechanics, vol. 35, no. 1, pp. 49–62, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar
 J. Wang, X. Wang, and Y. Zhou, “A wavelet approach for activepassive vibration control of laminated plates,” Acta Mechanica Sinica, vol. 28, no. 2, pp. 520–531, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Ü. Lepik, “Solving fractional integral equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 214, no. 2, pp. 468–478, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Google Scholar  MathSciNet
 J. Cang, Y. Tan, H. Xu, and S. Liao, “Series solutions of nonlinear Riccati differential equations with fractional order,” Chaos, Solitons and Fractals, vol. 40, no. 1, pp. 1–9, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
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.