- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2013 (2013), Article ID 198926, 11 pages

http://dx.doi.org/10.1155/2013/198926

## Fourier Operational Matrices of Differentiation and Transmission: Introduction and Applications

^{1}Department of Applied Mathematics, School of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad, Iran^{2}Center of Excellence on Modelling and Control Systems, Ferdowsi University of Mashhad, Mashhad, Iran^{3}Department of Mathematics, Universiti Putra Malaysia (UPM), 43400 Serdang, Selangor, Malaysia

Received 1 January 2013; Revised 4 March 2013; Accepted 4 March 2013

Academic Editor: Carlos Vazquez

Copyright © 2013 F. Toutounian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper introduces Fourier operational matrices of differentiation and transmission for solving high-order linear differential and difference equations with constant coefficients. Moreover, we extend our methods for generalized Pantograph equations with variable coefficients by using Legendre Gauss collocation nodes. In the case of numerical solution of Pantograph equation, an error problem is constructed by means of the residual function and this error problem is solved by using the mentioned collocation scheme. When the exact solution of the problem is not known, the absolute errors can be computed approximately by the numerical solution of the error problem. The reliability and efficiency of the presented approaches are demonstrated by several numerical examples, and also the results are compared with different methods.

#### 1. Introduction

In the last four decades, numerical methods which are based on the operational matrices of integration (especially for orthogonal polynomials and functions) have received considerable attention for dealing with a huge size of applied mathematics problems such as *system identification*, *state space analysis*, *optimal control*, and *senstivity analysis*. The key idea of these methods is based on the integral expression
where is an arbitrary basis vector and is a constant matrix, called the operational matrix of integration. The matrix has already been determined for many types of orthogonal (or nonorthogonal) bases such as Walsh functions [5–7], block-pulse functions [8], Laguerre polynomials [9], Chebyshev polynomials [10], Legendre polynomials [11], Hermite polynomials [12], Fourier series [13], Bernstein polynomials [14], and Bessel series [15]. As a primary research work which was based on the operational matrices of integration, one can refer to the work of Corrington [5]. In [5], the author proposed a method of solving nonlinear differential and integral equations using a set of Walsh functions as the basis. His method is aimed at obtaining piecewise constant solutions of dynamic equations and requires previously prepared tables of coefficients for integrating Walsh functions. To alleviate the need for such tables, Chen and Hsiao [6, 7] introduced an operational matrix to perform integration of Walsh functions. This operational matrix approach has been applied to various problems such as time-domain analysis and synthesis of linear systems, piecewise constant-feedback-gain determination for optimal control of linear systems and for inverting irrational Laplace transforms.

On the other hand, since the beginning of 1994, the Bernoulli, Chebyshev, Laguerre, Bernstein, Legendre, Taylor, Hermite, and Bessel matrix methods have been used in the works [4, 16–28] to solve high-order linear and nonlinear Fredholm Volterra integro-differential difference equations and their systems. The main characteristic of these approaches is based on the operational matrices of differentiation instead of integration. The best advantage of these techniques with respect to the integration methods is that, in the fundamental matrix relations, there is not any approximation symbol, meanwhile in the integration forms such as (1), the approximation symbol could be seen obviously. In other words, where is the operational matrix of differentiation for any selected basis such as the above-mentioned polynomials, functions, and truncated series. The readers can see that there is no approximation symbol in (2), meanwhile this can be seen in (1) by using operational matrices of integration. For justifying this expression, one can refer to this subject that after the differentiating of an th degree polynomial we usually reach to a polynomial which has less than th degree. However, in the integration processes the degree of polynomials would be increased.

To the best of our knowledge, this is the first work concerning the Fourier matrix method for solving high-order linear ordinary differential equations (ODEs) in a differentiation form of view. This partially motivated our interest in such methods. In this paper, in the light of the methods [4, 16–20, 24, 25] and by means of the matrix relations between the truncated Fourier series and their derivatives, we propose a new method, namely, the Fourier matrix method for solving th order linear ODEs with constant coefficients in the form where is a known function and the initial conditions , for are given.

Also, we present a new method by using matrix relations between the truncated Fourier series and their transmissions for solving the following th order difference equation: Moreover, by means of the matrix relations between the truncated Fourier series and their derivatives together with their transmissions and by using suitable Gaussian collocation nodes, we obtain the numerical solution of generalized Pantograph equation with variable coefficients in the form subject to the initial conditions , for .

