复Morlet小波的谱图(从数据中重新创建文章结果)

时间:2013-10-14 06:53:54

标签: matlab wavelet spectrogram

我正在尝试重新创建以下结果:

Spectrogram

来自以下数据:

https://www.dropbox.com/s/mi3szqgzgku29rn/FS40.dat

时间以毫秒为单位(频率为40000 Hz)

文章指出他们使用Complex Morlet小波来创建谱图: “来自平均LFP的功率估计值来自于 通过将信号与形式的复Morlet小波进行卷积,得到1-88 Hz数据的时频谱图

w(t,f0)= Aexp(-t ^ 2/2 *σ^ 2)exp(2 * i π f0 * t)

对于每个感兴趣的频率f0,其中σ= m /2πf0,并且i是虚数单位。规范化 因子为A = 1 /(σ(2π)^ 0.5),定义时间和频率分辨率之间折衷的常数m为7.

我只是在matlab中使用频谱图功能获得了一些“好”的结果。 但我不太了解如何使用morlet复小波。 尝试将cwt与'morl'窗口一起使用时,我得到了错误的结果

谢谢。

P.S。

我正在尝试重新创建这篇文章:

由细胞型选择性光遗传驱动驱动的不同新皮层振荡的计算模型:由低阈值尖峰和快速尖峰中间神经元控制的可分离谐振电路。

1 个答案:

答案 0 :(得分:0)

我设法找到了这个MATLAB包

https://sites.google.com/site/rwfwuwx/Home/mfeeg

它与文章中的操作几乎完全相同..

包中的相关文件是: mf_tfcm.m& mf_cmorlet.m