matlab代码对一eeg单通道数据进行去眼电

以下是一个简单的Matlab代码，用于去除EEG单通道数据中的眼电：

% 导入EEG数据
data = load('eeg_data.mat');
eeg_data = data.eeg_data;

% 设定参数
fs = 256; % 采样频率
fc = 30; % 截止频率
order = 4; % 滤波器阶数

% 设计滤波器
[b, a] = butter(order, fc/(fs/2), 'high');

% 应用滤波器
filtered_data = filter(b, a, eeg_data);

% 提取眼电信号
eye_data = filtered_data - eeg_data;

% 对眼电信号进行去除
window_size = 256; % 窗口大小
step_size = 128; % 步长

for i = 1:step_size:length(eye_data)-window_size
    window = eye_data(i:i+window_size-1);
    % 求均值
    mean_window = mean(window);
    % 减去均值
    eye_data(i:i+window_size-1) = window - mean_window;
end

% 重新组合EEG和眼电信号
cleaned_data = filtered_data - eye_data;

该代码首先导入EEG数据，然后设定一些参数，如采样频率、截止频率和滤波器阶数。然后，它使用Butterworth滤波器对EEG数据进行高通滤波，以去除低频噪声和基线漂移。接下来，它从滤波后的EEG数据中提取眼电信号，并在每个窗口中减去窗口的均值，以去除眼动信号。最后，它重新组合滤波后的EEG数据和清理后的眼电信号，以获得去眼电的EEG数据。