Definito N (a soli qualche mesi dalla creazione del file), sistemato problema edges e dt

%bin che dura dt

Tseg=256; %segments' length in seconds

M=fix(t(end)-t(1))/Tseg; %number of segments

M=fix((t(end)-t(1))/Tseg); %number of segments

%con fix prendo la parte intera, scarto l'ultimo segmento che tanto non

%sarà mai di lunghezza Tseg (molto improbabile)

% Time matrix with bin midpoints for each segment tm(m,j) ---------------

% tmid(m) is the midpoint in time for the m-th segment

dt = Tseg/N;

% dt = Tseg/N; %

dt = dt_psd;

N = fix(Tseg/dt);

tm = zeros(M,N);

tmid = zeros(M,1);

for m = 1:M

    tm(m,1) = t(1)+(m-1)*Tseg + dt/2;

    tm(m,1) = t(1)+(m-1)*Tseg; %t(j) è già centrato

    for j = 2:N

        tm(m,j)= tm(m,1)+(j-1)*dt; 

    end

%for each segment (lavoro su tm)

for m=1:M

    [Cm,edges]=(histcounts(x(m,:),round((tm(m,end)-tm(m,1))/dt_psd))); % R - convincitene

    edges=edges(end)-edges(2); %mi dà il tempo preciso di tutta la TdF, che sarà leggermente diversa da length(C)*dt per come è definito histcounts 

    %% La riga 137 non può essere giusta, si perde sicuramente un pezzo 

    %% di tm(1) che è preso come primo edge da histcounts: va sostituito 

    %% edges(2) con edges(1)

    edges=edges(end)-edges(1); %mi dà il tempo preciso di tutta la TdF, che sarà leggermente diversa da length(C)*dt per come è definito histcounts 

    Y=fft(Cm).'; clear Cm

    F=((0:length(Y)-1)./edges).'; clear edges

    L=length(F); %lunghezza iniziale, servirà per lo zero-padding