`International Journal of Mathematics and Mathematical SciencesVolume 2014 (2014), Article ID 304745, 10 pageshttp://dx.doi.org/10.1155/2014/304745`
Research Article

## Legendre Wavelet Operational Matrix Method for Solution of Riccati Differential Equation

Department of Mathematics, SASTRA University, Thanjavur 613 401, India

Received 12 February 2014; Revised 20 May 2014; Accepted 5 June 2014; Published 24 June 2014

Copyright © 2014 S. Balaji. 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 Legendre wavelet operational matrix method (LWM) is presented for the solution of nonlinear fractional-order Riccati differential equations, having variety of applications in quantum chemistry and quantum mechanics. The fractional-order Riccati differential equations converted into a system of algebraic equations using Legendre wavelet operational matrix. Solutions given by the proposed scheme are more accurate and reliable and they are compared with recently developed numerical, analytical, and stochastic approaches. Comparison shows that the proposed LWM approach has a greater performance and less computational effort for getting accurate solutions. Further existence and uniqueness of the proposed problem are given and moreover the condition of convergence is verified.

#### 1. Introduction

In recent years, use of fractional-order derivative goes very strongly in engineering and life sciences and also in other areas of science. One of the best advantages of use of fractional differential equation is modeling and control of many dynamic systems. Fractional-order derivatives are used in fruitful way to model many remarkable developments in those areas of science such as quantum chemistry, quantum mechanics, damping laws, rheology, and diffusion processes [15] described through the models of fractional differential equations (FDEs). Modeling of a physical phenomenon depends on two parameters such as the time instant and the prior time history; because of this reason, reasonable modeling through fractional calculus was successfully achieved. The abovementioned advantages and applications of FDEs attracted researchers to develop efficient methods to solve FDEs in order to get accurate solutions to such problems and more active research is still going on in those areas. Most of the FDEs are complicated in their structure; hence finding exact solutions for them cannot be simple. Therefore, one can approach the best accurate solution of FDEs through analytical and numerical methods. Designing accurate or best solution to FDEs, many methods are developed in recent years; each method has its own advantages and limitations. This paper aims to solve a FDE called fractional-order Riccati differential equation, one of the important equations in the family of FDEs. The Riccati equations play an important role in engineering and applied science [6], especially in quantum mechanics [7] and quantum chemistry [8, 9]. Therefore, solutions to the Riccati differential equations are important to scientists and engineers. Solving fractional-order Riccati differential equation, the most significant methods are Adomian decomposition method [10], homotopy perturbation method [1114], homotopy analysis method [15, 16], Taylor matrix method [17], Haar wavelet method [18], and combination of Laplace, Adomian decomposition, and Padé approximation [19] methods.

Several numerical methods for approximating the solution of nonlinear fractional-order Riccati differential equations are known. Raja et al. [20] developed a stochastic technique based on particle swarm optimization and simulated annealing. They were used as a tool for rapid global search method and simulated annealing for efficient local search method. A fractional variational iteration method described in the Riemann-Liouville derivative has been applied in [21], to give an analytical approximate solution to nonlinear fractional Riccati differential equation. A combination of finite difference method and Padé-variational iteration numerical scheme was proposed by Sweilam et al. [22]. Moreover, an analytical scheme comprising the Laplace transform, the Adomian decomposition method (ADM), and the Padé approximation is given in [19].

However, the abovementioned methods have some restrictions and disadvantages in their performance. For example, very complicated and toughest Adomian polynomials are constructed in the Adomian decomposition method. In the variational iteration method identification Lagrange multiplier yields an underlying accuracy. The homotopy perturbation method needs a linear functional equation in each iteration to solve nonlinear equations; forming these functional equations is very difficult. The performance of the homotopy analysis method very much depends on the choice of the auxiliary parameter of the zero-order deformation equation. Moreover, the convergence region and implementation of these results are very small.

