We obtain the approximate analytical solution for the fractional quadratic Riccati differential equation with the Riemann-Liouville derivative by using the Bernstein polynomials (BPs) operational matrices. In this method, we use the operational matrix for fractional integration in the Riemann-Liouville sense. Then by using this matrix and operational matrix of product, we reduce the problem to a system of algebraic equations that can be solved easily. The efficiency and accuracy of the proposed method are illustrated by several examples.
1. Introduction
The Riccati differential equation is named after the Italian Nobleman Count Jacopo Francesco Riccati (1676–1754). The book of Reid [1] contains the fundamental theories of the Riccati equation, with applications to random processes, optimal control, and diffusion problems. Moreover, it is well known that the one-dimensional static Schrödinger equation is closely related to a Riccati differential equation [2]. Solitary wave solution of a nonlinear partial differential equation can be represented as a polynomial in two elementary functions satisfying a projective Riccati equation [3].
In this paper, we are dealing with the fractional quadratic Riccati differential equation as follows:
(1)Dtαay(t)=f(t)+b(t)y(t)+a(t)y2(t),a<t≤b,
subject to the initial conditions
(2)y(k)(a)=yak,k=0,1,…,⌈α⌉-1,
where yak(k = 0,1,…,⌈α⌉-1) are constants and a(t), b(t), and f(t) are known functions. We can see [4] to guarantee the existence and uniqueness of the solution of (1) with initial conditions (2).
The general response expression (1) contains a parameter α, the order of the fractional derivative that can be varied to obtain various responses. In the case that α is integer, then (1) is reduced to the classical Riccati differential equation.
This problem to develop the analytical and numerical method to solve the Riccati differential equation with standard derivative, the Caputo fractional derivative and the Riemann-Liouville fractional derivative, has attracted much attention and has been studied by many authors [5–18].
The aim of this work is using the Bernstein polynomials for solving the problem (1) and (2). We notice that the problem presented [18] was in the Caputo sense but in our work, the problem is with the Riemann-Liouville derivative; therefore we considered a more general space of functions. Also, in [18], the authors used the polynomials in the form of Bi,mα(x) = (mi)xαi(1-xα)m-i(i = 0,1,…,m) that is different from the standard Bernstein polynomials. So, the operational matrices in this work are different from those in [18].
The organization of this paper is as follows. In Section 2, the Bernstein polynomials are introduced. Some basic definitions and properties of the fractional calculus and also the BPs operational matrix for the Riemann-Liouville fractional integration are presented in Section 3. In Section 4, by BPs operational matrices, we solve the fractional quadratic Riccati differential equation. In Section 5, we discuss the convergence of the proposed method. In Section 6, several examples are considered to evaluate the power and effectiveness of the presented method. Some conclusions are summarized in the last section.
2. The Bernstein Polynomials and Their Properties
On the interval [0,1] we define the Bernstein polynomials (BPs) of mth degree as follows [19]:
(3)Bi,m(x)=(mi)xi(1-x)m-i,i=0,1,…,m.
Set {B0,m(x),B1,m(x),…,Bm,m(x)} in the Hilbert space L2[0,1] is a complete basis. We can write Φm(x)=ATm(x), where A is a matrix upper triangular, Tm(x)=[1,x,…,xm]T, and Φm(x)=[B0,m(x),B1,m(x),…,Bm,m(x)]T [20].
As a result, any polynomial of degree m can be expanded in terms of linear combination of Bi,m(x)(i=0,1,…,m) as given below:
(4)P(x)=∑i=0mciBi,m(x)=cTΦm(x),
where
(5)c=(∫01Φ(x)Φ(x)Tdx)-1(∫01P(x)Φ(x)dx).
The approximation of functions within the Bernstein polynomials and convergence analysis can be found in [20, 21].
3. BPs Operational Matrix for the Riemann-Liouville Fractional Integration
In this section, firstly, we give some basic definitions and properties of the fractional calculus which are used further in this paper.
Definition 1 (see [4, 22–24]).
Let α≥0; the operator Itαa, defined on L1[a,b] by
(6)Itαaf(t)=1Γ(α)∫at(t-x)α-1f(x)dx,a≤t≤b,It0af(t)=f(t),
is called the Riemann-Liouville fractional integral operator of order α.
Definition 2 (see [4, 22–24]).
Let α≥0 and n=⌈α⌉; the operator Dtαa, defined by
(7)Dtαaf(t)=1Γ(n-α)(ddt)n∫at(t-x)n-α-1f(x)dx,a≤t≤b,Dt0af(t)=f(t),
is called the Riemann-Liouville fractional derivative operator of order α.
Definition 3 (see [4, 22–24]).
Let α≥0, n=⌈α⌉, and dnf(x)/dxn∈L1[a,b]. The operator Dtαac, defined by
(8)Dtαacf(t)=1Γ(n-α)∫at(t-x)n-α-1dnf(x)dxndx,a≤t≤b,Dt0acf(t)=f(t),
is called the Caputo fractional derivative operator of order α.
Lemma 4.
If α≥0, n=⌈α⌉, and a≤t≤b, then
(9)(1)DtαacItαaf(t)=f(t),(10)(2)ItαaDtaacf(t)=f(t)-∑k=0n-1f(k)(a)k!(x-a)k,(11)(3)Dtαacf(t)=Dtαaf(t)-∑k=0n-1f(k)(a)Γ(k-α+1)(x-a)k-α.
Proof.
See [22–24].
Theorem 5.
One can get BPs operational matrix Fα from order (m+1)×(m+1) for the Riemann-Liouville fractional integral as
(12)ItαaΦm(t)=1Γ(α)∫at(t-x)α-1Φm(x)dx≈FαΦm(t).
Proof.
See [21].
4. BPs for Solving the Fractional Quadratic Riccati Differential Equation
Firstly, we use the initial conditions to reduce a given initial-value problem to a problem with zero initial conditions.
So, we define
(13)y(t)=y^(t)+z(t),
where y^(t) is some known function that satisfied the initial conditions (2) and z(t) is a new unknown function.
Substituting (13) in (1) and (2), we have an initial-value problem as follows:
(14)Dtαaz(t)=f~(t)+b~(t)z(t)+a~(t)z2(t),a<t≤b,
subject to the initial conditions
(15)z(k)(a)=0,k=0,1,…,⌈α⌉-1.
On the other hand, by (11) in Lemma 4 we can write
(16)Dtαacz(t)=Dtαaz(t).
Also, by using Lemma 3.3 in [20] the inputs f~(t), a~(t), b~(t), and Dtαaz(t) can be approximated as follows:
(17)f~(t)≈FTΦm(t),(18)Dtαaz(t)≈CTΦm(t),(19)a~(t)≈ATΦm(t),(20)b~(t)≈BTΦm(t),
where F, A, and B are known (m+1) column vectors and C is an unknown (m+1) column vector.
From (10), (15), (16), (18), and (12), we have
(21)z(t)=(10)ItαaDtαacz(t)+∑k=0n-1z(k)(a)k!(x-a)k=(15)ItαaDtαacz(t)=(16)ItαaDtαaz(t)≈(18)Itαa(CTΦm(t))=CTItαaΦm(t)≈(12)CTFαΦm(t)=CαTΦm(t),
where CαT=CTFα.
Now, by substituting (17)–(21) into (14), we obtain
(22)CTΦm(t)=FTΦm(t)+BTΦm(t)Φm(t)TFαTC+ATΦm(t)×(CTFαΦm(t)Φm(t)TFαTC)=FTΦm(t)+BTΦm(t)Φm(t)TCα+ATΦm(t)×(CαTΦm(t)Φm(t)TCα).
Then, from Lemma 3.5 in [20] we have
(23)CαTΦm(x)Φm(x)T≈Φm(x)TC~α,(24)ATΦm(x)Φm(x)T≈Φm(x)TA^,(25)BTΦm(x)Φm(x)T≈Φm(x)TB^.
Therefore we can reduce (22) by (23)–(25) as
(26)CTΦm(t)=FTΦm(t)+Φm(t)TB^Cα+Φm(t)TA^C~αCα.
Finally, we obtain the following nonlinear system of algebraic equation:
(27)C=F+B^Cα+A^C~αCα,
such that by solving this system we can obtain the vector C. Then, we can get
(28)y(t)≈y^(t)+CTFαΦm(t).
5. Convergence Analysis
In this section, we investigate the convergence analysis for the method presented in Section 4.
The problem (14) changes to the following problems
(29)Dtαaz(t)=f~(t)+b~(t)(ItαaDtαaz(t))+a~(t)(ItαaDtαaz(t))2,a<t≤b,
since
(30)z(t)=(10)ItαaDtαacz(t)+∑k=0n-1z(k)(a)k!(x-a)k=(15)ItαaDtαacz(t)=(16)ItαaDtαaz(t).
By taking u(t)=Dtαaz(t) we obtain the following fractional integral equation:
(31)u(t)=f~(t)+b~(t)Itαau(t)+a~(t)(Itαau(t))2,a<t≤b.
If we use the approximation u(t)≈CTΦm(t), then the problem (31) from space C1[0,1] reduces to the following problem in space Sm=Span{B0,m(t),B1,m(t),…,Bm,m(t)}:
(32)CTΦm(t)=f~(t)+b~(t)Itαa(CTΦm(t))+a~(t)(Itαa(CTΦm(t)))2,a<t≤b.
Now, similar to Theorem 6.1 in [21], we propose the next theorem.
Theorem 6.
Suppose that u*(t)∈C1[0,1] is the exact solution of (31) and μm=J[um]=Minu∈SmJ[u], where
(33)J[u]=∥u(t)-f~(t)-b~(t)Itαau(t)-a~(t)(Itαau(t))2∥.
Then one has μm→0 asm→∞ (i.e., um(t)→u*(t) as m→∞).
Proof.
By substituting (33) instead of (41) in [21], we can use the proof of Theorem 6.1 in [21].
Theorem 7.
Suppose that u*(t)∈C1[0,1] is the exact solution of (31) and um(t)∈Sm is the obtained solution of (26). Then one has um(t)→u*(t) asm→∞.
Proof.
Substituting (17), (19), (20), and (12) in (32) we have
(34)CTΦm(t)=(FTΦM(t)+ef~M)+(BTΦM(t)+eb~M)(CT(FαΦM(t)+EIM))+(ATΦM(t)+ea~M)(CT(FαΦM(t)+EIM))2.
From Lemmas 2.3 and 4.1 in [21], we have ef~M,eb~M,ea~M,EIM→0 as M→∞. So we can observe that as M increases, (34) gets close to (32). Now, by taking M=m we propose the following problem that gets close to (32) as m increases:
(35)CTΦm(t)=FTΦm(t)+BTΦm(t)Φm(t)TCα+ATΦm(t)(CαTΦm(t)Φm(t)TCα).
Then by (23) and (25), (35) reduces to the following equation:
(36)CTΦm(t)=FTΦm(t)+(ΦM(t)TB^+E^BM)Cα+ATΦm(t)((ΦM(t)TC~α+E~CαM)Cα).
Equation (36) gets to (35) as M→∞, because from Lemma 3.1 in [21] E^BM,E~CαM→0 as M→∞. Then by deleting E^BM, E~CαM, taking M=m, and using (24) in (36), we have
(37)CTΦm(t)=FTΦm(t)+Φm(t)TB^Cα+(ΦM(t)TA^+E^AM)C~αCα,
where, from Lemma 3.1 in [21], E^AM→0 as M→∞. Now, by taking M=m and deleting E^AM in (37), we get (26). Obviously, if u~m(t) is solution of (32), then we have u~m-um→0 as m→∞.
On the other hand, from Theorem 6 we obtained u~m→u* as m→∞. Therefore we can write um→u* as m→∞ and the proof is complete.
6. Illustrative Numerical Examples
In this section, we apply our method with m=10 (BPs of degree m=10) to solve the following examples. We define ym(t) and y(t) for the approximate solution and the exact solution, respectively.
Example 1.
Consider the nonlinear Riccati differential equation [14]:
(38)Dtα0y(t)=-y(t)2+1,0<t≤1,
subject to the initial condition as y(0)=0. The exact solution of the equation for α=1 is given as
(39)y(t)=e2t-1e2t+1.
Numerical results compared to [14] are given in Table 1 and also Figure 1 shows the absolute error for our method for α=1 and Figure 2 shows behavior y10(t) for different values of α.
Numerical results for α=1 and m=10 in Example 1 with comparison to exact solution and [14].
t
Exact
Present methodm=10
Reference [14]
0.1
0.099668
0.099668
0.099668
0.2
0.197375
0.197375
0.197375
0.3
0.291313
0.291313
0.291313
0.4
0.379949
0.379949
0.379944
0.5
0.462117
0.462117
0.462078
0.6
0.537050
0.537050
0.536857
0.7
0.604368
0.604368
0.603631
0.8
0.664037
0.664037
0.661706
0.9
0.716298
0.716298
0.709919
1
0.761594
0.761594
0.746032
Behavior of the absolute error function in Example 1 for α=1 and m=10.
Behavior of y10(t) in Example 1 for different α and exact solution for α=1.
Example 2.
Consider the following quadratic Riccati differential equation of fractional order [14]
(40)Dtα0y(t)=2y(t)-y(t)2+1,0<t≤1,
subject to the initial condition as y(0)=0. The exact solution of the equation for α=1 is given as
(41)y(t)=1+2tanh(2t+12log(2-12+1)).
Numerical results compared to [14] are given in Table 2 and also, Figure 3 shows the absolute error for our method for α=1 and Figure 4 shows behavior y10(t) for different values of α.
Numerical results for α=1 and m=10 in Example 2 with comparison to exact solution and [14].
t
Exact
Present methodm=10
Reference [14]
0.1
0.110295
0.110295
0.110294
0.2
0.241977
0.241977
0.241965
0.3
0.395105
0.395105
0.395106
0.4
0.567812
0.567812
0.568115
0.5
0.756014
0.756014
0.757564
0.6
0.953566
0.953566
0.958259
0.7
1.152949
1.152949
1.163459
0.8
1.346364
1.346364
1.365240
0.9
1.526911
1.526911
1.554960
1
1.689499
1.689499
1.723810
Behavior of the absolute error function in Example 2 for α=1 and m=10.
Behavior of y10(t) in Example 2 for different α and exact solution for α=1.
Example 3.
Consider the nonlinear fractional differential equation:
(42)Dtα0y(t)=y2(t)+1,0<t≤1,n-1<α≤n,
with the initial conditions
(43)y(k)(0)=0,k=0,1,…,n-1.
This problem has been studied by using ADM [25], FDTM [26], and BPFs [27]. Our results with α=1.5, α=2.5 are compared to [25–27] in Tables 3 and 4. Therefore, we see that our method is very effective and obtained solutions that are in good agreement with the results in [25–27]. Also, Figure 5 shows behavior ym(t) for different values of α.
Numerical results for α=1.5 in Example 3 with comparison to [25–27].
t
Present methodm=10
ADM [25]
FDTM [26]
BPFs [27]
0.1
0.023779
0.023790
0.023790
0.023800
0.2
0.067336
0.067330
0.067330
0.067335
0.3
0.123886
0.123896
0.123896
0.123900
0.4
0.191373
0.191362
0.191362
0.191368
0.5
0.268851
0.268856
0.268856
0.268862
0.6
0.356235
0.356238
0.356238
0.356244
0.7
0.453958
0.453950
0.453950
0.453956
0.8
0.562999
0.563007
0.563007
0.563014
0.9
0.685066
0.685056
0.685056
0.685067
1
0.822540
0.822511
0.822509
0.822525
Numerical results for α=2.5 in Example 3 with comparison to [25, 26].
t
Present methodm=10
ADM [25]
FDTM [26]
0.1
0.000952
0.000952
0.000952
0.2
0.005383
0.005383
0.005383
0.3
0.014833
0.014833
0.014833
0.4
0.030450
0.030450
0.030450
0.5
0.053197
0.053197
0.053197
0.6
0.083925
0.083925
0.083925
0.7
0.123412
0.123412
0.123412
0.8
0.172391
0.172391
0.172391
0.9
0.231574
0.231574
0.231574
1
0.301676
0.301676
0.301676
Plot of y10(t) for different α in Example 3.
7. Conclusion
In this paper, we proposed a numerical method for solving the fractional quadratic Riccati differential equations by the operational matrices of the Bernstein polynomials. We applied operational matrix for fractional integration in the Riemann-Liouville sense. Then by using this matrix and operational matrix of product, we reduced the fractional quadratic Riccati differential equation to a system of algebraic equations that can be solved easily. Finally, examples have been simulated to demonstrate the high performance of the proposed method. We saw that the results were in good agreement with the analytical solutions and the solutions in the open literatures. Also, we observed that the solutions approach to classical solutions as the order of the fractional derivatives approaches 1, for fixed m.
ReidW. T.197286New York, NY, USAAcademic PressMathematics in Science and EngineeringMR0357936KravchenkoV. V.2009Basel, SwitzerlandBirkhäuserFrontiers in Mathematics10.1007/978-3-0346-0004-0MR2510077ConteR.MusetteM.Link between solitary waves and projective Riccati equations1992252156095623MR1192020DiethelmK.20102004Berlin, GermanySpringerLecture Notes in Mathematics10.1007/978-3-642-14574-2MR2680847AbbasbandyS.Iterated He's homotopy perturbation method for quadratic Riccati differential equation2006175158158910.1016/j.amc.2005.07.035MR2216365AbbasbandyS.Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian's decomposition method2006172148549010.1016/j.amc.2005.02.014MR2197917AbbasbandyS.A new application of He's variational iteration method for quadratic Riccati differential equation by using Adomian's polynomials20072071596310.1016/j.cam.2006.07.012MR2332947TanY.AbbasbandyS.Homotopy analysis method for quadratic Riccati differential equation20081335395462-s2.0-3484888079810.1016/j.cnsns.2006.06.006GengF.LinY.CuiM.A piecewise variational iteration method for Riccati differential equations20095811-122518252210.1016/j.camwa.2009.03.063MR2557400GengF. Z.LiX. M.A new method for Riccati differential equations based on reproducing kernel and quasilinearization methods20122012860374810.1155/2012/603748MR2898051KhanN. A.AraA.JamilM.An efficient approach for solving the Riccati equation with fractional orders20116192683268910.1016/j.camwa.2011.03.017MR2795392MomaniS.ShawagfehN.Decomposition method for solving fractional Riccati differential equations200618221083109210.1016/j.amc.2006.05.008MR2282552CangJ.TanY.XuH.LiaoS.-J.Series solutions of non-linear Riccati differential equations with fractional order20094011910.1016/j.chaos.2007.04.018MR2517912OdibatZ.MomaniS.Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order200836116717410.1016/j.chaos.2006.06.041MR2372056LiC.WangY.Numerical algorithm based on Adomian decomposition for fractional differential equations200957101672168110.1016/j.camwa.2009.03.079MR2528411JafariH.TajadodiH.He's variational iteration method for solving fractional Riccati differential equation20102010876473810.1155/2010/764738MR2607724MerdanM.On the solutions fractional Riccati differential equation with modified Riemann-Liouville derivative201220121734608910.1155/2012/346089MR2925289YüzbaşıS.Numerical solutions of fractional Riccati type differential equations by means of the Bernstein polynomials2013219116328634310.1016/j.amc.2012.12.006MR3018474CheneyE. W.19982ndProvidence, RI, USAAMS ChelseaMR1656150AlipourM.RostamyD.BaleanuD.Solving multi-dimensional FOCPs with inequality constraint by BPs operational matrices201210.1177/1077546312458308RostamyD.AlipourM.JafariH.BaleanuD.Solving multi-term orders fractional differential equations by operational matrices of BPs with convergence analysis2013652334349PodlubnyI.1999198New York, NY, USAAcademic PressMR1658022KilbasA. A.SrivastavaH. M.TrujilloJ. J.2006San Diego, Calif, USAElsevierMR2218073BaleanuD.DiethelmK.ScalasE.TrujilloJ. J.20123Hackensack, NJ, USAWorld ScientificSeries on Complexity, Nonlinearity and Chaos10.1142/9789814355216MR2894576ShawagfehN. T.Analytical approximate solutions for nonlinear fractional differential equations20021312-351752910.1016/S0096-3003(01)00167-9MR1920243ArikogluA.OzkolI.Solution of fractional differential equations by using differential transform method20073451473148110.1016/j.chaos.2006.09.004MR2335397LiY.SunN.Numerical solution of fractional differential equations using the generalized block pulse operational matrix20116231046105410.1016/j.camwa.2011.03.032MR2824692