function [X,freq]=positiveFFT_zero_padding(x,Fs,N)
k=0:N-1;     %create a vector from 0 to N-1
T=N/Fs;      %get the frequency interval
freq=k/T;    %create the frequency range
X=fft(x,N)/length(x); % normalize the data
cutOff = ceil(N/2); 
X = X(1:cutOff);
freq = freq(1:cutOff);