In recent years, wavelets theory is one of the growing and predominantly new methods in the area of mathematical and engineering research. It has been applied in vast range of engineering sciences; particularly, they are used very successfully for waveform representation and segmentations in signal analysis and time-frequency analysis and in the mathematical sciences it is used in thriving manner for solving variety of linear and nonlinear differential and partial differential equations and fast algorithms for easy implementation [23]. Moreover, wavelets build a connection with fast numerical algorithms [24]; this is due to the fact that wavelets admit the exact representation of a variety of functions and operators. The application of Legendre wavelet and its operational matrix for solving differential, integral, and fractional-order differential equations is thoroughly considered in [25, 26].

In this work, the nonlinear Riccati differential equations of fractional-order are approached analytically by using Legendre wavelets method. The operational matrix of Legendre wavelet is generalized for fractional calculus in order to solve fractional and classical Riccati differential equations. The Legendre wavelet method (LWM) is illustrated by application, and obtained results are compared with recently proposed method for the fractional-order Riccati differential equation. We have adopted Legendre wavelet method to solve Riccati differential equations not only due to its emerging application but also due to its greater convergence region.

The rest of the paper is as follows. In Section 2, definitions and mathematical preliminaries of fractional calculus are presented. In Section 3, Legendre wavelet, its properties, function approximations, and generalized Legendre wavelet operational matrix fractional calculus are discussed. Section 4 establishes application of proposed method in the solution of Riccati differential equations, existence and uniqueness solution of the proposed problem, and convergence analyses of the proposed approach. Section 5 deals with the illustrative examples and their solutions by the proposed approach. Section 6 ends with our conclusion.

#### 2. Preliminaries and Notations

The notations, definitions, and preliminary facts present in this section will be used in forthcoming sections of this work. Several definitions of fractional integrals and derivatives have been proposed after the logical definition given by Liouville. Important and few of these definitions include the Riemann-Liouville, the Caputo, the Weyl, the Hadamard, the Marchaud, the Riesz, the Grunwald-Letnikov, and the Erdelyi-Kober. As stated in [25], the Caputo fractional derivative uses initial and boundary conditions of integer order derivatives having some physical interpretations. Because of this specific reason, in this work, we will use the Caputo fractional derivative proposed by Caputo [27] in the theory of viscoelasticity.

The Caputo fractional derivative of order and is continuous and is defined by where is the Riemann-Liouville fractional integral operator of order and is the gamma function.

The fractional integral of  is given as Properties of fractional integrals and derivatives are as follows [28], for .

The fractional-order integral satisfies the semigroup property The integer order derivative and fractional-order derivative commute with each other: The fractional integral operator and fractional derivative operator do not satisfy the commutative property. In general, But in the reverse way we have

#### 3. Generalized Legendre Wavelet Operational Matrix to Fractional Integration

##### 3.1. Legendre Wavelets

A family of functions were constituted by wavelets and constructed from dilation and translation of a single function called mother wavelet. When the parameters of dilation and of translation vary continuously, following are the family of continuous wavelets [29]: If the parameters and are restricted to discrete values as ,,,, and and are positive integers, following are the family of discrete wavelets: where form a wavelet basis for . In particular, when and form an orthonormal basis [29].

Legendre wavelets have four arguments; can assume any positive integer, is the order for Legendre polynomials, and is the normalized time. They are defined on the interval [0,1) as [30, 31] where and . The coefficient is for orthonormality, the dilation parameter is , and translation parameter is . are the well-known Legendre polynomials of order m defined on the interval and can be determined with the aid of the following recurrence formulae: The Legendre wavelet series representation of the function defined over is given by where , in which denotes the inner product. If the infinite series in (12) is truncated, then (12) can be written as where  and   are  matrices given by Taking suitable collocation points as we defined the -square Legendre matrix where  ; correspondingly, we have The Legendre matrix is an invertible matrix, and the coefficient vector is obtained by .

##### 3.2. Operational Matrix of the Fractional Integration

The integration of the defined in (14) can be approximated by Legendre wavelet series with Legendre wavelet coefficient matrix : where the -square matrix is called Legendre wavelet operational matrix of integration.

The -set of block-pulse functions is defined on [0, l) as follows: where .

The functions are disjoint and orthogonal. That is, The orthogonality property of block-pulse function is obtained from the disjointness property.

