我正在尝试重新创建以下结果:
来自以下数据:
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。
我正在尝试重新创建这篇文章:
由细胞型选择性光遗传驱动驱动的不同新皮层振荡的计算模型:由低阈值尖峰和快速尖峰中间神经元控制的可分离谐振电路。
答案 0 :(得分:0)
我设法找到了这个MATLAB包
https://sites.google.com/site/rwfwuwx/Home/mfeeg
它与文章中的操作几乎完全相同..
包中的相关文件是: mf_tfcm.m& mf_cmorlet.m