% ------------------------------
% Matlab csv數(shù)據(jù)批量預(yù)處理輸出至Excel
% 時(shí)間：2022年8月28日
% 所需軟件代碼：Matlab R2020a
% 本代碼版本：3.00
% 程序員：wishing
% ------------------------------
%% 初始化
% 清理運(yùn)行環(huán)境
clearvars
clc
%% 單個(gè)文件處理
%% 目標(biāo)：
% 1.篩選出每個(gè)被試4種速度條件下的數(shù)據(jù)（無(wú)應(yīng)答率还惠，正確率泞莉，反應(yīng)時(shí)，最快10%平均反應(yīng)時(shí)）,需計(jì)算出結(jié)果寞秃，最終每個(gè)被試有四中速度下的結(jié)果瑰枫，每種速度下有四個(gè)值踱葛。
% 2.整理成RW狀態(tài)下speed1（+-0.5）/speed2（+-1）/speed3(+-1.5)/speed4(+-2)下的所有被試的數(shù)據(jù)。（RW下四個(gè)Sheet,SD下四個(gè)sheet光坝，每個(gè)sheet下每個(gè)被試有四個(gè)值)

%% 循環(huán)處理
for No = 10 : 37
        filename_RW = string(append('D:\desk\speed\data_1\behavior\RW','\','1',string(No),'.csv'));%文件位置尸诽，制作NO循環(huán)
        filename_SD = string(append('D:\desk\speed\data_1\behavior\SD','\','2',string(No),'.csv'));%文件位置，制作NO循環(huán)
       %% 數(shù)據(jù)準(zhǔn)備
        data_RW = readmatrix(filename_RW); % 導(dǎo)入成矩陣
        data_SD = readmatrix(filename_SD);
        A = data_RW(:,[2 18 19 20]);%提取2,18,19,20列組成新數(shù)組盯另，分別代表speed性含，res_keys,res_corr,res_rt
        B = data_SD(:,[2 18 19 20]);
        R_speed_1 = A(A(:,1)==7.5|A(:,1)==8.5,:);% 篩選出speed1,2,3,4的數(shù)據(jù)
        R_speed_2 = A(A(:,1)==7|A(:,1)==9,:);
        R_speed_3 = A(A(:,1)==6.5|A(:,1)==9.5,:);
        R_speed_4 = A(A(:,1)==6|A(:,1)==10,:);
        S_speed_1 = B(B(:,1)==7.5|B(:,1)==8.5,:);
        S_speed_2 = B(B(:,1)==7|B(:,1)==9,:);
        S_speed_3 = B(B(:,1)==6.5|B(:,1)==9.5,:);
        S_speed_4 = B(B(:,1)==6|B(:,1)==10,:);
         %% 計(jì)算指標(biāo)
         %RW
         R_speed_1_norate = numel(find(isnan(R_speed_1(:,2))))/14;%計(jì)算無(wú)應(yīng)答率，nan個(gè)數(shù)/14
         R_speed_1_corate = sum(R_speed_1(:,3))/14;%計(jì)算正確率鸳惯，1求和/14
         R_speed_1_meanrt = mean(R_speed_1(:,4),'omitnan');%計(jì)算均值商蕴，忽略nan
         R_speed_1_fastrt = min(R_speed_1(:,4));%最快10%，14個(gè)中就是最小值
         R_speed_2_norate = numel(find(isnan(R_speed_2(:,2))))/14;
         R_speed_2_corate = sum(R_speed_2(:,3))/14;
         R_speed_2_meanrt = mean(R_speed_2(:,4),'omitnan');
         R_speed_2_fastrt = min(R_speed_2(:,4));
         R_speed_3_norate = numel(find(isnan(R_speed_3(:,2))))/14;
         R_speed_3_corate = sum(R_speed_3(:,3))/14;
         R_speed_3_meanrt = mean(R_speed_3(:,4),'omitnan');
         R_speed_3_fastrt = min(R_speed_3(:,4));
         R_speed_4_norate = numel(find(isnan(R_speed_4(:,2))))/14;
         R_speed_4_corate = sum(R_speed_4(:,3))/14;
         R_speed_4_meanrt = mean(R_speed_4(:,4),'omitnan');
         R_speed_4_fastrt = min(R_speed_4(:,4));
         %SD
         S_speed_1_norate = numel(find(isnan(S_speed_1(:,2))))/14;
         S_speed_1_corate = sum(S_speed_1(:,3))/14;
         S_speed_1_meanrt = mean(S_speed_1(:,4),'omitnan'); 
         S_speed_1_fastrt = min(S_speed_1(:,4));
         S_speed_2_norate = numel(find(isnan(S_speed_2(:,2))))/14;
         S_speed_2_corate = sum(S_speed_2(:,3))/14;
         S_speed_2_meanrt = mean(S_speed_2(:,4),'omitnan');
         S_speed_2_fastrt = min(S_speed_2(:,4));
         S_speed_3_norate = numel(find(isnan(S_speed_3(:,2))))/14;
         S_speed_3_corate = sum(S_speed_3(:,3))/14;
         S_speed_3_meanrt = mean(S_speed_3(:,4),'omitnan');
         S_speed_3_fastrt = min(S_speed_3(:,4));
         S_speed_4_norate = numel(find(isnan(S_speed_4(:,2))))/14;
         S_speed_4_corate = sum(S_speed_4(:,3))/14;
         S_speed_4_meanrt = mean(S_speed_4(:,4),'omitnan');
         S_speed_4_fastrt = min(S_speed_4(:,4));     
        %% 輸出結(jié)果
     % RW
      writematrix([R_speed_1_norate R_speed_1_corate R_speed_1_meanrt R_speed_1_fastrt],'RW_speed.xlsx','sheet',1,'Range',string(append('A',string(No),':D',string(No))));%循環(huán)寫入表格
      writematrix([R_speed_2_norate R_speed_2_corate R_speed_2_meanrt R_speed_2_fastrt],'RW_speed.xlsx','sheet',2,'Range',string(append('A',string(No),':D',string(No))));
      writematrix([R_speed_3_norate R_speed_3_corate R_speed_3_meanrt R_speed_3_fastrt],'RW_speed.xlsx','sheet',3,'Range',string(append('A',string(No),':D',string(No))));
      writematrix([R_speed_4_norate R_speed_4_corate R_speed_4_meanrt R_speed_4_fastrt],'RW_speed.xlsx','sheet',4,'Range',string(append('A',string(No),':D',string(No))));
     % SD
      writematrix([S_speed_1_norate S_speed_1_corate S_speed_1_meanrt S_speed_1_fastrt],'SD_speed.xlsx','sheet',1,'Range',string(append('A',string(No),':D',string(No))));
      writematrix([S_speed_2_norate S_speed_2_corate S_speed_2_meanrt S_speed_2_fastrt],'SD_speed.xlsx','sheet',2,'Range',string(append('A',string(No),':D',string(No))));
      writematrix([S_speed_3_norate S_speed_3_corate S_speed_3_meanrt S_speed_3_fastrt],'SD_speed.xlsx','sheet',3,'Range',string(append('A',string(No),':D',string(No))));
      writematrix([S_speed_4_norate S_speed_4_corate S_speed_4_meanrt S_speed_4_fastrt],'SD_speed.xlsx','sheet',4,'Range',string(append('A',string(No),':D',string(No))));
end