An arbitrary function can be expanded into block-pulse functions as where are the coefficients of the block-pulse function, given by The Legendre wavelets can be expanded into -set of block-pulse functions as where .

The fractional integral of block-pulse function vector can be written as where is given in [32].

Now, we introduce the derivation process of the Legendre wavelet operational matrix of the fractional integration: where the -square matrix is called Legendre wavelet operational matrix of the fractional integration.

Using (23) and (25), we have From (25) and (26), we get and by (23), (27) becomes Then, the Legendre wavelet operational matrix  of fractional integration is given by Following is the Legendre wavelet operational matrix  of fractional-order integration, for the particular values of = 2, = 3, and = 0.5:

#### 4. Application to Fractional Riccati Differential Equation

In this section, we use the generalized Legendre wavelet operational matrix to solve nonlinear Riccati differential equation and we discuss the existence and uniqueness of solutions with initial conditions and convergence criteria of the proposed LWM approach.

Consider the fractional-order Riccati differential equation of the form subject to the initial condition Let us suppose that the functions , , , and are approximated using Legendre wavelet as follows: where , , , , and are given in (14).

Using (6), we can write By (25) and (32), (34) leads to where .

Substituting (33) and (35) into (31), we have Substituting (23) into (36), we have where , , , and are known. Equation (37) represents a system of nonlinear equations with unknown vector . This system of nonlinear equations can be solved by Newton method for the unknown vector and we can get the approximation solution by including into (35).

##### 4.1. Existence and Uniqueness of Solutions

Consider the fractional-order Riccati differential equation of the forms (31) and (32). The nonlinear term in (31) is and , are known functions. For , the fractional-order Riccati converts into the classical Riccati differential equation.

Definition 1. Let , , and be the class of all continuous functions defined on , with the norm which is equivalent to the sup norm of . That is, .

Remark. Assume that solution of fractional-order Riccati differential equations (31) and (32) belongs to the space , in order to study the existence and uniqueness of the initial value problem.

Definition 2. The space of integrable functions in the interval is defined as

Theorem 3. The initial value problem given by (31) and (32) has a unique solution:

Proof. By (1), the fractional differential equation (31) can be written as and becomes Now we define the operator and then hence, we have which implies that the operator given by (43) has a unique fixed point and consequently the given integral equation has a unique solution . Also we can see that Now, from (42), we have from which we can deduce that and .
Now, again from (42), (43), and (46), we get which implies that the integral equation (46) is equivalent to the initial value problem (32) and the theorem is proved.

##### 4.2. Convergence Analyses

Let where form a wavelet basis for . In particular, when form an orthonormal basis [29].

By (14), let be the solution of (31) where in which denotes the inner product: Let where .

Let be a sequence of partial sums. Then, Further, As , from Bessel’s inequality, we have is convergent.

It implies that is a Cauchy sequence and it converges to (say).

Also, which is possible only if . That is, both and converge to the same value, which indeed give the guarantee of convergence of LWM.

#### 5. Numerical Examples

In order to show the effectiveness of the Legendre wavelets method (LWM), we implement LWM to the nonlinear fractional Riccati differential equations. All the numerical experiments were carried out on a personal computer with some MATLAB codes. The specifications of PC are Intel core i5 processor and with Turbo boost up to 3.1 GHz and 4 GB of DDR3 memory. The following problems of nonlinear Riccati differential equations are solved with real coefficients.

Example 1. Consider the following nonlinear fractional Riccati differential equation: with initial condition Exact solution for was found to be The integral representation of (54) and (55) is given by Let and then By substituting (59) and (60) into (58), we get the following system of algebraic equations: By solving the above system of linear equations, we can find the vector . Numerical results are obtained for different values of , , and . Solution obtained by the proposed LWM approach for = 1, = 1, and = 3 is given in Figure 1 and for different values of = 0.6, 0.7, 0.8, and 0.9 and for = 2 and = 5 is graphically given in Figure 2. It can be seen from Figure 1 that the solution obtained by the proposed LWM approach is more close to the exact solution. Table 1 describes the efficiency of the proposed method by comparing with the methods in [20, 22] through their absolute error. The following is used for the errors of the approximation ; that is, . Table 1 shows that very high accuracies are obtained for = 3 and = 5 by the present method.

