批量预处理实战：把流程变成脚本

% 在你做完一遍 GUI 预处理后，命令行输入：
eegh
% 会输出类似：
%   EEG = pop_eegfilt(EEG, 0.1, 45, [], []);
%   EEG = pop_runica(EEG, 'extended', 1);
%   EEG = pop_epoch(EEG, {'stimulus'}, [-0.2 0.8], 'newname', ...);
% 这就是你的脚本原料。

%% ============ 参数区：只改这里 ============
clear; clc; close all;
dataDir   = 'D:\project\raw\';      % 原始数据目录（英文路径！）
outDir    = 'D:\project\processed\';% 输出目录
qcDir     = 'D:\project\qc\';       % 质控图目录
subjects  = {'sub01','sub02','sub03','sub04'};   % 被试列表
hp = 0.1; lp = 45;                  % 滤波参数
epochWin = [-0.2 0.8];              % 分段窗口 (s)
baseline = [-200 0];                % 基线 (ms)

if ~exist(outDir, 'dir'), mkdir(outDir); end
if ~exist(qcDir, 'dir'),  mkdir(qcDir);  end

%% ============ 循环区：所有被试同一流程 ============
for s = 1:length(subjects)
    sub = subjects{s};
    fprintf('==== [%d/%d] %s ====\n', s, length(subjects), sub);
    try
        % 1. 读入
        EEG = pop_loadset('filename', [sub '_raw.set'], 'filepath', dataDir);

        % 2. 滤波（零相位）
        EEG = pop_eegfiltnew(EEG, 'fcutoff', hp, 'ftype', 'highpass');
        EEG = pop_eegfiltnew(EEG, 'fcutoff', lp, 'ftype', 'lowpass');

        % 3. ICA（长数据时耗时，可在此跳过调试）
        EEG = pop_runica(EEG, 'extended', 1);
        % TODO: ICA 成分剔除 —— 建议人工查看后把成分编号写进一个表

        % 4. 分段 + 基线
        EEG = pop_epoch(EEG, {'stimulus'}, epochWin, 'newname', [sub '_ep']);
        EEG = pop_rmbase(EEG, baseline);

        % 5. 试次级剔除（±100 uV）
        EEG = pop_rejtrend(EEG, 1, EEG.nbchan, -100, 100);

        % 6. 保存
        pop_saveset(EEG, 'filename', [sub '_pp.set'], 'filepath', outDir);

        % 7. 质控图：ERP 总平均 + 剩余试次数写进日志
        figure('Visible','off');
        plot(EEG.times, mean(EEG.data, 2));       % 通道平均的粗略 ERP
        title(sprintf('%s  trials=%d', sub, EEG.trials));
        print('-dpng', '-r150', [qcDir sub '_erp.png']);
        close;

        fprintf('    OK, %d trials kept\n', EEG.trials);
    catch err
        % 单个被试失败不中断整批
        fprintf(2, '    FAILED: %s\n', err.message);
        fid = fopen([qcDir 'error_log.txt'], 'a');
        fprintf(fid, '%s\t%s\n', sub, err.message);
        fclose(fid);
    end
end
disp('全部完成，请检查 qc 目录。');