From the multivariable Nyquist plot, if 1/k < or 1/k > 5.36, or 0.19 < k < 2.7, then the closed loop system is stable.
|
|
- Florence Rice
- 6 years ago
- Views:
Transcription
1 1. (a) Multivariable Nyquist Plot (diagonal uncertainty case) dash=neg freq solid=pos freq Im Re From the multivariable Nyquist plot, if 1/k < -.37 or 1/k > 5.36, or.19 < k < 2.7, then the closed loop system is stable. (b)-(i) 6 singular value plot frequency (rad/sec) max singular value = bound on Delta =.1839 max (k 1, k 2 ) < 1/ =.184 (b)-(ii) bounds on (kl,k2): ++,+-,-+,-- ++: , +-: , -+: , --:
2 Plots of Hermitian part of diag(k 1 +/-,k 2 +/- ) G(jw). (k 1 >,k 2 >) (k 1 >,k 2 <) minimum eigenvalue of (G(jw)+G T (-jw))/2: k1>,k2> minimum eigenvalue of (G(jw)+G T (-jw))/2: k1>,k2< frequency (rad/sec) frequency (rad/sec) (k 1 <,k 2 >) (k 1 <,k 2 <) minimum eigenvalue of (G(jw)+G T (-jw))/2: k1<,k2> minimum eigenvalue of (G(jw)+G T (-jw))/2: k1<,k2< (b)-(iii) frequency (rad/sec) frequency (rad/sec) With k 1 in -.18,.18 and k 2 in -.18,.18, the following set of LMI s P (A+B diag(k 1 +,k 2 + )C) + (A+B diag(k 1 +,k 2 + )C) T P< P (A+B diag(k 1 +,k 2 - )C) + (A+B diag(k 1 +,k 2 - )C) T P< P (A+B diag(k 1 -,k 2 + )C) + (A+B diag(k 1 -,k 2 + )C) T P< P (A+B diag(k 1 -,k 2 - )C) + (A+B diag(k 1 -,k 2 - )C) T P< where (A,B,C,) is a realization of G(s), has a solution. Output of LMI solver: Iteration : Best value of t so far Result: best value of t: f-radius saturation:. of R = 1.e *** feasible! *** stability bound: k1: -.18,.18 k2: -.18,.18
3 (b)-(iv) ********* real mu *********... Mu upper bound (peak value): 5.372e+ Robust stability margin : 1.861e-1 ********* complex mu *********... Mu upper bound (peak value): 5.372e+ Robust stability margin : 1.861e-1 ********* Delta = unstructured (should be 1/H-infty norm of F) *********... Mu upper bound (peak value): 5.465e+ Robust stability margin : 1.83e-1 maximum real mu = maximum real bound for diag(k1,k2) = maximum complex mu = maximum complex bound for diag(k1,k2) = maximum bound for unstructured Delta = maximum spectral radius of F = maximum bound for k when Delta = k I (k complex) = (b)-(v) Closed loop characteristic polynomial: 2 s + (2 + k1 + 4 k2) s + 4 k2 + k k1 k2 Assume -k1max < k1 < k1max, -k2max < k2 < k2max, then the coefficients can be bounded by: 2 - k1max - 4 k2max < (2 + k1 + 4 k2) < 2 + k1max + 4 k2max - 4 k2max - k1max + 1 < 4 k2 + k1 + 1 < 4 k2max + k1max k1max k2max < 1-2 k1 k2 < 1 + k1max k2max all Kharitonov polynomials are stable when k1 <.18 k2 <.18 Solution Script (fin1.m) solution of problem 1 (Multivariable Final 22) clear all; close all; open loop plant
4 G=1 2;3 4*tf(1,1 1); aa,bb,cc,dd=ssdata(ss(g)); multivariable nyquist w=logspace(-4,4,5); sys=ss(aa,bb,cc,dd); sigmah1=mvarnyq(sys,w); diagonal uncertainty case k ; k -1/k<-.88-1/k> <k<11.4 figure(1); plot(real(sigmah1(1,:)),imag(sigmah1(1,:)),'-r',... real(sigmah1(1,:)),-imag(sigmah1(1,:)),'--r',... real(sigmah1(2,:)),imag(sigmah1(2,:)),'-b',... real(sigmah1(2,:)),-imag(sigmah1(2,:)),'--b'); title('multivariable Nyquist Plot (diagonal uncertainty case)',... ' dash=neg freq solid=pos freq') xlabel('re');ylabel('im');grid; axis equal small gain w=logspace(-1,2,2); sys=ss(aa,bb,cc,dd); sv=sigma(sys,w); figure(3);semilogx(w,sv);title('singular value plot'); xlabel('frequency (rad/sec)'); disp('max singular value = ',num2str(max(sv(1,:)))); disp('bound on Delta = ',num2str(1/max(sv(1,:)))); print -f3 -depsc smallgainexample.eps positive real realpartpp=pr(sys,w); figure(4);semilogx(w,realpartpp); title('minimum eigenvalue of (G(jw)+G^T(-jw))/2: k1>,k2>'); xlabel('frequency (rad/sec)'); realpartpn=pr(1 ; -1*sys,w); figure(5);semilogx(w,realpartpn); title('minimum eigenvalue of (G(jw)+G^T(-jw))/2: k1>,k2<'); xlabel('frequency (rad/sec)'); realpartnp=pr(-1 ; 1*sys,w); figure(6);semilogx(w,realpartnp); title('minimum eigenvalue of (G(jw)+G^T(-jw))/2: k1<,k2>'); xlabel('frequency (rad/sec)'); realpartnn=pr(-1 ; -1*sys,w); figure(7);semilogx(w,realpartnn); title('minimum eigenvalue of (G(jw)+G^T(-jw))/2: k1<,k2<'); xlabel('frequency (rad/sec)');
5 print -f4 -depsc prexample1.eps print -f5 -depsc prexample2.eps print -f6 -depsc prexample3.eps print -f7 -depsc prexample4.eps bpp=min(realpartpp); bpn=min(realpartpn); bnp=min(realpartnp); bnn=min(realpartnn); disp('bounds on (kl,k2): ++,+-,-+,--'); disp(sprintf('++:.5g, +-:.5g, -+:.5g, --:.5g\n',... 1/bpp,1/bpn,1/bnp,1/bnn)); LMI n=size(aa,1); L1=1 ; ; L2= ; 1; a1=-bb*l1*cc; a2=-bb*l2*cc; define lmi variable to be solved setlmis(); P,ndec,Pdec=lmivar(1,n 1); set up uncertain variable range k1=-1,1*.18 range of k1 k2=-1,1*.18; range of k2 k1a=k1(1);k1b=k1(2); k2a=k2(1);k2b=k2(2); set up LMI to be solved left hand side of first lmi: m1 1,1th block of the lmi: l1,l2 first argument = m1,l1,l2,p second/third a,b 's' to specify a*p*b+b'*p*a' lmiterm(1,1,1,p,(aa+k1a*a1+k2a*a2)',1,'s'); left hand side of second lmi: m1 1,1th block of the lmi: l1,l2 a*p*b first argument = m1,l1,l2,p second/third a,b lmiterm(2,1,1,p,(aa+k1b*a1+k2a*a2)',1,'s'); left hand side of third lmi: m1 1,1th block of the lmi: l1,l2 a*p*b
6 first argument = m1,l1,l2,p second/third a,b lmiterm(3,1,1,p,(aa+k1a*a1+k2b*a2)',1,'s'); left hand side of fourth lmi: m1 1,1th block of the lmi: l1,l2 a*p*b first argument = m1,l1,l2,p second/third a,b lmiterm(4,1,1,p,(aa+k1b*a1+k2b*a2)',1,'s'); example=getlmis; options=,1,-1,-1,; tmin,pfeas=feasp(example,options); disp(tmin) if tmin<; disp('*** feasible! ***'); disp('stability bound: '); disp('k1: ',num2str(k1a),',',num2str(k1b)); disp('k2: ',num2str(k2a),',',num2str(k2b)); else disp('*** no feasible solution (reduce uncertainty bound)! ***'); end mu sys1=ltisys(aa,bb,cc,dd,1); delta1=ublock(1,1,'ltisr'); delta2=ublock(1,1,'ltisr'); deltar=udiag(delta1,delta2); disp('********* real mu ********* '); marginr,peakf,fs,ds,gs = mustab(sys1,deltar); delta1=ublock(1,1,'ltisc'); delta2=ublock(1,1,'ltisc'); deltac=udiag(delta1,delta2); disp('********* complex mu ********* '); marginc,peakf,fs,ds,gs = mustab(sys1,deltac); deltaf=ublock(2,1,'ltifc'); disp('********* Delta = unstructured ',... '(should be 1/H-infty norm of F) ********* '); marginf,peakf,fs,ds,gs = mustab(sys1,deltaf); should be same as H- inf disp('maximum real mu = ',num2str(1/marginr)); disp('maximum real bound for diag(k1,k2) = ',num2str(marginr)); disp('maximum complex mu = ',num2str(1/marginc)); disp('maximum complex bound for diag(k1,k2) = ',num2str(marginc)); disp('maximum bound for unstructured Delta = ',num2str(marginc)); spectral radius rho_c=(max(abs(sigmah1))); disp('maximum spectral radius of F = ',... num2str(max(rho_c))); disp('maximum bound for k when Delta = k I (k complex) = ',... num2str(1/max(rho_c)));
7 Kharitonov syms s k1 k2 k1mx=.18;k2mx=.18; kk=k1 ; k2; charpoly=collect(simplify(det(s*eye(n,n)-(aa-bb*kk*cc))),s); pretty(charpoly); a1u=2+k1mx+4*k2mx; a1l=2-k1mx-4*k2mx; au=4*k2mx + k1mx *k1mx *k2mx; al=-4*k2mx - k1mx *k1mx *k2mx; poly1=1 a1l au; poly2=1 a1u al; poly3=1 a1l al; poly4=1 a1u au; r1=max(real(roots(poly1))); r2=max(real(roots(poly2))); r3=max(real(roots(poly3))); r4=max(real(roots(poly4))); if r1<&r2<&r3<&r4< disp('all Kharitonov polynomials are stable'); disp(' k1 <',num2str(k1mx),' ',' k2 <',num2str(k2mx)); else disp('not all Kharitonov polynomials are stable, reduce bounds!'); end
8 Problem 2. Given the state space representation of the plant Parameterize an n K -th order controller as ẋ = Ax + Bu y = Cx. ẋ K = A K x K + B K y = A K x K + B K Cx u = C K x K + D K y = C K x K + D K Cx. Substituting back into the plant state equation we have ẋ =(A + BD K C)x + BC K x K. Putting x and x K together, the closed loop system is described by ẋ A + BDK C BC = K x ẋ K B K C A K x K = A B AK B + K I x I C K D K C }{{}}{{}}{{}}{{} J A The closed loop system is stable if and only if there exists X cl > such that B C x K. X cl (A + BJC)+(A + BJC) T X cl <. (1) Define H = X cl A + A T X cl Q = B T X cl P = C. Then (1) becomes H + Q T JP + P T J T Q. (2) Note that with loss of generality we can assume B and C full rank, then Q and P are also full rank. Now apply Lemma 7.2, then we have (2) equivalent to N T QHN Q <, N T P HN P < (3) where R(N Q )=N (Q) and R(N P )=N (P ). In our case, N Q = X 1 cl N T c N P = No where R(N o )=N (C) and R(N T c )=N (B T ). 1
9 Let and write X X2 X cl = X2 T X 3 X 1 cl = Y Y2 Y2 T Y 3 By Lemma 7.9, X cl > if and only if X I X I and rank I Y I Y lasubsituting (4) (5) into (3) and simplify, we obtain Together with (6), this proves Theorem 5.8. No T (A T X + XA)N o < N c (AY + YA T )N c <. (4). (5) n + n K. (6) 2
10 (a) Let the input of G be denoted z and output w. Then from the block diagram we have z = K(Gz + w) = (I + KG) 1 Kw = K(I + GK) 1 w. Given a general interconnection between a system P (with input (w, u) and output (z,y)) and K (feedback from y to u), the closed loop transfer function is T zw = P 11 P 12 K(I + P 22 K) 1 P 21. Matching terms, it s clear that P 11 =, P 12 = I, P 21 = I, P 22 = G,or I P =. I G From the robust stability theorem, the closed loop system is stable for all G ν if and only if T zw H < 1/ν. (b) The condition T zw H < 1/ν is equivalent to T νz,w H < 1. A state space realization for T νz,w is ẋ = Ax + w + Bu = Ax + w +(B/ν)νu νz = x + w + I νu y = Cx + I w + νu. The H optimization requires the solution of X and Y.ForX, we have the Riccati Equation: A T X + X A + C T 1 C 1 + X (B 1 B T 1 B 2 B T 2 )X =. Since B 1 =, C 1 =, B 2 = B/ν,wehave Similarly, for Y (note that C 2 = C): A T X + X A ν 2 X BB T X =. = AY + Y A T + B1 T B 1 + Y (C1 T C 1 C2 T C 2 )Y = AY + Y A T Y C T CY. Let W c be the controllability grammian of ( A +,B + ) and W o be the observability grammian of (C +, A + ): W c ( A + ) T +( A + )W c + B+B T + = W o ( A + )+( A + ) T W o + C + C+ T =. Then it follows that X and Y given below satisfy the Riccati equations: W 1 c ν X = 2 W 1 o, Y =. 1
11 The condition ρ(x Y ) < 1 (ρ is the spectral radius) is equivalent to ρ(wc 1 This is further equivalent to ρ(w c W o ) >ν 2, which is the stated condition. The H optimal controller is of the form: ˆx = (A B 2 B T 2 X L C)ˆx + L y = (A B + B T +W 1 c L C)ˆx + L y L = (I Y X ) 1 Y C T =(I W 1 o Wc 1 ) 1 Wo 1 C+ T νu = B2 T X ˆx = BT + ν W c 1 ν 2ˆx = B+W T c 1 νˆx. The actual control signal is then u = B+W T c 1 ˆx. (c) Given g(s) =1/(s 1) 2. Choose, for example, the realization 1 1 A =, B =, C = 1, D =. 1 1 Then W c = , W o = W 1 o )) < 1/nu 2. The supremum of ν is.136. For the controller design, choose something slightly smaller, e.g., at 95 value, ν =.984. The H controller is then s.423 K(s) = (s )(s +4.76). The closed loop K(I + GK) 1 H is 1.2 which is 1/ν. The closed loop poles are ( 33.2, 1, 1,.42). MATLAB Code fin3c.m: solution of problem 3(c) (Multivariable Final 22) A=1 1; 1;B=;1;C=1 ;D=; sys=ss(-a,b,c,d); Wc=gram(sys, c ); Wo=gram(sys, o ); numax=sqrt(min(abs(eig(wc*wo)))); nu=.95*numax; disp(sprintf( numax = 5.3g,numax)); disp(sprintf( Choose nu = 5.3g,nu)); disp(sprintf( Theoretical achivable H-infinity norm = 1/nu = 5.3g,1/nu)); B2=B/nu; 2
12 Xinf=inv(Wc)*nuˆ2; Yinf=inv(Wo); Linf=inv(eye(2,2)-Yinf*Xinf)*Yinf*C ; AK=A-B2*B2 *Xinf-Linf*C; BK=Linf; CK=B2 *Xinf; DK=; K=ss(AK,BK,CK,DK)/nu; G=ss(A,B,C,D); Gcl=feedback(-K,G,1); a,b,c,d=ssdata(gcl); sysgcl=pck(a,b,c,d); out=hinfnorm(sysgcl); disp(sprintf( closed loop H-infinity norm in 5.3g,5.3g,... out(1),out(2))); disp(sprintf( tolerable additive uncertainty bound: 4.3g,4.3g,... 1/out(1),1/out(2))); figure(1);bode(gcl) disp( closed loop poles: ); disp(pole(gcl)); 3
13 (a) Let k(s) be stabilizing and k = y/x a stable coprime factorization such that xm yn =1. If k(s) is stable, then x(s) cannot have zeros in the open right half plane. Since x(s) varies continuously in s, x(z 1 ) and x(z 2 ) must have the same sign for any real and non-negative z 1 and z 2 (this is necessary but not sufficient). Let z 1 and z 2 be two consecutive open RHP zeros of g(s) (i.e., n(z 1 )=n(z 2 )=). Then from Bezout Identity, we have x(z 1 )m(z 1 )=x(z 2 )m(z 2 ). Since x(z 1 ) and x(z 2 ) have the same signs, m(z 1 ) and m(z 2 ) must also have the same signs. Each real pole of g(s) (i.e., zeros of m(s)) between z 1 and z 2 will cause a sign change in m(s) as s varies from z 1 to z 2. Therefore, there must be an even number of poles between z 1 and z 2 in order for m(z 1 ) and m(z 2 ) to have the same signs. (b) Consider g(s) =s/(s 2 1). Between the zeros at and +, there is one real pole, therefore, g does not satisfy PIP. Let n = s/(s +2) 2, m =(s 2 1)/(s +2) 2. Use the nominal controller k(s) = (s +4)/(s.5), then we can choose y = (s +4)/(s +2)and x =(s.5)/(s +2). Now calculate We can set δ(s) =m(s)x(s) n(s)y(s) = 1 2s 3 + s 2 +6s (s +2) 3 2(s +4)(s +2)2 y 1 = y/δ = 2s 3 + s 2 +6s +1 (2s 1)(s +2)2 x 1 = x/δ = 2s 3 + s 2 +6s +1 then (y 1,x 1 ) would satisfy the Bezout Identity. The complete parameterization of stabilizing controllers is given by k = y 1 mq x 1 nq where q is any RH function. The stability of k is determined by the zeros of d k = x 1 nq. Evaluate d k at and, we get d k () = 4 and d k ( ) =1. Since they are of different signs, d k (s) must cross zero for some positive and real s. 1
From the multivariable Nyquist plot, if 1/k < or 1/k > 5.36, or 0.19 < k < 2.7, then the closed loop system is stable.
1. (a) Multivariable Nyquist Plot (diagonal uncertainty case) dash=neg freq solid=pos freq 2.5 2 1.5 1.5 Im -.5.5-2 -2.5 1 2 3 4 5 Re From the multivariable Nyquist plot, if 1/k < -.37 or 1/k > 5.36, or.19
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
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 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 informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
More information8 A First Glimpse on Design with LMIs
8 A First Glimpse on Design with LMIs 8.1 Conceptual Design Problem Given a linear time invariant system design a linear time invariant controller or filter so as to guarantee some closed loop indices
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 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 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 informationToday (10/23/01) Today. Reading Assignment: 6.3. Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10
Today Today (10/23/01) Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10 Reading Assignment: 6.3 Last Time In the last lecture, we discussed control design through shaping of the loop gain GK:
More informationRobust Multivariable Control
Lecture 1 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Addresses email: anders.helmerson@liu.se mobile: 0734278419 http://users.isy.liu.se/rt/andersh/teaching/robkurs.html
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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationRobust Control, H, ν and Glover-McFarlane
Robust Control, H, ν and Glover-McFarlane Department of Automatic Control LTH, Lund University Robust Control 1 MIMO performance 2 Robustness and the H -norm 3 H -control 4 ν-gap metric 5 Glover-MacFarlane
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 informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 21: Optimal Output Feedback Control connection is called the (lower) star-product of P and Optimal Output Feedback ansformation (LFT).
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationProblem Set 4 Solution 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 4 Solution Problem 4. For the SISO feedback
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
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 informationTopic # Feedback Control
Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationLecture 15: H Control Synthesis
c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control
More informationHuman Friendly Control : an application to Drive by Wire
Human Friendly Control : an application to Drive by Wire By Emmanuel Witrant Master Thesis, 2001-2002 Laboratoire d Automatique de Grenoble with Carlos Canudas-de-Wit 2 Juillet 2002 Typeset by FoilTEX
More informationStability and Robustness 1
Lecture 2 Stability and Robustness This lecture discusses the role of stability in feedback design. The emphasis is notonyes/notestsforstability,butratheronhowtomeasurethedistanceto instability. The small
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 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 informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
More informationChapter 8 SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH
Chapter 8 SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH 8.1 INTRODUCTION Practical systems are difficult to model, as this
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 informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
More informationStability of CL System
Stability of CL System Consider an open loop stable system that becomes unstable with large gain: At the point of instability, K( j) G( j) = 1 0dB K( j) G( j) K( j) G( j) K( j) G( j) =± 180 o 180 o Closed
More informationStationary trajectories, singular Hamiltonian systems and ill-posed Interconnection
Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute
More informationLecture 10: Linear Matrix Inequalities Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Matrix Inequalities A linear matrix inequality (LMI)
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 informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
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 informationRobustness Analysis and Controller Synthesis with Non-Normalized Coprime Factor Uncertainty Characterisation
211 5th IEEE onference on Decision and ontrol and European ontrol onference (D-E) Orlando, FL, USA, December 12-15, 211 Robustness Analysis and ontroller Synthesis with Non-Normalized oprime Factor Uncertainty
More informationProblem Set 8 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6245: MULTIVARIABLE CONTROL SYSTEMS by A Megretski Problem Set 8 Solutions 1 The problem set deals with H-Infinity
More informationLecture 2. FRTN10 Multivariable Control. Automatic Control LTH, 2018
Lecture 2 FRTN10 Multivariable Control Course Outline L1 L5 Specifications, models and loop-shaping by hand 1 Introduction 2 Stability and robustness 3 Specifications and disturbance models 4 Control synthesis
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More information3 Stabilization of MIMO Feedback Systems
3 Stabilization of MIMO Feedback Systems 3.1 Notation The sets R and S are as before. We will use the notation M (R) to denote the set of matrices with elements in R. The dimensions are not explicitly
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationLinear Matrix Inequality Approach
1 MEM800-007 Chapter 4a Linear Matrix Inequality Approach Reference: Linear matrix inequalities in system and control theory, Stephen Boyd et al. SIAM 1994. Linear Matrix Inequality (LMI) approach have
More informationRobust Controller Design By Convex Optimization. Alireza Karimi Laboratoire d Automatique, EPFL
Robust Controller Design By Convex Optimization Alireza Karimi Laboratoire d Automatique, EPFL References 1. Linear Matrix Inequalities in System and Control Theory, S. Boyd, L. El Ghaoui, E. Feron and
More informationChapter 3. State Feedback - Pole Placement. Motivation
Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based
More informationEEE 184: Introduction to feedback systems
EEE 84: Introduction to feedback systems Summary 6 8 8 x 7 7 6 Level() 6 5 4 4 5 5 time(s) 4 6 8 Time (seconds) Fig.. Illustration of BIBO stability: stable system (the input is a unit step) Fig.. step)
More informationHomework Assignment 3
ECE382/ME482 Fall 2008 Homework 3 Solution October 20, 2008 1 Homework Assignment 3 Assigned September 30, 2008. Due in lecture October 7, 2008. Note that you must include all of your work to obtain full
More informationChapter 5. Standard LTI Feedback Optimization Setup. 5.1 The Canonical Setup
Chapter 5 Standard LTI Feedback Optimization Setup Efficient LTI feedback optimization algorithms comprise a major component of modern feedback design approach: application problems involving complex models
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 informationOutline. Classical Control. Lecture 1
Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction
More informationLecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 3 Chapter 4: Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 3 p. 1/77 3.1 System Descriptions [4.1] Let f(u) be a liner operator, u 1 and u
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
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 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 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 informationRobust control of MIMO systems
Robust control of MIMO systems Olivier Sename Grenoble INP / GIPSA-lab September 2017 Olivier Sename (Grenoble INP / GIPSA-lab) Robust control September 2017 1 / 162 O. Sename [GIPSA-lab] 2/162 1. Some
More information9 The LQR Problem Revisited
9 he LQR Problem Revisited Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := lim E [ z(t) z(t) ] t for the LI system Assumptions: a) D
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationOutline. Linear Matrix Inequalities in Control. Outline. System Interconnection. j _jst. ]Bt Bjj. Generalized plant framework
Outline Linear Matrix Inequalities in Control Carsten Scherer and Siep Weiland 7th Elgersburg School on Mathematical Systems heory Class 3 1 Single-Objective Synthesis Setup State-Feedback Output-Feedback
More informationS I X SOLUTIONS TO CASE STUDIES CHALLENGES. Antenna Control: Stability Design via Gain K s s s+76.39K. T(s) =
S I X Stability SOLUTIONS TO CASE STUDIES CHALLENGES Antenna Control: Stability Design via Gain From the antenna control challenge of Chapter 5, Make a Routh table: 76.39K s 3 +151.32s 2 +198s+76.39K s
More informationK(s +2) s +20 K (s + 10)(s +1) 2. (c) KG(s) = K(s + 10)(s +1) (s + 100)(s +5) 3. Solution : (a) KG(s) = s +20 = K s s
321 16. Determine the range of K for which each of the following systems is stable by making a Bode plot for K = 1 and imagining the magnitude plot sliding up or down until instability results. Verify
More informationClosed-Loop Structure of Discrete Time H Controller
Closed-Loop Structure of Discrete Time H Controller Waree Kongprawechnon 1,Shun Ushida 2, Hidenori Kimura 2 Abstract This paper is concerned with the investigation of the closed-loop structure of a discrete
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
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 informationEntering Matrices
MATLAB Basic Entering Matrices Transpose Subscripts The Colon Operator Operators Generating Matrices Element-by-element operation Creating a Plot Multiple Data Sets in One Graph Line Styles and Colors
More informationH 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process
More informationMEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation
More informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
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.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationAnti-windup synthesis using Riccati equations
Anti-windup synthesis using Riccati equations Jorge Sofrony, Matthew C. Turner and Ian Postlethwaite Control and Instrumentation Research Group, Department of Engineering, University of Leicester, Leicester,
More informationTopic # Feedback Control Systems
Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control
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 informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) D. Linear Matrix Inequality D.1 Convex Optimization D.2 Linear Matrix Inequality(LMI) D.3 Control Design and LMI Formulation
More informationStructured Uncertainty. Block diagonal
Structured Uncertainty Block diagonal Example 5.8 0.1 δ 0.1, 1 δ 1 1 2 error : N yw ( s) Example 5.9 The Structured Singular Value and Stability Robustness Assume nominal closed-loop system N(s) is
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 informationAn Internal Stability Example
An Internal Stability Example Roy Smith 26 April 2015 To illustrate the concept of internal stability we will look at an example where there are several pole-zero cancellations between the controller and
More informationEngraving Machine Example
Engraving Machine Example MCE44 - Fall 8 Dr. Richter November 24, 28 Basic Design The X-axis of the engraving machine has the transfer function G(s) = s(s + )(s + 2) In this basic example, we use a proportional
More informationA Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method
A Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method Kin Cheong Sou and Anders Rantzer Abstract In this paper a controller reduction method which preserves
More informationIMPACT OF PLANT COUPLING ON DECENTRALIZED PI CONTROL LAW. Date: 08/30/2016
IMPACT OF PLANT COUPLING ON DECENTRALIZED PI CONTROL LAW Date: 8/3/216 1 / 37 Standard Negative Feedback System With Pre-Filter d i input disturbance d o output disturbance r reference command W Pre-Filter
More informationTopic # Feedback Control Systems
Topic #19 16.31 Feedback Control Systems Stengel Chapter 6 Question: how well do the large gain and phase margins discussed for LQR map over to DOFB using LQR and LQE (called LQG)? Fall 2010 16.30/31 19
More informationOQ4867. Let ABC be a triangle and AA 1 BB 1 CC 1 = {M} where A 1 BC, B 1 CA, C 1 AB. Determine all points M for which ana 1...
764 Octogon Mathematical Magazine, Vol. 24, No.2, October 206 Open questions OQ4867. Let ABC be a triangle and AA BB CC = {M} where A BC, B CA, C AB. Determine all points M for which 4 s 2 3r 2 2Rr AA
More informationRobust control for a multi-stage evaporation plant in the presence of uncertainties
Preprint 11th IFAC Symposium on Dynamics and Control of Process Systems including Biosystems June 6-8 16. NTNU Trondheim Norway Robust control for a multi-stage evaporation plant in the presence of uncertainties
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More information9. Two-Degrees-of-Freedom Design
9. Two-Degrees-of-Freedom Design In some feedback schemes we have additional degrees-offreedom outside the feedback path. For example, feed forwarding known disturbance signals or reference signals. In
More informationControl Systems. Design of State Feedback Control.
Control Systems Design of State Feedback Control chibum@seoultech.ac.kr Outline Design of State feedback control Dominant pole design Symmetric root locus (linear quadratic regulation) 2 Selection of closed-loop
More informationSYNTHESIS OF LOW ORDER MULTI-OBJECTIVE CONTROLLERS FOR A VSC HVDC TERMINAL USING LMIs
SYNTHESIS OF LOW ORDER MULTI-OBJECTIVE CONTROLLERS FOR A VSC HVDC TERMINAL USING LMIs Martyn Durrant, Herbert Werner, Keith Abbott Control Institute, TUHH, Hamburg Germany; m.durrant@tu-harburg.de; Fax:
More informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
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 informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
More information