1 Introduction

Time delays exist in a wide range of Engineering systems.Since time delay may improve the performance of the system,there are various applications of time delay in active vibration control and noise control[1–5].On the other hand,time delay may influence the stability and the dynamic behavior of the system[6–9],which deserves attention.Delay differential equations(DDEs)are employed to describe time-delay systems.They are a special class of differential equations,called functional differential equations.Thus,the approach to solving DDEs differs from that to solving ordinary differential equations.The main difficulty lies in the transcendental character of DDEs[10].

Extensive investigations have been carried out concerning the solution to DDE[10–16].Delay problems always bring about an infinite spectrum of frequencies.An approximation method of this type is the well-known Padé approximation,which results in a shortened repeating fraction for the approximation of the characteristic equation of the delay τ[11,12].The work of Lam[11]gave an algorithm for forming a balanced realization of the all-pass Padé approximants of e−τs.The algorithm generated the coefficients of the Padé denominators based on a scaled recurrence formula.The problem of finding reduced-order models for a class of delay systems was considered in Ref.[12].The delay element e−τs was suggested to be replaced by Padé approximants.

Another approach to solving DDE is to find the response through analyzing the entire delay spectrum instead of approximating the time delay.Several studies have been accomplished to find a solution to DDE by solving its characteristic equation under different conditions[13–15].A related study on analytic solution to linear DDEs can be found in Ref.[13].In their work,Wright[14]studied a Fourier-like analysis of the existence of the solution and its properties for the nonlinear DDEs.Similar approaches to linear DDEs were also reported by Bellman and Cooke[15].The uniqueness of the solution and its properties for the linear DDEs with varying coefficients were studied by W right[14].

Although different methods have been proposed to solve DDEs,they still have a few shortcomings.The solution to DDE can be investigated analytically by using the classical method of steps[16].However,the calculations will quickly become inconvenient without revealing essential properties of the solution.The readers can refer to Ref.[17]for the example of using the method of steps to solve DDEs.Since a few frequencies are actually dominant in the appearance of oscillations[10],it can be envisioned that an approximation in the frequency domain may be feasible.Therefore,it is of necessity to present a frequency method that guarantees the simplicity of calculation without sacrificing accuracy,which is the motivation of this paper.Jin et al.[18]studied the identification of time delay dynamic system in experiment,while this work can help to better understand the properties of delayed system.

This paper concerns solving a class of linear DDEs with constant coefficients.Firstly,the solution in the frequency domain is deduced and expressed by the transfer function,system parameters and equivalent initial conditions.The theoretical transfer function is represented by infinite number of poles and residues,which can be acquired by fitting the curve of the theoretical frequency response function(FRF).While the transfer function is the sum of infinite terms,it can be truncated by selecting the dominate poles and residues.Moreover,the guidance of selecting these dominant poles and residues is provided.

两组比较，学生对OSCE模式满意度较高，认为OSCE模式更符合教学大纲，在评分标准、难易程度、时间安排均比传统考核方式更合理，差异具有统计学意义（P＜0.05）（见表2）。

It is worth noting that the method presented by Qin[6]may lead to the opposite judgement of the stability of the delayed system when the term with time delay was expanded by Maclaurin series.To be more specific,after the delayed term was expanded,one may get the conclusion that the system was unstable while the system was actually stable.However,the method proposed in this paper will not result in the wrong judge of the stability.In other words,the truncation of the transfer function does not affect the stability of the original delayed system so that one will not make mistakes when judging the stability of the system.In addition,since the method is developed in the frequency domain,there will be fewer cumulative errors appearing in other iteration methods,which contributes to the accuracy of calculation.Another feature of the method is that usually a few truncation terms are sufficient for getting a reasonably accurate solution,leading to high efficiency.

The rest of this paper is organized as follows.Section 2 describes the category of the DDEs concerned in this paper.The frequency method to solve the DDE is proposed in Sect.3.Numerical examples of both stable and unstable systems are presented to illustrate the method in Sect.4.Conclusions are provided in Sect.5.

