Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 403137, 11 pages

http://dx.doi.org/10.1155/2014/403137

Research Article

Design and Stability Analysis of Uncertain Networked Predictive Control Systems with Multiple Forward Channels

Hongbo Song

College of Information, Zhejiang University of Finance and Economics, Hangzhou 310018, China

Received 2 January 2014; Accepted 24 February 2014; Published 2 April 2014

Academic Editor: Huaicheng Yan

Copyright © 2014 Hongbo Song. 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 is concerned with the design and stability of networked predictive control for uncertain systems with multiple forward channels. The delays and packet dropouts are distributed such that the classic networked predictive control (NPC) needs modifications to be implemented. An improved control signal selection scheme with distributed prediction length is proposed to increase the prediction accuracy and hence achieve better control performance. Moreover, stability analysis results are obtained for both constant and random cases. Interestingly, it is shown that the stability of the closed-loop NPC system is not related to the distributed delays when they are constant and the system model is accurate. Finally, a two-axis milling machine example is given to illustrate the effectiveness of the proposed method.

1. Introduction

As modern control systems become more and more complex, traditional point-to-point control architecture is no longer suitable under certain circumstances. Meanwhile, the technology of computer network improved significantly in the past decades and a new networked architecture emerges and attracts increasing attention. This kind of control systems are called networked control systems (NCSs), in which the control loops are closed via a network, for example, Field bus, Ethernet, and Internet [1–3]. Major advantages of NCSs include reduced cost, easy installation and maintenance, and high efficiency. However, the insertion of a network into the control loops introduces some challenging problems for NCSs such as network-induced delay, packet dropout, and quantization [4–7]. Moreover, beside the control problem for NCSs, state estimation or filtering for NCSs also attracts much attention [8–10]. Thus, NCSs have become a hot research area in the control and signal processing communities.

Among these problems, network-induced delay and packet dropout are two major issues that degrade the system performance or even cause instability [1]. Up to date, both network-induced delay and packet dropout issues have received extensive research attention and fruitful results were obtained. To name a few, time delay method [11, 12], stochastic system method [13], switched system method [14], and robust control method [15] were proposed to deal with the delay. Switched system method [16] and stochastic system method [17] were presented to handle the packet dropout. These methods commonly model the delay or/and packet dropout into the closed-loop NCS and the corresponding controller were then designed based on conditions that make the closed-loop NCS stable. Most of the conditions are only sufficient and thus the design is conservative. Moreover, essentially, they all passively compensate for the negative effects of the delay and packet dropout after accepting them. Alternatively, a novel actively compensating method called networked predictive control (NPC) was presented in [18, 19]. In NPC, future control sequences are generated and transmitted in a packet through the network. At the plant side, appropriate control signals from the sequences are selected to control the plant based on delay measurements.

With accurate system model and delay measurements, NPC was shown to have the ability to achieve desired control performance and the stability of the closed-loop NPC system is not related to the delay and packet dropout in the constant case [20]. In [21], an event driven NPC was presented to avoid the practical problem that delay measurements are inaccurate. In [22], a switched controller structure was proposed and the controllers were designed according to the delay based on switched system method. For more research results on NPC, see [23–28] and the references therein. However, two issues of NPC are not fully considered in the existing literature, which are the model uncertainty and multiple communication channels. The NPC for uncertain system with multiple distributed delays and packet dropouts in the feedback channels was studied in [29]. In [30], stability analysis was carried out for NPC systems with model uncertainty. In [31], a data reconstruction method was presented for NPC system with distributed delays in the feedback channels.

For a practical system, it is very difficult to obtain an accurate model. A nonlinear system is always linearized to be simplified. Thus, model uncertainty is an important issue to be investigated. On the other hand, when the plant to be controlled is spatially distributed, the corresponding NCSs have to be multiple communication channels. Obviously, for a NPC system with model uncertainty, larger input delay results in longer prediction length and thus leads to larger prediction error. Hence, the control performance becomes worse. In the case of multiple communication channels, delays and packet dropouts are distributed, which means that delays and packet dropout process may be different for each channel, while such case was not taken into account in classic NPC such as [18–20]. Introducing the queuing method in [32] can make the input delays the same for each channel, which makes the classic NPC applicable for the case without modification. However, input delays are actually enlarged in this way and the control performance will hence be degraded. Thus, this paper presents a modified control signal selection scheme such that the control signals which are predicted with distributed prediction length can be applied to control the plant. Such treatment essentially uses the most recent data and hence reduces the prediction length and improves the system performance. Only multiple forward channels are considered since the compensation scheme in the forward channels is different from the one in the feedback channels.

