%close all;
clear all;
%Procura o ponto mais proximo ja sabendo qual o arquivo ou segmente

NomeDoArquivo = 'C://Matlab//TOOLBOX//SVD_data//a2.sig';
NomeDoArquivop = 'C://Matlab//TOOLBOX//SVD_data//a2.pos';
NomeDoArquivoout = 'C://Matlab//TOOLBOX//SVD_data//centralcontact.txt';

NPoint = 50;
NumeroDeIteracaoNoArq = 200;
%NumeroDeIteracaoNoArq = 1200;

%shapes=zeros(NumeroDeIteracaoNoArq,3*NPoint);
%Abri o arquivo
fid  = fopen(NomeDoArquivo,'r');
fidp  = fopen(NomeDoArquivop,'r');
fidout = fopen(NomeDoArquivoout,'w');

%Shaper para o sinal onde ocorreu a interacao
SizeDelay = 20;

for Iteracao=1:NumeroDeIteracaoNoArq
	count = 18*50;

	[a count] = fscanf(fid,'%f\n', count);
%	[b connt2] = fscanf(fidp, '%f %f %f', 3);
	%organizo o Dado
	for k=0:(count/50)-1
		Dado(k+1,:) = a(1+(k*50):50+(k*50))';
	end

	%Shaper para o sinal onde ocorreu a interacao	
	DadoModificado = zeros(1,NPoint);
	DadoModificado(1:SizeDelay) = Dado(8,1:SizeDelay);
	DadoModificado(SizeDelay+1:NPoint) = Dado(8,SizeDelay+1:NPoint)-Dado(8,1:NPoint-SizeDelay);
    for ponto = 1:50
        fprintf(fidout,'%f ',Dado(9,ponto));
        %fprintf(fidout,'%f ',DadoModificado(ponto));
    end

n=[40 -40];
d=[41 -39];
%Dado(8,:)=dlsim(n,d,[Dado(8,:) ones(1,25)]);
shapes(Iteracao,:)=[Dado(7,:) zeros(1,25) dlsim(n,d,[Dado(8,:) ones(1,25)]) Dado(9,:) zeros(1,25)];

%shapes(Iteracao,:)=[Dado(7,:) zeros(1,25) DadoModificado(1,:) zeros(1,25)  Dado(9,:)  zeros(1,25) ];
%shapes(Iteracao,:)=[Dado(7,:) DadoModificado(1,:) Dado(9,:)];
    fprintf(fidout,'\n');
end

%fecha os arquivos
fclose(fid);
fclose(fidp);
fclose(fidout);

%remove funny shapes
shape(1:7,:)=shapes(1:7,:);
shape(8:28,:)=shapes(9:29,:);
shape(29:107,:)=shapes(33:111,:);
shape(108:189,:)=shapes(119:200,:);


clear Dado ans SizeDelay         
clear DadoModificado                count                         
clear Iteracao                      fid                           
clear NPoint                        fidout                        
clear NomeDoArquivo                 fidp                          
clear NomeDoArquivoout              k                             
clear NomeDoArquivop                ponto 
clear NumeroDeIteracaoNoA a shapes d n