采样点的布置根据不同的工作比例尺，遵循相应的化探规范，要求每个采样点应最大限度地控制上游汇水域。在面积较小的残山区选择在残山脚部低洼处布置采样点，布置1～2个点能有效控制残山面积即可。在这样的景观区由于受残山的面积大小和密集程度的制约，因此采样点布设不遵循“均匀性和不能连续出现3个空白小格”的原则。

2 Problem formulation

The equation of motion of a time delay system with n degrees of-freedom(DOF)is given as

with initial conditions

where M,C,K ∈ Rn×n are the mass matrix,damping matrix,stiffness matrix,respectively.u(t)∈Rn is the vector of external excitation.G ∈ Rn×n is the feedback gain matrix,and τ is the time-delay.A ll the system parameters and feedback parameters are real constants.x(t)∈Rn is the vector of the displacement response of the time delay system.

The characteristic equation of Eq.(1)is

is the transfer function of the closed-loop system.

In this paper,we discuss the systems where poles are complex conjugate pairs as

whereαi is the real part of the pole is the imaginary unit.fi is the frequency in Hz.

The solution to Eq.(1)consists of two parts.The first part is the response of the system over the time interval of t ∈ [0,τ),while the second part is the response over the time interval of t∈ [τ,∞).In this paper,the system over the time interval of t∈ [0,τ)and t∈ [τ,∞)w ill be termed as open-loop system and closed-loop system,respectively.

The derivation of the proposed method can be developed for a general class of inputs.In this research we specifically focus on transient response to highlight the essential aspect of the method proposed.Response under external excitation can be derived in a similar manner.Therefore,this paper presents the derivation of the following DDE

3 Solution to time delay system

where x op(t)∈Rn is the vector of the displacement response of the open-loop system.The subscript•op refers to “open loop”.

with initial conditions

The equation of motion of the open-loop system is given as

Since it is straightforward to obtain the solution to Eq.(6),the process is omitted here[19].The displacement and velocity responses can be acquired and denoted as x op(t)and,respectively.Substituting t=τ into x op(t)andyields the responses at the instant t=τ,which are given as

andw ill be used to calculate the equivalent initial conditions later.So far,the solution to the open-loop system has been obtained.

3.1 Solution to closed-loop system

This section presents the process of solving for the solution to the closed-loop system in detail.The equation of motion of the closed-loop system is a DDE,which can be expressed as

where x cl(t)∈Rn is the vector of the displacement response of the closed-loop system.The subscript•cl refers to“closed loop”.

The equivalent initial conditions are denoted as

It is worth noting that x cl0 and are two unknown vectors and will be determined with continuity conditions later.Equivalent initial conditions are the responses of the closed loop system at the instant t=0.However,the domain of Eq.(9)is t ∈ [τ,∞),which means there is no response of the closed-loop system over the time interval t∈ [0,τ).Therefore,the equivalent initial conditions have no physical meaning.This is the reason why they are termed as equivalent initial conditions.

The solution to the closed-loop system can be expressed as

where L−1(·)is inverse Laplace transform and

where s represents the pole of the system.

The velocitycan be obtained by differentiating the displacement.Substituting t=τ into the displacement and velocity yields the responses at the instant t=τ,namelyandHowever,it is worth noting that the two vectors are expressions containing the unknown equivalent initial conditions given in Eq.(10).We need to calculate the equivalent initial conditions.

3.2 Equivalent initial condition

As stated in the Sect.2,the solution to the time delay system is calculated over two time intervals,which means the responses of the open-loop and closed-loop systems are respectively acquired.However,since the displacement and velocity of the two systems are continuous at the instant t=τ,we have the following continuity conditions

The left-hand side of Eq.(13)is the responses of the closed-loop system containing the unknown equivalent initial conditions x cl0 and,while the right-hand side is the responses of the open-loop system and have already been obtained in Eq.(8).

2.2.7 3组小鼠最大呼气中期时间比较 对照组和脂多糖组在6、18、36 h后的最大呼气中期时间比较，差异无统计学意义(P>0.05)。18、36 h后，甲强龙组最大呼气中期时间较脂多糖组下降，差异有统计学意义(P<0.05)；6 h后，两组最大呼气中期时间比较，差异无统计学意义(P>0.05)。