2. Design of NPC for Systems with Multiple Forward Channels

2.1. Structure of the NPC System

The structure of NPC systems with multiple forward channels is shown in Figure 1. The feedback channels are assumed to be ideal such that the plant outputs are transmitted to the predictive controller (PC) side without any delays or packet dropouts. The task of PC is to generate future control prediction sequences and transmit them to the control signal selector (CSS). The sequences are subjected to distributed delays and packet dropouts effects since the system has multiple forward channels. CSS receives the distributed prediction sequences and select proper control signal from them according to the measured input delays. The selected control signal is then used to control the plant.

Figure 1: Structure of NPC systems with multiple forward channels.

The control performance of NPC systems highly depends on the prediction accuracy, which is related to both the discrepancy between the adopted model and real system dynamic and the prediction length. On the other hand, larger input delay results in longer prediction length. Based on the NPC in [19, 30], it is known that larger modeling error and larger delay lead to larger prediction errors and hence degrade the control performance. Thus, intuitively, it can be inferred that two ways to improve the control performance of NPC systems are to reduce the modeling error and to make the prediction length shorter.

For the NPC systems with multiple forward channels, the delays and packet dropouts process are distributed. If we use the classic NPC, which assumes single communication channel, the input delays of all the channels should be the same and equal to a value which is the largest one of the distributed input delays. This can be done by the queuing method. However, with this treatment it can be seen that the prediction length is enlarged for the channels except the one with largest input delays. Based on the thought that shorter prediction length leads to better prediction accuracy, another way is to use an improved CSS with distributed prediction length, which makes the prediction length shorter than using classic NPC. This is the main idea of the paper.

The plant is represented by the following uncertain linear discrete-time state space system model:

where , , and are the state, input, and output vectors of the plant, respectively, and , , and are the nominal system matrices with compatible dimensions. and represent the system uncertainties and it is assumed that and satisfy the following structure:

where , and , and are matrices with compatible dimensions.

The prediction sequences are designed based on the following standard state observer:

where is the state of the observer, and is the observer gain matrix which can be designed by standard methods such as pole placement. The state feedback controller is as follows:

where is the controller gain matrix which can be designed by Lyapunov method and so on.

2.2. Assumptions and Notations

Before proceeding further, we need to introduce some reasonable assumptions and some notations.

(1)Without loss of generality, it is assumed that there are forward channels in the system and they are denoted by channel to channel , respectively.   are transmitted through channel , , respectively.(2)The distributed delays are bounded. The upper bound for the delays is denoted by and the lower bound is assumed to be step without loss of generality.(3)The number of consecutive packet dropouts for channels to is bounded and the upper bound is denoted by . The lower bound is assumed to be steps, which means no packet dropout.(4)The transmitted packets are time-stamped and the clocks of the receiving sides and the transmitting sides are synchronized such that the input delays can be calculated.

At any time instant, it is probable that more than one packet are received or no packet is received because of packet dropouts or packet disorder. In both cases, the data in the most recent packet will be used. The input delays then can always be calculated by the time stamp.

Since the case random delays and packet dropouts are considered, the input delay for channel at instant is denoted by and it is clear that the lower and upper bound for the input delays are and , respectively, where . It can be seen that input delays may be of the same value for certain channels. That is to say, , . Thus, we denote the number of different values for the input delays as . The corresponding values for the input delays are denoted by , respectively, where for any , . The set of the indexes of the channels in which is denoted by . Introduce selection matrices , , which is defined by

where represents a diagonal matrix and

An example is given to make these notations clear. Assume that the system has forward channels, which means . At instant , it is measured that and . Then, it can be seen that there are two different values for the distributed delays; that is, and , . Furthermore, we have and , respectively.

For the constant delays case, the notations can be simply used by removing the index of the ones for the random delays and packet dropouts case.

