Convolution as spectral multiplication

  Рет қаралды 16,501

Mike X Cohen

Mike X Cohen

Күн бұрын

This video lesson is part of a complete course on neuroscience time series analyses.
The full course includes
- over 47 hours of video instruction
- lots and lots of MATLAB exercises and problem sets
- access to a dedicated Q&A forum.
You can find out more here:
www.udemy.com/course/solved-c...
For more online courses about programming, data analysis, linear algebra, and statistics, see
sincxpress.com/

Пікірлер: 23
@sumanmighty
@sumanmighty 3 жыл бұрын
Wow . Just wow. This is such a great explanation.
@mikexcohen1
@mikexcohen1 3 жыл бұрын
Glad it was helpful!
@addyad_yt
@addyad_yt 2 жыл бұрын
Hi, Thanks, nice illustrations and examples. Can you make a video about deconvolution of 1D signals? For example, if we have the observed signal and the input kernel, how to identify the system function (reflectivity series). It is really hard to understand this concept because it is an inverse problem.
@shubhamnayak5398
@shubhamnayak5398 4 жыл бұрын
How to make the signal that you were using to explain morlet wavelet as spectral filter or similar random signal that are combination of different frequencies and are not stationary in matlab?
@jiongwang7645
@jiongwang7645 2 жыл бұрын
such an elegant explanation! thanks a lot!
@mikexcohen1
@mikexcohen1 2 жыл бұрын
Glad you liked it!
@romanvereb7144
@romanvereb7144 4 жыл бұрын
Spectacular
@saeedmirzaee8315
@saeedmirzaee8315 2 жыл бұрын
Thank you sir ,keep going and make your powerful and fantastic totourial
@mikexcohen1
@mikexcohen1 2 жыл бұрын
Thank you, I will :)
@mandychan8282
@mandychan8282 3 жыл бұрын
Thank you very much for your wonderful illustration. I have a question about IFFT, hoping that you could explain further. As you mentioned the transformation of frequency spectrum back into time-domain using IFFT can result in the same time series pattern. How is the sequence conveyed in the frequency spectrum? Could it be possible that, after IFFT, the frequencies are all there but they are in different order? Thank you.
@mikexcohen1
@mikexcohen1 3 жыл бұрын
Hi Mandy. As long as the Fourier coefficients aren't shuffled, the reconstruction in the time domain will be accurate. I have some videos on IFFT (e.g., #12 in NEW ANTS2, and probably one in the OLD ANTS series) that I hope will help clarify this point.
@mandychan8282
@mandychan8282 3 жыл бұрын
@@mikexcohen1 I see. Thanks for the guide!! Your tutorials are fantastic and I'm looking forward to go through them!
@wwmheat
@wwmheat Жыл бұрын
Thank you Mike for the great video! Im struggling to wrap my head around one thing though.. At this panel that starts at around 11:00 we see the signal, kernel, result of convolution and the image of spectral multiplication. When I go through it from top to bottom, through convolution, it seems perfectly intuitive. We have this signal changing with time, and the convolutions with the wavelet of certain frequency helps us to clearly uncover, let's say, this "blob" of that frequency starting after approx. 2/3 of the orange time series (Result). So, there is some temporal dynamics, and the convolution helps to emphasize it. That's okay. However, I can't get how the inverse fourier transform (when mentally going from 'frequency' timeseries up to 'Result') could recover the same blob at the exact time? I mean, how the frequency spectrum, multiplied by that gaussian image of wavelet in frequency domain - which makes it very localized to short range of frequencies - how this frequency spectrum after it is deciphered back to time domain can still hold the information about that exact timing of that blob on the orange line? Thanks
@mikexcohen1
@mikexcohen1 Жыл бұрын
Hi Valeriy. The images I show are simplified. I only show the amplitude spectrum for simplicity and interpretability. But actually, the multiplication is with the complex-valued spectrum, which means it includes the phases. The phase information is where the temporal dynamics are stored.
@bokkieyeung504
@bokkieyeung504 3 жыл бұрын
Thanks for the video! at 17:00, what do you mean by saying "that length is n+m-1"? I think n+m-1 is the length of the convolution result. but for spectra multiplication, I think the frequency resolution matters? although the length of signal (n)/kernel (m) will affect frequency resolution. and I can't understand why IFFT will get a result of length n+m-1? I think n+m-1 is originated from the "time-series sliding with zero-padding" procedure, which is the convolution in time domain, i.e., the top panel of your comparison plot, but 5-steps here are related to the bottom panel, right? I know they are equivalent, but I expect with the procedure done in frequency domain, there is no need for zero-padding (thus no need to cut of "wings") as we did in the time-series sliding procedure? but seems that you also mentioned if the signal and kernel have different lengths, we will zero-padding when we do FFTs for them (separately). the purpose for this is to get the same frequency resolution in order to multiply spectra? I know you emphasized that zero-padding in the time-series sliding vs. FFT are different things, but I get confused... sorry for the long question list
@mikexcohen1
@mikexcohen1 3 жыл бұрын
Yes, that's how it works. The zero-padding is necessary to get the spectra to match perfectly in frequency resolution. Imagine if the signal were 1000 points long and the kernel were 20 points long -- you can't multiply their spectra unless you zero-pad. The n+m-1 number is to get it to match with the time-domain convolution. If you only zero-pad to max(m,n), you get what's called "circular convolution" (n+m-1 is linear convolution). That's also a valid operation, but the edge effects wrap around the signal, so artifacts from the beginning of time appear at the end. Linear convolution pushes those edge effects into the "wings," so they get cut off. That's why it's more commonly done.
@bokkieyeung504
@bokkieyeung504 3 жыл бұрын
thanks for the detailed explanation. You mentioned "match" between time and frequency domain. So (n+m-1) will be valid for both the "time-series sliding with zero-padding" procedure and the "separate FFT-multiplication-IFFT" procedure.
@brendawilliams8062
@brendawilliams8062 Жыл бұрын
Thankyou
@addyad_yt
@addyad_yt 2 жыл бұрын
what is the optimal amplitude of the kernel should be? after convolution, my signal's amplitude becomes very high, in my case more than my kernel.
@mikexcohen1
@mikexcohen1 2 жыл бұрын
Normalization is always tricky in signal processing. There are often different normalizations that achieve different goals. What I like to do is max-value normalize the wavelet in the frequency domain. So you scale the wavelet spectrum such that the maximum amplitude gain is 1. That normalizes for unit gain at the peak frequency.
@addyad_yt
@addyad_yt 2 жыл бұрын
@@mikexcohen1 thanks!
@tobi3497
@tobi3497 4 жыл бұрын
Why do you 0-pad m on boths sides, instead of m/2 on both sides. If you 0 padded only m/2 on both sides, wouldn't this mean you wouldn't require step 5 (cutting off the wings)
@mikexcohen1
@mikexcohen1 4 жыл бұрын
Good observation. When doing convolution in the time domain, you can indeed start the loop at m/2. But to get the results to match with spectral multiplication, you need the full resolution.
Complex Morlet wavelet convolution
12:44
Mike X Cohen
Рет қаралды 18 М.
How AI Discovered a Faster Matrix Multiplication Algorithm
13:00
Quanta Magazine
Рет қаралды 1,4 МЛН
The day of the sea 🌊 🤣❤️ #demariki
00:22
Demariki
Рет қаралды 100 МЛН
The filter-Hilbert method
23:07
Mike X Cohen
Рет қаралды 22 М.
Convolution in the time domain
23:37
Mike X Cohen
Рет қаралды 25 М.
Parameters of Morlet wavelet (time-frequency trade-off)
18:19
Mike X Cohen
Рет қаралды 14 М.
The Man Who Solved the World’s Hardest Math Problem
11:14
Newsthink
Рет қаралды 292 М.
But what is the Fourier Transform?  A visual introduction.
20:57
3Blue1Brown
Рет қаралды 10 МЛН
A simple guide to chaos theory - BBC World Service
5:10
BBC World Service
Рет қаралды 39 М.
The multi-taper method
11:04
Mike X Cohen
Рет қаралды 10 М.
Is graphene starting to live up to its hype?
28:03
RAZOR Science Show
Рет қаралды 168 М.
Wavelets: a mathematical microscope
34:29
Artem Kirsanov
Рет қаралды 609 М.
Main filter..
0:15
CikoYt
Рет қаралды 14 МЛН
Собери ПК и Получи 10,000₽
1:00
build monsters
Рет қаралды 2,1 МЛН
Мой инст: denkiselef. Как забрать телефон через экран.
0:54
После ввода кода - протирайте панель
0:18
Up Your Brains
Рет қаралды 1 МЛН
Худший продукт Apple
0:53
Rozetked
Рет қаралды 146 М.