The Fast Fourier Transform (FFT): Most Ingenious Algorithm Ever?

  Переглядів 1,830,102

Reducible

Reducible

День тому

In this video, we take a look at one of the most beautiful algorithms ever created: the Fast Fourier Transform (FFT). This is a tricky algorithm to understand so we take a look at it in a context that we are all familiar with: polynomial multiplication. You will see how the core ideas of the FFT can be "discovered" through asking the right questions. The key insights that are presented in this video is that polynomial multiplication can be improved significantly by multiplying polynomials in a special value representation. The challenge that presents itself is the problem of converting a polynomial from a standard coefficient representation to value representation.
We see that the FFT is an incredibly efficient recursive algorithm that performs this task, and we also discover that a slightly tweaked FFT (Inverse FFT) can also solve the reverse problem of interpolation. If this video doesn't blow your mind, I don't know what will.
0:00 Introduction
2:19 Polynomial Multiplication
3:36 Polynomial Representation
6:06 Value Representation Advantages
7:07 Polynomial Multiplication Flowchart
8:04 Polynomial Evaluation
13:49 Which Evaluation Points?
16:30 Why Nth Roots of Unity?
18:28 FFT Implementation
22:47 Interpolation and Inverse FFT
26:49 Recap
Also a subtle mistake that a commenter made me aware of -- at 26:40 instead of replacing w with (1/n * e^{-2 * pi i/ n}), the actual right way to do this is by taking the final output of the IFFT at the end of the recursion and dividing by n.
So the full change is w = e^{-2 pi i / n}
And then somewhere outside the scope of the IFFT function ifft_result = 1/n * IFFT(values)
The treatment of the FFT in this video is inspired by several well known references, mainly Introduction to Algorithms (Cormen et al.) and Algorithms (Papadimitriou et al.).
Support: / reducible
This video wouldn't be possible without the open source manim library created by 3blue1brown: github.com/3b1b/manim
Here is link to the repository that contains the code used to generate the animations in this video: github.com/nipunramk/Reducible
Elegant proof that the matrix used in the proof that (d + 1) points uniquely define a degree d polynomial is invertible: math.stackexchange.com/questi...
Music:
Lift Motif by Kevin MacLeod is licensed under a Creative Commons Attribution license (creativecommons.org/licenses/...)
Source: incompetech.com/music/royalty-...
Artist: incompetech.com/
All other music by Aakash Gandhi
SVG Attributions:
Earth: Designed by Flat Icons from www.flaticon.com, CC BY 4.0 creativecommons.org/licenses/..., via Wikimedia Commons
GPS: Icons made by www.flaticon.com/authors/pause08 from www.flaticon.com/
Wireless Comms: Design inspired by svgsilh.com/image/297434.html

