Advances in Numerical Analysis
Volume 2016 (2016), Article ID 8903184, 9 pages
http://dx.doi.org/10.1155/2016/8903184

Research Article

Analysis of Linear Piecewise Constant Delay Systems Using a Hybrid Numerical Scheme

H. R. Marzban and S. M. Hoseini

Department of Mathematical Sciences, Isfahan University of Technology, Isfahan 8415683111, Iran

Received 30 April 2016; Accepted 3 July 2016

Academic Editor: Marc Gerritsma

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.

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 block-pulse 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 time-delay 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 well-known 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, time-delay 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 block-pulse 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 block-pulse functions and Legendre polynomials. Second, the current paper is an extension of our previous ones. More specifically, the time-delay 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 block-pulse functions or by Legendre polynomials. Combining block-pulse 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 block-pulse 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 block-pulse 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 block-pulse functions and Legendre polynomials, respectively. Here, are the well-known Legendre polynomials of order which are orthogonal on the interval and satisfy the following recursive formula [20]:Since consists of block-pulse functions and Legendre polynomials, which are both complete and orthogonal, the set of the hybrid of block-pulse 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 time-varying 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 well-known 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 time-varying 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.

Table 1: Results of maximum error for Example 2.

Present method

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 block-pulse 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 block-pulse 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. Malek-Zavarei and M. Jamshidi, Time Delay Systems: Analysis, Optimization and Applications, North-Holland Systems and Control Series, Elsevier Science, New York, NY, USA, 1987.

H. Smith, An Introduction to Delay Differential Equations with Applications to the Life Sciences, Texts in Applied Mathematics, Springer, Berlin, Germany, 2010.

W. L. Chen and Y. P. Shih, “Shift Walsh matrix and delay-differential equations,” IEEE Transactions on Automatic Control, vol. 23, no. 6, pp. 1023–1028, 1978.

W. L. Chen and B. S. Jeng, “Analysis of piecewise constant delay systems via block-pulse functions,” International Journal of Systems Science, vol. 12, no. 5, pp. 625–633, 1981.

A. Bellen, “One-step collocation for delay differential equations,” Journal of Computational and Applied Mathematics, vol. 10, no. 3, pp. 275–283, 1984.

A. Bellen and M. Zennaro, “Numerical solution of delay differential equations by uniform corrections to an implicit Runge-Kutta method,” Numerische Mathematik, vol. 47, no. 2, pp. 301–316, 1985.

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.

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.

F.-C. Kung and H. Lee, “Solution and parameter estimation in linear time-invariant delayed systems using laguerre polynomial expansion,” Journal of Dynamic Systems, Measurement and Control, vol. 105, no. 4, pp. 297–301, 1983.

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.

I. R. Horng and J. H. Chou, “Analysis, parameter estimation and optimal control of time-delay systems via Chebyshev series,” International Journal of Control, vol. 41, no. 5, pp. 1221–1234, 1985.

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.

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.

H. R. Marzban and M. Razzaghi, “Solution of time-varying delay systems by hybrid functions,” Mathematics and Computers in Simulation, vol. 64, no. 6, pp. 597–607, 2004.

H. R. Marzban and M. Razzaghi, “Optimal control of linear delay systems via hybrid of block-pulse and Legendre polynomials,” Journal of the Franklin Institute, vol. 341, no. 3, pp. 279–293, 2004.

H. R. Marzban and M. Razzaghi, “Analysis of time-delay systems via hybrid of block-pulse functions and Taylor series,” Journal of Vibration and Control, vol. 11, no. 12, pp. 1455–1468, 2005.

H. R. Marzban and M. Razzaghi, “Solution of multi-delay systems using hybrid of block-pulse functions and Taylor series,” Journal of Sound and Vibration, vol. 292, no. 3–5, pp. 954–963, 2006.

H. R. Marzban and M. Shahsiah, “Solution of piecewise constant delay systems using hybrid of block-pulse and Chebyshev polynomials,” Optimal Control Applications & Methods, vol. 32, no. 6, pp. 647–659, 2011.

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 non-linear optimal control problems involving piecewise constant delay,” IMA Journal of Mathematical Control and Information, 2015.

P. Lancaster, Theory of Matrices, Academic Press, New York, NY, USA, 1969.

E. L. Ortiz, “The tau method,” SIAM Journal on Numerical Analysis, vol. 6, no. 3, pp. 480–492, 1969.