2.3. Predictive Controller

At the PC side, the state of the observer is obtained at instant according to (3). Since future system outputs and are not available at instant , the prediction sequences are generated as follows:

The prediction is up to steps to meet the worst case of delays and packet dropouts, that is, steps delays and consecutive packet dropouts. For the prediction procedure (6), it can be seen that the prediction accuracy at each instant depends on two factors. The first is how accurate represents the system dynamic. The second is how close to ,  . It is interesting to see that are fundamentally different for the random and constant cases.

For the constant delays case, even though the future control input is not applied, it can be inferred what control signal will be applied at instant and the control signal is already available at the PC side. That is to say, PC can use instead of in (6). For the random case, PC cannot determine what control signal will be used at instant ,  . Thus, can be approximated by (4) such that

This will also lead to different techniques used for the stability analysis of random and constant cases.

In classic NPC, the future control sequence is transmitted in a packet to the CSS as follows:

However, this paper considers the multiple forward channels case and thus the packet that contains ,  , is transmitted separately via channels. The packets transmitted through channel contains the following information, respectively. Consider

2.4. Control Signal Selector

Subjected to the distributed delays and packet dropouts effects, the packets are received by the CSS in a distributed delayed manner. Specifically, the following data are received by CSS at instant :

As mentioned, classic NPC with using queuing method can be used and in this case the selected control signal is

where .

Alternatively, a modified CSS scheme is presented here. We can see that for channel ,  , the predicted control signal with minimal prediction length is

Then, it follows by using the notations in the second subsection of this section that

It can be seen that the modified CSS selects control signal with distributed prediction length. Then, we can see that the following control signal to control the plant:

As a special case, if , , , which means that the length of the delays are all the same for each channel, then it can be seen that the control input (14) equals the control input (11).

3. Stability Analysis of the Closed-Loop NPC System

The last section designed a modified NPC strategy to control uncertain systems with multiple forward channels. Since stability is crucial for a control system, the stability of the closed-loop NPC system is studied in this section. Two cases are considered, which are the random distributed delays and packet dropouts case and the constant distributed delays case. As mentioned in the above section, the analysis techniques are different for the two cases.

3.1. Model of the Closed-Loop NPC System: Random Case

By (6) and (7), it follows that

Then, by (14) and (15) we have

Define an augmented vector to be

then the closed-loop NPC systems with random distributed delays and packet dropouts can be obtained by (1), (3), and (16) as follows:

where

Note that the dynamics of the closed-loop NPC system (18) is related to the input delays vector ; denote the set of all the possible values of as . When takes a value in , system (18) resides in the corresponding subsystem, which means that system (18) is a switched system. For simplicity of notation, introduce a one-to-one mapping that maps the set of to a set with numbers. For example, is mapped to 1, is mapped to 2, and so on. Denote the set of the numbers by . Moreover, the system has model uncertainties. Thus, system (18) can be transformed into the following uncertain switched system:

where

and and are and with is mapped by , .

3.2. Stability Results: Random Case

The closed-loop NPC system with random distributed delays and packet dropouts in the multiple forward channels is modeled as the switched uncertain system (20). A sufficient condition for the stability of the closed-loop NPC system (20) is presented in this subsection. A lemma is first given as follows.

Lemma 1. For given appropriate matrices , , and , with ,

holds for all if and only if there exists a scalar such that

Theorem 2. For the NPC system (20) with random distributed delays and packet dropouts, given all possible values of , and controller gain matrix and observer gain matrix , if there exists a matrix and a scalar such that the following matrix inequalities

hold, then the closed-loop NPC system (20) is stable.

Proof. It can be obtained by (2) that has the structure , where

By the stability result of switched system with arbitrary switching in [10], if there exists a matrix such that the following

holds, then the closed-loop NPC system (20) is stable. Then, by Schur complement and some matrix operations, we have

Then, by Lemma 1 it can be obtained that

Pre- and postmultiplying of (28) by and letting lead to (24). The proof is completed.

It should be pointed out that (24) is linear in matrix and scalar and thus can be conveniently solved by LMI toolbox in Matlab for example. This means that the stability of the closed-loop NPC system with random distributed delays and packet dropouts can be readily checked. However, a common Lyapunov function is used in the proof, which leads to some conservatism. One possible way to reduce conservatism is to use multiple Lyapunov function method.

