clear all
close all
x=dlmread('sw00426438000brtms_t1.txt',';');
tf=x(:,1);
fs=1/(tf(2,1)-tf(1,1));
%input trgtime
trgtime=299725433.152  ;
t=tf-trgtime;
idx0=find(t==0);
tt=t(idx0-500:idx0+500);
s_i=x(idx0-500:idx0+500, 2);
figure; plot(tt, s_i);
L=floor(length(s_i)/2); %%%%%%%%%%% Lunghezza finestra per ssa : pari alla metà della lunghezza della serie temporale
%
I=[1:12,L]; %%%%%%%%%%%%%%%Punto delicato per ottenere curva di luce filtrata da noise
[y,sr,vr]=ssa(s_i,L,I);