% Apos verificar todos os sujeitos, observei que o menor t inicial ocorreu em 150 pontos, por isso o offset foi estabelecido como a media dos
100 primeiros pontos.
% Dada Freq. Coleta 1000 Hz e Tempo Coleta 4 seg
z(:, t) = z(:, t) - z_off; x(:, t) = x(:, t) - x_off;
% [a, b] = butter(ORDEM, Freq_Nyquist = Freq_Corte/(Freq_Amostragem/2) ); %Freq_Corte=200HZ; Freq_Amostragem=1000HZ
[a, b] = butter(2, fc/(freq/2));
xf = filtfilt( a, b, x(:, t) );
zf = filtfilt( a, b, z(:, t) );
% =================
% GRAFICOS INICIAIS
% =================
scr = get(0,'ScreenSize'); % Vetor com resolucao da tela [1 1 1024 768]
h1 = figure('Position',[1 1 scr(3) 0.9*scr(4)]); % Maximiza janela do grafico
h1 = plot( time, z(:, t), 'r.'); hold on; h1=plot(time, z(:, t), 'k-');
axis([min(time) max(time) min(z(:, t)) max(z(:, t))]);
xlabel('tempo (s)');ylabel('FZ (PC)');
title([label ' - Trial: ' char(tentativa(t))]);
grid;
clc;
disp('Marcacao dos ptos para zoom (Com MOUSE)');
disp('--------------------------------------');
disp('1o pto: antes do inicio do apoio');
disp('2o pto: logo apos fim do apoio');
disp(' ' );
% definicao do zoom
[xg1 yg1] = ginput(2);
%1o pto - antes do inicio do apoio
%2o pto - logo apos fim do apoio
axis([xg1(1) xg1(2) min(z(:,t)) yg1(2)]);
% Determina INICIO APOIO
ok = 1;
% while ok == 1;
while not(isempty(ok));
[tg2 zg2] = ginput(1);
i_find = find(time > tg2); i_1 = i_find(1);
i = (i_1 - 1) + find( z (i_1:length(z), t ) > zg2 );
inic_apoio = i(1) - 1 ;
h1 = plot( time(inic_apoio), z(inic_apoio, t), 'ro');
set(h1, 'markersize',6,'linewidth',3);
disp('Confirmacao da Marcacao do Inicio do Apoio');
disp(' ');
ok = input('Ponto do Inicio do Apoio ok? [S] ou N = 1 > ');
%if ok == 1; h1 = plot( time(inic_apoio), z(inic_apoio, t), 'ko'); set(h1, 'markersize',6,'linewidth',3); end;
if not(isempty(ok)); h1 = plot( time(inic_apoio), z(inic_apoio, t), 'wo'); set(h1, 'markersize',6,'linewidth',3); end;
end;
disp(' ' );
disp('Confirmacao da Marcacao do FIM do Apoio');
disp(' ');
% Determina FIM APOIO
ok = 1;
while not(isempty(ok));
[tg2 zg2] = ginput(1);
i_find = find(time > tg2); i_1 = i_find(1);
i = (i_1 - 1) + find( z(i_1:length(z), t) < zg2 );
fim_apoio = i(1);
h1 = plot( time(fim_apoio), z(fim_apoio, t), 'ro');
set(h1, 'markersize',6,'linewidth',3);
ok = input('Ponto do Fim do Apoio ok? [S] ou N = 1 > ');
%if ok == 1; h1 = plot( time(fim_apoio), z(fim_apoio, t), 'ko'); set(h1, 'markersize',6,'linewidth',3); end;
if not(isempty(ok)); h1 = plot( time(fim_apoio), z(fim_apoio, t), 'wo'); set(h1, 'markersize',6,'linewidth',3); end;
end;
inic_fim(t, 1) =inic_apoio;
inic_fim(t, 2) =fim_apoio;
t_apoio(t,1) = (inic_fim(t, 2) - inic_fim(t, 1))/1000;
%%%%% inicio da interplacao inicial %%%%%%%%%%%%%%%%
%Interpolacao do Apoio
t101 = 0:100; t101 = t101'; % tempo de apoio normalizado (0% a 100% passo 1%)