%% 0.!!!需手动调整参数
data_path = 'C:\Users\Amax\Desktop\basetest_flod\raw_data_flod\';
svae_path = 'C:\Users\Amax\Desktop\basetest_flod\save_fold\';
file_name_target = 'ADHD';
file_name_nontarget = 'Normal';
stuct_name =  'Data';

raw_temp_data = load ([data_path ,file_name_target,'1']);
raw_temp_data = raw_temp_data.Data.data;
fs_raw = 500;
fs_down = 250; %降采样后采样率
subject_num = [1;10];
temp_abandon = 20; %去掉的暂态，以秒为单位
epoch_length = 2;
perocess_channel = [1 ; 128];

经验建议：

采样率降到300Hz以下，目前主流研究认为，静息态脑电有效频率<100Hz，根据脑奎斯特定理，采样率到200Hz以上即可。采样率是最影响计算时间的参数之一，尤其影响后期的非线性计算，绘图、特征分析等。
建议每批处理人数<10人，如有大量被试应分为多组。一是由于多被试处理时间长，中间运算和调试周期较长；二是多被试计算结果存储量大，MATLAB对于2GB以上的数据存储较慢，并且下一阶段读取也缓慢。
去暂态20s，temp_abandon = 20，是由于一般脑电设备开机后需要稳定波形，这个数值要根据设备、现场环境修改。

3.1 数据导入

导入原始数据，并输出显示基本数据信息：

%% 1.Standard_input
% 0.1 rest_data 输入数据格式为：通道数*采样点数
standard_time = floor(size(raw_temp_data,2)/fs_raw- temp_abandon);
remain_raw_time = standard_time - temp_abandon;
remain_trial = floor(remain_raw_time/epoch_length);
remain_confrim_time = remain_trial * epoch_length;

disp(['数据总时长: ' , num2str(standard_time) , '||暂态时长: ' , num2str(temp_abandon), '||修正剩余时长: ' , num2str(remain_confrim_time)]);
disp(['采样率: ' , num2str(fs_raw), '||每段时长: ' , num2str(epoch_length), '||试次数量: ' , num2str(remain_trial)]);

3.2 前置准备-主函数

主要功能依靠调用 Standard_cut_rest 函数实现：

%% 1.切分数据
disp(['数据分段中...']);
Standard_input_target = [];
Standard_input_nontarget = [];

Standard_input_target.fs = fs_down;
Standard_input_nontarget.fs = fs_down;

Standard_input_target.subject_num = subject_num;
Standard_input_nontarget.subject_num = subject_num;

[Standard_input_target.data,Standard_input_nontarget.data] = Standard_cut_rest(data_path,stuct_name,file_name_target,file_name_nontarget,fs_raw,fs_down,subject_num,remain_trial,epoch_length,temp_abandon,perocess_channel);
Standard_input_target.file = file_name_target;
Standard_input_nontarget.file = file_name_nontarget;
disp(['||已完成标准分段||']);

3.3 前置准备-功能函数

主功能函数 Standard_cut_rest ：

function[Standard_input_target,Standard_input_nontarget] = Standard_cut_rest(data_path,stuct_name,file_name_target,file_name_nontarget,fs_raw,fs_down,subject_num,remain_trial,epoch_length,temp_abandon,perocess_channel)
%% 0.参数说明
% data_path 数据路径
% file_name_target 目标数据统一名
% file_name_nontarget 非目标数据统一名
% fs_raw 原始数据采样率
% subject_num 被试起止编号
% remain_trial 剩余试次数
% epoch_length 单试次时长
% temp_abandon 暂态去除时长
% perocess_channel 有效数据通道的首末
Standard_input_target = cell(remain_trial,subject_num(2,1)-subject_num(1,1)+1);
Standard_input_nontarget = cell(remain_trial,subject_num(2,1)-subject_num(1,1)+1);

