Chapter Robust Performance and Introduction to the Structured Singular Value Function Introduction As discussed in Lecture 0, a process is better desc
|
|
- Cameron Harrington
- 6 years ago
- Views:
Transcription
1 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
2 Chapter Robust Performance and Introduction to the Structured Singular Value Function Introduction As discussed in Lecture 0, a process is better described in terms of a set of plants centered around a nominal model The robust stabilization problem is concerned with nding non conservative conditions on the stable nominal closed loop system that guarantee the stability of all possible closed loop systems An equally important problem is the robust performance problem which is concerned with nding non conservative conditions on the nominal closed loop system that guarnatee that the performance is met for all possible closed loop systems Robust Disturbance Rejection We will focus our discussion on one prototype problem, namely, the robust disturbance rejection problem shown in Figure This motivates the following problem: Robust Disturbance Rejection Problem (RP) Find conditions on the nominal closed-loop system (P o K) such that K robustly stabilizes all P, where fp j P (I W )P o kk < g: k(i P K) ; W k for all P From Lecture 0, a performance objective in terms of the H -norm of some closed loop map between some exogenous input w, to a regulated variable z, is mathematically equivalent to a robust stabilization problem with a perturbation block mapping the regulated output z to the exogenous input w Obviously, the new perturbed system is stable if and only if kt zw k, which is the performance
3 - - - m to k () m e? - W - W m - -?? P 0 - ; m - - K Figure : Uncertain Plant with Disturbance objective Notice that if the performance objective consists of several closed loop maps, then several perturbation blocks can be introduced in exactly the same fashion Proceeding for RP, we can \wrap" a frequency-weighted perturbation from the output to the input of interest, which results in the model of Figure Next, we can re-arrange the system into the - W z - w w z W?? m - - m P o ; K Figure : Robust Performance Model M- feedback form (a nominal stable M in feedback with the perturbation ) as in Figure In this case, however, there are multiple inputs and outputs to consider We use the following procedure to generate M and : Dene w i z i to be the output and input, respectively, of the perturbation i For a total of m perturbations, compute the matrix transfer function M as the map from w : w z z w m z m In other words, all the blocks are removed, and the transfer functions \seen" by the blocks from each input w j to each output z i are calculated and used as the (i j) th element of M The perturbation matrix will have the structure m i k < : () For a SISO system, each i is a scalar, so that becomes a diagonal matrix with complex entries In the MIMO case, is block-diagonal
4 Example (Robust Disturbance Rejection) Applying the robust performance procedure to Figure yields: ;W (I P 0 K) ; P 0 K ;W (I P 0 K) ; P 0 K M : () W (I P 0 K) ; W (I P 0 K) ; The transfer functions on the diagonal are identical to those in the single-block robust stability and disturbance-rejection problems, respectively, while the o-diagonal terms account for the interaction between the two constraints Having found the appropriate M and, we have thereby reduced the robust performance problem to a stability problem for the system of Figure - l - M - m! Figure : M- Feedback Form A sucient condition for robust stability is given by the small gain theorem, namely, max [M(jw)] max [(jw)] < for all w: Since is norm bounded by one, this condition translates to kmk This condition, however, is far from necessary since has a block diagonal structure The Structured Singular Value For an unstructured perturbation, the supremum of the maximum singular value of M (ie kmk ) provides a clean and numerically tractable method for evaluating robust stability Recall that, for the standard M- loop, the system fails to be robustly stable if there exists an admissible such that (I ; M) is singular What distinguishes the current situation from the unstructured case is that we have placed constraints on the set Given this more limited set of admissible perturbations, we desire a measure of robust stability similar to kmk This can be derived from the structured singular value (M) Denition The structured singular value of a complex matrix M with respect to a class of perturbations is given by (M) : () inff max () j det(i ; M) 0g If det(i ; M) 0 for all, then (M) 0
5 Theorem The M- System is stable for all with kk < if and only if sup (M) :! Proof: Immediate, from the denition Clearly, if, then the norm of the smallest allowable destabilizing perturbation must by denition be greater than Properties of the Structured Singular Value It is important to note that is a function that depends on the perturbation class (sometimes, this function is denoted by to indicate this dependence) The following are useful properties of such a fucntion (M) 0 If fi j C g, then (M) (M), the spectral radius of M (which is equal to the magnitude of the eigenvalue of M with maximum magnitude) If f j is an arbitrary complex matrixg then max (M), from which sup! kmk Property shows that the spectral radius function is a particular function with respect to a perturbation class consisting of matrices of the form of scaled identity Property shows that the maximum singular value function is a particular function with respect to a perturbation class consisting of arbitrary norm bounded perturbations (no structural constraints) If fdiag( : : : n ) j i complexg, then (M) (D ; MD) for any D diag(d : : : d n ) jd i j > 0 The set of such scales is denoted D This can be seen by noting that det(i ; AB) det(i ; BA), so that det(i ; D ; MD) det(i ; MDD ; ) det(i ; M) The last equality arises since the diagonal matrices and D commute If diag( : : : n ) i complex, then (M) (M) max (M) This property follows from the following observation: If, then It is clear that the class of perturbations consisting of scaled identity matrices is a subset of which is a subset of the class of all unstructured perturbations From and we have that (M) (D ; MD) inf DD max (D ; MD) Computation of In general, there is no closed-form method for computing Upper and lower bounds may be computed and rened, however In these notes we will only be concerned with computing the upper bound If diag( : : : n ), then the upper bound on is something that is easy to calculate Furthermore, property above suggests that by inmizing max (D ; MD) over all possible diagonal scaling matrices, we obtain a better approximation of This turns out to be a convex optimization problem at each
6 () frequency, so that by inmizing over D at each frequency, the tightest upper bound over the set of D may be found for We may then ask when (if ever) this bound is tight In other words, when is it truly a least upper bound The answer is that for three or fewer 's, the bound is tight The proof of this is involved, and is beyond the scope of this class Unfortunately, for four or more perturbations, the bound is not tight, and there is no known method for computing exactly for more than three perturbations Robust Disturbance Rejection (SISO) As shown earlier, the disturbance rejection requirement could be converted to a robust stability problem with two blocks of uncertainty, as in Figure, where and are SISO stable systems Hence is the set of diagonal complex matrices (which result from evaluating at each frequency) Now, since this is a two-block problem, it should be possible to nd by inmizing max (D ; MD) We have D diag(d d ), so that (M) inf d d >0 8 >< >: W ; P 0 K d ; W K > P 0 K d P 0 K max () d W P 0 W d P 0 K P 0 K {z A() with the \pure" robust stability requirement occup ying the upper left diagonal, and the nominal performance requirement on the lower right Setting d d and xing!, and taking the denition of A() from (), we have Now, for nominal performance, we require that For robust stability, we nee 9 } > (M) inf f (A ()A())g: () max jj>0 W P 0 K : W P 0 K : P 0 K For robust performance, the necessary and sucient condition is A bit of algebra yields from which we have max (A A) jj () (8) (M) : (9) W K W P 0 K P 0 K W KP 0 P 0 K jj W P 0 P 0 K W P 0 K W inf max (A A) : P 0 K P 0 K (0) ()
7 W (j ω ) - ω o L(j ω )W (j ω ) L(j ω) Figure : Robust Performance/Nyquist Criterion This minimum occurs at jj jw P 0 j jw Kj () which is not equal to in general, so that sup! kmk In other words, is a less conservative measure than kk in this case Once again, there is a graphical interpretation of the SISO robust disturbance rejection problem, in terms of the Nyquist criterion From (), we have (M) () W P 0 K P 0 K W : () P 0 K Letting L represent the nominal loop gain P 0 K, this can be rewritten as: jw Lj jw j j Lj: () Graphically, we can represent this at each frequency! as a circle centered at ; of radius jw j, and a second circle centered at L of radius jw Lj Robust performance will be achieved as long as the two circles never intersect Loop-shaping Revisited Loop-shaping is a well-established method of control design that concentrates on the frequency-domain characteristics of the open-loop transfer function L P 0 K Based primarily on design experience, there are certain characteristics of the loop transfer function that translate into desirable control performance Other open-loop characteristics are known by experience to result in undesirable or unpredictable behavior This method diers from -synthesis and H methods, which concentrate on optimizing the characteristics of the closed-loop transfer function Since, presumably, a controller with good behavior designed by loop-shaping should be similar in some way to a controller designed by more recent methods, it is of interest to look for parallels in the heuristic rules of loop-shaping and the more methodical methods of -synthesis and H
8 L Identifying the sensitivity and complementary sensitivity functions from (), we can write the RP requirement as jw T j jw Sj : () Model uncertainty typically increases with frequency, so it is important that the complementary sensitivity function decreases with increasing frequency For disturbance rejection, which is typically most critical over a low frequency range, we require that S remain small The weighting functions W and W are designed to reect this, and so might take on the form of Figure Normally, at low W W Figure : Typical Weighting Functions frequency, L >> and at high frequency, L << Now, L T 0 S 0 () L L so that at low frequency, T 0 and S 0 L Thus we can approximate the RP requirement at the low end as: jw j W jw j (8) ) jlj ; jw j At high frequency, the approximation is T 0 L and S 0, which leads to: jw Lj jw j ) jlj ; jw j : (9) jw j These constraints are summarized in Figure, which also notes another design rule, which is that the 0 db crossing should occur at a slope no more negative than -0 db per decade If W and W do not overlap signicantly in frequency, then the upper and lower bounds reduce to jw j and jw j, respectively Example (Loop Shaping) Assume P 0 is minimum phase stable with relative degree Designing a controller by shaping the loop gain L P 0 K is not aected by P 0 just the relative degree is needed
9 W - W 0 db L(j ) -0 db/decade max - W W Figure : Typical Loop-shaping Problem Suppose the multiplicative uncertainty is described by s W 0(0:0s ) ie, the multiplicative perturbations of the plant are upper bounded by W at each frequency The objective is to track sinusoidal signals at the reference input in the frequency range [0 ] rads We would like to make the tracking error small however, we do not know yet by how much Let W have the following frequency response jw j a 0! 0 otherwise Note that this may not correspond to a stable W (s) however, this does not aect the resulting loop shape We are going to exhibit the design by trial and error Let At high frequency,! 0, b L(s) : cs ; jw j L! 0: jw j jw j If we pick c, then the largest value for b such that the above is satised is b 0 Hence 0 L(s) : s At low frequency,!, jw j a jlj : ; jw j ; jw j
10 which implies that a : Checking the RP condition Since jlj is decreasing and jw j is increasing in the range [0 ], the largest a can be solved for: a jl(j)j ; jw (j)j jw Sj jw T j 0:9 8! which implies RP is achieved and the tracking error is smaller than : in the range [0 ] If a better performance is desired, a possibly more complicated L needs to be used The discussion in this chapter has focused on perturbations that are arbitrary dynamic systems This alowed us to think of any class of structured perurbations as sets of arbitrary (structured) matrices at each frequency point These matrices correspond to evaluating the dynamic system at a given frequency In practical applications, some perturbations may be static and not dynamic These arise in problems with real parameter uncertainties We can still proceed as before and transform such problems to the general M- diagram In this case, will have a combination of both static and dynamic perturbations for such a class can be dened as before, and it will provide a necessary and sucient condition for robust stability The main issue here is computing a good upper bound for Of course, we can always embed this class of perturbations in a larger class containing dynamic perturbations and use D-scaling to obtain an upper bound This, however, gives conservative conditions Computing non-conservative upper bounds of for such perturbations remains an active area of research Rank-One Although we do not have methods for computing exactly, there is one particular situation where this is possible This situation occurs if M has rank, ie where a b C is given by However, M ab n Then it follows that with respect to containing complex diagonal perturbations inf f max () j det(i ; M) 0g: (M) det(i ; M) det(i ; ab ) det(i ; b a) 0 det I ; [ n ] ; [ n ] b a b a b na n b a b a b na n C A
11 m (M) N (M) jb m - - D m : and max () max i j i j Hence, inf ] ::: n 8 >< >: max j i j i [ n b a b a b na n : 9 > > Optimizing the RHS, it follows that (verify) P n $ (M) ia i j i X n i jb ia i j Notice that the SISO robust disturbance rejection problem is a rank-one problem This follows since Then M (M) which is the condition we derived before Coprime Factor Perturbations Consider the class of SISO systems N(s) D(s) ;W K P 0 [ ]: P 0 K P 0 K W W P 0 K P 0 K W P 0 K N 0 W D D 0 W k i k < where the nominal plant is N 0 D 0 with the property that both N 0 and D 0 are stable with no common zeros in the RHP Assume that K stabilizes N 0 D 0 This block diagram is shown in Figure - W z - w w z W ;? u? - - y - N 0 ; 0 K Figure : Coprime Factor Perturbation Model
12 M ; : The closed loop block diagram can be mapped to the M- diagram where ; W K W ; K D 0 N 0 K D 0 N 0 K Hence, M has rank and (M) W D 0 N 0 K ; W K D 0 N 0 K W D 0 N 0 K W K D 0 N 0 K W D 0 N 0 K [ ]: W D 0 N 0 K Robust Hurwitz Stability of Polynomials with Complex Perturbations Another application of the structured singular value with rank one matrices is the robust stabil- ity of a family of polynomials with complex perturbations of the coecients In this case let T n; n; : : : 0 and consider the polynomial family P (s ) s n (a n; n; n; )s n; : : : (a ) where a i, i, and i C and j i j We want to obtain a condition that is both necessary and sucient for the Hurwitz stability of the entire family of polynomials P (s ) We can write the polynomials in this family as P (s ) P (s 0) (0) P~(s ) ; n n; n; s a n; n; n; s : s : : : a 0 : : 0 0 () which can also be rewritten as n; 0 0 : : : 0 n; s n; 0 n; 0 : : : 0 n; s n; P (s ) P (s 0) : : : : 0 s 0 0 : : : We assume that the center polynomial P (s 0) is Hurwitz stable This implies that the stability of the entire family P (s ) is equivalent to the condition that : : : P (j! 0) n; 0 0 : : : 0 n; n; 0 n; 0 : : : 0 n; n; : : :
13 @ for all! R and j i j This is equivalent to the condition that 0 B n; n; n; n; for all! R and with kk Now using the concept of the structured singular value we arrive at the following condition which is both necessary and sucient for the Hurwitz stability of the entire family (M) < C A det I : : : 0 P (j! 0) 0 for all! R, where M P (j! 0) n; (j! ) n; n; (j! ) n; 0 : : : : Clearly this is a rank one matrix and by our previous discussion the structured singular value can be computed analytically resulting in the following test for all! R jp (j! 0)j X n i j n;i jj!j n;i <
14 Exercises Exercise In decentralized control, the plant is assumed to be diagonal and controllers are designed independently for each diagonal element If however, the real process is not completely decoupled, the interactions between these separate subsystems can drive the system to instability Consider the plant P P P (s) : P P Assume that P and P are stable and relatively small in comparison to the diagonal elements, and only a bound on their frequency response is available Suppose a controller K diag(k K ) is designed to stabilize the system P 0 diag(p P ) Set-up the problem as a stability robustness problem, ie, put the problem in the M ; form Derive a non-conservative condition (necessary and sucient) that guarantees the stability robustness of the above system Assume the o-diagonal elements are perturbed independently Reduce the result to the simplest form (an answer like (M) < is not acceptable this problem has an exact solution which is computable) How does your answer change if the o-diagonal elements are perturbed simultaneously with the same Exercise Consider the rank problem Suppose, contains only real perturbations Compute the exact expression of (M) Exercise Consider the set of plants characterized by the following sets of numerators and denominators of the transfer function: N(s) N 0 (s) N (s) D(s) D 0 (s) D (s) Where both N 0 and D 0 are polynomials in s, R n, and N D are polynomial row vectors The set of all plants is then given by: Let K be a con troller that stabilizes such that the system is stable N(s) f j R n j i j g D(s) N 0 Compute the exact stability margin ie, compute the largest D 0
15 MIT OpenCourseWare J / 8J Dynamic Systems and Control Spring 0 For information about citing these materials or our Terms of Use, visit:
Chapter 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 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 informationChapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o
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 7 Interconnected
More informationChapter 30 Minimality and Stability of Interconnected Systems 30.1 Introduction: Relating I/O and State-Space Properties We have already seen in Chapt
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
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 informationROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ]
ROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ] General control configuration with uncertainty [8.1] For our robustness analysis we use a system representation in which the uncertain perturbations are
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.24: Dynamic Systems Spring 20 Homework 9 Solutions Exercise 2. We can use additive perturbation model with
More informationChapter 6 Balanced Realization 6. Introduction One popular approach for obtaining a minimal realization is known as Balanced Realization. In this appr
Lectures on ynamic Systems and Control Mohammed ahleh Munther A. ahleh George Verghese epartment of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 6 Balanced
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 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 informationStructured Uncertainty and Robust Performance
Structured Uncertainty and Robust Performance ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Structured uncertainty: motivating example. Structured
More informationStructured singular value and µ-synthesis
Structured singular value and µ-synthesis Robust Control Course Department of Automatic Control, LTH Autumn 2011 LFT and General Framework z P w z M w K z = F u (F l (P,K), )w = F u (M, )w. - Last week
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 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 informationEECE 460. Decentralized Control of MIMO Systems. Guy A. Dumont. Department of Electrical and Computer Engineering University of British Columbia
EECE 460 Decentralized Control of MIMO Systems Guy A. Dumont Department of Electrical and Computer Engineering University of British Columbia January 2011 Guy A. Dumont (UBC EECE) EECE 460 - Decentralized
More informationrad/sec
In this chapter we present Kharitonov's Theorem on robust Hurwitz stability of interval polynomials, dealing with both the real and complex cases This elegant result forms the basis for many of the results,
More informationChapter 3 Least Squares Solution of y = A x 3.1 Introduction We turn to a problem that is dual to the overconstrained estimation problems considered s
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
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 informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
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 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 informationChapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin
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 Observers,
More informationChapter 6 Balanced Realization 6. Introduction One popular approach for obtaining a minimal realization is known as Balanced Realization. In this appr
Lectures on ynamic Systems and Control Mohammed ahleh Munther A. ahleh George Verghese epartment of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 6 Balanced
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 informationhapter 8 Simulation/Realization 8 Introduction Given an nth-order state-space description of the form x_ (t) = f (x(t) u(t) t) (state evolution equati
Lectures on Dynamic Systems and ontrol Mohammed Dahleh Munther Dahleh George Verghese Department of Electrical Engineering and omputer Science Massachuasetts Institute of Technology c hapter 8 Simulation/Realization
More informationCHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER
114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers
More informationRobustness under Bounded Uncertainty with Phase Information A. L. Tits Dept. of Elec. Eng. and Inst. for Syst. Res. University of Maryland College Par
Robustness under Bounded Uncertainty with Phase Information A L Tits Dept of Elec Eng and Inst for Syst Res University of Maryland College Park, MD 74 USA V Balakrishnan School of Elec and Comp Eng Purdue
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 informationDidier HENRION henrion
POLYNOMIAL METHODS FOR ROBUST CONTROL PART I.1 ROBUST STABILITY ANALYSIS: SINGLE PARAMETER UNCERTAINTY Didier HENRION www.laas.fr/ henrion henrion@laas.fr Pont Neuf over river Garonne in Toulouse June
More informationLecture 7 (Weeks 13-14)
Lecture 7 (Weeks 13-14) Introduction to Multivariable Control (SP - Chapters 3 & 4) Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 7 (Weeks 13-14) p.
More informationClosed-loop system 2/1/2016. Generally MIMO case. Two-degrees-of-freedom (2 DOF) control structure. (2 DOF structure) The closed loop equations become
Closed-loop system enerally MIMO case Two-degrees-of-freedom (2 DOF) control structure (2 DOF structure) 2 The closed loop equations become solving for z gives where is the closed loop transfer function
More informationMIMO analysis: loop-at-a-time
MIMO robustness MIMO analysis: loop-at-a-time y 1 y 2 P (s) + + K 2 (s) r 1 r 2 K 1 (s) Plant: P (s) = 1 s 2 + α 2 s α 2 α(s + 1) α(s + 1) s α 2. (take α = 10 in the following numerical analysis) Controller:
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 information3.1 Overview 3.2 Process and control-loop interactions
3. Multivariable 3.1 Overview 3.2 and control-loop interactions 3.2.1 Interaction analysis 3.2.2 Closed-loop stability 3.3 Decoupling control 3.3.1 Basic design principle 3.3.2 Complete decoupling 3.3.3
More informationDesign Methods for Control Systems
Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9
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 informationCDS 101/110a: Lecture 10-1 Robust Performance
CDS 11/11a: Lecture 1-1 Robust Performance Richard M. Murray 1 December 28 Goals: Describe how to represent uncertainty in process dynamics Describe how to analyze a system in the presence of uncertainty
More informationSynthesis of Controllers for Non-minimum Phase and Unstable Systems Using Non-sequential MIMO Quantitative Feedback Theory
Synthesis of Controllers for Non-minimum Phase and Unstable Systems Using Non-sequential MIMO Quantitative Feedback Theory Chen-yang Lan Murray L Kerr Suhada Jayasuriya * Graduate Student Graduate Student
More informationA brief introduction to robust H control
A brief introduction to robust H control Jean-Marc Biannic System Control and Flight Dynamics Department ONERA, Toulouse. http://www.onera.fr/staff/jean-marc-biannic/ http://jm.biannic.free.fr/ European
More informationTHIS paper deals with robust control in the setup associated
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 50, NO 10, OCTOBER 2005 1501 Control-Oriented Model Validation and Errors Quantification in the `1 Setup V F Sokolov Abstract A priori information required for
More informationStability Margin Based Design of Multivariable Controllers
Stability Margin Based Design of Multivariable Controllers Iván D. Díaz-Rodríguez Sangjin Han Shankar P. Bhattacharyya Dept. of Electrical and Computer Engineering Texas A&M University College Station,
More informationClassify a transfer function to see which order or ramp it can follow and with which expected error.
Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,
More informationAn Overview on Robust Control
Advanced Control An Overview on Robust Control P C Scope Keywords Prerequisites allow the student to assess the potential of different methods in robust control without entering deep into theory. Sensitize
More informationControl of Electromechanical Systems
Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance
More informationChapter 8 Stabilization: State Feedback 8. Introduction: Stabilization One reason feedback control systems are designed is to stabilize systems that m
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of echnology c Chapter 8 Stabilization:
More informationStability radii and spectral value sets for real matrix perturbations Introduction D. Hinrichsen and B. Kelb Institut fur Dynamische Systeme Universitat Bremen D-859 Bremen, F. R. G. dh@mathematik.uni-bremen.de
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 informationA Worst-Case Estimate of Stability Probability for Engineering Drive
A Worst-Case Estimate of Stability Probability for Polynomials with Multilinear Uncertainty Structure Sheila R. Ross and B. Ross Barmish Department of Electrical and Computer Engineering University of
More informationRecitation 11: Time delays
Recitation : Time delays Emilio Frazzoli Laboratory for Information and Decision Systems Massachusetts Institute of Technology November, 00. Introduction and motivation. Delays are incurred when the controller
More informationLecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30
289 Upcoming labs: Lecture 12 Lab 20: Internal model control (finish up) Lab 22: Force or Torque control experiments [Integrative] (2-3 sessions) Final Exam on 12/21/2015 (Monday)10:30-12:30 Today: Recap
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
More informationLecture plan: Control Systems II, IDSC, 2017
Control Systems II MAVT, IDSC, Lecture 8 28/04/2017 G. Ducard Lecture plan: Control Systems II, IDSC, 2017 SISO Control Design 24.02 Lecture 1 Recalls, Introductory case study 03.03 Lecture 2 Cascaded
More informationThe Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals
The Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals Timo Hämäläinen Seppo Pohjolainen Tampere University of Technology Department of Mathematics
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More informationEECE 460 : Control System Design
EECE 460 : Control System Design SISO Pole Placement Guy A. Dumont UBC EECE January 2011 Guy A. Dumont (UBC EECE) EECE 460: Pole Placement January 2011 1 / 29 Contents 1 Preview 2 Polynomial Pole Placement
More informationA Comparative Study on Automatic Flight Control for small UAV
Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic
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 informationCDS 101/110a: Lecture 8-1 Frequency Domain Design
CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve
More informationCalculating determinants for larger matrices
Day 26 Calculating determinants for larger matrices We now proceed to define det A for n n matrices A As before, we are looking for a function of A that satisfies the product formula det(ab) = det A det
More informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition adings: Guzzella, Chapter 9.4 6 Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Emilio Frazzoli Institute
More informationControl System Design
ELEC ENG 4CL4: Control System Design Notes for Lecture #11 Wednesday, January 28, 2004 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Relative Stability: Stability
More informationData Based Design of 3Term Controllers. Data Based Design of 3 Term Controllers p. 1/10
Data Based Design of 3Term Controllers Data Based Design of 3 Term Controllers p. 1/10 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed
More informationAn Optimization-based Approach to Decentralized Assignability
2016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 2016 Boston, MA, USA An Optimization-based Approach to Decentralized Assignability Alborz Alavian and Michael Rotkowitz Abstract
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationRobust Control 9 Design Summary & Examples
Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via
More informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationFrequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability
Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods
More informationDidier HENRION henrion
GRADUATE COURSE ON POLYNOMIAL METHODS FOR ROBUST CONTROL PART II.1 ROBUST STABILITY ANALYSIS: POLYTOPIC UNCERTAINTY Didier HENRION www.laas.fr/ henrion henrion@laas.fr Cité de l espace with Ariane launcher
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationControl System Design
ELEC ENG 4CL4: Control System Design Notes for Lecture #14 Wednesday, February 5, 2003 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Chapter 7 Synthesis of SISO Controllers
More informationEconomics 472. Lecture 10. where we will refer to y t as a m-vector of endogenous variables, x t as a q-vector of exogenous variables,
University of Illinois Fall 998 Department of Economics Roger Koenker Economics 472 Lecture Introduction to Dynamic Simultaneous Equation Models In this lecture we will introduce some simple dynamic simultaneous
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 64: Dynamic Systems Spring Homework 4 Solutions Exercise 47 Given a complex square matrix A, the definition
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 informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationThe main aim of this exercise is to stress the structure of matrix arithmetic. a. If we had been given the equation
The main aim of this exercise is to stress the structure of matrix arithmetic. a. If we had been given the equation where XI A, B, and C were numbers, we would have at once concluded that X = 3 (AB + C).
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationRobust linear optimization under general norms
Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn
More informationUncertainty and Robustness for SISO Systems
Uncertainty and Robustness for SISO Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Outline Nature of uncertainty (models and signals). Physical sources of model uncertainty. Mathematical
More information12. Perturbed Matrices
MAT334 : Applied Linear Algebra Mike Newman, winter 208 2. Perturbed Matrices motivation We want to solve a system Ax = b in a context where A and b are not known exactly. There might be experimental errors,
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
More informationChapter 7 - Solved Problems
Chapter 7 - Solved Problems Solved Problem 7.1. A continuous time system has transfer function G o (s) given by G o (s) = B o(s) A o (s) = 2 (s 1)(s + 2) = 2 s 2 + s 2 (1) Find a controller of minimal
More informationSingular Value Decomposition Analysis
Singular Value Decomposition Analysis Singular Value Decomposition Analysis Introduction Introduce a linear algebra tool: singular values of a matrix Motivation Why do we need singular values in MIMO control
More informationTheory of Robust Control
Theory of Robust Control Carsten Scherer Mathematical Systems Theory Department of Mathematics University of Stuttgart Germany Contents 1 Introduction to Basic Concepts 6 1.1 Systems and Signals..............................
More information1 Matrices and Systems of Linear Equations
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 207, v 260) Contents Matrices and Systems of Linear Equations Systems of Linear Equations Elimination, Matrix Formulation
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework ypes of Uncertainty In this Lecture, we will cover Unstructured,
More informationROBUST CONTROL THROUGH ROBUSTNESS ENHANCEMENT
ROBUST CONTROL THROUGH ROBUSTNESS ENHANCEMENT Control Configurations And Two-Step Design Approaches Carles Pedret i Ferré Submitted in partial fulfillment of the requirements for the degree of Doctor Engineer
More informationLecture 1: Feedback Control Loop
Lecture : Feedback Control Loop Loop Transfer function The standard feedback control system structure is depicted in Figure. This represend(t) n(t) r(t) e(t) u(t) v(t) η(t) y(t) F (s) C(s) P (s) Figure
More informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
More informationLinear Algebra (part 1) : Matrices and Systems of Linear Equations (by Evan Dummit, 2016, v. 2.02)
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 206, v 202) Contents 2 Matrices and Systems of Linear Equations 2 Systems of Linear Equations 2 Elimination, Matrix Formulation
More informationECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27
1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system
More informationRICE UNIVERSITY. System Identication for Robust Control. Huipin Zhang. A Thesis Submitted. in Partial Fulfillment of the. Requirements for the Degree
RICE UNIVERSITY System Identication for Robust Control by Huipin Zhang A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Master of Science Approved, Thesis Committee: Athanasios
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 5: Matrix Perturbations Readings: DDV, Chapter 5 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology February 16, 2011 E. Frazzoli
More information16.30/31, Fall 2010 Recitation # 2
16.30/31, Fall 2010 Recitation # 2 September 22, 2010 In this recitation, we will consider two problems from Chapter 8 of the Van de Vegte book. R + - E G c (s) G(s) C Figure 1: The standard block diagram
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 informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More information