Difference Equations, Article ID 359093, 8 pages http://dx.doi.org/0.55/04/359093 Research Article Hermite Wavelet Method for Fractional Delay Differential Equations Umer Saeed and Mujeeb ur Rehman School of Natural Sciences, National University of Sciences and Technology, Sector H-, Islamabad, Pakistan Correspondence should be addressed to Umer Saeed; umer.math@gmail.com Received 4 April 04; Accepted June 04; Published July 04 Academic Editor: Honglei Xu Copyright 04 U. Saeed and M. ur Rehman. 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. We proposed a method by utilizing method of steps and Hermite wavelet method, for solving the fractional delay differential equations. This technique first converts the fractional delay differential equation to a fractional nondelay differential equation and then applies the Hermite wavelet method on the obtained fractional nondelay differential equation to find the solution. Several numerical examples are solved to show the applicability of the proposed method.. Introduction The future state of a physical system depends not only on the present state but also on its past history. Functional differential equations provide a mathematical model for such physical systems in which the rate of change of the system may depend on the influence of its hereditary effects. Delay differential equations have numerous applications in mathematical modeling [], for example, physiological and pharmaceutical kinetics, chemical kinetics, the navigational control of ships and aircrafts, population dynamics, and infectious diseases. Delay differential equation is a generalization of the ordinary differential equation, which is suitable for physical system that also depends on the past data. During the last decade, several papers have been devoted to the study of the numerical solution of delay differential equations. Therefore different numerical methods [ 7] have been developed and applied for providing approximate solutions. Method of steps is easy to understand and implement. In the method of steps [8], we convert the delay differential equation to a nondelay differential equation. The method of steps is utilized in [9] for solving integer order delay differential equations. Hermite wavelet method [0] is implemented for finding the numerical solution of the boundary value problems and compares the obtained solutions with exact solution. In [], authors utilized the physicists Hermite wavelet method for solving linear singular differential equations. According to our information, Hermite wavelet method has not been implemented for delay differential equations. Inthepresentwork,weestablishedatechniquebycombining both the method of steps and the Hermite wavelets method for solving the fractional delay differential equation. We also implemented the Hermite wavelet method for solving fractional delay differential equation, as described in Example 6, which was not implemented before. Shifted Chebyshev nodes are used as the collocation points. Comparison of solutions by these two methods, proposed method and Hermite wavelet method, with each other and with exact solution are also presented.. Preliminaries.. Hermite Wavelets. The Hermite polynomials H m (x), of order m, aredefinedontheinterval[, ] and given by the following recurrence formulae: H 0 (x) =, H (x) =x, H m+ (x) =xh m (x) mh m (x), m=,,3,... ()
Difference Equations The polynomials H m (x) are orthogonal with respect to the weight function e x ;thatis, 0, m e x H m (x) H n (x) dx = { =n; n! n π, m = n. The discrete wavelets transform is defined as () ψ j,k (x) = j/ ψ ( j x k). (3) The set ψ j,k forms an orthogonal basis of L (R).Thatis, ψ j,k (x),ψ l,m (x) =δ jl δ km. (4) The Hermite wavelets are defined on interval [0, ) by ψ n,m (x) { k/ = n! { n π H m ( k n x n), x< n+, k k { 0, elsewhere, (5) where k =,, 3,..., is the level of resolution, n =,, 3,..., k, n = n, is the translation parameter, m =,,...,M is the order of the Hermite polynomials, M>0... Function Approximations. We can expand any function L [0, ) into truncated Hermite wavelet series as y (x) a nm ψ n,m (x) = a T Ψ (x), (6) where a and Ψ are m ( m = k M) matrices, given by a =[a 0,a,...,a M,a 0,a,...,a M,..., a k 0,a k,...,a ]T, Ψ (x) =[ψ,0 (x),ψ, (x),...,ψ M (x),ψ,0 (x), ψ, (x),...,ψ,m (x),..., ψ k,0 (x),ψ k, (x),...,ψ k,m (x)] T. 3. Convergence Analysis Let L ([0, ]) be a Hilbert space for which ψ n,m (x) form an orthonormal sequence in L ([0, ]).Let L ([0, ]);we have y (x) (7) a nm ψ n,m (x), (8) where a nm =,ψ n,m (x) is an inner product of and ψ n,m (x).equation(8)canbewrittenas y (x) k n= M y (x),ψ n,m (x) ψ n,m (x). (9) m=0 For simplicity, let j = M(n ) + m + ;wecanwrite(9)as y (x) = m j= y (x),ψ j (x) ψ j (x) = m j= a j ψ j (x) = a T Ψ (x), (0) where a j = a mn, ψ j (x) = ψ m,n (x), m = k M,anda = [a,a,...,a m ] T ; Ψ(x) = [ψ (x), ψ (x),...,ψ m (x)] T. By following the procedure in [], we obtained the convergenceofalltheorthogonalwaveletmethodsforall levels of resolution k;thatis, m j= a jψ j (x) converges to, as m. Since m = k M and method converges if m ;that is, when we use higher order Hermite polynomials M or use higher level of resolution k or use both higher M and k, we get more accurate results. 4. Procedure for Implementation of Proposed Scheme The method of steps [8] isusedtoconvertthediscrete delay differential equations to nondelay differential equations on a given interval. Consider the following fractional delay differential equation with discrete delay: c D α y (x) =g(x) +f(,y (x),y(qx τ),y (qx τ)), a x b, <α, y (x) =φ(x), b x a, () where g(x) is a source function and f is a continuous linear or nonlinear function. Also q is constant, τ is delay, and qx τis called delay argument. The delay τ(x, ) is called constant delay, time dependent delay, and state dependent delay if the delay τ(x, ) is constant, function of time x, and function of time x and,respectively. Proposed method consists of two methods, method of stepsandhermitewaveletmethod.wefirstimplementthe method of steps to the fractional delay differential equation () and get the fractional nondelay differential equation by utilizing initial function, φ(x), andthenweutilizethe Hermite wavelet method for solving the obtained fractional nondelay differential equation. 4.. Method of Steps. In the fractional delay differential equation the solution is known on [ b, a],sayφ(x),and call this solution y 0 (x);thatis,y 0 (qx τ)=φ(qx τ),which is known. Now the fractional delay differential equation on [a, b] takes the form c D α y (x) =g(x) +f(,y (x),y 0 (qx τ), y 0 (qx τ)), a x b, <α, subject to the initial conditions y(a) = φ(a), y (a) = φ (a). ()
Difference Equations 3 It is a fractional nondelay differential equation because y 0 (qx τ)and y 0 (qx τ) are known. 4.. Hermite Wavelet Method for Fractional Nondelay Differential Equation. We solve the obtained fractional nondelay differential equation () on[a, b] by using the Hermite wavelet method. The procedure for implementation of Hermite wavelet method for fractional differential equation is as follows. Step. Approximate the unknown function of () by the Hermite wavelet method as y (x) a nm ψ n,m (x). (3) Step. Substitute (3)in() to get the residual a c nm Dα ψ n,m (x) =g(x) +f( a nm ψ n,m (x), a nm ψ n,m (x), y 0 (qx τ),y 0 (qx τ)). (4) Step 3. Set the residual (4) tobeequaltozeroatthesetof Chebyshev nodes, x i = ((b a)/) cos((i + )π/ k M) + (a + b)/, i=0,,,..., M,oninterval[a, b]. Consider a c nm Dα ψ n,m (x i ) g(x i ) f( a nm ψ n,m (x i ), a nm ψ n,m (x i), y 0 (qx i τ),y 0 (qx i τ))=0. (5) We obtain k M pequations, where p is the number of conditions of the delay equation. According to (), two conditions are given, p=,soweget k M equations from (5) by using Chebyshev nodes x i.twomoreequations are obtained from the conditions of (); that is, y (a) =φ(a), y (a) =φ (a), a nm ψ n,m (a) =φ(a), a nm ψ n,m (a) =φ (a). (6) We obtained k M equationseitherlinearornonlinear along with k M unknown coefficients a nm,whichissolved by Newton iterative method to get a nm s and use it in (3) to get the approximate solution. Denote the obtained solution as y (x), which is defined on [a, b]. Continue the procedure for the subsequent interval; delay differential equation on [b, b] becomes c D α y (x) =h(x) +f(,y (x),y (qx τ),y (qx τ)), b x b, <α, (7) subject to the initial conditions y(b) = y (b), y (b) = y (b), which is again a fractional nondelay differential equation and solve it by the Hermite wavelet method to get y (x) on [b, b]. This procedure may be continued for subsequent intervals. 5. Numerical Solutions In this section, we utilize the proposed scheme for finding the numerical solution of linear and nonlinear fractional delay differential equations. The notations y pro, y exact,ande abs represent the solution by proposed method, exact solution, and their absolute error, respectively. We use the results up to 00 decimal places. Through this work we use Caputo derivatives. For the details of fractional derivatives and integrals we refer the readers to [3]. 5.. Linear Delayed Fractional Differential Equations Example. Consider the fractional delay differential equation c D α y (x) + y(x τ) = Γ (3 α) x α Γ ( α) x α +τx τ τ, x > 0, 0 < α <, y (x) =0, x 0. (8) The exact solution [], when α=,is = x x. The results obtained by the proposed method, by taking α = 0.5, k =, M = 3,andτ = 0.0e x,areshown in Figure alongwiththeexactsolution.table indicates that results obtained from the proposed method are closer to exact solution and better than the method []. E abs and error represent the absolute error by proposed method and method [], respectively. Example. Consider the following fractional delay differential equation: c D α y (x) = y(x 0.3) +e x+0.3, (9) 0 x, <α 3, subject to the initial conditions y(0) =, y (0) =, and y (0) =. The exact solution, when α=,is = e x.
4 Difference Equations Table : Comparison of the solution by proposed method at M=3, with Exact solution and method []. x τ = 0., α = 0. τ = 0.0e x, α = 0.9 τ = 0.0e x, α = 0. E abs Error [] E abs Error [] E abs Error [] 5.9000E 97 0.006 4.0000E 98 0.000.7000E 97 0.0074 0.7000E 97 0.034.000E 97 4.75E 4 4.0000E 97 0.008 50 3.0000E 96 0.0690 0.0000E + 00 7.0303E 4 5.0000E 96 0.005 Table : Comparisonof the solution byproposed method,m =5, with Exact solution and Adomian decomposition method [3], at α=3. x y pro y exact E abs E ADM [3] 0.0.0000.0000 0.00E + 00 8.5E 4 0. 0.887 0.887 9.E 4 3.83E 4 0.4 0.6703 0.6703 6.90E 4.68E 3 0.6 0.5488 0.5488.38E 40 6.00E 4 0.8 0.4493 0.4493 6.0E 4 6.66E 5.0 0.3679 0.3679 3.58E 36 4.57E 4 500 000 0.9 0.8 500 0.7 000 0.6 0.5 500 0.4 0 0 0 0 30 40 50 at α = 0.5 Exact solution at α = 0.5 Figure : Solutions by the proposed method, τ = 0.0e x, and exact solutions, at α = 0.5. According to Table, proposed method provides more accurate results as compared to Adomian decomposition method [3]. These results are obtained by fixing M = 5 and k=,atα=3.solutionsbytheproposedmethodat different values of α are shown in Figure,whichshowsthat solutions by proposed method at different α converge to the exact solution at α=3,whenα approaches to 3. Example 3. Consider the fractional pantograph equation c D α y (x) = ex/ y( x )+ y (x), 0<α,0 x, y (0) =. The exact solution, when α=,is = e x. (0) 0 0 0. 0.4 0.6 0.8 at α =.8 at α =.7 at α =.6 at α =.5 at α=3 Exact solution at α=3 Figure : Solutions by the proposed method at different α, M=7, k=, and exact solution at α=3. By fixing k = and M = 4, we plot the solutions by proposed method at different values of α and exact solution at α =,asshowninfigure3. Itshowsthat proposed solution approaches to the exact solution while α approaches to. Comparison of proposed solution y pro at α=, M=30, and k =, with Adomian decomposition method [3] and the spline function technique [4],isshowninTable3. The notations E abs, E ADM,andE spline represent the absolute error by proposed method, Adomian decomposition method, and the spline function technique, respectively. We can get more accurate results while increasing M.
Difference Equations 5 Table 3: for M=30, k=. α= x y pro y exact E abs E ADM [3] E spline [4] 0..40.40.9659E 5 0.00 3.0E 5 0.4.498.498.5474E 5.E 6 7.54E 5 0.6.8.8.5683E 5.E 6.39E 4 0.8.554.554.03409E 50.33E 5.3E 4.0.788.788 8.866E 49 4.88E 5 3.9E 4 Table 4: Absolute errors using proposed method at α=and M=7. α= SLC method [5] RKHSM[6] RKTmethod[7] M=7 x E abs 0. 4.7E 7.4E 4 8.68E 4.3E 4 0..70E 7.7E 4.49E 3 3.4E 5 0.3 5.94E 7 9.45E 4.90E 3.74E 4 0.4 8.0E 7 7.59E 4.6E 3.79E 4 0.5 8.7E 7 6.03E 4.8E 3 3.33E 4 0.6.95E 6 4.73E 4.3E 3 3.94E 4 0.7.56E 6 3.64E 4.7E 3 3.56E 4 0.8 8.80E 7.75E 4.7E 3.35E 5 0.9.03E 6.03E 4.03E 3.0E 3.0.3E 6.43E 4.86E 3.58E Example 4. Consider the following fractional neutral functional differential equation with proportional delay: c D α y (x) 5 4.5 = + 0.y (0.8x) + 0.5 c D α y (0.8x) + (0.3x 0.5) e 0.8x +e x, x 0,0<α, y (0) =0, () 4 3.5 3.5 which has the exact solution, when α=,whichisxe x. We implement the proposed method by fixing k =, M = 5, at different values of α, the results are shown in Figure 4 along with the exact solution at α =. Table 4 shows that proposed method gives more accurate results as compared to the spectral shifted Legendre-Gauss collocation (SLC)method [5], the reproducing kernel Hilbert space method (RKHSM) [6],and a Runge-Kutta-type (RKT) method [7]. 5.. Nonlinear Delayed Fractional Differential Equations Example 5. Consider the fractional nonlinear delay differential equation c D α y (x) = y ( x ), 0 x,<α, () subject to the initial conditions, y(0) = and y (0) = 0.The exact solution [4], when α=, is = cos(x)..5 0 0. 0.4 0.6 0.8 at α = 0.5 at α = 0.6 at α = 0.7 at α = 0.8 at α = 0.9 at α= Exact solution at α= Figure 3: Solutions by the proposed method at different α, M=4, k=, and exact solutions at α=. at different values of α is plotted in Figure 5, whichshowsthatproposedsolution converges to the exact solution when α approaches to. According to Table 5, absolute error reduces while increasing M.
6 Difference Equations 0.4 0.35 0.3 0.5 0. 0.5 0. 0.05 subject to the initial conditions, y(0) =. The exact solution [5], when α=,is = e x. 5.4. Hermite Wavelet Method for Fractional Delay Differential Equation. We can approximate the solution of (3) by the Hermite wavelet method as y (x) a nm ψ n,m (x). (4) In the delay equations, we also have to approximate the delay unknown function y(x/) in terms of the Hermite wavelet series at delay time as 0 0 0. 0.4 0.6 0.8 at α = 0.7 at α = 0.8 at α = 0.9 at α = 0.95 at α= Exact solution at α= Figure 4: Solution by the proposed method at different α, M=5 and exact solution at α=. 0.9 0.8 0.7 0.6 0.5 0.4 0.3 at α =.65 at α =.75 at α =.9 0. 0. 0 0 0. 0.4 0.6 0.8 at α= Exact solution at α= Figure 5: Solution by the proposed method at different α, M=5 and exact solution at α=. 5.3. Comparison of Proposed Method and Hermite Wavelet Method Example 6. Consider the fractional nonlinear neutral delay differential equation c D α y (x) = y (x) + y(x ) c D α y( x ), x 0, 0<α, (3) y( x ) k M a nm ψ n,m ( x ). (5) We call this series as the delay Hermite wavelet series. Substituting (4)and(5)in(3), we get the residual as a c nm Dα ψ n,m (x) k M = nm ψ n,m (x) + a n= m=0 ( a c nm Dα ψ n,m ( x )). a nm ψ n,m ( x ) (6) Set the residual (6) to be equal to zero at the set of Chebyshev nodes, x i = (/) cos((i + )π/ k M) + (/), k=0,,,...,m,oninterval[0, ],weget a c nm Dα ψ n,m (x i ) n= ( m=0 a nm ψ n,m ( x i ) a c nm Dα ψ n,m ( x i a nm ψ n,m (x i ) )) = 0. (7) We get equations from (7) byusingchebyshev nodes x i. One more equation is obtained from the condition of (3); that is, y (0) =, a nm ψ n,m (0) =. (8) We obtained k M nonlinear equations along with k M unknown coefficients a nm,whichissolvedbynewtoniterative method to get a nm s and used in (4) to get the approximate solution by Hermite wavelet method.
Difference Equations 7 Table 5: Absolute errors by using proposed method at different M and α=. α= M=5 M=5 M=30 M=40 x E abs E abs E abs E abs 0..7077E 3.96989E 4 8.0964E 5 5.4668E 60 0. 4.49566E 3.8499E 4 7.79780E 5.568E 59 0.3.6700E 6.69E 4.4878E 5.5773E 59 0.4 3.6080E 3.5484E 4.38398E 5 5.4980E 60 0.5.8700E.5484E 4 4.954E 5.09577E 58 0.6 4.60073E 4 7.94684E 4 6.7556E 5 9.58443E 59 0.7 4.895E.8456E 40.E 50 3.69647E 57 0.8.358E.3608E 40.58769E 50.69940E 55 0.9 5.88496E 4.947E 40 5.3588E 50 5.9603E 54.0 5.3053E 9.539E 37.9993E 47 3.70755E 53 Table6:ComparisonofproposedandHermitewaveletmethodsM=9,k=. α= x y HWM y pro y exact E HWM E pro 0.0.0000000000.0000000000.0000000000 0.00000E + 00 0.00000E + 00 0..057098.057098.057098.3388E 5.37737E 0..40758.40758.40758 3.933E.4859E 0.3.3498588076.3498588076.3498588076.83949E.3830E 0.4.49846977.49846976.49846976.6536E 6.3057E 0.5.6487707.6487707.6487707 8.4958E 3.08447E 0.6.888004.888004.888004 5.79353E.57380E 0.7.03757076.03757075.03757075 8.9374E 4.348E 0.8.5540985.5540984.5540985 3.39446E 6.0808E 0.9.459603.459603.459603 6.779E 6.7473E.0.788876.788876.788885 8.49E 0 8.35594E 0 We fix M=9, k=, and implement the Hermite wavelet andproposedmethodto(3). The results are shown in Table 6 along with the absolute errors. According to Table 6, Hermite wavelet and proposed methods give good results; that is, for some points proposed method is more accurate as compared tohermitewaveletmethodandviceversa. For the problem (3), run time of proposed method and Hermite wavelet method is. and 6. seconds, respectively. Proposed method is more efficient than the Hermite wavelet method. y HWM and E HWM represent the solution by Hermite wavelet method and their absolute error, respectively. For this purpose, we use Maple 3 in system with Core Duo CPU.00 GHz and RAM.50 GB. 6. Conclusion It is shown that proposed method gives excellent results when applied to different fractional linear and nonlinear delay differential equations. The results obtained from the proposed method are more accurate and better than the results obtained from other methods, as shown in Tables 4. The solution of the fractional delay differential equation converges to the solution of integer delay differential equation, as shown in Figures 5. According to the convergence analysis, error by the proposed method reduces while increasing M,as shown in Table 5.Table 6 indicates that both Hermite wavelet method and proposed method give good results. Proposed method is more efficient than the Hermite wavelet method. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgment The authors are grateful to the anonymous reviewers for their valuable comments which led to the improvement of the paper. References [] V. Lakshmikantham and S. Leela, Differential and Integral Inequalities, Academic Press, New York, NY, USA, 969. [] Z. Wang, A numerical method for delayed fractional-order differential equations, JournalofAppliedMathematics,vol. 03, Article ID 5607, 7 pages, 03.
8 Difference Equations [3] D. J. Evans and K. R. Raslan, The Adomian decomposition method for solving delay differential equation, International JournalofComputerMathematics,vol.8,no.,pp.49 54,005. [4]A.El-Safty,M.S.Salim,andM.A.El-Khatib, Convergence of the spline function for delay dynamic system, International Computer Mathematics, vol.80,no.4,pp.509 58, 003. [5] A. H. Bhrawy, L. M. Assas, E. Tohidi, and M. A. Alghamdi, A Legendre-Gauss collocation method for neutral functionaldifferential equations with proportional delays, Advances in Difference Equations, vol. 03, article 63, 03. [6] X.LvandY.Gao, TheRKHSMforsolvingneutralfunctionaldifferential equations with proportional delays, Mathematical Methods in the Applied Sciences,vol.36,no.6,pp.64 649,03. [7] W. Wang, Y. Zhang, and S. Li, Stability of continuous Runge- Kutta-type methods for nonlinear neutral delay-differential equations, Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 33,no.8,pp.339 339,009. [8] H. Smith, An Introduction to Delay Differential Equations with Applications to the Life Sciences, Springer, New York, NY, USA, 0. [9] C. E. Falbo, Some Elementary Methods for Solving Functional Differential Equations, http://www.mathfile.net/hicstat FDE.pdf. [0] A. Ali, M. A. Iqbal, and S. T. Mohyud-Din, Hermite wavelets method for boundary value problems, International Modern Applied Physics,vol.3,no.,pp.38 47,03. [] M. Usman and S. T. Mohyud-Din, Physicists Hermite wavelet method for singular differential equations, International Journal of Advances in Applied Mathematics and Mechanics, vol., no., pp. 6 9, 03. [] S. G. Venkatesh, S. K. Ayyaswamy, and S. Raja Balachandar, The Legendre wavelet method for solving initial value problems of Bratu-type, Computers & Mathematics with Applications, vol. 63, no. 8, pp. 87 95, 0. [3] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 999. [4] F. Karakoç andh.bereketoğlu, Solutions of delay differential equations by using differential transform method, International Computer Mathematics,vol.86,no.5,pp.94 93, 009. [5] M. M. Khader and S. T. Mohamed, Numerical treatment for first order neutral delay differential equations using spline functions, Engineering Mathematics Letters,vol.,no.,pp.3 43, 0.
Advances in Operations Research Advances in Decision Sciences Applied Mathematics Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Discrete Mathematics Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis International Stochastic Analysis Optimization