陶粒/碎石复合骨料透水混凝土自身具有吸滤水中污染物有害离子的作用，使混凝土对氮磷等离子的吸滤功能进一步提高。同时，两种混凝土通过水泥水化粘结为一体，增加结构整体性。养护完后的预制件灌注混有草籽、灌木籽、长效缓释肥料的种植土泥浆，定期维护，保证预制件现场装配前的草、体量较小的湿生灌木正常生长。

Once Eq.(13)is solved,complete solution to the closed loop system can be completely obtained by substituting the specific values of the equivalent initial conditions into Eq.(11).

本文主要将西药与中药瓜萎薤白桂枝汤加味联合治疗心肌缺血，结果显示，治疗组与对照组相比，临床有效率分别为90.0%和72.5%；总胆固醇、甘油三酯和低密度胆固醇水平降低，而高密度胆固醇水平升高。差异有统计学意义（P＜0.05）。

3.3 Truncation of transfer function

The transfer function given in Eq.(12)can be written in the form of a matrix as

Substituting Eq.(14)into Eq.(11)yields the solution to the closed-loop system,given by

where m j and cq j are the elements of the mass matrix and damping matrix,respectively.

It is obvious that with the increase of the truncation order,the accuracy of the approximate solution is improved.The convergence of the curves shown in Fig.12 reveals the truncation order does not have to be very high to get a quite satisfactory approximate solution.Therefore,the proposed method makes the calculation efficient without sacrificing accuracy.

where sr and are the pole and its complex conjugate pair.r Rij and r Rij∗are the residue and its complex conjugate pair.

Thus,Eq.(16)can be further rewritten in detail as

It can be seen from Eq.(17)that the transfer function is the sum of infinite terms.From the view of feasibility,the transfer function should be truncated and written in the form of finite items as

where N is the truncation order.

The truncated transfer function in Eq.(19)is the sum of finite terms expressed in terms of the poles and residues.Theoretically,there are numerous choices we can make to select the N pairs of poles to truncate the transfer function.Therefore,it is necessary to provide a rule of selecting the poles.

选取我院2016年1月～2017年12月医院中入住的82名急性心肌梗死（AMI）行急诊急诊冠状动脉介入治疗患者，男66名，女16名，年龄33～92岁。

As can be seen from the expressions of the solution(Eqs.(15)and(18)),the real parts of the poles are the powers of the exponents.This means the larger a real part is,the greater influence its corresponding pole will ultimately have on the solution.In other words,the poles with larger real parts dominate the solution.Consequently,the rule of selecting poles is based on the values of their real parts.

全面深化改革开放中“全面”是状语，“深化”成为动词了。a译文将“全面深化”翻译成定语是仓促草率之举。对于全面深化改革开放，这一标准译文应给我们的启示是，状语仍译成状语，是保险的策略。

Specifically,if the system is stable,we first sort the poles according to the values of their real parts(see Eq.(4)).The sequence rule is shown in the following inequality as

It is easy to understand that the approximate solution with higher truncation order will have higher accuracy.Nevertheless,in practical problems usually a small number of terms would be sufficient.

Considering an unstable system with U unstable poles,we still sort the poles according to the magnitudes of their real parts.The sequence rule is presented as

It can be seen that even if there is only one pair of poles chosen in the truncated transfer function,the stability of the original system will not be changed since the selected poles have positive real parts.Therefore,this method of truncation will not lead to the wrong judgment of the stability of the delayed system.

When the truncation order N is determined,the poles that satisfy the inequality(20)or(21)and their corresponding residues will be chosen to truncate the transfer function.Then,substituting the poles and the residues into Eqs.(15)and(18),the approximate solution to the closed-loop system can be obtained with the infinite signs∞replaced by the truncation order N,given by

where

Choosing truncation orders ranging from 1 to 10 and repeating the above procedures,we can obtain the corresponding approximate solutions.Here,the figures of the results are omitted due to the limit of space.Figure 8 shows the errors of the FFTs of these approximate solutions with different truncation orders(Figs.9,10).

So far,the solution to the original time delay system has been completely obtained since the solution to the open loop system and the one to the closed-loop system have been acquired.

