IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 52, NO. 6, JUNE 2005
1179
Hierarchical Identification of Lifted State-Space Models for General Dual-Rate Systems Feng Ding and Tongwen Chen
Abstract—This paper is motivated by practical consideration that the input updating and output sampling rates are often limited due to sensor and actuator speed constraints. In particular, for general dual-rate systems with different updating and sampling periods, we derive the lifted state-space models (mapping relations between available dual-rate input-output data), and, by using a hierarchical identification principle, present combined parameter and state estimation algorithms for identifying the canonical lifted models based on the given dual-rate input-output data, taking into account the causality constraints of the lifted systems. Finally, we give an illustrative example to indicate that the proposed algorithm is effective. Index Terms—Hierarchical identification principle, Kalman filtering, least squares, multirate systems, parameter estimation, state-space model, stochastic approximation, system identification.
I. INTRODUCTION
M
ANY conventional algorithms are used to estimate system parameters, see, e.g., [1]–[4], but most of them assume that input-output data are available at every sampling instant. For multirate sampled-data systems, these algorithms cannot be applied directly. For decades, the study of multirate optimal control systems has been very active, including [5]–[8], adaptive control [9]–[11], and so on; application examples include distillation columns [9], polymer reactors [12], fermentation processes [13], and bioreactors [14], [15]. In the process control area, Sheng et al. proposed the generalized predictive control scheme for multirate systems [16]; Li et al. applied dual-rate modeling to Octane quality inferential control and did the performance analysis [17], [18]. In the process identification literature, Lu and Fisher used projection and least-squares algorithms for estimating intersample outputs [19], [20], but their algorithms handle only noise-free dual-rate systems; Li et al. assumed that the system states were available and used them and multirate input-output data to estimate the parameters of lifted state-space models for multirate systems [21]. (However, the system states are usually unavailable in Manuscript received April 15, 2004; revised October 16, 2004. This work was supported by the Natural Sciences and Engineering Research Council of Canada and the National Natural Science Foundation of China. This paper was recommended by Associate Editor C. Hadjicostis. F. Ding is with the Control Science and Engineering Research Center, Southern Yangtze University, Wuxi 214122, China, and also with the University of Alberta, Edmonton, AB T6G 2V4, Canada (e-mail:
[email protected];
[email protected]). T. Chen is with the Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB T6G 2V4, Canada (e-mail:
[email protected]). Digital Object Identifier 10.1109/TCSI.2005.849144
practice.) This lifting conversion from a time-varying system into a time-invariant one is similar to the construction presented in [22]. Recently, subspace state-space identification (4SID for short) methods have been developed for state-space models [23]–[31], the basic idea is to determine the extended observability matrix from the singular value decomposition or RQ factorization of an information matrix consisting of given input-output data, and then to compute the system parameter matrices. But, as the size of the information matrix grows, the difficulty and complexity in computation increase. Though these 4SID methods are popular and effective, the main disadvantage is that they cannot deal with the causality constraints generated by the lifted multirate systems [5]–[7], [16], [21]. For related work, see recent studies on identification of dual-rate systems by using the polynomial transformation technique [32], [33] and by using an auxiliary model identification method [34], [35], and also a survey on various identification schemes for multirate systems [36]. This paper focuses on the modeling and identification problem for general dual-rate systems. The key idea is the so-called hierarchical identification, and is inspired by the hierarchical control based on the decomposition-coordination principle for large-scale systems [37]–[39]. Hierarchical identification uses subsystem decomposition in identification, and is also called bootstrap identification. We tackle the causality problem by using the hierarchical identification principle [40]–[43], and directly identify the parameters of the canonical lifted models of dual-rate systems. The approach here differs from the one in [44] which identified the parameters of the noncanonical lifted models of dual-rate systems. Earlier work on conventional multiple-input–multiple-output (MIMO) system identification exists: some methods estimated the parameters of the difference equation models obtained from state-space models, e.g., [45]–[49]; others estimated the parameters of state-space models by using the two-stage bootstrap identification technique [3], [4]. They all assumed that the system observability indexes were known. In this paper, we extend this two-stage approach to identify the parameters of dual-rate sampled-data systems with causality constraints. The paper is organized as follows. Section II describes the modeling issues related to dual-rate systems. Section III derives the lifted state-space models of dual-rate systems. Section IV discusses the identification problem of the lifted models. Section V presents an illustrative example for the proposed algorithms in this paper. Finally, concluding remarks are given in Section VI.
1057-7122/$20.00 © 2005 IEEE
1180
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 52, NO. 6, JUNE 2005
with the sampling period to get , and assume that state-space model of the form (see, e.g., [50])
has a
(3) Fig. 1. General dual-rate systems (T =
ph; T
=
qh)
.
where the two matrices
and
are as follows:
II. PROBLEM FORMULATION The focus of this paper is a class of multirate systems—the general dual-rate systems with additive disturbance [21]—as is a continuous-time process, the input depicted in Fig. 1, to is produced by a zero-order hold with period , processing a discrete-time signal ; the output of , which is corrupted by the noise , is sampled by with period , yielding a discrete-time signal a sampler . If , then we obtain a conventional discrete-time system, i.e., single-rate system. Without loss of generality, we and and are two coprime assume that integers, for otherwise, we can absorb any common factor of and into is a positive real number called the base period. s, s, then , For example, if and s. For such a dual-rate sampled-data system, the input–output . Since data available are the zero-order hold is used, we have for
For tractability, the lifting technique is used to transform a periodic time-varying system into a time-invariant one with a larger (i.e., the least common multiple of and ), period known as the frame period. However, if we group every input values and every output samples together, respectively, then we obtain a LTI system operating over period prois ; this is the idea of lifting [17], [18], vided that and [21]. to in Fig. 1, we For the dual-rate system from intend to lift the input and output to get a lifted system, which is ; and a single-rate system with underlying frame period and by , where represents hence, we need to lift by the -fold lifting operator [21]. The lifted process is denoted by from to , i.e., . For convenience, we introduce the following notation:
(1)
is a linear time-inThroughout the paper, we assume variant (LTI) continuous-time process with the following statespace representation: (2) is the state vector, the control input where the output vector, a stochastic vector, , and are matrices noise vector with zero mean, and of appropriate dimensions. is LTI, the system from to Although in Fig. 1 is linear periodically time-varying due to different updating and sampling periods [5]–[7], [21]. For such a timevarying dual-rate system, our objective is two-fold. • First, by using the lifting technique, to establish time-invariant models of the system from to ; That is, to find the mapping relationship (lifted state-space models) between the available input and output data. • Second, by means of the hierarchical identification principle, to develop recursive algorithms for combined parameter and state estimation based on the given dual-rate , measurement data taking causality constraints into consideration. III. LIFTED STATE-SPACE MODELS In order to derive the state-space model for the dual-rate system, we discretize via the zero-order hold method
Define the lifted input vector and lifted noise vector as
, lifted output vector
.. .
.. .
.. . Theorem 1: For the dual-rate system in Fig. 1, the coprime, ness of integers and implies that for every and such that there exist integers
DING AND CHEN: HIERARCHICAL IDENTIFICATION OF LIFTED STATE-SPACE MODELS
Then, the lifted state-space model is given by
1181
Thus (4)
where
.. .
or
the equation at the bottom of the page, and the following equations are true:
(5) Similarly, we have
is
is of dimensions . Proof: Replacing , we have
is in (3) with
is and noting that
Hence, it is not difficult to get
or
Using (1), we have
for
.. .
(6)
..
.
..
.
.. . .. . .. .
1182
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 52, NO. 6, JUNE 2005
Combining (5) and (6) leads directly to the results of Theorem 1.
.. .
IV. IDENTIFICATION ALGORITHMS For convenience, let
, and , then (4) can be rewritten as .. .
(7) Noting that is a (block) lower triangular matrix, i.e., the are zero; this represents the upper triangular blocks in so-call causality constraint. Since the models in (7) are derived from the continuous-time models, the causality constraint is automatically satisfied. However, we need to deal with the causality constraint for the lifted system identification. Identifiability of a system depends on its controllability and observability. Therefore, it is important if the lifted state-space models in (7) are controllable and observable. What conditions are required for the lifted models in (7) to be controllable and observable? To answer this question, we assume that the stateis conspace model in (3) is a minimal realization, i.e., trollable and is observable. Note that this is guaranteed is controllable as long as the continuous-time process and is observable and if the sampling period is nonpathological [50]. In the same way, for the nonpathological samis controllable and observpling period able. Controllability and observability of (7) can be achieved under a mild condition. Lemma 1: points 1) If for every eigenvalue of , none of the is an eigenimplies that value of , then controllability of of each , which in turn implies controllability of the model in (7). of , none of the 2) If for every eigenvalue points is an eigenvalue of , controllability of implies that of each , which in turn implies controllability of the model in (7); observability of implies that of each , which in turn implies observability of the model in (7). The proof is omitted here, but can be obtained in a similar way as that given in [51]. The following is to apply the hierarchical identification principle to study the identification problem of the lifted model. Let
.. . Any controllable and observable state-space model is equivalent to an observability canonical form that has the least number of parameters. Assume that (3) is observable and the conditions of Lemma 1 hold. Then, there exists a nonsingular matrix such that the transform makes the model in (7) become the following observability canonical form: (8) where
.. .
.. .
.. .
.. .
..
.. .
.. .
.. .
.. . .. .
.. .
.. .
.. . .
DING AND CHEN: HIERARCHICAL IDENTIFICATION OF LIFTED STATE-SPACE MODELS
The transform matrix is given by ( the matrix )
represents the th row of
1183
Then, (9) is decomposed into
susbsystems:
or
.. . .. .
.. .
.. .
..
.. . Here, the observability structure indexes are assumed to . (If not, they can be known and satisfy be determined by using a similar approach as in [46].) Our goal is to estimate the unknown parameters in the system and unknown state parameter matrices of the canonical model in (8) by utilizing the lifted vector . input-output data From the structure of , and , a question arises: whether the number of parameters in the lifted system is less or or not if the original continuous-time system the discrete-time system is assumed to be canonical? The answer is negative. The assumption of or being canonical would not help, because there is no guarantee that the lifted system would be anonical; moreis canonical does not imply is canonical, and over, that is an visa versa. For example, if one assumes that observability canonical form, i.e., is a companion matrix and is a matrix consisting of unit row vectors, the lifted is not a companion matrix; hence, the transformation is also required. In order to estimate the unknown parameters in using the idea of decomposition identification [3], [4], [45], [47]–[49] as in the bootstrap algorithm, we divide (8) into the following two systems:
.. . .
.. .
.. .
.. .
(11)
.. .
.. .
(12) th row of Define the parameter matrix
and information vector
as
From (11), we have
(13)
(9) (10) According to the hierarchical identification principle, we first estimate the parameter matrices and state vector of the system in (9) from the available multirate inputoutput data, and then the parameter matrices of the system in (10) from the obtained state estimate and lifted input-output data. Let
.. .
.. .
Let be a forward shift operator: the th equation of (13) by summation is
, multiplying , and their
(14) From (12) and (14), we have (15) Here, we can see that this identification expression contains both the unknown state vector in and the unknown parameter vector . In order to identify them, according to the hierarchical identification principle, when recursive estimating the parameter vector , the unknown state vector in is replaced by its corresponding estimate , and by ; in the same way, when estimating the state vector , the unknown parameter vectors are also replaced by their estimates at time . This lead to the following hierarchical
1184
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 52, NO. 6, JUNE 2005
identification algorithms consisting of the parameter estimation and state estimation one. Parameter Estimation Algorithm: (30) (16) (17)
(18) are the estimate of where given later. of State Estimation Algorithm:
at time
is the estimate
where the estimates of are formed by using the entries of the obtained or . Here, the parameter estimator is derived by using the state estimates based on the least squares principle, and the state estimator is derived by using the obtained parameter estimates based on the Kalman filtering principle. But the state estimation algorithm in (19)–(21) based on the Kalman filtering involves heavy computational efforts; for computational efficiency, a stochastic approximation algorithm may be used for state estimation
(19) (20)
or
(21) .. .
.. .
.. .
.. .
..
(22)
.. .
(31) where
For example, we may take (23)
.
is the convergence factor satisfying
or
(32)
The following is to discuss the identification problem of the of the system in (10). Let parameter matrices .. .
.. .
(24)
So the system in (10) may be equivalently written as
.. .
(25)
From this, all the entries in the matrix the simple least squares algorithm
.. .
(26)
.. .
(27)
can be estimated by
or (28)
and (29)
for large However, since the upper triangular blocks in the zero entries of this upper triangular block in
are zero, i.e., need not be
DING AND CHEN: HIERARCHICAL IDENTIFICATION OF LIFTED STATE-SPACE MODELS
identified, in order to tackle the so-called causality constraint for the lifted model, we decompose the system in (10) into subsystems
1185
V. EXAMPLE For a system depicted in Fig. 1, take the process model to be
and s, , hence s, s, and s. The identification procedure is summarized as follows: First, we generate a random signal sequence with zero mean and unit ; Second, we lift by and variance as the input signal to form the lifted signals lift by
where
and then the corresponding lifted state-space model is given by
.. . th block row of
Then, for large , the estimate of
is obtained from
is taken as a white noise sequence with zero mean and variance . Next, take , based on the lifted signals, we apply the proposed algorithm in (16)–(18) and (31)–(34) to estimate the parameters of the lifted model of the canonical form
.. .
The estimate
may be also recursively computed by
(33)
(34) The combined parameter and state estimation algorithms in (16)–(30) and (33)–(34) form the hierarchical identification algorithms of the lifted state-space model for dual-rate systems based on the least squares and Kalman filtering principles; and (16)–(18) and (31)–(34) form the hierarchical identification algorithms based on the least squares and stochastic approximation principles. To obtain consistent parameter estimation, the input signals are required to be persistently exciting from the formulation of the information vectors above. To initialize the above algorithms, we take , and with normally a large positive number (e.g., ), and and , some small real vector, e.g., , and so on, with being an -dimensional vector whose elements are all 1.
The true parameters and their estimates with different data lengths are shown in Table I, and the parameter versus is shown in Fig. 2, where estimation error represents the true parameter vector, the estimate of . The estimated lifted model is given by
1186
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 52, NO. 6, JUNE 2005
TABLE I PARAMETERS AND THEIR ESTIMATES ( = 0:10 )
Fig. 3 compares step responses of the actual system and the estimated lifed model. From Table I and Figs. 2 and 3, we can see that the parameter estimation error is becoming smaller (in general) as increases, and the step response of the estimated . This inmodel is very close to that of the actual system dicates that the estimated lifted model can capture the process dynamics very well, and can achieve satisfactory results. VI. CONCLUSION
Fig. 2. Parameter estimation error versus k ( = 0:10 ).
This paper presents hierarchical state-space model identification algorithms for general dual-rate stochastic systems. Although the algorithms are developed for output error models, namely, output measurement containing additive white noise disturbances, the method adopted can be easily extended to other cases, for example, dual-rate stochastic systems with colored noises and general multirate multivariable systems with each of input and output channels having different sampling periods, and with colored noises. Under what conditions the parameter estimation by the proposed algorithms is convergent still requires further investigation. REFERENCES
Fig. 3. Step responses of the actual system and the estimated model.
The estimated lifted model is very close to the true lifted canonical model
[1] T. Söderström and P. Stoica, System Identification. Englewood-Cliffs, NJ: Prentice-Hall, 1988. [2] L. Ljung, System Identification: Theory for the User, 2nd ed. Englewood-Cliffs, NJ: Prentice-Hall, 1999. [3] H. El-Sherief and N. Sinha, “Bootstrap estimation of parameters and states of linear multivariable systems,” IEEE Trans. Autom. Contr., vol. AC-24, no. 2, pp. 340–343, Feb. 1979. [4] R. M. Prasad and A. K. Sinha, “On bootstrap identification using stochastic approximation,” IEEE Trans. Autom. Contr., vol. 22, no. 4, pp. 671–672, Apr. 1977. [5] T. Chen and L. Qiu, “ design of general multirate sampled-data control systems,” Automatica, vol. 30, no. 7, pp. 1139–1152, 1994. [6] L. Qiu and T. Chen, “ -optimal design of multirate sampled-data systems,” IEEE Trans. Autom. Contr., vol. 39, no. 12, pp. 2506–2511, Dec. 1994. [7] , “Multirate sampled-data systems: All suboptimal controllers and the minimum entropy controllers,” IEEE Trans. Autom. Contr., vol. 44, no. 3, pp. 537–550, Mar. 1999. [8] M. F. Sagförs, H. T. Toivonen, and B. Lennartson, “State-space solution to the periodic multirate H control problem: A lifting approach,” IEEE Trans. Autom. Contr., vol. 45, no. 12, pp. 2345–2350, Dec. 2002. [9] R. Scattolini, “Self-tuning control of systems with infrequent and delayed output sampling,” in Proc. Inst. Elect Eng., Contr. Theory Appl., vol. 135, 1988, pp. 213–221. [10] C. Zhang, R. H. Middleton, and R. J. Evans, “An algorithm for multirate sampling adaptive control,” IEEE Trans. Autom. Contr., vol. 34, no. 7, pp. 792–795, Jul. 1989.
H H
H
DING AND CHEN: HIERARCHICAL IDENTIFICATION OF LIFTED STATE-SPACE MODELS
1187
[41]
, “Iterative least squares solutions of coupled Sylvester matrix equations,” Syst. Contr. Lett., vol. 54, no. 2, pp. 95–107, 2005. F. Ding, J. B. Yang, and Y. M. Xu, “Convergence of hierarchical stochastic gradient identification for transfer function matrix model,” Contr. Theory Appl., vol. 18, no. 6, pp. 949–953, 2001. , “Hierarchical identification and its convergence for the transfer function matrix,” in Proc. IFAC Symp. System Identification, Santa Barbara, CA, Jun. 21–23, 2000. F. Ding, T. Chen, and D. Xiao, “State-space modeling and identification for general dual-rate stochastic systems,” Acta Automatica Sinica, vol. 30, no. 5, pp. 652–663, 2004. R. Guidorzi, “Canonical structures in the identification of multivariable systems,” Automatica, vol. 11, no. 4, pp. 361–374, 1975. R. Guidorzi, M. Losito, and T. Muratori, “The range error test in the structural identification of linear multivariable systems,” IEEE Trans. Autom. Contr., vol. AC-27, no. 5, pp. 1044–1054, Oct. 1982. A. Gauthier and I. D. Landau, “On the recursive identification of multiinput multioutput systems,” Automatica, vol. 14, no. 6, pp. 609–614, 1978. , “Parametric identification of a state-space model of multivariable systems using the extended least-squares method,” IEEE Trans. Syst. Man Cybern., vol. 11, no. 3, pp. 223–227, Mar. 1981. N. K. Sinha and Y. H. Kwong, “Recursive estimation of the parameters of linear multivariable systems,” Automatica, vol. 15, no. 4, pp. 471–475, 1979. T. Chen and B. Francis, Optimal Sampled-Data Control Systems. London, U.K.: Springer-Verlag, 1995. B. A. Francis and T. T. Georgiou, “Stability theory for linear time-invariant plants with periodic digital controllers,” IEEE Trans. Autom. Contr., vol. 33, no. 9, pp. 820–832, Sep. 1988.
[11] P. Albertos, J. Salt, and J. Tormero, “Dual-rate adaptive control,” Automatica, vol. 32, no. 7, pp. 1027–1030, 1996. [12] M. Ohshima, I. Hashimoto, M. Takeda, T. Yoneyama, and F. Goto, “Multirate multivariable model predictive control and its application to a semi-commercial polymerization reactor,” in Proc. ACC, vol. 2, 1992, pp. 1576–1581. [13] R. D. Gudi, S. L. Shah, and M. R. Gray, “Multirate state and parameter estimation in an antibiotic fermentation with delayed measurements,” Biotechnol. Bioeng., vol. 44, no. 11, pp. 1271–1278, 1994. , “Adaptive multirate state and parameter estimation strategies with [14] application to a bioreactor,” AIChE J., vol. 41, no. 11, pp. 2451–2464, 1995. [15] S. Tatiraju, M. Soroush, and R. Mutharasan, “Multirate nonlinear state and parameter estimation in a bioreactor,” Biotechnol. Bioeng., vol. 63, no. 1, pp. 22–32, 1999. [16] J. Sheng, T. Chen, and S. L. Shah, “Generalized redictive control for nonuniformly sampled systems,” J. Process Contr., vol. 12, no. 8, pp. 875–885, 2002. [17] D. Li, S. L. Shah, and T. Chen, “Analysis of dual-rate inferential control systems,” Automatica, vol. 38, no. 6, pp. 1053–1059, 2002. [18] D. Li, S. L. Shah, T. Chen, and K. Z. Qi, “Application of dual-rate modeling to CCR octane quality inferential control,” IEEE Trans. Contr. Syst. Technol., vol. 11, no. 1, pp. 43–51, Jan. 2003. [19] W. P. Lu and D. G. Fisher, “Output estimation with multirate sampling,” Int. J. Contr., vol. 48, no. 1, pp. 149–160, 1988. , “Least-squares output estimation with multirate sampling,” IEEE [20] Trans. Autom. Contr., vol. 34, no. 6, pp. 669–672, Jun. 1989. [21] D. Li, S. L. Shah, and T. Chen, “Identification of fast-rate models from multirate data,” Int. J. Contr., vol. 74, no. 7, pp. 680–689, 2001. [22] D. Hristu and K. Morgansen, “Limited communication control,” Syst. Contr. Lett., vol. 37, no. 4, pp. 193–205, 1999. [23] M. Verhaegen and P. Dewilde, “Subspace model identification part 1: The output-error state space model identification class of algorithms,” Int. J. Contr., vol. 56, no. 5, pp. 1187–1210, 1992. [24] , “Subspace model identification part 2: Analysis of the elementary output-error state space model identification algorithm,” Int. J. Contr., vol. 56, no. 5, pp. 1211–1241, 1992. [25] M. Verhaegen, “Subspace model identification part 3: Analysis of the ordinary output-error state space model identification algorithm,” Int. J. Contr., vol. 58, no. 3, pp. 555–586, 1993. [26] , “Identification of the deterministic part of MIMO state space models given in innovations form from input-output data,” Automatica, vol. 30, no. 1, pp. 61–71, 1994. [27] P. Van Overshee and B. De Moor, “N4SID: Subspace algorithm for the identification of combined determine-stochastic systems,” Automatica, vol. 30, no. 1, pp. 75–93, 1994. [28] C. T. Chou and M. Verhaegen, “Subspace algorithms for the identification of multivariable dynamic error-in-variable models,” Automatica, vol. 33, no. 10, pp. 1857–1869, 1997. [29] T. Gustafsson, “Instrumental variable subspace tracking using projection approximation,” IEEE Trans. Signal Process., vol. 46, no. 3, pp. 669–681, Mar. 1998. [30] M. Lovera, T. Gustafsson, and M. Verhaegen, “Recursive subspace identification of linear and nonlinear Wiener state-space models,” Automatica, vol. 36, no. 11, pp. 1651–1658, 2000. [31] H. Oku and H. Kimura, “Recursive 4 SID algorithms using gradient type subspace tracking,” Automatica, vol. 38, no. 6, pp. 1035–1043, 2002. [32] F. Ding and T. Chen, “Parameter identification and intersample output estimation of a class of dual-rate systems,” in Proc. 42nd IEEE Conf. Decision and Control (CDC), Maui, HI, Dec. 9–12, 2003, pp. 5555–5560. [33] , “Parameter estimation for dual-rate systems with finite measurement data,” Dynam. Contr. Discr. Imp. Syst., vol. 11, no. 1, pp. 101–121, 2004. [34] , “Combined parameter and output estimation of dual-rate systems using an auxiliary model,” Automatica, vol. 40, no. 10, pp. 1739–1748, 2004. [35] , “Identification of dual-rate systems based on finite impulse response models,” Int. J. Adapt. Contr. Signal Process., vol. 18, no. 7, pp. 589–598, 2004. [36] , “Modeling and identification for multirate systems,” Acta Automatica Sinica, vol. 31, no. 1, pp. 105–122, 2005. [37] M. G. Singh, Dynamical Hierarchical Control. New York: North-Holland , 1980. [38] H. Tamura and T. Yoshikawa, Large-Scale Systems Control and Decision Making. New York: Marcel Dekker, 1990. [39] M. Drouin, H. Abou-Kandil, and M. Mariton, Control of Complex Systems, M. G. Singh, Ed. New York: Plenum, 1991. [40] F. Ding and T. Chen, “Hierarchical gradient-based identification of multivariable discrete-time systems,” Automatica, vol. 41, no. 2, pp. 315–325, 2005.
[42] [43] [44] [45] [46] [47] [48] [49] [50] [51]
Feng Ding was born in Guangshui, China. He received the B.Sc. degree in electrical engineering from the Hubei University of Technology, Wuhan, China, in 1984, and the M.Sc. and Ph.D. degrees in automatic control both from Tsinghua University, Beijing, China, in 1991 and 1994, respectively. From 1984 to 1988, he was an Electrical Engineer at the Hubei Pharmaceutical Factory, Hubei, China. He was with Department of Automation, Tsinghua University. He is now a Professor in the Control Science and Engineering Research Center, Southern Yangtze University, Wuxi, China, and has been a Research Associate at the University of Alberta, Edmonton, AB, Canada, since 2002. His current research interests include model identification and adaptive control. He co-authored the book Adaptive Control Systems (Beijing, China: Tsinghua University Press, 2002), and published over 80 papers on modeling and identification as the first author.
Tongwen Chen received the B.Sc. degree from Tsinghua University, Beijing, China, in 1984, and the M.A.Sc. and Ph.D. degrees from the University of Toronto in 1988 and 1991, respectively, all in electrical engineering. From 1991 to 1997, he was an Assistant/Associate Professor in the Department of Electrical and Computer Engineering, University of Calgary, Calgary, AB, Canada. Since 1997, he has been with the Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB, Canada, and is presently a Professor of Electrical Engineering. He held visiting positions at the Hong Kong University of Science and Technology, Tsinghua University, and Kumamoto University. His current research interests include process control, multirate systems, robust control, network based control, digital signal processing, and their applications to industrial problems. He co-authored with B.A. Francis the book Optimal Sampled-Data Control Systems (New York: Springer, 1995). Dr. Chen received a University of Alberta McCalla Professorship for 2000/2001, and a Fellowship from the Japan Society for the Promotion of Science for 2004. He was an Associate Editor of the IEEE TRANSACTIONS ON AUTOMATIC CONTROL during 1998–2000. Currently he is an Associate Editor of Automatica, Systems and Control Letters, and DCDIS Series B. He is a registered Professional Engineer in Alberta, Canada.