NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
|
|
- Frank Hancock
- 5 years ago
- Views:
Transcription
1 TWMS J Pure Appl Math V5, N2, 24, pp NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay differential equations (DDEs) A discretized form of DDEs at collocation points based on some useful properties of Haar wavelets transforms original problem into a nonlinear algebraic equations Finally, the numerical experiments are given to demonstrate the conclusions Keywords: Delay differential equation, haar wavelet, collocation method, numerical solution AMS Subject Classification: 34K28, 42C5, 65L3 Introduction Delay differential equations (DDEs) arise in many areas of mathematical modelling, as physiological and pharmaceutical kinetics, chemical kinetics, population dynamics, the navigational control of ships and aircraft, infectious diseases, and more general control problems Given the importance of applications of DDEs, many researchers have addresssed this issue [, 4, 2, 7] Since analytical solutions of the DDEs may be obtained only in very restricted cases, many methods have been proposed for numerical approximation of them Some of these techniques include Hermite interpolation [3], θ-methods [], Runge-Kutta methods [, 9], Parallel continuous Runge-Kutta methods [2], Spline collocation methods [5], Collocation methods [6], linear multistep methods [9] Recently, Haar wavelets as a useful mathematical tool, have been applied extensively for signal processing in communications, physics researches and optimal control problems [3] Haar wavelets have the simplest orthogonal series with compact support and this characteristic introduces Haar wavelets as a good candidate for application to differential equations In this paper, using the benefits of Haar wavelets, a novel collocation approach is introduced to find approximate solutions of the following delay differential equation ẋ(t) = f ( t, x(t), x(t τ ), x(t τ 2 ),, x(t τ s ) ) < t t f, x(t) = ξ(t) t, () where x(t) R n, f : R (s+)n+ R n is a nonlinear smooth function and τ i >, i =,, s are constant delays In the next sections, first Haar wavelets and their properties are introduced Then the approximation of a function by Haar wavelets is discussed By introducing operational integration matrix and delay operational matrix, a discretization method is established Finally, by some numerical examples the proficiency of the given approach is examined School of Mathematics and Computer Science, Damghan University, Damghan, Iran soasmath84@gmailcom, borzabadi@duacir Manuscript received January 23 22
2 222 TWMS J PURE APPL MATH, V5, N2, 24 2 Haar wavelets and its properties The orthogonal set of Haar wavelets φ i (t) is a group of square waves with magnitude + or - in some intervals and zeros elsewhere, φ (t) =, t <, (2) { if t < φ (t) = 2, if 2 t <, (3) if k t < k +, φ i (t) = φ (2 j 2 j 2 j 2 j+ t k) = if k + t < k +, (4) 2 j 2 j+ 2 j 2 j otherwise i = 2 j + k, j =,, M, k =,, 2 j, integer m = 2 j, (j =,,, J), indicates the level of the wavelet, k =,,, m is the translation parameter, maximal level of resolution is J and the maximal value of i is M = 2 J+ A simple calculation shows that { φ i (t)φ l (t)dt = M if i = l = 2j + k, if i l Consequently, the functions φ i (t) are orthogonal This allows us to transform any function square integrable on the interval time [,] into Haar wavelets series 2 Function approximation We just pointed out that a square integrable function can be expressed in terms of Haar orthogonal basis on interval l [, ] However, before the procession to this transfer, it is necessary to unify the time interval Using a linear transformation, the actual time t can be expressed as a function of l via t = [(t f t )l + t ], where t is the initial time and t f is the final time in a square integrable function f(t) Any function f(t) which is square integrable in the interval [, ] can be expanded in a Haar series with an infinite number of terms f(t) = a i φ i (t), i = 2 j + k, j, k < 2 j, t [, ), (6) i= where the Haar coefficients a i = 2 j are determined in such a way that the integral square error ε = (5) f(t)φ i (t)dt, (7) M (f(t) a i φ i (t)) 2 dt (8) is minimum Here ε is vanished when M tends to infinity Usually, the series expansion of (5) contains an infinite number of terms for smooth f(t) If f(t) is a piece wise constant or may be approximated as a piecewise constant, then the summation (6) will be terminated after M terms, that is, f(t) M i= i= a i φ i (t) = A T Ψ M (t), (9) where the coefficient vector A = [a, a,, a M ] T and Ψ M (t) = [φ, φ,, φ M ] T Let us define the collocation points t k = (k 5)/M, (k =,, M) With these chosen
3 S ASADI, AH BORZABADI: NUMERICAL SOLUTION OF DELAY 223 collocation points, the function is discretized into a series of nodes with equivalent distances Let the Haar matrix H M M be the combination of Ψ M (t) at all the collocation points Thus, φ (t ) φ (t ) φ (t M ) φ (t ) φ (t ) φ (t M ) H M M = [Ψ M (t ),, Ψ M (t M )] = φ M (t ) φ M (t ) φ M (t M ) M M () For example, H 2 2 = [ Ψ 2 (t ) Ψ 2 (t ) ] [ ] = Therefore, the function f(t) may be approximated as 2 2 f(t k ) c T MH M M, () where c T M = [c, c 2,, c M ] are the wavelet coefficients The integration of the vector Ψ M (t) defined in (9) can be approximated by t Ψ M (t )dt = P M M Ψ M (t), (2) where P M M is M M operational integration matrix which satisfies the following recursive formula [8], [ ] P M M = P M/2 M/2 2M H 2M H M/2 M/2 M/2 M/2 where (M/2) (M/2) is a null matrix of order(m/2) (M/2), P = [ ], (3) 2 22 The delay operational matrix The delay function Ψ M (t τ) is the shift of Ψ M (t) in (9) along the time axis by τ and to estimate them respect to Ψ M (t) the delay operational matrix D(τ) is defined as Ψ M (t τ) D(τ)Ψ M (t), t > τ, t <, (4) where τ [, ] is the delay constant Here is shown how to generate the delay operational matrix D(τ) = [d ij ] for τ 2 with four basis functions These basis functions are given by φ, φ, φ 2, φ 3 By (4) and (4) we have φ (t τ) φ (t) φ (t τ) φ 3 (t τ) = [d ij(τ)] φ (t) φ 3 (t), i, j =, 2, 3, 4, φ 4 (t τ) φ 4 (t) where φ (t τ) =, τ t <, and φ (t τ) = { ifτ t < 2 + τ, if 2 + τ t <, and φ 2 (t τ) = φ 3 (t τ) = { ifτ t < 4 + τ, if 4 + τ t < 2 + τ, { if 2 + τ t < τ, if τ t <,
4 224 TWMS J PURE APPL MATH, V5, N2, 24 To find the entries d ij (τ), i, j =, 2, 3, 4, we use the inner product For example if τ =, we have d =< φ (t τ), φ (t) >= φ (t τ)φ (t)dt = 9, d 3 =< φ 3 (t τ), φ (t) >= φ 3 (t τ)φ (t)dt = If we calculate all d ij (τ) as d and d 3, the 4 4 operational matrix D(τ) will be obtained as D 4 4 () = In a similar manner if we use the vector function φ(t) with dimension 2 n+, then 2 n+ 2 n+ delay matrix D(τ) with τ 2 can be obtained Note that for any dimension if τ = then n matrix is diagonal and we have d il = 4 4 { 2 j if i = l = 2 j + k, φ i (t)φ l (t)dt = ifi l, 3 The collocation approach We discretize the functions φ i (t) in (4) by dividing the interval l [, ], to M equidistance intervals with distance parameter t = /M and introduce the collocation points l k = (k 5)/M, k =,, M, where M is the number of nodes used in the discretization and also is the maximum wavelet index number Also we approximate state variables ẋ(l) by Haar wavelets with M collocation points, ie, ẋ(l) c T x Ψ M (l), (5) where Haar coefficient vector c T x are defined as c T x = [c x,, c xm ] Using the operational integration matrix P M M defined in (2), x(l) = l ẋ(l )dl + x() = l c T x Ψ M (l )dl + x() = c T x P M M Ψ M (l) + x() (6) As stated in (3), the expansion of the matrix Ψ M (l) at the M collocation points will yield the M M Haar matrix H M M and by (6) it can be concluded that ẋ(l k ) = c T x Ψ M (l k ), x(l k ) = c T x P M M Ψ M (l k ) + x(), k =,, M, (7) Now we focus on the analysis of time-delayed systems Choose N i as following manner, N i = Mτ i + 5, (8) i =, 2,, s and let N N 2 N s, (where denotes the bracket function) Using (2), (4) and (), { ξ(l τi ) if l < τ x(l τ i ) = i Cx T i =, 2,, s (9) P M M D(τ i )Ψ M (l) + x() if τ i l < t f By substituting x(l k ), ẋ(l k ) and x(l k τ i ), (i =, 2,, s) in () and using (4)-(9), we have c T x Ψ M (τ k ) = (t f t ) { f(l k, c T x P M M Ψ M (l k ) + x(), ξ(l k τ ), ξ(l k τ 2 ),, ξ(l k τ s )) }, k =,, N, (2),
5 S ASADI, AH BORZABADI: NUMERICAL SOLUTION OF DELAY 225 c T x Ψ M (l k ) = (t f t ) { f ( l k, c T x P M M Ψ M (l k ) + x(), c T x P M M D(τ )Ψ M (l k ) + x(), ξ(l k τ 2 ),, ξ(l k τ s ) )}, k = N +,, N 2, (2) c T x Ψ M (l k ) = (t f t ) { f(l k, c T x P M M Ψ M (l k ) + x(), c T x P M M D(τ )Ψ M (l k ) + x(), c T x P M M D(τ 2 )Ψ M (l k ) + x(),, c T x P M M D(τ s )Ψ M (l k ) + x()) }, k = N s +,, M, (22) In this way, the DDE problem is transformed into nonlinear algebraic equations for the coefficients c T x From the equations (2)-(22), the vector unknown c T x using an iterative approach, eg Newton iterative method, may be evaluated Since the first and last collocation points are not set as the initial and final time, the initial and final variables are calculated according to x() = ξ(), (23) x f = x(m) + ẋ(m)/2m (24) 4 Numerical results In this section, two DDE problems and the results of applying the discussed collocation approach for them is considered Also the obtained approximate solutions have been compared with exact ones Example We consider the equation ẋ(t) = x(t /2) + sin(t /2) + cos(t), < t, (25) ξ(t) = sin(t), t (26) Using the Haar wavelets collocation method with M = 6 collocation points we have τ = 2, N = 8, (27) c T x h k = sin(l k /2) + sin(l k /2) + cos(l k ), k =, 2,, N, (28) c T x h k = c T x P M M D(/2)Ψ M (l k ) x() + sin(l k /2) + cos(l k ), k = N +,, M, (29) and boundary constraints x() = ξ() =, (3) x(t f ) = c T x P M M Ψ M (l M ) + x() + c T x Ψ M (l M )/2M, (3) where c x are all variables which must be obtained at the collocation points The approximate solution with respect to the number of nodes with exact solution of equation are shown in Fig Example 2 In the second example the two-dimensional DDE equation ẋ (t) = x 2 (t), ẋ 2 (t) = x 2 (t) + exp () x 2(t /) + x (t), (32) with delay functions ξ (t) = exp ( t), ξ 2 (t) = exp ( t), for t is considered The approximate solution with using 6 collocation points by applying the collocation approach, in comparison with the exact solution are shown in Figs 2 and 3
6 226 TWMS J PURE APPL MATH, V5, N2, 24 5 Conclusions In this paper an collocation method using the profits of Haar wavelets is presented Applying the given approach is simple and expected that along with the increase of collocation points the approximate solution converges to exact solution of DDE Of course the approximate solution depends on the applied method of solving the nonlinear equation, which will be generated in the process of yielding the mentioned approach and so, like any iterative scheme the nonlinearity of the functional f in () has direct impact on the precision of the approximate solution References [] Baker CTH, Paul CAH, (994), Computing stability regions Runge-Kutta methods for delay differential equations, IMA J Numer Anal, V(4), pp [2] Baker CTH, Paul CAH, (993), Parallel continuous Runge-Kutta methods and vanishing lag delay differential equations, Adv Comp Math, V(), pp [3] Dai R, JECochran, Jr, (29), Wavelet Collocation Method for Optimal control Problems, J Optim Theory Appl, V(43), pp [4] Driver RD, (977), Ordinary and Delay Differential Equations, Applied Mathematical Sciences, Springer- Verlag, New York-Heidelberg, 2p [5] El-Hawary HM, Mahmoud SM, (23), Spline collocation methods for solving delay-differential equations, Appl Math Comput, V(46), pp [6] Engelborghs K, Luzyanina T, Int l Hout KJ, Roose, D, (2), Collocation methods for the computation of periodic solutions of delay differential equations, SIAM J Sci Comput, V(22), pp [7] Gopalsamy K, (992), Stability and Oscillations in Delay Differential Equations of Population Dynamics, Mathematics and its Applications, Kluwer Academic Publishers Group, Dordrecht, 74p [8] Gu JS, Jian WS, (996), The Haar Wavelet operation matrix of integration, Int J Syst Sci, V(27)(7), pp [9] Houwen PJ, Sommelier BP, (984), Stability in linear multistep methods for pure delay equations, J Comp Appl Math, V(), pp55-63 [] Kolmanovskii VB, Nosov VR, (986), Stability of Functional-Differential Equations, Mathematics in Science and Engineering, Academic Press, Inc [Harcourt Brace Jovanovich, Publishers], London, 8 p [] Koto T, (23), Stability of θ-methods for delay integro-differential equations, J Comput Appl Math, V(6), pp [2] Kuang K, (993), Delay Differential Equations with Applications in Population Dynamics, Mathematics in Science and Engineering, Academic Press, Inc, Boston, MA, 9p [3] Oberle HJ, Pesch HJ, (98), Numerical treatment of delay differential equations by Hermite interpolation, Numer Math, V(37), pp [4] Oppelstrup, J, (978), The RKFHB4 method for delay differential equations, Lecture Notes in Mathematics, V(63), pp33-46
7 S ASADI, AH BORZABADI: NUMERICAL SOLUTION OF DELAY Haar wavelet with 6 nodes Haar wavelet with 32 nodes 7 Numeric solution for x time Figure Approximate and exact solution for Example 4 2 Haar wavelet with 6 nodes exact solution 8 6 x Time Figure 2 Approximate and exact solution x for Example x Haar wavelet with 6 nodes exact solution Time Figure 3 Approximate and exact solution x 2 for Example 4
8 228 TWMS J PURE APPL MATH, V5, N2, 24 Solayman Asadi was born in Sari, Iran, 987 He received his BSc in pure mathematics in Ferdowsi University of Mashhad in 2, MSc in applied mathematics in 22 from Damghan University in IRAN His research areas include optimal control, optimization and analysis Akbar Hashemi Borzabadi received his BSc from Birjand University, Birjand, Iran and MSc and PhD from Ferdowsi University of Mashhad, Mashhad, Iran Now he is an Associate Professor at the School of Mathematics and Computer Science, Damghan University, Iran and his research areas are optimal control, optimization, approximation theory and numerical analysis
Numerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationSolving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted Legendre Functions
International Journal of Mathematical Modelling & Computations Vol. 5, No. 3, Summer 215, 219-23 Solving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More informationModified Milne Simpson Method for Solving Differential Equations
Modified Milne Simpson Method for Solving Differential Equations R. K. Devkate 1, R. M. Dhaigude 2 Research Student, Department of Mathematics, Dr. Babasaheb Ambedkar Marathwada University, Aurangabad,
More informationResearch Article Hermite Wavelet Method for Fractional Delay Differential Equations
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
More informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
More informationPAijpam.eu NUMERICAL SOLUTION OF WAVE EQUATION USING HAAR WAVELET Inderdeep Singh 1, Sangeeta Arora 2, Sheo Kumar 3
International Journal of Pure and Applied Mathematics Volume 98 No. 4 25, 457-469 ISSN: 3-88 (printed version); ISSN: 34-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/.2732/ijpam.v98i4.4
More informationON THE SOLUTION OF DIFFERENTIAL EQUATIONS WITH DELAYED AND ADVANCED ARGUMENTS
2003 Colloquium on Differential Equations and Applications, Maracaibo, Venezuela. Electronic Journal of Differential Equations, Conference 13, 2005, pp. 57 63. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationSuperconvergence analysis of multistep collocation method for delay Volterra integral equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 216, pp. 25-216 Superconvergence analysis of multistep collocation method for delay Volterra integral equations
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationNumerical solution of linear time delay systems using Chebyshev-tau spectral method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp. 445 469 Numerical solution of linear time
More informationProperties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation
Applied Mathematical Sciences, Vol. 6, 212, no. 32, 1563-1569 Properties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation Ahmad Shahsavaran 1 and Abar Shahsavaran
More informationSOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 3(2) (214), 281 288 Palestine Polytechnic University-PPU 214 SOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX Naresh Berwal, Dinesh Panchal and
More informationAssignment #09 - Solution Manual
Assignment #09 - Solution Manual 1. Choose the correct statements about representation of a continuous signal using Haar wavelets. 1.5 points The signal is approximated using sin and cos functions. The
More informationA NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS. Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan Przeradzki
Opuscula Math. 33, no. 1 213, 5 17 http://dx.doi.org/1.7494/opmath.213.33.1.5 Opuscula Mathematica A NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan
More informationA convolution test equation for double delay integral equations
Journal of Computational and Applied Mathematics 228 (2009) 589 599 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationConverse Lyapunov-Krasovskii Theorems for Systems Described by Neutral Functional Differential Equation in Hale s Form
Converse Lyapunov-Krasovskii Theorems for Systems Described by Neutral Functional Differential Equation in Hale s Form arxiv:1206.3504v1 [math.ds] 15 Jun 2012 P. Pepe I. Karafyllis Abstract In this paper
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationAsymptotic stability of solutions of a class of neutral differential equations with multiple deviating arguments
Bull. Math. Soc. Sci. Math. Roumanie Tome 57(15) No. 1, 14, 11 13 Asymptotic stability of solutions of a class of neutral differential equations with multiple deviating arguments by Cemil Tunç Abstract
More informationRESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials
Journal of Modern Methods in Numerical Mathematics ISSN: 29-477 online Vol, No, 2, c 2 Modern Science Publishers wwwm-sciencescom RESEARCH ARTICLE Spectral Method for Solving The General Form Linear Fredholm
More informationSolving Delay Differential Equations By Elzaki Transform Method
SCITECH Volume 3, Issue 1 RESEARCH ORGANISATION February 22, 2017 Boson Journal of Modern Physics www.scitecresearch.com Solving Delay Differential Equations By Elzaki Transform Method Ebimene, J. Mamadu
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
More informationGeneralized B-spline functions method for solving optimal control problems
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 2, No. 4, 24, pp. 243-255 Generalized B-spline functions method for solving optimal control problems Yousef Edrisi Tabriz
More informationNumerical Solution of Volterra-Fredholm Integro-Differential Equation by Block Pulse Functions and Operational Matrices
Gen. Math. Notes, Vol. 4, No. 2, June 211, pp. 37-48 ISSN 2219-7184; Copyright c ICSRS Publication, 211 www.i-csrs.org Available free online at http://www.gean.in Nuerical Solution of Volterra-Fredhol
More informationSolving Delay Differential Equations (DDEs) using Nakashima s 2 Stages 4 th Order Pseudo-Runge-Kutta Method
World Applied Sciences Journal (Special Issue of Applied Math): 8-86, 3 ISSN 88-495; IDOSI Publications, 3 DOI:.589/idosi.wasj.3..am.43 Solving Delay Differential Equations (DDEs) using Naashima s Stages
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationALMOST PERIODIC SOLUTIONS OF NONLINEAR DISCRETE VOLTERRA EQUATIONS WITH UNBOUNDED DELAY. 1. Almost periodic sequences and difference equations
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 10, Number 2, 2008, pages 27 32 2008 International Workshop on Dynamical Systems and Related Topics c 2008 ICMS in
More informationTHE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands
THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay
More informationPERSISTENCE AND PERMANENCE OF DELAY DIFFERENTIAL EQUATIONS IN BIOMATHEMATICS
PERSISTENCE AND PERMANENCE OF DELAY DIFFERENTIAL EQUATIONS IN BIOMATHEMATICS PhD Thesis by Nahed Abdelfattah Mohamady Abdelfattah Supervisors: Prof. István Győri Prof. Ferenc Hartung UNIVERSITY OF PANNONIA
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationTHE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS
THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. BY DR. LOTHAR COLLATZ
More informationarxiv: v2 [math.fa] 27 Sep 2016
Decomposition of Integral Self-Affine Multi-Tiles Xiaoye Fu and Jean-Pierre Gabardo arxiv:43.335v2 [math.fa] 27 Sep 26 Abstract. In this paper we propose a method to decompose an integral self-affine Z
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationNumerical Analysis: Approximation of Functions
Numerical Analysis: Approximation of Functions Mirko Navara http://cmp.felk.cvut.cz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office 104a
More informationSOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD
SOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD B. SEPEHRIAN AND M. RAZZAGHI Received 5 November 24 Single-term Walsh series are developed to approximate
More informationApplication of Semiorthogonal B-Spline Wavelets for the Solutions of Linear Second Kind Fredholm Integral Equations
Appl Math Inf Sci 8, No, 79-84 (4) 79 Applied Mathematics & Information Sciences An International Journal http://dxdoiorg/78/amis/8 Application of Semiorthogonal B-Spline Wavelets for the Solutions of
More informationNumerical Investigation of the Time Invariant Optimal Control of Singular Systems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol. 8, 24, no. 2, 6-68 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ams.24.4863 Numerical Investigation of the Time Invariant Optimal Control of Singular Systems
More informationNumerical solution of delay differential equations via operational matrices of hybrid of block-pulse functions and Bernstein polynomials
Computational Methods for Differential Equations http://cmdetabrizuacir Vol, No, 3, pp 78-95 Numerical solution of delay differential equations via operational matrices of hybrid of bloc-pulse functions
More information(Received 10 December 2011, accepted 15 February 2012) x x 2 B(x) u, (1.2) A(x) +
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol13(212) No4,pp387-395 Numerical Solution of Fokker-Planck Equation Using the Flatlet Oblique Multiwavelets Mir Vahid
More informationBLOCK-PULSE FUNCTIONS AND THEIR APPLICATIONS TO SOLVING SYSTEMS OF HIGHER-ORDER NONLINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 214 (214), No. 54, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu BLOCK-PULSE FUNCTIONS
More informationNumerical solution of optimal control problems by using a new second kind Chebyshev wavelet
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 2, 2016, pp. 162-169 Numerical solution of optimal control problems by using a new second kind Chebyshev wavelet Mehdi
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationSolving Linear Time Varying Systems by Orthonormal Bernstein Polynomials
Science Journal of Applied Mathematics and Statistics 2015; 3(4): 194-198 Published online July 27, 2015 (http://www.sciencepublishinggroup.com/j/sjams) doi: 10.11648/j.sjams.20150304.15 ISSN: 2376-9491
More informationNUMERICAL TREATMENTS FOR VOLTERRA DELAY INTEGRO-DIFFERENTIAL EQUATIONS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol92009, No3, pp292 308 c 2009 Institute of Mathematics of the National Academy of Sciences of Belarus NUMERICAL TREATMENTS FOR VOLTERRA DELAY INTEGRO-DIFFERENTIAL
More informationASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS. Janusz Migda and Małgorzata Migda
Opuscula Math. 36, no. 2 (2016), 265 278 http://dx.doi.org/10.7494/opmath.2016.36.2.265 Opuscula Mathematica ASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS Janusz Migda and Małgorzata
More informationSteepest descent method on a Riemannian manifold: the convex case
Steepest descent method on a Riemannian manifold: the convex case Julien Munier Abstract. In this paper we are interested in the asymptotic behavior of the trajectories of the famous steepest descent evolution
More informationIdentification of Parameters in Neutral Functional Differential Equations with State-Dependent Delays
To appear in the proceedings of 44th IEEE Conference on Decision and Control and European Control Conference ECC 5, Seville, Spain. -5 December 5. Identification of Parameters in Neutral Functional Differential
More informationStability of Adjointable Mappings in Hilbert
Stability of Adjointable Mappings in Hilbert arxiv:math/0501139v2 [math.fa] 1 Aug 2005 C -Modules M. S. Moslehian Abstract The generalized Hyers Ulam Rassias stability of adjointable mappings on Hilbert
More informationAn Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley Equation
Applied Mathematical Sciences, Vol. 11, 2017, no. 30, 1467-1479 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.7141 An Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley
More informationOperational Matrices for Solving Burgers' Equation by Using Block-Pulse Functions with Error Analysis
Australian Journal of Basic and Applied Sciences, 5(12): 602-609, 2011 ISSN 1991-8178 Operational Matrices for Solving Burgers' Equation by Using Block-Pulse Functions with Error Analysis 1 K. Maleknejad,
More informationA METHOD FOR SOLVING DUAL FUZZY GENERAL LINEAR SYSTEMS*
Appl Comput Math 7 (2008) no2 pp235-241 A METHOD FOR SOLVING DUAL FUZZY GENERAL LINEAR SYSTEMS* REZA EZZATI Abstract In this paper the main aim is to develop a method for solving an arbitrary m n dual
More informationThe exponential asymptotic stability of singularly perturbed delay differential equations with a bounded lag
J. Math. Anal. Appl. 270 (2002) 143 149 www.academicpress.com The exponential asymptotic stability of singularly perturbed delay differential equations with a bounded lag Hongjiong Tian Department of Mathematics,
More informationLecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms
ECE785: Hybrid Systems:Theory and Applications Lecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms Wei Zhang Assistant Professor Department of Electrical
More informationSemiorthogonal Quadratic B-Spline Wavelet Approximation for Integral Equations
Mathematical Sciences Vol. 3, No. (29) 99- Semiorthogonal Quadratic B-Spline Wavelet Approximation for Integral Equations Mohsen Rabbani a,, Nasser Aghazadeh b a Department of Mathematics, Islamic Azad
More informationDecentralized control synthesis for bilinear systems using orthogonal functions
Cent Eur J Eng 41 214 47-53 DOI: 12478/s13531-13-146-1 Central European Journal of Engineering Decentralized control synthesis for bilinear systems using orthogonal functions Research Article Mohamed Sadok
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationLecture 2: Haar Multiresolution analysis
WAVELES AND MULIRAE DIGIAL SIGNAL PROCESSING Lecture 2: Haar Multiresolution analysis Prof.V. M. Gadre, EE, II Bombay 1 Introduction HAAR was a mathematician, who has given an idea that any continuous
More informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationLOCALLY POSITIVE NONLINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci. 3 Vol. 3 No. 4 55 59 LOCALLY POSITIVE NONLINEAR SYSTEMS TADEUSZ KACZOREK Institute of Control Industrial Electronics Warsaw University of Technology ul. Koszykowa 75 66
More informationDirect method for variational problems by using hybrid of block-pulse and Bernoulli polynomials
Direct method for variational problems by using hybrid of block-pulse and Bernoulli polynomials M Razzaghi, Y Ordokhani, N Haddadi Abstract In this paper, a numerical method for solving variational problems
More informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
More informationHomotopy method for solving fuzzy nonlinear equations
Homotopy method for solving fuzzy nonlinear equations S. Abbasbandy and R. Ezzati Abstract. In this paper, we introduce the numerical solution for a fuzzy nonlinear systems by homotopy method. The fuzzy
More informationOSCILLATION AND ASYMPTOTIC STABILITY OF A DELAY DIFFERENTIAL EQUATION WITH RICHARD S NONLINEARITY
2004 Conference on Diff. Eqns. and Appl. in Math. Biology, Nanaimo, BC, Canada. Electronic Journal of Differential Equations, Conference 12, 2005, pp. 21 27. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT. A. Khani
Acta Universitatis Apulensis ISSN: 1582-5329 No 38/214 pp 11-22 ADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT A Khani Abstract In this
More informationTime-dependent order and distribution policies in supply networks
Time-dependent order and distribution policies in supply networks S. Göttlich 1, M. Herty 2, and Ch. Ringhofer 3 1 Department of Mathematics, TU Kaiserslautern, Postfach 349, 67653 Kaiserslautern, Germany
More informationExponentially Fitted Error Correction Methods for Solving Initial Value Problems
KYUNGPOOK Math. J. 52(2012), 167-177 http://dx.doi.org/10.5666/kmj.2012.52.2.167 Exponentially Fitted Error Correction Methods for Solving Initial Value Problems Sangdong Kim and Philsu Kim Department
More informationNumerical bifurcation analysis of delay differential equations
Numerical bifurcation analysis of delay differential equations Dirk Roose Dept. of Computer Science K.U.Leuven Dirk.Roose@cs.kuleuven.be Acknowledgements Many thanks to Koen Engelborghs Tatyana Luzyanina
More informationA Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations
Acta Appl Math (21) 19: 861 873 DOI 1.17/s144-8-9351-y A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations M. Rasty M. Hadizadeh Received:
More informationMLISP: Machine Learning in Signal Processing Spring Lecture 8-9 May 4-7
MLISP: Machine Learning in Signal Processing Spring 2018 Prof. Veniamin Morgenshtern Lecture 8-9 May 4-7 Scribe: Mohamed Solomon Agenda 1. Wavelets: beyond smoothness 2. A problem with Fourier transform
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationDisconjugate operators and related differential equations
Disconjugate operators and related differential equations Mariella Cecchi, Zuzana Došlá and Mauro Marini Dedicated to J. Vosmanský on occasion of his 65 th birthday Abstract: There are studied asymptotic
More informationApplied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations
Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationApplication of modified Leja sequences to polynomial interpolation
Special Issue for the years of the Padua points, Volume 8 5 Pages 66 74 Application of modified Leja sequences to polynomial interpolation L. P. Bos a M. Caliari a Abstract We discuss several applications
More informationTwo Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. (8) 7-77 ISSN - Two Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation T.V. Triantafyllidis,
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics ON A HYBRID FAMILY OF SUMMATION INTEGRAL TYPE OPERATORS VIJAY GUPTA AND ESRA ERKUŞ School of Applied Sciences Netaji Subhas Institute of Technology
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationNEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS
NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS Nikolaos D. Atreas Department of Mathematics, Aristotle University of Thessaloniki, 54006, Greece, e-mail:natreas@auth.gr Abstract We
More informationLegendre Wavelets Based Approximation Method for Cauchy Problems
Applied Mathematical Sciences, Vol. 6, 212, no. 126, 6281-6286 Legendre Wavelets Based Approximation Method for Cauchy Problems S.G. Venkatesh a corresponding author venkamaths@gmail.com S.K. Ayyaswamy
More informationNumerical Solution of Fredholm Integro-differential Equations By Using Hybrid Function Operational Matrix of Differentiation
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISS 8-56) Vol. 9, o. 4, 7 Article ID IJIM-8, pages Research Article umerical Solution of Fredholm Integro-differential Equations
More informationNon-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems
Non-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems Jalil Rashidinia Reza Jalilian Abstract In this paper we use quintic non-polynomial spline functions to develop numerical
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationMatrix power converters: spectra and stability
Matrix power converters: spectra and stability Stephen Cox School of Mathematical Sciences, University of Nottingham supported by EPSRC grant number EP/E018580/1 Making It Real Seminar, Bristol 2009 Stephen
More informationSHADOWING PROPERTY FOR INDUCED SET-VALUED DYNAMICAL SYSTEMS OF SOME EXPANSIVE MAPS
Dynamic Systems and Applications 19 (2010) 405-414 SHADOWING PROPERTY FOR INDUCED SET-VALUED DYNAMICAL SYSTEMS OF SOME EXPANSIVE MAPS YUHU WU 1,2 AND XIAOPING XUE 1 1 Department of Mathematics, Harbin
More informationDelay Differential Equations Part II: Time and State dependent Lags
Delay Differential Equations Part II: Time and State dependent Lags L.F. Shampine Department of Mathematics Southern Methodist University Dallas, Texas 75275 shampine@smu.edu www.faculty.smu.edu/shampine
More informationThe use of cubic spline and far-side boundary condition for the collocation solution of a transient convection-diffusion problem
Korean J Chem Eng, 24(2), 204-208 (2007) SHORT COMMUNICATION The use of cubic spline and far-side boundary condition for the collocation solution of a transient convection-diffusion problem Woohyun Yun
More information