Hi Erick,

a package that comprises both, beat detection from PPG and HRV analysis function, should be this one: https://physionet.org/physiotools/physionet-cardiovascular-signal-toolbox/

Greetings, Sebastian

Hi Sebastian,

Thank you very much!

Hi

The methodology (code) is not different from HRV analyses having RR-interval from ECG.

fs = 4;

time = (Rtimes(1):1/fs:Rtimes(end));

IBIint = interp1(Rtimes,IBIs,time','spline')';

y = detrend(IBIint, 'linear'); % you can apply different detrending methods

window = min(300*fs, length(y));

noverlap = window/2;

nfft = max(256,2^nextpow2(window));

if rem(nfft, 2) == 0

DFT = (nfft/2)+1;

else

DFT = (nfft+1)/2;

end

% you can decide which DFT you want for example DFT = 0:0.001:2;

% for AR metod

order = 20; % you can use different order

[Pyy, f] = pburg(y,order,DFT,fs, 'onesided');

% for FFT method

[Pyy, f] = pwelch(y,window,noverlap,DFT,fs,'PSD','onesided');

VLF_ab = sum(Pyy(f >= 0.003 & f < 0.04));

LF_ab = sum(Pyy(f >= 0.04 & f < 0.15));

HF_ab = sum(Pyy(f >= 0.15 & f < 0.4));