MATLAB实现PCA去除眼电信号
一、去除伪迹的讨论
查阅文献[1]可知,所获取的低频脑电信号中,主要受心电(ECG)和眼电(EOG)的干扰较大。
但一般实验会忽略掉心电伪迹(为了简化,俺也是),若非要去除,可采用:在事件相关电位(ERP)研究中,通过设置刺激信号的随机出现,再对脑电信号叠加平均即可去除心电干扰。另外也可以通过成分分离的方法去除。
眼电的具体去除见后文,此处仅科普产生机理…
→ 眼电的产生机理:由于跨膜间存在静息电位,角膜侧为正,巩膜侧为负,这可以看作是在那里有一定大小的电偶极子存在。眼球转动过程中,偶极子的方向相对于记录电极发生改变。眼球初始位置与终止位置时的偶极子所产生的电位差,构成了眼电伪迹,眼球运动方向的电极记录到正的眼电伪迹,眼球运动相反方向的电极记录得到的眼电伪迹为负
除了上述两种信号伪迹,还有:
肌电信号EMG(解决办法:受试者尽量减少头颈部、手臂运动。加上肌电信号主要频率在100Hz以上,大部分已被之前的带通滤波去除)
呼吸(导致电极相对于头皮出现滑动)、皮肤电(被试出汗或者心理紧张所引起皮肤阻抗发生变化)
这类信号的频率较低,一般小于1Hz,脑电采集过程中将直接予以滤除。因此伪迹往往被忽略,不予考虑
非生理信号伪迹(在这些噪声中,50Hz的工频干扰尤为严重, 频谱不明显