Table 1: Numerical results of Example 1 for .
Figure 1: Numerical results of Example 1 by LWM for .
Figure 2: Numerical results of Example 1 by LWM for different values of .

Example 2. Consider another fractional-order Riccati differential equation: with initial condition Exact solution for the above equation was found to be The integral representation of (62) and (63) is given by Let and then By substituting (66) and (67) in (62), we get the following system of algebraic equations: By solving the above system of linear equations, we can find the vector . Numerical results are obtained for different values of , , and . Results obtained by LWM for = 1, = 2, and = 3 are shown in Figure 3 and it can be seen from the figure that solution given by the LWM merely coincides with the exact solution. Figure 4 shows the obtained results of (62) and (63) by LWM for different values of and for = 2 and = 5. Table 2 describes the efficiency of the proposed method by comparing with the methods in [20, 22] through their absolute error. Table 1 shows that very high accuracies are obtained for = 3 and = 5 by the present method and from these results we can identify that guarantee of convergence of the proposed LWM approach is very high.

Table 2: Numerical results of Example 2 for .
Figure 3: Numerical results of Example 2 by LWM for .
Figure 4: Numerical results of Example 2 by LWM for different values of .

Example 3. Let us consider another problem of nonlinear Riccati differential equation: with initial condition When = 1, its exact solution is given by where is the Bessel function of first kind.
The integral representation of (69) and (70) is given by Let and then By substituting (73) and (74) into (69), we get the following system of algebraic equations: By solving the above system of linear equations, we can find the vector . Numerical results are obtained for different values of ,   , and . Obtained results for (69) and (70) are shown in Figures 5 and 6. Figure 5 shows the solutions obtained by LWM for different values of and for = 2 and = 4. Figure 6 compares the solution obtained by LWM with the exact solution of (69) and (70) when = 1, = 1, and = 2. So far, there are no published results of absolute error for this problem and hence we are unable to compare absolute error of our method with the existing methods. From these results, we can see that the proposed LWM approach gives the solution which is very close to the exact solution and outperformed recently developed approaches for the nonlinear fractional-order Riccati differential equations in terms of solution quality and convergence criteria.

Figure 5: Numerical results of Example 3 by LWM for different values of .
Figure 6: Numerical results of Example 3 by LWM for .

#### 6. Conclusions

Nonlinear fractional-order Riccati differential equations play an important role in the modeling of many biological, physical, chemical, and real life problems. Therefore, it is necessary to develop a method which would give more accurate solutions to such type of problems with greater convergence criteria. In this work, a Legendre’s wavelet operational matrix method called LWM was proposed for solving nonlinear fractional-order Riccati differential equations. Comparison was made for the solutions obtained by the proposed method and with the other recent approaches developed for the same problem through their error analysis; obtained results show that the proposed LWM yields more accurate and reliable solutions even for small values of and , which assures the best approximate solution in less computational effort. Further, we have discussed the convergence criteria of proposed scheme, which indeed provides the guarantee of consistency and stability of the proposed LWM scheme for the solutions of nonlinear fractional Riccati differential equations.

#### Conflict of Interests

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

#### References

