-
Notifications
You must be signed in to change notification settings - Fork 0
/
GetFDandPowerSpectrum.m
34 lines (25 loc) · 1.04 KB
/
GetFDandPowerSpectrum.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
% Get Jenkinson's FD and Power Spectrum
% Define variables
Fs = 0.754; % This is the TR
T = 1/Fs;
numVols = 616;
subject = 'sub-023';
cd ~/Desktop
mov = dlmread('prefiltered_func_data_mcf.par');
mov = mov(:,[4:6,1:3]);
% ------------------------------------------------------------------------------
% Compute fd (Jenkinson2002) (used by Satterthwaite2012)
% ------------------------------------------------------------------------------
% Get FD
fdJenk = GetFDJenk(mov, 50);
% ------------------------------------------------------------------------------
% Get frequency
% ------------------------------------------------------------------------------
[f P1 P2 freqA psdxA freqB psdxB] = PowerSpect(Fs,fdJenk);
% ------------------------------------------------------------------------------
% Get plots
% ------------------------------------------------------------------------------
PlotPowerSpect(subject,numVols,fdJenk,f,P1,P2,freqA,psdxA,freqB,psdxB);
fig = gcf;
set(fig,'PaperPositionMode','Auto')
print(fig,[subject,'_powerspect.bmp'],'-dbmp')