Research Article  Open Access
H. R. Marzban, S. M. Hoseini, "Analysis of Linear Piecewise Constant Delay Systems Using a Hybrid Numerical Scheme", Advances in Numerical Analysis, vol. 2016, Article ID 8903184, 9 pages, 2016. https://doi.org/10.1155/2016/8903184
Analysis of Linear Piecewise Constant Delay Systems Using a Hybrid Numerical Scheme
Abstract
An efficient computational technique for solving linear delay differential equations with a piecewise constant delay function is presented. The new approach is based on a hybrid of blockpulse functions and Legendre polynomials. A key feature of the proposed framework is the excellent representation of smooth and especially piecewise smooth functions. The operational matrices of delay, derivative, and product corresponding to the mentioned hybrid functions are implemented to transform the original problem into a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the proposed numerical scheme.
1. Introduction
Delay differential equations (DDEs) naturally arise in diverse areas of science and engineering such as transmission lines, communication networks, biological models, population dynamics, and transportation systems [1, 2]. So far, a large body of literature has been devoted to the theoretical aspects and numerical treatments of DDEs with constant delays [3–18]. It is known that, except for some simple cases, it is either extremely difficult or impossible to analytically solve DDEs. Accordingly, a numerical algorithm has to be adopted in most cases. The situation becomes more complicated when the timedelay is a piecewise constant function. Owing to the nature of DDEs, none of the smooth basis functions is able to properly model the inherent behavior of this class of systems. This is due to the lack of smoothness of analytical solution associated with DDEs. It should be pointed out that the approximation of a piecewise smooth function by a finite number of smooth functions often fails to converge because of the existence of the wellknown Gibbs phenomenon. Consequently, a suitable basis is required to accurately model the true locations of the switching points that occur in the exact solution of a delay differential equation. It is generally assumed that the delay function is either constant or continuous. However, in some situations, timedelay is a piecewise constant function. To the best of our knowledge, few research works have been dedicated to the development of computational algorithms for solving DDEs involving piecewise constant delay [4, 19]. Recently, Marzban and Shahsiah proposed an efficient numerical scheme for solving DDEs containing piecewise constant delay. Their method is based on a hybrid of blockpulse functions and Chebyshev polynomials. It has been demonstrated that the method implemented in [19] is effective and produces very accurate results.
In what follows, we describe some similarities and differences between our method and the procedures developed in [15, 16]. First, the foundations of the proposed framework and those used in the previous works are based on a hybrid of blockpulse functions and Legendre polynomials. Second, the current paper is an extension of our previous ones. More specifically, the timedelay systems considered in [15, 16] involve constant delay, while here we investigate linear piecewise constant delay systems. Obviously, the latter systems are a general class of constant DDEs. Third, the approach employed here is based on the derivative matrix corresponding to the mentioned hybrid functions while the method implemented in our earlier works is based on the operational matrix of integration. Fourth, the operational matrix of delay associated with the piecewise constant delay systems is constructed. To derive this matrix, we use the same approach as that of [19]. The purpose of this paper is to introduce an efficient numerical technique for solving DDEs with a piecewise constant delay. It should be emphasized that the analytical solutions of DDEs cannot be obtained solely either by blockpulse functions or by Legendre polynomials. Combining blockpulse functions and Legendre polynomials allows one to simultaneously make use of the best properties of the two mentioned bases. It is worth noting that the value of , the order of blockpulse functions, plays an essential role in modelling of the problem under consideration. Indeed, by selecting the suitable value of , we are able to correctly determine the exact locations of the switching points that occur in the solution associated with a piecewise constant delay system. The excellent properties of hybrid functions together with the operational matrices of delay, derivative, and product are then utilized to transform the delay differential equation under investigation into a system of algebraic equations whose solution is much easier than the original one.
The rest of the paper is organized as follows. In Section 2, the basic properties of hybrid of blockpulse functions and Legendre polynomials are presented. In Section 3, the operational matrices of derivative, product, and delay are presented. Section 4 is devoted to the problem statement and its approximation. In Section 5, three examples are tested to show the efficiency and accuracy of the proposed numerical scheme.
2. Hybrid Functions
Hybrid functions , , , are defined on the interval as [15]where and are the orders of blockpulse functions and Legendre polynomials, respectively. Here, are the wellknown Legendre polynomials of order which are orthogonal on the interval and satisfy the following recursive formula [20]:Since consists of blockpulse functions and Legendre polynomials, which are both complete and orthogonal, the set of the hybrid of blockpulse functions and Legendre polynomials is a complete orthogonal set in the Hilbert space .
3. Operational Matrices
In this section, we first present the operational matrix of derivative based on the weak representation of the derivative operator. We then state the operational matrices of delay and product corresponding to the proposed hybrid functions.
A function can be approximated by the hybrid functions as follows:whereis the vector of coefficients andThe coefficients , , , are obtained by the following formula:
3.1. The Operational Matrix of Derivative
We approximate the derivative of by the derivative of (3); that is,The right hand side of the preceding equation can be represented in terms of hybrid functions aswhere
The relationship between the two vectors and is expressed bywhere matrix is the operational matrix of derivative.
The structure of the operational matrix of derivative is given by [21].
In the case ,in the case ,and, in the case ,where, for , we set
3.2. The Operational Matrix of Product
Let be an arbitrary vector of order . Then, the expression can be expanded in terms of hybrid functions as follows:in which is a matrix of order . This matrix is called the operational matrix of product and has the following structure:where , , are matrices given in [15].
3.3. The Operational Matrix of Delay
The goal of this subsection is to determine the operational matrix of delay associated with the developed hybrid functions. For this purpose, letin which is defined bywhere and , , are known constants. Furthermore, assume that is a nonempty set, , , are nonnegative rational numbers, and for .
To obtain the operational matrix of delay corresponding to the proposed hybrid functions, we apply an approach analogous to the one devised in [19]. To do this, we first divide the time interval into subintervals of the same length where the value of is obtained in the following manner.
Define as the smallest positive integer number in such a way thatSuppose that . Next, we choose as the greatest common divisor of the integers and , , ; that is, Letwhere denotes the greatest integer value less than or equal to . It should be noted that is chosen in such a way that the number of subintervals is minimized. As a consequence, we obtain the following subintervals:where . Now, from (20) and (22) we can find the integer numbers and , such thatNote that , for , . Therefore, we can rewrite as follows:For convenience, we define as
Therefore, the problem is reduced to find the operational matrix of delay for the following delay function:where
In order to find the matrix , we first derive the matrix for , so that the following relation holds:
With the use of (1), it is worth noting that, in the case of , the only functions with nonzero values are , for . Since by expanding in terms of , we get the identity matrix as the corresponding coefficients. Therefore, we have where is the dimensional identity matrix and is an matrix in which the only nonzero entry is equal to one and located at the th row and th column.
Remark 1. If , then is a zero matrix of order .
As a consequence, if we expand in terms of hybrid functions , we deduce that
To illustrate the derivation process of the operational matrix of delay associated with the mentioned hybrid functions, we present an example. DefineIt is easily verified that . As a result, we get For and , we conclude that . Therefore, and are zero matrices of order . If , it follows that . Consequently, is a matrix in which the only nonzero entry is equal to one and located at the 1st row and 2nd column. If , then we deduce that . Hence, is a matrix in which the only nonzero entry is equal to one and located at the 1st row and 4th column. A similar argument can be used for and to specify matrices , and . Using the above comments, we obtain
Consequently
4. Problem Statement and Its Approximation
4.1. Piecewise Constant Delay Systems
Consider the linear timevarying piecewise constant delay system described bywhere , , and , , and are matrices with appropriate dimensions. The problem is to find , , satisfying (37)–(39).
4.2. Approximation Using Hybrid Functions
We approximate the system dynamics by hybrid functions as follows. LetUsing (3), each and each , , , can be written in terms of hybrid functions as follows:From (40), it follows thatwhere and are the and identity matrices, respectively, and denotes the Kronecker product [22]. It should be noted that and are vectors of orders and , respectively, given by Alsowhere is a vector of order defined by Similarly where , and are matrices of orders , , and , respectively. The delay vector can also be expanded in terms of hybrid functions aswhere is the operational matrix of delay described by (18).
Now, using (16), we havewhere and can be calculated in a way similar to the construction method of matrix given by (17). Moreover, ConsequentlyNow, using the operational matrix of derivative and substituting (44)–(50) in (37), we get Because the elements of are linearly independent functions over the interval , it follows that
Therefore, the original problem is transformed into a system of linear algebraic equations. The aforementioned system with the associated initial condition can be easily solved by the wellknown Tau method [20, 23] for the unknown vector .
5. Illustrative Examples
In this section, three examples are investigated to evaluate the performance of the method.
Example 1. Consider the timevarying piecewise constant delay systemThe exact solution to this problem is given by [19]To solve this problem by the method developed in the present paper, we first determine the value of , the required number of subintervals, with the use of (22). For this problem, we choose . Also, we select . Let By expanding and in terms of hybrid functions, we get Substituting (58)–(60) in (53) yieldswith the following initial condition:where is the delay operational matrix described by (18). By solving the resulting algebraic system, we can find the vector asWith the use of (58), we would obtain the same value as the exact value of .
Example 2. Consider the following piecewise constant delay system [13]:The analytical solution to this problem is presented by [13]To solve this problem by the proposed approach, we first determine the value of . With the use of (22), we take . In addition, let be an arbitrary positive integer number. Definewhere and denote the approximate solution obtained by the present method and the exact solution, respectively. In Table 1, the maximum absolute error of corresponding to and different values of are summarized. This table shows that there is an excellent agreement between the analytical and approximate solutions. After choosing , small values for are needed to achieve a satisfactory approximation.

