Robust Nyquist array analysis based on uncertainty descriptions from system identication
|
|
- Geraldine Townsend
- 6 years ago
- Views:
Transcription
1 Automatica 38 (2002) Brief Paper Robust Nyquist array analysis based on uncertainty descriptions from system identication Dan Chen, Dale E. Seborg Department of Chemical Engineering, University of California, Santa Barbara, CA 93106, USA Received 25 April 2000; revised 18 July 2001; received in nal form 21 August 2001 Abstract A robust Nyquist array analysis for MIMO systems is proposed based on uncertainty descriptions obtained from system identication. Two types of statistical-based uncertainty error bounds for the frequency response are obtained: element bounds and column bounds. Gershgorin s theorem and the concepts of diagonal dominance and Gershgorin bands are extended to include model uncertainty. Robust stability theorems are developed based on ellipsoidal uncertainty descriptions obtained from system identication. An example is given to illustrate the robust Nyquist array analysis.? 2002 Elsevier Science Ltd. All rights reserved. Keywords: Statistical model uncertainty; Gershgorin theorem; Robust Nyquist array analysis 1. Introduction Robust control based on frequency domain analysis such as H methods and -synthesis has been an active research eld during the past 20 years (Skogestad & Postlethwaite, 1996; Zhou, Doyle, & Glover, 1996). These techniques provide very powerful, complex approaches, but the resulting control systems usually have a complicated structure and are high order. Thus, these design techniques are not directly applicable to the decentralized PI=PID control systems that are widely used in the process industries. Furthermore, a typical starting point for contemporary robust control algorithms is that hard bounds are used to characterize model uncertainties. This formulation disregards the statistical model uncertainty descriptions that are readily available from system identication (Goodwin, 1999). Eorts to obtain a better match between robust control and system identication have received attention recently. For example, Braatz and Crisalle (1998) developed a robustness analysis for SISO systems with ellipsoidal parametric uncertainty This paper was not presented at any IFAC meeting. This paper was recommended for publication in revised form by Associate Editor Tor Arne Johansen under the direction of Editor Sigurd Skogestad. Corresponding author. Tel.: ; fax: address: seborg@engineering.ucsb.edu (D.E. Seborg). that is naturally obtained from system identication. A robust stability analysis directly based on the condence ellipses of system frequency response from identication has been developed by Cooley and Lee (1998). Some general conditions for the robustness of MIMO systems with structured uncertainty such as ellipsoidal uncertainty have been developed in the framework of -synthesis (Khatri & Parrilo, 1998; Chellaboina, Haddad, & Bernstein, 1998), but just as for other -synthesis methods these robustness analyses cannot be applied for xed structure control system design. In the present paper, the research objective is to develop a robust control technique based on statistical uncertainty information obtained from system identication, which can be applied not only as a robust stability analysis method for MIMO systems but also as a simple robust design approach for decentralized control systems. Frequency domain techniques, such as Nyquist array analysis (Rosenbrock, 1974) and quantitative feedback theory (QFT) (Horowitz, 1982), provide powerful tools for the analysis and design of decentralized control systems. QFT is a generalized loop-shaping method based on Nichols plots that addresses model uncertainty for single loop and decentralized control systems (Horowitz, 1992); thus, it is primarily a robust design methodology. Modications of QFT have been developed that address robustness to parametric and unstructured uncertainties (Braatz, 1994; Jayasuriya & Zhao, 1994; Chait, Chen, & Hollot, /02/$ - see front matter? 2002 Elsevier Science Ltd. All rights reserved. PII: S (01)
2 468 D. Chen, D.E. Seborg / Automatica 38 (2002) ; Lee, Chait, & Steinbuch, 2000). A disadvantage of the QFT design procedure is that it is rather complicated for MIMO systems. Design approaches based on combinations of the inverse Nyquist array (INA) analysis and QFT have been proposed (Nwokah, Nordgren, & Grewal, 1995; Chen, Wang, & Wang, 1991), but do not consider ellipsoidal uncertainty descriptions. The Nyquist array method for MIMO systems is much more straightforward than the QFT design procedure and is easier to implement, especially for xed structure controllers such as PID control (Maciejowski, 1989). Nyquist array analysis provides a useful extension of the SISO Nyquist method to MIMO systems. Both Gershgorin bands, which address loop interactions, and Nyquist loci have been used to design controllers. Unlike the QFT method, model uncertainty is not considered directly in the original Nyquist array method. Based on a large amount of numerical calculation, Gao and Zhang (1993) derived empirical robust stability conditions for the inverse Nyquist array method and 2 2 and 3 3 systems. More general and theoretical robust Nyquist array methods have been developed to accommodate system uncertainty that can be described by deterministic error bounds on individual transfer functions (Arkun, Manousiouthakis, & Putz, 1984), or be expressed as parametric uncertainty in the rational polynomial functions that describe transfer functions (De Santis & Vicino, 1996; Kontogiannis & Munro, 1997). Since these methods are restricted to transfer function models with rational forms, systems with time delays cannot be analyzed without approximating each time delay as a rational function (e.g., a Pade approximation). Furthermore, these approaches are not general enough to utilize the statistical uncertainty information, especially correlation information that is available from system identication. Thus, unnecessary conservativeness can be introduced by these approaches. In the present paper, Nyquist array analysis is extended to more general system descriptions and to more general model uncertainty descriptions, based on statistical uncertainty bounds obtained from system identication. Two types of uncertainty bounds are obtained for the frequency response matrix: bounds for individual elements and column bounds for each column which are structured error bounds and include correlation information between transfer functions. A new robust Gershgorin theorem is derived for these error bounds descriptions. The new results provide a less conservative robust stability analysis and a simple robust design approach for decentralized control systems. 2. Uncertainty description n y Consider a MIMO linear system with n u inputs and outputs. A discrete-time dynamic model can be written as y(t)=g(q 1 ; )u(t)+h(q 1 ; )e(t); (1) where y(t) R ny, e(t) R ny, u(t) R nu and =[ 1 ;:::; p ] T ; (2) G(q 1 ; )= G k ()q k ; (3) H(q 1 ; )=I + H k ()q k : (4) Vector denotes the model parameters to be identied. In Eqs. (1) (4), q 1 denotes the backward shift operator and e(t) is a sequence of independent random vectors with zero mean values and covariance matrices e. The subsequent analysis is not restricted to the model form given in Eqs. (1) (4). In fact, if the system can be described by any linear model, the subsequent analysis is still valid. Assume that ˆ =[ˆ 1 ;:::; ˆ p ] T is an asymptotically normal and unbiased estimate of. The estimated covariance matrix of the parameter estimation error is assumed to be obtained during system identication. Therefore, an asymptotic limit exists, ( ˆ ) T 1 ( ˆ ) 2 p; (5) where p is the number of model parameters and 2 p is the chi-square distribution with p degrees of freedom. This result enables the construction of condence bounds P[( ˆ ) T 1 ( ˆ ) 6 2 ;p]=1 ; (6) where 1 is the probability that the true parameters lie within the specied ellipsoidal condence region about ˆ and 2 ;p is the upper 100 percentage point of the 2 probability distribution with p degrees of freedom. Let g kl (e j! ; ) be the frequency response from input l to output k. The corresponding estimated frequency response is ĝ kl (e j! ), g kl (e j! ; ˆ). Express the frequency response in terms of real and imaginary parts: g kl (e j! ; )=a kl (!; )+jb kl (!; ); ĝ kl (e j! )=â kl (!)+jˆb kl (!); (7) where â kl (!) and ˆb kl (!) denote the corresponding estimates. The gradient of a kl (!; ) with respect to the parameters is represented by a T a kl(!; ) (8)
3 D. Chen, D.E. Seborg / Automatica 38 (2002) which is a 1 p row vector. Denote the gradient of b kl (!; ) as b b kl(!; ): T (9) Dene akl; (!) g kl (!), b kl; (!) ; (10) ˆ a2 p matrix which is evaluated at = ˆ. Therefore, the estimation error for a kl (!; ) and b kl (!; ) can be approximated by a truncated Taylor series expansion, âkl (!) a kl (!; ) g kl (!), ˆb kl (!) b kl (!; ) g kl (!)( ˆ ): (11) For the remainder of the paper, we use the simpler notation, a kl (!)=a kl (!; ) and b kl (!)=b kl (!; ). The subsequent discussion utilizes the following lemma. Lemma 1 (Wahlberg & Ljung; 1992). Let x R n ; R n n be positive denite; and assume that x T 1 x 6 1: (12) Consider y = Ax; where y R p ;p6 n; A R p n and A has full row rank. Then y T (AA T ) 1 y 6 1: (13) 2.1. Element bounds for the frequency response Condence bounds for each frequency response, g kl, can be derived using Lemma 1 with Eqs. (6) and (11) g kl (!) T [ g kl (!) g kl (!) T ] 1 g kl (!) 6 2 ;p: (14) Let kl (!)= g kl (!) g kl (!) T, which represents the 2 2 covariance matrix for â kl (!) and ˆb kl (!) (Ljung, 1985). The singular value decomposition (SVD) for the symmetric kl (!) matrix is 1;kl (!) 0 kl (!)=U kl (!) U kl (!) T ; (15) 0 2;kl (!) where 1;kl (!) 2;kl (!) and U kl (!) isa2 2 unitary matrix. The following theorem and proof is a formal statement of a result in Cooley and Lee (1998) that was stated without proof. Theorem 1 (Element bounds). An error bound for each frequency response; g kl (e j! ; ); is given by g kl (e j! ; ) ĝ kl (e j! ) = g kl (!) 2 6 ;p 2 1;kl (!): (16) Proof. Substituting Eq. (15) into Eq. (14) gives, 1 0 g kl (!) T U kl (!) 1;kl (!) 1 0 2;kl (!) U kl (!) T g kl (!) 6 2 ;p: (17) Let y =[y 1 ;y 2 ] T = U kl (!) T g kl (!), therefore the above equation can be written as y 2 1 1;kl (!) + y2 2 2;kl (!) 6 2 ;p: (18) This implies that y 2 1 1;kl (!) + y2 2 1;kl (!) 6 2 ;p (19) because 1;kl (!) 2;kl (!). Hence, y ;p 1;kl (!): (20) Because U kl (!) is a unitary matrix, y 2 2 = g kl(!) 2 2, which implies that inequality (16) holds. Theorem 1 provides a circular bound for each frequency response element, g kl (e j!; ; ). This theorem can also be proven by using the Rayleigh Ritz inequality (Rugh, 1996). The uncertainty bounds given by Theorem 1 include a degree of conservatism that is introduced by Lemma 1 and Eq. (19) Column bounds for the frequency response The error bounds in Eq. (16) provide independent bounds on each element of the frequency response matrix. However, they ignore the covariance information between dierent elements. In multivariable design problems, the interaction between elements is an important factor. In Nyquist array analysis, the radius of each Gershgorin circle depends on the entire column or row of the transfer function matrix. Therefore, ignoring the uncertainty interaction between dierent elements can produce very conservative stability results. In this section, a new formulation of the error bound for each column of the frequency response matrix is derived from statistical information provided by system identication.
4 470 D. Chen, D.E. Seborg / Automatica 38 (2002) Dene g l (!)=col([ g 1l (!);:::; g nyl(!)]); (21) where col( ) is an operator that stacks the columns on top of each other and yields a column vector. Thus, g l (!), col([g 1l (!);:::;g nyl(!)]) From Lemma 1, g l (!)( ˆ ): (22) g l (!) T [ g l (!) g l (!) T ] 1 g l (!) 6 2 ;p: (23) Let l (!)= g l (!) g l (!) T. Therefore, l (!) isa 2n y 2n y symmetric matrix representing the covariance matrix of [â 1l (!); ˆb 1l (!);:::;â nyl(!); ˆb nyl(!)] T (Ljung, 1985). The SVD for this matrix is l (!)=U l (!) diag[ 1;l (!);:::; 2ny;l(!)]U l (!) T ; (24) where 1;l (!) 2;l (!) 2ny;l(!), U l (!) is a 2n y 2n y unitary matrix, and diag[ ] denotes a diagonal matrix with elements 1;l (!);:::; 2ny;l(!). Therefore, a new bound for the estimation error of each column of the system frequency response matrix can be derived based on covariance information and the largest singular value. Theorem 2 (Column bounds). The bound for g l (!) 2 ; the error for the lth column of frequency response matrix G(e j! ; ); is given by g l (!) 2 = 6 ( ny ) 1=2 g kl (e j! ; ) ĝ kl (e j! ) 2 2 ;p 1;l (!): (25) Proof. Similar to Theorem 1. In analogy with Theorem 1, the uncertainty bounds given by Theorem 2 may also be conservative due to Lemma 1 and Eq. (19). Eq. (25) gives the error bound for a column of the estimated frequency response matrix while Eq. (16) gives the error bound for each element. Thus, it is inappropriate to compare these two bounds directly. But if the error bound for an entire column of the estimated frequency response matrix is obtained by adding the error bounds in Eq. (16) for each element in this column, this error bound will be larger than the error bound obtained directly from Eq. (25), unless this column has only one nonzero element. The element bound (16) is more conservative because the correlation between different elements in the column is ignored while the covariance information of the entire column is included in the error bound in Eq. (25). 3. Robust Gershgorin bands The key to Nyquist array analysis is Gershgorin s theorem. Next, a new robustness theorem is developed for uncertain complex matrices based on Gershgorin s theorem. Theorem 3 (Robust Gershgorin s theorem). For an n n uncertain complex matrix Z =[z kl ] n n =[a kl + jb kl ] n n with the nominal value Ẑ =[ẑ kl ] n n =[â kl + j ˆb kl ] n n ; the eigenvalues of Z lie in the union of the n circles (Z) ẑ ll 6 R l ; ; 2;:::;n; (26) where R l = ˆR l +R l and ˆR l = min ẑ kl ; 1=2 n 1 ẑ kl 2 ; (27) R l = min z kl ẑ kl ; ( ) 1=2 n z kl ẑ kl 2 : (28) The eigenvalues of Z also lie in the union of the circles; each centered on ẑ kk with radius R k = ˆR k +R k ; where ˆR k = min ẑ kl ; 1=2 n 1 ẑ kl 2 ;l k ;l k ; (29) R k = min z kl ẑ kl ; ( ) 1=2 n z kl ẑ kl 2 : (30) Proof. See the appendix. Theorem 3 is a robust version of Gershgorin s theorem for uncertain matrices. According to this theorem, the eigenvalues of an uncertain matrix must be located inside the union of circles dened by the estimated values and the uncertainty bounds. Thus, for an MIMO control problem, the eigenvalues of the frequency response matrix lie within the union of circles centered on the estimated frequency response. Therefore, Theorem 3 can be used for a Nyquist array analysis of uncertain models. Denition (Diagonal dominance). For an n n uncertain complex matrix Z which has a nominal value Ẑ, if ẑ ll R l = ˆR l +R l for l =1; 2;:::;n (31) then, Z is said to have robust column diagonal dominance for ˆR l and R l in Eqs. (27) and (28). If ˆR k and R k are
5 D. Chen, D.E. Seborg / Automatica 38 (2002) dened as in Eqs. (29) and (20), Z is said to have robust row diagonal dominance. Denition (Robust Gershgorin bands). Consider an n n transfer function matrix Z(s) with the nominal value Ẑ(s). At each ẑ ll (j!), superimpose a circle of radius ˆR l (!) +R l (!) dened as Eqs. (27), (28) or (29), (30) (making the same choice for all the diagonal elements at each frequency). The bands obtained in this way are called robust Gershgorin bands; each is composed of robust Gershgorin circles. 4. Robust Nyquist array analysis Theorem 3 and robust diagonal dominance provide the basis for a robust Nyquist stability criterion for uncertain process models. First, a preliminary result is presented in Lemma 2. Lemma 2. Suppose that Z(s) is an n n uncertain matrix with the nominal value Ẑ(s) and that Z(s) is robust diagonally dominant for l =1; 2;:::;nand for all s on the Nyquist contour. If the lth robust Gershgorin band of Z(s) encircles the origin N l times counterclockwise and det[z(s)] encircles the origin N times counterclockwise; then N = N l : (32) Proof. The robust diagonal dominance of Z(s) ensures that each of its robust Gershgorin bands does not cover the origin. According to Theorem 3, the union of these bands contains the characteristic loci of Z(s) which are denoted by l (s)(l =1;:::;n). Therefore by the principle of the argument (Maciejowski, 1989), 2N =arg det[z(s)] = arg l (s) = 2N l ; (33) where arg denotes the change in the argument as s traverses the Nyquist contour. Hence, Eq. (32) holds. Based on Theorem 3 and Lemma 2, a new robust stability theorem can be obtained for the closed-loop system consisting of a square transfer function G(s) connected with a diagonal controller gain matrix K = diag{k 1 ;:::;k n }. Each k i is a real, nonzero constant in a negative-feedback loop. Theorem 4 (Robust Nyquist stability). Suppose that G(s) is an n n system and has an estimated frequency response matrix Ĝ(s)=[ĝ kl (s)] n n ; that K = diag[k 1 ;:::;k n ] and that the matrix; K 1 + G(s); has robust column diagonal dominance on the Nyquist contour; i.e.; 1 +ĝ k ll (s) l R l(s)= ˆR l (s)+r l (s); (34) where ˆR l (s)=min ĝ kl (s) ; 1=2 n 1 ĝ kl (s) 2 ; (35) R l (s)=min g kl (s) ĝ kl (s) ; n ( ) 1=2 g kl (s) ĝ kl (s) 2 (36) for each l and for all s on the Nyquist contour. Let the lth robust Gershgorin band of G(s); which is composed of circles centered at ĝ ll (s) with radius R l (s); encircle the critical point ( 1=k l ; 0); N l times counterclockwise. Then the negative feedback system with return ratio, G(s)K, is stable if and only if N l = P 0 ; (37) where P 0 is the number of unstable poles of G(s). Proof. Let N and N be the number of counterclockwise encirclements of the origin made by det[i + G(s)K] and by det[k 1 + G(s)], respectively. We must have N = N (Maciejowski, 1989). Because the lth robust Gershgorin band of G(s) encircles the point ( 1=k l ; 0); N l times counterclockwise, this means that the lth robust Gershgorin band of K 1 + G(s) encircles the origin N l times counterclockwise. Condition (34) ensures that K 1 + G(s) is robust diagonally dominant, so that the assumptions of Lemma 2 are satised for Z(s)=K 1 + G(s). Therefore, N = N l N = N l : (38) From the corollary of Theorem 4:1 on page 141 in Rosenbrock (1974), a necessary and sucient condition for stability is that N = P 0. Therefore, Eq. (37) is the necessary and sucient condition for stability. Remark. For a continuous-time system, G(s) gives the frequency response of the system as s traverses the Nyquist contour. Correspondingly, for a discrete-time system, the matrix G(e j! ) for all! [ ; ] provides
6 472 D. Chen, D.E. Seborg / Automatica 38 (2002) the same system frequency response information. Therefore, the results in Lemma 2 and Theorem 4 are also applicable to discrete-time systems if G(s) is replaced by G(e j! ). Corollary. Consider an n n system G(e j! ) with the estimated frequency response matrix Ĝ(e j! )= [ĝ kl (e j! )] n n and covariance information as described in Section 2. Thus; the robust Gershgorin bands at the (1 )100% condence level are composed of the circles centered at ĝ ll (e j! ) and having radii; R l (!)= ˆR l (!)+ R l (!); dened as ˆR l (!)=min ĝ kl (e j! ) ; n 1 1=2 ĝ kl (e j! ) 2 ; (39) R l (!)=min{r e l; R c l}; (40) R e l(!)= 2 ;p 1;kl (!); (41) R c l(!)= n ;p 2 1;l (!): (42) Let P 0 be the number of unstable open-loop poles which is assumed to remain constant for G(e j! ) and Ĝ(e j! ). A decentralized controller; K = diag[k 1 ;:::;k n ]; can be designed so that the robust Gershgorin bands exclude the critical point ( 1=k l ; 0) but encircle it counterclockwise N l times for all l =1; 2;:::;n. Then the closed-loop system is stable at the (1 )100% condence level if Eq. (37) holds. Proof. Substituting the uncertainty descriptions of Eqs. (16) and (25) in Section 2 into Eqs. (35) and (36), then the radius of each robust Gershgorin circle at the (1 )100% condence level is obtained as Eqs. (39) (42). Therefore, according to Theorem 4 the closed-loop system is stable at the (1 )100% condence level if Eq. (37) holds. Theorem 4 and the corollary provide a simple method for analyzing robust stability for a multivariable system which has an estimated model and statistical uncertainty description obtained from system identication. The upper bound on the magnitude of the controller gain k l for robust stability can be easily obtained from this theorem. Furthermore, robust design methods for decentralized control systems can be developed based on the proposed robust stability theorems. 5. Example Wood and Berry (1973) developed the following empirical model of a pilot-scale distillation column that is used to separate a methanol water mixture, XD (s) = X B (s) 12:8e s 16:7s +1 6:6e 7s 10:9s +1 18:9e 3s 21s +1 19:4e 3s 14:4s +1 R(s) ; (43) S(s) where X D and X B are the overhead and bottom compositions of methanol, respectively, R is the reux ow rate, and S is the steam ow rate to the reboiler. A MIMO PRBS signal was designed and used as the excitation signal for process identication. Gaussian measurement noise was added to provide a signal-to-noise ratio of three. The sampling period was chosen as t =1 min. A subspace identication technique, canonical variate analysis (CVA) (Larimore, 1996), and the ADAPTX software were used to identify a state space model from a 1200-point data set. The frequency response and the corresponding covariance for each element are also calculated by the software. Fig. 1 compares the estimated model with the true model and 95% condence intervals. Based on the estimated frequency response and covariance information, the uncertainty bounds for the radii of the Gershgorin bands were calculated. The 95% condence bounds for the Gershgorin band radii are shown in Fig. 2. Notice that the element bounds are more conservative than the column bounds. By applying Theorem 4, a proportional-only, decentralized controller, K = diag[0:3; 0:16], has been obtained that guarantees robust stability for all the possible systems that are within the 95% condence bands of the estimated model. Fig. 3 shows that the Gershgorin bands for this controller do not quite cover the critical point. Thus, it can be concluded that any controller with smaller gains (in absolute value) can guarantee robust stability, if this control structure is used. 6. Conclusions A robust Nyquist array analysis for MIMO systems has been developed based on uncertainty descriptions that are readily available from system identication techniques. Two types of uncertainty bounds on the system frequency response are derived: an uncertainty bound for each element of the system frequency response and an uncertainty bound for each column of the system frequency response. Gershgorin s theorem is extended to include model uncertainty. Robust diagonal dominance and robust Gershgorin bands are dened for uncertain system
7 D. Chen, D.E. Seborg / Automatica 38 (2002) Fig. 1. The frequency response with 95% condence bands. Fig. 2. The uncertainty bounds for the Gershgorin band radii. Fig. 3. Gershgorin bands with proportional-only, decentralized control. models. Based on these results, new robust stability theorems are developed for systems with model uncertainty. The robust Nyquist array analysis is illustrated in an application to the widely used, Wood Berry distillation column model. Acknowledgements The authors acknowledge the UCSB Process Control Consortium for nancial support and Dr. Wallace Larimore (Adaptics, Inc.) for his assistance and for providing
8 474 D. Chen, D.E. Seborg / Automatica 38 (2002) the ADAPTx software. The reviewers suggested a number of useful modications and references which have been included. Appendix Proof of Theorem 3. Let be any eigenvalue of Z. Then there is at least one nonnull x such that x T Z = x T : (A.1) Suppose the lth component of x has the largest modulus, so x can be normalized as [x 1 ;:::;x l 1 ; 1;x l+1 ;:::;x n ] T with x k 6 1 for k l. Equating the lth element on each side of (A.1) gives x k z kl = x l = z ll = x k z kl : (A.2) Hence ẑ ll (z ll ẑ ll ) = ẑ kl x k + ẑ ll 6 (z kl ẑ kl )x k (A.3) ẑ kl x k + (z kl ẑ kl )x k : (A.4) By using the triangle inequality and the Cauchy Schwarz inequality (Anton, 1987), two inequalities can be obtained for the rst term on the right-hand side, ẑ kl x k 6 ẑ kl x k = ẑ kl x k 6 ẑ kl x k 6 6 n 1 ẑ kl ; (A.5) ẑ kl 2 x k 2 ẑ kl 2 1=2 1=2 : (A.6) Therefore, the following inequality must be held: ẑ kl x k 6 min ẑ kl ; n 1 1=2 ẑ kl 2 : (A.7) Similarly, the following inequality can be obtained for the second term on the right-hand side: (z kl ẑ kl )x k ( ) 1=2 6 min z kl ẑ kl ; n z kl ẑ kl 2 : (A.8) By combining Eqs. (A.4), (A.7) and (A.8), Eq. (26) can be obtained. Thus Theorem 3 is proved. References Anton, H. (1987). Elementary linear algebra (5th ed.). New York: Wiley. Arkun, Y., Manousiouthakis, B., & Putz, P. (1984). Robust Nyquist array methodology: A new theoretical framework for analysis and design of robust multivariable feedback systems. International Journal of Control, 40, Braatz, R. D. (1994). A reconciliation between quantitative feedback theory and robust multivariable control. Proceedings of 1994 American control conference, Baltimore, MD (pp ). Braatz, R. D., & Crisalle, O. D. (1998). Robustness analysis for systems with ellipsoidal uncertainty. International Journal of Robust and Nonlinear Control, 8, Chait, Y., Chen, Q., & Hollot, C. V. (1999). Automatic loop-shaping of QFT controllers via linear programming. Journal of Dynamic Systems, Measurement and Control-Transactions of the ASME, 121, Chellaboina, V.-S., Haddad, W. M., & Bernstein, D. S. (1998). Structured matrix norms for robust stability and performance with block-structured uncertainty. International Journal of Control, 71, Chen, C. W., Wang, G. G., & Wang, S. H. (1991). New interpretation of a MIMO quantitative feedback theory. Proceedings of 1991 American control conference, Boston, MA (pp ). Cooley, B. L., & Lee, J. H. (1998). Integrated identication and robust control. Journal of Process Control, 8, De Santis, E., & Vicino, A. (1996). Diagonal dominance for robust stability of MIMO interval systems. IEEE Transactions on Automatic Control, 41, Gao, D., & Zhang, L. (1993). The robust inverse Nyquist array (RINA) method for the design of multivariable control systems. Proceedings of IEEE TENCON 93, Beijing, China (pp ). Goodwin, G. C. (1999). Identication and robust control: Bridging the gap. Seventh Mediterranean conference, Israel (pp. 1 22). Horowitz, I. M. (1982). Quantitative feedback theory. IEE Proceedings D, 129, Horowitz, I. M. (1992). Quantitative feedback design theory (QFT). Boulder, CO: QFT Publications. Jayasuriya, S., & Zhao, Y. (1994). Robust stability of plants with mixed uncertainties and quantitative feedback theory. Journal of Dynamic Systems, Measurement and Control-Transactions of the ASME, 116, Khatri, S., & Parrilo, P. A. (1998). Spherical. Proceedings of 1998 American control conference, Philadelphia, PA (pp ). Kontogiannis, E., & Munro, N. (1997). Extreme point solutions to the diagonal dominance problem and stability analysis of uncertain systems. Proceedings of 1997 American control conference, Albuquerque, NM (pp ).
9 D. Chen, D.E. Seborg / Automatica 38 (2002) Larimore, W. E. (1996). ADAPTX users manual. Mclean, VA: Adaptics Inc. Lee, J. W., Chait, Y., & Steinbuch, M. (2000). On QFT tuning of multivariable controllers. Automatica, 36, Ljung, L. (1985). Asymptotic variance expressions for identied black-box transfer function models. IEEE Transactions on Automatic Control, AC-30, Maciejowski, J. M. (1989). Multivariable feedback design. Reading, MA: Addison-Wesley. Nwokah, O. D. I., Nordgren, R. E., & Grewal, G. S. (1995). Inverse Nyquist array: A quantitative theory. IEE Proceedings D, 142, Rosenbrock, H. H. (1974). Computer-aided control system design. New York: Academic Press. Rugh, W. J. (1996). Linear system theory (2nd ed.). Upper Saddle River, NJ: Prentice-Hall. Skogestad, S., & Postlethwaite, I. (1996). Multivariable feedback control: Analysis and design. New York: Wiley. Wahlberg, B., & Ljung, L. (1992). Hard frequency-domain model error bounds from least-squares like identication techniques. IEEE Transactions on Automatic Control, 37, Wood, R. K., & Berry, M. W. (1973). Terminal composition control of a binary distillation column. Chemical Engineering Science, 28, Zhou, K., Doyle, J. C., & Glover, K. (1996). Robust and optimal control. Upper Saddle River, NJ: Prentice-Hall. Dale E. Seborg received his B.S. degree from the University of Wisconsin and his Ph.D. degree from Princeton University, both in chemical engineering. From 1968 to 1977, he was a faculty member at the University of Alberta in Edmonton, Canada. He joined UCSB in 1977 and served as the department chair for three years. Dr. Seborg s research interests are in the areas of process control and monitoring. He is the co-author of a widely used textbook, Process Dynamics and Control (1989), with Prof. Duncan Mellichamp (UCSB) and Prof. Tom Edgar (UT-Austin). This textbook has been translated into Korean and Japanese. Dr. Seborg is the co-editor of three books and has taught a variety of short courses around the world. He is an active industrial consultant. Dr. Seborg is the recipient, or co-recipient, of several national awards that include the American Statistical Association s Statistics in Chemistry Award (1994), the American Automatic Control Council s Education Award (1993), the American Society of Engineering Education s Meriam-Wiley Distinguished Author Award (1990), and the Joint Automatic Control Conference Best Paper Award (1973). Within IFAC, he was the NOC Co-Chair for SYSID-2000 and the IPC Chair for ADCHEM He also served as the General Chair for the 1992 American Control Conference and co-organized the 1981 Chemical Process Control (CPC-2) Conference. He has been a director of the American Automatic Control Council and the AIChE CAST Division. He currently serves on the editorial board of the IEE Proc. on Control Theory and Applications. Previously, he was an Associate Editor at Large for the IEEE Trans. on Automatic Control, and a member of the editorial boards for the Int. J. of HVAC & R Research and the Int. J. of Adaptive Control & Signal Processing. Dan Chen received her B.S. degree in Chemical Engineering and her M.S. degree in Industrial Process Control from Zhejiang University, China, in 1993 and She is currently a Ph.D. student in Department of Chemical Engineering at the University of California, Santa Barbara. Her research interests include robustness analysis by utilizing system identication and statistical uncertainty descriptions, decentralized control system design and analysis.
CONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. VIII - Design Techniques in the Frequency Domain - Edmunds, J.M. and Munro, N.
DESIGN TECHNIQUES IN THE FREQUENCY DOMAIN Edmunds, Control Systems Center, UMIST, UK Keywords: Multivariable control, frequency response design, Nyquist, scaling, diagonal dominance Contents 1. Frequency
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More information80DB A 40DB 0DB DB
Stability Analysis On the Nichols Chart and Its Application in QFT Wenhua Chen and Donald J. Ballance Centre for Systems & Control Department of Mechanical Engineering University of Glasgow Glasgow G12
More informationThe Generalized Nyquist Criterion and Robustness Margins with Applications
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA The Generalized Nyquist Criterion and Robustness Margins with Applications Abbas Emami-Naeini and Robert L. Kosut Abstract
More informationAnalysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems
Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA
More informationFilter Design for Linear Time Delay Systems
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 11, NOVEMBER 2001 2839 ANewH Filter Design for Linear Time Delay Systems E. Fridman Uri Shaked, Fellow, IEEE Abstract A new delay-dependent filtering
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 5: Uncertainty and Robustness in SISO Systems [Ch.7-(8)] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 5:Uncertainty and Robustness () FEL3210 MIMO
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 6: Robust stability and performance in MIMO systems [Ch.8] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 6: Robust Stability and Performance () FEL3210
More informationDESPITE SEVERAL decades of industrial PID control
1270 IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 17, NO. 6, NOVEMBER 2009 Tuning of Multivariable Decentralized Controllers Through the Ultimate-Point Method Lucíola Campestrini, Luiz Carlos
More informationChapter Stability Robustness Introduction Last chapter showed how the Nyquist stability criterion provides conditions for the stability robustness of
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter Stability
More informationComputation of Stabilizing PI and PID parameters for multivariable system with time delays
Computation of Stabilizing PI and PID parameters for multivariable system with time delays Nour El Houda Mansour, Sami Hafsi, Kaouther Laabidi Laboratoire d Analyse, Conception et Commande des Systèmes
More informationMIMO Toolbox for Matlab
MIMO Toolbox for Matlab Oskar Vivero School of Electric and Electronic Engineering University of Manchester, Manchester, UK M60 1QD Email: oskar.vivero@postgrad.manchester.ac.uk Jesús Liceaga-Castro Departamento
More informationNull controllable region of LTI discrete-time systems with input saturation
Automatica 38 (2002) 2009 2013 www.elsevier.com/locate/automatica Technical Communique Null controllable region of LTI discrete-time systems with input saturation Tingshu Hu a;, Daniel E. Miller b,liqiu
More informationOptimization based robust control
Optimization based robust control Didier Henrion 1,2 Draft of March 27, 2014 Prepared for possible inclusion into The Encyclopedia of Systems and Control edited by John Baillieul and Tariq Samad and published
More informationAsymptotic behaviour of Toeplitz matrix in multi-input multi-output model predictive control
23 European Control Conference ECC) July 7-9, 23, Zürich, Switzerland Asymptotic behaviour of Toeplitz matrix in multi-input multi-output model predictive control Quang Tran, Leyla Özkan, Jobert Ludlage
More informationEects of small delays on stability of singularly perturbed systems
Automatica 38 (2002) 897 902 www.elsevier.com/locate/automatica Technical Communique Eects of small delays on stability of singularly perturbed systems Emilia Fridman Department of Electrical Engineering
More informationExploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability
Exploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability Nuntanut Raksasri Jitkomut Songsiri Department of Electrical Engineering, Faculty of Engineering,
More informationNew Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems
Systems & Control Letters 43 (21 39 319 www.elsevier.com/locate/sysconle New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems E. Fridman Department of Electrical
More informationChapter Robust Performance and Introduction to the Structured Singular Value Function Introduction As discussed in Lecture 0, a process is better desc
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter Robust
More informationNP-hardness of the stable matrix in unit interval family problem in discrete time
Systems & Control Letters 42 21 261 265 www.elsevier.com/locate/sysconle NP-hardness of the stable matrix in unit interval family problem in discrete time Alejandra Mercado, K.J. Ray Liu Electrical and
More informationA New Subspace Identification Method for Open and Closed Loop Data
A New Subspace Identification Method for Open and Closed Loop Data Magnus Jansson July 2005 IR S3 SB 0524 IFAC World Congress 2005 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems
More informationFrequency-Domain Robust Control Toolbox
Frequency-Domain Robust Control Toolbox Alireza Karimi Abstract A new frequency-domain robust control toolbox is introduced and compared with some features of the robust control toolbox of Matlab. A summary
More informationNorm invariant discretization for sampled-data fault detection
Automatica 41 (25 1633 1637 www.elsevier.com/locate/automatica Technical communique Norm invariant discretization for sampled-data fault detection Iman Izadi, Tongwen Chen, Qing Zhao Department of Electrical
More informationRobust cross-directional control of paper making machines with saturating actuators
Robust cross-directional control of paper making machines with saturating actuators J. A. G. Akkermans, A. G. Wills and W. P. Heath March 2, 2004 1 Introduction Preface Technical Report EE03044 In this
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (2011) 797 802 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: wwwelseviercom/locate/aml Model order determination using the Hankel
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationRobust fixed-order H Controller Design for Spectral Models by Convex Optimization
Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by
More informationRobust Performance Example #1
Robust Performance Example # The transfer function for a nominal system (plant) is given, along with the transfer function for one extreme system. These two transfer functions define a family of plants
More informationDECENTRALIZED PI CONTROLLER DESIGN FOR NON LINEAR MULTIVARIABLE SYSTEMS BASED ON IDEAL DECOUPLER
th June 4. Vol. 64 No. 5-4 JATIT & LLS. All rights reserved. ISSN: 99-8645 www.jatit.org E-ISSN: 87-395 DECENTRALIZED PI CONTROLLER DESIGN FOR NON LINEAR MULTIVARIABLE SYSTEMS BASED ON IDEAL DECOUPLER
More informationMRAGPC Control of MIMO Processes with Input Constraints and Disturbance
Proceedings of the World Congress on Engineering and Computer Science 9 Vol II WCECS 9, October -, 9, San Francisco, USA MRAGPC Control of MIMO Processes with Input Constraints and Disturbance A. S. Osunleke,
More informationMULTILOOP PI CONTROLLER FOR ACHIEVING SIMULTANEOUS TIME AND FREQUENCY DOMAIN SPECIFICATIONS
Journal of Engineering Science and Technology Vol. 1, No. 8 (215) 113-1115 School of Engineering, Taylor s University MULTILOOP PI CONTROLLER FOR ACHIEVING SIMULTANEOUS TIME AND FREQUENCY DOMAIN SPECIFICATIONS
More informationA characterization of consistency of model weights given partial information in normal linear models
Statistics & Probability Letters ( ) A characterization of consistency of model weights given partial information in normal linear models Hubert Wong a;, Bertrand Clare b;1 a Department of Health Care
More informationRobust Loop Shaping Controller Design for Spectral Models by Quadratic Programming
Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Gorka Galdos, Alireza Karimi and Roland Longchamp Abstract A quadratic programming approach is proposed to tune fixed-order
More informationAPPLICATION OF D-K ITERATION TECHNIQUE BASED ON H ROBUST CONTROL THEORY FOR POWER SYSTEM STABILIZER DESIGN
APPLICATION OF D-K ITERATION TECHNIQUE BASED ON H ROBUST CONTROL THEORY FOR POWER SYSTEM STABILIZER DESIGN Amitava Sil 1 and S Paul 2 1 Department of Electrical & Electronics Engineering, Neotia Institute
More information1 Introduction 198; Dugard et al, 198; Dugard et al, 198) A delay matrix in such a lower triangular form is called an interactor matrix, and almost co
Multivariable Receding-Horizon Predictive Control for Adaptive Applications Tae-Woong Yoon and C M Chow y Department of Electrical Engineering, Korea University 1, -a, Anam-dong, Sungbu-u, Seoul 1-1, Korea
More informationLecture 6. Chapter 8: Robust Stability and Performance Analysis for MIMO Systems. Eugenio Schuster.
Lecture 6 Chapter 8: Robust Stability and Performance Analysis for MIMO Systems Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 6 p. 1/73 6.1 General
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationChapter 9 Robust Stability in SISO Systems 9. Introduction There are many reasons to use feedback control. As we have seen earlier, with the help of a
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Robust
More informationreality is complex process
ISS0080 Automation and Process Control Lecture 5 1 Process models the desire to describe reality Model of the process, model simplication, identication. model reality is complex process Replaces the original;
More informationMANY adaptive control methods rely on parameter estimation
610 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 52, NO 4, APRIL 2007 Direct Adaptive Dynamic Compensation for Minimum Phase Systems With Unknown Relative Degree Jesse B Hoagg and Dennis S Bernstein Abstract
More informationRobust Gain Scheduling Synchronization Method for Quadratic Chaotic Systems With Channel Time Delay Yu Liang and Horacio J.
604 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 56, NO. 3, MARCH 2009 Robust Gain Scheduling Synchronization Method for Quadratic Chaotic Systems With Channel Time Delay Yu Liang
More informationA Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 12, NO. 5, SEPTEMBER 2001 1215 A Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing Da-Zheng Feng, Zheng Bao, Xian-Da Zhang
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationTWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN INTERVAL SYSTEM USING GENETIC ALGORITHM
International Journal of Electrical and Electronics Engineering Research (IJEEER) ISSN:2250-155X Vol.2, Issue 2 June 2012 27-38 TJPRC Pvt. Ltd., TWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN
More informationAdvanced Aerospace Control. Marco Lovera Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano
Advanced Aerospace Control Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano ICT for control systems engineering School of Industrial and Information Engineering Aeronautical Engineering
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationA FEEDBACK STRUCTURE WITH HIGHER ORDER DERIVATIVES IN REGULATOR. Ryszard Gessing
A FEEDBACK STRUCTURE WITH HIGHER ORDER DERIVATIVES IN REGULATOR Ryszard Gessing Politechnika Śl aska Instytut Automatyki, ul. Akademicka 16, 44-101 Gliwice, Poland, fax: +4832 372127, email: gessing@ia.gliwice.edu.pl
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationIan G. Horn, Jeffery R. Arulandu, Christopher J. Gombas, Jeremy G. VanAntwerp, and Richard D. Braatz*
Ind. Eng. Chem. Res. 996, 35, 3437-344 3437 PROCESS DESIGN AND CONTROL Improved Filter Design in Internal Model Control Ian G. Horn, Jeffery R. Arulandu, Christopher J. Gombas, Jeremy G. VanAntwerp, and
More informationHANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING
Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain HANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING Björn Wittenmark Department of Automatic Control Lund Institute of Technology
More informationRobust Control. 2nd class. Spring, 2018 Instructor: Prof. Masayuki Fujita (S5-303B) Tue., 17th April, 2018, 10:45~12:15, S423 Lecture Room
Robust Control Spring, 2018 Instructor: Prof. Masayuki Fujita (S5-303B) 2nd class Tue., 17th April, 2018, 10:45~12:15, S423 Lecture Room 2. Nominal Performance 2.1 Weighted Sensitivity [SP05, Sec. 2.8,
More informationAn LQ R weight selection approach to the discrete generalized H 2 control problem
INT. J. CONTROL, 1998, VOL. 71, NO. 1, 93± 11 An LQ R weight selection approach to the discrete generalized H 2 control problem D. A. WILSON², M. A. NEKOUI² and G. D. HALIKIAS² It is known that a generalized
More informationPROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS
PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix
More informationFeedback Control of Dynamic Systems
THIRD EDITION Feedback Control of Dynamic Systems Gene F. Franklin Stanford University J. David Powell Stanford University Abbas Emami-Naeini Integrated Systems, Inc. TT Addison-Wesley Publishing Company
More informationMULTILOOP CONTROL APPLIED TO INTEGRATOR MIMO. PROCESSES. A Preliminary Study
MULTILOOP CONTROL APPLIED TO INTEGRATOR MIMO PROCESSES. A Preliminary Study Eduardo J. Adam 1,2*, Carlos J. Valsecchi 2 1 Instituto de Desarrollo Tecnológico para la Industria Química (INTEC) (Universidad
More informationA new robust delay-dependent stability criterion for a class of uncertain systems with delay
A new robust delay-dependent stability criterion for a class of uncertain systems with delay Fei Hao Long Wang and Tianguang Chu Abstract A new robust delay-dependent stability criterion for a class of
More informationQuantitative Feedback Theory based Controller Design of an Unstable System
Quantitative Feedback Theory based Controller Design of an Unstable System Chandrima Roy Department of E.C.E, Assistant Professor Heritage Institute of Technology, Kolkata, WB Kalyankumar Datta Department
More informationIntroduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31
Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured
More informationFair and Square Computation of Inverse -Transforms of Rational Functions
IEEE TRANSACTIONS ON EDUCATION, VOL. 55, NO. 2, MAY 2012 285 Fair and Square Computation of Inverse -Transforms of Rational Functions Marcos Vicente Moreira and João Carlos Basilio Abstract All methods
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
More informationControl Configuration Selection for Multivariable Descriptor Systems
Control Configuration Selection for Multivariable Descriptor Systems Hamid Reza Shaker and Jakob Stoustrup Abstract Control configuration selection is the procedure of choosing the appropriate input and
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationOptimal triangular approximation for linear stable multivariable systems
Proceedings of the 007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 007 Optimal triangular approximation for linear stable multivariable systems Diego
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationH State-Feedback Controller Design for Discrete-Time Fuzzy Systems Using Fuzzy Weighting-Dependent Lyapunov Functions
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL 11, NO 2, APRIL 2003 271 H State-Feedback Controller Design for Discrete-Time Fuzzy Systems Using Fuzzy Weighting-Dependent Lyapunov Functions Doo Jin Choi and PooGyeon
More informationExperimental evidence showing that stochastic subspace identication methods may fail 1
Systems & Control Letters 34 (1998) 303 312 Experimental evidence showing that stochastic subspace identication methods may fail 1 Anders Dahlen, Anders Lindquist, Jorge Mari Division of Optimization and
More informationCompensatorTuning for Didturbance Rejection Associated with Delayed Double Integrating Processes, Part II: Feedback Lag-lead First-order Compensator
CompensatorTuning for Didturbance Rejection Associated with Delayed Double Integrating Processes, Part II: Feedback Lag-lead First-order Compensator Galal Ali Hassaan Department of Mechanical Design &
More information1.1 Notations We dene X (s) =X T (;s), X T denotes the transpose of X X>()0 a symmetric, positive denite (semidenite) matrix diag [X 1 X ] a block-dia
Applications of mixed -synthesis using the passivity approach A. Helmersson Department of Electrical Engineering Linkoping University S-581 83 Linkoping, Sweden tel: +46 13 816 fax: +46 13 86 email: andersh@isy.liu.se
More informationWhile using the input and output data fu(t)g and fy(t)g, by the methods in system identification, we can get a black-box model like (In the case where
ESTIMATE PHYSICAL PARAMETERS BY BLACK-BOX MODELING Liang-Liang Xie Λ;1 and Lennart Ljung ΛΛ Λ Institute of Systems Science, Chinese Academy of Sciences, 100080, Beijing, China ΛΛ Department of Electrical
More informationRobust Control. 1st class. Spring, 2017 Instructor: Prof. Masayuki Fujita (S5-303B) Tue., 11th April, 2017, 10:45~12:15, S423 Lecture Room
Robust Control Spring, 2017 Instructor: Prof. Masayuki Fujita (S5-303B) 1st class Tue., 11th April, 2017, 10:45~12:15, S423 Lecture Room Reference: [H95] R.A. Hyde, Aerospace Control Design: A VSTOL Flight
More informationNON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS
Control 4, University of Bath, UK, September 4 ID-83 NON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS H. Yue, H. Wang Control Systems Centre, University of Manchester
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationModel Uncertainty and Robust Stability for Multivariable Systems
Model Uncertainty and Robust Stability for Multivariable Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Representing model uncertainty.
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationDecoupled Feedforward Control for an Air-Conditioning and Refrigeration System
American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, FrB1.4 Decoupled Feedforward Control for an Air-Conditioning and Refrigeration System Neera Jain, Member, IEEE, Richard
More informationRobust Tuning of Power System Stabilizers Using Coefficient Diagram Method
International Journal of Electrical Engineering. ISSN 0974-2158 Volume 7, Number 2 (2014), pp. 257-270 International Research Publication House http://www.irphouse.com Robust Tuning of Power System Stabilizers
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationAvailable online at ScienceDirect. Procedia Engineering 100 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 100 (015 ) 345 349 5th DAAAM International Symposium on Intelligent Manufacturing and Automation, DAAAM 014 Control of Airflow
More informationRobust sampled-data stabilization of linear systems: an input delay approach
Automatica 40 (004) 44 446 www.elsevier.com/locate/automatica Technical Communique Robust sampled-data stabilization of linear systems: an input delay approach Emilia Fridman a;, Alexandre Seuret b, Jean-Pierre
More informationDiscrete-Time H Gaussian Filter
Proceedings of the 17th World Congress The International Federation of Automatic Control Discrete-Time H Gaussian Filter Ali Tahmasebi and Xiang Chen Department of Electrical and Computer Engineering,
More informationA Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case
A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of
More informationAN OPTIMIZATION-BASED APPROACH FOR QUASI-NONINTERACTING CONTROL. Jose M. Araujo, Alexandre C. Castro and Eduardo T. F. Santos
ICIC Express Letters ICIC International c 2008 ISSN 1881-803X Volume 2, Number 4, December 2008 pp. 395 399 AN OPTIMIZATION-BASED APPROACH FOR QUASI-NONINTERACTING CONTROL Jose M. Araujo, Alexandre C.
More information(Continued on next page)
(Continued on next page) 18.2 Roots of Stability Nyquist Criterion 87 e(s) 1 S(s) = =, r(s) 1 + P (s)c(s) where P (s) represents the plant transfer function, and C(s) the compensator. The closedloop characteristic
More informationLocal Modelling with A Priori Known Bounds Using Direct Weight Optimization
Local Modelling with A Priori Known Bounds Using Direct Weight Optimization Jacob Roll, Alexander azin, Lennart Ljung Division of Automatic Control Department of Electrical Engineering Linköpings universitet,
More informationThe Uniformity Principle: A New Tool for. Probabilistic Robustness Analysis. B. R. Barmish and C. M. Lagoa. further discussion.
The Uniformity Principle A New Tool for Probabilistic Robustness Analysis B. R. Barmish and C. M. Lagoa Department of Electrical and Computer Engineering University of Wisconsin-Madison, Madison, WI 53706
More informationFloor Control (kn) Time (sec) Floor 5. Displacement (mm) Time (sec) Floor 5.
DECENTRALIZED ROBUST H CONTROL OF MECHANICAL STRUCTURES. Introduction L. Bakule and J. Böhm Institute of Information Theory and Automation Academy of Sciences of the Czech Republic The results contributed
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XII - Lyapunov Stability - Hassan K. Khalil
LYAPUNO STABILITY Hassan K. Khalil Department of Electrical and Computer Enigneering, Michigan State University, USA. Keywords: Asymptotic stability, Autonomous systems, Exponential stability, Global asymptotic
More informationNumerical Solution for Multivariable Idle Speed Control of a Lean Burn Natural Gas Engine
Numerical Solution for Multivariable Idle Speed Control of a Lean Burn Natural Gas Engine Arvind Sivasubramanian, Peter H. Meckl Abstract This paper proposes a numerical solution to the idle speed and
More informationQFT design for uncertain non-minimum phase and unstable plants revisited
Loughborough University Institutional Repository QFT design for uncertain non-minimum phase and unstable plants revisited This item was submitted to Loughborough University's Institutional Repository by
More informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationA NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN
Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain A NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN Chyi Hwang,1 Chun-Yen Hsiao Department of Chemical Engineering National
More informationAnalysis of SISO Control Loops
Chapter 5 Analysis of SISO Control Loops Topics to be covered For a given controller and plant connected in feedback we ask and answer the following questions: Is the loop stable? What are the sensitivities
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 17: Robust Stability Readings: DDV, Chapters 19, 20 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology April 6, 2011 E. Frazzoli
More informationAn Iteration-Domain Filter for Controlling Transient Growth in Iterative Learning Control
21 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 21 WeC14.1 An Iteration-Domain Filter for Controlling Transient Growth in Iterative Learning Control Qing Liu and Douglas
More informationFEEDFORWARD CONTROLLER DESIGN BASED ON H ANALYSIS
271 FEEDFORWARD CONTROLLER DESIGN BASED ON H ANALYSIS Eduardo J. Adam * and Jacinto L. Marchetti Instituto de Desarrollo Tecnológico para la Industria Química (Universidad Nacional del Litoral - CONICET)
More informationAdditional Closed-Loop Frequency Response Material (Second edition, Chapter 14)
Appendix J Additional Closed-Loop Frequency Response Material (Second edition, Chapter 4) APPENDIX CONTENTS J. Closed-Loop Behavior J.2 Bode Stability Criterion J.3 Nyquist Stability Criterion J.4 Gain
More informationThe model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho
Model Reduction from an H 1 /LMI perspective A. Helmersson Department of Electrical Engineering Linkoping University S-581 8 Linkoping, Sweden tel: +6 1 816 fax: +6 1 86 email: andersh@isy.liu.se September
More information