oren_a
Newbie level 1
- Joined
- Oct 14, 2013
- Messages
- 1
- Helped
- 0
- Reputation
- 0
- Reaction score
- 0
- Trophy points
- 1
- Activity points
- 13
I am trying to recreate the following results:
from the following data:
https://www.dropbox.com/s/mi3szqgzgku29rn/FS40.dat
the time is in milliseconds (frequency is 40000 Hz)
The article state that they used Complex Morlet wavelet to create the spectrogram:
" Power estimates from the averaged LFPs were calculated from
time–frequency spectrograms of the data from 1–88 Hz by convolving the signals with a complex Morlet wavelet of the form
w(t,f0 ) = Aexp(−t^2 / 2*σ^2 )exp(2*i*π*f0*t)
for each frequency of interest f0, where σ = m/2πf0, and i is the imaginary unit. The normalization
factor was A = 1 /(σ (2π)^0.5 ), and the constant m defining the compromise between time and frequency resolution was 7.
I only managed to get some "good" results using spectrogram function in matlab.
But I dont have much idea of how to use the morlet complex wavelet.
I got bad result when trying to use cwt with 'morl' window
Thank You.
P.S.
I'm trying to recreate this article:
Computational modeling of distinct neocortical oscillations driven by cell-type selective optogenetic drive: separable resonant circuits controlled by low-threshold spiking and fast-spiking interneurons.
from the following data:
https://www.dropbox.com/s/mi3szqgzgku29rn/FS40.dat
the time is in milliseconds (frequency is 40000 Hz)
The article state that they used Complex Morlet wavelet to create the spectrogram:
" Power estimates from the averaged LFPs were calculated from
time–frequency spectrograms of the data from 1–88 Hz by convolving the signals with a complex Morlet wavelet of the form
w(t,f0 ) = Aexp(−t^2 / 2*σ^2 )exp(2*i*π*f0*t)
for each frequency of interest f0, where σ = m/2πf0, and i is the imaginary unit. The normalization
factor was A = 1 /(σ (2π)^0.5 ), and the constant m defining the compromise between time and frequency resolution was 7.
I only managed to get some "good" results using spectrogram function in matlab.
But I dont have much idea of how to use the morlet complex wavelet.
I got bad result when trying to use cwt with 'morl' window
Thank You.
P.S.
I'm trying to recreate this article:
Computational modeling of distinct neocortical oscillations driven by cell-type selective optogenetic drive: separable resonant circuits controlled by low-threshold spiking and fast-spiking interneurons.