Linear Systems of Equations, Least Squares Regression, Pseudoinverse

  Рет қаралды 88,225

Steve Brunton

Steve Brunton

Күн бұрын

This video describes how the SVD can be used to solve linear systems of equations. In particular, it is possible to solve nonsquare systems (overdetermined or underdetermined) via least squares regression and the pseudoinverse.
These lectures follow Chapter 1 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
Book PDF: databookuw.com/databook.pdf
Brunton Website: eigensteve.com
This video was produced at the University of Washington

Пікірлер: 90
@dragoncurveenthusiast
@dragoncurveenthusiast 3 жыл бұрын
This is so cool! When I started this lecture series in order to understand PCA better, I had no idea it would also relate to least squares regression! This blew my mind! Thank you so much for making these. They must be a lot of work, but they are so appreciated!
@neoblackcyptron
@neoblackcyptron 3 жыл бұрын
Thank you for explaining hard to grasp concepts in a filtered simple manner for us to understand. Your lectures are a great complement to prof strangs both high quality content.
@philmccavity
@philmccavity 6 ай бұрын
Exceptionally clear explanation, crisp hand-written notes, wonderful!
@leo1fun
@leo1fun 3 жыл бұрын
I'm always impressed by how clean the board is, so it looks like there's nothing at all.
@pranjalsahu
@pranjalsahu 3 жыл бұрын
Excellent Lecture ! So clear to understand! Thank You !
@chengkeattan6571
@chengkeattan6571 2 жыл бұрын
Thank you very much for the clear explanation of pseudo-inverse.
@macmos1
@macmos1 4 жыл бұрын
please keep going with the numerical linear algebra/numerical analysis/scientific computation/applied math stuff thanks :)
@idkravitz
@idkravitz 4 жыл бұрын
Am I right, that you write on real glass in front of camera and the image is just mirrored by editing? If so its brilliant.
@bastianian2939
@bastianian2939 3 жыл бұрын
Damn i just concluded he was an expert at writing mirrored
@sihebi973
@sihebi973 3 жыл бұрын
The real video will look like this: www.mirrorthevideo.com/watch?v=PjeOmOz9jSY He is also using his left hand.
@reihanehvafadar6139
@reihanehvafadar6139 4 жыл бұрын
can't thank you enough for sharing your knowledge with entire world
@PunmasterSTP
@PunmasterSTP Жыл бұрын
Pseudoinverse? More like "Super videos for us!" Thank you so much for making all of them.
@amodamatya
@amodamatya 3 жыл бұрын
Thank you Professor for this valuable lecture
@mauzaomin3872
@mauzaomin3872 4 жыл бұрын
Thank you Steve for video. We make the assumption that it is an economy SVD at time 6:28. Then, how can we guarantee that V multiplies with V* will become identity matrix, especially for the under-determined system?
@Atlas-ck9vm
@Atlas-ck9vm 4 жыл бұрын
Absolutely Great Content
@zachfang4424
@zachfang4424 3 жыл бұрын
Thank you for the awesome materials and everything is well explained! One question I have is how we calculate the inverse of the singular matrix which is a non-square matrix? Isn't that back to the problem, i.e. inverse the non-square A matrix, we had at the first place?
@u2coldplay844
@u2coldplay844 3 жыл бұрын
I loves your lectures! it is so clear and save us from mist of information. I noticed that you wrote Sigma matrix is a invertible, if A is not invertible, shouldn't the Sigma matrix also is not invertible, hence Pseudo inverse of Sigma matrix? Thank you for clarification
@amaniarman460
@amaniarman460 3 жыл бұрын
am I just missing something, why is n
@utatistics9293
@utatistics9293 4 жыл бұрын
having second thoughts about doing master cuz your videos are just too helpful
@patf9770
@patf9770 3 жыл бұрын
with all the amazing resources on the internet, it seems like higher ed is turning into mostly gatekeeping
@mikatshow3932
@mikatshow3932 2 жыл бұрын
Great explanation, sir!
@ehkim1977
@ehkim1977 3 жыл бұрын
This is really nice lecture I've ever seen. I want to recommend this lecture for engineering graduate student basic class!! :-) Thank you so much~ I'll buy the book! ^-^
@Eigensteve
@Eigensteve 3 жыл бұрын
Awesome, thanks so much, and hope you like the book!
@JohnJTraston
@JohnJTraston Жыл бұрын
Yeah. This is a very good and useful lecture.
@drditup
@drditup 3 жыл бұрын
I'm amazed. This is so clearly explained!!
@1985lama
@1985lama 2 жыл бұрын
In the case of over-determined matrix X, why VVT is equal to identity since we are using economy matrices?
@FelipeCondo
@FelipeCondo 3 жыл бұрын
Thanks for the video Professor. Could you help me with something please? I am trying to fit a sinusoidal surface the depends on (x,y,t) but in my case b is not a vector but a matrix, what can I do in this case? Thank you
@julianschmitz9525
@julianschmitz9525 Ай бұрын
great video
@annflintoft9222
@annflintoft9222 Жыл бұрын
Wonder if Steve Brunton can cover for under/over determined systems the nontrivial solutions of Ax = 0, using SVD.
@EngineeringChampion
@EngineeringChampion 3 жыл бұрын
Since A has fewer rows than columns, then A is an m×n matrix with m
@ahsanahmed2505
@ahsanahmed2505 2 жыл бұрын
so, he made a mistake?
@FortranCastle
@FortranCastle 9 ай бұрын
@@ahsanahmed2505 He is using n by m in the videos and in his book instead of the usual convention of m by n which is quite confusing and contrary to most linear algebra resources.
@jay89boy
@jay89boy 3 жыл бұрын
Great content !
@Tyokok
@Tyokok 2 жыл бұрын
Thanks for the great video! One question: 7:50, if you have zero singular value, how do you compute sigma inverse? Thank you!
@diogenescruz-figueroa2719
@diogenescruz-figueroa2719 Жыл бұрын
This might be a year too late, but the Moore-Penrose Invert satisfies (AB)+ = B+A+. So I think there was a mistake, and it should have been Sigma+, not Sigma^{-1}, so you separate it from U and V. Since U and V are unitary, they are invertible, and thus their Moore-Penrose inverse is the regular inverse. As for Sigma, since it has only elements on the diagonal, its Moore-Penrose inverse is just the transpose, and then have the reciprocals for each of its elements (an "d" in the diagonal becomes a "1/d").
@eddiechen6389
@eddiechen6389 2 жыл бұрын
Wonderful lecture
@giorgoschristopoulos9087
@giorgoschristopoulos9087 3 жыл бұрын
In an underdetermined case, if you use the economy SVD, is V*V' equal to an identity?
@ahsanahmed2505
@ahsanahmed2505 2 жыл бұрын
I am confused with the dimension of A matrix. in overperformed, shouldn't overdetermined have m>n
@ankitsinha4124
@ankitsinha4124 3 жыл бұрын
Man. It was really awesome..
@3d_chip
@3d_chip 7 күн бұрын
My god, you just explained what my professor is trying to explain for 5 lectures
@maettu102
@maettu102 3 жыл бұрын
Thanks for the video was really nice! There are 2 points which seem to be important to me. The invers of Σ is actually not always computable (if there exist a single value =0) so the more nice Expression would also be Σ+ . Where Σ+ is the matrix where every non zero single value is inverted but the zeroes are left as they are. And why not follow the convention of naming matrizes? Normally a matrix is called a mxn matrix. It seems you use a nxm matrix here which i think is a bit confusing at first glance.
@drickquantum5216
@drickquantum5216 2 жыл бұрын
Thank you for the Single-Value-Matrix-Plus thing. Now I undertand the video.
@sexwax4191
@sexwax4191 2 жыл бұрын
The inverse is not defined by inverting every element, but that its multiplication with the matrix yields the identity matrix. In other words if S_inverse is the inverse of S, then S*S_inverse = 1, where 1 is the identity matrix. Therefore the inverse is computable even with zeros in the diagonal.
@subhadeepbej3241
@subhadeepbej3241 3 жыл бұрын
Excellent
@georgeyu7987
@georgeyu7987 4 жыл бұрын
i know this is just about the notation, but i think the majority of linear algebra text use m by n, rather than n by m. It's sometimes a little confusing here...
@ukaszorpik3913
@ukaszorpik3913 2 жыл бұрын
Damn, thank you, I finally got it thanks to you after 3h research xD
@deathtrick
@deathtrick 2 жыл бұрын
How does this board work? Are you writing in opposite direction?
@kangningwei934
@kangningwei934 Жыл бұрын
Thanks a lot
@surajchess3114
@surajchess3114 2 жыл бұрын
If we have x and b vectors how to find A matrix ?
@wojtekskaba9757
@wojtekskaba9757 3 жыл бұрын
Why in definition of A-dagger you put Sigma-inverse rather than Sigma-dagger? Sigma is non-square and has only pseudo-inverse rather than inverse.
@noahbarrow7979
@noahbarrow7979 2 жыл бұрын
knowing that the pseudoinverse exists makes me feel really powerful
@chengeng6472
@chengeng6472 4 жыл бұрын
Thanks for your video. However, I'm confused about the proof of the theorem(why the norm is minimized), can you give a simple proof? Thanks.
@sergiohuaman6084
@sergiohuaman6084 3 жыл бұрын
I believe it has to do with Eckart-Young theorem. you can check it out on KZfaq.
@amirkhazama7464
@amirkhazama7464 10 ай бұрын
what if we had a singular value of 0 ? ...cant that happen ? when A has dependent columns? ....in that case we wouldnt have Sigma inverse correct ? ....do we first idk,... drop some columns of A so they are all independent /.?
@uralmutlu4320
@uralmutlu4320 7 ай бұрын
I have a set of 3D positions and vectors. I try to find their intersection, so I linearized the equations like written in academic papers, but I dont get the results. I dont have a clue what could be wrong. I have Ax=b I tried: 1) (AA')^-1*A'*b 2) pinv(A)*b 3) A\b 4) I''ve just tried the SVD method, I get the same results with all. I checked my values 100 times. The vectors and directions are correct, I manually calculated them and confirmed with the matbal code.
@lukf8347
@lukf8347 3 жыл бұрын
Hello, thank you for this nice video series. This is so helpful and I use it with your book for my masters thesis. But while going through the equations, one question popped to my head: at 7:32 you use the inverse of Sigma, but for the SVD Sigma is not a quadratic matrix, rather than a nxm - Matrix and as such not invertable (in the "classic" sense of invertable matrices). So while I understand that if I use the economic SVD, this Matrix would be a mxm - Matrix, but I don't understand it for the case of nxm. Is there a video or a page in the book, where this case is discussed? Other than that, thank you very much for saving my masters degree :D
@sexwax4191
@sexwax4191 2 жыл бұрын
Sigma might be a nxm Matrix but not all rows and columns are non zero. If you remove all zero columns and rows you will get a quadratic matrix. So its essentially a quadratic matrix padded with zeros.
@Catwomen4512
@Catwomen4512 2 жыл бұрын
He replied to another person saying "Usually we will invert the first "m x m" sub-block, which is square, and then only use the first "m" columns of "U". Or, we could be even more aggressive and only invert the first "r x r" sub-block of Sigma, and only use the first "r" columns of "U" and "V", where r is much less than m." However, if you just consider when nm then we would use economy SVD (seen in his next video) so sigma would also be a square matrix.
@DiegoAndresAlvarezMarin
@DiegoAndresAlvarezMarin Жыл бұрын
I had the same question
@JoaoBarbosa-pq5pv
@JoaoBarbosa-pq5pv 3 жыл бұрын
Extremely nice lectures Steve Brunton, thank you very much for all the effort of creating and sharing them! do any of you - or anyone that reads this :) - know of any reference that explores the math of why you get min |x|2 in the underdetermined case? thank you in advance!
@Andres186000
@Andres186000 Жыл бұрын
This might be late. The reason for this min |x|2 is that any other solution would be this x hat plus something in the null-space of A. That addition would be orthogonal to this x hat and thus only be able to increase the magnitude of x hat. I am essentially reading this right out of pages 404 to 405 of Gilbert Strang's "Introduction to Linear Algebra" fourth edition.
@michaellewis7861
@michaellewis7861 4 жыл бұрын
Wouldn’t the error be ||X-x dagger||_2 not just ||x dagger||_2 in itself?
@mmarva3597
@mmarva3597 3 жыл бұрын
waw!! Merci beaucoup
@michalwisler9616
@michalwisler9616 4 жыл бұрын
Hi Sir, I am quite curious about the name of the 'transparent' blackboard, I want to buy one, where I can get it? Thank you.
@ahmaddarawshi91
@ahmaddarawshi91 4 жыл бұрын
I think this is not a "transparent blackboard". It is just glass and behind him is a black wall.
@hanyingjiang6864
@hanyingjiang6864 3 жыл бұрын
@@ahmaddarawshi91 I'm always curious how he writes. He must write all stuff in a reverse direction - I mean, b would look like d from his perspective.
@ahmaddarawshi91
@ahmaddarawshi91 3 жыл бұрын
@@hanyingjiang6864 he writes normally as he would write on a whiteboard but then the video is flipped (reflected) digitally.
@drickquantum5216
@drickquantum5216 2 жыл бұрын
I get a bit lost when he constructs the left pseudo-inverse. I cannot grasp why the single value matrix is guaranteed to have an inverse.
@julijangrajfoner1730
@julijangrajfoner1730 Ай бұрын
How can sigma be invertible when it's a rectangular matrix?
@mrgledjan3091
@mrgledjan3091 5 ай бұрын
It is not Sigma^-1 but Sigma^-T so the traspose of the inverse ?????
@Trubripes
@Trubripes 4 ай бұрын
LOL So moor-penrose pseudoinverse is just inverting one SVD term at a time, and since they U and V are orthogonal they are just transposed. Well that saved me a lot of effort looking into where it comes from.
@somethingnew7538
@somethingnew7538 3 жыл бұрын
Great!!
@puneetgupta22
@puneetgupta22 Жыл бұрын
Forget about everything what visual do you use to write and record....
@harsha123144
@harsha123144 Жыл бұрын
Hi sir, Can you solve one linear regression problem using svd and upload in KZfaq.? Please it helps me lot and others as well.!
@Eigensteve
@Eigensteve Жыл бұрын
Check the playlist, I already have an example
@harsha123144
@harsha123144 Жыл бұрын
@@Eigensteve Respected sir, can you send me the link? Sorry, for troubling...!
@michaelvandeborne9382
@michaelvandeborne9382 4 жыл бұрын
Sir, by solving for x, I understand that A dagger cannot be a true inverse of A, but since the SVD is an exact egality of A, and since multiplying by U transpose S-1 and V is again an exact operation, I don't understand where the approximation step was introduced in the calculus.
@Eigensteve
@Eigensteve 4 жыл бұрын
Great question. I address this exact question in the first 2 minutes of the next video: kzfaq.info/get/bejne/ZpiBdtd5spfVkmg.html
@michaelvandeborne9382
@michaelvandeborne9382 4 жыл бұрын
Indeed, thanks :)
@nami1540
@nami1540 2 жыл бұрын
Thanks, just had this in mind as well
@DiegoAndresAlvarezMarin
@DiegoAndresAlvarezMarin Жыл бұрын
Are you using the economy version of the SVD? otherwise you are not able to take the inverse of Sigma, which is a nxm matrix. EDIT YES >> see next video in this list
@alexyang6755
@alexyang6755 3 жыл бұрын
I FOUND GOLD.
@user-or7ji5hv8y
@user-or7ji5hv8y 3 жыл бұрын
I don’t know why such was never explained when I took econometrics.
@periklisdrakousis6537
@periklisdrakousis6537 2 жыл бұрын
why U trans U cancel to identity. Don't they result in a square matrix???
@fununterhaltung6556
@fununterhaltung6556 Жыл бұрын
Because U is a orthogonal Matrix
@arjunkc3227
@arjunkc3227 6 ай бұрын
This is cool but computational cost to determine svd is nightmare
@tanjamikovic2739
@tanjamikovic2739 Жыл бұрын
does this guy actually writes backward?
@forooghfarajzade8206
@forooghfarajzade8206 Жыл бұрын
when you say thank you i emphasize thank you!!!! (not me) :)))
@hyperduality2838
@hyperduality2838 4 жыл бұрын
Under is dual to over, left is dual to right, up is dual to down, in is dual to out. Thesis is dual to anti-thesis -- The Generalized or time independent Hegelian dialectic. Alive is dual to not alive -- Schrodinger's/Hegel's cat. Duality creates reality.
Least Squares Regression and the SVD
5:43
Steve Brunton
Рет қаралды 57 М.
Singular Value Decomposition (SVD): Mathematical Overview
12:51
Steve Brunton
Рет қаралды 382 М.
УГАДАЙ ГДЕ ПРАВИЛЬНЫЙ ЦВЕТ?😱
00:14
МЯТНАЯ ФАНТА
Рет қаралды 3,8 МЛН
9. Four Ways to Solve Least Squares Problems
49:51
MIT OpenCourseWare
Рет қаралды 117 М.
Linear Systems of Equations, Ax=b
10:04
Steve Brunton
Рет қаралды 35 М.
Fourier Series: Part 1
12:16
Steve Brunton
Рет қаралды 181 М.
Singular Value Decomposition (the SVD)
14:11
MIT OpenCourseWare
Рет қаралды 599 М.
Hardware Control with Collimator Part 1
11:29
Collimator
Рет қаралды 708
Телефон-електрошокер
0:43
RICARDO 2.0
Рет қаралды 1,3 МЛН
Cheapest gaming phone? 🤭 #miniphone #smartphone #iphone #fy
0:19
Pockify™
Рет қаралды 4,3 МЛН
$1 vs $100,000 Slow Motion Camera!
0:44
Hafu Go
Рет қаралды 27 МЛН
iPhone, Galaxy или Pixel? 😎
0:16
serg1us
Рет қаралды 1,1 МЛН