Fa2004 16.3339-5 Find the closed-loop dynamics Ax+B(r-Ka =(A- BK)a+Br ld.C+B 4 C Objective: Pick K so that Ac has the desired properties, e. g A unstable, want A stable Put2 poles at-2±2 o Note that there are n parameters in k and n eigenvalues in a, so it looks promising, but what can we achieve Example #1: Consider 0 Then det(s-A)=(s-1)(s-2) 3s+1=0 so the system is unstable -Define u=-[ki k2]=-Ka,then 1-k11-k2 Ad=a-BK= [k1k2] So then we have that det(I-Aa)=s2+(k1-3)s+(1-2k1+k2)=0
� � � � � � � � Fall 2004 16.333 9–5 • Find the closedloop dynamics: x˙ = Ax + B(r − Kx) = (A − BK)x + Br = Aclx + Br y = Cx • Objective: Pick K so that Acl has the desired properties, e.g., – A unstable, want Acl stable – Put 2 poles at −2 ± 2j • Note that there are n parameters in K and n eigenvalues in A, so it looks promising, but what can we achieve? • Example #1: Consider: 1 1 x˙ = 1 2 1 x + u 0 – Then det(sI − A) = (s − 1)(s − 2) − 1 = s2 − 3s + 1 = 0 so the system is unstable. – Define u = − k1 k2 x = −Kx, then � � � � 1 1 1 � � 1 − k1 1 − k2 Acl = A−BK = k1 k2 = 1 2 − 0 1 2 – So then we have that det(sI − Acl) = s2 + (k1 − 3)s + (1 − 2k1 + k2) = 0
Fa2004 16.3339-6 Thus, by choosing ki and k2, we can put Ai(Ac) anywhere in the complex plane(assuming complex conjugate pairs of poles) To put the poles ats=-5,-6, compare the desired characteristic equation (s+5)(s+6)=s2+118+30=0 with the closed-loop one s2+(k1-3)x+(1-2k+k2)=0 to conclude that k1-3=11 k1=14 1-2k1+k2=30 so that K=[14 57, which is called Pole Placement Of course, it is not always this easy, as the issue of controllability must be addressed Example #2: Consider this system 11 0 with the same control approach 1-k11-k2 A=A-Bk so that det(si- ad)=(s-1+k1(s-2)=0 The feedback control can modify the pole at s but it cannot move the pole at s= 2 This system cannot be stabilized with full-state feedback control o What is the reason for this problem?
� � � � � � � � � � � � � Fall 2004 16.333 9–6 – Thus, by choosing k1 and k2, we can put λi(Acl) anywhere in the complex plane (assuming complex conjugate pairs of poles). • To put the poles at s = −5, −6, compare the desired characteristic equation (s + 5)(s + 6) = s2 + 11s + 30 = 0 with the closedloop one 2 s + (k1 − 3)x + (1 − 2k1 + k2) = 0 to conclude that k1 − 3 = 11 k1 = 14 1 − 2k1 + k2 = 30 k2 = 57 so that K = 14 57 , which is called Pole Placement. • Of course, it is not always this easy, as the issue of controllability must be addressed. • Example #2: Consider this system: 1 1 x˙ = 0 2 1 x + u 0 with the same control approach 1 1 1 � � 1 − k1 1 − k2 Acl = A − BK = k1 k2 = 0 2 − 0 0 2 so that det(sI − Acl) = (s − 1 + k1)(s − 2) = 0 The feedback control can modify the pole at s = 1, but it cannot move the pole at s = 2. • This system cannot be stabilized with fullstate feedback control. • What is the reason for this problem?
Fa2004 16.3339-7 It is associated with loss of controllability of the e mode Basic test for controllability: rank Mc=n BlAB So that rank M=1<2 Must assume that the pair(A, B) are controllable
� � � � � Fall 2004 16.333 9–7 – It is associated with loss of controllability of the e2t mode. • Basic test for controllability: rank Mc = n � � � � � B AB 1 1 = 1 0 1 0 Mc = 0 2 So that rank Mc = 1 < 2. • Must assume that the pair (A, B) are controllable
Fa2004 16.3339-8 Ackermann's Formula The previous outlined a design procedure and showed how to do it by hand for second-order systems - Extends to higher order(controllable) systems, but tedious Ackermann's Formula gives us a method of doing this entire design process is one easy step K=[0…01]Ma重 [BAB…A-B] dpa(s is the characteristic equation for the closed-loop poles. which we then evaluate for s= a It is explicit that the system must be controllable because we are inverting the controllability matrix · Revisit Example#1:重(Ss)=s2+115+30 11 BlAB 0 So 11 11 01 +301 01 4314 =1457 Automated in Matlab: place. m acker m(see polyvalmm too
� � � � � � � � �� �� Fall 2004 16.333 9–8 Ackermann’s Formula • The previous outlined a design procedure and showed how to do it by hand for secondorder systems. – Extends to higher order (controllable) systems, but tedious. • Ackermann’s Formula gives us a method of doing this entire design process is one easy step. K = 0 . . . 0 1 M−1 Φd(A) c – Mc = B AB . . . An−1B – Φd(s) is the characteristic equation for the closedloop poles, which we then evaluate for s = A. – It is explicit that the system must be controllable because we are inverting the controllability matrix. • Revisit Example #1: Φd(s) = s2 + 11s + 30 � � � � � � � 1 1 1 � 1 1 1 Mc = B AB = = 0 1 2 0 0 1 So ⎞ � � 1 1 �−1 ⎛� �2 � � � 1 1 1 1 K = 0 1 ⎝ + 11 + 30I⎠ 0 1 1 2 1 2 � � 43 14 � � = 0 1 = 14 57 14 57 • Automated in Matlab: place.m & acker.m (see polyvalm.m too)
Fa2004 16.3339-9 Origins? For simplicity, consider a third-order system(case #2),but this extends to any order ai-a 100B=0C=[bb2b3] This form is useful because the characteristic equation for the system is obvious => det(sI-A)=s+a1S+a2S+a3=0 Can show that a A=A-BK 100 0[k1 k2 ka 1-k1-a2-k2 so that the characteristic equation for the system is still obvious 重a(s)=de-Aa)=s3+(a1+k)s2+(a2+k2)s+(a3+k3)=0 We then compare this with the desired characteristic equation devel oped from the desired closed -loop pole locations 重(s)=s32+(a1)s2+(a2)s+(a3)=0 to get that +hl +k kn Pole placement is a very powerful tool and we will be using it for most of our state space work
� � � � Fall 2004 16.333 9–9 • Origins? For simplicity, consider a thirdorder system (case #2), but this extends to any order. ⎡ ⎤ ⎡ ⎤ −a1 −a2 −a3 1 A = ⎣ 1 0 0 ⎦ B = ⎣ 0 ⎦ C = b1 b2 b3 0 1 0 0 – This form is useful because the characteristic equation for the 2 system is obvious ⇒ det(sI − A) = s3 + a1s + a2s + a3 = 0 • Can show that ⎡ ⎤ ⎡ ⎤ −a1 −a2 −a3 1 Acl = A − BK = ⎣ 1 0 0 ⎦ − ⎣ 0 ⎦ k1 k2 k3 0 1 0 0 ⎡ ⎤ −a1 − k1 −a2 − k2 −a3 − k3 = ⎣ 1 0 0 ⎦ 0 1 0 so that the characteristic equation for the system is still obvious: 2 Φcl(s) = det(sI −Acl) = s3+(a1+k1)s +(a2+k2)s+(a3+k3) = 0 • We then compare this with the desired characteristic equation developed from the desired closedloop pole locations: 2 Φd(s) = s3 + (α1)s + (α2)s + (α3) = 0 to get that ⎫ a1 + k1 = α1 ⎬ k1 = α . 1 − a1 . . . . . ⎭ an + kn = αn kn = αn − an • Pole placement is a very powerful tool and we will be using it for most of our state space work