Continue to Site

Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

Spectrogram from Complex Morlet wavelet ( recreating article result from data)

Status
Not open for further replies.

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:

muMxX.jpg


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.
 

Status
Not open for further replies.

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Back
Top