Research Article  Open Access
Controller Design for Asynchronous Multirate SampledData Systems
Abstract
This paper considers the analysis and synthesis of a linear discrete asynchronous multirate sampleddata system. An controller based on an observer is proposed, which guarantees the stability of the closed system and makes the norm of the closed system less than a given attenuation level. To improve the performance further, a tradeoff strategy is applied. That is, the exogenous signals sampled at different rates are lifted to an appropriate signal rate, while the endogenous signals are not lifted for avoiding the causal constraint and the dimension multiplied again. The controller is obtained by solving the corresponding matrix inequality, which can be calculated by Matlab. Finally, an example is presented to demonstrate the validity of these methods.
1. Introduction
Sampleddata control systems, which provide digital control techniques, and emerged with the development of computers, have been long studied by many researchers [1, 2]. Multirate sampleddata control systems are required when the signals of interest are sampled at different rates compared with others. There are two primary reasons for this requirement. First, it has been shown that multirate sampling can improve the closedloop control performance significantly. Second, technological or economic constraints in some real applications necessitate the use of control schemes where sensor measurements and control inputs have to be performed at different sampling rates. Figure 1 shows an ideal model of a standard multirate sampleddata control system. For example, in hybrid electric vehicles, a hierarchical controller is designed to decide the torque demands of motor, generator, internal combustion engine, and mechanical brake according to the driver’s torque demand, speed of vehicle, and battery’s state of charge (SOC), where the SOC is estimated by battery management system (BMS) and the speed of vehicle is fed by sensor [3]. As a largescale and complex system, the sensors and control inputs of subsystems are difficult to be performed at the same sampling rates; thus a multirate sampleddata system will be a better choice for modeling.
Multirate systems have been studied since the late 1950s, resulting in the development and application of various sampling methods. Classical approaches include the lifting method, control [4, 5], LMI method [6], prediction control, LQG control, iterative control [7], and parameter estimation [8]. Among these, the lifting method [9] is a vital tool for studying multirate systems, whereby such systems can be changed into singlerate systems. This method can improve the performance of a plant to some extent [10]. However, during the lifting process the dimensions of the system will be multiplied and it is increasingly difficult to analyze the new system effectively, especially systems with very different sampling rates, which can lead to the “model explosion” problem. control is a useful method for analyzing systems and for synthesizing control problems based on the lifting technique, and LMI tools can be used to solve the inequalities [11]. By the way in real applications, fastinput and slowoutput is a good way of improving the performance of the system, and this mechanism and its results have been discussed in previous studies [12, 13].
Most of lifting methods require that all signals should be sampled and held during the same time slot [1, 2, 9, 10, 14, 15]. That is, at a specific , where , , are unit matrices and , , and , will be defined in Section 2. However, in most actual situations, this condition cannot be satisfied because the overall system lacks a slot when the entire signal can be sampled or held, and there is a class of multirate systems that each output of them has its own frequency of measurement and each input has its own frequency of updating [16–19]. A system with this type of sampling mechanism is referred to as asynchronous system in this paper. The traditional lifting technique is not applicable for the sampling frequencies of the sampling and holding elements are incommensurate [4, 20].
LQG control based on an observer is an effective method for reducing the dimensions of a multirate system without using the lifting technique [14, 21]. An LQG optimal controller derived is periodic and can guarantee the overall closedloop stability. To compensate for the skewed input signal caused by slowoutput feedback, only the input signal is supplemented by the observer. However, in the previous results, the controller and observer are designed separately, which may result in conservativeness.
Based on the methods mentioned above, the design of a controller for a multirate sampleddata system that does not multiply the dimensions necessitates the synthesis of various appropriate measures. In this paper, we propose an observerbased method for analysis and synthesis of linear discrete asynchronous multirate sampleddata systems. In contrast to the traditional lifting method, the control inputs and measured outputs are not lifted to avoid the dimension multiplication and causal constraints. However, the exogenous system signals are lifted rationally to change the system into a fastinputslowoutput system. The controller can be obtained by solving an inequality with LMI tools of Matlab. Although the inequality is nonlinear, it can be calculated linearly under a constraint [22]. The approach proposed in this paper is different from the observerbased LQG method in [21]. Here, the controller and observer are obtained simultaneously by solving a matrix inequality under the same optimal performance constraint, thus reducing the conservativeness and improving the system performance. In addition, the concept of asynchronism in this paper is different from those in [23, 24], and so forth, where the sampling at sensors is assumed to be asynchronized with timevarying delay and it falls into the framework of timedelay systems, not the multirate sampleddata systems.
The remainder of this paper is organized as follows. Section 2 provides a description of the multirate system and its background, including lemmas and definitions. In Section 3, the main results obtained using the synthesis methods are presented. In Section 4, a simulated numerical example is provided to demonstrate the effectiveness of the method. Finally, Section 5 concludes the paper.
2. System Descriptions
Definition 1. A multirate sampleddata system is called an asynchronous system if the signals from its sampler or holder cannot be obtained simultaneously during any time slot.
For example, as shown in Figure 2, Signal 1 and Signal 2 are the signals of the sampler and the holder with period 3, respectively, while Signal 3 is the signal with period 2. Thus, they cannot be obtained simultaneously during any time slot because the odd period signals appear alternately. In Figure 2, the filled circles indicate that the signal is sampled or held in the slot, whereas the hollow cylinders indicate that the signal is not sampled or held. Therefore, this system cannot satisfy the assumption mentioned in most previous studies where all signals can be sampled or held during the same time slot.
We consider a linear discrete asynchronous multirate sampleddata system: where , , , , , , , is a sampling device, is the output measured vector, and is the initial state.
Assume that the th component , , of the input vector can be modified every time instant . Then, is the output of the following discretetime period system, which is defined as the inputholding mechanism: where where is a new state variable and is a new input variable. The integer () describes the skew of the inputsholding mechanism.
Then, matrix and system (2) have a period of , and
Similarly, the output sampling mechanism is given as where where () indicates the skew of the outputs mechanism.
Matrix has a period, and
After combining (2) and (3), we obtain a new samplingholding system: where
Based on the of matrix and the of matrix , system (9) has the period of
The state observer of system (2) is where asymptotically stabilizes system (12).
The observerbased controller is chosen as where .
The state error is defined as
According to (2), (12) and (14), the following equation can be obtained:
Combining (9), (13), and (15) yields where
Lemma 2 (see [21, 25] ). If(i)the pair is stabilizable,(ii)the pair is stabilizable,(iii)the pair is detectable,(iv)for any pair of distinct eigenvalues of , and , , , it follows that ,(v)there is no eigenvalue of , , such that ,then the pairs and are controllable, and the pair is detectable.
Lemma 3 (see [10] (the Bounded Real Lemma) ). For a given scalar , if there is a matrix that satisfies the matrix inequality then the system (16) is asynchronously stable and it achieves a specified attenuation level , such that .
3. Main Results
For the linear discrete asynchronous multirate sampleddata system (9), the exogenous parts of the system are lifted to improve its performance. To avoid multiplying the dimensions of the control inputs and measured outputs, we do not lift them, thereby avoiding causal constraints of the traditional lifting method.
If we assume that there is an underlying clock with a base period of , then system (9) represents a fast discrete system with a subperiod of for the discrete process described in [9]. The sampler samples the continuous signal for a discrete period and the holder holds the output signal during period . The rates of the sampler and the holder can be set at different flexible rates; that is, . If we select (13) as the new input, is the controller of system (9).
These concepts are described in Figure 3, where
Theorem 4. Consider the sampleddata system (9) with the lifting process shown in Figure 3. For a given attenuation level , if there exist matrices , , , and ( denotes the iteration time), which are suboptimal solutions of the inequalities below, then the closed system of (9) with the controller (13) is asynchronously stable and satisfies the norm constraint :
Proof. After lifting system (9) according to Figure 3, the transform function of a completely discrete system with a period of , which is partly contained within the dotdash line, is
Based on the above, we can derive minimal statespace representations for as follows based on the transfer function theory given in [9, 15]. However, in contrast to the theory, and are not blocks of unit matrices because sampling and holding are not synchronous; that is, the entire signal cannot be sampled or held at the same time:
Next, we calculate the minimal statespace realization of (22) as follows.
(a) Transfer Function for . Note that , which comes directly from the theory in [8], and its corresponding state model is
where
(b) Transfer Function for . As in is not an blocks unit matrix, the theory cannot be applied directly. However, note that
Then,
The corresponding state model can be obtained as
where and have been defined previously and
(c) Transfer Function for . Similarly, because is also not an blocks unit matrix, can be changed to
Thus, the state model can be derived as
where and were defined previously and
(d) Transfer Function for . Like to parts and , the state description of is
where , , and were defined previously, while
Finally, formula (23) is represented using the statespace form as follows:
where the underlining below and indicates that the signals are lifted.
According to Lemma 2, the controller exists and can stabilize the system.
Let , and after combining (13), (15), and (34), we obtain
where
In the light of [5], we set
By substituting in Lemma 3 with , we can obtainThen, according to [26], after multiplying the left and righthand sides of inequality (38) by , we can obtain inequality (20).
According to quadratic stability theory, the system should be stable at each step . Thus, for a discrete linear timeinvariant system, the eigenvalues of the closed system should lie within a unit circle.
From Theorem 4, if matrix inequalities (20) can be solved, there exists a state feedback controller that can asymptotically stabilize closed system (35) of enlarged system (34), which satisfies the performance index. Since system (34) is lifted from the original system (9) and we know that the lifting operator is an isometric isomorphism, that is, the feedback stability and system norm can be preserved, then system (16), which is the closed system of (9), is asynchronously stable and it satisfies the norm constraint .
The synchronous case can be regarded as a special example of the asynchronous case. In this case, the elements of sampling and holder are unit matrices with the assumption that the whole signals of sampling and holding can be obtained at the beginning of every period; then the system (9) turns into a linear discrete synchronous multirate sampleddata system with period as follows:
As and are unit matrices, that is, constant now, , , , and related to and are also constant. Theorem 4 can be applied to system (39) straightforwardly.
Corollary 5. Consider the sampleddata system (39) with the lifting process shown in Figure 3. For a given attenuation level , if there exist matrices , , , and , which are suboptimal solutions of the inequalities below, then the closed system of (39) with the controller (13) is synchronously uniformly stable and it satisfies the norm constraint :
Remark 6. Let us compare the designed observer (12) with that in [21], where the observer (predictor) is The sum of outputs is used to avoid the singularity but may cause more noise. The observer applied in this paper is simpler; that is, just the current output is introduced in the observer. This will result in better tracking performance and smaller state ripple than [21], where the past outputs other than current output are used for feedback.
Remark 7. In Theorem 4 and Corollary 5, due to the existence of inverse matrices and , inequalities (20) and (40) are not LMIs and cannot be directly solved by the LMI tool. However, notice that the matrices and their inverse matrices appear in pairs; we can solve the inequalities by iteration on LMIs according to Cone Complementarity Linearization Algorithm [22].
4. Numerical Example with Simulation
Consider the original plant described by where
If the fast discretization formulae follow , , , , and , then system (2) can be obtained. The initial state is .
If , the inputholding mechanism can be described as
If and , the outputsampling matrices are assumed to be
Therefore, this is a linear discrete asynchronous dualrate sampleddata system [18], which applies a fastinputslowoutput mechanism for .
For comparison, we introduce a previous method [21], where the system is not lifted according to the pure LQG theory, and also let be Gaussian white noise; amplitude is increased to 10, which is much bigger to the system. The results obtained using this method are shown in Figures 4, 5, 6, 7, and 8, where s.
Next, the synthesis method described in Theorem 4 is applied to the same plant. The simulations obtain the transients of the system shown in Figures 9, 10, 11, 12, and 13, where , , and , while (monorate).
A comparison of Figures 4–8 with Figures 9–13 shows clearly that the method performs better than LQG; the controlled states are consistent with the observed states with smaller ripples, especially state , the biggest ripple of which is almost half of the LQG method. And the output of method has better performance of disturbance attenuation than LQG method.
When s and (monorate), applying Theorem 4 again, we can get Figures 14, 15, 16, 17, and 18. Then the lifting number is modified to (multirate); the results shown in Figures 19, 20, 21, 22, and 23 are obtained.
Comparing Figures 14–18 with Figures 19–23, when the lifting amplitude increases from 1 to 2, the controlled states are almost consistent with the observed states with smaller ripples; here the biggest ripple of state is almost half of the monorate. And the output of the multirate case has better performance of disturbance attenuation than the monorate case.
Furthermore, we change the lifting amplitude to , adjust the underlying period of system (2) to and can draw the similar pictures; here we omit them but only list their minimum attenuation levels to compare them with different cases; see Table 1.