КОМЕНТАРІ: 1 400
@Reducible
@Reducible 3 роки тому
For those of you who want to see me break down further how the FFT Implementation works: checkout out this video: ukposts.info/have/v-deo/jKlhfpmJbXqnyqc.html Also a subtle mistake that a commenter made me aware of -- at 26:40 instead of replacing w with (1/n * e^{-2 * pi i/ n}), the actual right way to do this is by taking the final output of the IFFT at the end of the recursion and dividing by n. Sorry about that error, and in general all bugs like this one will be updated in the description of the video. So the full change is w = e^{-2 pi i / n} And then somewhere outside the scope of the IFFT function ifft_result = 1/n * IFFT() As one additional side note, there seems to be a little confusion about the point by point multiplication at 6:31. Remember, all that's happening here is we are sampling points from polynomials A(x) and B(x) and using these points to find a set of points for C(x) = A(x)B(x). Every corresponding y coordinate for C(x) will be the y coordinate of A(x) multiplied by the y coordinate of B(x). Therefore, point by point multiplication really means "let's just keep the x-coordinates and multiply the y-coordinates."
@ashtonsmith1730
@ashtonsmith1730 3 роки тому
:) first to reply to u
@Reducible
@Reducible 3 роки тому
@@electronicmusicwave7459 Great question and the answer is that multiplication of polynomials is an indirect way to get to the DFT. The key idea is that one of the fastest ways to do polynomial multiplication is take coefficient representations of each polynomial, evaluate them at the appropriate number of roots of unity, multiply the y-coordinates in this new representation, and convert back to coefficient form. The DFT shows up in the evaluation step because evaluation of polynomials is just a matrix vector product where the matrix is the DFT. Remember, at its core, the FFT algorithm is just a fast way to multiply a vector by the DFT matrix -- polynomial multiplication just happens to be one case where the DFT shows up, so we can use the FFT. The DFT also shows up in all sorts of other applications such as taking a time domain signal and converting it to a frequency domain representation, image compression algorithms, acoustics, etc. So there's nothing special about polynomial multiplication other than the fact that the fastest way to do it involved polynomial evaluation which has the DFT.
@jameshogge
@jameshogge 3 роки тому
Was just about to comment the error but saw you already had it! Great video and really clear too!
@LeviATallaksen
@LeviATallaksen 3 роки тому
@@Reducible Divide by 2, rather than by n, I think? Otherwise, in the recursions, we would divide by n/2, while calling on recursions that divide by n/4 themselves, and so on... Really like the video, by the way!
@wormius51
@wormius51 3 роки тому
Great video. Now I got it to work. Maybe you could add a correction in the video though?
@fuckYTIDontWantToUseMyRealName
@fuckYTIDontWantToUseMyRealName 3 роки тому
"If your mind isn't blown, you weren't paying attention." I promise, I tried.
@nedisawegoyogya
@nedisawegoyogya 3 роки тому
it's too much for my 3 am brain lol
@michaelmeichtry316
@michaelmeichtry316 3 роки тому
Tried to blow your mind, or tried to pay attention? I tried both, still working on it.
@nickwallette6201
@nickwallette6201 2 роки тому
Me too! I knew I was in trouble when he started with the part "everybody" understands -- polynomial math. All the sudden I got the feeling like I had boarded a plane, and it was starting to taxi to the runway, and the flight attendant just announced that we're headed to ... wait.. where?!? Uh oh!
@scottysantos377
@scottysantos377 Рік тому
@@nickwallette6201 Same, I barely remembered what polynomials were, I last studied them was like 7 years ago when I was in Middle school.
@rafaelmenna8384
@rafaelmenna8384 Рік тому
Oh I paid attention that’s when my mind started trying to blow not because it understood the concept but because it tried so hard to understand. I think I just wrote a complicated sentence🤔
@mikemarsh5433
@mikemarsh5433 2 роки тому
Back in 1972 when I was a young physics grad working in an optics lab, I built an FFT optical spectrometer using a Michaelson interferometer using a laser reference beam as a ruler to trigger my computer program as to when to take samples of the sample beam, using an appropriate set of A/D converters. I scanned the interferometer by advancing one leg with a fine screw and we took 16,384 samples, which made for a nice resolution. Written in Fortran for a Data General Nova with a 5 pass compiler, then hand optimized the assembly code to make it efficient enough to run in real time, displayed on a long retention phosphor oscilloscope as it scanned. I even made a hardware assisted program instruction to do the bit reversal and shift to speed it up.
@TrapAddict
@TrapAddict Рік тому
That sounds so cool!
@jyun3102
@jyun3102 Рік тому
yes, further optimization by using photonics !
@7654321220
@7654321220 Рік тому
sounds like a nice experiment for undergrad students
@headyshotta5777
@headyshotta5777 Рік тому
I had to read the last sentence just to make sure this wasn't one of those overly complex "andddd I have no idea what I'm talking about. I just made all that up." Memes 🤣 don't mind me, smooth brain over here 😭
@metroboominauditorybellow563
@metroboominauditorybellow563 Рік тому
As a physics undergrad with some knowledge about FORTRAN and Assembly, and made in '72.... WOAH
@timotejbernat462
@timotejbernat462 3 роки тому
6:58 not going to lie, I spent an inordinate amount of time trying to figure out how it could possibly be O(n)! before realizing the exclamation was only punctuation, not a factorial
@edtsch
@edtsch 3 роки тому
I was thinking as I saw that: "he shouldn't have used the exclamation point there."
@moegreen3870
@moegreen3870 3 роки тому
LOL... i wonder, i think in the spanish fonts they have an exclamation that is flipped upside down :p
@Cyberplayer5
@Cyberplayer5 3 роки тому
Exclamation ! math is hard....XD
@DarthZackTheFirstI
@DarthZackTheFirstI 3 роки тому
spent a crazy amount at 17:30 because w^(n-1) (that would mean 15 points - 1 = 14/2=7; but that lacks one point to reach the opposite site) doesnt really work out when the first is 1 ... . suppose its 16 points divided by 2. then it works just fine :-p . or did i overlook some math context from before? has he uploaded a correct python code by now?
@wuxi8773
@wuxi8773 3 роки тому
Exclamation ! well deserved
@CodeParade
@CodeParade 3 роки тому
Just got recommended your channel, really like it! If you need any topic suggestions, I'd love to see a video about the disjoint set and union-find algorithm. That's one of my favorites, but it's not as well known. It even uses the inverse Ackermann function for the big-O analysis too which is super weird and interesting.
@Reducible
@Reducible 3 роки тому
Hey, didn't expect this! I absolutely LOVE the work you do CodeParade! Thanks for leaving this comment! I also love the union find algorithm. The ideas that lead up to the ultimate optimization of path compression and the Ackermann function is also quite beautiful and definitely something I'll keep in mind!
@jaikumar848
@jaikumar848 3 роки тому
@@Zcon18 which post ?
@Zcon18
@Zcon18 3 роки тому
@@jaikumar848 the community post from a couple of days ago
@jamesbalajan3850
@jamesbalajan3850 3 роки тому
I'd love to see a video on union-find. One of the most underrated data structures. Its so incredibly useful for graph problems.
@RudraSingh-pb5ls
@RudraSingh-pb5ls 3 роки тому
@@jamesbalajan3850 absolutely
@1wsm1
@1wsm1 3 роки тому
Who came here from 3blue1brown and fell in love with this channel?
@Reducible
@Reducible 3 роки тому
That shoutout was awesome!
@shadowthehedgehog2727
@shadowthehedgehog2727 3 роки тому
I didn’t know they are from 3blue1brown, but now that you mention it they remind me of one another
@nagoshi01
@nagoshi01 3 роки тому
Me! This video wasn't quite as easy to follow as Grant's videos, it was a bit too fast paced, but this guy seems to just be starting out, and he's doing a great job so far!
@JSBACH-yv9sx
@JSBACH-yv9sx 3 роки тому
Me too!
@WaluigiisthekingASmith
@WaluigiisthekingASmith 3 роки тому
@@shadowthehedgehog2727 they are different but using the same tools
@blitzkringe
@blitzkringe 2 роки тому
I think the most impressive thing about this Fourier Transform related video is that it never mentions words "sine" or "frequency".
@P-nk-m-na
@P-nk-m-na Рік тому
i mean they do mention frequency at the start but i see what you mean
@owenpenning1597
@owenpenning1597 11 місяців тому
Fourier transform and fast Fourier transform are different
@JeanBaptisteEmanuelZorg
@JeanBaptisteEmanuelZorg 11 місяців тому
or integrals hell
@ibraheemsony9947
@ibraheemsony9947 10 місяців тому
​​@@owenpenning1597 but they gave similar results
@peterhemmings2929
@peterhemmings2929 7 місяців тому
Often seems like e^jt is the more fundanmental concept, cos and sin are just the real and imaginary parts of that
@geeshta
@geeshta 2 роки тому
Years ago when I got to sound design, I was so excited that you can represent sound as both a waveform and a set of frequencies. I heard about this mysterious "Fourrier transform" that was responsible for that and the famous "FFT" algorithm that made it possible for computers to do things such as equalization quickly. I saw some videos on UKposts that explained it but lacking the proper mathematical background, I couldn't wrap my head around it. Now this video got to my recommended, with your excellent step-by-step explanation with actual examples. I still can't wrap my head around it.
@markuscwatson
@markuscwatson 2 роки тому
You need to learn the FFT from an engineer's point of view. This is more math/algorithm than you likely need.
@imlxh7126
@imlxh7126 2 роки тому
@@markuscwatson Even better, learn it from a music producer's point of view. Once you get a grip for what filters and EQ do to audio, thinking about signals in the frequency domain makes a LOT more sense.
@nickwallette6201
@nickwallette6201 2 роки тому
As an idiot, I feel like a blind man at a demonstration of optical illusions watching this video. I can tell there are really cool things happening, I just can't see them. OTOH, the 3Blue1Brown videos on FFT were inspirational. I really need to take some math classes. I love everything about audio. I like building speaker enclosures (physics!), I like coding (math!), I like building electronics, and I would love to do some DSP development. But when I start trying to grok the math behind all of that stuff, it feels like somebody took the numbers and four-function buttons off my calculator and now it's all Greek to me.
@ricardoraymond9037
@ricardoraymond9037 2 роки тому
This guy has got to be joking if he seriously thinks we got any mileage in knowledge from this video. He was just talking to himself
@charlesmrader
@charlesmrader 2 роки тому
@@imlxh7126 The reason music relates so closely to frequency is that our ears are sensitive to frequency, but not to phase. A sine wave of a given frequency sounds almost exactly the same if it is delayed by a few milliseconds. Two sinewaves of different frequency sound the same if one is shifted relative to the other. Computers that recognize speech begin with frequency analysis.
@benjaminv3748
@benjaminv3748 3 роки тому
This is a great video. I love the close interconnection of serious mathematics and efficient programming. This didn't just help me better understand tricks for efficient algorithms (the log n trick is awesome), but also understand why the fourier transform is defined as it is. I mean, I understood the mathematics since before, but this gave so much intuition. Some background knowledge defenitly helps here, I think, but I love the level. Please keep it up! Simpler concepts are covered pretty well elsewhere too, but this might be the best video on FFT on all of UKposts by quite a bit (imo). Keep it up! See you when you've totally blown up :)
@alexarnold8461
@alexarnold8461 3 роки тому
Holy crap! How has it already been 30 mins!?! You know it's a good video when you completely loose track of time!
@Reducible
@Reducible 3 роки тому
Thank you so much for this comment! I love hearing comments like this one.
@ccgarciab
@ccgarciab 3 роки тому
I didn't realize, that's amazing!
@nickikallman390
@nickikallman390 2 роки тому
My reaction was instead "Holy crap! How has it already been 3 hours!?!" because I like to pause and think things through and just marvel at the beauty of it
@1fareast14
@1fareast14 3 роки тому
I smell Grant's python library at work (oh, its credited in description, never mind)
@bimbumbamdolievori
@bimbumbamdolievori 3 роки тому
I thought soo!!!
@felixpaulus9728
@felixpaulus9728 3 роки тому
Astonishing how the quality side by side with the level of complexity of your videos keeps improving further and further. Thx.
@tienthanh8490
@tienthanh8490 2 роки тому
I’ve been curious about FFT for a while and has some rough ideas about its applications, but how it works is still a black box to me. Thanks a lot for making such useful educational content - this is the best explanation on FFT I’ve seen so far, oh and great video editing as well. I’m sure this will stick with me, and even if I need some recap I’ll go back to this video. Keep up the good work!
@dragonfyre1589
@dragonfyre1589 3 роки тому
This is a fantastic video! Feels very much like the style of 3b1b- not just bc of the animation, but also bc you have a similar way of carefully yet thoroughly explaining complex ideas. Well done!
@intrinsicload
@intrinsicload 3 роки тому
WOW. I'm literally blown away. This video is EXCELLENT. PLEASE keep doing what you're doing! Instant share.
@Reducible
@Reducible 3 роки тому
Thank you! Don't plan on stopping any time soon :)
@MadScientist267
@MadScientist267 2 роки тому
Literally?
@anatheistsopinion9974
@anatheistsopinion9974 3 роки тому
Beyond excellent! Now I must spend my day binging your entire content. Can't wait for more already!
@Superwatson28
@Superwatson28 2 роки тому
Incredible the way you explain such a complex topic. I could follow all of the details. Great job!
@Jellylamps
@Jellylamps 3 роки тому
I’m just getting into programming and i love learning why things work the way they do. I feel like this channel is going to be very helpful to me
@skoky76
@skoky76 Рік тому
It won’t - watching this as an common programmer is like if as a worker you would watch how shovel is made in a forge. Nice to watch and know how it is done, but won’t do the shoveling 😀
@petergilliam4005
@petergilliam4005 3 роки тому
Honestly among my favorite youtube videos ever created. I've never felt such an urge to like a video. Every time I'd come back, I'd be disappointed that I've already liked the video and can't do it again.
@hellNo116
@hellNo116 3 роки тому
ok i remember being blown away the first time i learned about fft and i only saw it only the context of signal processing. this is truly amazing. i would watch the big O notation if i didn't just have a month doing so already this semester. your video was amazing and i am really glad i came across it. really inspiring keep up the good work!!!
@CrypticPulsar
@CrypticPulsar Місяць тому
ive parked on this topic for the past 9 weeks, and only now thanks to this incredible video the pieces are finally coming to fall into place! Thank you so much!
@abc-by1kb
@abc-by1kb 3 роки тому
13:00 for those of you who are wondering why he is plugging in the square of x_i. It's because when we use the coefficient form, if we just plug in x_i, it would be interpreted as powers 0,1,2,...., so if we plug in x^2 it would be interpreted as (x^2)^0. (x^2)^1, (x^2)^2, (x^2)^3, which are powers 0,2,4,6...which are even terms.
@vylbird8014
@vylbird8014 3 роки тому
An aside: When you mention that any polynomial of degree N can be uniquely determined by N+1 points? That's how many forms of forward error correction work, suitably modified for modular arithmetic. If you want to send your message you break it up into N parts, each of which becomes a point on a polynomial of degree N-1. Then generate as many extra points on the polynomial as you need to generate a list of K points where K>N, and transmit those K points through your unreliable communications channel. As long as at least N of the K make it through to the other end, you can recover the original message.... and thus your phone call can keep on going even with a little radio interference taking place. You can tune K to be optimal for the reliability of your communications channel and the acceptable level of overhead.
@Reducible
@Reducible 3 роки тому
You've just discovered a future video I'm planning on error correcting codes (Shh don't tell anyone :P).
@AJMansfield1
@AJMansfield1 3 роки тому
@@Reducible Convolutional error correcting codes too, pls?
@Reducible
@Reducible 3 роки тому
​@@AJMansfield1 the ones that I was thinking about going through initially are Reed-Solomon codes and the Berlekamp Welsh algorithm. Both of those techniques rely on similar ideas with polynomials but over a finite field. Pretty cool stuff. Convolutional error correction codes are also pretty interesting so may think about addressing that as well. I have an insanely long list of video ideas that I'm just adding to one by one based on request and my interest so I can't make any promises :)
@edwardfanboy
@edwardfanboy 3 роки тому
Shamir's secret sharing algorithm is based on that same property.
@gabelluc9573
@gabelluc9573 3 роки тому
This is absolutely brilliant !
@abc3631
@abc3631 3 роки тому
What an awesome presentation .. also love the fact that you share your code for people to gain the knowledge you have. So hats off 👍
@xXKILLBOT27Xx
@xXKILLBOT27Xx Рік тому
I like the fact you actually laid out the pseudocode for the algorithm at the end of each section, most channels that cover topics like these usually leave it at the block diagram. Shows you really know your stuff.
@BernardWu97
@BernardWu97 3 роки тому
I absolutely love this. Complex math being presented in an intuitive way makes it so much more accessible. Keep up the good work!
@Maltanx
@Maltanx 3 роки тому
Dude, you uploaded this video like two minutes before I searched for a video on how FFT works. Not even subbed but this was one of the first videos to pop up. That was a ridiculously perfect timing! Also, I'm getting strong 3blue1brown vibes from your channel, subscribed right away!
@Reducible
@Reducible 3 роки тому
Haha that's awesome! Always a great time to learn about the FFT!
@l.f.3835
@l.f.3835 3 роки тому
I never knew Applejack's brother was into algorithm theory!! Awesome :^)
@avi7278
@avi7278 2 роки тому
I think that has something to do with polynomial multiplication or something, ehhhemmm...
@viharivemuri7202
@viharivemuri7202 11 місяців тому
I just read about the FFT from CLRS, where the authors directly went into choice of evaluation points as complex roots of unity. But your intuition and the way you deduced that is brilliant! Thanks man!
@aloshe1000
@aloshe1000 2 роки тому
This is the most elegant way to explain FFT I have ever encountered. Brilliant!
@fdmeneses
@fdmeneses 3 роки тому
It's awesome to see 3b1b inspiring great new channels like Reducible. Fantastic material thus far and I'm truly excited for what's to come. Well done!
@GameBoyzCZ
@GameBoyzCZ 3 роки тому
2 minutes in and i love it already
@histravelstories8458
@histravelstories8458 Місяць тому
This might be one of the first videos other than travel vlogs where I fell in love with it and subscribed. FFT had been a nightmare, everytime I started learning it mind would just wander off. this video is soooooo good. Thank you for making this. And being beautiful ... you have added beauty to the FFT algorithm.
@renatobritto4096
@renatobritto4096 2 роки тому
Wow just wow, this has gotta be the best video of the year for me. Absolutely astounding work!
@shishka4022
@shishka4022 3 роки тому
Man, this is one of the best videos I’ve ever seen! You are very clear and keep it very interesting. Serious props to you, great job!
@MrKyltpzyxm
@MrKyltpzyxm 3 роки тому
This popped up in my recommendations. I have a bad habit of watching long, intricately detailed descriptions of concepts that I do not have the foundational education to understand. This is probably why every time I see a video like this that assumes I have some prior knowledge, I wonder how I ended up here. The video itself, however, is spectacular. Very easy to follow, even with my lack of understanding. (I don't think I am describing this very well, but you'll have to trust me.) I work as a cable technician and have a very VERY basic understanding of RF signal processing theory. It's always interesting to peer into the depths of what keeps it all tied together. I keep telling myself that one day I'll work up the nerve to pursue a more structured, formal education in these topics, but in the mean time I am glad that such accessible presentations exist out here in the wild.
@Reducible
@Reducible 3 роки тому
Thank you so much for taking the time to write out such a nice comment! I try to make the videos as accessible as possible, but it can be tough especially when we are talking about something as advanced as algorithms to calculate fourier transforms. But I'm glad you were able to get something out of it!
@savinolupo2379
@savinolupo2379 2 роки тому
Thank you so much for the beautiful arguments you used to show the magnificence of the FFT algorithm!
@MurakDurak
@MurakDurak 3 роки тому
i had this exact algorithm in a numerical programming class in university and my professor did a fairly poor job at explaining what FFT and IFFT really are - he just showed us the algorithm and how its recursive call works. this video makes me truly appreciate how ingenious this algorithm actually is, and knowing where the algorithm comes from makes understanding it a lot easier. thank you for this amazing explanation!
@bghnv123
@bghnv123 3 роки тому
dude, this channel is really well made! im definitely subbing for more
@MarekKnapek
@MarekKnapek 3 роки тому
26:44 "So, if your mind isn't blown, you haven't been paying attention." Not gonna lie, you had me in the first half.
@zhenhaogu629
@zhenhaogu629 3 роки тому
So clearly explained and well demonstrated! Thank you so much for your work.
@dhruvakashyap
@dhruvakashyap 3 роки тому
I'm glad I found this channel, it's really beautiful. Good luck!!
@AndrewAkaHrun
@AndrewAkaHrun 3 роки тому
Amazing! I somehow missed this part of signal processing course so was using FFT as a black wonderbox almost everyday. Now I have a clear view on how it's acually working. Thanks a lot, man! Your work is highly appreciated. Special thanks for using Python for code representation.
@NEMountainG
@NEMountainG 3 роки тому
One of the first videos on UKposts I’ve seen that intuitively implements a complex mathematical idea in code. I appreciate how must respect you have for your viewers by showing them the steps to solve this problem without “baby talk”. Absolutely fantastic work! Also, I am aware of 3B1B, Primer, etc. Those channels are also fantastic but I enjoyed how this video actually showed the code implemented from scratch.
@Reducible
@Reducible 3 роки тому
Thank you so much for this awesome comment!
@NEMountainG
@NEMountainG 3 роки тому
You deserve all the compliments you’re reading. I’m fairly familiar with this topic and you did a great job explaining it, easily better than how it was first explained to me.
@eric3813
@eric3813 Рік тому
Damn, i just found that Video and got SOO MUCH MORE insight in the magic behind FFT! Thank you so much!!
@chapmag6578
@chapmag6578 3 роки тому
Great video. Nicolet launched the 440 FFT dual channel analyser in the mid 70’s whilst I was half way through my eng degree. It was lust at first sire from me, and lead me to a wonderfull career in signal analysis and machine condition monitoring.
@mohammedbelgoumri
@mohammedbelgoumri 3 роки тому
One of the absolute best math videos I ever watched. Only problem is I couldn't find the code for this particular video in the repo.
@Reducible
@Reducible 3 роки тому
Apologies, I made a rookie mistake and forgot to do that before uploading this video -- just pushed the code for this video.
@mohammedbelgoumri
@mohammedbelgoumri 3 роки тому
Thank you kindly. Keep up the amazing work!
@jajssblue
@jajssblue 3 роки тому
Agreed. Felt like a masterclass. I wish I had this video back in college when I covered the material.
@shantanurahman5084
@shantanurahman5084 3 роки тому
You can find the code (with explanation too) over here cp-algorithms.com/algebra/fft.html
@mohammedbelgoumri
@mohammedbelgoumri 3 роки тому
@@shantanurahman5084 Thanks a lot. I meant the code for the animation, not that of the FFT.
@muskyoxes
@muskyoxes 3 роки тому
this is a ridiculously awesome amount of work, possibly hundreds of times as long as the video. it hurts to think this will make less money than someone doing a makeup tutorial
@RotBaron
@RotBaron 3 роки тому
the creator of these animations has a super loaded patreon page and other stuff. he is doing great!
@davidgustavsson4000
@davidgustavsson4000 2 роки тому
No reason to shit on other creators though. A lot of experience and planning can go into a makeup tutorial too.
@matthewlee6922
@matthewlee6922 Рік тому
Insanely brilliant idea separated into simple concepts! Great video.
@mukhamediyar
@mukhamediyar Рік тому
This is single-handedly the best explanation on the topic I have ever seen. Perfect for college students!!
@dombowombo3076
@dombowombo3076 3 роки тому
I think 3Blue1Brown is very happy to see that his library is put to a good use. You even addopted his style of presentation. :) Nice work.
@benprice9036
@benprice9036 3 роки тому
Thank you so much for an excellent presentation on a beautiful and complex topic. I was going to ask what inspired the dissection using polynomials, but you already left references in the description! Absolutely fantastic!
@farrukhzamir
@farrukhzamir 3 місяці тому
Beautifully described and explained. This is the Best video on FFT on youtube which describes the underlying concepts of FFT.
@STD282
@STD282 3 роки тому
Absolutely beautiful, I always wanted to understand as to why FFT works! The FFT has always been presented to me as a thing that can map vector convolution to vector multiplication element-wise, and now it looks astonishingly clear why. The trick is, the convolution (or as I sometimes call it, long multiplication without carrying) is *exactly* what you should do in order to multiply the polynomials in coefficient representation, therefore, you can apply it anywhere where convolution takes place. Thank you for this video!
@desmonddulaney7131
@desmonddulaney7131 Рік тому
Very helpful video. However, I'd like to point out that @10:50, the decomposed functions should be P_e(x) and P_o(x) not P_e(x^2) and P_o(x^2) because once we do substitute in x^2, it should get us back to the original function. No one has pointed out this mistake so I figured I would because I was confused for a little bit.
@jingyuanxu7272
@jingyuanxu7272 3 місяці тому
For anyone has the same confusion - It's actually correct. You start with a function P(x), which in the example, is something to the x^5. You x is input, the function P(x) is something of x^5. To decompose, you split this x^5 into two parts, in the video example, it's x^2, x^4 and odd part is x^3 and x^5, except x^3 and x^5 can be expressed as x*x^2 and x*x^4. The two decomposed function, P_e and P_o, is a function of x, or you can think it as function of y, where y=x^2. So essentially, 2x^4+7x^2+1 is P_e(y) = 2y^2 + 7y + 1. So it's not a mistake, x^2 just mean P_e(y) will take input where the input is x^2 from the original function. FYI - All the video content basically came from the book called "Algorithms". I think combining the book with the video makes a lot more sense...coming from somebody reading the book, don't understand a thing, had to watch the video, then video makes a high level sense but omits a lot of details so had to go back to book again.
@dannyshaw6979
@dannyshaw6979 3 роки тому
Wow, I am very pleased to see such substantial inspiring video showing up here, thank you.
@Number_Cruncher
@Number_Cruncher 3 роки тому
Thank you very much. Very impressive!!! I've been working with FFT and did my own implementation. But I never thought, that it could be that simple to understand. And thank you for sharing the code for the video, this will be helpful to try the manim library on my own.
@hannesk2119
@hannesk2119 3 роки тому
I am glad 3b1b gave you a shoutout so that I got to see this! very nice video, I especially liked seeing the actual code implementation and of course the smooth animations :)
@JoJoModding
@JoJoModding 3 роки тому
Did 3b1b? Seems like I missed that. Where did he give that shoutout?
@sainathsingineedi2922
@sainathsingineedi2922 3 роки тому
@@JoJoModding community post
@tsunghan_yu
@tsunghan_yu 2 роки тому
@@sainathsingineedi2922 Patreon?
@tsunghan_yu
@tsunghan_yu 2 роки тому
@@JoJoModding ukposts.info/have/v-deo/n2iDn3epn6p0xGg.html
@miklov
@miklov 3 роки тому
My mind wasn't blown even once but I didn't arm it first. I'll rewatch this when I have read up on roots of unity and DFT matrices. I did enjoy the content even though it was currently way over my head ^^
@sathviknallapuri4667
@sathviknallapuri4667 Рік тому
Surely, one of the best explanations of FFT I have ever listened to.
@davidvose2475
@davidvose2475 2 роки тому
Beautifully constructed tutorial, well done. I learned a lot!
@sdsa007
@sdsa007 2 роки тому
This is such an amazingly lucid visual explanation of Fourier! I think Laplace, himself would almost be proud, if not infuriated. So much to visualize, but so little time! I mean frequency…
@Abrissbirne66
@Abrissbirne66 Рік тому
infourierted
@user-ki3sb8st5p
@user-ki3sb8st5p 3 роки тому
Mind-blowing... Never do I think of FFT when speaking about polynomial multiplication, nor do I ever see so clear the core of the algorithm. My first encounter with FFT is horrible, with the overwhelming notations making me dizzy. This channel absolutely deserves more views.
@Patrick-li8nh
@Patrick-li8nh 3 роки тому
Love it! Thank you for this amazing presentation, really enjoyed it!
@fizqialfairuz5744
@fizqialfairuz5744 3 роки тому
Im so glad that I found this channel, keep making the good stuff!.
@TrapAddict
@TrapAddict Рік тому
The best way to understand the FFT algorithm is to do it several times for different signal types by hand, starting off with simple signals. Get the equations, draw the diagrams, and then get the outputs and even write a script to do it if you're inclined. Sometimes you just have to follow the rules and steps, and trust that it will all make sense over time!
@arno.claude
@arno.claude Рік тому
I really wanted to understand FFT, so I took the time to thoroughly watch 3Blue1Brown's, Veritasium's, and your video on the topic, and your approach helped me the most! Maybe because I watched it last, but still, you did a phenomenal job. Well done!
@justarandomcatwithmoustache
@justarandomcatwithmoustache Рік тому
same man.
@zihengtong4816
@zihengtong4816 3 роки тому
Sir, you've explained this topic really clearly! And luckily I've seen your video 3days before my exam. Thank you sooo much!!
@andresquesnel1350
@andresquesnel1350 4 місяці тому
This is one of the best explanations I have ever seen. Thank you
@PowerhouseCell
@PowerhouseCell 3 роки тому
Really neat video! Speaking from experience, this must have taken forever to make haha. Brilliantly explained and produced :)
@Reducible
@Reducible 3 роки тому
I just watched some of your videos! Man, you are going to explode -- very well done. I can definitely learn a few things from you :) Keep at it!
@jakobv8
@jakobv8 3 роки тому
I don't even have a layman's understanding of complex mathematics, but I was very pleased to be along for the ride, just out of curiosity. Beautiful presentation and I hope to be able to revise and understand more one day.
@valthorhalldorsson9300
@valthorhalldorsson9300 2 роки тому
Loved this video - had to pause a few times but the explanation really clicked
@marcoronzani7197
@marcoronzani7197 2 роки тому
You managed to make me fall in love with the FFT ahahaha! I already knew it’s use cases, theoretical bases and all, but the explanation you provided is simply perfect!
@elael2
@elael2 3 роки тому
Great video!! I learned FFT in the classical context, but this is amazing! (Thanks 3b1b for suggesting 😄) Just one detail, I feel we cannot apply the 1/n directly to omega (if the inverse is right). Because 1/n is applied to the 1 values on the matrix as well as the omega is raised to some power that would affect the 1/n too. =)
@Reducible
@Reducible 3 роки тому
Yeah another astute commenter like yourself noticed this, I updated my pinned comment with this note to address this. Sorry about that!
@JohnWick-xd5zu
@JohnWick-xd5zu 3 роки тому
This guy was my TA at Cal ❤️❤️
@Reducible
@Reducible 3 роки тому
Go bears! :)
@electronicsandroboticsclub750
@electronicsandroboticsclub750 3 роки тому
John Wick is at caltech now?
@isoEH
@isoEH 3 роки тому
Brings back memories. Thank you! Nice tribute and explanation..
@catalinmihit
@catalinmihit 3 роки тому
This was not only interesting, it was really well explained AND pleasantly animated. You did more in 30 mins than my teacher did in a year, subbed xp
@anywallsocket
@anywallsocket 3 роки тому
This is insanely good work. I wish I had a better sense of how Euler’s identity, complex numbers, polynomials, and the cycle of differentiation and integration all relate, because everything seems to boil down to these ideas at the heart of mathematics.
@peterfireflylund
@peterfireflylund Рік тому
Add linear algebra to the list! (And Generating Functions - I think there is a 3blue1brown video about them.)
@fredoo6627
@fredoo6627 3 роки тому
This video is both useful and beautiful.
@gushinets
@gushinets 3 роки тому
Just brilliant explaination!! For the past several months i have been working on FFT vectorized implementation and digged into this algo, but still i couldn`t quite understand it. Thank you for this video!
@joshgorton3843
@joshgorton3843 3 роки тому
I'm a software developer with a background in sound engineering... so I'm use to peering into "black boxes" only to find the math is way over my head... I just trust the magic, LOL. THIS video was about to be more of the same were it not for the final wrap up at the end - WOW! Well done, sir. Well done!
@capableText3
@capableText3 3 роки тому
This is awesome ❤️
@natefidalgo7625
@natefidalgo7625 2 роки тому
Most videos or how people would approach introducing FFT would be by means of Fourier transforms and complex numbers. Including myself and i liked this way of explaining it since it was independent of knowing Fourier transforms or complex analysis to any significant amount. Very very neat way to expose people to FFT thru the problem of fast polynomial multiplication. Just so happened i came across your video and realized it was so close to what i was trying to do with Schönhage-Strassen algorithm and Fürer's_algorithm to get a Harvey and van der Hoeven algorithm O(nlogn). Still got some work to do but was a delight to see somebody creating a video so well on a similar problem. Keep it up great channel!
@thomasdeloof9720
@thomasdeloof9720 10 місяців тому
Great video ! i just discovered your chanel and i am very impressed by the quality of your content. Thanks !
@MasamuneX
@MasamuneX 3 роки тому
this is way over my head but I appreciate things being talked about in terms of computability not just in the abstract
@DJTimeLock
@DJTimeLock 3 роки тому
"We'll take a look at something you're all familiar with. Polynomial multiplication." 1:53 minutes in and you already lost me lmao
@communistpropagandist4608
@communistpropagandist4608 2 роки тому
bruh
@leptok3736
@leptok3736 2 роки тому
We've gone too deep 😂
@lorenpearson1230
@lorenpearson1230 2 роки тому
I used to teach this, and I always found a concrete example was helpful to students. While FFT applies anywhere, we were working with time based signals and extracting the frequencies. While you showed that at the start, I feel you should finish with that. The classic teaching approach of telling the class where you will be taking them, taking them there, and summing up by revisiting where they've been.
@chrisb198
@chrisb198 3 роки тому
Pretty much the best explanation out there for this topic... Thank you!
@brendanmitchell4581
@brendanmitchell4581 2 роки тому
Really awesome video. Right on my level, understandable but challenging enough to make me think. My mind was kinda blown a couple times! Really well done! For me, the one extra thing I'd love to see is an explanation of the connection between this algorithm and the arithmetic or discrete Fourier Transform. I think of a Fourier transform as a method for viewing data or functions from the time domain in the frequency domain. This video outlined a method for evaluating polynomials at the roots of unity, but it isn't obvious to me that that is equivalent to a transformation to frequency space. At first glance, they seem like completely different problems. After some work I was able to convince myself that they were equivalent (and experimentation confirms it) but I would love to see a follow up video to solidify that connection in my mind! Again, awesome work and thanks for a great video!
@batner
@batner 3 роки тому
13:30 Thats it i'm lost. Good video man! Will have to review from top.
@Leon-be4lx
@Leon-be4lx 3 роки тому
batner Same Here
@MrNerdHair
@MrNerdHair 2 місяці тому
Excellent. I was looking for info on the NTT and iNTT over finite fields, and this is a perfect explanation that doesn't incorporate any of the "motivating" nonsense about waves or rotations that aren't relevant for my particular situation!
@nachopuma13
@nachopuma13 3 роки тому
Immediate sub! This is AWESOME, keep it up. Watching you from Argentina 🇦🇷🇦🇷
@icicle8334
@icicle8334 3 роки тому
You can never tell how UKposts search algorithm leads you to the right place
@benmaghsoodi2067
@benmaghsoodi2067 2 роки тому
It uses FFT
@amicloud_yt
@amicloud_yt 3 роки тому
holy crap wow this video really made things click for me i've never gotten goosebumps from understanding something before but i just did
@airxperimentboom
@airxperimentboom 3 роки тому
This is beautiful. Thank you for this astounding work :D
@zacontraption
@zacontraption 3 роки тому
Thanks for taking the time to make this video
@Renslay
@Renslay 3 роки тому
I never went into understanding FFT, because it always seemed so scary to me - even after I understood the Fourier transformation itself. But holy math, this explanation is brilliant! The animations are so smooth (3B1B did a really amazing job with his module!), and everything is clean and understandable. The length of the final code is completely mind blowing! I got so hyped I wrote my own Python implementation, which contains the entire polynomial multiplication, using the FFT implementation you provided (with a few minor modifications). So here it is! Notes: - in FFTCore(), sign decides whenever we use FFT (sign=1), or IFFT(sign=-1) - using list comprehensions instead of for cycles - expanding the polynomials so we always operate on length of powers of 2 - pretty printing the result (the limitations of float values bring near-zero but non-zero imaginary components after IFFT) import math, cmath def FFTCore(P, sign): n = len(P) if n == 1: return P ye, yo = FFTCore(P[::2], sign), FFTCore(P[1::2], sign) w = [cmath.exp(sign*2*cmath.pi*1j*i/n) for i in range(n//2)] y1 = [ye[j] + w[j]*yo[j] for j in range(n//2)] y2 = [ye[j] - w[j]*yo[j] for j in range(n//2)] return y1 + y2 # concatenation def FFT(P, inverse=False): n = len(P) assert((n & (n-1) == 0) and n != 0), "Length of P must be a power of 2" if inverse: return list(map(lambda x: x/n, FFTCore(P, -1))) return FFTCore(P, 1) def polyMul(A, B): lenA, lenB = len(A), len(B) lenC = pow(2, math.ceil(math.log2(lenA+lenB))) A, B = A + [0]*(lenC-lenA), B + [0]*(lenC-lenB) Af, Bf = FFT(A), FFT(B) Cf = [Af[i] * Bf[i] for i in range(lenC)] return FFT(Cf, inverse=True) def prettyPrint(C, lenC): res = [f"{x.real:0.2f}" if cmath.isclose(x.imag, 0.0, abs_tol=0.0001) else f"{x.real:0.2f}{['+',''][x.imag < 0]}{x.imag:0.2f}j" for x in C] print(', '.join(res[:lenC-1])) A = [2, 3, 1] B = [1, 0, 2] C = polyMul(A, B) prettyPrint(C, len(A)+len(B))
@diegoG21
@diegoG21 3 роки тому
I was blown by the hability to make me lose track of the topic faster than any teacher, engineer and mathematician combined.
@nedisawegoyogya
@nedisawegoyogya 3 роки тому
sarcasm?
@abdullahsy7072
@abdullahsy7072 3 роки тому
I really enjoyed the explanation and how elegant it is. Great video :)
@lyhed5234
@lyhed5234 Місяць тому
Incredible video and explanation! My Mind is blown and it'll probably be that again when I actually understand each step
@GeofreySanders
@GeofreySanders 3 роки тому
This... this is a perfect example of FM technology.
The Discrete Fourier Transform: Most Important Algorithm Ever?
29:24
But what is the Fourier Transform?  A visual introduction.
20:57
3Blue1Brown
Переглядів 10 млн
Підставка для яєць
00:37
Afinka
Переглядів 118 тис.
I PUT MY ARMOR ON (Creeper) (PG Version)
00:19
Sam Green
Переглядів 6 млн
The Fourier Series and Fourier Transform Demystified
14:48
Up and Atom
Переглядів 758 тис.
FFT Example: Unraveling the Recursion
7:39
Reducible
Переглядів 79 тис.
3. Divide & Conquer: FFT
1:20:52
MIT OpenCourseWare
Переглядів 304 тис.
¿Qué es la Transformada de Fourier? Una introducción visual
19:48
3Blue1Brown Español
Переглядів 225 тис.
Understanding the Discrete Fourier Transform and the FFT
19:20
MATLAB
Переглядів 100 тис.
Denoising Data with FFT [Python]
10:03
Steve Brunton
Переглядів 164 тис.
The Discrete Fourier Transform (DFT)
17:36
Steve Brunton
Переглядів 322 тис.
But what is a convolution?
23:01
3Blue1Brown
Переглядів 2,4 млн
La destacable historia detrás del algoritmo más importante de todos los tiempos
24:20