本文介绍了EM算法在MATLAB环境下的具体实现方法,并应用于SSMT(二次短时傅里叶变换)时频分析中,展示了该算法的有效性和应用价值。
该存储库包含发表于PNAS的论文“状态空间multitaper时频分析”中的算法实现代码。如果这些代码对您的研究有所帮助,请引用以下论文:S.-E. Kim,M. Behr,D. Ba 和 E.N.Brown,“状态空间多锥时频分析”,美国国家科学院院刊,第1卷。115,No. 1,pp.E5-E14,2018年。
代码文件包括:
- main.m:主要代码
- EM_parameters.m:使用EM算法计算噪声和状态方差的函数
- periodogram.m:用于计算周期图的函数
- multitaper.m:用于计算多峰频谱图的函数
- SS_ST.m:用于计算SS周期图的函数
- SS_MT.m:用于计算SS多谱图的函数
说明:将所有代码下载到同一目录中,并运行main.m文件。该代码会显示在全身麻醉下记录的EEG数据中的四个频谱图。