matlab去eog,使用ICA去除32通道脑电数据中的眼电

使用ICA去除32通道脑电数据中的眼电

% 去除听音乐1的脑电数据眼电

% Method: ICA

% 1.导入数据 32通道 1000hz采样率

% 2. 数据预处理(50hz陷波+0.3hz-250hz带通，可以自己改参数调整)

% 3. ICA(下载的FastICA_25)

% 4. 计算ICA的各ICs分量与原水平、垂直眼电信号相关性

% 5. 将眼电噪声的ICs置零后重建成没有眼电的脑电信号

% 注意代码中设置断点

% copyright：seu-ww

% 2019-05-07

clc;clear all;close all;

MUSIC1PATH=‘D:\Share\music_breath_2018\1-EEG\MUSIC2’;

SAVEPATH1 = ‘D:\Share\music_breath_2018\1-EEG\CLEAN_DATA\MUSIC2’;

data_list = dir(MUSIC1PATH);

eeg_number = length(data_list)-2;

fs = 1000;

for i =3:48 %45名被试的数据

tic;

filename = data_list(i).name;

datapath = [MUSIC1PATH,filename];

num_med = regexp(data_list(i).name,‘D’,‘split’);

subjectID = str2num(cell2mat(num_med(1)));

load(datapath);

data_org = data_music2