dtau unico

%Fourier transform on original time-series --------------------------------

%per mantenere l'informazione di fase, non faccio il valore assoluto al quadrato della fft

%for each segment (lavoro su tm)

dtau = zeros(1,N);

taul = zeros(N+1,1);

% dtau = zeros(1,N);

% taul = zeros(N+1,1);

taul = zeros(2,1);

disp(M);

for m=1:M

            % tic

            tau = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttempdif(1:N)).^((1:s_s).').'),2);

            taul(1:N) = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttempdifl(1:N)).^((1:s_s).').'),2);

            taul(N+1) = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttemp(N)-tmid(m)+0.5*dt).^((1:s_s).').'));

            for j=1:N

                dtau(j) = ((taul(j+1)-taul(j))); 

            end

            % taul(1:N) = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttempdifl(1:N)).^((1:s_s).').'),2);

            % taul(N+1) = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttemp(N)-tmid(m)+0.5*dt).^((1:s_s).').'));

            % for j=1:N

            %     dtau(j) = ((taul(j+1)-taul(j))); 

            % end

            taul(1:2) = tmid(m) + sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((ttempdifl(1:2)).^((1:s_s).').'),2);

            dtau = ((taul(2)-taul(1))); 

            % toc

            % Y1=(2./sum(X1).*abs(fft(X1)).^2).'; %normalizzazione Leahy, giusto?????

            % tic

            % Y1 = fft(interp1(ttemp(:),xtemp(:),tau,'linear',0)).';

            Y1 = fft((interp1(tau(:),(xtemp(:)./dtau(:)),ttemp(:),'linear',0)).*dt).';

            % Y1 = fft((interp1(tau(:),(xtemp(:)./dtau(:)),ttemp(:),'linear',0)).*dt).';

            Y1 = fft((interp1(tau(:),(xtemp(:)./dtau),ttemp(:),'linear',0)).*dt).';

            % clear X1

            %F1=((0:length(Y1)-1)./(ttemp(end)-ttemp(1))).'; 

            F1=((0:length(Y1)-1)./(length(Y1)*dt_psd)).';