3.3. Stability Results: Constant Case

In this section, the constant case is considered. Different from the random case, a necessary and sufficient condition is obtained. Define the state error as

Subtracting (3) from (1) leads to the following state error equation:

For an integer , it can be seen by (6) that the following two equalities hold:

By shifting (3) backward for steps, it can be obtained that

Then, by (33) and subtracting (31) from (32) we have the following:

Applying (34) recursively results in the following equation that holds for any integer :

We now have by (14) and (35)

For convenience of the representation of the closed-loop NPC system, is transformed as follows:

where ,  . Substituting in (1) and (30), respectively, by (26) yields

Define an augmented vector as

Then, the closed-loop NPC system with constant distributed delays is as follows:

where

It is clear that system (40) is a standard uncertain system and the stability result is as follows.

Theorem 3. For the NPC system (40) with constant distributed delays, given , and controller and observer , if and only if there exists a matrix and a scalar such that the following matrix inequality

holds, then the closed-loop NPC system (40) is stable, where

Proof. It can be obtained by (2) that has the structure . Following the robust control system results such as [33], it can be obtained that system (40) is stable if and only if there exists a matrix such that the following inequality holds:

Then, similar to the procedure in Theorem 2, the stability result can be obtained readily. The rest of the proof is thus omitted.

Assume that there are no model uncertainties; that is, and , and then closed-loop NPC system (40) with constant distributed delays can be represented in a more concise form as follows:

where

Clearly, a necessary and sufficient condition for the stability of the system (45) is the eigenvalues of and that are within the unit circle, which means that the stability of the closed-loop NPC system (45) is not related to the distributed delays. This extends the results in [20] to the multiple forward channels case.

4. Illustrative Example

To illustrate the effectiveness of the modified NPC with multiple forward channels, we consider the two-axis example of a three-axis milling machine. More details about the example can be referred to [34]. The parameters of the nominal model and the uncertain parts for the system are chosen as follows:

where ,  , , and . The sampling period is chosen as  s. Clearly, . The observer gain matrix and state feedback control gain matrix are calculated by pole placement method and Lyapunov method, respectively. They are given as follows:

In the simulation, the system model for the plant and the initial condition are chosen to be

We can see from the structure of the system matrices that and are independent of each other. It is assumed that there are forward channels for the system and that and are transmitted via channels and , respectively.

First, the effects of NPC for uncertain NCSs with delays in the forward channels are considered. Let the delays in both channel 1 and channel 2 be equal and constant and see the system performance. Take for example, as shown in Figure 2; the system performance becomes worse with larger delay and the system is unstable with more than 6 steps delays. The trajectory of with NPC is shown in Figure 3, from which we can see that the delays are effectively compensated. By applying Theorem 3, the closed-loop NPC system is still stable with 10 steps delay.

Figure 2: Trajectory of without NPC compensation.

Figure 3: Trajectory of with NPC compensation.

Then, we consider the effectiveness of the modified NPC with distributed delays and packet dropouts. In this simulation of constant case, the delays in channel and channel are chosen to be steps and steps, respectively. By Theorem 3 it can be obtained that the corresponding closed-loop system is stable. From the analysis, it can be inferred that the trajectory of will be the same by using modified and classic NPC since the delay in channel 1 is always the largest. While the performance of will be better by using modified NPC than classic one. Figure 4 show the trajectory of with modified and classic NPC, supporting the theory.

Figure 4: Comparison of modified and classic NPC with distributed constant delays.

Finally, the case of random distributed delays and packet dropouts is simulated. and are chosen to be and , respectively. That is to say, the delays in channels and are both random between step and steps, and the number of consecutive packet dropouts in both channels is up to steps. By Theorem 2, it follows that the closed-loop system is stable. The simulated distributed random delays and packet dropouts process are shown in Figures 5 and 6, respectively. The corresponding trajectories of and are shown in Figures 7 and 8, respectively. Clearly, we can see that the modified NPC achieves better performance than classic NPC method for this example.

Figure 5: Random distributed delays.

Figure 6: Random distributed packet dropouts process.

