%% fast furier transofrmation %%
function [freqf,ampf,metadata]=furiertrans(cas,amp);
  
ns=length(cas);
dt=(cas(end)-cas(1))/(ns-1); 
fs=1/dt;



spektrum=fft(2*amp/ns,ns);  %furierova transformace
freqf=(fs*(0:ns/2-1)/ns)'; %% osa x - prepocet z ns na frekvenci
ampf=abs(spektrum(1:ns/2)); %% vyjadreni v absolutních hodnotách, zároven jen polovina


%% pro reverzni fft
metadata=spektrum(1:ns/2);

end