It should be noted that delay differential equations have a wide range of application in science and engineering. Functional-differential equations with proportional delays are usually refereed as Pantograph equations or generalized Pantograph equations. The applications of Pantograph equation are in different fields such as number theory, economy, biology, control, electrodynamics, nonlinear dynamical systems, quantum mechanics, probability theory, astrophysics, cell growth, and other industrial applications. For some applications of this equation, we refer the interested reader to [21]. Pantograph equation was studied by many authors numerically and analytically. A complete survey about numerical and analytical methods about generalized Pantograph equations has been provided in [21].

The rest of the paper is organized as follows. In Section 2, we introduce some mathematical preliminaries of Fourier series together with their operational matrices of differentiation and transmission. Section 3 is devoted to apply the Fourier matrix methods for solving (3), (4), and (5) using the Fourier operational matrices of differentiation and transmission. An error estimation of the collocation scheme for solving the generalized Pantograph equation (5) is provided in Section 4. In Section 5, the proposed method is applied to several numerical examples. Also conclusions and future works are given in Section 6.

#### 2. Some Properties of Fourier Series

Fourier series decomposes periodic functions or periodic signals into the sum of a (possibly infinite) set of simple oscillating functions, namely, sines and cosines. The Fourier series has many applications in electrical engineering, vibration analysis, acoustics, optics, signal processing, image processing, quantum mechanics, econometrics, thin-walled shell theory, and so forth. A square integrable function can be expanded in terms of Fourier series as follows: where the coefficients and are given by In practice, we use only the first terms of the Fourier series. Therefore, our aim is to approximate the solution of (3), (4), and (5) as follows: The vector representation of the above-mentioned truncated series can be written in the following form: According to the property , we have , where denotes the first derivative of . By repeating this procedure, we have , where denotes the th derivative of . On the other hand, the matrix relation between and can be constructed by using the elementary calculus as follows: where is the Fourier operational matrix of differentiation. The above relation can be written in the form . Accordingly, the th derivative of can be given by Therefore, has the following form: We now construct the operational matrices of transmission in both cases of forward and backward by using the following simple relations from the elementary calculus: Thus, the relation between and is as follows: Moreover, the relation between and can be written in the following form: In other words, and . Therefore, It should be noted that (12) and (16) are the fundamental matrix relations of this section. In other words, the operational matrices of differentiation and transmission would be applied in our methods in the next section abundantly.

#### 3. Matrix and Collocation Methods Based on Fourier Operational Matrices

In this section, we use the matrix relations (12) and (16) for solving (3), (4), and (5). In the solution approximation procedure of (3), we use only the fundamental relation (12), and for approximating the solution of (4), we use only the fundamental relation (16), meanwhile both of relations (12) and (16) together with the collocation method based on Legendre Gauss points are applied to solve numerically the generalized Pantograph equation (5).

##### 3.1. Fourier Matrix Method for Solving High-Order Linear ODEs

We now derive an algorithm for solving (3). To do this, let the solution of (3) be approximated by the first -terms of the truncated Fourier series. Hence, if we write where the unknown Fourier coefficient vector and the Fourier vector are given in (9), then the th derivative of can be expressed in the matrix form by for . Moreover, we need to approximate the given function which exists in (3) by using (7). In other words, we assume that where is a known vector with the aid of (7). By substituting the approximations (12) and (18) into (3), we get According to the completeness property of Fourier series, the matrix equation corresponding to the basic equation (3) can be written as or in the form of augmented matrix , where Therefore, the basic equation (3) is changed into a system of linear algebraic equations with unknown coefficients which can be written in the form or in augmented matrix form Now we should impose the initial conditions for to the associated system of algebraic equations in a proper way. For this reason and for clarity of presentation, we can remove the last rows of and the last entries of without loss of generality. However, we can do this work by the first rows of the and the first entries of . Thus, we have According to (12), we can write and hence the matrix-vector representation of the initial conditions is as follows: Therefore, the matrix form of the initial conditions can be written as where Finally, from (23) and (25), (3) subject to the considered initial conditions reduces to the following system of algebraic equations: where The matrix coefficients of the above system is sparse, and for solving system of linear equations , we can use some iterative Krylov subspace methods and determine the vector and hence the approximate solution is obtained.

##### 3.2. Fourier Matrix Method for Solving High-Order Linear Difference Equations

In this subsection, we propose an idea for approximating the solution of (4) by using Fourier operational matrix of transmission. For this purpose, we assume that the solution of (4) can be approximated by truncated Fourier series in the form . According to (16), we have for . Similar to the previous method, the known function can be approximated by truncated Fourier series in the form where the entries of would be calculated from (7). By substituting these matrix forms in (4), we get Hence, the matrix equation corresponding to the basic equation (4) can be written as or in the form of augmented matrix , where Therefore, the basic equation (4) is changed into a system of linear algebraic equations with unknown coefficients which can be written in the form or in augmented matrix form According to the structure of Fourier operational matrices of transmissions, the matrix coefficients of the above system is tridiagonal (see for instance (14)), and for solving system of linear equations , we can use some iterative krylov subspace methods and determine the vector and hence the approximated solution is obtained.

