Numerical solution of linear time delay systems using Chebyshev-tau spectral method
|
|
- Bertram Hunt
- 5 years ago
- Views:
Transcription
1 Available at Appl. Appl. Math. ISSN: Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp Numerical solution of linear time delay systems using Chebyshev-tau spectral method M. Mousa-Abadian, S. H. Momeni-Masuleh Corresponding Author Department of Mathematics Shahed University P.O. Box: Tehran, Iran Received: May 3, 06; Accepted: January 0, 07 Abstract In this paper, a hybrid method based on method of steps and a Chebyshev-tau spectral method for solving linear time delay systems of differential equations is proposed. The method first converts the time delay system to a system of ordinary differential equations by the method of steps and then employs Chebyshev polynomials to construct an approximate solution for the system. In fact, the solution of the system is expanded in terms of orthogonal Chebyshev polynomials which reduces the solution of the system to the solution of a system of algebraic equations. Also, we transform the coefficient matrix of the algebraic system to a block quasi upper triangular matrix and the latter system can be solved more efficiently than the first one. Furthermore, using orthogonal Chebyshev polynomials enables us to apply fast Fourier transform for calculating matrix-vector multiplications which makes the proposed method to be more efficient. Consistency, stability and convergence analysis of the method are provided. Numerous numerical examples are given to demonstrate efficiency and accuracy of the method. Comparisons are made with available literature. Keywords: Linear delay systems; commensurate delays; method of steps; Chebyshev-tau spectral method; stability analysis; convergence analysis MSC 00 No.: 65L03, 65L06, 65D99 445
2 446 M. Mousa-Abadian & S. H. Momeni-Masuleh. Introduction Delay differential equations (DDEs) have been applied to model real phenomena. These equations often appear in many real-world problems. They occur in a wide variety of physical, chemical, engineering, economic and biological systems and their networks. One can cite many examples where delay plays an important role. As an example, in biological mathematics that delay plays a crucial role is drug therapy for the disease of human immunodeficiency virus (HIV) infection (Kirschner (996), Nelson and Perelson (00)). These equations also known as difference-differential equations, are a special class of differential equations called functional differential equations (Hale and Lunel (03), Elsgolts and Norkin (973)). Delay differential equations were initially introduced in the 8th century by Laplace and Condorcet (Gorecki et al. (989)). Numerically approximating a solution to a DDE has been studied by many authors (see, for example, Bellen and Zennaro (03), Bellman and Cooke (965) and Feldstein and Neves (984)). Currently, most DDE solver are generated using adapted continuous numerical methods for solving ODE. To obtain a numerical solution of a DDE by means of a discrete method often we need to know the solution on a set of points that differ from the grid. To calculate the approximate solutions on non-grid points, a local interpolation is needed. In such situation, stability of discrete method may be lost. Bellen and Zennaro (988), introduced the concept of stable interpolant for Runge-Kutta methods. In other words, a stable interpolant maintains the stability properties of the Runge-Kutta method itself. Hermite interpolation can be also used to solve a DDE numerically. Oberle and Pesch (98), combined Runge- Kutta methods and Hermite interpolation for the numerical solution of a DDE. In fact, they have been shown that the combination has high order of accuracy and it is reliable. Nevertheless, the most popular discrete methods are Runge-Kutta methods that use interpolation to obtain the numerical solution of a DDE. Wherever the step size must be reduced often to produce answers at specified points, a Runge-Kutta formula becomes inefficient. Bogacki and Shampine (990), have been presented a structural way to handle this issue. When a high order method for solving a DDE is applied, the amount of computations for the interpolation process increase considerably. To overcome this difficulty, Zennaro (985a), presented a one-step collocation method. Also, for super convergence rate, Zennaro (985b), considered one-step collocation method at Gaussian points. Wang and Wang (00), proposed the single-step and multiple-domain Legendre-Gauss collocation integration processes for nonlinear DDEs. One of the most important advantages of their method is its accuracy and efficiency. Another class of numerical methods that can be considered are Spline methods. Spline collocation methods are used by Banks and Kappel (979), Kappel and Kunisch (98) and Kemper (975). So far, all of the mentioned methods for solving a DDE have finite order of accuracy. To get infinite order of accuracy for a smooth problem usually spectral methods are applied. For more details about these methods, see for example Gottlieb and Orszag (977) and Canuto et al. (987). Ito et al. (99), have been employed shifted Legendre polynomials to construct a spectral approximate solution for a DDE. More precisely,
3 AAM: Intern. J., Vol., Issue (June 07) 447 they considered Legendre-tau method to construct an approximate solution for a DDE with one constant delay. Sedaghat et al. (0) proposed a numerical scheme to solve delay differential equations of pantograph type. Their method consists of expanding the required approximate solution as the elements of the shifted Chebyshev polynomials. Aziz and Amin (06) applied Haar wavelet collocation method to obtain the numerical solution of DDEs. Behroozifar and Yousefi (03) introduced hybrid of block-pulse functions and Bernstein polynomials to construct an approximate solution for DDEs. Ghasemi and Tavassoli Kajani (0) employed Chebyshev wavelets to obtain a numerical scheme for solving DDEs. They introduced operational matrix of delay and utilize it to reduce the solution of time-varying delay systems to the solution of algebraic equations. In this paper, we propose a method for solving linear time delay systems of differential equations. At the first step, the proposed method uses the method of steps (Bellen and Zennaro (03)) to convert the linear system of DDE on a given interval to an ODE over that interval, by employing the known history function, say φ(t), for that interval. Then, the approximate solution for this system is obtained by tau-galerkin method (Gottlieb and Orszag (977), Canuto et al. (987)). We construct an approximate solution by Chebyshev polynomials expansion (Fox and Parker (968)). More precisely, we consider the solution of the ODE system as a truncated series of Chebyshev polynomials. This leads us to solve a linear system of algebraic equations. Furthermore, we transform the coefficient matrix of the algebraic system to a block quasi upper triangular form. Therefore, the algebraic system can be solved more efficiently than the original one. In addition, properties of the Chebyshev polynomials, enriches the proposed method to apply fast Fourier transform for computing matrix-vector multiplications. At the next step, the newly found approximate solution is used as the history function for the next interval and the above process is carried out iteratively to obtain the approximate solution for the whole time interval. Also, the consistency, stability and convergence of the proposed method are provided. Numerical results are presented to demonstrate the accuracy and efficiency of the proposed method. In the following subsection we explain notations, preliminaries and the problem statement... Notations, preliminaries and problem statement We adopt the following notations throughout this paper. R n is Euclidean space with corresponding norm, R n n is space of real square matrices, L ([, ]; R n ) or for short L (, ) is the Hilbert space of square integrable functions that map [, ] to a member of R n with inner product (f, g) w = f(x)g(x)w(x)dx, where
4 448 M. Mousa-Abadian & S. H. Momeni-Masuleh w(x) = ( x ) / (that for simplicity we may drop the subscript w), and corresponding norm L w, H ([, ]; R n ) or for convenience H (, ) is the Sobolov space of absolutely continuous functions with square integrable derivative with norm,w. Let Ω be an open bounded domain in R, with piecewise smooth boundary Γ. We assume that we want to approximate the boundary-value problem Lx = f in Ω, () Bx = 0 on Γ, () where L is a linear differential operator in Ω, and B is a set of linear boundary differential operators on Γ. Let X to be a space of real functions defined in Ω that are square integrable with respect to w(x). The domain of definition of L, i.e., the subset D(L) of those functions x of X for which Lx is still an element of X, is supposed to be a dense subspace of X. Thus, L is a linear operator from D(L) to X. Prescribing the boundary conditions () amounts to restricting the domain of L to the subspace D B (L) of D(L) defined by D B (L) = {x D(L) Bx = 0 on Γ}. Suppose that X N and Y N are finite-dimensional subspaces of X having the same dimension. In the tau method, we have x N X N, (Lx N, v) = (f, v) v Y N. Let P N be the orthogonal projection of L w onto the space spanned by {p k } N k=0, in which p k is a polynomial of degree k and we denote this set by P N. Lemma.. Canuto et al. (987) Suppose that N P N x(t) = â k T k (t), k=0 where T k (t) is the Chebyshev polynomial of degree k. Then, for all x(t) H (, ) we have x(t) P N x(t) L w CN x(t),w, (3)
5 AAM: Intern. J., Vol., Issue (June 07) 449 where the positive constant C is independent of the function x, the integer N, and the diameter of the domain. Theorem.. Lax-Richtmyer equivalence theorem, Gottlieb and Orszag (977) A consistent approximation to a well-posed linear problem is stable if and only if it is convergent. Lemma.3. Babu ska inf-sup condition, Canuto et al. (987) Let X N F and Y N G for all positive N. If there exists a constant α > 0 independent of N such that then the following estimate holds (Lx, v) α x F sup x X N, v Y N v G v 0 x N F C α f, where the positive constant C is independent of N. One important class of delay differential equations is linear time invariant systems (LTI) or briefly time-delayed LTI systems. General representation of a time-delayed LTI systems can be given by ẋ(t) = A 0 x(t) + N k= A kx(t τ k ) + u(t), x(t) = φ(t), t [ τ, 0), x(0) = x 0, (4) where 0 τ τ τ N = τ are positive delays, x(t), u(t) R n are the state and input vectors respectively, A k (k = 0,,,, N) R n n are system matrices which are given and φ L ([ τ, 0]; R n ). Theorem.4. Delfour and Mitter (975) Let (x 0, φ) M = R n L ([ τ, 0]; R n ). Then, the system (4) for any T 0 admits a unique solution x(t) L ([ τ, T ]; R n ) H ([ τ, T ]; R n ), which depends continuously on the initial data (x 0, φ) M.
6 450 M. Mousa-Abadian & S. H. Momeni-Masuleh In this paper, we consider the system (4) with two discrete commensurate delays. However the method can be applied to systems with more than two commensurate delays. It can be also extended to solve nonlinear time delay problems. The remaining sections of this paper is organized as follows. In Section, we recall some properties of Chebyshev polynomials, while Section 3 describes the Chebyshev-tau spectral approach for solving the system (4) with two commensurate delays. Section 4 provides the stability and convergence analysis of the method while Section 5 presents illustrative examples. Finally, Section 6 contains some concluding remarks.. Properties of Chebyshev polynomials Here, we review some properties of Chebyshev polynomials (Fox and Parker (968)). The Chebyshev polynomials of the first kind, T k (x), k = 0,,..., x, can be defined as the solution of the following ordinary differential equation d ( x dt ) k dx dx (x) k + T x k(x) = 0, (5) which satisfy T k () =. Thus, T 0 (x) =, T (x) = x, and so on. Chebyshev polynomials {T k (x)} k=0 satisfy the orthogonality relation (T k, T l ) w = 0, for k l, (6) and they form a basis for L (, ). More precisely, any u L (, ) can be written as u = where û k T k, k=0 and û k = u(x)t k (x)w(x)dx, k 0, πc k c k = {, k = 0,, k, (7) with
7 AAM: Intern. J., Vol., Issue (June 07) 45 u L w = k=0 c k π û k. These polynomials also satisfy the recursion relation If u is represented as T k+ (x) = xt k (x) T k (x). (8) u = then N û k T k (x), k=0 where u = du N dx = û () k T k(x), (9) k=0 û () k = c k N j=k+ j+k odd jû j, k = 0,,, N. (0) One important property that we will use it for efficient implementation of the method is the following relation (Canuto et al. (987)) From (), we have T k (x) = k + T k+ (x) k T k (x), k. () and so kû k = c k û () k û() k+, k, c k û () k = û () k+ + (k + )û k+, 0 k N. () The generalization of () is given by (Canuto et al. (987))
8 45 M. Mousa-Abadian & S. H. Momeni-Masuleh c k û (q) k = û (q) k+ + (k + )û(q ) k+, k 0. (3) 3. Chebyshev-tau approximation Now, we focus on the Chebyshev approximation of solutions to the system (4) with two discrete point delays which are commensurate. The proposed method uses the method of steps (Bellen and Zennaro (03)) and Chebyshev-tau spectral method (Gottlieb and Orszag (977), Canuto et al. (987)). We begin our discussion by considering the following system ẋ(t) = A 0 x(t) + A x(t τ ) + A x(t τ ) + u(t) x(t) = φ(t), t [ τ, 0) x(0) = x 0, (4) where 0 < τ τ, τ τ = k N (N stands for natural number set). It is clear that for t [0, τ ], the linear delay system (4) is equivalent to the following initial value problem { ẋ(t) = A 0 x(t) + A φ(t τ ) + A φ(t τ ) + u(t), t [0, τ ], x(0) = x 0. (5) It is assumed that the approximate solution x N (t) to be expanded in terms of Chebyshev polynomials, i.e., x N (t) = N k=0 â k T k ( t τ τ ), t [0, τ ], (6) where T k (θ) are the Chebyshev polynomials, â k R n, k = 0,,,..., N, are the unknowns to be computed. On the other hand, we also expand φ(t τ ), φ(t τ ) and u in terms of a truncated Chebyshev series as follows where φ N (t τ ) = φ N (t τ ) = u N (t) = N k=0 N k=0 N k=0 ˆφ () k T k( t τ τ ), (7) ˆφ () k T k( t τ τ ), (8) û k T k ( t τ τ ), t [0, τ ], (9)
9 AAM: Intern. J., Vol., Issue (June 07) 453 ˆφ () k = πc k ˆφ () k = πc k and the c k are given by (7). û k = πc k φ ( τ (s )) T k (s) ds, s φ ( τ (s + ) τ ) Tk (s) ds, s u ( τ (s + )) T k (s) ds, s Now, by using the method of weighted residuals (Canuto et al. (987)), we obtain (ẋn (t) A 0 x N (t) A φ N (t τ ) A φ N (t τ ) u N (t), T k ( t τ ) ) = 0, (0) τ ω for 0 k N. From (6), (7), (8), (9), the orthogonality of the Chebyshev polynomials and ẋ N (t) = N τ k=0 b k T k ( t τ τ ), where b k are defined to be the same as û () k in (0), we obtain b k = A 0 â k + A τ ˆφ() k + A ˆφ() k + û k, 0 k N. () Linear algebraic system () has N + unknowns but only N equations. To obtain a unique solution for the system, we need an additional equation that deals with the initial condition. We note that the important difference between the tau approximation (Gottlieb and Orszag (977), Canuto et al. (987)), and the Galerkin approximation is that the basis functions T k (θ) are not required to satisfy the initial condition in (5). The additional equation is x N (0) = N â k T k ( ) = k=0 N â k ( ) k = x 0. () Rewriting () and () in matrix form, leads to the following algebraic system for the N + unknowns â k k=0 where V N a N = W N a N + b N, t [0, τ ], (3)
10 454 M. Mousa-Abadian & S. H. Momeni-Masuleh a N = (â 0, â,..., â N ) T, b N = (A ˆφ() 0 + A ˆφ() 0 + û 0, A ˆφ() + A ˆφ() + û,..., A ˆφ() N + A () ˆφ N + û N, x 0 ) T, and W N = diag(a 0, A 0,..., A 0, 0). For N even (for N odd, only the last column of V N is different), V N is given by V N = τ N N (N ) I, (4) (N ) N τ τ τ τ τ where I is n n identity matrix and denotes the Kronecker product. We can repeat the above procedure to solve (4) on [τ, τ ], [τ, 3τ ], and so on. We note that for each time interval [iτ, (i + )τ ], i = 0,,, only the right hand side of Equation (3) changes, but the coefficient matrix V N W N remains unchange. Now, we discuss some implementation issues. First, computation of the right hand side of () () Equation (3) needs to compute the Chebyshev coefficients ˆφ k, ˆφ k of the initial function and Chebyshev coefficients û k of the forcing function on each time interval [iτ, (i + )τ ], i = 0,,. In fact we use the matrix C jk = cos πjk N, k =,,..., N, to transform Chebyshev space to physical space and inverse transform, i.e., transform from physical space to Chebyshev space that done by where c j = (C ) jk = N c j c k cos πjk N, {, j = 0, N,, j =,..., N. Both transforms may be evaluated by the fast Fourier transform (Canuto et al. (987)). Second, the linear system of Equation (3) can be reduced to a block quasi upper triangular form with two off-diagonal band. If ordinary Gauss elimination were used to solve the linear system of Equations (3) with n(n + ) equations and unknowns, it would be require about (n(n + )) 3 multiplications for the L N U N decomposition of V N W N.
11 AAM: Intern. J., Vol., Issue (June 07) 455 However, by taking into consideration the spectral structures of matrix V N and W N the LU-factorization can be carried out in only c n(n + ) multiplications and about c n(n + ) multiplications for each right hand side, where c < c and c j N, j =,. From (), we have If we invoke the recursion relation(3) with q =, then b k = τ ( ) A0 â k + A ˆφ() k + A ˆφ() k + û k. (5) Now, from (5) we obtain â k+ = ( ) ck b k b k+. (6) (k + ) â k+ = ( (τ c k (k + ) τ By simplifying, for 0 k N we obtain + A ˆφ() k + û k ) ) (A 0â k + A ˆφ() k ( A0 â k+ + A ˆφ() k+ + A () ˆφ k+ + û k+) ). (7) If we put (k + )â k+ τ c k A 0â k + τ A 0â k+ = τ c k A () ˆφ k + τ c k A () ˆφ k g k = τ + τ c k ûk τ A ˆφ () k+ τ A ( ) c k A ˆφ() k + c k A ˆφ() k + c k û k û k+ A ˆφ() k+ A () ˆφ k+, we arrive at a linear system of equations with the following coefficient matrix ˆφ () k+ τ ûk+. τ τ A 0 A τ A τ 0 4 A τ A N A τ 0 6 A = , (8) τ A 0 (N ) τ A τ A 0 N τ A 0 τ I τ I τ I τ I τ I for the right hand side g = (g 0, g,..., g N, x 0 ) T and the unknown vector a N = (â 0, â,..., â N ) T. Therefore, the following linear system of equations should be solved
12 456 M. Mousa-Abadian & S. H. Momeni-Masuleh A N a N = g. (9) The linear system of Equations (3) and (9) are equivalent. However, the system matrix in (9) is of block quasi upper triangular form and therefore, Equation (9) can be solved more efficiently than Equation (3). One way is to apply Gaussian elimination procedure to a block at a time (Isaacson and Keller (996)). An operation count reveals that block elimination for a block quasi upper triangular matrix of the form (9) requires about n 3 N multiplications and about as many additions (or subtractions) for the LU-factorization plus about 3n N for solving (9) for each right-hand side (Ito et al. (99)). 4. Consistency, stability and convergence analysis In this section, we show that the approximate solution x N converges to the exact solution x of the system (4). In the following discussions, we put x N,i (s) = x N ((i )τ + s), s [0, τ ], i, (30) to denote the approximate solution x N on the interval [(i )τ, iτ ]. First, we prove the consistency and then the stability of the method will be proved. For the proof of the consistency let us to define the operator L as following where Lx i ẋ i A 0 x i, (3) In this case, x i (s) = x((i )τ + s), s [0, τ ], i. and D(L) = {x i C (0, τ )} Then, thanks to Lemma., we have D B (L) = {x i D(L) x i (0) = x i (τ )}. x i P N x i L w 0, as N, for all x i D B (L). (3)
13 AAM: Intern. J., Vol., Issue (June 07) 457 Therefore, the method is consistent. Second, we show that the method is stable, i.e., estimates of the form x N,i C x N,i, i =,, 3,, M, (33) can be obtained in some appropriately chosen norm and the constant C is independent of N. For doing this job, let us to define spaces and X N = {x i P N x i (0) = x i (τ )} Note that for all x i X N, we have Y N = P N. Therefore, P N Lx i = Lx i A 0 (P N x i x i ). (Lx i, P N Lx i ) w = Lx i L w (Lxi, A 0 (P N x i x i )) w According to Lemma., we have Lx i L w Lxi L w A 0 x i P N x i L w. (34) A 0 x i P N x i L w C 0 N x i,w. Moreover, it is possible to prove the following a priori estimate x i,w C Lx i L w, for a suitable constant C. Now, using (34), we get (Lx i, P N Lx i ) w ( C 0 C N ) Lx i L w (C ) x i,w, provided N is so large that
14 458 M. Mousa-Abadian & S. H. Momeni-Masuleh Since ( C 0 C N ). we conclude the following estimate P N Lx i L w C 3 x i,w, where α = C C 3. (Lx i, P N Lx i ) w P N Lx i L w α x i,w, (35) If we define and F = {x i H (0, τ ) x i (0) = x i (τ )} G = L w(0, τ ), then the stability of the method (i.e., inequality (33)) results from Lemma.3. Now, the convergence of the method can be established as a consequence of Lax-Richtmyer equivalence Theorem.. 5. Numerical results In this section, we provide some numerical examples and compare the proposed method with the exact and Runge-Kutta solutions. Runge-Kutta solutions are obtained by using dde3 function of Matlab 04a which is adapted for solving delay differential equations. The computational codes were conducted on an Intel(R) Core(TM) i7-6700k processor, equipped with 8 GB of RAM. Example. In this example (Banks and Kappel (979)), the equation for a damped oscillator with delay restorting force and constant external force is considered
15 AAM: Intern. J., Vol., Issue (June 07) 459 { ẍ(t) + ẋ(t) + x(t ) = 0, x(t) = cos t, ẋ(t) = sin t, t [, 0]. (36) We can write the above equation in terms of a linear first order system as follows: [ ] ẋ (t) = ẋ (t) [ ] [ 0 0 x (t) x (t) ] [ ] [ x (t ) x (t ) where x (t) = x(t) and x (t) = ẋ(t). The exact solution for 0 t is { x (t) = sin(t ) cos(t ) e t a + 0t 9 + sin() 9, x (t) = cos(t ) + sin(t ) + e t a + 0, where a = 0 cos() + sin(). The exact solution for t is where ] [ ] 0 +, (37) 0 x (t) = 5t 9t 0.5 sin(t ) e t+ (0 + 0t + 0.5) cos() 0.5 cos()te t+ + sin() ( 0.5e t te t+ + t ).36e t b, x (t) = 0t cos(t ) + 0te t cos()te t+ 0.5 sin()te t+ + sin() +.36e t b, b = 9 sin() cos() cosh() + cos() sinh() 0 cosh() +0 sinh() + sin() cosh() sin() sinh(). Numerical and analytical solutions are plotted in Figure and the maximum error between Rune-Kutta and the proposed method solutions is given in Table. As it was expected, the proposed method is more accurate than the Runge-Kutta method. In this example, the spectral method only uses 8 nodes to get the results, while Runge-Kutta method gives the results by taking 5 nodes. Note that if one try to obtain the accuracy of the spectral method, needs to refine the mesh and this leads to a heavy computation cost of Runge-Kutta method. Table : Maximum error between the exact solution, Runge-Kutta solution and the Chebyshev-tau spectral solution for 0 t of Example. Example. x i Runge-Kutta Current work x x Consider the scalar equation (Bellen and Zennaro (03))
16 460 M. Mousa-Abadian & S. H. Momeni-Masuleh 4 x (t) Runge Kutta solution Exact solution Spectral solution Runge Kutta solution Exact solution Spectral solution x (t) t Figure : Comparison of state vector (x (t) and x (t)) in the exact solution, Runge-Kutta solution and Chebyshev-tau spectral solution (N = 8) of Example. with the initial function ẋ(t) = x(t ), x(t) = 0.5t, t [, 0]. The exact solution for 0 t can be expressed as { x(t) = 4 t + t, 0 t, (t )3 4 t + t, t. In Figure numerical and exact solutions are demonstrated. For a longer time interval, the spectral and Runge-Kutta solutions are plotted in Figure 3 which confirms that the spectral method produces a smoother solution compared to Runge-Kutta method. In addition, the maximum error for Runge-Kutta method was with 37 nodes, while for the proposed method was which demonstrates the accuracy of the method as expected.
17 AAM: Intern. J., Vol., Issue (June 07) x(t) Spectral solution Exact solution Runge Kutta solution t Figure : Comparison between the exact solution, Runge-Kutta solution and Chebyshev-tau spectral solution (N = 8) of Example. Example 3. Consider the following system (Banks and Kappel (979)) ẋ(t) = 0 0 x(t) x(t ), which is equipped with the initial function x(t) =, t [, 0]. The analytical solution for 0 t is x (t) = 3 t3 + t +, x (t) = t +, x 3 (t) = t +,
18 46 M. Mousa-Abadian & S. H. Momeni-Masuleh Spectral solution Runge Kutta solution x(t) t Figure 3: Comparison between Runge-Kutta method and Chebyshev-tau spectral method (N = 8) of Example. and for t is x (t) = t + 4t + 3, x (t) = t +, x 3 (t) = 3 t3 + t Numerical simulations for N = 8 are demonstrated in Figure 4. In this case, the maximum error between exact and numerical simulations is given in Table. The number of nodes for Runge-Kutta method was 39. Table : Maximum error between the exact solution, Runge-Kutta solution and the Chebyshev-tau spectral solution for 0 t of Example 3. x i Runge-Kutta Current work x x x
19 AAM: Intern. J., Vol., Issue (June 07) 463 Example 4. In this example, we consider the following scalar equation with rational delay (Bellen and Zennaro (03)) with the following initial function ẋ(t) = x(t 0.5), The exact solution for 0 t is x(t) = 0.5t, t [ 0.5, 0]. { 0.5(t 0.5) + x(t) = 6, 0 t, (t )3 6 t , t. Approximate solution obtained by the Chebyshev-tau spectral method for N = 8 is illustrated in Figure 5. Moreover, the maximum error for Runge-Kutta method was , whereas for Chebyshev-tau spectral method was If the time interval is taken.5 x (t) Runge Kutta solution Exact solution Spectral solution x (t) 0 Runge Kutta solution Exact solution Spectral solution x 3 (t) 4 3 Runge Kutta solution Exact solution Spectral solution t Figure 4: Comparison of state vector (x (t), x (t) and x 3 (t)) in the exact solution, Runge- Kutta solution and Chebyshev-tau spectral solution (N = 8) of Example 3.
20 464 M. Mousa-Abadian & S. H. Momeni-Masuleh more longer, one can see that the Chebyshev-tau spectral method has a smooth transient response compared to Runge-Kutta method by using 37 nodes. This event can be observed in Figure x(t) Spectral solution Exact solution Runge Kutta solution t Figure 5: Comparison between the exact solution, Runge-Kutta solution and Chebyshev-tau spectral solution (N = 8) of Example 4. Example 5. In this example, we consider the following inhomogeneous system which constrained to the following history ẋ(t) = x(t ) + t, The exact solution for 0 t is x(t) = t, t [, 0]. { x(t) = 3 t3 + t t, 0 t, t4 + 6 t3 t t 3, t.
21 AAM: Intern. J., Vol., Issue (June 07) Spectral solution Runge Kutta solution x(t) t Figure 6: Comparison between Runge-Kutta method and Chebyshev-tau spectral method (N = 8) of Example 4. Numerical and analytical solutions are illustrated in Figure 7. Furthermore, the maximum error for Runge-Kutta method was by using nodes, whereas for Chebyshev-tau spectral method was Example 6. As a final example, consider the following system with two discrete point delays ẋ(t) = x(t 0.5) + x(t ), with the history function as x(t) = 0.5t, t [, 0]. The analytical solution for 0 t.5 is
22 466 M. Mousa-Abadian & S. H. Momeni-Masuleh.5 Runge-Kutta solution Exact solution Spectral solution x(t) t Figure 7: Comparison between the exact solution, Runge-Kutta method and Chebyshev-tau spectral method (N = 8) of Example t 0.75t, 0 t 0.5, x(t) = 6 (t 0.5)3 3 8 (t 0.5) + 0.5(t ) 5 6, 0.5 t, 4 t(t )3 + (t 3 )3 5 6 t 3 8 (t ) 6 96, t.5. Numerical and analytical solutions are depicted in Figure 8. The maximum error for Runge- Kutta method was by employing 37 nodes, while for the proposed method was which confirms the accuracy of the method as expected. The CPU time of the above examples are given in Table 3. As one observe, the CPU time required by the proposed method using 8 nodes is about one half of that used by Runge-Kutta method. Table 3: CPU time (s) for the proposed (N = 8) and Runge-Kutta method Example Current work Runge-Kutta No. nodes
23 AAM: Intern. J., Vol., Issue (June 07) Spectral solution Exact solution Runge Kutta solution 0.3 x(t) t Figure 8: Comparison between the exact solution, Runge-Kutta method and Chebyshev-tau spectral method (N = 8) of Example Conclusion In this article, the method of steps and Chebyshev-tau spectral method are combined to solve linear time delay systems of differential equations. Usually, Runge-Kutta method is used to solve ODE systems. Here, we used Chebyshev-tau spectral method to solve ODE systems which improve the accuracy of approximate solution dramatically and decrease the computational costs. In fact, we construct the approximate solution using small number of nodes. Implementation of the method leads to a linear system of algebraic equations which has, in general, a dense coefficient matrix. Here we convert the coefficient matrix to a block quasi upper triangular form using Chebyshev polynomials which can be solve more efficiently. Moreover, employing Chebyshev-tau spectral method enables us to apply FFT which reduces the CPU time of calculations. The consistency and stability analysis of the proposed method is provided, and consequently its convergence is established. Numerical results show a very good agreement with available literature. The proposed method can be extended to nonlinear problems and systems with more than two commensurate delays. This will be the subject of further studies. Acknowledgment: The authors would like to thank the anonymous referees for their valuable comments and suggestions which have helped to improve the manuscript. REFERENCES
24 468 M. Mousa-Abadian & S. H. Momeni-Masuleh Aziz, I. and Amin, R. (06). Numerical solution of a class of delay differential and delay partial differential equations via Haar wavelet, Appl. Math. Model., Vol. 40, pp Banks, H. T. and Kappel, F. (979). Spline approximations for functional differential equations, J. Differ. Equations, Vol. 3, pp Behroozifar, M. and Yousefi, S. A. (03). Numerical solution of delay differential equations via operational matrices of hybrid of block-pulse functions and Bernstein polynomials, C. M. D. E., Vol., pp Bellen, A. and Zennaro, M. (988). Stability properties of interpolants for Runge-Kutta methods, SIAM J. Numer. Anal., Vol. 5, pp Bellen, A. and Zennaro, M. (03). Numerical Methods for Delay Differential Equations, Oxford University Press. Bellman, R. and Cooke, K. L. (965). On the computational solution of a class of functional differential equations, J. Math. Anal. Appl., Vol., pp Bogacki, P. and Shampine, L. F. (990). Interpolating high-order Runge-Kutta formulas, Comput. Math. Appl., Vol. 0, pp Canuto, C., Hussaini, H. Y., Quarteroni, A. and Zang, T. A. (987). Spectral Methods in Fluid Dynamics, Springer. Delfour, M. C. and Mitter, S. K. (975). Hereditary differential systems with constant delays. II. A class of affine systems and the adjoint problem, J. Differ. Equations, Vol. 8, pp Elsgolts, L. E. and Norkin, S. B. (973). Introduction to the Theory and Application of Differential Equations with Deviating Arguments, Elsevier. Feldstein, A. and Neves, K. W. (984). High order methods for state-dependent delay differential equations with nonsmooth solutions, SIAM J. Numer. Anal., Vol., pp Fox, L. and Parker, I. B. (968). Chebyshev Polynomials in Numerical Analysis, Oxford University Press. Ghasemi, M. and Tavassoli Kajani, M. (0). Numerical solution of time-varying delay systems by Chebyshev wavelets, Appl. Math. Model., Vol. 35, pp Gorecki, H., Grabowski, P., Fuksa, S. and Korytowski, A. (989). Analysis and Synthesis of Time Delay Systems, John Wiley and Sons. Gottlieb, D. and Orszag, S. A. (977).Numerical Analysis of Spectral Methods: Theory and Applications, SIAM Press. Hale, J. K. and Lunel, S. M. V. (03). Introduction to Functional Differential Equations, Springer Science and Business Media. Isaacson, E. and Keller, H. B. (996). Analysis of Numerical Methods, John Wiley and Sons. Ito, K., Tran, H. T. and Manitius, A. (99). A fully-discrete spectral method for delaydifferential equations, SIAM J. Numer. Anal., Vol. 8, pp. 40. Kappel, F. and Kunisch, K. (98). Spline approximations for neutral functional differential equations, SIAM J. Numer. Anal., Vol. 8, pp Kemper, G. A. (975). Spline function approximation for solutions of functional differential equations, SIAM J. Numer. Anal., Vol., pp
25 AAM: Intern. J., Vol., Issue (June 07) 469 Kirschner, D. (996). Using mathematics to understand HIV immune dynamics, Notices Amer. Math. Soc., Vol. 43, pp Nelson, P. W. and Perelson, A. S. (00). Mathematical analysis of delay differential equation models of HIV- infection, Math. Biosci., Vol. 79, pp Oberle, H. J. and Pesch, H. J. (98). Numerical treatment of delay differential equations by Hermite interpolation, Numer. Math., Vol. 37, pp Sedaghat, S., Ordokhani, Y. and Dehghan, M. (0). Numerical solution of the delay differential equations of pantograph type via Chebyshev polynomials, Commun. Non. Sci. Numer. Simulat., Vol. 7, pp Wang, Z. Q. and Wang, L. L (00). A Legendre-Gauss collocation method for nonlinear delay differential equations, Discret Contin. Dyn. S., Vol. 3, pp Zennaro, M. (985a). On the P-stability of one-step collocation for delay differential equations, International Series of Numerical Mathematics, Vol., 74 pp Zennaro, M. (985b). One-step collocation: Uniform superconvergence, predictor-corrector method, local error estimate, SIAM J. Numer. Anal., Vol., pp
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 informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 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
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 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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
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 informationA Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind
Journal of Mathematical Extension Vol. 8, No. 1, (214), 69-86 A Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind Sh. Javadi
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
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 informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
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 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 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 information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
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 informationCollocation approximation of the monodromy operator of periodic, linear DDEs
p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,
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 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 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 informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationNumerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method
Punjab University Journal of Mathematics (ISSN 116-2526) Vol.48(1)(216) pp. 47-53 Numerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method Kobra Karimi Department
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
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 informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationRATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages 72 83 c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR
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 informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
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 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 informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
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 informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationA new class of highly accurate differentiation schemes based on the prolate spheroidal wave functions
We introduce a new class of numerical differentiation schemes constructed for the efficient solution of time-dependent PDEs that arise in wave phenomena. The schemes are constructed via the prolate spheroidal
More informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
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 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 information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationA Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals
arxiv:1503.00037v2 [math.na] 19 Mar 2015 A Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals Riccardo Fazio and Alessandra Jannelli Department
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
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 informationBanach Journal of Mathematical Analysis ISSN: (electronic)
Banach J. Math. Anal. 6 (2012), no. 1, 139 146 Banach Journal of Mathematical Analysis ISSN: 1735-8787 (electronic) www.emis.de/journals/bjma/ AN EXTENSION OF KY FAN S DOMINANCE THEOREM RAHIM ALIZADEH
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationUsing Hybrid Functions to solve a Coupled System of Fredholm Integro-Differential Equations of the Second Kind
Advances in Dynamical Systems and Applications ISSN 973-532, Volume 9, Number, pp 5 (24) http://campusmstedu/adsa Using Hybrid Functions to solve a Coupled System of Fredholm Integro-Differential Euations
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 informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationA Numerical Solution of Volterra s Population Growth Model Based on Hybrid Function
IT. J. BIOAUTOMATIO, 27, 2(), 9-2 A umerical Solution of Volterra s Population Growth Model Based on Hybrid Function Saeid Jahangiri, Khosrow Maleknejad *, Majid Tavassoli Kajani 2 Department of Mathematics
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
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 Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationIntroduction to Signal Spaces
Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationIterated Defect Correction with B-Splines for a Class of Strongly Nonlinear Two-Point Boundary Value Problems
American Review of Mathematics and Statistics June 2016, Vol. 4, No. 1, pp. 31-44 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). All Rights Reserved. Published by American Research
More informationConstrained controllability of semilinear systems with delayed controls
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Constrained controllability of semilinear systems with delayed controls J. KLAMKA Institute of Control Engineering, Silesian
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationApproximate solution of linear integro-differential equations by using modified Taylor expansion method
ISSN 1 746-7233, Engl, UK World Journal of Modelling Simulation Vol 9 (213) No 4, pp 289-31 Approximate solution of linear integro-differential equations by using modified Taylor expansion method Jalil
More informationApplication of the Bernstein Polynomials. for Solving the Nonlinear Fredholm. Integro-Differential Equations
Journal of Applied Mathematics & Bioinformatics, vol., no.,, 3-3 ISSN: 79-66 (print), 79-6939 (online) International Scientific Press, Application of the Bernstein Polynomials for Solving the Nonlinear
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 234 (21) 1364 1373 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationS. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA
RECURSIVE MULTIPLE WIENER INTEGRAL EXPANSION FOR NONLINEAR FILTERING OF DIFFUSION PROCESSES Published in: J. A. Goldstein, N. E. Gretsky, and J. J. Uhl (editors), Stochastic Processes and Functional Analysis,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationHausdorff operators in H p spaces, 0 < p < 1
Hausdorff operators in H p spaces, 0 < p < 1 Elijah Liflyand joint work with Akihiko Miyachi Bar-Ilan University June, 2018 Elijah Liflyand joint work with Akihiko Miyachi (Bar-Ilan University) Hausdorff
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationThe variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
Cent. Eur. J. Eng. 4 24 64-7 DOI:.2478/s353-3-4-6 Central European Journal of Engineering The variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
More informationApplication of Taylor-Padé technique for obtaining approximate solution for system of linear Fredholm integro-differential equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 1 (June 2017), pp. 392 404 Applications and Applied Mathematics: An International Journal (AAM) Application of Taylor-Padé
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
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 informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction 1 Introduction In this module, we develop solution techniques for numerically solving ordinary
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
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 informationNumerical Methods for Two Point Boundary Value Problems
Numerical Methods for Two Point Boundary Value Problems Graeme Fairweather and Ian Gladwell 1 Finite Difference Methods 1.1 Introduction Consider the second order linear two point boundary value problem
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationBASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA
1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.
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 informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More information