Figure 7: Trajectory of with modified and classic NPC with distributed random delays and packet dropouts.

Figure 8: Trajectory of of modified and classic NPC with distributed random delays and packet dropouts.

5. Conclusions

This paper studied the design and stability analysis of uncertain networked predictive control systems with distributed delays and packet dropouts in the forward channels. A modified NPC was proposed, in which the key point is an improved control signal selection scheme. The CSS with distributed prediction length uses the most recent data and hence can make modified NPC achieve better control performance. Stability analysis results are obtained for both constant and random cases. They are formulated as linear matrix inequalities and can be readily checked. Moreover, it is shown that the stability of the closed-loop NPC system is not related to the distributed delays when they are constant and system model is accurate. An example was given to show that the modified NPC method achieves better performance than classic NPC in the case of distributed delays and packet dropouts in the forward channels.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research work was supported by the National Nature Science Foundation of China (Grant no. 61304040) and the Nature Science Foundation of Zhejiang Province (Grant no. LQ13F030005).

References

W. Zhang, M. S. Branicky, and S. M. Philips, “Stability of networked control systems,” IEEE Control Systems Magazine, vol. 21, no. 1, pp. 84–99, 2001.

L. X. Zhang, H. J. Gao, and O. Kaynak, “Network-induced constraints in networked control systems—a survey,” IEEE Transactions on Industrial Informatics, vol. 9, no. 1, pp. 403–416, 2013.

H. Gao, T. Chen, and J. Lam, “A new delay system approach to network-based control,” Automatica, vol. 44, no. 1, pp. 39–52, 2008.

H. Yan, H. Shi, H. Zhang, and F. Yang, “Quantized control for networked systems with communication constraints,” Asian Journal of Control, vol. 15, no. 5, pp. 1468–1476, 2013.

J. P. Hespanha, P. Naghshtabrizi, and Y. Xu, “A survey of recent results in networked control systems,” Proceedings of the IEEE, vol. 95, no. 1, pp. 138–172, 2007.

Y.-G. Sun and Q.-Z. Gao, “Stability and stabilization of networked control system with forward and backward random time delays,” Mathematical Problems in Engineering, vol. 2012, Article ID 834643, 15 pages, 2012.

H. Zhang, H. Yan, F. Yang, and Q. Chen, “Quantized control design for impulsive fuzzy networked systems,” IEEE Transactions on Fuzzy Systems, vol. 19, no. 6, pp. 1153–1162, 2011.

D. Zhang, Q. G. Wang, L. Yu, and Q. K. Shao, “ filtering for networked systems with multiple time-varying transmissions and random packet dropouts,” IEEE Transactions on Industrial Informatics, vol. 9, no. 3, pp. 1705–1716, 2013.

B. Chen, L. Yu, and W.-A. Zhang, “Robust Kalman filtering for uncertain state delay systems with random observation delays and missing measurements,” IET Control Theory & Applications, vol. 5, no. 17, pp. 1945–1954, 2011.

H. Zhang, Q. Chen, H. Yan, and J. Liu, “Robust filtering for switched stochastic system with missing measurements,” IEEE Transactions on Signal Processing, vol. 57, no. 9, pp. 3466–3474, 2009.

D. Yue, E. Tian, and Q.-L. Han, “A delay system method for designing event-triggered controllers of networked control systems,” IEEE Transactions on Automatic Control, vol. 58, no. 2, pp. 475–481, 2013.

H. Yan, Z. Su, H. Zhang, and F. Yang, “Observer-based control for discrete-time stochastic systems with quantisation and random communication delays,” IET Control Theory & Applications, vol. 7, no. 3, pp. 372–379, 2013.

L. Zhang, Y. Shi, T. Chen, and B. Huang, “A new method for stabilization of networked control systems with random delays,” IEEE Transactions on Automatic Control, vol. 50, no. 8, pp. 1177–1181, 2005.

M. C. F. Donkers, W. P. M. H. Heemels, N. van de Wouw, and L. Hetel, “Stability analysis of networked control systems using a switched linear systems approach,” IEEE Transactions on Automatic Control, vol. 56, no. 9, pp. 2101–2115, 2011.