##### 3.3. Fourier Collocation Method for Solving Generalized Pantograph Equations

In this subsection, we use the collocation method (via Legendre Gauss points) based on Fourier operational matrices of differentiation and transmission to solve numerically the generalized Pantograph equation (5). To do this, let the solution of (5) be approximated by the first -terms truncated Fourier series in the form , where the Fourier coefficient vector is unknown (should be determined) and the Fourier vector is given in (9). According to (12), the th derivative of can be expressed in the matrix form .

By substituting the above approximations into (5), we get To find the solution , we first collocate (33) at the nodes as the roots of (where denotes the shifted Legendre polynomial of degree in the interval ) yields where The matrix system (34) can be rewritten as follows: where Hence, the matrix equation (36) corresponding to the basic equation (5) can be written as or in the form of augmented matrix , where Therefore, the basic equation (5) is changed into a system of linear algebraic equations with unknown coefficients , which can be written in the form or in augmented matrix form

Now we should impose the initial conditions to the associated algebraic system of equations in a proper manner. These initial conditions can be written as Therefore, the matrix form of the initial conditions can be written as where Finally, from (40) and (42), the generalized Pantograph equation (5) subject to the initial conditions reduces to the following system of algebraic equations: where The matrix coefficients of the above system may be sparse, and for solving system of linear equations , we can use some iterative krylov subspace methods and determine the vector and hence the approximated solution is obtained.

*Remark 1. *For more information about iterative krylov subspace methods, one can point out to [29]. In this book, several iterative methods have been introduced for solving large sparse linear systems.

#### 4. Error Estimation

In this section, we will give an error estimation (which was previously proposed in the works [27, 28]) for the collocation scheme that is proposed for solving Pantograph equation (5) with the aid of residual function. This idea may help us to estimate the error especially if the exact solution of (5) does not exist. For this purpose, we can rewrite (5) in the form , where with the initial conditions . On the other hand, our collocation scheme finds the approximate solution of the problem , where the residual function vanishes at the collocation nodes and with the initial conditions . Now, we define the error function as . Trivially one can write with the initial conditions . In other words, we construct the error problem as follows: with the initial conditions . Similar to the previous collocation method, we now solve the above error problem and approximate by the aid of truncated Fourier series (with the terms) in the following form: If the exact solution of (5) is unknown, then the absolute errors are not found. However, the absolute errors can be computed approximately with the aid of the estimated absolute error function .

#### 5. Numerical Examples

In this section, several numerical examples are given to illustrate the accuracy and effectiveness of the proposed methods, and all of them are performed on a computer using programs written in MAPLE 13. In this regard, we have reported in the tables and figures the values of the exact solution , the approximate solution , the absolute error function , and estimate absolute error function (or the absolute residual functions associated to (3), (4), and (5)) at any selected points of the given computational interval. It should be noted that in the first two examples we provide equations in which our results are exact, meanwhile other approaches which were based on the operational matrices of integration and differentiation such as [5, 11, 14, 16–18] could not obtain the exact solutions.

*Example 2. *As the first example we consider the following ODE with the exact solution
subject to the initial condition .

According to (8), one can approximate the solution of the above ODE with assumption as follows:
Our aim is to determine the vector components with the aid of Fourier operational matrix of differentiation. By using (10) (i.e., ), we have
Moreover, we can approximate the known function by the aid of the truncated Fourier series in the form
By substituting the matrix forms of (51)–(53) into (50) we have
Since , then . Thus, the above matrix equation transforms into the following matrix equation:
The solution of the above system is and hence
which is the exact solution.

*Example 3. *As the second example we consider the following linear difference equation with the exact solution

According to (8), one can approximate the solution of the above difference equation with assumption as follows:
Our aim is to determine the vector components with the aid of Fourier operational matrix of transmission. By using (16) (i.e., , ), we have
Moreover, we can approximate the known function by the aid of the truncated Fourier series in the form
By substituting the matrix forms of (59)–(61) into (57), we have
The solution of the above system is and hence
which is the exact solution.

*Example 4. *As the third example, we consider the following Pantograph equation with the exact solution
Since the computational interval of this problem is , we need the roots of the shifted Legendre polynomial for any arbitrary values of . Now, we solve this Pantograph equation by means of the proposed method in Section 3 for different values of such as and 3. These approximate solutions are as follows:

