1 / 55

Lecture 40: Review of Controls

Lecture 40: Review of Controls. Converting a problem to mathematics. Linearization. Stabilization (making x go to zero). Tracking (making x go to x r — the hardest). Observers (stabilizing the system when you can’t measure the entire state). Converting a problem to mathematics.

ornice
Download Presentation

Lecture 40: Review of Controls

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Lecture 40: Review of Controls Converting a problem to mathematics Linearization Stabilization (making x go to zero) Tracking (making x go to xr— the hardest) Observers (stabilizing the system when you can’t measure the entire state)

  2. Converting a problem to mathematics You can write the physics using free body diagrams or the Euler-Lagrange method There may be an auxiliary electrical equation Let’s recall magnetic suspension; it has a number of useful features.

  3. Converting a problem to mathematics yis positive down Numerical parameters for later

  4. Converting a problem to mathematics We’ll need to convert the second order dynamical equation to a pair of first order equations We’ll want to rearrange the electrical equation to look like a state equation Our state vector will be where v denotes the rate of change of y

  5. Converting a problem to mathematics Manipulate the equations combine these input

  6. Converting a problem to mathematics and I have replaced e by the generic input variable u This is a nonlinear problem, which brings us to the next step

  7. Linearization We linearize about some reference state, some point in state space that is stationary This is frequently, but not always, zero For the magnetic suspension system is it not zero The reference state may include a reference input, as it does for the magnetic suspension

  8. Linearization Find the reference state We have Clearly v0 = 0. y0 cannot be zero, so we have Finally

  9. Linearization We are at liberty to select one of these variables. y0 makes the most sense. The reference state is then We get the linear equations by writing x = x0 + x’, e = e0 + e’

  10. Linearization One way, the brute force and awkwardness way, is to substitute for x and e and expand all the terms, keeping only those linear in the primed variables This works if the problem is simple, and it can be done if the problem is hard The better way is to differentiate the nonlinear force F

  11. Linearization Linearization by differentiation Think of expanding F in a three variable Taylor series around the equilibrium

  12. Linearization substitute the equilibrium values

  13. Linearization Remember that So that the expressions become And these are the columns of the linearized matrix

  14. Linearization Let’s look at another example of linearization by differentiation Find an equilibrium such that this vanishes and linearize about that equilibrium

  15. Linearization Take the simple equilibrium, for which x = 0

  16. Linearization The more complicated equilibrium is In this case we have

  17. QUESTIONS?

  18. Stabilization We have arrived at the canonical problem where our goal is to choose u in such a way that x will go to zero We do that using full state feedback, setting u = -gx We need to findg, and we have a ritual to do that The ritual only works if the system is controllable and the first step in the ritual is figuring out if the system is controllable. It’s probably worth noting that: A is a kxk square matrix gis a kx 1 column vector bis a 1 xk row vector

  19. Stabilization Step one is to form the controllability matrix and calculate its determinant. If the determinant is nonzero the system is controllable and we can move forward Find the inverse ofQand select its last row This is a row vector. We calculate k – 1 more rows to make a square matrix T.

  20. Stabilization And we have the square matrix which is invertible We use T to do a change of variable (which just amounts to a change of basis for the mathematically inclined of you) We see what happens on the next slide

  21. Stabilization We can multiply this last equation by T to get We agree to rename the matrices and vectors to make this equation A1 and B2 are in companion form; here’s a k = 6 example

  22. Stabilization Now we are ready to talk about gains, and we have to start in terms of z x will go to zero if z goes to zero. z will go to zero if the eigenvalues of Ac all have negative real parts. We call these eigenvalues poles. Because of the companion form, it is easy to choose gains to put the poles where we want them

  23. Stabilization The k = 6 example The characteristic polynomial Remember this ritual for finding characteristic polynomials You don’t need to calculate determinants

  24. Stabilization We can form the characteristic polynomial for the poles we want Expand that in powers of s and equate the coefficients to determine the gains (in z space) Once we have those, we need to go back tox space, and let’s see how that goes

  25. Stabilization Were there but worlds enough and time I would now have you work a stabilization problem But there ain’t, so I’ll work one

  26. Stabilization The inverse: swap the diagonal terms change the signs of the off-diagonal terms divide by the determinant Build T

  27. Stabilization The inverse: swap the diagonal terms change the signs of the off-diagonal terms divide by the determinant The companion forms

  28. Stabilization Now we can form Ac and work on its eigenvalues (the poles) Characteristic polynomials

  29. Stabilization The choice of poles is up to you Critical damping: s1 = -r = s2 Conjugate pair: s1 = -r+ jq , s2 = -r- jq

  30. QUESTIONS?

  31. Tracking This is probably the hardest one, and we haven’t spent much time on it Find some ur such that Find some Arsuch that Two options present themselves to us: We found the second method to work better in the one problem we did with both

  32. Tracking We compare with Friedland equation (6.37)

  33. Tracking We can ask that the derivative of the error goes to zero, from which is stable, hence invertible and we can write the latter form being equivalent to equation 6.40 We find g in the same manner we just rehearsed

  34. Tracking It may happen that we can actually make e vanish in its entirety If that is the case, then we will be very happy It is more generally true that we cannot that all we can make vanish is some subset of e typically the part that looks like an error in y This is what Friedland does, multiplies the error by the matrix/vector C For SISO systems, it is the vector c, I’ll leave it a vector for now

  35. Tracking rearrange we want this to work for any xr so . . .

  36. Tracking Let’s look at another 2 x 2 system and see if we can make it track

  37. Tracking From the earlier part of this lecture, we can find the gains given the poles I will take them to be at -1, giving me a critically damped system

  38. Tracking We can build the error matrix and in this case we can actually kill all the error by choosing We select the feedback as we had before

  39. Tracking We can plug all this in and integrate numerically. Here’s the result for Y = 1, w = 2

  40. QUESTIONS?

  41. Observers This last bit depends on the Lecture 38 (and I’m sticking to a single output) Suppose that all we know is an output y, a scalar for an SISO system

  42. Observers We use y, the control and an artificial dynamical system to find an estimate for x and we select the carat matrices to minimize the difference between x and its estimate.

  43. Observers As usual we want to construct a system for the evolution of the error and, of course, make the error go to zero Note that k is a column vector We eliminate the explicit appearance of the estimate

  44. Observers e will go to zero if I get rid of the contributions from x and u and make sure the new matrix is stable This becomes a pole placement exercise

  45. Observers The substitutions give and the error will vanish of the eigenvalues have negative real parts. We need to find k such that this is so. We have two questions: Is this possible? How?

  46. Observers For our SISO system the K vector is a column and the system looks like We can make this look like something we already know how to do by taking the transpose The eigenvalues of the transpose of a matrix are the same as those of the original matrix

  47. Observers cT looks like b and kT looks like g We have the analog of the controllability theorem The observability theorem We can build an observer if det(N) ≠ 0.

  48. Observers

  49. Observers Let’s apply this to a simple model problem I will skip the calculation of the regular gain matrix, because we’ve just gone through one The result is

  50. Observers

More Related