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);