In Table 1, we provide the absolute values of errors for and 3. Moreover, we estimate the mentioned errors by solving the associated error problem which is constructed in Section 4 by using , and 5. This table shows the agreement of the estimate errors with the actual errors.

*Example 5. * As the final example we consider the linear Pantograph differential equation of first order
Since the problem (66) has a complicated exact solution, we compare the results of collocation method based on presented method (PM) with the other four methods in Table 2. These methods that are used to obtain the approximate solutions of this Example are Walsh series approach [1], Laguerre series technique [2], Taylor series method [3], and Hermit series collocation approach [4]. Moreover, it seems that the Walsh series approach [1] and Laguerre series technique [2] are not efficient as other Taylor, Hermite, and presented Fourier methods. Therefore, for illustrating the efficiency of the proposed method, we depict the history of the residual of the obtained solution of our scheme together with the associated residuals of Taylor and Hermite methods in Figure 1. From this figure, one can see that the presented method is more applicable and efficient for solving such Pantograph equations.

#### 6. Conclusions and Future Works

The aim of this paper is to introduce two new operational matrices of Fourier series for solving a huge class of high-order linear differential equations. It should be noted that operational matrices of differentiation have more accuracy with respect to the integration ones usually. Also Fourier series approximation is very effective and useful for ODEs with periodic solutions. These two advantages together with the superior results of the proposed scheme (in the numerical illustrations) encourage the authors to extend this method for solving high-order linear and nonlinear Fredholm integro-differential equations [16], systems of Volterra integral equations [30], high-order linear partial differential equations [20, 21], and optimal control problems [31].

#### Conflict of Interests

The authors declare that they do not have any conflict of interests in their submitted paper.

#### Acknowledgment

The authors thank the referee for his or her valuable comments and helpful suggestions which led to the improved version of the paper.

#### References

