Changed C to Cm, few comments

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

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

%for each segments (lavoro su tm)

%for each segment (lavoro su tm)

for m=1:M

    [C,edges]=(histcounts(tm(m,:),round((tm(m,end)-tm(m,1))/dt_psd))); 

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

    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 

    Y=fft(C).'; clear C

    Y=fft(Cm).'; clear Cm

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

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

    %il secondo è un vettore colonna di zeri con lunghezza pari al numero di elementi che mi serve per tornare alla risoluzione originale

    %Fourier transform -----------------------------------------------------

    [C,edges]=(histcounts(X1,round((X1(end)-X1(1))/dt_psd))); 

    [Cm,edges]=(histcounts(X1,round((X1(end)-X1(1))/dt_psd))); 

    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 

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

    clear C

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

    clear Cm

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

    F1=F1(1:round(length(F1)/2));

    Y1=Y1(1:round(length(Y1)/2));