4 Numerical simulation

Numerical solutions to DDEs can be obtained by XPPAUT[20],a common software used in the calculation of differential equations.Since there is no analytical solution to DDE,the solutions obtained by XPPAUT will be regarded as“accurate solutions”,which the approximate solutions will be compared with.

In this section,numerical simulations of both stable and unstable time delay systems are given to show the efficiency and accuracy of the proposed method.The equations of motion of both stable and unstable time delay systems are the same as Eq.(1).

4.1 Stable time delay system

In this case,the system parameters are

The feedback parameters are set as

在本研究中，选择了90例患者为主要研究对象，之后以患者最后的病理结果为基准，分析这三类诊断方式在宫颈癌治疗当中的价值后发现，经阴道超声联合弹性成像检查方式的灵敏度、特异度以及准确度均高于其他两组，P＜0.05差异具有统计学意义。

The initial conditions are given as

巷道岩爆控制，往往在掘进断面的辅助眼中，沿正中心及两侧接近拱腰部位呈三角形地各选取眼间距不超过2 m的1个辅助眼，实施掘进进尺2倍的超深凿岩，并装满药与辅助眼一同爆破[8]；而且在巷道拱腰部位的两帮，每个掘进进尺各凿一个深约2～2.5 m的震动爆破孔，仅孔底装药1/3卷(约40～50 g)并紧密堵塞20 cm长的黄泥，实施震动爆破[9]。一般巷道帮墙震动爆破的孔深为巷道宽度的70%～80%，巷道宽度较大时震动孔深度取下线，否则取上线；为了避免爆破震动损坏巷道壁面，一般震动孔深度不小于2 m[9]。

According to Eq.(12),the theoretical FRF matrix of the stable delayed system and that of the corresponding uncontrolled system can be obtained.Take one FRF element,say H1,1,as an example.The two FRF curves of the stable delayed system and the uncontrolled system are shown in Fig.1.

After using the vector fitting method[21]to fit the FRF curve of the delayed system,we can obtain as many poles as we want.The first 30 poles(15 pairs)are shown on the complex plane in Fig.2.

It can be seen from Fig.2 that the real parts of the poles are all negative since the delayed system is stable.Also,the poles are sorted by their imaginary parts.However,according to sequence(20),the poles should be sorted by their real parts.Table 1 shows the first 10 pairs of the poles sorted according to sequence(20).The circles with squares in Fig.2 represented the 10 pairs of poles listed in Table 1.

Fig.1 Curves of FRF H1,1 of the stable delayed and uncontrolled system

Fig.2 Distribution of the poles of the stable delayed system

Table1 The sorted poles according to real parts

Sequence Real part α Imaginary part f(Hz)Pole s 1−0.068 10.440 −0.068±65.593i 2−0.689 3.789 −0.689±23.810i 3−1.680 8.928 −1.680±56.097i 4−4.198 11.548 −4.198±72.561i 5−4.278 9.492 −4.278±59.641i 6−4.671 7.556 −4.671±47.474i 7−5.474 5.869 −5.474±36.878i 8−5.646 13.260 −5.646±83.314i 9−5.866 4.211 −5.866±26.458i 10 −5.924 2.519 −5.924±15.825i

It can be seen from Table 1 that the poles at the beginning of the sequence have real parts with larger values.As is shown in Eq.(23),the real parts of the poles are the powers of the exponents.Thus,the first three frequencies dominate in the approximate solution will be around 10.5,3.8,and 9.0Hz.

Fig.3 Accurate and approximate solutions of DOF 1(N=1)

4.1.1 Free decay

We consider the free decay of the stable time delay system,which means there is no external excitation.Choosing different truncation orders,we can get several groups of approximate solutions by the mentioned method.Three groups of results with truncation order N being 1,3,and 10 are picked out and analysed to show the validity of the method.

晶闸管的额定电压通常选取断态重复峰值电压UDRM和反向重复峰值电压URRM中较小的标值作为该器件的额定电压。晶闸管体的额定电流值通常采用其普通状态时一般电流的2-3倍。在桥式体整流的电压电路中晶闸管内两边承受的最高正反向电压的额定值都是晶闸管正常工作下的额定电压通常采取它的最高正反向平均电压的2-2.5倍。