- G. P. Rao and K. R. Palaisamy, “Walsh stretch matrices and functional differential equations,”
*IEEE Transactions on Automatic Control*, vol. 27, no. 1, pp. 272–276, 1982. View at Google Scholar · View at Scopus - C. Hwang and Y. P. Shih, “Laguerre series solution of a functional differential equation,”
*International Journal of Systems Science*, vol. 13, no. 7, pp. 783–788, 1982. View at Google Scholar - M. Sezer, S. Yalçinbaş, and M. Gülsu, “A Taylor polynomial approach for solving generalized pantograph equations with nonhomogenous term,”
*International Journal of Computer Mathematics*, vol. 85, no. 7, pp. 1055–1063, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Yalçinbaç, M. Aynigül, and M. Sezer, “A collocation method using Hermite polynomials for approximate solution of pantograph equations,”
*Journal of the Franklin Institute*, vol. 348, no. 6, pp. 1128–1139, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. S. Corrington, “Solution of differential and integral equations with walsh functions,”
*IEEE Transactions on Circuit Theory*, vol. 20, no. 5, pp. 470–476, 1973. View at Google Scholar · View at Scopus - C. F. Chen and C. H. Hsiao, “Time-domain synthesis via walsh functions,”
*Proceedings of the Institution of Electrical Engineers*, vol. 122, no. 5, pp. 565–570, 1975. View at Google Scholar · View at Scopus - C. F. Chen and C. H. Hsiao, “Walsh series analysis in optimal control,”
*International Journal of Control*, vol. 21, no. 6, pp. 881–897, 1975. View at Google Scholar · View at Scopus - N. S. Hsu and B. Cheng, “Analysis and optimal control of time-varying linear systems via block-pulse functions,”
*International Journal of Control*, vol. 33, no. 6, pp. 1107–1122, 1981. View at Google Scholar · View at Scopus - C. Hwang and Y. P. Shih, “Parameter identification via laguerre polynomials,”
*International Journal of Systems Science*, vol. 13, no. 2, pp. 209–217, 1982. View at Google Scholar · View at Scopus - I. R. Horng and J. H. Chou, “Shifted chebyshev direct method for solving variational problems,”
*International Journal of Systems Science*, vol. 16, no. 7, pp. 855–861, 1985. View at Google Scholar · View at Scopus - R. Y. Chang and M. L. Wang, “Shifted Legendre direct method for variational problems,”
*Journal of Optimization Theory and Applications*, vol. 39, no. 2, pp. 299–307, 1983. View at Publisher · View at Google Scholar · View at Scopus - G. T. Kekkeris and P. N. Paraskevopoulos, “Hermite series approach to optimal control,”
*International Journal of Control*, vol. 47, no. 2, pp. 557–567, 1988. View at Google Scholar · View at Scopus - M. Razzaghi and M. Razzaghi, “Fourier series direct method for variational problems,”
*International Journal of Control*, vol. 48, no. 3, pp. 887–895, 1988. View at Google Scholar · View at Scopus - E. H. Doha, A. H. Bhrawy, and M. A. Saker, “Integrals of Bernstein polynomials: an application for the solution of high even-order differential equations,”
*Applied Mathematics Letters*, vol. 24, no. 4, pp. 559–565, 2011. View at Publisher · View at Google Scholar · View at Scopus - P. N. Paraskevopoulos, P. G. Sklavounos, and G. C. Georgiou, “The operational matrix of integration for Bessel functions,”
*Journal of the Franklin Institute*, vol. 327, no. 2, pp. 329–341, 1990. View at Google Scholar · View at Scopus - A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals,”
*Applied Mathematics and Computation*, vol. 219, no. 2, pp. 482–497, 2012. View at Publisher · View at Google Scholar - A. Akyüz-Dascioglu, “Chebyshev polynomial approximation for high-order partial differential equations with complicated conditions,”
*Numerical Methods for Partial Differential Equations*, vol. 25, no. 3, pp. 610–621, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Gülsu, B. Gürbüz, Y. Öztürk, and M. Sezer, “Laguerre polynomial approach for solving linear delay difference equations,”
*Applied Mathematics and Computation*, vol. 217, no. 15, pp. 6765–6776, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. Sezer, S. Yalçinbaş, and N. Şahin, “Approximate solution of multi-pantograph equation with variable coefficients,”
*Journal of Computational and Applied Mathematics*, vol. 214, no. 2, pp. 406–416, 2008. View at Publisher · View at Google Scholar · View at Scopus - E. Tohidi, “Legendre approximation for solving linear HPDEs and comparison with taylor and bernoulli matrix methods,”
*Applied Mathematics*, vol. 3, pp. 410–416, 2012. View at Google Scholar - E. Tohidi, A. H. Bhrawy, and Kh. Erfani, “A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation,”
*Applied Mathematical Modelling*, vol. 37, no. 6, pp. 4283–4294, 2013. View at Publisher · View at Google Scholar - E. Tohidi, “Bernoulli matrix approach for solving two dimensional linear hyperbolic partial differential equations with constant coefficients,”
*American Journal of Computational and Applied Mathematics*, vol. 2, no. 4, pp. 136–139, 2012. View at Publisher · View at Google Scholar - F. Toutounian, E. Tohidi, and S. Shateyi, “A collocation method based on Bernoulli operational matrix for solving high order linear complex differential equations in a rectangular domain,”
*Abstract, Applied, Analysis*, Article ID 823098, 2013. View at Publisher · View at Google Scholar - S. A. Yousefi and M. Behroozifar, “Operational matrices of Bernstein polynomials and their applications,”
*International Journal of Systems Science*, vol. 41, no. 6, pp. 709–716, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. Yuzbasi,
*Bessel polynomial solutions of linear differential, integral and integro-differential equations [M.S. thesis]*, Graduate School of Natural and Applied Sciences, Mugla University, 2009. - S. Yuzbasi, “A numerical approximation based on the Bessel functions of first kind for solutions of Riccati type differen- tialdifference equations,”
*Computers & Mathematics with Applications*, vol. 64, no. 6, pp. 1691–1705, 2012. View at Publisher · View at Google Scholar - S. Yuzbasi and M. Sezer, “An improved Bessel collocation method with a residual error function to solve a class of LaneEmden differential equations,”
*Computer Modeling*, vol. 57, pp. 1298–1311, 2013. View at Google Scholar - S. Yuzbasi, M. Sezer, and B. Kemanci, “Numerical solutions of integro-differential equations and application of a population model with an improved Legendre method,”
*Applied Mathematical Modelling*, vol. 37, pp. 2086–2101, 2013. View at Google Scholar - Y. Saad,
*Iterative Methods for Sparse Linear Systems*, SIAM, Philadelphia, Pa, USA, 2nd edition, 2003. - O. R. N. Samadi and E. Tohidi, “The spectral method for solving systems of Volterra integral equations,”
*Journal of Applied Mathematics and Computing*, vol. 40, no. 1-2, pp. 477–497, 2012. View at Publisher · View at Google Scholar - E. Tohidi and O. R. N. Samadi, “Optimal control of nonlinear Volterra integral equations via Legendre polynomials,”
*IMA Journal of Mathematical Control and Information*, vol. 30, no. 1, pp. 67–83, 2013. View at Publisher · View at Google Scholar