Added missing data preprocessing script

OptimizedPSOCodes/segdatacond.m

@@ -0,0 +1,46 @@
+function [whtndfiltdata, whtndstd, TFtotal]=segdatacond(segdata, PSD, sampFreq, tidxs)
+
+fmin = 30;
+% sampfreq = 4096;
+
+filtdata = segdata;
+
+datalen = length(filtdata)/sampFreq;
+
+N = length(filtdata);
+Tsig = N/sampFreq;
+timeVec = (0:N-1)*(1/sampFreq);
+
+%Highpass Filter the data and handle NaN-regions
+% [filtdata, PSD, nanchunk_start_idxs, nanchunk_end_idxs] = nanhandle(data, sampFreq, fmin, tstart, seglen, winlen);
+
+%Create Transfer Function
+TF = 1./sqrt(PSD);
+
+%Set Values of TF below fmin to 0
+TF(1:fmin*datalen) = 0;
+
+%Create entire Transfer Function vector
+negFStrt = 1-mod(N,2);
+kNyq = floor(N/2)+1;
+
+TFtotal = [TF, TF((kNyq-negFStrt):-1:2)];
+
+%Whiten Filter Data with Transfer Function
+%Window data before FFT with a Tukey-window with a 0.5 seconds rolloff
+
+rolloff = 0.5; %Roll-off in seconds
+winfiltdata = filtdata.*tukeywin(length(filtdata), rolloff*sampFreq/N)';
+
+fftfiltdata = fft(winfiltdata);
+
+whtndfftfiltdata = fftfiltdata.*TFtotal;
+
+whtndfiltdata = ifft(whtndfftfiltdata);
+
+%Divide by variance of whitened strain so that final whitened vector has
+%unit strain
+tstart = tidxs(1);
+tend = tidxs(2);
+whtndstd = std(whtndfiltdata(tstart: tend));
+whtndfiltdata = whtndfiltdata/whtndstd;