%% 1.切分数据
print_count=0;
sub_count = 1;
for sub_loop = subject_num(1,1):subject_num(2,1)
raw_data_target = load ([data_path ,file_name_target , num2str(sub_loop)]);
raw_data_target = raw_data_target.(stuct_name).data;
raw_data_target = downsample(raw_data_target',fs_raw/fs_down)';


raw_data_nontarget = load ([data_path ,file_name_nontarget , num2str(sub_loop)]);
raw_data_nontarget = raw_data_nontarget.(stuct_name).data;
raw_data_nontarget = downsample(raw_data_nontarget',fs_raw/fs_down)';

temp_abandon_point = temp_abandon * fs_down;
for trial_loop = 1:remain_trial
    cut_temp = [];
    cut_temp = raw_data_target(perocess_channel(1,1):perocess_channel(2,1), temp_abandon_point+ (trial_loop-1)*fs_down*epoch_length+1 : temp_abandon_point + trial_loop*fs_down*epoch_length);
    Standard_input_target{trial_loop,sub_count} = cut_temp;
    
    cut_temp = [];
    cut_temp = raw_data_nontarget(:, temp_abandon_point+ (trial_loop-1)*fs_down*epoch_length+1 : temp_abandon_point + trial_loop*fs_down*epoch_length);
    Standard_input_nontarget{trial_loop,sub_count} = cut_temp;
end
% disp(['输入标准化分段: ' , num2str(sub_loop/subject_num)]);

fprintf(repmat('\b',1,print_count));   
print_count=fprintf('计算标准化分段进度 : %f',sub_count/(subject_num(2,1) - subject_num(1,1)+1));
sub_count = sub_count+1;
end
fprintf('\n');
end

3.4 前置准备-结果保存

最终，结果保存：

%% 2.保存标准输入文件
disp(['标准分段保存中...']);
save([ svae_path , 'Standard_input_target_',num2str(subject_num(1,1)),'_',num2str(subject_num(2,1))],'Standard_input_target');
save([ svae_path , 'Standard_input_nontarget_',num2str(subject_num(1,1)),'_',num2str(subject_num(2,1))],'Standard_input_nontarget');
disp(['||已完成标准分段保存||']);
t_Standard_input_cost = toc;
disp(['标准输入格式调整完毕，耗时： ',num2str(t_Standard_input_cost)]);

四、前置准备整体代码

静息态信号整体代码：

disp(['||0.静息态类数据处理-标准输入格式||']);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
%% 0.!!!需手动调整参数
data_path = 'C:\Users\Amax\Desktop\basetest_flod\raw_data_flod\';
svae_path = 'C:\Users\Amax\Desktop\basetest_flod\save_fold\';
file_name_target = 'ADHD';
file_name_nontarget = 'Normal';
stuct_name =  'Data';

raw_temp_data = load ([data_path ,file_name_target,'1']);
raw_temp_data = raw_temp_data.Data.data;
fs_raw = 500;
fs_down = 250; %降采样后采样率
subject_num = [1;10];
temp_abandon = 20; %去掉的暂态，以秒为单位
epoch_length = 2;
perocess_channel = [1 ; 128];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
tic;
%% 1.Standard_input
% 0.1 rest_data 输入数据格式为：通道数*采样点数
standard_time = floor(size(raw_temp_data,2)/fs_raw- temp_abandon);
remain_raw_time = standard_time - temp_abandon;
remain_trial = floor(remain_raw_time/epoch_length);
remain_confrim_time = remain_trial * epoch_length;

disp(['数据总时长: ' , num2str(standard_time) , '||暂态时长: ' , num2str(temp_abandon), '||修正剩余时长: ' , num2str(remain_confrim_time)]);
disp(['采样率: ' , num2str(fs_raw), '||每段时长: ' , num2str(epoch_length), '||试次数量: ' , num2str(remain_trial)]);

%% 1.切分数据
disp(['数据分段中...']);
Standard_input_target = [];
Standard_input_nontarget = [];

Standard_input_target.fs = fs_down;
Standard_input_nontarget.fs = fs_down;

Standard_input_target.subject_num = subject_num;
Standard_input_nontarget.subject_num = subject_num;

[Standard_input_target.data,Standard_input_nontarget.data] = Standard_cut_rest(data_path,stuct_name,file_name_target,file_name_nontarget,fs_raw,fs_down,subject_num,remain_trial,epoch_length,temp_abandon,perocess_channel);
Standard_input_target.file = file_name_target;
Standard_input_nontarget.file = file_name_nontarget;
disp(['||已完成标准分段||']);

%% 2.保存标准输入文件
disp(['标准分段保存中...']);
save([ svae_path , 'Standard_input_target_',num2str(subject_num(1,1)),'_',num2str(subject_num(2,1))],'Standard_input_target');
save([ svae_path , 'Standard_input_nontarget_',num2str(subject_num(1,1)),'_',num2str(subject_num(2,1))],'Standard_input_nontarget');
disp(['||已完成标准分段保存||']);
t_Standard_input_cost = toc;
disp(['标准输入格式调整完毕，耗时： ',num2str(t_Standard_input_cost)]);