W.-A. Zhang and L. Yu, “A robust control approach to stabilization of networked control systems with time-varying delays,” Automatica, vol. 45, no. 10, pp. 2440–2445, 2009.

M. Yu, L. Wang, T. Chu, and G. Xie, “Stabilization of networked control systems with data packet dropout and network delays via switching system approach,” in Proceedings of the 43rd IEEE Conference on Decision and Control (CDC '04), pp. 3539–3544, Paradise Island, Bahamas, December 2004.

P. Seiler and R. Sengupta, “An approach to networked control,” IEEE Transactions on Automatic Control, vol. 50, no. 3, pp. 356–364, 2005.

G. P. Liu, J. X. Mu, D. Rees, and S. C. Chai, “Design and stability analysis of networked control systems with random communication time delay using the modified MPC,” International Journal of Control, vol. 79, no. 4, pp. 288–297, 2006.

G.-P. Liu, Y. Xia, J. Chen, D. Rees, and W. Hu, “Networked predictive control of systems with random network delays in both forward and feedback channels,” IEEE Transactions on Industrial Electronics, vol. 54, no. 3, pp. 1282–1297, 2007.

G. P. Liu, “Predictive controller design of networked systems with communication delays and data loss,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 6, pp. 481–485, 2010.

W. Hu, G. Liu, and D. Rees, “Event-driven networked predictive control,” IEEE Transactions on Industrial Electronics, vol. 54, no. 3, pp. 1603–1613, 2007.

R. Wang, G.-P. Liu, W. Wang, D. Rees, and Y.-B. Zhao, “ control for networked predictive control systems based on the switched Lyapunov function method,” IEEE Transactions on Industrial Electronics, vol. 57, no. 10, pp. 3565–3571, 2010.

J. H. Zhang, Y. Q. Xia, and P. Shi, “Design and stability analysis of networked predictive control systems,” IEEE Transactions on Control Systems and Technology, vol. 21, no. 4, pp. 1495–1501, 2013.

Y. Xia, W. Xie, B. Liu, and X. Wang, “Data-driven predictive control for networked control systems,” Information Sciences, vol. 235, pp. 45–54, 2013.

G. Pin and T. Parisini, “Networked predictive control of uncertain constrained nonlinear systems: recursive feasibility and input-to-state stability analysis,” IEEE Transactions on Automatic Control, vol. 56, no. 1, pp. 72–87, 2011.

A. Onat, T. Naskali, E. Parlakay, and O. Mutluer, “Control over imperfect networks: model-based predictive networked control systems,” IEEE Transactions on Industrial Electronics, vol. 58, no. 3, pp. 905–913, 2011.

Y.-B. Zhao, J. Kim, and G.-P. Liu, “Error bounded sensing for packet-based networked control systems,” IEEE Transactions on Industrial Electronics, vol. 58, no. 5, pp. 1980–1989, 2011.

Y.-B. Zhao, J. Kim, X.-M. Sun, and G.-P. Liu, “On the delay effects of different channels in Internet-based networked control systems,” International Journal of Systems Science, vol. 44, no. 11, pp. 2119–2129, 2013.

H. B. Song, G. P. Liu, and L. Yu, “Networked predictive control of uncertain systems with multiple feedback channels,” IEEE Transactions on Industrial Electronics, vol. 60, no. 11, pp. 5228–5238, 2013.

G.-P. Liu, “Analysis of networked predictive control systems with uncertainties,” in Proceedings of the IEEE International Conference on Systems, Man and Cybernetics (SMC '09), pp. 4819–4824, San Antonio, Tex, USA, October 2009.

Y. B. Zhao, J. Kim, G. H. Yang, and G. P. Liu, “Model-based compensation for multi-packet transmission in networked control systems,” in Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference, pp. 3136–3141, Orlando, Fla, USA, December 2011.

R. Luck and A. Ray, “An observer-based compensator for distributed delays,” Automatica, vol. 26, no. 5, pp. 903–908, 1990.

S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, vol. 15 of SIAM Studies in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1994.

F.-L. Lian, J. Moyne, and D. Tilbury, “Modelling and optimal controller design of networked control systems with multiple delays,” International Journal of Control, vol. 76, no. 6, pp. 591–606, 2003.