Topic # Feedback Control Systems
|
|
- Audra Carr
- 5 years ago
- Views:
Transcription
1 Topic # Feedback Control Systems State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality
2 Fall Combined Estimators and Regulators Now evaluate stability and/or performance of a controller when K designed assuming u(t) = Kx(t) but implemented as u(t) = K ˆx(t) Assume we have designed a closed-loop estimator with gain L ˆx(t) = (A LC)ˆx(t) + Bu(t) + Ly(t) ŷ(t) = C ˆx(t) The closed-loop system dynamics are given by: ẋ(t) = Ax(t) + Bu(t) ˆx(t) = (A LC)ˆx(t) + Bu(t) + Ly(t) y(t) = Cx(t) ŷ(t) = C ˆx(t) u(t) = K ˆx(t) Which can be compactly written as: [ ] [ ] [ ] ẋ(t) A BK x(t) = ˆx(t) LC A BK LC ˆx(t) ẋ cl (t) = A cl x cl (t) This does not look too good at this point not even obvious that the closed-system is stable. λ i (A cl ) =??
3 Fall Can fix this problem by introducing a new variable x = x ˆx and then converting the closed-loop system dynamics using the similarity transformation T x cl (t) [ x(t) x(t) Note that T = T 1 ] = [ I 0 I I ] [ x(t) ˆx(t) ] = T x cl (t) Now rewrite the system dynamics in terms of the state x cl (t) A cl T A cl T 1 Ācl Since similarity transformations preserve the eigenvalues we are guaranteed that λ i (A cl ) λ i (Ācl) Work through the math: [ ] [ ] [ I 0 A BK I 0 Ā cl = I I LC A BK LC I I [ ] [ ] A BK I 0 = A LC A + LC I I [ ] A BK BK = 0 A LC ] Because Ācl is block upper triangular, we know that the closedloop poles of the system are given by det(si Ā cl ) det(si (A BK)) det(si (A LC)) = 0
4 Fall Observation: Closed-loop poles for this system consist of the union of the regulator poles and estimator poles. So we can just design the estimator/regulator separately and combine them at the end. Called the Separation Principle. Keep in mind that the pole locations you are picking for these two sub-problems will also be the closed-loop pole locations. Note: Separation principle means that there will be no ambiguity or uncertainty about the stability and/or performance of the closedloop system. The closed-loop poles will be exactly where you put them!! And we have not even said what compensator does this amazing accomplishment!!!
5 Fall The Compensator Dynamic Output Feedback Compensator is the combination of the regulator and estimator using u(t) = K ˆx(t) ˆx(t) = (A LC)ˆx(t) + Bu(t) + Ly(t) ˆx(t) = (A BK LC)ˆx(t) + Ly(t) u(t) = K ˆx(t) Rewrite with new state x c (t) ˆx(t) ẋ c (t) = A c x c (t) + B c y(t) u(t) = C c x c (t) where the compensator dynamics are given by: A c A BK LC, B c L, C c K Note that the compensator maps sensor measurements to actuator commands, as expected. Closed-loop system stable if regulator/estimator poles placed in the LHP, but compensator dynamics do not need to be stable. λ i (A BK LC) =??
6 Fall For consistency with the implementation of classical controllers, define the compensator transfer function as G c (s) = C c (si A c ) 1 B c so that u(t) = G c (s)y(t) Note that it is often very easy to provide classical interpretations (such as lead/lag) for the compensator G c (s). Can then implement this compensator with a reference command r(t) using feedback on e(t) = r(t) y(t) rather than just y(t) r(t) e(t) u(t) y(t) G c (s) G(s) u(t) = G c (s)e(t) = G c (s)(r(t) y(t)) Still have u(t) = G c (s)y(t) if r(t) = 0. So now write state space model of the compensator as ẋ c (t) = A c x c (t) + B c e(t) u(t) = C c x c (t) Intuitively appealing setup because it is the same approach used for the classical control, but it turns out not to be the best approach More on this later.
7 Fall Mechanics Basics: e(t) = r(t) y(t), u(t) = G c e(t), y(t) = Gu(t) G c (s) : ẋ c (t) = A c x c (t) + B c e(t), u(t) = C c x c (t) G(s) : ẋ(t) = Ax(t) + Bu(t), y(t) = Cx(t) Loop dynamics L(s) = G(s)G c (s) y(t) = L(s)e(t) ẋ(t) = Ax(t) +BC c x c (t) ẋ c (t) = +A c x c (t) +B c e(t) L(s) : [ ẋ(t) ẋ c (t) ] = [ A BCc 0 A c ] [ x(t) x c (t) ] y(t) = [ C 0 ] [ x(t) x c (t) ] [ 0 + B c ] e(t) To close the loop, note that e(t) = r(t) y(t), then [ ] [ ] [ ] [ ] ( ẋ(t) A BCc x(t) 0 = + r(t) [ C 0 ] [ x(t) ẋ c (t) 0 A c x c (t) B c x c (t) [ ] [ ] [ ] A BC c x(t) 0 = + r(t) B c C A c x c (t) B c y(t) = [ C 0 ] [ ] x(t) x c (t) A cl is not exactly the same as on page 16-1 because we have re-arranged where the negative sign enters into the problem. Same result though. ])
8 Fall Simple Compensator Example Let G(s) = 1/(s 2 + s + 1) with state-space model given by: [ ] [ ] A =, B =, C = [ 1 0 ], D = Design the regulator to place the poles at s = 4 ± 4j λ i (A BK) = 4 ± 4j K = [ 31 7 ] Regulator pole time constant τ c = 1/ζω n 1/4 = 0.25 sec Put estimator poles with faster time constant τ e 1/10 Use real poles, so Φ e (s) = (s + 10) 2 [ ] 1 [ ] C 0 L = Φ e (A) CA 1 [ ] 2 [ ] [ ] [ = [ ] [ ] [ ] = = ] 1 [ 0 1 ]
9 Fall Compensator: A c = A [ BK ] LC [ = [ ] 19 1 = B c = L = [ ] C c = K = [ 31 7 ] ] [ 31 7 ] [ ] [ 1 0 ] Compensator transfer function: G c (s) = C c (si A c ) 1 B c U E s = 1149 s s Note that the compensator has a low frequency real zero and two higher frequency poles. Thus it looks like a lead compensator.
10 Fall Figure 1: Plant is pretty simple and the compensator looks like a lead 2 10 rads/sec.
11 Fall Figure 2: Loop transfer function L(s) shows the slope change near ω c = 5 rad/sec. Note that we have a large PM and GM.
12 Fall Figure 3: Freeze compensator poles and zeros; look at root locus of closed-loop poles versus an additional loop gain α (nominally α = 1.) Figure 4: Closed-loop transfer function.
13 Fall Dynamic Output Feedback Example (reg est.m) 1 % Combined estimator/regulator design for a simple system 2 % G= 1/(s^2+s+1) 3 % 4 % Jonathan How 5 % Fall % 7 close all;clear all 8 for ii=1:5 9 figure(ii);clf;set(gcf, DefaultLineLineWidth,2);set(gcf, DefaultlineMarkerSize,10) 10 end a=[0 1;-1-1];b=[0 1] ;c=[1 0];d=0; 13 k=acker(a,b,[-4+4*j;-4-4*j]); 14 l=acker(a,c,[-10-10] ) ; 15 % 16 % For state space for G_c(s) 17 % 18 ac=a-b*k-l*c;bc=l;cc=k;dc=0; G=ss(a,b,c,d); 21 Gc=ss(ac,bc,cc,dc); f=logspace(-1,2,400); 24 g=freqresp(g,f*j);g=squeeze(g); 25 gc=freqresp(gc,f*j);gc=squeeze(gc); figure(1);clf 28 subplot(211) 29 loglog(f,abs(g),f,abs(gc), -- );axis([.1 1e2.02 1e2]) 30 xlabel( Freq (rad/sec) );ylabel( Mag ) 31 %legend( Plant G, Compensator Gc ); 32 grid 33 subplot(212) 34 semilogx(f,180/pi*angle(g),f,180/pi*angle(gc), -- ); 35 axis([.1 1e ]) 36 xlabel( Freq (rad/sec) );ylabel( Phase (deg) );grid 37 legend( Plant G, Compensator Gc ) L=g.*gc; figure(2);clf 42 subplot(211) 43 loglog(f,abs(l),[.1 1e2],[1 1]);axis([.1 1e2.2 1e2]) 44 xlabel( Freq (rad/sec) );ylabel( Mag ) 45 legend( Loop L ); 46 grid 47 subplot(212) 48 semilogx(f,180/pi*phase(l. ),[.1 1e2],-180*[1 1]); 49 axis([.1 1e ]) 50 xlabel( Freq (rad/sec) );ylabel( Phase (deg) );grid 51 % 52 % loop dynamics L = G Gc 53 % 54 al=[a b*cc;zeros(2) ac]; 55 bl=[zeros(2,1);bc]; 56 cl=[c zeros(1,2)]; 57 dl=0; 58 figure(3) 59 rlocus(al,bl,cl,dl) 60 % 61 % closed-loop dynamics 62 % unity gain wrapped around loop L 63 % 64 acl=al-bl*cl;bcl=bl;ccl=cl;dcl=d; 65 % 66 % scale closed-loop to get zero SS error 67 N=inv(ccl*inv(-acl)*bcl) hold on;plot(eig(acl), d, MarkerFaceColor, b );hold off 70 grid on
14 Fall axis([ ]) 72 % 73 % closed-loop freq response 74 % 75 Gcl=ss(acl,bcl*N,ccl,dcl*N); 76 gcl=freqresp(gcl,f*j);gcl=squeeze(gcl); figure(4);clf 79 loglog(f,abs(g),f,abs(gcl), -- ); 80 axis([.1 1e2.01 1e2]) 81 xlabel( Freq (rad/sec) );ylabel( Mag ) 82 legend( Plant G, closed-loop G_{cl} );grid 83 title([ Factor of N=,num2str(N), applied to Closed loop ]) figure(5);clf 86 margin(al,bl,cl,dl) figure(1);orient tall;print -dpng -r300 reg_est1 89 figure(2);orient tall;print -dpng -r300 reg_est2 90 figure(3);print -dpng -r300 reg_est3 91 figure(4);print -dpng -r300 reg_est4 92 figure(5);print -dpng -r300 reg_est5
15 Fall Figure 5: Example #1: G(s) = (s+8)(s+14)(s+20)
16 Fall Figure 6: Example #1: G(s) = (s+8)(s+14)(s+20) closed-loop poles, open-loop poles, Compensator poles, Compensator zeros Two compensator zeros at ±6.63j draw the two lower frequency plant poles further into the LHP. Compensator poles are at much higher frequency. Looks like a lead compensator.
17 Fall Figure 7: Example #2: G(s) = 0.94 s
18 Fall Figure 8: Example #2: G(s) = 0.94 s closed-loop poles, open-loop poles, Compensator poles, Compensator zeros Compensator zero at draws the two lower frequency plant poles further into the LHP. Compensator poles are at much higher frequency. Looks like a lead compensator.
19 Fall Figure 9: Example #3: G(s) = (s 8)(s 14)(s 20)
20 Fall Figure 10: Example #3: G(s) = (s 8)(s 14)(s 20) closed-loop poles, open-loop poles, Compensator poles, Compensator zeros Compensator zeros at 3.72±8.03j draw the two higher frequency plant poles further into the LHP. Lowest frequency one heads into the LHP on its own. Compensator poles are at much higher frequency. Note sure what this looks like.
21 Fall Figure 11: Example #4: G(s) = (s 1) (s+1)(s 3)
22 Fall Figure 12: Example #4: G(s) = (s 1) (s+1)(s 3) closed-loop poles, open-loop poles, Compensator poles, Compensator zeros Compensator zero at -1 cancels the plant pole. Note the very unstable compensator pole at s = 9!! Needed to get the RHP plant pole to branch off the real line and head into the LHP. Other compensator pole is at much higher frequency. Note sure what this looks like.
23 Fall Separation principle gives a very powerful and simple way to develop a dynamic output feedback controller Note that the designer now focuses on selecting the appropriate regulator and estimator pole locations. Once those are set, the closed-loop response is specified. Can almost consider the compensator to be a by-product. These examples show that the design process is extremely simple.
24 Fall Reference Input - II On page 16-5, compensator implemented with a reference command by changing to feedback on e(t) = r(t) y(t) rather than y(t) r(t) e(t) u(t) y(t) G c (s) G(s) So u(t) = G c (s)e(t) = G c (s)(r(t) y(t)) Intuitively appealing because same approach used for classical control, but it turns out not to be the best. Can improve the implementation by using a more general form: ẋ c (t) = A c x c (t) + B c y(t) + Gr(t) u(t) = C c x c (t) + Nr(t) Now explicitly have two inputs to the controller (y(t) and r(t)) N performs the same role that we used it for previously. Introduce G as an extra degree of freedom in the problem. Turns out that setting G = B N is a particularly good choice. Following presents some observations on the impact of G
25 Fall First: this generalization does not change the closed-loop poles of the system, regardless of the selection of G and N, since ẋ(t) = Ax(t) + Bu(t), y(t) = Cx(t) ẋ c (t) = A c x c (t) + B c y(t) + Gr(t) [ ẋ(t) ẋ c (t) u(t) = C c x c (t) + Nr(t) ] = [ ] [ ] A BCc x(t) + B c C A c x c (t) ] y(t) = [ C 0 ] [ x(t) x c (t) So the closed-loop poles are the eigenvalues of [ ] A BCc [ B N G ] r(t) B c C A c (same as 16 6 except in a different place, gives same closed-loop eigenvalues) regardless of the choice of G and N G and N impact the forward path, not the feedback path Second: if N = 0 and G = L = Bc, then we recover the original implementation, on 16 5 since the controller reduces to: ẋ c (t) = A c x c (t) + B c (y(t) r(t)) = A c x c (t) + B c ( e(t)) u(t) = C c x c (t) With G c (s) = C c (si A c ) 1 B c, then this compensator can be written as u(t) = G c (s)e(t) as before (since the negative signs cancel).
26 Fall Third: Given this extra freedom, what is the best way to use it? One good objective is to select G and N so that the state estimation error is independent of r(t). With this choice, changes in r(t) do not tend to cause such large transients in x(t) Note that for this analysis, take x(t) = x(t) x c (t) since x c (t) ˆx(t) x(t) = ẋ(t) ẋ c (t) = Ax(t) + Bu(t) (A c x c (t) + B c y(t) + Gr(t)) = Ax(t) + B( C c x c (t) + Nr(t)) ({A BC c B c C}x c (t) + B c Cx(t) + Gr(t)) x(t) = Ax(t) + B( Nr(t)) ({A B c C}x c (t) + B c Cx(t) + Gr(t)) = (A B c C)x(t) + B Nr(t) ({A B c C}x c (t) + Gr(t)) = (A B c C) x(t) + B Nr(t) Gr(t) = (A B c C) x(t) + (B N G)r(t) Thus we can eliminate the effect of r(t) on x(t) by setting G B N
27 Fall Fourth: if this generalization does not change the closed-loop poles of the system, then what does it change? Zeros of SISO y/r transfer function are given by: si A BC c B N general det B c C si A c G = 0 C 0 0 original det si A BC c 0 B c C si A c B c C 0 0 = 0 new det si A BC c B N B c C si A c B N C 0 0 = 0 Hard to see how this helps, but consider the scalar new case: C( BC c B N + (si A c )B N) = 0 CB N(BC c (si [A BC c B c C])) = 0 CB N(sI [A B c C]) = 0 So zero of the y/r path is the root of si [A B c C] = 0, which is the pole of the estimator. So by setting G = B N, the estimator dynamics are canceled out of the response of the system to a reference command. Cancelation does not occur with original implementation.
28 Fall So in summary, if the SISO system is G(s) = b(s)/a(s), then with DOFB control, the closed-loop transfer function will be of the form: T (s) = Y (s) R(s) = K γ(s)b(s) g Φ e (s)φ c (s) Designer determines Φ c (s) and Φ e (s) Plant zeroes b(s) remain, unless canceled out Many choices for γ(s), one is to set it to Φ e (s) to cancel out the estimator dynamics. Fifth: select N to ensure that the steady-state error is zero. As before, this can be done by selecting N so that the DC gain of the closed-loop y/r transfer function is 1. y [ C 0 ] [ ] 1 [ ] A BC c B N = I r DC B c C B A c The new implementation of the controller is ẋ c (t) = A c x c (t) + B c y(t) + B Nr(t) u(t) = C c x c (t) + Nr(t) Which has two separate inputs y(t) and r(t) Selection of N ensure that the steady-state performance is good The new implementation gives better transient performance.
29 Fall Figure 13: Example #1: G(s) = (s+8)(s+14)(s+20). Method #1: original implementation. Method #2: original, with the reference input scaled to ensure that the DC gain of y/r DC = 1. Method #3: new implementation with both G = B N and N selected. New method should show less of a transient impact of applying r(t).
30 Fall Figure 14: Example #2: G(s) = 0.94 s Method #1: original implementation. Method #2: original, with the reference input scaled to ensure that the DC gain of y/r DC = 1. Method #3: new implementation with both G = B N and N selected.
31 Fall Figure 15: Example #3: G(s) = (s 8)(s 14)(s 20). Method #1: original implementation. Method #2: original, with the reference input scaled to ensure that the DC gain of y/r DC = 1. Method #3: new implementation with both G = B N and N selected.
32 Fall Dynamic Output Feedback Codes (dofb examp2.m) 1 % Code for Topic 16 2 % Examples of dynamic output feedback 3 % 4 % Jonathan How 5 % 6 fig=1; 7 % system 8 switch model 9 case 1 10 G=tf(8*14*20,conv([1 8],conv([1 14],[1 20])));name= nexamp1 11 tt=[ ]*2.5; 12 t=[0:.0025:1]; 13 case 2 14 G=tf(8*14*20,conv([1-8],conv([1-14],[1-20])));name= nexamp2 ; 15 tt=[ ]*2.5; 16 t=[0:.0025:1]; 17 case 3 18 G=tf(.94,[ ]);name= nexamp3 ; 19 tt=[ ]/5; 20 t=[0:.01:5]; 21 case 4 22 G=tf([1-1],conv([1 1],[1-3]));name= nexamp4 ; 23 tt=[ ]; 24 t=[0:.01:10]; 25 case 5 26 G=tf(conv([1-2],[1-4]),conv(conv([1-1],[1-3]),[1 2*.2*2 2^2 0 0])); 27 name= examp5 ; 28 case 6 29 G=tf(8*14*20,conv([1 8],conv([1 14],[1 20])));name= nexamp1h ; 30 tt=[ ]*2.5; 31 t=[0:.0001:1]; 32 otherwise 33 return 34 end %%% 37 [a,b,c,d]=ssdata(g);na=length(a); 38 % 39 figure(fig);clf;fig=fig+1; 40 % controller gains to place poles at specified locations -- see if model == 6 42 R=.00001; 43 else 44 R=.01; 45 end 46 % choose the regular poles using LQR [k,p,reg_poles]=lqr(a,b,c *c,r); 48 %hold on 49 %plot(reg_poles+j*eps, md, MarkerSize,12, MarkerFaceColor, m ) 50 %hold off % design the estimator by doubling the real part of the regulator poles 53 PP=2*real(reg_poles)+imag(reg_poles)*j; 54 % see ke=place(a,c,pp);l=ke 56 % see ac=a-b*k-l*c;bc=l;cc=k;dc=0; 58 % see acl=[a -b*cc;bc*c ac];ccl=[c,zeros(1,na)];dcl=0; 60 % standard case 61 bcl1=[zeros(na,1);-bc]; 62 %scale up the ref path so that at least the Gcl(s) has unity DC gain 63 bcl2=[zeros(na,1);-bc]; % original implementation on scale=ccl*inv(-acl)*bcl2; 65 bcl2=bcl2/scale; 66 % 67 bcl3=[b;b]; % new implementation on scale=ccl*inv(-acl)*bcl3; 69 bcl3=bcl3/scale; 70 %
33 Fall Gcl1=ss(acl,bcl1,ccl,dcl); % original implementation (16-5) unscaled 72 Gcl2=ss(acl,bcl2,ccl,dcl); % original implementation (16-5) scaled 73 Gcl3=ss(acl,bcl3,ccl,dcl); % new implementation on scaled 74 Gc=ss(ac,bc,cc,dc); f=logspace(-2,3,400); 77 g=freqresp(g,f*j);g=squeeze(g); 78 gc=freqresp(gc,f*j);gc=squeeze(gc); 79 gcl1=freqresp(gcl1,f*j);gcl1=squeeze(gcl1); 80 gcl2=freqresp(gcl2,f*j);gcl2=squeeze(gcl2); 81 gcl3=freqresp(gcl3,f*j);gcl3=squeeze(gcl3); figure(fig);fig=fig+1;clf 84 orient( tall ) 85 subplot(211) 86 loglog(f,abs(g),f,abs(gc), --, LineWidth,2);axis([.1 1e3 1e-2 500]) 87 xlabel( Freq (rad/sec) );ylabel( Mag ) 88 legend( Plant G, Compensator Gc );grid 89 subplot(212) 90 semilogx(f,180/pi*unwrap(angle(g)),f,180/pi*unwrap(angle(gc)), --, LineWidth,2); 91 axis([.1 1e ]) 92 xlabel( Freq (rad/sec) );ylabel( Phase (deg) );grid 93 legend( Plant G, Compensator Gc ) L=g.*gc; figure(fig);fig=fig+1;clf 98 orient( tall ) 99 subplot(211) 100 loglog(f,abs(l),[.1 1e3],[1 1], LineWidth,2);axis([.1 1e3 1e-2 10]) 101 xlabel( Freq (rad/sec) );ylabel( Mag ) 102 legend( Loop L ); 103 grid 104 subplot(212) 105 semilogx(f,180/pi*phase(l. ),[.1 1e3],-180*[1 1], LineWidth,2); 106 axis([.1 1e ]) 107 xlabel( Freq (rad/sec) );ylabel( Phase (deg) );grid 108 % 109 % loop dynamics L = G Gc 110 % see al=[a b*cc;zeros(na) ac];bl=[zeros(na,1);bc];cl=[c zeros(1,na)];dl=0; figure(fig);fig=fig+1;clf 114 margin(al,bl,cl,dl); 115 figure(fig);fig=fig+1;clf 116 rlocus(al,bl,cl,dl) 117 hold on; 118 plot(eig(acl)+eps*j, bd, MarkerFaceColor, b ) 119 plot(eig(a)+eps*j, mv, MarkerFaceColor, m ) 120 plot(eig(ac)+eps*j, rs, MarkerSize,9, MarkerFaceColor, r ) 121 plot(tzero(ac,bc,cc,dc)+eps*j, ro, MarkerFaceColor, r ) 122 hold off;grid on 123 % 124 % closed-loop freq response 125 % 126 figure(15);clf 127 loglog(f,abs(g),f,abs(gcl1),f,abs(gcl2), --, LineWidth,2); 128 axis([.1 1e3.01 1e2]) 129 xlabel( Freq (rad/sec) );ylabel( Mag ) 130 legend( Plant G, Gcl unscaled, Gcl scaled );grid figure(fig);fig=fig+1;clf 133 loglog(f,abs(g),f,abs(gcl1),f,abs(gcl2),f,abs(gcl3), --, LineWidth,2); 134 axis([.1 1e3.01 1e2]) 135 xlabel( Freq (rad/sec) );ylabel( Mag ) 136 legend( Plant G, Gcl1, Gcl2, Gcl3 );grid ZZ=-1+.1*exp(j*[0:.01:1]*2*pi); 139 figure(fig);fig=fig+1;clf 140 semilogy(unwrap(angle(l))*180/pi,abs(l)) 141 hold on; 142 semilogy(unwrap(angle(zz))*180/pi-360,abs(zz), g- ) 143 hold off 144 axis([ ])
34 Fall hold on;plot(-180,1, rx );hold off 146 if max(real(eig(a))) > title( Nichols: Unstable Open-loop System ) 148 else 149 title( Nichols: Stable Open-loop System ) 150 end 151 ylabel( Mag );xlabel( Phase (deg) ) figure(fig);fig=fig+1;clf 154 plot(unwrap(angle(l))*180/pi,abs(l), LineWidth,1.5) 155 hold on 156 plot(unwrap(angle(l))*180/pi,.95*abs(l), r-., LineWidth,1.5) 157 plot(unwrap(angle(l))*180/pi,1.05*abs(l), m:, LineWidth,1.5) 158 plot(unwrap(angle(zz))*180/pi,abs(zz), g- ) 159 plot(-180,1, rx );hold off 160 hold on; 161 semilogy(unwrap(angle(zz))*180/pi-360,abs(zz), g- ) 162 hold off 163 legend( 1, 0.95, 1.05 ) 164 axis([ ]) 165 if max(real(eig(a))) > title( Nichols: Unstable Open-loop System ) 167 else 168 title( Nichols: Stable Open-loop System ) 169 end 170 ylabel( Mag );xlabel( Phase (deg) ) figure(fig);fig=fig+1;clf 173 loglog(f,abs(1./(1+l)),f,abs(l), --, LineWidth,1.5);grid 174 title( Sensitivity Plot ) 175 legend( S, L ) 176 xlabel( Freq (rad/sec) );ylabel( Mag ) 177 axis([.1 1e3 1e-2 100]) figure(fig);fig=fig+1;clf 180 [y1,t]=step(gcl1,t); 181 [y2,t]=step(gcl2,t); 182 [y3,t]=step(gcl3,t); 183 plot(t,y1,t,y2,t,y3,t,ones(size(t)), --, LineWidth,1.5) 184 setlines(2) 185 legend( meth1, meth2, meth3 ) 186 title( Step Response ) for ii=[1:gcf 15] 189 eval([ print -f,num2str(ii), -dpng -r300,name, _,num2str(ii),.png ]) 190 if ii==5; 191 figure(ii);axis(tt) 192 eval([ print -dpng -r300,name, _,num2str(ii), a.png ]) 193 % axis(tt) 194 end 195 end eval([ save,name, R G Gc Gcl1 Gcl2 Gcl3 k l P PP ]) return 200
Topic # /31 Feedback Control Systems
Topic #16 16.30/31 Feedback Control Systems Add reference inputs for the DOFB case Reading: FPE 7.8, 7.9 Fall 2010 16.30/31 16 2 Reference Input - II On page 15-6, compensator implemented with reference
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 informationTopic # /31 Feedback Control Systems. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #15 16.3/31 Feedback Control Systems State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Reading: FPE 7.6 Fall 21 16.3/31 15 2 Combined
More informationTopic # Feedback Control Systems
Topic #15 16.31 Feedback Control Systems State-Space Systems Open-loop Estimators Closed-loop Estimators Observer Theory (no noise) Luenberger IEEE TAC Vol 16, No. 6, pp. 596 602, December 1971. Estimation
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 informationTopic # Feedback Control Systems
Topic #14 16.31 Feedback Control Systems State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or, even if we can change the pole locations. Where do we
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 informationTopic # Feedback Control
Topic #4 16.31 Feedback Control Stability in the Frequency Domain Nyquist Stability Theorem Examples Appendix (details) This is the basis of future robustness tests. Fall 2007 16.31 4 2 Frequency Stability
More informationTopic # /31 Feedback Control Systems
Topic #17 16.30/31 Feedback Control Systems Improving the transient performance of the LQ Servo Feedforward Control Architecture Design DOFB Servo Handling saturations in DOFB Fall 2010 16.30/31 17 1 LQ
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationTopic # /31 Feedback Control Systems
Topic #12 16.30/31 Feedback Control Systems State-Space Systems Full-state Feedback Control How do we change location of state-space eigenvalues/poles? Or, if we can change the pole locations where do
More informationPole placement control: state space and polynomial approaches Lecture 2
: state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state
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 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 informationProfessor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C8 / ME C34 Problem Set 7 Solution Fall Jansen Sheng and Wenjie Chen, UC Berkeley. 35 pts Lag compensation. For open loop plant Gs ss+5s+8 a Find compensator gain Ds k such that the
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
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 informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
More informationIntroduction to Feedback Control
Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System
More informationELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques
CAIRO UNIVERSITY FACULTY OF ENGINEERING ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 00/0 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques [] For the following system, Design a compensator such
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 informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
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 informationSynthesis via State Space Methods
Chapter 18 Synthesis via State Space Methods Here, we will give a state space interpretation to many of the results described earlier. In a sense, this will duplicate the earlier work. Our reason for doing
More informationMODERN CONTROL DESIGN
CHAPTER 8 MODERN CONTROL DESIGN The classical design techniques of Chapters 6 and 7 are based on the root-locus and frequency response that utilize only the plant output for feedback with a dynamic controller
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 information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More informationEECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
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 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 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 informationReturn Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems
Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations
More informationMAE 143B - Homework 7
MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u
More informationModule 08 Observability and State Estimator Design of Dynamical LTI Systems
Module 08 Observability and State Estimator Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha November
More informationECSE 4962 Control Systems Design. A Brief Tutorial on Control Design
ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got
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 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 informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationReview: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system
Plan of the Lecture Review: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system Plan of the Lecture Review: transient and steady-state
More informationOverview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros)
Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros) J. McNames Portland State University ECE 222 Bode Plots Ver.
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 informationEECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C128/ ME C134 Final Wed. Dec. 14, 211 81-11 am Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth 1 points total. Problem Points Score 1 16 2 12
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 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 informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Steady-state error Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling Analysis Design Laplace
More informationFREQUENCY-RESPONSE DESIGN
ECE45/55: Feedback Control Systems. 9 FREQUENCY-RESPONSE DESIGN 9.: PD and lead compensation networks The frequency-response methods we have seen so far largely tell us about stability and stability margins
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 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 informationINTRODUCTION TO DIGITAL CONTROL
ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant
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 informationIntro. Computer Control Systems: F9
Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se
More informationAPPLICATIONS FOR ROBOTICS
Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table
More informationOutline. Classical Control. Lecture 5
Outline Outline Outline 1 What is 2 Outline What is Why use? Sketching a 1 What is Why use? Sketching a 2 Gain Controller Lead Compensation Lag Compensation What is Properties of a General System Why use?
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
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 information] [ 200. ] 3 [ 10 4 s. [ ] s + 10 [ P = s [ 10 8 ] 3. s s (s 1)(s 2) series compensator ] 2. s command pre-filter [ 0.
EEE480 Exam 2, Spring 204 A.A. Rodriguez Rules: Calculators permitted, One 8.5 sheet, closed notes/books, open minds GWC 352, 965-372 Problem (Analysis of a Feedback System) Consider the feedback system
More information6.1 Sketch the z-domain root locus and find the critical gain for the following systems K., the closed-loop characteristic equation is K + z 0.
6. Sketch the z-domain root locus and find the critical gain for the following systems K (i) Gz () z 4. (ii) Gz K () ( z+ 9. )( z 9. ) (iii) Gz () Kz ( z. )( z ) (iv) Gz () Kz ( + 9. ) ( z. )( z 8. ) (i)
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 informationDynamic Compensation using root locus method
CAIRO UNIVERSITY FACULTY OF ENGINEERING ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 00/0 CONTROL ENGINEERING SHEET 9 Dynamic Compensation using root locus method [] (Final00)For the system shown in the
More informationME 132, Fall 2015, Quiz # 2
ME 132, Fall 2015, Quiz # 2 # 1 # 2 # 3 # 4 # 5 # 6 Total NAME 14 10 8 6 14 8 60 Rules: 1. 2 sheets of notes allowed, 8.5 11 inches. Both sides can be used. 2. Calculator is allowed. Keep it in plain view
More informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Routh-Hurwitz stability criterion Examples Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
More informationTRACKING AND DISTURBANCE REJECTION
TRACKING AND DISTURBANCE REJECTION Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 13 General objective: The output to track a reference
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 informationTopic # Feedback Control
Topic #11 16.31 Feedback Control State-Space Systems State-space model features Observability Controllability Minimal Realizations Copyright 21 by Jonathan How. 1 Fall 21 16.31 11 1 State-Space Model Features
More informationRichiami di Controlli Automatici
Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici
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 informationECE382/ME482 Spring 2005 Homework 7 Solution April 17, K(s + 0.2) s 2 (s + 2)(s + 5) G(s) =
ECE382/ME482 Spring 25 Homework 7 Solution April 17, 25 1 Solution to HW7 AP9.5 We are given a system with open loop transfer function G(s) = K(s +.2) s 2 (s + 2)(s + 5) (1) and unity negative feedback.
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
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 informationBangladesh University of Engineering and Technology. EEE 402: Control System I Laboratory
Bangladesh University of Engineering and Technology Electrical and Electronic Engineering Department EEE 402: Control System I Laboratory Experiment No. 4 a) Effect of input waveform, loop gain, and system
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 informationDr Ian R. Manchester
Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign
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 informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
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 informationInternal Model Principle
Internal Model Principle If the reference signal, or disturbance d(t) satisfy some differential equation: e.g. d n d dt n d(t)+γ n d d 1 dt n d 1 d(t)+...γ d 1 dt d(t)+γ 0d(t) =0 d n d 1 then, taking Laplace
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 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 information6 OUTPUT FEEDBACK DESIGN
6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method
.. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to
More informationSchool of Mechanical Engineering Purdue University. ME375 Feedback Control - 1
Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System
More informationAutomatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback
More informationPrashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides. Department of Chemical Engineering University of California, Los Angeles
HYBRID PREDICTIVE OUTPUT FEEDBACK STABILIZATION OF CONSTRAINED LINEAR SYSTEMS Prashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides Department of Chemical Engineering University of California,
More informationStochastic Optimal Control!
Stochastic Control! Robert Stengel! Robotics and Intelligent Systems, MAE 345, Princeton University, 2015 Learning Objectives Overview of the Linear-Quadratic-Gaussian (LQG) Regulator Introduction to Stochastic
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More informationTopic # /31 Feedback Control Systems
Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationDue Wednesday, February 6th EE/MFS 599 HW #5
Due Wednesday, February 6th EE/MFS 599 HW #5 You may use Matlab/Simulink wherever applicable. Consider the standard, unity-feedback closed loop control system shown below where G(s) = /[s q (s+)(s+9)]
More informationMechanical Systems Part A: State-Space Systems Lecture AL12
AL: 436-433 Mechanical Systems Part A: State-Space Systems Lecture AL Case study Case study AL: Design of a satellite attitude control system see Franklin, Powell & Emami-Naeini, Ch. 9. Requirements: accurate
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More information(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:
1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.
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 informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC426 SC426 Fall 2, dr A Abate, DCSC, TU Delft Lecture 5 Controllable Canonical and Observable Canonical Forms Stabilization by State Feedback State Estimation, Observer Design
More information