q_rattio_erds.m

用于脑电分析的各类源码

function coefficient=q_rattio_erds(data,fs,ref_t,rep_t)
% 命令格式:coefficient=q_rattio_erds(data,fs,ref_t,rep_t)
% 用于在事件相关去同步化（同步化）分析中计算data序列的ERD/ERS系数
% 输入参数： data，输入数据，必须为一维数组
%           fs,采样频率
%           ref_t，参考时段时间坐标，如：[0,2],即0s到2秒
%           rep_t，响应时段时间坐标，如：[3,5]，即3秒到5秒
% 输出参数： coefficient，ERD/ERS系数，为单个数值


if min(size(data))>1 
   error('input data must be one dimension');
end

ref_index=[1+floor(fs*ref_t(1)):floor(fs*ref_t(2))];
rep_index=[1+floor(fs*rep_t(1)):floor(fs*rep_t(2))];

ref=mean(data(ref_index));
rep=mean(data(rep_index));

coefficient=(rep-ref)/abs(ref);