1. N. A. Khan, M. Jamil, A. Ara, and S. Das, “Explicit solution for time-fractional batch reactor system,” International Journal of Chemical Reactor Engineering, vol. 9, article A91, 2011.
2. V. Feliu-Batlle, R. R. Pérez, and L. S. Rodríguez, “Fractional robust control of main irrigation canals with variable dynamic parameters,” Control Engineering Practice, vol. 15, no. 6, pp. 673–686, 2007.
3. I. Podlubny, “Fractional-order systems and $PI\sp\lambda D\sp\mu$-controllers,” IEEE Transactions on Automatic Control, vol. 44, no. 1, pp. 208–214, 1999.
4. R. Garrappa, “On some explicit Adams multistep methods for fractional differential equations,” Journal of Computational and Applied Mathematics, vol. 229, no. 2, pp. 392–399, 2009.
5. M. Jamil and N. A. Khan, “Slip effects on fractional viscoelastic fluids,” International Journal of Differential Equations, vol. 2011, Article ID 193813, 19 pages, 2011.
6. 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.
7. R. Shankar, Principles of Quantum Mechanics, Plenum Press, New York, NY, USA, 1980.
8. S. Fraga, M. J. García de la Vega, and E. S. Fraga, The Schrödinger and Riccati Equations, vol. 70 of Lect. Notes. Chem., 1999.
9. L. B. Burrows and M. Cohen, “Schrodinger’s wave equation-A lie algebra treatment,” in Fundamental World of Quantum Chemistry: A Tribute to the Memory of Per-Olov Löwdin, E. J. Brändas and E. S. Kryachko, Eds., Kluwer, Dordrecht, The Netherlands, 2004.
10. 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.
11. 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.
12. N. A. Khan, A. Ara, and M. Jamil, “An efficient approach for solving the Riccati equation with fractional orders,” Computers & Mathematics with Applications, vol. 61, no. 9, pp. 2683–2689, 2011.
13. 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.
14. 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.
15. 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.
16. 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.
17. M. Gülsu and M. Sezer, “On the solution of the Riccati equation by the Taylor matrix method,” Applied Mathematics and Computation, vol. 176, no. 2, pp. 414–421, 2006.
18. Y. Li and L. Hu, “Solving fractional Riccati differential equations using Haar wavelet,” in Proceedings of the 3rd International Conference on Information and Computing (ICIC '10), pp. 314–317, Wuxi, China, June 2010.
19. N. A. Khan and A. Ara, “Fractional-order Riccati differential equation: analytical approximation and numerical results,” Advances in Difference Equations, vol. 2013, article 185, 2013.
20. M. A. Z. Raja, J. A. Khan, and I. M. Qureshi, “A new stochastic approach for solution of Riccati differential equation of fractional order,” Annals of Mathematics and Artificial Intelligence, vol. 60, no. 3-4, pp. 229–250, 2010.
21. M. Merdan, “On the solutions fractional Riccati differential equation with modified RIEmann-Liouville derivative,” International Journal of Differential Equations, vol. 2012, Article ID 346089, 17 pages, 2012.
22. N. H. Sweilam, M. M. Khader, and A. M. S. Mahdy, “Numerical studies for solving fractional Riccati differential equation,” Applications and Applied Mathematics, vol. 7, no. 2, pp. 595–608, 2012.
23. C. K. Chui, Wavelets: A Mathematical Tool for Signal Processing, SIAM, Philadelphia, Pa, USA, 1997.
24. G. Beylkin, R. Coifman, and V. Rokhlin, “Fast wavelet transforms and numerical algorithms, I,” Communications on Pure and Applied Mathematics, vol. 44, no. 2, pp. 141–183, 1991.
25. M. ur Rehman and R. Ali Khan, “The Legendre wavelet method for solving fractional differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4163–4173, 2011.
26. S. Balaji, “A new approach for solving Duffing equations involving both integral and non-integral forcing terms,” Ain Shams Engineering Journal, 2014.
27. M. Caputo, “Linear models of dissipation whose Q is almost frequency independent II,” Geophysical Journal of the Royal Astronomical Society, vol. 13, pp. 529–539, 1967.
28. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of North-Holland Mathematics Studies, Elsevier, 2006.
29. J. S. Gu and W. S. Jiang, “The Haar wavelets operational matrix of integration,” International Journal of Systems Science, vol. 27, no. 7, pp. 623–628, 1996.
30. M. Razzaghi and S. Yousefi, “Legendre wavelets direct method for variational problems,” Mathematics and Computers in Simulation, vol. 53, no. 3, pp. 185–192, 2000.
31. M. Razzaghi and S. Yousefi, “Legendre wavelets method for constrained optimal control problems,” Mathematical Methods in the Applied Sciences, vol. 25, no. 7, pp. 529–539, 2002.
32. A. Kilicman and Z. A. A. Al Zhour, “Kronecker operational matrices for fractional calculus and some applications,” Applied Mathematics and Computation, vol. 187, no. 1, pp. 250–265, 2007.