As can be seen from the enveloping lines,the accurate solution(blue line)includes multiple frequencies,while the approximate solution(red line)includes only one frequency,which is especially obvious in the zoom-in view of the responses from 0.6 to 1.2s.During the last second,the tendency the red curve,which represents approximate solution coincides well with the blue one standing for accurate solution.However,the amplitudes of the red curve are larger than those of the blue one.Similar phenomena can be observed for DOF 2 and DOF 3.Therefore,the figures of comparisons are omitted here.

Now the truncation order is chosen as N=3 and Fig.4 shows the comparison between the accurate solution and the approximate solution of DOF 1.

According to Fig.4,the tendency of the approximate solution(red line)is consistent with that of the accurate solution(blue line).The zooms of the beginning and the last second show that the curve of the approximate solution coincides quite well with that of the accurate solution.The consistency of the two curves indicates approximation solution is capable of describing the dominant components in the accurate solution.Combining Figs.3and4,it can be obviously noticed that the approximate solution with N=3 is far more satisfactory than that with N=1.

With the truncation order N=1,Fig.3 shows the comparison between the accurate solution and the approximate solution of DOF 1.

Fig.4 Accurate and approximate solutions of DOF 1(N=3)

Now the truncation order is selected as N=10.The approximate solutions of DOF 1,2,and 3 are obtained and presented in the following.Besides,the fast Fourier transforms(FFTs)of the solutions are conducted to indicate the frequency components,as is shown in Figs.5,6,and 7.

On one hand,Figs.5a,6a,and 7a indicate that the approximation solutions of the three degrees of freedom are highly close to the accurate solutions.On the other hand,Figs.5b,6b,and 7b show the curves of the FFTs of the approximate solutions are in good agreement with those of the accurate ones.Take Fig.7b as an example.The three peaks in Fig.7b have the frequencies of 3.784,9.094,and 10.437Hz,respectively,which are consistent with the first three frequencies shown in Table 1.

The velocity can be obtained by differentiating the displacement.Then,the equivalent initial conditions can be determined by continuity conditions given in Eq.(13).Lastly,substituting the specific values of the equivalent initial conditions into Eq.(22)yields the approximation solution to the closed-loop system.

渠道的选择，应该严格按照霍童古镇规划的要求为前提，在整体品牌形象下进行渠道拓展工作。霍童古镇采用统一的品牌形象和VI识别系统：景区的标志、小品塑造、指示牌、宣传手册，在式样、颜色的选择上精心地加以区别，并且在购物、娱乐、交通、卫生、安全等方面形成统一的高质量服务保障体系。加强对古镇的宣传，建立多渠道销售模式，结合宁德市其他旅游资源的优势和影响力，借助“清新福建”品牌让霍童古镇旅游产品和服务走向全国。

It is obvious that with the increase of truncation order,the accuracy of the approximate solution is improved.Also,the error reaches a plateau when the truncation order increases to 3,which means that approximate results with the truncation orders higher than 3 are quite similar to those with N=3.The convergence of the curves reveals that the truncation order does not have to be very high to get an approximate solution with satisfactory accuracy.Based on this,the proposed method is proved to have convenience and efficiency in calculation without the sacrifice of accuracy.

4.1.2 Forced response

Now we consider the forced response of the stable time delay system.A sinusoidal excitation u1(t)=200sin(6πt)is applied on DOF 1 and no external excitation is applied on the other two degrees of freedom.Therefore,the frequency of the only external excitation is 3 Hz.

民族地区乡镇文化人才的培养，需要各级政府提高意识，加强认识，在人力、财力、政策上倾斜，将民族地区乡镇文化工作真正作为民生工程来打造，结合当地民族特色，人民群众的需求，完成特色的乡镇文化顺利传承。