Example 3. As a more complicated problem, consider the following nonlinear piecewise constant delay system:The exact solution to this problem is described byAlthough the above problem is a nonlinear delay differential equation, the method developed in the current paper is applicable. To employ the procedure described in Section 4, we first select . Moreover, we choose . LetBy expressing , , and in terms of hybrid functions, we obtainSubstituting (73)–(75) in (67) implieswith the following initial condition:in which and is the operational matrix of delay given by (18). Moreover, and can be determined in a way similar to the construction method of the matrix presented by (17). By solving (77)(78), it follows that Using (73), we would obtain the same value as the exact value of .
6. Conclusion
An efficient procedure has been successfully developed for solving delay differential equations with a piecewise constant delay function. The method is based upon a hybrid of blockpulse functions and Legendre polynomials. The nice properties of the hybrid functions together with the associated operational matrices were used to convert the original problem into a system of algebraic equations. The proposed framework allows one to simultaneously make use of the best advantages of the two mentioned bases. After determining the appropriate value of , small values for are required to obtain an admissible approximation. It is worth noting that the correct choice of has a fundamental effect on the solution accuracy. It should be emphasized that the analytical solutions of Examples 1 and 3 cannot be derived solely either by blockpulse functions or by Legendre polynomials. The simulation results demonstrate the reliability and effectiveness of the proposed approximation scheme.
Competing Interests
The authors declare that they have no competing interests.
References
 M. MalekZavarei and M. Jamshidi, Time Delay Systems: Analysis, Optimization and Applications, NorthHolland Systems and Control Series, Elsevier Science, New York, NY, USA, 1987. View at: MathSciNet
 H. Smith, An Introduction to Delay Differential Equations with Applications to the Life Sciences, Texts in Applied Mathematics, Springer, Berlin, Germany, 2010. View at: Publisher Site  MathSciNet
 W. L. Chen and Y. P. Shih, “Shift Walsh matrix and delaydifferential equations,” IEEE Transactions on Automatic Control, vol. 23, no. 6, pp. 1023–1028, 1978. View at: Publisher Site  Google Scholar  MathSciNet
 W. L. Chen and B. S. Jeng, “Analysis of piecewise constant delay systems via blockpulse functions,” International Journal of Systems Science, vol. 12, no. 5, pp. 625–633, 1981. View at: Publisher Site  Google Scholar
 A. Bellen, “Onestep collocation for delay differential equations,” Journal of Computational and Applied Mathematics, vol. 10, no. 3, pp. 275–283, 1984. View at: Publisher Site  Google Scholar  MathSciNet
 A. Bellen and M. Zennaro, “Numerical solution of delay differential equations by uniform corrections to an implicit RungeKutta method,” Numerische Mathematik, vol. 47, no. 2, pp. 301–316, 1985. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Bellen and S. Maset, “Numerical solution of constant coefficient linear delay differential equations as abstract Cauchy problems,” Numerische Mathematik, vol. 84, no. 3, pp. 351–374, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. P. Shih, C. Hwang, and W. K. Chia, “Parameter estimation of delay systems via block pulse functions,” Journal of Dynamic Systems, Measurement and Control, vol. 102, no. 3, pp. 159–162, 1980. View at: Publisher Site  Google Scholar  MathSciNet
 F.C. Kung and H. Lee, “Solution and parameter estimation in linear timeinvariant delayed systems using laguerre polynomial expansion,” Journal of Dynamic Systems, Measurement and Control, vol. 105, no. 4, pp. 297–301, 1983. View at: Publisher Site  Google Scholar
 L. Lee and F.C. Kung, “Shifted Legendre series solution and parameter estimation of linear delayed systems,” International Journal of Systems Science, vol. 16, no. 10, pp. 1249–1256, 1985. View at: Publisher Site  Google Scholar
 I. R. Horng and J. H. Chou, “Analysis, parameter estimation and optimal control of timedelay systems via Chebyshev series,” International Journal of Control, vol. 41, no. 5, pp. 1221–1234, 1985. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. T. H. Baker, C. A. H. Paul, and D. R. Willé, “Issues in the numerical solution of evolutionary delay differential equations,” Advances in Computational Mathematics, vol. 3, no. 3, pp. 171–196, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. B. Datta and B. M. Mohan, Orthogonal Functions in Systems and Control, vol. 9 of Advanced Series in Electrical and Computer Engineering, World Scientific, Singapore, 1995.
 M. Razzaghi and H. R. Marzban, “A hybrid domain analysis for systems with delays in state and control,” Mathematical Problems in Engineering, vol. 7, no. 4, pp. 337–353, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Marzban and M. Razzaghi, “Solution of timevarying delay systems by hybrid functions,” Mathematics and Computers in Simulation, vol. 64, no. 6, pp. 597–607, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. R. Marzban and M. Razzaghi, “Optimal control of linear delay systems via hybrid of blockpulse and Legendre polynomials,” Journal of the Franklin Institute, vol. 341, no. 3, pp. 279–293, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Marzban and M. Razzaghi, “Analysis of timedelay systems via hybrid of blockpulse functions and Taylor series,” Journal of Vibration and Control, vol. 11, no. 12, pp. 1455–1468, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Marzban and M. Razzaghi, “Solution of multidelay systems using hybrid of blockpulse functions and Taylor series,” Journal of Sound and Vibration, vol. 292, no. 3–5, pp. 954–963, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Marzban and M. Shahsiah, “Solution of piecewise constant delay systems using hybrid of blockpulse and Chebyshev polynomials,” Optimal Control Applications & Methods, vol. 32, no. 6, pp. 647–659, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Scientific Computation, Springer, Berlin, Germany, 2006.
 H. R. Marzban and S. M. Hoseini, “Numerical treatment of nonlinear optimal control problems involving piecewise constant delay,” IMA Journal of Mathematical Control and Information, 2015. View at: Publisher Site  Google Scholar
 P. Lancaster, Theory of Matrices, Academic Press, New York, NY, USA, 1969. View at: MathSciNet
 E. L. Ortiz, “The tau method,” SIAM Journal on Numerical Analysis, vol. 6, no. 3, pp. 480–492, 1969. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2016 H. R. Marzban and S. M. Hoseini. 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.