The inequality acquired using the theorem is almost linear, except in a few terms. Although these terms are just the inverses of their corresponding terms, they could easily be calculated by LMI tool of Matlab under a constraint. Thus, it is far easier to solve the controller and observer under the same constraint comparing with most of the existing methods, which use highly complex formulae.
5. Conclusion
Motivated by a desire to provide a better solution for the design of a controller for a discretetime asynchronous multirate sampleddata system, we developed synthesis methods in this paper to improve performance and guarantee closedloop transient behavior via observerbased control. Furthermore, an extended lifting technique was adopted to change the original system into a fastinputslowoutput discrete system, thereby improving its performance. Finally, we presented an example that demonstrates the effectiveness of these methods in controlling a discretetime asynchronous multirate sampleddata system, whose performance is thus improved considerably.
In this paper, we focus on the standard control of asynchronous multirate sampleddata system, and the input noise discussed here is Gaussian white for simplicity. As is well known, practical input signals often have finite frequency (FF), so it is reasonable to consider the performance in finite frequency range for sampleddata systems. As in [27–31], where KalmanYakubovicPopov (KYP) lemma is used to establish the equivalence between a frequency domain inequality (FDI) and a linear matrix inequality, it is a valuable work to extend the results in this paper to the FF framework based on KYP lemma. This will be the direction for future development of this paper.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by Zhejiang Provincial Natural Science Foundation of China under Grant no. LR12F03002 and National Natural Science Foundation of China (NSFC) under Grant no. 61074045.
References
 M. F. Sågfors, H. T. Toivonen, and B. Lennartson, “Statespace solution to the period multirate ${H}_{\infty}$ control problem: a lifting approach,” IEEE Transaction on Automatic Control, vol. 37, no. 5, pp. 2345–2350, 2000. View at: Google Scholar
 X. Zhao, Y. Yao, and J. Ma, “Overview of multirate sampleddata control system,” Journal of Harbin Institute of Technology, vol. 9, no. 4, pp. 405–410, 2002. View at: Google Scholar
 Y. Fukada, “Slipangle estimation for vehicle stability control,” Vehicle System Dynamics, vol. 32, no. 4, pp. 375–388, 1999. View at: Google Scholar
 T. Chen and L. Qiu, “${H}_{\infty}$ design of general multirate sampleddata control systems,” Automatica, vol. 30, no. 7, pp. 1139–1152, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Wang, B. Huang, and H. Unbehauen, “Robust ${H}_{\infty}$ observer design of linear state delayed systems with parametric uncertainty: the discretetime case,” Automatica, vol. 35, no. 6, pp. 1161–1167, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in Systems and Control Theory, vol. 15 of SIAM Studies in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1994. View at: Publisher Site  MathSciNet
 S. Patra and A. Lanzon, “A closedloop data based test for robust performance improvement in iterative identification and control redesigns,” Automatica, vol. 48, no. 10, pp. 2710–2716, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Ding, Y. Shi, H. Wang, and F. Ding, “A modified stochastic gradient based parameter estimation algorithm for dualrate sampleddata systems,” Digital Signal Processing, vol. 20, no. 4, pp. 1238–1247, 2010. View at: Publisher Site  Google Scholar
 T. W. Chen and B. Francis, Optimal SampledData Control Systems, Springer, London, UK, 1995.
 L. Shen and M. J. Er, “Multiobjectives design of a multirate output controller,” in Proceedings of the IEEE International Conference on Control Applications (CCA '99), pp. 193–198, August 1999. View at: Google Scholar
 M. Chilali and P. Gahinet, “${H}_{\infty}$ design with pole placement constraints: an LMI approach,” IEEE Transactions on Automatic Control, vol. 41, no. 3, pp. 358–367, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Mo, X. Chen, J. Zhao, J. Qian, and Z. Shao, “A twostage method for identification of dualrate systems with fast input and very slow output,” Industrial and Engineering Chemistry Research, vol. 48, no. 4, pp. 1980–1988, 2009. View at: Publisher Site  Google Scholar
 M. De la Sen and S. AlonsoQuesada, “Model matching via multirate sampling with fast sampled input guaranteeing the stability of the plant zeros: extensions to adaptive control,” IET Control Theory & Applications, vol. 1, no. 1, pp. 210–225, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. F. Sågfors and H. T. Toivonen, “${H}_{\infty}$ and LQG control of asynchronous sampleddata systems,” Automatica, vol. 33, no. 9, pp. 1663–1668, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 S. Lópezlópez, Optimal H_{∞} Design of Causal Multirate Controllers and Filters, Mechanical and Aerospace Engineering, University of California, Irvine, Calif, USA, 2010.
 U. Borison, “Selftuning regulators for a class of multivariable systems,” Automatica, vol. 15, no. 2, pp. 209–215, 1979. View at: Publisher Site  Google Scholar
 V. S. Ritchey and G. F. Franklin, “A stability criterion for asynchronous multirate linear systems,” IEEE Transactions on Automatic Control, vol. 34, no. 5, pp. 529–535, 1989. View at: Publisher Site  Google Scholar  MathSciNet
 B. Shen, Z. Wang, and X. Liu, “Sampleddata synchronization control of dynamical networks with stochastic sampling,” IEEE Transactions on Automatic Control, vol. 57, no. 10, pp. 2644–2650, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Longo, G. Herrmann, and P. Barber, “Robust scheduling of sampleddata networked control systems,” IEEE Transactions on Control Systems Technology, vol. 20, no. 6, pp. 1613–1621, 2011. View at: Publisher Site  Google Scholar
 P. Voulgaris, “Control of asynchronous sampled data systems,” IEEE Transactions on Automatic Control, vol. 39, no. 7, pp. 1451–1455, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Colaneri, R. Scattolini, and N. Schiavoni, “LQG optimal control of multirate sampleddata systems,” IEEE Transactions on Automatic Control, vol. 37, no. 5, pp. 675–682, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. El Ghaoui, F. Oustry, and M. AitRami, “A cone complementarity linearization algorithm for static outputfeedback and related problems,” IEEE Transactions on Automatic Control, vol. 42, no. 8, pp. 1171–1176, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Seuret, “Stability analysis of networked control systems with asynchronous sampling and input delay,” in Proceedings of the American Control Conference (ACC '11), pp. 533–538, July 2011. View at: Google Scholar
 B. Tavassoli, “Stability of nonlinear networked control systems over multiple communication links with asynchronous sampling,” IEEE Transaction on Automatic Control, 2013. View at: Publisher Site  Google Scholar
 L. Yang and J. M. Li, “Sufficient and necessary conditions of controllability and observability of a class of linear switching system,” Systems Engineering and Electronics, vol. 25, no. 5, pp. 588–590, 2003. View at: Google Scholar
 W.J. Mao, “Observerbased energy decoupling of linear timedelay systems,” Journal of Zhejiang University, vol. 37, no. 5, pp. 499–503, 2003. View at: Google Scholar
 T. Iwasaki and S. Hara, “Generalized KYP lemma: unified frequency domain inequalities with design applications,” IEEE Transactions on Automatic Control, vol. 50, no. 1, pp. 41–59, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 H. Gao and X. Li, “${H}_{\infty}$ filtering for discretetime statedelayed systems with finite frequency specifications,” IEEE Transactions on Automatic Control, vol. 56, no. 12, pp. 2935–2941, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 X. Li, H. Gao, and C. Wang, “Generalized KalmanYakubovichPopov lemma for 2D FM LSS model,” IEEE Transactions on Automatic Control, vol. 57, no. 12, pp. 3090–3103, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 X. Li and H. Gao, “Robust finite frequency ${H}_{\infty}$ filtering for uncertain 2D systems: the FM model case,” Automatica, vol. 49, no. 8, pp. 2446–2452, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 X. W. Li and H. J. Gao, “A heuristic approach to static outputfeedback controller synthesis with restricted frequencydomain specifications,” IEEE Transaction on Automatic Control, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Xiaoqiang Sun and Weijie Mao. 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.