Fig.5 Accurate and approximate solutions of DOF 1(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.6 Accurate and approximate solutions of DOF 2(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.7 Accurate and approximate solutions of DOF 3(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.8 Errors of the FFTs with different truncation orders

The process of obtaining approximate solution of the forced response is similar to that of the free response.The following shows the approximation results with the truncation order N=10.

It can be seen from the above figures that the approximation solutions of the three degrees of freedom are highly close to the accurate solutions both in time domain and in frequency domain.The three peaks in Fig.11b with the frequencies of 3.792,9.071,and 10.445Hz are consistent with the first three frequencies shown in Table 1,respectively.The peak with the frequency of 3.0Hz is caused by the only external sinusoidal excitation.

Choosing truncation orders ranging from 1 to 10 and repeating the above procedures,we can obtain the approximate solutions to the DDE.Figure 12 shows the errors of the FFTs of these approximate solutions with different truncation orders.

也不怪他疑心。刚才讲价钱的时候太爽快了也是一个原因。她匆匆下楼，那店员见她也下来了，顿了顿没说什么。她在门口却听见里面楼上楼下喊话。

The element of the transfer function matrix can be expressed as the sum of terms involving the poles and their residues as

4.2 Unstable time-delay system

In this case,we study an unstable time-delay system whose system parameters and initial conditions are the same as those of the stable delayed system(see Sect.4.1).The differences are the feedback parameters,which are set as

Since the system is unstable,no external force is applied.

According to Eq.(12),the theoretical FRF matrix of the unstable delayed system and that of the corresponding uncontrolled system can be obtained.Take one FRF element,say H1,1,as an example.The two FRF curves are shown in Fig.13.

After using the vector fitting method to fit the FRF curve of the delayed system,we can obtain as many poles as we want.The first 30 poles are shown on the complex plane in Fig.14.

It can be seen from Fig.14 that four pairs of poles(marked by crosses)have positive real parts which means the delayed system is unstable.According to sequence(21),the poles should be sorted by their real parts.Table 2 shows the first 10 pairs of the poles sorted according to sequence(21).The poles with squares in Fig.14 represented the10 pairs of poles listed in Table 2.

Fig.9 Accurate and approximate solutions of DOF 1(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.10 Accurate and approximate solutions of DOF 2(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.11 Accurate and approximate solutions of DOF 3(N=10).a Solutions to the DDE.b FFTs of the solutions

Fig.12 Errors of the FFTs with different truncation orders

Fig.13 The curves of FRF H1,1 of the unstable delayed and uncontrolled system

Fig.14 The distribution of the poles of the unstable delayed system

Table2 Sorted poles according to real parts

Sequence Real part α Imaginary part f(Hz) Pole s 1 2.307 10.416 2.307±65.448i 2 1.177 8.945 1.177±56.203i 3 0.724 11.715 0.724±73.609i 4 0.274 7.467 0.274±46.919i 5 −0.357 3.818 −0.357±23.992i 6 −0.485 5.826 −0.485±36.604i 7 −0.644 13.338 −0.644±83.800i 8 −0.889 2.500 −0.889±15.705i 9 −1.042 0.833 −1.042±5.236i 10 −1.215 4.181 −1.215±26.270i

Table 2 shows that the poles in the front of the sequence have the real parts with larger values.After choosing the finite poles with larger real parts,the frequencies of these poles will dominant in the appearance of oscillations.

The process of obtaining approximate solution of the unstable delayed system similar to that of the stable delayed system.Using the proposed method and choosing the truncation order as N=15,we can get the approximate solutions,shown in Fig.15.

Since this delayed system is unstable and will never be put into practice,it is unnecessary to compare the specific values of the solutions.Figure 15 illustrates that the proposed method will not mistakenly treat an unstable delayed system as a stable one.This means the stability of the delayed system will not be changed due to the truncation of the transfer function,enabling the method to ensure the correctness of stability.

Fig.15 Comparisons between accurate and approximate solutions with N=15.a Solutions of DOF 1.b Solutions of DOF 2.c Solutions of DOF 3

5 Conclusions

The aim of this paper is to provide a frequency method to obtain the solutions of a class of linear DDEs with constant coefficients.

Firstly,the equation of motion is given,which clearly describes the form of the DDE concerned in this paper.The solution to the DDE consists of two parts.One is the solution over t ∈ [0,τ),which is called an open-loop system.The other is the solution to the system over t ∈ [τ,∞),called a closed-loop system.

Secondly,the solution to the closed-loop system is expressed in terms of the transfer function.Since the theoretical transfer function consists of infinite poles,it can be truncated,which is the basic idea of this method.While there are infinite poles,only a few of them dominate the solution.Therefore,the guideline of selecting poles is provided according to their real parts.Poles with larger real parts and their corresponding residues are selected to truncate the transfer function.

Finally,numerical simulations are given where the approximate solutions obtained by the presented method are compared with the solutions obtained by a well-known software.As is shown in the results,only a few truncation items are needed to get a satisfactory approximate solution,which makes the method quite efficient.

Acknowledgements This work was supported by the National Natural Science Foundation of China(11272235).

References

1.Cai,G.P.,Lim,C.W.:Optimal tracking control of a flexible hubbeam system with time delay.Multibody Syst.Dyn.16,331–350(2006)

2.Xu,J.,Sun,Y.X.:Experimental studies on active control of a dynamic system via a time-delayed absorber.Acta Mech.Sinica 31,229–247(2015)

3.Olgac,N.,Holmhansen,B.T.:A novel active vibration absorption technique:delayed resonator.J.Sound Vib.176,93–104(1994)

4.Wang,Q.,Wang,Z.H.:Optimal feedback gains of a delayed proportional-derivative(PD)control for balancing an inverted pendulum.Acta Mech.Sinica 33,635–645(2017)

5.Sun,X.T.,Xu,J.,Fu,J.S.:The effect and design of time delay in feedback control for a nonlinear isolation system.Mech.Syst.Signal Process.87,206–217(2017)

6.Qin,Y.X.:The Motion Stability of Dynamical Systems with Time delay,2nd edn.Science Press,Beijing(1984).(In Chinese)

7.Jiang,S.Y.,Xu,J.,Yan,Y.:Stability and oscillations in a slow–fast flexible joint system with transformation delay.Acta Mech.Sinica 30,727–738(2014)

8.Wang,Z.H.,Hu,H.Y.:A modified averaging scheme with application to the secondary hopf bifurcation of a delayed van der pol oscillator.Acta Mech.Sinica 24,449–454(2008)

9.Chanda,A.,Fischer,A.,Eberhard,P.,et al.:Stability analysis of a thin-walled cylinder in turning operation using the semi discretization method.Acta Mech.Sinica 30,214–222(2014)

10.Asl,F.M.,Ulsoy,A.G.:Analysis of a system of linear delay differential equations.J.Dyn.Syst.T.ASME 125,215–223(2003)

11.Lam,J.:Balanced realization of pade approximants of e-st.IEEE Trans.Autom.Control 36,1096–1100(1991)

12.Lam,J.:Model reduction of delay systems using pade approximants.Int.J.Control 57,377–391(1993)

13.Falbo,C.E.:Analytic and numerical solutions to the delay differential equations.In:Joint Meeting of the Northern and Southern California Sections of the MAA,San Luis Obispo(1995)

14.W right,E.M.:The non-linear difference-differential equation.Q.J.Math.17,245–252(1946)

15.Bellman,R.,Cooke,K.L.:Differential-Difference Equations.Academic Press,New York(1963)

16.Zwillinger,D.:Handbook of Differential Equations.Academic Press,Boston(1989)

17.Erneux,T.:Applied Delay Differential Equations.Springer,New York(2009)

18.Jin,M.S.,Sun,Y.Q.,Song,H.W.,et al.:Experiment-based identification of time delays in linear systems.Acta Mech.Sinica 33,429–439(2017)

19.Caughey,T.K.:Classical normal modes in damped linear dynamic systems.J.Appl.Mech.T.ASME 32,583–588(1965)

20.Ermentrout,B.:Simulating,Analyzing,and Animating Dynamical Systems:A Guide to XPPAUT for Researchers and Students.Siam,Philadelphia(2002)

21.Gustavsen,B.,Semlyen,A.:Rational approximation of frequency domain responses by vector fitting.IEEE Trans.Power Deliv.14,1052–1061(1999)