generatePowerSpectrav2

get time vectors

0001 function [dataFreqRoach]=generatePowerSpectrav2(time,data)
0002 
0003 %get time vectors
0004 timeroach=time;
0005 lengthvec = length(data)
0006 Fs=1./(time(10)-time(9))
0007 
0008 datamean=data-mean(data);
0009 
0010 
0011 
0012 
0013 
0014 
0015 %Fs=100; %sampling rate of the data
0016 %FsTemps = 1; %sampling rate of the temps
0017 T = 1/Fs;                     % Sample time
0018 L = length(data);                     % Length of signal
0019 t = (0:L-1)*T; 
0020 
0021 
0022 NFFT = 2^nextpow2(L); % Next power of 2 from length of y
0023 
0024 f = Fs/2*linspace(0,1,NFFT/2+1);
0025 
0026 
0027 Yd = fft(datamean,NFFT)/L;
0028 
0029 
0030 
0031 Ydb = 2*abs(Yd(1:NFFT/2+1));
0032 
0033 
0034 bin_size=100;
0035 %YRRbin=bin(YRRb,bin_size);
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 dataFreqRoach=[f' Ydb];