【附代码】SSVEP解码算法 - 多变量同步指数（MSI）

1 算法来源

1. 该算法来自电子科技大学张杨松博士，针对该算法的计算在张博士的博士论文中有详细介绍，有兴趣的读者可以下载阅读，本文重点在对该方法的代码实现。
2. Zhang, Yangsong, et al. “Multivariate synchronization index for frequency recognition of SSVEP-based brain–computer interface.” Journal of neuroscience methods 221 (2014): 32-40.
3. DOI：10.1016/j.jneumeth.2013.07.018

2 算法原理

1. 准备好一段大小为N×M的脑电信号X，参考信号Y的大小为2Nh×M。其中N是EEG信号的通道数，M是信号的采样点数。Nh是参考信号中的正弦和余弦的谐波数。
   
   其中Fs为采样率，M是信号采样点数（采样率×时间(秒)）。

2. 对X和Y进行标准化处理，目的是为了不失一般性。但本人测试了一下，是否做标准化的差别不太大。

3. 计算X和Y之间的联合相关矩阵：
   
   

   <