pca滤波 matlab,MATLAB实现PCA去除眼电信号

MATLAB实现PCA去除眼电信号

一、去除伪迹的讨论

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

但一般实验会忽略掉心电伪迹(为了简化，俺也是)，若非要去除，可采用：在事件相关电位(ERP)研究中，通过设置刺激信号的随机出现，再对脑电信号叠加平均即可去除心电干扰。另外也可以通过成分分离的方法去除。

眼电的具体去除见后文，此处仅科普产生机理…

→ 眼电的产生机理：由于跨膜间存在静息电位，角膜侧为正，巩膜侧为负，这可以看作是在那里有一定大小的电偶极子存在。眼球转动过程中，偶极子的方向相对于记录电极发生改变。眼球初始位置与终止位置时的偶极子所产生的电位差，构成了眼电伪迹，眼球运动方向的电极记录到正的眼电伪迹，眼球运动相反方向的电极记录得到的眼电伪迹为负

除了上述两种信号伪迹，还有：

肌电信号EMG(解决办法：受试者尽量减少头颈部、手臂运动。加上肌电信号主要频率在100Hz以上，大部分已被之前的带通滤波去除)

呼吸(导致电极相对于头皮出现滑动)、皮肤电(被试出汗或者心理紧张所引起皮肤阻抗发生变化)

这类信号的频率较低，一般小于1Hz，脑电采集过程中将直接予以滤除。因此伪迹往往被忽略，不予考虑

非生理信号伪迹(在这些噪声中，50Hz的工频干扰尤为严重, 频谱不明显࿰