MATLAB实现PCA去除眼电信号

MATLAB实现PCA去除眼电信号

一、去除伪迹的讨论

查阅文献[1]可知，所获取的低频脑电信号中，主要受心电（ECG）和眼电（EOG）的干扰较大。

• 但一般实验会忽略掉心电伪迹（为了简化，俺也是），若非要去除，可采用：在事件相关电位（ERP）研究中，通过设置刺激信号的随机出现，再对脑电信号叠加平均即可去除心电干扰。另外也可以通过成分分离的方法去除。
• 眼电的具体去除见后文，此处仅科普产生机理…
  → 眼电的产生机理：由于跨膜间存在静息电位，角膜侧为正，巩膜侧为负，这可以看作是在那里有一定大小的电偶极子存在。眼球转动过程中，偶极子的方向相对于记录电极发生改变。眼球初始位置与终止位置时的偶极子所产生的电位差，构成了眼电伪迹，眼球运动方向的电极记录到正的眼电伪迹，眼球运动相反方向的电极记录得到的眼电伪迹为负
  
  除了上述两种信号伪迹，还有：
• 肌电信号EMG（解决办法：受试者尽量减少头颈部、手臂运动。加上肌电信号主要频率在100Hz以上，大部分已被之前的带通滤波去除）
• 呼吸（导致电极相对于头皮出现滑动）、皮肤电（被试出汗或者心理紧张所引起皮肤阻抗发生变化）