Linearizing Around a Fixed Point [Control Bootcamp]
Ғылым және технология
This lecture describes how to obtain linear system of equations for a nonlinear system by linearizing about a fixed point. This is worked out for the simple pendulum "by-hand" and in Matlab.
Chapters available at: databookuw.com/databook.pdf
These lectures follow Chapter 8 from:
"Data-Driven Science and Engineering: Machine Learning, Dynamical Systems, and Control" by Brunton and Kutz
Amazon: www.amazon.com/Data-Driven-Sc...
Book Website: databookuw.com
Brunton Website: eigensteve.com
This video was produced at the University of Washington
Пікірлер: 95
This is an "honest-to-goodness" fantastic series. Thanks !
@architektwo
2 ай бұрын
I want to like this comment but the like count is too nice. Also, is anyone else reminded of Mr. Mackey mkay?
A friend of mine's daughter waited till the last minute to ask for help from me for her graduate ODE class. Happens. But I haven't taken ODE for decades and when I did, well, let's say the results were less than impressive. My friend's daughter's professor seemed, from the murky lecture notes I got, to not be a great communicator, so I needed to get a quick warm-up on this topic. Thank goodness I ran into your clear-as-a-bell lecture series. College professors, beware! There is some serious online talent who will teach your kids if you're too lazy to write coherent notes and give coherent lectures. In short, bravo! You were a life-saver!!!
I joined my master's course after 7 years being away from college. and your videos are my saviour for my control systems class! Thank you so much for such brilliant inductive way of introducing and explaining concepts!
I just wanna call this professor and let him know he is awesome. I got what exactly I needed. Learn something and see how it can be applied in real world. You are just amazing.
Thank you Sir....I have seen the whole playlist and it cleared a lot of my concepts about control theory. Your videos are just great and your way of teaching complex things in simple manner is appreciable. Thanks Again.
just finished about 16 lectures from ur ME564 course and wanted to start getting into the control playlist .. i gotta say having examples with numbers after evey part of the lecture + the matlab examples was a luxury :D
You have some great instructional video for dynamical systems, stability, and control (linear and nonlinear). This series of videos would be very helpful to new graduate level control engineering students, and I’d recommend they watch them all.
Thanks a lot, I've just started the series and im in love with them. Really good work!
Your manner of teaching is very amazing!
Bro you are honestly the best. You make everything easier to understand.
@Eigensteve
3 жыл бұрын
That is really nice of you!
I did not know I needed this series
Thank you Sir, this is another life-saving Series. Now I feel much better and more clear about what I'm dealing with in Modern Control
@Eigensteve
7 ай бұрын
You're welcome!
Hi, I want to say that your videos are just AMAZING!!! no words to say, this series also amazing like your other series. but why you have stopped posting more such videos on this bootcamp playlist? I am eagerly waiting for more videos on this bootcamp, apecially MPC controller and controlling of Non-linear systems. Please keep this playlist running and post more such contents please....
Thank you so much for the series. It really gives me a chance to pass the test with only one day studying. The cool part is, during the semester I wasn't interested at all in this subject (hence the one day study) but seeing your videos, and understanding these concepts so easily, now I find this area really fascinating, and will learn more about it. So thanks again for the solid explainers and for showing me how interesting this subject really is.
Thanks prof. for sharing this piece of art with public
Thank you so much for sharing your knowledge. This is truly remarkable and amazing content
Thank you for posting this excellent series!
@Eigensteve
3 жыл бұрын
Glad you enjoy it!
@kubigaming3606
3 жыл бұрын
@@Eigensteve if i could have lecturers like You :0. Thanks for video, it was very helpful :]
This is so good. Tanks a lot for this lectures.
Multiplying by delta(x) - top tip! Makes life helluva easier
Seriously, you are awesome!!. Thank you very much for all this information.
I have never been the brightest student, so even I like and find fascinating control it was alwasy dificult for me. Thank you this was really easy to understand.
This is an absolutely fantastic video series. Thanks for making these!
@messiahjustin1964
2 жыл бұрын
i dont mean to be so off topic but does someone know a trick to get back into an instagram account? I was dumb forgot my account password. I would appreciate any tips you can offer me!
@jaycemarcelo2634
2 жыл бұрын
@Messiah Justin Instablaster ;)
@messiahjustin1964
2 жыл бұрын
@Jayce Marcelo Thanks for your reply. I got to the site through google and I'm in the hacking process now. Seems to take quite some time so I will get back to you later with my results.
@messiahjustin1964
2 жыл бұрын
@Jayce Marcelo it did the trick and I finally got access to my account again. I'm so happy! Thank you so much, you really help me out!
@jaycemarcelo2634
2 жыл бұрын
@Messiah Justin glad I could help :)
Thanks for this amazing series! Question: if the pendulum had no dampening, could we still linearize the system around the fixed points? The eigenvalues of the Jacobian would be +/- i, which the Hartman-Grobman theorem (IIUC) says would make it unsuitable for simple linearization.
Wonderful .. Highest quality of presentation
3blue1brown has a really good video about Jacobian and how it represents linearization. Definitely worth checking out.
@vikeshbubbles205
3 жыл бұрын
Could you point me to this? I haven't been able to find it
@darshangadginmath838
3 жыл бұрын
@@vikeshbubbles205 kzread.info/head/PL2QY0xcsWhz6ghPUixeNQm47xapSKw4SM
@jay89boy
2 жыл бұрын
@@vikeshbubbles205 kzread.info/dash/bejne/qXeOk7h-nbnQhtI.html
@humzahkhan6299
Жыл бұрын
@@vikeshbubbles205 kzread.info/dash/bejne/lKOcrptqaM6-gc4.html&ab_channel=KhanAcademy
Jacibain is a covector (if $f$ just a function of several vars), means vector-row. So, \partial f_i are in the rows
Thank you Professor! 13:35. Hartman-Grobman theorem
@ManuelMenzella04
3 жыл бұрын
Does this imply that an undampened pendulum cannot be linearized around the fixed points? The eigenvalues of the Jacobian would be +/- i, with no real part.
Fantastic video!! What were to happen if the equilibrium points were 0 for both? Would there be any reason to linearize the problem?
very nice, can anyone upload some more practice problems where we need to determine System Matrix(A) for given Differential equations ?
Love this series. Are we able to linearize around non-fixed points? For many real world mechanical problems, we're not trying to control to an equilibrium position, so I guess this must already be a solved problem.
@Eigensteve
2 жыл бұрын
Awesome, glad you like it! Yes, we can definitely linearize about non-fixed points. The vector field will flow through this point though, so it is not an equilibrium, and there will be a "constant" vector field plus the Ax term
Amazing as usual
Hello Steve. Thanks for the lecture. I cannot help but think that we cheated a bit because we have some previous knowledge of the physics and we know about those two fixed points. But I’m wondering about the sort of systems where we “honest to goodness” don’t know if there are fixed points and where do they could be although we expect them to exist because we have this intuition it can become stable even if we can tell nothing else. Thanks for supplement the lecture with Matlab, it makes it somewhat easier to grasp from an engineering perspective.
@Eigensteve
3 жыл бұрын
Really good question. You are right, usually we start with some knowledge of the system. When we truly have no knowledge, we can sometimes use system identification and then analyze the resulting system. We have a method (sparse identification of nonlinear dynamics; SINDy) that will give equations of motion from data, and then these can be analyzed. But it is extra tricky when the system is naturally unstable. (we have a video on the SINDy method too)
Excellent explanation. When you are obtaining the fixed points don't you have to consider the control signal/force?
The A matrix for down position of the inverted pendulum without friction has eigenvalues +-i which are purely imaginary. Is it correct to point out that this violates the Hartman-Grobman condition and hence linearization might not work in this case?
Great series, I'm learning a lot. Not sure if this comment is getting lost, may be posting multiple times, so apologies up front. On the dynamics of the pendulum, you have the friction/dissipation term as $-\delta \dot{\theta}$. Since friction works to oppose the acceleration regardless of which way it is swinging, shouldn't the term use the absolute value of theta-dot?: $-\delta \lvert \dot{\theta} vert$? I am new to this, pulling on old physics education, so I may be misunderstanding something. Or perhaps it doesn't affect the main thrust of the lesson, and doesn't materially affect the Jacobian at the fixed point? Thanks if you are able to address this question.
This is really cool video. Thank a lot.
Is there a reason why we linearize around fixed points, apart for the sake of simplicity ? Couldn't we, almost as easily, use an affine approximation around any point, even if its not an equilibrium, and then design a controller that stabilizes the system around that point ?
@alexandreblanc7253
3 жыл бұрын
Ok so I actually worked it out on a piece of paper, by linearizing around theta = pi/3 in a damped pendulum. As it turns out, the dynamics is stable. Crazy, right ? Well it is stable but the equilibrium point is not theta/3, which would be what we want. We linearize around points of equilibrium because if we are not around a point of equilibrium the approximation does not hold in the long run, and is thus not interesting to study stability. The general workflow is to first shift the equilibrium (in that case by adding a constant term to the control -f(x bar)), then stabilize around this new equilibrium. edit : see www.cds.caltech.edu/~murray/courses/cds101/fa02/faq/02-10-09_linearization.html
Why does a hyperbolic point have neighboring stable and unstable manifolds, but not a center manifold? Don't you need an "inflection point" between them, which would have eigenvalues of zero and be a center manifold?
just a one word 'Awesome'
I'm confused between this and the extended the kalman filter. In EKF we linearize about the mean and not a stationery point. Why is that any thoughts?
You are amazing! Thanks!
@Eigensteve
3 жыл бұрын
Thank you too!
This is awesome, but I'm in an application where I need to accurately control a pendulum throughout its entire range of motion, not just around the fixed points. Is there a solution for linearizing the dynamics around an arbitrary point?
This makes me want to do another masters
Coming back to this video, I have a question. Usually, the theoretical justification for representing a NL system close to an equilibrium with its linearized version is the Hartman-Grobman theorem, which states that for non-hyperbolic equilibria the phase portraits of the NL and linearized systems are topologically equivalent. However, for an undamped pendulum the Jacobian matrix at the stable equilibrium point has two purely imaginary eigenvalues, which actually are hyperbolic (they have zero real part). So how come the linearization does not fail? What piece of information am I missing here? Thanks to anyone who is willing to answer!
you are amazing ❤
“Honest to goodness” love this term😂
I have a Question at 29:47 min: What is the stable and the unstable direction? Are those "directions" represented physically somehow?
Sir is this book available in paperback format in India..?? Hard cover book is very very expensive & as a student I cannot afford..
I am using this method for the cart / inverted pendulum problem in video 12. I can linearize the A matrix with the Jacobian by setting u = 0, but then what do I do with the B matrix which also has nonlinear terms in it? Can I also take the Jacobian of the B matrix and evaluate that around the stability points? Or does even solving the A matrix require knowing what u(t) is prior, since the variables vary with u? Link to a text file with my code: drive.google.com/file/d/11aVEOvqi4T7-qTYThAk-zAFcJJEKwMWa/view?usp=sharing EDIT: I should clarify, my code is actually for a rod with uniform mass on a cart, so not exactly the same as the inverted pendulum.
@ryanshaw9301
4 жыл бұрын
I found this video here, which goes through how to linearize a non-linear system dx = Ax +Bu, answers my question (beginning around 9:30) m.kzread.info/dash/bejne/dKyjksF6ZafWY9o.html
I have many questions. How would I contact you?
Why there isnt a factorial term in the taylor series? Great lecture, btw.
@anantchopra1663
4 жыл бұрын
Yeah, he missed it in the third term (2!), but it didn't really matter since we neglected that term anyway.
15:30 sir your cleaning seems to have an adaptive control with awesome actuators, no marks left :). apart from joke this video and this whole lecture are unmatched.
Is the approximation at the fixed point is a Taylor expansion?
@anantchopra1663
4 жыл бұрын
Yes! In general, we can linearize about any point, but linearizing about a fixed point gives the advantage that the first term in the Taylor series i.e. f(x-bar) = 0 and so we do not get a constant in the linear equation.
Thanks.
Thanks
I'm just curious, how did he write oppositly on board?
The pleasure is mine:)
What if you want to control to a non stable point?
Fantastic class except for the shrill noises ...
where is 2nd part?
Great teacher thank you
=) loved it.
Man oh Man, amazing content! But if I may give one humble opinion, it would be a lot better if you could just present your lectures the usual way rs, Just a good old, run of the mill black board without the fancy refletive of shooting it. But that might be just me, keep up the great work though.
👌
Ugh still above my head
is he writing in reverse?
@suyashmuley9096
3 жыл бұрын
@@mouseShyu He is writing on the mirror!
Did we already cover how to find the fixed point? It is not clear to me how to solve f(x) = 0 here
@abhinavraghunandankannan3546
3 жыл бұрын
He goes over how to calculate the fixed point in this video. Here f(x) is the following vector ( x2, -sin(x1) - delta*x2)^T Since f(x) = 0, x2 = 0 and sin(x1) = 0 sin(x1 ) is 0 when x1 = 0, n*pi which are fixed points.
Since you missed the 1/2 in the Taylor series you must mean it when you say you don’t use nonlinear approximations often.
squeak squeak squeak squeak
@amitra2692
2 ай бұрын
Teri Ammi ka halala
As I expected, the guy forgot about 1/2 in the third term of the Taylor expansion (see 9:33). Of course, higher orders are later negclected, but still it does not look fine.