Discussion: “Analysis of a System of Linear Delay Differential Equations” (, and , 2003, ASME J. Dyn. Syst., Meas., Control, 125, pp. 215–223) PUBLIC ACCESS

[+] Author and Article Information
Naci Zafer

Department of Mechanical Engineering, Eskisehir Osmangazi University, 26480 Eskisehir, Turkeynzafer@ogu.edu.tr

J. Dyn. Sys., Meas., Control 129(1), 121-122 (Jan 01, 2007) (2 pages) doi:10.1115/1.2428282 History:

This work presents a commentary of the article published by Asl and Ulsoy (2003, ASME J. Dyn. Syst., Meas., Control, 125, pp. 215–223). We show by an example that their method leads to inaccurate results and is therefore erroneous.


Consider a system of first order delay differential equations of the formDisplay Formula

where A and B are n×n matrices, and y is an n×1 vector.

The exact solution of Eq. 1 cannot, in general, be obtained. Special cases of Eq. 1 were considered by Chen et al. (1) and exact analytical solutions were obtained. In (2), Asl and Ulsoy offer a closed form solution to the general case in matrix form and compute the stability lobes numerically. Both studies in (1-2) are based on a solution of a transcendental equation expressed in terms of the Lambert function, which was first derived by Briggs in (3). Because the solution involves the Lambert function and it is given as a series, reaching an exact analytical stability bound for Eq. 1 is still a problem.

Employing the Laplace transform method, Eq. 1 leads toDisplay Formula

where Y(s) is the Laplace transform of y(t). The authors of (2) call the relationDisplay Formula
as the characteristic equation in matrix form of Eq. 1. Clearly, this equation can be written in the formDisplay Formula
From 4, the authors of (2) deduce thatDisplay Formula
where W() is the well-known Lambert function. For further properties of the Lambert function, see (3-4) and the references therein.

If the equality in 4 were a matrix identity, then passing from 4 to 5 would be trivial. However, the equality in 4 means thatDisplay Formula

which is derived from 2, with Y(s)=TePZ(s). In this case, one must have the rank of M less than n in order to have a nontrivial solution, which is equivalent to detM=0. In fact, this condition leads to a correct characteristic equation of 1 asDisplay Formula
Unfortunately, M is set equal to zero as a condition in (2). We note that the solution space of 6 becomes all of Rn in this case, i.e., every Y(s) is a solution. Therefore, the use of 5 may not give the correct answer in general, as illustrated numerically in the next section.

Let us reconsider the case study in (2) with


As stated in (2), the s values are computed as eigenvalues of (1T)W(Q)B, i.e.,Display Formula


In Fig. 1, the data points N=1T,KcKm, represented by a diamond sign, are obtained from 8 numerically for the principal branch of the Lambert function. An important observation is that these data points do not match with Fig. 9 of (2). Figure 1 also displays, by the solid curve, the actual stability lobe diagram predicted analytically (5).

Using certain data points from Fig. 1, one can easily verify that 8 is satisfied but 7 is not. It is easy to check that (19.9,0.193993) at s=180.520834i is such a point. On the other hand, it can also be shown that 7 holds but 8 does not for some data points from Fig. 1. For instance, this can be verified by using (14.245976,0.105) at s=157.321327i.

As a result, the method given by Asl and Ulsoy (2) results in s values which do not satisfy 7. Thus, the closed form solution of Eq. 1 given in (2) is not correct.

I would like to thank the anonymous referee for the valuable suggestions.

Copyright © 2007 by American Society of Mechanical Engineers
View article in PDF format.



Grahic Jump Location
Figure 1

Stability lobes diagram




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In